Relation between Mass Sensitivity and Complex Power Flow in Love Wave Sensors

In this paper, we investigate the connection between average power flows in Love wave waveguides with the mass sensitivity of Love wave sensors. In fact, loading with a Newtonian liquid gives rise to two extra power flows, in the transverse direction towards the loading Newtonian liquid. The first is an active power flow feeding viscous losses in the Newtonian liquid and the second is a reactive power flow that is responsible for the phase delay of the Love wave and consequently for the changes in phase velocity of the Love wave. Since loading with a lossless mass also leads to changes in the phase velocity, we assert that mass sensitivity Sσvp of Love wave sensors is connected to the average reactive power flow, in the transverse direction x2, bouncing back and forth, between the interior of the waveguide and the loading Newtonian liquid. Subsequently, we found the thickness of the effective surface layer of mass that is equivalent to loading with a semi-infinite Newtonian liquid. The analytical formulas developed in this paper are illustrated by the results of numerical calculations performed for an exemplary Love wave waveguide composed of a PMMA surface layer deposited on an ST-Quartz substrate.


Introduction
Bulk and surface elastic waves [1][2][3][4] are widely used in sensors of physical properties of materials as well as in biosensors and chemosensors [5][6][7][8][9][10][11][12][13][14]. The elastic surface waves of the Love type are especially attractive for use in sensors working in a liquid environment, since Love surface waves can propagate long distances when loaded by low viscosity liquids, such as water. Since most chemosensors and biosensors [15][16][17][18] are used in a water-like environment, the Love wave sensors offer a significant advantage over other types of elastic sensors based on bulk elastic waves or Rayleigh-type elastic surface waves.
In this paper, we will try to connect the average power flow occurring in Love wave waveguides loaded with a Newtonian liquid with some engineering parameters of Love wave sensors, such as their coefficient of mass sensitivity S v p σ . At first glance, the connection between the average power flow and the coefficient S v p σ is not obvious. However, the results obtained in this paper suggest that in fact the average power flow and the mass sensitivity of the sensor are intimately related, i.e., they display the same qualitative dependencies (maxima) as a function of wave frequency f and thickness h 1 of the surface layer of the waveguide.
In Love wave waveguides, loaded with a Newtonian liquid, the average (total) power flows P 1 = ∞ −∞ P 1 (x 2 )dx 2 and P 2 = ∞ −∞ P 2 (x 2 )dx 2 in both directions, i.e., in the direction of propagation x 1 and in the transverse direction x 2 , can be of the active or reactive type. The corresponding Poynting vectors P 1 (x 2 ) and P 2 (x 2 ), were defined in (Auld, 1990). In this paper we will focus on the average power flow P 2 in the transverse direction x 2 as well as on the corresponding complex Poynting vector P 2 (x 2 ) evaluated at the interface x 2 = 0, between the loading Newtonian liquid and the elastic surface layer of the waveguide. Indeed, we will argue that the reactive part ImP 2 of the average power flow in the transverse In other words, the active power flow ReP 2 in the transverse directions x 2 feeds viscous losses in the loading Newtonian liquid and is ultimately dissipated to heat. On the other hand, the reactive power flow ImP 2 in the transverse directions x 2 gives rise to delay of the Love wave and therefore contributes to changes in its phase velocity.
The key parameter of the Love wave sensor is its coefficient of mass sensitivity S v p σ , i.e., the relative change ∆v p /v p in the phase velocity v p due to the loading of the waveguide with the surface mass density σ kg/m 2 [5]. The mass sensitivity S v p σ depends on the surface layer thickness h 1 , frequency f , and other material parameters of the Love wave waveguide [19,20].
Despite the fact that the coefficient of mass sensitivity S v p σ was initially defined for waveguides loaded with an infinitesimally thin layer of lossless mass, in this paper we extend the notion of the mass sensitivity S v p σ on waveguides loaded with a semi-infinite lossy Newtonian liquid of a density ρ 0 and viscosity η 0 . To this end, we introduce the notion of the effective layer of thickness h = δ 0 /2 that provides the same phase delay as a semi-infinite Newtonian liquid loading the waveguide, where δ 0 is the penetration depth of the Love wave into the Newtonian liquid. As a result, the equivalent surface mass density σ of a semi-infinite Newtonian liquid loading Love wave waveguide equals σ = ρ 0 (δ 0 /2) [kg/m 2 ], see Section 4.2.
The analytical formulas for the mass sensitivity S v p σ and the reactive power flow ImP 2 in the transverse direction x 2 , were supported by the numerical calculations performed for a Love wave waveguide composed of a PMMA-Poly(methyl methacrylate) surface layer deposited on ST-Quartz substrate. In fact, Quartz is very attractive as a material for substrates in Love wave sensors, since it is the only common piezoelectric material that supports pure shear SH bulk waves that can be generated and received via a piezoelectric effect, using interdigital transducers in a wide range of frequencies, namely from ∼1 MHz to ∼1 GHz . On the other hand, PMMA is a very good candidate as a material for surface layers in Love wave waveguides, due to its low phase velocity (∼1100 m/s) of bulk SH waves that promotes strong energy trapping in the surface layer itself.
To the best of our knowledge, the relationship between the mass sensitivity S v p σ and complex power flow P 2 in Love wave sensors was not yet analyzed in the literature and therefore can be considered as an original contribution to the state-of-the art.
The results of the theoretical analysis and numerical calculations presented in this paper provide a deeper insight into the physical phenomena occurring in Love wave sensors and can serve as a basis for better design and optimization of Love wave sensors, biosensors, and chemosensors, working in a liquid environment.

Physical Model
In this paper, we investigate Love wave waveguides with a single elastic surface layer rigidly bonded to a semi-infinite elastic substrate. The top surface of the waveguide is loaded with a Newtonian liquid of a semi-infinite extent (see Figure 1).
The waveguide structure shown in Figure 1 represents the main components of the simplest Love wave sensor working in a liquid environment.
Part of the energy of the Love wave penetrates into the Newtonian liquid, where it is dissipated into heat. Therefore, wavenumber k of the Love wave is a complex quantity given by the following formula: where k p = ω/v p determines the phase velocity v p of the Love wave, α is the attenuation, ω stands for the angular frequency and "j" is the imaginary unit. Part of the energy of the Love wave penetrates into the Newtonian liquid, where it is dissipated into heat. Therefore, wavenumber of the Love wave is a complex quantity given by the following formula: where = / determines the phase velocity of the Love wave, is the attenuation, stands for the angular frequency and " " is the imaginary unit.

Mathematical Model
Love surface waves, analyzed in this paper, are time-harmonic, propagate in the direction , and are uniform along the transverse direction . Therefore, their mechanical displacement ( ) and shear stresses ( ) in the Newtonian liquid ( = 0), elastic surface layer ( = 1) and in the elastic substrate ( = 2), see Figure 1, will be sought in the following generic form: where ( ) ( ), ( ) , and ( ) express variations of the mechanical displacement and shear stress in the transverse direction , is the wavenumber of the Love surface wave, its angular frequency and index = 0, 1 and 2. By definition, shear stresses ( ) , ( ) of the Love wave are connected to the mechanical displacement ( ) by the following formulas: where ( ) is the shear modulus of elasticity of the constituting medium number .

Mathematical Model
Love surface waves, analyzed in this paper, are time-harmonic, propagate in the direction x 1 , and are uniform along the transverse direction x 3 . Therefore, their mechanical displacement u (i) 3 and shear stresses τ (i) 23 in the Newtonian liquid (i = 0), elastic surface layer (i = 1) and in the elastic substrate (i = 2), see Figure 1, will be sought in the following generic form: where u 44 is the shear modulus of elasticity of the constituting medium number i.

Differential Equations
Mechanical displacement u (i) 3 of Love surface waves is governed by the following wave equation [21] resulting from the second Newton's law of motion: where v i = c (i) is the phase velocity of SH bulk waves in medium number i, c (i) 44 is its shear modulus of elasticity and ρ i is the density.
Substituting Equation (2) for the mechanical displacement u (i) 3 into wave equation Equation (4) one obtains the following ordinary differential equation of the Helmholtz type: where k i = ω/v i is the wavenumber of bulk SH waves in medium number i = 0, 1 and 2.

Mechanical Displacement u (i)
3 (x 2 ) and Shear Stresses τ Since the mechanical displacement u 3 (x 2 ) of the Love wave must vanish for x 2 → −∞ , we will seek the solution to Equation (5) in the following form: where q 0 = k 2 − k 2 0 1/2 is the complex transverse wave number of the Love wave in the Newtonian liquid and C 0 is an arbitrary constant. In order to fulfill the condition 3 (x 2 ) → 0 for x 2 → −∞ the real part of the transverse wavenumber q 0 must be positive, Req 0 > 0.
By virtue of Equation (3a), shear stress τ (0) 23 of the Love wave in the Newtonian liquid is given by: where c  Similarly, using Equation (3b), shear stress τ (0) 13 of the Love wave in the Newtonian liquid is given by: where c

Elastic Surface Layer
Since the elastic surface layer is of a finite thickness h 1 the solution to Equation (5) can be sought in the following form: where C 1 and C 2 are arbitrary constants and the transverse wave number of the Love wave in the elastic surface layer q 1 = k 2 1 − k 2 1/2 . Since the elastic surface layer is of a finite thickness h 1 , the real part of q 1 can be either positive or negative. By virtue of Equation (3a), shear stress τ (1) 23 of the Love wave in the elastic surface layer is given by: The shear stress τ (1) 23 will enter into the appropriate boundary conditions at two interfaces x 2 = 0 and x 2 = h 1 of the waveguide as well as into equations for the Poynting vector P (1) 2 (x 2 ) of the Love wave in the transverse direction x 2 . Analogously, employing Equation (3b), shear stress τ (1) 13 of the Love wave in the elastic surface layer is given by: The shear stress τ (1) 13 . will enter into equations for the Poynting vector P 1 (x 2 ) of the Love wave in the direction of propagation x 1 .

Elastic Substrate
Since the amplitude of the Love surface wave in the elastic substrate must tend to zero for x 2 → ∞ , as a solution to Equation (5) we choose the following expression: where the transverse wave number of the Love wave in the substrate 44 is its modulus of elasticity and ρ 2 is its density. By virtue of Equation (3a), shear stress τ (2) 23 of the Love wave in the elastic substrate is given by: The shear stress τ (2) 23 will enter into the appropriate boundary conditions at the interface x 2 = h 1 of the waveguide as well as into equations for the Poynting vector P On the other hand, using Equation (3b), shear stress τ (2) 13 of the Love wave in the elastic substrate is given by: The shear stress τ (2) 13 will enter into equations for the Poynting vector P 1 (x 2 ) of the Love wave in the direction of propagation x 1 .

Boundary Conditions
Boundary conditions at two interfaces x 2 = 0 and x 2 = h 1 of the waveguide shown in Figure 1 require the continuity of the mechanical displacement u   [1]. Consequently, at the interface x 2 = 0 between the Newtonian liquid, and the elastic surface layer we can write: Similarly, at the interface x 2 = h 1 between the elastic surface layer and the elastic substrate we have: u Substituting Equations (6), (7), (9), (10), (12) and (13) into boundary conditions Equations (15)- (18), one obtains the following set of homogeneous linear algebraic equations for the unknown coefficients C 0 , C 1 , C 2 and C 3 : Since the set of 4 linear algebraic equations (Equation (19)) is homogeneous we can determine only 3 independent coefficients C i in function of the remaining one, say C 0 . Consequently, the coefficients C 1 , C 2 , C 3 can be expressed in terms of the coefficient C 0 as:

Dispersion Equation
The set of homogeneous linear algebraic equations given by Formula (19) in a matrix form has a non-trivial solution if the determinant of its left-hand matrix equals zero. This condition leads to the following dispersion relation for the phase velocity and attenuation of the Love wave propagating in the investigated waveguide shown in Figure 1: Formula (21) is a complex dispersion equation of the Love wave, propagating in the waveguide structure shown in Figure 1. Equation (21) can be split into its real and imaginary parts, providing therefore a set of two nonlinear algebraic equations for unknown k p and α [21]. The resulting set of two nonlinear transcendental algebraic equations can be solved numerically, using for instance a two-dimensional Newton-Raphson procedure.
Indeed, solving numerically the dispersion equation is a prerequisite in the analysis of Love surface waves, propagating in any waveguide structure, since without the knowledge of the complex wavenumber k = k p + jα, we cannot evaluate numerical values of other field quantities of the Love wave, such as the mechanical displacement u 3 , shear stress τ 23 or Poynting vector P 2 .

Complex Poynting Vector
The power flux in the investigated Love wave waveguide is represented by the complex Poynting vector with two components, i.e., P 1 (x 2 ) in the direction of propagation x 1 and P 2 (x 2 ) in the transverse direction x 2 .
In this section, we will develop analytical formulas for the complex Poynting vector P 2 (x 2 ) in the transverse direction x 2 of the waveguide that is loaded with a Newtonian liquid at its surface (x 2 = 0).
In free waveguides not loaded with a Newtonian liquid, the power flux P 2 (x 2 ) across the free interface x 2 = 0 is obviously zero, P 2 (x 2 = 0) = 0, since no power can be transmitted from the waveguide into vacuum.
This situation changes drastically in waveguides loaded with a Newtonian liquid, since the acoustic power flux P 2 (x 2 ) across the interface x 2 = 0 is now clearly non-zero, P 2 (x 2 = 0) = 0. As a result, an additional complex power flow, between the interior of the waveguide and the bulk of the loading Newtonian liquid, occurs.
Thus, the complex Poynting vector P 2 (x 2 ) = ReP 2 (x 2 ), evaluated at the interface x 2 = 0, contains unique information about viscous η 0 and inertial ρ 0 properties of the loading Newtonian liquid.
By definition, the complex Poynting vector P 2 (x 2 ) W/m 2 in the transverse direction x 2 can be expressed as [2]: where the asterisk " * " stands for complex conjugation and "j" is the imaginary unit. Due to the interaction with the viscous Newtonian liquid, Love surface waves propagating in the direction x 1 in the waveguide structure shown in Figure 1, will be gradually attenuated, what can be expressed analytically as P 2 (x 1 , Substituting Equations (6) and (7), for the mechanical displacement u   (20) for the coefficients C 1 , C 2 , C 3 into Equation (22) one obtains the following expression for the complex Poynting vector P where Re(q 0 ) > 0.

Elastic Surface Layer
Substituting Equations (9) and (10), for the mechanical displacement u 3 (x 2 ) and shear stress τ (1) 23 (x 2 ) of the Love wave in the elastic surface layer, as well as Formula (20) for the coefficients C 1 , C 2 , C 3 into Equation (22), one obtains the following expression for the complex Poynting vectors P (1) 2 (x 2 ) in the elastic surface, along the transverse direction x 2 : where the auxiliary functions Substituting Equations (12) and (13), for the mechanical displacement u 3 (x 2 ) and shear stress τ (2) 23 (x 2 ) of the Love wave in the elastic substrate, as well as Formula (20) for the coefficients C 1 , C 2 , C 3 into Equation (22), one obtains the following expression for the complex Poynting vectors P 2 (x 2 ) in the elastic substrate, along the transverse direction x 2 : where Re(q 2 ) > 0 and h 1 stands for thickness of the elastic surface layer. As it was stated before in this section, the complex Poynting vector P 2 (x 2 ) evaluated at the interface x 2 = 0 contains unique information about viscous η 0 and inertial ρ 0 properties of the loading Newtonian liquid.
According to Equation (23) the complex Poynting vector P (0) 2 (x 2 ) evaluated at the interface x 2 = 0 is given by: The transverse wavenumber q 0 of the Love wave in the Newtonian liquid is a complex quantity; therefore, it can be represented as q 0 = a 0 + jb 0 . On the other hand, the elastic modulus of the Newtonian liquid c (0) 44 = −jωη 0 . Therefore, Equation (26) can be written as: Since for the Newtonian liquid q 0 = k 2 − k 2 0 1/2 and k 2 0 = ω 2 ρ 0 /(−jωη 0 ) as well as the wavenumber of the Love wave k = k p + jα, we can write the following equation: Thus, real a 0 and imaginary b 0 part of the transverse wavenumber q 0 of the Love wave in the Newtonian liquid equal: where k p = ω/v p Consequently, from Equation (27) it follows that: where a 0 , b 0 are given by Equation (29). Equation (30) shows that the real and imaginary parts of the complex Poynting vector 2 (x 2 ), evaluated at the interface x 2 = 0, are involved functions of v p , α, ω, η 0 and ρ 0 . It should be noticed that v p and α depend in turn on the material and geometrical parameters of the waveguide and Newtonian liquid.  By definition, the complex Poynting vector P 1 (x 2 ) W/m 2 in the direction of propagation x 1 is given by [2]:

Since in general ReP
where the asterisk " * " stands for complex conjugation and "j" is the imaginary unit.
The integrals in Equation (33) can be evaluated analytically. However, due to their excessive length and complexity they will be not reproduced here. Since Poynting vectors 2 (x 2 ) and P (2) 2 (x 2 ) are complex, the average power flow P 2 in the transverse direction x 2 is a complex-valued quantity, therefore it can be written as P 2 = ReP 2 + jImP 2 , where in general ReP 2 = 0 and ImP 2 = 0.
The reactive part ImP 2 of the average power flow P 2 , in the transverse direction x 2 , in the Newtonian liquid P 2 (x 2 )dx 2 , will be plotted in Section 5.1, as a function of frequency f and thickness of the elastic surface layer h 1 .
It should be noticed that the reactive parts of the average power flow in the elastic surface layer P 2 (x 2 )dx 2 and in the elastic substrate P 2 (x 2 )dx 2 are always non-zero in free waveguides, not loaded with a Newtonian liquid, since they are responsible for the energy storage occurring in Love wave waveguides in the transverse direction x 2 . However, the active part of the average power flow in the elastic surface layer P     The average (total) power flow P 1 in the direction of propagation x 1 is defined as: where P (0) 1 (x 2 ) and P 1 (x 2 ) are, respectively, complex Poynting vectors in the Newtonian liquid, elastic surface layer and in the elastic substrate, given by Equation (32).
The integrals in Equation (34) can be evaluated analytically. However, due to their excessive length and complexity they will be not reproduced here. Since Poynting vectors 1 (x 2 ) and P (2) 1 (x 2 ) are complex, the average power flow P 1 in the direction of propagation x 1 is a complex-valued quantity as well and therefore can be written as P 1 = ReP 1 + jImP 1 , where in general ReP 1 = 0 and ImP 1 = 0. The active part ReP 1 of the average power flow P 1 in the direction of propagation x 1 will serve as a normalization factor in Section 5.1.

Coefficient of Mass Sensitivity S v p σ
One of the most important parameters characterizing the quality of Love wave sensors is the coefficient of mass sensitivity. Until now, this coefficient of mass sensitivity was evaluated using the perturbation theory and/or numerical methods such as the finite element method, etc. [21][22][23][24][25][26]. In this paper, to evaluate the mass sensitivity, we used analytical methods that provide a deeper insight into the physical background of the operation of Love wave sensors [19,20].
The sensitivity of Love wave sensors to mass loading was initially defined in lossless waveguides, where a semi-infinite Newtonian liquid in Figure 1 is replaced by an infinitesimal layer of lossless mass that only changes the phase velocity v p of the Love wave, without introducing any extra attenuation, see Figure 2.

Average (Total) Power Flow = ( ) in the Direction of Propagation
The average (total) power flow in the direction of propagation is defined as: where ( ) ( ), ( ) ( ) and ( ) ( ) are, respectively, complex Poynting vectors in the Newtonian liquid, elastic surface layer and in the elastic substrate, given by Equation (32). The integrals in Equation (34) can be evaluated analytically. However, due to their excessive length and complexity they will be not reproduced here. Since Poynting vectors ( ) ( ), ( ) ( ) and ( ) ( ) are complex, the average power flow in the direction of propagation is a complex-valued quantity as well and therefore can be written as = + , where in general ≠ 0 and ≠ 0. The active part of the average power flow in the direction of propagation will serve as a normalization factor in in Section 5.1.

Coefficient of Mass Sensitivity
One of the most important parameters characterizing the quality of Love wave sensors is the coefficient of mass sensitivity. Until now, this coefficient of mass sensitivity was evaluated using the perturbation theory and/or numerical methods such as the finite element method, etc. [21][22][23][24][25][26]. In this paper, to evaluate the mass sensitivity, we used analytical methods that provide a deeper insight into the physical background of the operation of Love wave sensors [19,20].
The sensitivity of Love wave sensors to mass loading was initially defined in lossless waveguides, where a semi-infinite Newtonian liquid in Figure 1 is replaced by an infinitesimal layer of mass is lossless it changes only phase velocity of the Love wave, without introducing any extra attenuation, see Figure 2.  The coefficient of mass sensitivity S v p σ of Love wave sensors is defined as: It can be shown [20] that the dispersion equation for Love surface waves propagating in lossless waveguides loaded with an infinitesimal layer of lossless mass density σ is formally identical to Equation (21) developed in Section 3.4 in this paper, providing that the term c As a result, the modified dispersion equation is an implicit function of the phase velocity v p and surface mass density σ, what can be symbolically written as F v p , σ = 0. The derivative dv p /dσ in Equation (35) can be calculated analytically from the modified dispersion equation, using the rules of differentiation of implicit functions. In fact, the differentiation of the modified dispersion equation F v p , σ = 0 with respect to v p and σ leads to the following relation dv p /dσ ∂F/∂v p + (dσ/dσ)∂F/∂σ = 0. Consequently, the derivative dv p /dσ can be written as: By virtue of Equations (35) and (36) the coefficient of mass sensitivity S v p σ is given by the following explicit formula [19,20]: where h 1 is the thickness of the guiding surface layer, q 1 and q 2 are, respectively, transverse wavenumbers of the Love wave in the guiding surface layer and in the substrate and It has to be stressed that Equation (37)  One may wonder what on earth is common between the two waveguide configurations shown in Figures 1 and 2. First waveguide is loaded with a semi-infinite lossy Newtonian liquid ( Figure 1) and second with an infinitesimally thin layer of a lossless mass (Figure 2). The first waveguide is lossy with a non-zero amplitude of the Love wave extending from x 2 = −∞ to x 2 = +∞ and the second is lossless with the amplitude of the Love wave limited to the positive half-space x 2 = (0, +∞).
Our answer to the above question is the following: the common factor shared by these two seemingly dissimilar waveguide configurations is the fact that in both of them the phase velocity v p of the Love wave is affected by the loading medium, i.e., by a semi-infinite lossy Newtonian liquid or an infinitesimal layer of lossless mass. Now, we are in a position to go one step further and identify the suspect physical phenomenon, occurring in both waveguides, which is responsible for the phase delay (velocity changes) of the Love wave propagating in both structures. Since loading with a Newtonian liquid gives rise to a non-zero complex Poynting vector P

Equivalent Inertial Properties of the Newtonian Liquid Seen by Love Surface Waves
The amplitude of Love surface waves, propagating in waveguides loaded with a Newtonian liquid, decays very rapidly in the Newtonian liquid, as a function of the transverse direction x 2 (see Figure 1). The penetration depth δ 0 of the Love wave into a Newtonian liquid equals δ 0 = 1/Req 0 , where the real part Req 0 = a 0 of the transverse wave number q 0 of the Love wave in the Newtonian liquid is given by Equation (29). The penetration depth δ 0 of the Love wave in the Newtonian liquid is given approximately by [27]: where η 0 is the viscosity and ρ 0 the density of the Newtonian liquid. Now, we will answer another question. Can semi-infinite Newtonian liquid, with an exponentially decaying mechanical displacement of the Love wave, be replaced by an equivalent thin layer of a thickness h , with a constant mechanical displacement of the Love wave equal to u 3 (x 2 = 0) that provides the same changes in phase velocity v p of the Love wave?
To answer this question we will determine the average power flow P 2 (x 2 = 0)e 2Re(q 0 )x 2 , see Equation (26), the average power flow in the Newtonian liquid equals: By definition, the penetration depth δ 0 of the Love wave into the loading Newtonian liquid equals δ 0 = 1/Req 0 . Therefore, we can write that: As a result, the sought equivalent thickness h is given by: Since, for thin layers of a lossless mass, the surface mass density of the load can be written as σ = ρ 0 h , we postulate that the effective mass loading exerted by a semi-infinite Newtonian liquid is equivalent to mass loading of a thin layer of lossless mass of thickness δ 0 /2 and density ρ 0 . As a result, we can write: where σ given by Equation (42) represents the equivalent surface mass density [kg/m 2 ] of a semi-infinite Newtonian liquid loading Love wave waveguide. Note that, the equivalent surface mass density σ depends on density ρ 0 and viscosity η 0 of the Newtonian liquid as well as on angular frequency ω of the Love wave. Now, we are able to determine the coefficient of mass sensitivity S v p σ of Love wave sensors, loaded with a semi-infinite Newtonian liquid, using Equation (37) with the effective surface mass density σ given by Equation (42).
The coefficient of mass sensitivity S v p σ will be plotted in Section 5.2, as a function of frequency f and thickness of the elastic surface layer h 1 .

Numerical Results
The analytical formulas obtained in this paper will be illustrated by numerical results performed for the Love wave waveguide structure showed in Figure 1 (Section 2) and Figure 2 (Section 4.1) with the waveguide parameters enclosed in Table 1.
The necessary condition for Love surface waves to exist is as follows: v 1 < v 2 , i.e., the elastic surface layer must be slower than the elastic substrate (see Table 1). In the PMMA elastic surface layer the phase velocity of bulk SH waves v 1 = 1100 m/s and in the ST-Quartz substrate v 2 = 5060 m/s. Both materials are assumed to be lossless and frequency independent. The piezoelectric effect in the ST-Quartz substrate was neglected except for the phase velocity v 2 of bulk SH waves in Quartz (see Appendix A for more details). Since the piezoelectric effect in the Quartz is relatively weak (∼0.3%) the above assumption can be considered as a good first order approximation.  (39)) in the transverse direction x 2 in the Newtonian liquid is a complex quantity; therefore, it can be written as P

Average Power Flow P
The active part ReP In Figures 3 and 4 we present the reactive part ImP   The analytical formulas obtained in this paper will be illustrated by numerical results performed for the Love wave waveguide structure showed in Figure 1 (Section 2) and Figure 2 (Section 4.1) with the waveguide parameters enclosed in Table 1.
The necessary condition for Love surface waves to exist is as follows: < , i.e., the elastic surface layer must be slower than the elastic substrate (see Table 1). In the PMMA elastic surface layer the phase velocity of bulk SH waves = 1100 m/s and in the ST-Quartz substrate = 5060 m/s. Both materials are assumed to be lossless and frequency independent. The piezoelectric effect in the ST-Quartz substrate was neglected except for the phase velocity of bulk SH waves in Quartz (see Appendix A for more details). Since the piezoelectric effect in the Quartz is relatively weak (~0.3%) the above assumption can be considered as a good first order approximation. The average power flow ( ) (Equation (39)) in the transverse direction in the Newtonian liquid is a complex quantity; therefore, it can be written as

Coefficient of Mass Sensitivity
By contrast to the average power flow ( ) , presented in Section 5.1, the coefficient of mass sensitivity is a real-valued quantity, since by definition (Equation (35)) it is related only to changes in phase velocity of the Love wave.
The coefficient of mass sensitivity of Love wave sensors, loaded with a semiinfinite Newtonian liquid, was determined using Equation (37) with the effective surface mass density of the Newtonian liquid given by Equation (42).
The coefficient of mass sensitivity is plotted in Figures 5 and 6, as a function of frequency and thickness of the elastic surface layer ℎ .

Coefficient of Mass Sensitivity S v p σ
By contrast to the average power flow P (0) 2 , presented in Section 5.1, the coefficient of mass sensitivity S v p σ is a real-valued quantity, since by definition (Equation (35)) it is related only to changes in phase velocity of the Love wave.
The coefficient of mass sensitivity S v p σ of Love wave sensors, loaded with a semiinfinite Newtonian liquid, was determined using Equation (37) with the effective surface mass density σ of the Newtonian liquid given by Equation (42).
The coefficient of mass sensitivity S v p σ is plotted in Figures 5 and 6, as a function of frequency f and thickness of the elastic surface layer h 1 .

Discussion
In this paper, we investigated two phenomena occurring in Love wave waveguides loaded with a semi-infinite Newtonian liquid of a viscosity and density . First, the average power flow = ( ) , in the transverse direction , and the corresponding complex Poynting vector ( ) ( = 0), evaluated at the interface = 0, between the waveguide and the Newtonian liquid. Second, the change ∆ ⁄ in phase velocity of the Love wave, represented by the mass coefficient of sensitivity . At first glance, these two phenomena may seem to be completely unrelated. However, in this paper we showed that the first phenomenon provides physical basis for the

Discussion
In this paper, we investigated two phenomena occurring in Love wave waveguides loaded with a semi-infinite Newtonian liquid of a viscosity η 0 and density ρ 0 . First, the average power flow P 2 = ∞ −∞ P 2 (x 2 )dx 2 , in the transverse direction x 2 , and the corresponding complex Poynting vector P (0) 2 (x 2 = 0), evaluated at the interface x 2 = 0, between the waveguide and the Newtonian liquid. Second, the change ∆v p /v p in phase velocity of the Love wave, represented by the mass coefficient of sensitivity S v p σ . At first glance, these two phenomena may seem to be completely unrelated. However, in this paper we showed that the first phenomenon provides physical basis for the second one. In fact, if the average reactive power flow across the interface between the waveguide and the Newtonian liquid, represented by the imaginary part of the Poynting vector