METHOD FOR MEASURING VISCOELASTIC MODULUS OF SUBSTANCE, AND APPARATUS FOR MEASURING VISCOELASTIC MODULUS OF SUBSTANCE

Abstract
[Problem] A method for measuring a viscoelastic modulus of a substance and an apparatus for measuring the viscoelastic modulus of the substance are provided for allowing information on viscoelasticity of an adsorption substance to be expressed by moduli G′ and G″ which are generally used when expressing viscoelasticity, and for further allowing calculation of the viscoelastic modulus in real time.
Description
TECHNICAL FIELD

The present invention relates to a method for measuring a viscoelastic modulus of a substance and an apparatus for measuring a viscoelastic modulus of a substance, using a sensor for adsorption measurement and physical property evaluation of a substance in a solution in the field of chemistry, physics, biochemistry, pharmacy, material, or the like.


BACKGROUND ART

A relationship between frequency changes and mass loads caused by adsorption of QCM is expressed by Sauerbrey equation indicated by the following equation 14.










Δ





Fs

=


-


2


Nf
0
2





ρ
Q



μ
Q








Δ





m

A






Equation





14







However, in the case of measurement in the field of chemistry, biochemistry, pharmacy, or the like in which the measurement is mainly made in a solution, since the solution is contained in an adsorption substance unlike in the atmosphere, the adsorption substance is known to be a film with viscoelastic properties. A frequency change of a resonance frequency Fs measured in conventional QCM has therefore indicated a value that includes a mass load due to adsorption, a viscous load of the solution, and a viscoelastic effect of the adsorption substance itself, and these three elements could not be separated.


Accordingly, the present applicant found from the invention described in Patent Document 1 that the three elements included in a frequency variation of the above Fs can be separated and respectively calculated as frequency variations.


However, even if conversion into mass can be made using the above Sauerbrey equation from a frequency change calculated by a separated mass load, a viscoelastic modulus of an adsorption substance cannot be calculated from a separated viscoelastic element.


Specifically, it was found from the invention described in Patent Document 1 that the three elements included in a frequency variation of the above Fs can be separated and respectively calculated as frequency variations. However, the mutual separation of amass load term, a viscous load term, and a viscoelastic term cannot be made in all of measurements, and there was a case that a calculated viscous change term was not zero even if there was no viscous change. Furthermore, each separated term was a value of a frequency, and a viscoelastic modulus could not be calculated.


CITATION LIST
Patent Documents
Patent Document 1: JP-A-2007-10519
SUMMARY OF THE INVENTION
Problems to be Solved by the Invention

Accordingly, an object of the present invention is to provide a method for measuring a viscoelastic modulus of a substance that allows information on viscoelasticity of an adsorption substance to be expressed by moduli G′ and G″ which are generally used when expressing viscoelasticity, and further allows calculation of the viscoelastic modulus in real time. The present invention also has an object to calculate the viscoelasticity loss coefficient, rigidity coefficient, and viscosity coefficient of the adsorption substance from the obtained viscoelastic modulus.


Means of Solving the Problems

In order to solve the above problems, a first solution means in the present invention is characterized in that, in a system for forming a film by adsorbing a substance to the surface of a piezoelectric element or to a film fixed onto the piezoelectric element in a solution, using a sensor with the piezoelectric element submerged on both sides or one side in the solution, at least two of N-th waves (N=1, 3, 5 . . . (N=2n+1)) of the piezoelectric element are used, and further at least two of a resonance frequency Fs, and half-value frequencies F1 and F2 (F2>F1) having half conductance values of a conductance value of the resonance frequency in each N-th wave are used to calculate a mass load term, a viscoelastic term (1), a viscoelastic term (2), and a viscoelastic term (3) represented by the following equations 1 to 3:











Δ






F
2


=



-


f
0


π






z
q





ω






ρ
1



h
1


+



f
0


π






z
q






(


G


+

G



)




G


2




ω
2



ρ
2



η
2



h
1














Mass





load





Viscoelastic





term






(
1
)






Equation





1








Δ





Fw

=


-


f
0


π






z
q







G





G


2




ω
2



ρ
2



η
2



h
1













Viscoelastic





term






(
2
)






Equation





2








Δ





Fs

=



-


f
0


π






z
q






ωρ
1



h
1


+



f
0


π






z
q






G





G


2




ω
2



ρ
2



η
2



h
1














Mass





load





Viscoelastic





term






(
3
)






Equation





3







Viscoelastic moduli G′ (storage elastic modulus) and G″ (loss elastic modulus) of the film are calculated by the following equation 4 and a relational expression G″=C*G′,










G


=



ω
·

[

mass





load

]



[

viscoelastic





term






(
2
)


]


·



ρ
2



η
2




(

1
+

C
2


)

·

ρ
1








Equation





4







wherein ΔFn=Δ(F1−F2)/2, G is a complex elastic modulus (MPa), G′ is a storage elastic modulus (dynamic elastic modulus) (MPa), G″ is a loss elastic modulus (dynamic loss) (MPa), ω is an angular frequency, ρ2 is the density (g/cm3) of the solution, η2 is the viscosity (Pa s) of the solution, h1 is the thickness (nm) of the formed film, ρ1 is the density (g/cm3) of the formed film, f3 is a fundamental frequency (Hz), and Zq is the shear mode acoustic impedance (gm/sec/cm2) of crystal.


A second solution means in the present invention is characterized in that the viscoelasticity loss coefficient tanδ of the adsorbed substance, the rigidity coefficient μ of the adsorption substance, and the viscosity coefficient η of the adsorbed substance are calculated using the viscoelastic modulus G′ in the first solution means.


A third solution means in the present invention is characterized in that the viscoelastic modulus, the loss coefficient, the rigidity coefficient, or the viscosity coefficient of the adsorption substance are calculated along with the measurement of the frequency variation in the first solution means.


A fourth solution means in the present invention is characterized in that the piezoelectric element is a quartz crystal oscillator, an APM (acoustic plate mode sensor), an FPW (flexural plate-wave sensor), or a SAW (sourface acoustic-wave sensor) in the first solution means.


A measuring apparatus in the present invention, which is an apparatus for measuring the viscoelastic modulus of the substance in a system for forming a film by adsorbing a substance to the surface of a piezoelectric element or to a film fixed onto the piezoelectric element in a solution, using a sensor with the piezoelectric element submerged on both sides or one side in the solution,


includes a measuring means for measuring at least two of a resonance frequency Fs, and half-value frequencies F1 and F2 (F2>F1) having half conductance values of a conductance value of the resonance frequency in the piezoelectric element, and


includes an arithmetic means for using the values measured by the measuring means to calculate amass load term, a viscoelastic term (1), a viscoelastic term (2), and a viscoelastic term (3) represented by the following equations 1 to 3:











Δ






F
2


=



-


f
0


π






z
q





ω






ρ
1



h
1


+



f
0


π






z
q






(


G


+

G



)




G


2




ω
2



ρ
2



η
2



h
1














Mass





load





Viscoelastic





term






(
1
)






Equation





1








Δ





Fw

=


-


f
0


π






z
q







G





G


2




ω
2



ρ
2



η
2



h
1













Viscoelastic





term






(
2
)






Equation





2








Δ





Fs

=



-


f
0


π






z
q






ωρ
1



h
1


+



f
0


π






z
q






G





G


2




ω
2



ρ
2



η
2



h
1














Mass





load





Viscoelastic





term






(
3
)






Equation





3







and further to calculate viscoelastic moduli G′ (storage elastic modulus) and G″ (loss elastic modulus) of the film by the following equation 4 and a relational expression G″=C*G′,










G


=



ω
·

[

mass





load

]



[

viscoelastic





term






(
2
)


]


·



ρ
2



η
2




(

1
+

C
2


)

·

ρ
1








Equation





4







wherein ΔFs=Δ(F1−F2)/2, G is a complex elastic modulus (MPa), G′ is a storage elastic modulus (dynamic elastic modulus) (MPa), G″ is a loss elastic modulus (dynamic loss) (MPa), ω is an angular frequency, ρ2 is the density (g/cm3) of the solution, η2 is the viscosity (Pa s) of the solution, h1 is the thickness (nm) of the formed film, ρ1 is the density (g/cm3) of the formed film, f0 is a fundamental frequency (Hz), and Zq is the shear mode acoustic impedance (gm/sec/cm2) of crystal.


Effects of the Invention

Accurate information on viscoelasticity has not been conventionally obtained because of attaining the measurement of only the resonance frequency Fs including an element of a mass load due to adsorption and an element of viscoelasticity of an adsorption substance. Further, although elements of mass, viscosity, and viscoelasticity could be separately obtained by the invention proposed in Patent Document 1, they could be only represented as frequency variations.


However, the present invention attained a method for calculating general G′ and G″ values as parameters representing viscoelasticity from the obtained frequency variations of the respective elements. Consequently, viscoelastic moduli G′ and G″ can be obtained in real time during the measurement of the frequency variations, and physical information (a loss coefficient, an elasticity coefficient, and a viscosity coefficient) other than the frequencies can be acquired in an instant. A more accurate physical property evaluation of the adsorption substance can thereby be made.





BRIEF DESCRIPTION OF DRAWINGS


FIG. 1 is an explanatory diagram of apparatus configuration of the embodiment of the present invention.



FIG. 2 is a graph showing the measured result of one example of the present invention.



FIG. 3 is a graph showing changes in the value calculated based on the measured result of the same example.



FIG. 4 is a graph showing changes in the viscoelastic coefficient and the viscosity coefficient in the same example.



FIG. 5 is an explanatory diagram of the measured frequencies of N-th waves etc. of one embodiment in the present invention.



FIG. 6 is a graph showing frequency change values for use in other examples of the present invention.



FIG. 7 is a graph showing the result of an analysis made using a measuring method of the present invention using the frequency change values in FIG. 6.



FIG. 8 is a graph showing the result of an analysis made using a measuring method of the present invention using the frequency change values in FIG. 6.





DESCRIPTION OF REFERENCE NUMERAL




  • 1 Cell


  • 2 Measuring means (network analyzer)


  • 3 Control means


  • 4 Temperature control means


  • 5 Temperature regulating means


  • 6 First arithmetic part


  • 7 Second arithmetic part


  • 8 Display part



MODE FOR CARRYING OUT THE INVENTION

The measurement principle of the present invention is explained as follows.


The change of impedance Z in the case of adsorbing a viscoelastic film to a quartz crystal oscillator in a solution by a transmission theory of Martin et al. (V. E. Granstaff, S. J. Martin, J. Appl. Phys. 1994, 75, 1319) is expressed by the following equation (1).









Equation





5











Z
=




(


ωρ
2




η
2

/
2


)


1
/
2




(

1
+
j

)



jωρ
1



h
1


+



(


G


-

j






G




)




G


2




ω
2



ρ
2



η
2



h
1







(
1
)







From the equation (1), the change of a resonance frequency Fs is expressed by the equation (2), and the variation of half (F1−F2)/2 (=Fw) of a half-value frequency is expressed by the equation (3).














Equation











6













Δ





Fs

=


Im


(
Z
)


=



-


f
0


π






z
q







(

ω






ρ
2




η
2

/
2


)


1
/
2



-



f
0


π






z
q





ωρ
1



h
1


+



f
0


π






z
q






G





G


2




ω
2



ρ
2



η
2



h
1








(
2
)











Equation





7













Δ



(


F
1

-

F
2


)

2


=


Δ





Fw

=


Re


(
Z
)


=



-


f
0


π






z
q







(


ωρ
2




η
2

/
2


)


1
/
2



-



f
0


π






z
q






G





G


2




ω
2



ρ
2



η
2



h
1









(
3
)







From the equations (2) and (3), the variation of a frequency F2 can be obtained as shown in the equation (4).









Equation





8












Δ






F
2


=



-


f
0


π






z
q






ωρ
1



h
1


+



f
0


π






z
q






(


G


-

G



)




G


2




ω
2



ρ
2



η
2



h
1







(
4
)







In this specification, G is a complex elastic modulus (MPa), G′ is a storage elastic modulus (dynamic elastic modulus) (MPa), G″ is a loss elastic modulus (dynamic loss) (MPa), ω is an angular frequency, ρ2 is the density (g/cm3) of the solution, η2 is the viscosity (Pa s) of the solution, h1 is the thickness (nm) of the formed film, ρ1 is the density (g/cm3) of the formed film, η1 is a viscosity coefficient (Pa s) of the formed film, f0 is a fundamental frequency (Hz), and Zq is the shear mode acoustic impedance (gm/sec/cm2) of crystal.


A Voight model frequently used as a viscoelasticity model of a film is applied to G′ and G″.


A model in which a spring G of an elastic element and a dashpot η are connected in parallel is expressed by the following expression:






G=G′+jG″=μ+jωη


Herein, in the case of ωη=Cμ (C: constant, μ: rigidity coefficient (MPa)), the equation (2) can be transformed as follows:













Δ





Fw

=


-


f
0


π






z
q







G





G


2




ω
2



ρ
2



η
2



h
1








=


-


f
0


π






z
q






μ

(


μ
1
2

+


ω
2



η
1
2



)




ω
2



ρ
2



η
2



h
1









Equation





9







In the case of a N-th wave, the equation is as follows:










Δ






Fw
N


=


-


f
0


π






z
q






1

(

1
+


C
2



N
2



)






N
2



ω
2



ρ
2



η
2



h
1



μ
1







Equation





10







An equation including a constant C is obtained by making division of respective ΔFw values of two frequencies.











Δ






Fw
3



Δ






Fw
1



=

A
=




-

1

(

1
+

9


C
2



)






9


ω
2



ρ
2



η
2



h
1



μ
1





-

1

(

1
+

C
2


)







ω
2



ρ
2



η
2



h
1



μ
1




=


9


(

1
+

C
2


)



1
+

9


C
2










Equation





11







The above equation is an equation in the case of using the frequency variation Fw of a fundamental wave (N=1) and a threefold wave (N=3).


A value obtained by dividing the frequency variation Fw3 obtained by measurement, by Fw1 is made A, and the value A is substituted in the left side of the above equation to calculate the following constant C.









C
=


1
3

*



9
-
A


A
-
1








Equation





12







The frequency changes of the respective terms indicated in the following equations 1 to 3 are obtained from the calculated value C and the measured frequency variations ΔF2, ΔFw, and ΔFs values of the fundamental wave. This measurement is made on condition that there is no occurrence of the viscous load of the solution.


The change of the viscous load occurs when a buffer solution and a sample solution greatly differ in viscosity, for example, when a sample containing a glycerol solution of high concentration is added to the buffer solution for measurement. Inmost cases, the sample solution and the buffer solution are almost equal in viscosity, and it is therefore considered that no viscous change occurs.











Δ






F
2


=



-


f
0


π






z
q





ω






ρ
1



h
1


+



f
0


π






z
q






(


G


+

G



)




G


2




ω
2



ρ
2



η
2



h
1














Mass





load





Viscoelastic





term






(
1
)






Equation





1








Δ





Fw

=


-


f
0


π






z
q







G





G


2




ω
2



ρ
2



η
2



h
1













Viscoelastic





term






(
2
)






Equation





2








Δ





Fs

=



-


f
0


π






z
q






ωρ
1



h
1


+



f
0


π






z
q






G





G


2




ω
2



ρ
2



η
2



h
1














Mass





load





Viscoelastic





term






(
3
)






Equation





3







First, since ΔFw is established only with the viscoelastic term (2),





Viscoelastic term (2)=ΔFw


is satisfied.


Secondly, the viscoelastic term (1) is obtained by the following expression. Since the constant C can be expressed by G″/G′, C=G″/G′ is substituted and multiplied by the viscoelastic term (2) to calculate the viscoelastic term (1).





Viscoelastic term (1)=−(1+C)*Viscoelastic term (2)


The mass load term is obtained by substituting the obtained value of the viscoelastic term (1) in an expression of ΔF2.





Mass load term=ΔF2−Viscoelastic term (1)


Finally, the viscoelastic term (3) is obtained by substituting the obtained value of the mass load term in an expression of ΔFs.





Viscoelastic term (3)=ΔFs−Mass load term


On this principle, it is proved that the respective terms can be separated as the frequency variations.


Further, the obtained frequency variations of the mass load term and the viscoelastic term (2) are used to calculate the viscoelastic moduli G′ and G″ by the following equation, wherein ω is an angular frequency, ρ2 is the density (g/cm3) of a solution, η2 is the viscosity (Pa s) of a solution, and pi is the density (g/cm3) of a film.











G


=



ω
·

[

mass





load

]



[

viscoelastic





term






(
2
)


]


·



ρ
2



η
2




(

1
+

C
2


)

·

ρ
1












G


=

C
*

G








Equation





4







The constant C also corresponds to a loss coefficient tans (=G″/G′) representing a loss portion of a viscoelastic body. The loss coefficient is one of the parameters indicating the hardness of the viscoelastic body by indicating an ideal solid when the loss coefficient is 0 and indicating an ideal liquid when the loss coefficient is ∞.


Further, since the Voight model G=G′+jG″=μ+jωη is applied to the obtained viscoelastic modulus, the rigidity coefficient μ (MPa) and viscosity coefficient η (Pa s) of the adsorption substance can be calculated from the G′ and G″ values.


The frequencies Fs, F1, and F2 used in the above measuring method can be measured by a method based on an oscillating circuit, a method obtained by frequency scanning from external equipment such as an impedance analyzer or a network analyzer, etc., and the measuring method is not limited as long as the resonance frequency Fs, and the half-value frequencies F1 and F2 (F2>F1) having half conductance values of a conductance value of the resonance frequency are measured.


As shown in FIG. 1, the apparatus of the present invention includes a control means 3 controlling a cell 1 having a piezoelectric element, through a network analyzer 2 which is a measuring means; a first arithmetic part 6 for computing the frequency components of frequencies measured by the network analyzer 2; a second arithmetic part 7 for computing film thickness and viscoelasticity using the result obtained in the first arithmetic part 6; and a display means 8 for alternately or simultaneously displaying the values of the frequencies, viscoelasticity, etc. obtained in the respective arithmetic parts. The control means 3, the arithmetic parts 6, 7, and the display part 8 can be composed of usual personal computers and monitors.


In the example shown in the figure, a temperature control means 4 such as a Peltier element is provided at the lower face of the cell 1 to regulate the temperature of the cell 1, and a temperature regulating means 5 for regulating the temperature control means 4 is controlled in the same way by the control means 3.


The fundamental wave and the threefold wave are used in the above explanation. However, the present invention can be used if at least two frequencies are used out of overtone frequencies including the fundamental wave as shown in FIG. 5.


The piezoelectric element used in the present invention is not limited as long as it can measure the target frequencies, and a quartz crystal oscillator, an APM (acoustic plate mode sensor), an FPW (flexural plate-wave sensor), or a SAW (sourface acoustic-wave sensor) can be also used.


EXAMPLES

The obtained result of viscoelastic moduli of a biomolecule film and the like is shown in the following table 1 using the apparatus explained referring to the FIG. 1, and using the fundamental wave and the threefold wave of 27 MHz according to the present invention.


Each sample used in the experiment is as follows.


1) 50 nm beads (the final concentration when measured in latex beads is 10 μg/ml)


2) Neutravidin (the final concentration when measured in globular protein is 10 μg/ml)


3) ssDNA 30 mer (single-stranded DNA with a length of 30 mer, and the final concentration when measured is 1 μg/ml)


4) ssDNA 60 mer (single-stranded DNA with a length of 60 mer, and the final concentration when measured is 1 μg/ml)


5) ssDNA 90 mer (single-stranded DNA with a length of 90 mer, and the final concentration when measured is 1 μg/ml)


6) dsDNA 90 mer (double-stranded DNA with a length of 90 mer, and the final concentration when measured is 1 μg/ml)












TABLE 1







G′ (Mpa)
G″ (Mpa)




















50 nm beads
2.090
1.264



Neutravidin
1.660
1.066



ssDNA 30 mer
0.360
0.315



ssDNA 60 mer
0.340
0.237



ssDNA 90 mer
0.301
0.201



dsDNA 90 mer
1.790
0.443










Next, the measured data of the fundamental wave (N1) and threefold wave (N3) of the ssDNA 90 mer are specifically shown in FIGS. 2 to 4.


The FIG. 2 is a graph showing the measured data of respective ΔF2 and ΔFw of the above two kinds of different frequencies, and a graph in FIG. 3 shows viscoelastic moduli G′ and G″ of a DNA film calculated in real time from the measured data.


A graph in FIG. 4 shows the result of an elasticity coefficient and a viscosity coefficient of the DNA film obtained in real time in the same way from the viscoelastic moduli obtained in FIG. 3.


From the above result, the viscoelastic moduli G′ and G″ could be obtained from the frequency variations by the method in the present example, and it was proved that the loss coefficient, the elasticity coefficient, and the viscosity coefficient of the object to be measured can be measured in real time at the time of measurement.


Other examples using the measuring method in the present invention are explained.



FIG. 6 shows the frequency change values obtained by measurement, and FIGS. 7 and 8 show the result of making an analysis using the measuring method of the present invention with the use of the frequency change values.


In each figure, a to c are as follows:


a. a region for modifying and stabilizing NeutrAvidin (protein) on a gold electrode in a solution


b. a change region when injecting 60 mer biotinylated DNA


c. a change region when injecting 60 mer DNA which is a complementary strand of the above biotinylated DNA to cause hybridization



FIGS. 6 to 8 proved the occurrence of an adsorption phenomenon from changes in ΔFs and ΔF2 in the region b or c in FIG. 1 in the conventional QCM measuring method, but the viscoelasticity of the adsorption film cannot be digitized. However, physical property values such as viscoelastic moduli, mass, and film thickness of the film formed on the electrode in the solution can be calculated by the method of the present invention.


The viscoelastic change of a ΔG′ value occurring from b to c in FIG. 7 shows an increase in hardness in the DNA film due to the occurrence of hybridization of the single-stranded DNA with the complementary strand. The tans called the loss coefficient shows a value 0 as a complete solid and a value ∞ as a complete solution. However, the film hardness was proved to be increased from the decrease of the tanδ.


Further, FIG. 8 proved that the adsorption of the single-stranded DNA and the adsorption of the complementary strand can be measured by the increase of mass.

Claims
  • 1. A method for measuring a viscoelastic modulus of a substance in a system for forming a film by adsorbing a substance to the surface of a piezoelectric element or to a film fixed onto the piezoelectric element in a solution, using a sensor with the piezoelectric element submerged on both sides or one side in the solution, the method comprising: calculating a mass load term, a viscoelastic term (1), a viscoelastic term (2), and a viscoelastic term (3) represented by the following equations 1 to 3:
  • 2. The method for measuring a viscoelastic modulus of a substance according to claim 1, wherein the viscoelasticity loss coefficient tanδ of the adsorbed substance, the rigidity coefficient μ of the adsorption substance, and the viscosity coefficient η of the adsorbed substance are calculated using the viscoelastic modulus G′.
  • 3. The method for measuring a viscoelastic modulus of a substance according to claim 1, wherein the viscoelastic modulus, the loss coefficient, the rigidity coefficient, or the viscosity coefficient of the adsorption substance are calculated along with the measurement of the frequency variation.
  • 4. The method for measuring a viscoelastic modulus of a substance according to claim 1, wherein the piezoelectric element is a quartz crystal oscillator, an APM (acoustic plate mode sensor), an FPW (flexural plate-wave sensor), or a SAW (sourface acoustic-wave sensor).
  • 5. An apparatus for measuring a viscoelastic modulus of a substance in a system for forming a film by adsorbing the substance to the surface of a piezoelectric element or to a film fixed onto the piezoelectric element in a solution, using a sensor with the piezoelectric element submerged on both sides or one side in the solution, the apparatus comprising: a measuring means for measuring at least two of a resonance frequency Fs, and half-value frequencies F1 and F2 (F2>F1) having half conductance values of a conductance value of the resonance frequency in the piezoelectric element; andan arithmetic means for using the values measured by the measuring means to calculate amass load term, a viscoelastic term (1), a viscoelastic term (2), and a viscoelastic term (3) represented by the following equations 1 to 3:
Priority Claims (1)
Number Date Country Kind
2010-235657 Oct 2010 JP national
PCT Information
Filing Document Filing Date Country Kind 371c Date
PCT/JP11/05810 10/18/2011 WO 00 4/16/2013