Next Article in Journal
The Development of a Digital Twin Framework for an Industrial Robotic Drilling Process
Next Article in Special Issue
Acoustoelectric Effect for Rayleigh Wave in ZnO Produced by an Inhomogeneous In-Depth Electrical Conductivity Profile
Previous Article in Journal
Near Infrared Emitting Semiconductor Polymer Dots for Bioimaging and Sensing
Previous Article in Special Issue
High-Frequency Surface Acoustic Wave Resonator with Diamond/AlN/IDT/AlN/Diamond Multilayer Structure
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multimode Design and Piezoelectric Substrate Anisotropy Use to Improve Performance of Acoustic Liquid Sensors

1
Kotelnikov Institute of Radio Engineering and Electronics of RAS, 125009 Moscow, Russia
2
Acoustoelectronic and Piezokeramic ELPA Corporation, 124460 Moscow, Russia
3
State Key Laboratory of Mechanics and Control of Mechanical Structures, College of Aerospace Engineering, Nanjing University of Aeronautic and Astronautic, Nanjing 210016, China
4
Piezoelectric Device Laboratory, School of Mechanical Engineering and Mechanics, Ningbo University, Ningbo 315211, China
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(19), 7231; https://doi.org/10.3390/s22197231
Submission received: 28 August 2022 / Revised: 16 September 2022 / Accepted: 20 September 2022 / Published: 23 September 2022
(This article belongs to the Special Issue Piezoelectric Resonator-Based Sensors)

Abstract

:
Using acoustic wave modes propagation in piezoelectric plates loaded with conductive liquids, peculiarities of the mode-liquid acoustoelectric interaction are studied. It is found that (i) in contrast to bulk and surface acoustic waves propagating in piezoelectric semiconductors, the acoustoelectric attenuation of the modes is not symmetric in respect to its maximum, (ii) a large increase in attenuation may be accompanied by a small decrease in phase velocity and vice versa, (iii) the peculiarities are valid for “pure” (without beam steering) and “not pure” (with beam steering) modes, as well as for modes of different orders and polarizations, and (iv) conductivity of test liquid increases electromagnetic leakage between input and output transducers, affecting results of the measurements. To decrease the leakage, the liquid should be localized between transducers, outside the zone over them. If so, the mode sensitivity may be as large as 8.6 dB/(S/m) for amplitude and 107°/(S/m) for phase. However, because of comparable cross-sensitivity towards viscosity and dielectric permittivity, modes with selective detection of liquid conductivity are not found.

1. Introduction

It is well-known that a change in electric conductivity of an adjacent medium located on a surface of a piezoelectric substrate produces a change in phase velocity and attenuation of any acoustic wave with nonzero coupling constant (Rayleigh [1,2,3], Love [3], Lamb [4,5], zero-order shear-horizontal [6,7], higher-order [8], slot [9,10], etc.). For most waves, the phase velocity decreases with conductivity monotonically, approaching its minimum value when tangential components of electric fields accompanying the wave are completely shorted [1,2,3,4,5,6,7,8,9,10]. On the other hand, for particular cases of slightly inhomogeneous piezoactive surface acoustic waves with shear-horizontal polarization (Bleustein-Gulyaev and/or Love waves), an increase in conductivity of a substrate or conductive layer deposited on the substrate produces, first, an increase in phase velocity of the wave until some maximal value and then, a decrease in the velocity to some minimal value [11,12]. This property originates from the fact that for small conductivities, a change in electrical boundary conditions enhances localization of the wave electric fields near the interface and, thereby, enlarges the wave velocity in a propagation medium. The effect of the field localization is larger than the effect of screening the tangential electric wave components up to a certain conductivity. However, when the conductivity is large enough, the screening effect becomes dominant and the phase velocity of the wave behaves in a common manner, monotonically decreasing [11,12].
The acoustoelectric attenuation behaves similar for all types of acoustic waves: it increases for small conductivities, approaches some maximum for intermediate conductivities, and decreases to zero for large conductivities of the substrate or layer [1,2,3,4,5,6,7,8,9,10,11,12]. This property is due to the fact that at low conductivities, the layer behaves like dielectric, while at high conductivities it is more like an ideal conductor. For intermediate conductivities, the structure is characterized by a certain Maxwellian relaxation frequency fM = σSeff, where σS is the surface conductivity of the layer, and εeff is the effective permittivity depending on σS, plate thickness, and wavelength. When acoustic wave frequency f becomes equal to the Maxwellian relaxation frequency fM (f = fM), a resonant interaction between an acoustic wave and a structure takes place and the acoustoelectric attenuation approaches maximum [4,7]. It should be noted that for a substrate coated with a conductive solid layer, the dependence of the attenuation on the surface conductivity is symmetrical. On the contrary, the same dependence for a zero-order shear-horizontal wave in a plate coated with a conductive liquid layer is turned to be asymmetrical [13].
The dependence of acoustic wave properties on a liquid conductivity made it possible to propose a large number of acoustoelectronic sensors. First, of them were based on shear-horizontal surface acoustic waves characterized by maximal value of the mechanical displacement component lying in the plane perpendicular to the wave vector [14,15,16,17]. Further, it was suggested to use for the same sensors the zero-order waves with shear-horizontal polarization propagating in the thin plate whose thickness is comparable with the acoustic wavelength. The waves of this type are characterized by the largest electromechanical coupling coefficients among all other waves in the same materials [7,13,18]. Furthermore, the spectrum of higher-order waves in piezoelectric plates with a thickness comparable to or greater than the wavelength is very wide and the number and properties of the plate waves are purposefully varied by the plate thickness h, wavelength λ, mode order n, and propagation direction [19,20,21,22,23,24,25,26]. Electric fields of the modes penetrate into contacting liquid and interact with free charge carriers. This interaction results in a change in the phase velocities v (phases ϕ) and amplitudes (attenuation α) of the waves. Although the acoustoelectronic interaction takes place only in a thin liquid layer (~10 μm) adjacent to the interface and only in a narrow range of liquid conductivities (σ ≈ 0.01 ÷ 10 S/m [7]), first-order plate modes have proven to be useful for applications such as the detection of yeasts and bacteria [27,28] and registration of the water-ice phase transition [20,29]. The propagation of shear-horizontal waves in the conductive liquid-loaded functionally graded porous piezoelectric media (FGPPM) has been investigated recently. The dispersion and the attenuation curves have been plotted to describe the impact of various parameters, namely gradient parameter of both the piezoelectric material, thickness of FGPPM plate, imperfect interface parameter, relative permittivity, and conductivity of the liquid medium [30]. A thickness shear solidly mounted resonator based on Yttrium-doped AlN has been suggested recently as a high-sensitive liquid sensor [31]. A lateral electric field-excited resonator based on piezoceramic has also been suggested as a sensor for liquid conductivity measurements [32].
Attractive property of the acoustic plate waves is also the change in the power flow direction through the change in the electrical boundary conditions and propagation direction. This property, originated from anisotropy of the plate, has already been demonstrated in [33,34].
In conclusion, expansion of the acoustic plate wave applications requires the study of acoustoelectronic interaction in a more general case for the whole spectrum of the waves including modes of high orders and modes with beam steering. It is also interesting to estimate achievable sensitivity of the wave towards electric conductivity and compare it with those for liquid viscosity and dielectric permittivity.
The purpose of this paper is the experimental studying of the effect of liquid conductivity on the characteristics of the higher-order waves in piezoelectric plates, as well as the use of piezoelectric plate anisotropy to improve the performance of the liquid sensors.

2. Materials and Methods

Schematic view of a test sample is shown in Figure 1. It consists of 4 delay lines with input and output interdigital transducers (IDTs) arranged on same piezoelectric plate at different angles (Θ = 0°, 30°, 60°, 90°) with respect to crystallographic axis X (Table 1). The number of acoustic plate modes generated in each line is about 10–15. The total amount of modes detected in the plate is 40–60. Properties of all modes are different. They are varied with mode order n, propagation direction (angle Θ), material of the plate, plate thickness h, and wavelength λ (period of the transducers).
The plate materials used in experiments are commercially available 128Y-LiNbO3 (Euler angles 0°, 37.86°, Θ), 41Y- LiNbO3 (0°, −49°, Θ), Y-LiNbO3 (0°, 90°, Θ), and 36Y-LiTaO3 (0°, −54°, Θ) with thickness h = 350 and 500 μm. The plates have one grinded and one polished surface. The grinded surface (optical class Δ10) has averaged horizontal and vertical roughness 0.16 μm and 0.8 μm, respectively. The polished surface (optical class Δ14) has averaged horizontal and vertical roughness 0.01 μm and 0.05 μm, respectively. The polished surface is coated with liquid cell (teflon) 18 mm in diameter L and four pairs of IDTs at (Θ = 0°, 30°, 60°, 90° off X axis), as shown in Figure 1. Volume of the liquid sufficient for the measurements is as small as 100 μL.
Each transducer comprises 20 finger electrodes followed with period λ = 200 μm patterned from 100-nanometer thick Cr and 1000-nanometer thick Al. Large number of electrodes pairs (20) provides narrow pass band of the transducers (5%) and good frequency resolution of the modes (±0.5 MHz) with close velocities vn (±0.2 × 106 m/s). The normalized plate thickness equal to h/λ = 350 μm/200 μm = 1.75 and 500 μm/200 μm = 2.5 ensure large number of the modes in each direction and the mode’s variety. Different orientations of the transducers allow to compare “pure” directions (Θ = 0°, 90°), where acoustic beams of all modes are parallel to propagation direction (no beam steering, the power flow angle Ψn = 0°), and “not pure” directions (Θ = 30°, 60°), where there is beam steering for all modes (Ψn ≠ 0°).
The angle Ψn for acoustic plate modes is calculated for other acoustic waves [35]:
Ψn = arctg(1/vn) (dvn/dΘ),
where dvn is a change in a mode velocity due to a change (dΘ) in propagation direction, and Ψn is an angle between a mode beam and a propagation direction on a plate face. Both Θ and Ψn are taken in radians in (1).
Calculations of Ψn are accomplished for the first 10 modes (n = 0–9), using well-approved software [36] and material constants [35]. In order to calculate orientation dependence Ψn(Θ) for each mode, the angle Θ between propagation direction and X axis is changed from Θ = 0° (+X axis) to Θ = 180° (–X axis) with the step 3°, holding other two Euler angles constant. In order to follow a given mode and avoid a skip from one mode to the other, the depth profiles of the modes are controlled. It allows for observing a smooth evolution of any investigated mode during the calculations for different Θ. Like for surface and bulk acoustic waves, the as-calculated dependence Ψn(Θ) is originated from anisotropy of the crystals, but unlike waves of other types, this dependence for acoustic plate modes is changed also with mode order n and plate thickness h/λ [37]: for a given mode n and angle Θ, the angle Ψn(Θ) is varied with plate thickness h/λ; for a given plate thickness h/λ and angle Θ, the angle Ψn(Θ) is varied with mode order n.
Acoustic delay line is placed in climatic chamber UC-20CE (Terchy Technology LTD, Taiwan, 20l in volume) fixed at T = 20 ± 0.1 °C (293.15 K). From the chamber, the sample is connected with KEYSIGHT 5061B network analyzer (Keysight, Santa Rosa, CA, USA) operating in amplitude (S21) or phase (ϕ21) formats. First, insertion loss (transfer function S21) of the delay line versus frequency f is measured without any liquid in the cell (S21AIR). This measurement provides the spectrum of acoustic modes existing in plate at relevant frequencies fn = vn/λ. Second, the measurements of S21 and ϕ21 are carried out, when reference liquid (distilled water, 200 mL) and test liquid (see below) are introduced into the cell one by one. It provides modes with best sensitivity towards a test liquid, i.e., the modes with large amplitude response ΔS12 = S12 LqS1H2O and/or large phase response Δϕ12 = ϕ12Lq − ϕ12H2O. Finally, the time variations in the amplitude ΔS12(t) and phase Δϕ12(t) responses are measured for the most sensitive modes when drops of test liquid (40 mL) are induced one by one into the cell with distilled water at regular time intervals Δt ≈ 100 sec. Knowing mass of the drop, drops number, mass of distilled water in the cell, and relevant acoustic responses, calibration curves for a given mode and a given liquid are plotted. Experimental results for different modes, plates, and propagation directions are compared with one another and best of them are tabulated. Precisions of the measurements are ±0.1 dB for S12 and ±0.1 degree for phase. The corresponding flow chart of the experiments and experimental setup is presented in Figure 2.
Absolute values of the phases are not useful for comparing different samples because electric contributions to the phases are changed from sample to sample unpredictably depending on length of wires, quality of contacts, location of sample in holder, etc. Comparison is useful only for phase responses Δϕ = ϕ Lq − ϕH2O, when one and the same electric contributions to ϕLq and ϕH2O mutually compensate one another for each sample.
Distilled water is used as reference liquid because of its low conductivity (σ < 10−3 S/m) and small viscosity (η = 1.03 cP).
Water solutions of NaCl are used as test liquids with variable conductivity changed from σ = 0 (distilled water) to σ = 10 S/m (7.6 weight % NaCl in water) by changing weight per cent of the solution components. Viscosity (<13%), density (<8%), and permittivity (<1%) of the solutions are almost constant [38] and do not contribute much to the conductivity responses.
The cross-sensitivity of the modes towards viscosity η is studied using water solutions of glycerin as test liquid. Viscosities η of the solutions are varied from 1 to 1490 cP, while conductivity σ < 10−3 S/m, density ρ, and permittivity εlq stay permanent within ±10% and give negligible contributions to acoustic responses for conductivity and viscosity.
Numerical data for σ, η, ρ, and εlq of the solutions are determined from anhydrous solute weight per cent of the solution components [38]. Fabrication of the solutions is accomplished by mixing partial components in forced vibrator for about 5 min. The error in weight concentrations is about ±1%.
Cross-sensitivity of the samples towards temperature is studied by heating climatic camera UC-20CE from T = 0 to 55 °C with the step ΔT = 5 °C. The measurements are accomplished for amplitude response ΔS12 whose temperature variations without liquid in the cell are small (~0.1 dB) [39]. On the contrary, when liquid is present, the temperature response depends on liquid viscosity η at relevant temperature T. Extracting 1st data set from the 2nd, the temperature dependence of viscosity η(T) for a given liquid is measured. Precision of the measurements is ± 20%.
Results of the measurements are compared with numerical data calculated from simple formulas applied before to shear-horizontal acoustic modes in quartz plates loaded with conductive liquids [3]:
α β = k 2 2 ε p l + ε 0 ε p l + ε l q σ ω ( ε p l + ε l q ) σ 2 + ω 2 ( ε p l + ε l q ) 2  
Δ v v 0 = k 2 2 ε p l + ε 0 ε p l + ε l q σ 2 σ 2 + ω 2 ( ε p l + ε l q ) 2  
where vo is the velocity of the wave at zero liquid conductivity, Δv is the change in the velocity arising from liquid conductivity, α is the acoustoelectric attenuation, β is the wave number, k2 is the coupling constant of the wave, ω = 2πf is the wave frequency, and εlq and εpl are dielectric constants of the liquid and the plate, respectively. In particular, for quartz plates with small k2 = 3.2 × 10−4, acoustoelectric interaction at ω = 109 s−1 produces quite measurable change in wave velocity or phase Δv/vo = Δϕ/ϕ ≈ 20 ppm, but negligible change in wave amplitude (acoustoelectric attenuation α = ΔS12/L ≈ 0.03 dB/mm). Because of that, the value of the attenuation for quartz plates was neither measured nor calculated [3]. However, substrates exhibiting higher coupling constant k2, such as lithium niobate, display grater attenuation, as shown in present paper.
Though formulas (2) and (3) directly relate acoustic (v, α) and electric (σ, εlq) characteristics, they do not account for the difference in the modes of different orders, while the mode-liquid interaction should depend on electric potentials of the modes within conductive liquid. Therefore, results of the calculations may be compared with experimental data only qualitatively. In present paper, it is made for one of the modes detected in lithium niobate plate, as an example (vo = 14.093 m/s, k2 = 116 × 10−4, ω =109 s−1).

3. Results and Discussion

Results obtained are presented in Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9 and Figure 10 and in Table 2 and Table 3.
According to formulas (2) and (3), both acoustic velocity v and attenuation α increase with conductivity σ of the liquid, while they decrease with liquid permittivity εlq (Figure 3). The changes produced by σ and εlq to v and α are comparable with each other.
In contrast to bulk and surface acoustic waves in piezoelectric semiconductors [40], acoustoelectric attenuation for acoustic plate modes and conductive liquids is not symmetric in respect to its maximum (Figure 3b): at the beginning, α increases with σ quickly, then approaches maximum at about σ ≈ 0.5 S/m, and finally, decreases slowly for σ > 0.5 S/m. The nature of the property is not clear yet, but very probable that it is originated from inhomogeneity of the propagation medium composed of a plate of one material and a liquid of another material and from inhomogeneity of the mode profiles through the depth of the propagation medium. Slow decrease in α after its maximum makes the range of wave sensitivity towards liquid conductivity σ wider.
Figure 4 and Figure 5 demonstrate an important peculiarity of acoustic sensing conductive liquids: results of the measurements depend on the location of the liquid on the plate. When the liquid is over the whole propagation path, including two regions opposite to IDTs (Figure 4b, insert), insertion loss S21(f) are damaged by large electromagnetic leakage (horizontal level in Figure 4a) and phase ϕ is decreased with conductivity σ (Figure 4b). As far as ϕ~1/v, this behavior of the phase indicates that partial shorting of the plate face produced by the conductive liquid results in an increase in mode velocity v that is not right [7,41]. On the other hand, when the conductive liquid is only between transducers, outside the zone over them (Figure 4b), the leakage is much lower (15 dB below acoustic signal), insertion loss S12(f) is not damaged (Figure 4a), and phase ϕ is increased with conductivity σ (Figure 5), indicating a correct decrease in the mode velocity v in accordance with appropriate boundary conditions [7,41]. Following this peculiarity, all measurements of the present paper are accomplished with liquids localized, as shown in Figure 5.
Figure 6 shows typical spectra of acoustic plate modes in crystal plates with different propagation directions (angles Θ). For Θ = 0° (X axis) and Θ = 90° (perpendicular to X axis), all modes are “pure”: the beams of the modes are parallel to the wave vector and the power flow angles Ψn are zero for all modes.
On the contrary, for propagation directions at Θ = 30° and Θ = 60° off the X axis, the beams are steering (“not pure” directions) and the power flow angles Ψn are not zero for most modes. The values of the angles Ψn are ranged from −13° to +13° depending on the mode order n (Figure 7, Table 2). As a result, each IDT radiates the modes not only at different frequencies (because fn = vn/λ, where λ is the same IDT period for all modes, while vn is different for all modes), but at various angles Ψn off the wave vector β, i.e., like a fan (Figure 8). In this case, some part of the beam’s energy misses output transducers, increasing the insertion loss S12 of the delay lines. For example, the mode n = 0 at Θ = 0° has Ψn = 0° and S12(AIR) = 30 dB (Figure 6a), while the same mode at Θ = 30° has Ψn = −13° and S12(AIR) = 56 dB (Figure 6b). Nevertheless, “not pure” modes may also be used for liquid sensing, in particular, for electric measurements. The modes have different responses Δϕ12 = ϕ12Lq − ϕ12H2O and ΔS12 = S12LqS12H2O for all propagation directions. The best of them for each direction is indicated in Figure 6 by bold arrows.
Figure 9 shows typical behavior of the mode responses versus step-by-step increase in liquid conductivity. At the beginning (t < 0 s), when a cell contains distilled water as a reference liquid with σ = 0 S/m, both the amplitude and phase of mode are permanent (S12, ϕ12 = constants). When the first drop of conductive liquid is introduced into the water (first vertical arrow), the conductivity of the water-drop solution σ is increased from 0 to 2.8 S/m and both responses (ΔS1, Δϕ1) are changed immediately and remarkably: the amplitude of the mode is decreased (the attenuation is increased) and the phase of the mode is increased (the velocity is decreased). Next drops of the conductive liquid (other four vertical arrows from left to right) result in a further increase in conductivity σ of the liquid solution from 2.8 to 6.8 S/m and decreases in mode attenuation (ΔS4/4 in average) and velocity (Δϕ4/4 in average). The same behavior is inherent for all modes, plates, and propagation directions studied in the paper. The most attractive of them is presented in Table 3.
Table 3 shows that (i) the modes with enhanced electric sensitivity exist both for “pure” (Θ = 0°, 90°) and “not pure” (Θ = 30°, 60°) directions, (ii) the amplitude responses of the modes approach ΔS1 = 24 dB or 8.6 dB/(S/m) for the first drop and ΔS4/4 = 2.1 dB or 0.5 dB/(S/m) for the next four drops on average, (iii) relevant phase responses are as large as Δϕ1/4 = 300° or 107°/(S/m) and (Δϕ4/4) = 21.3° or 5.3°/(S/m), (iv) large amplitude responses may be accompanied by small phase responses and vice versa, though according to formulas (2) and (3), both acoustic characteristics should increase and decrease synchronously following the coupling constant k2 of the mode, and (iv) there are few modes with a good combination of all four sensing parameters (ΔS1, ΔS4/4, Δϕ1/4, and Δϕ4/4). Two of them are found along the Θ = 60° and Θ = 90° propagation directions (bold) in one and the same Y-LiNbO3 plate.
Figure 10 presents results of a mode calibration: (i) the curves for viscosity η (Figure 10a,b) are typical for waves vibrating at ultrasonic frequencies [3,42]. They are almost linear for small η, when liquid behaves as ideal (Newtonian), and saturates for large η, when liquid behaves as a solid, (ii) the curves for conductivity σ (Figure 10c,d) are also typical, in general [7,13]. Velocity of the mode decreases with σ monotonously, while acoustic attenuation jumps up initially, approaches maximum, and finally drops down slowly to zero for large σ, (iii) at the same time, the acoustoelectric attenuation has an important peculiarity. In contrast to surface and bulk acoustic waves propagating in piezoelectric semiconductors where the curve ΔS12(σ) is symmetric with respect to its maximum [40], the same curve for acoustic plate modes and conductive liquid is asymmetric (Figure 10c). It makes the range of detectable values wider: for phase response, the range of the measurements is restricted by about 5 S/m (Figure 10d), while for amplitude response, it is prolonged to 10–15 S/m (Figure 10c). This property is in qualitative agreement with numerical data calculated with formulas (2) and (3) (Figure 10e,f), but the nature of the property is not clear yet as already mentioned, (iv) because of the particular dependence on liquid conductivity σ (Figure 10d), the phase measurement provides a simple value of σ because each ϕ corresponds to a single σ; on the other hand, the phase response is temperature-dependent and demands careful thermal stabilization. Vice versa, the amplitude response ΔS12 is almost temperature-independent, but because of particular dependence on liquid conductivity σ (Figure 10c), it gives two different values σ for each ΔS12 at the same time and, thereby, demands additional measurements to avoid the wrong value. The maximal range of acoustic measurements is 0–10 S/m for phase output and 0–20 S/m for amplitude output. The volume of the test probe is about 100 μL in both cases, and (v) the sensitivity of the mode towards liquid conductivity (Figure 10c,d) is comparable with that for viscosity (Figure 10a,b), in general. Modes providing selective detection of conductivity, but not viscosity are not found in the paper.

4. Conclusions

Three peculiarities of acoustoelectric interaction between acoustic plate modes and conductive liquids deposited on piezoelectric plates are found in the paper: electromagnetic leakage between input and output transducers increases remarkably because of liquid conductivity, large increase in acoustoelectric attenuation may be accompanied by small decrease in phase velocity, and vice versa, the modes attenuation versus liquid conductivity is not symmetrical with respect to its maximum, making the range of detectable values wider. All peculiarities of the interaction are inherent for modes of different orders, polarizations, and beam steering. The sensitivity of the modes towards liquid conductivity is as large as 8.6 dB/(S/m) for amplitude and 107°/(S/m) for phase. The range of the measurements is restricted by about 10–15 S/m.
Because of anisotropy of the plates, there is beam steering for all modes. The magnitude of the steering depends on the mode order. The beam of different modes generated at various frequencies forms a fan around propagation direction on the surface of the plate. For some modes, the power flow angles are as large as ±13°.
Acoustic measurements of liquid conductivity through the phase of the waves are preferable when temperature of the liquid is constant. On the other hand, when the liquid temperature is varied somehow, the measurements of the wave amplitudes are more suited, though they give two different conductivity values for any acoustic amplitude. In this case, additional efforts to avoid a wrong value and to find the correct counterpart are required.
Responses of the modes towards liquid conductivity, viscosity, and dielectric permittivity found in the paper are comparable with one another in general.

Author Contributions

Conceptualization, V.A. and I.K.; methodology, N.V.; software, Z.Q., P.L. and H.E.; validation, E.S., A.S., and I.K.; formal analysis, T.M.; investigation, V.A.; resources, I.K.; data curation, I.K.; writing—original draft preparation, V.A.; writing—review and editing, I.K.; visualization, A.S.; supervision, V.A.; project administration, I.K.; funding acquisition, I.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Russian Science Foundation grant #20-19-00708.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wixforth, A.; Scriba, J.; Wassermeier, J.; Kotthaus, G.; Weimann, G.; Schlapp, W. Interaction of surface acoustic waves with a two-dimensional electron system in LiNbO3-GaAs/AlGaAs sandwich structure. J. Appl. Phys. 1988, 64, 2213–2215. [Google Scholar] [CrossRef]
  2. Zaitsev, B.D.; Kuznetsova, I.E. Surface acoustic waves in a gallium arsenide-conducting layer structure. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 1996, 43, 1130–1132. [Google Scholar] [CrossRef]
  3. Ballantine, D.S.; White, R.M.; Martin, S.J.; Ricco, A.J.; Zellers, E.T.; Frye, G.C.; Wohltjen, H. Acoustic Wave Sensors; Academic Press: San Diego, CA, USA, 1997; pp. 83–88. [Google Scholar]
  4. Zaitsev, B.D.; Joshi, S.G.; Kuznetsova, I.E.; Borodina, I.A. Influence of conducting layer and conducting electrode on acoustic waves propagating in potassium niobate plates. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2001, 48, 624–626. [Google Scholar] [CrossRef] [PubMed]
  5. Zhu, F.; Wang, B.; Qian, Z.; Kuznetsova, I.; Ma, T. Influence of surface conductivity on dispersion curves, mode shapes, stress, and potential for Lamb waves propagating in piezoelectric plate. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2020, 67, 855–862. [Google Scholar] [CrossRef] [PubMed]
  6. Kuznetsova, I.; Nedospasov, I.A.; Kolesov, V.V.; Qian, Z.; Wang, B.; Zhu, F. Influence of electrical boundary conditions on profiles of acoustic field and electric potential of shear-horizontal acoustic waves in potassium niobate plates. Ultrasonics 2018, 86, 6–13. [Google Scholar] [CrossRef]
  7. Zaitsev, B.D.; Kuznetsova, I.E.; Joshi, S.G.; Borodina, I.A. Shear horizontal acoustic waves in piezoelectric plates bordered with conductive liquids. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2001, 48, 627–631. [Google Scholar] [CrossRef]
  8. Anisimkin, V.I.; Caliendo, C.; Verona, E. Acoustic Plate Mode sensing in liquids based on free and electrically shorted plate surfaces. Ultrasonics 2016, 68, 29–32. [Google Scholar] [CrossRef]
  9. Borodina, I.; Zaitsev, B.; Teplykh, A.; Burygin, G.; Guliy, O. Sensor based on PZT ceramic resonator with lateral electric field for immunodetectionof bacteria in the conducting aquatic environment. Sensors 2020, 20, 3003. [Google Scholar] [CrossRef]
  10. Guliy, O.I.; Zaitsev, B.D.; Burygin, G.L.; Karavaeva, O.A.; Fomin, A.S.; Staroverov, S.A.; Borodina, I.A. Prospects for the use of gold nanoparticles to increase the sensitivity of an acoustic sensor in the detection of microbial cells. Ultrasound Med. Biol. 2020, 46, 1727–1737. [Google Scholar] [CrossRef]
  11. Zaitsev, B.D.; Kuznetsova, I.E.; Joshi, S.G. Anomalous resisto-acoustic effect. J. Appl. Phys. 1999, 86, 6868–6874. [Google Scholar] [CrossRef]
  12. Kuznetsova, I.E.; Zaitsev, B.D. The peculiarities of the Bleustein-Gulyaev wave propagation in structures containing conductive layer. Ultrasonics 2015, 59, 45–49. [Google Scholar] [CrossRef] [PubMed]
  13. Zaitsev, B.D.; Kuznetsova, I.E.; Joshi, S.G.; Borodina, I.A. Acoustic waves in piezoelectric plates bordered with viscous and conductive liquids. Ultrasonics 2001, 39, 45–50. [Google Scholar] [CrossRef]
  14. Josse, F. Acoustic wave liquid-phase-based microsensors. Sens. Actuators A Phys. 1994, 44, 199–208. [Google Scholar] [CrossRef]
  15. Fisher, B.H.; Malocha, D.C. Study of the acoustoelectric effect for SAW sensors. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2010, 57, 698–706. [Google Scholar] [CrossRef]
  16. Caliendo, C.; Hamidullah, M. A theoretical study of love wave sensors based on ZnO-glass layered structures for application to liquid environments. Biosensors 2016, 6, 59. [Google Scholar] [CrossRef]
  17. Kondoh, J.; Nakayama, K.; Kuznetsova, I. Study of frequency dependence of shear horizontal surface acoustic wave sensor for engine oil measurements. Sens. Actuators A Phys. 2021, 325, 112503. [Google Scholar] [CrossRef]
  18. Yang, C.-H.; Shue, C.J. Guided waves propagating in a piezoelectric plate immersed in a conductive fluid. NDT E Int. 2001, 39, 199–206. [Google Scholar] [CrossRef]
  19. Lee, Y.-C.; Kuo, S.H. Leaky Lamb wave of a piezoelectric plate subjected to conductive fluid loading: Theoretical analysis and numerical calculation. J. Appl. Phys. 2006, 100, 073519. [Google Scholar] [CrossRef]
  20. Anisimkin, V.; Kolesov, V.; Kuznetsova, A.; Shamsutdinova, E.; Kuznetsova, I. An analysis of the water-to-ice phase transition using acoustic plate waves. Sensors 2021, 21, 919. [Google Scholar] [CrossRef]
  21. Anisimkin, I. New type of an acoustic plate modes: Quasi-longitudinal normal wave. Ultrasonics 2004, 42, 1095–1099. [Google Scholar] [CrossRef]
  22. Yantchev, V.; Katardjiev, I. Micromachined thin film plate acoustic resonators utilizing the lowest order symmetric Lamb wave mode. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2007, 54, 87–95. [Google Scholar] [CrossRef] [PubMed]
  23. Soluch, W.; Lysakowska, M. Properties of shear horizontal acoustic plate modes in BT-cut quartz. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2011, 58, 2239–2243. [Google Scholar] [CrossRef] [PubMed]
  24. Caliendo, C.; Lo Castro, F. Quasi-linear polarized modes in Y-rotated piezoelectric GaPO4 plates. Crystals 2014, 4, 228–240. [Google Scholar] [CrossRef]
  25. Chen, Z.; Fan, L.; Zhang, S.; Zhang, H. Theoretical research on ultrasonic sensors based on high-order Lamb waves. J. Appl. Phys. 2014, 115, 14–20. [Google Scholar] [CrossRef]
  26. Anisimkin, V.I.; Voronova, N.V. Features of normal higher-order acoustic wave generation in thin piezoelectric plates. Acoust. Phys. 2020, 66, 1–4. [Google Scholar] [CrossRef]
  27. Anisimkin, V.I.; Pokusaev, B.G.; Skladnev, D.A.; Sorokin, V.V.; Tyupa, D.V. Application of an acoustoelectronic technique to study ordered microstructured disperse systems with biological objects in a hydrogel. Acoust. Phys. 2016, 62, 754–759. [Google Scholar] [CrossRef]
  28. Anisimkin, V.I.; Kuznetsova, I.E.; Kolesov, V.V.; Pyataikin, I.I.; Sorokin, V.V.; Skladnev, D.A. Plate acoustic wave sensor for detection of small amounts of bacterial cells in micro-litre liquid samples. Ultrasonics 2015, 62, 156–159. [Google Scholar] [CrossRef]
  29. Anisimkin, V.I.; Voronova, N.V. New modification of the acoustic Lamb waves and its application for liquid and ice sensing. Ultrasonics 2021, 116, 106496. [Google Scholar] [CrossRef]
  30. Kumaru, S.; Sahu, S.A. SH wave vibration in functionally graded porous piezoelectric composite loaded with Newtonian conductive/non-conductive liquid. Waves Rand. Compl. Media 2022, 1–23. [Google Scholar] [CrossRef]
  31. Laidoudi, F.; Kanouni, F.; Assali, A.; Caliendo, C.; Amara, S.; Nezarri, H.; Boubenider, F. Thickness shear SMRresonator based on Yttrium-doped AlN for high sensitive liquid sensor. Sens. Actuators A Phys. 2022, 333, 113238. [Google Scholar] [CrossRef]
  32. Borodina, I.A.; Zaitsev, B.D.; Teplykh, A.A. Influence of the conductivity of a liquid contacting with a lateral electric field excited resonator based on PZT ceramics on its characteristics. Ultrasonics 2020, 102, 106059. [Google Scholar] [CrossRef] [PubMed]
  33. Anisimkin, V.I.; Voronova, N.V. Controlling Acoustic Wave Energy Fluxes in Piezoelectric Crystals. J. Commun. Technol. Electron. 2020, 65, 90–92. [Google Scholar] [CrossRef]
  34. Kuznetsova, I.E.; Zaitsev, B.D.; Teplykh, A.A.; Joshi, S.G.; Kuznetsova, A.S. Power flow angle (PFA) of acoustic waves in thin piezoelectric plates. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2008, 55, 1984–1991. [Google Scholar] [CrossRef]
  35. Slobodnik, A.J.; Conway, J.R.; Delmonico, E.D. Microwave Acoustic Handbook; V.1A. Surface Wave Velocities; Air Force Systems Command: Cambridge, MA, USA, 1973. [Google Scholar]
  36. Adler, E.L.; Slaboszewicz, J.K.; Farnell, G.W.; Jen, C.K. PC software for SAW propagation in anisotropic multi-layers. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 1990, 37, 215–220. [Google Scholar] [CrossRef]
  37. Anisimkin, V.I. Anisotropy of the acoustic plate modes in ST-Quartz and 128Y-LiNbO3. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2014, 61, 120–132. [Google Scholar] [CrossRef]
  38. Weast, R.C.; Astle, M.J. Chemical Rubber Company Handbook of Chemistry and Physics, 66th ed.; Chemical Rubber D: Boca Raton, FL, USA, 1985; 254p. [Google Scholar]
  39. Anisimkin, V.I.; Voronova, N.V.; Puchkov, Y.V. General properties of the acoustic plate modes at different temperatures. Ultrasonics 2015, 62, 46–49. [Google Scholar] [CrossRef] [PubMed]
  40. Gulyaev, Y.V.; Pustovoit, V.I. Amplification of surface waves in semiconductors. Sov. Phys. J. Exp. Theor. Phys. 1964, 47, 2251–2253. [Google Scholar]
  41. Zaitsev, B.D.; Kuznetsova, I.E.; Joshi, S.G. Anomalous Resisto-acoustic Effect in a Piezoelectric–Conducting Liquid Structure. Tech. Phys. 2001, 46, 767–769. [Google Scholar] [CrossRef]
  42. Anisimkin, I.V.; Anisimkin, V.I. Attenuation of acoustic normal modes in piezoelectric plates loaded by viscous liquids. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2006, 53, 1487–1492. [Google Scholar] [CrossRef]
Figure 1. Schematic view (a) and photo (b) of experimental sample with 4 acoustic delay lines arranged on same piezoelectric plate; 1—plate; 2—interdigital transducers (IDTs); 3—liquid cell. Acoustic plate modes propagate in the plate at angles 0°, 30°, 60°, and 90° off X axis.
Figure 1. Schematic view (a) and photo (b) of experimental sample with 4 acoustic delay lines arranged on same piezoelectric plate; 1—plate; 2—interdigital transducers (IDTs); 3—liquid cell. Acoustic plate modes propagate in the plate at angles 0°, 30°, 60°, and 90° off X axis.
Sensors 22 07231 g001
Figure 2. Flow chart of the experiments and schematic view of the experimental setup.
Figure 2. Flow chart of the experiments and schematic view of the experimental setup.
Sensors 22 07231 g002
Figure 3. Changes in phase velocity Δv/vo (a) and acoustoelectric attenuation α/β (b) versus liquid conductivity σ and dielectric permittivity εlq calculated with formulas (2) and (3). Liquids: εlq/ε0 = 79.3 (1), 52.5 (2), 32.9 (3), and 24.2 (4). Plate: 128°YX-LiNbO3, h = 500 μm, εpl = 34ε0. Wave: λ = 300 μm, v = 14.093 m/s, k2 = 116 × 10−4, ω =109 s−1, displacements on the plate faces are {u1, u2, u3,} = {1, 0.36, 0.34}. β is the wave number.
Figure 3. Changes in phase velocity Δv/vo (a) and acoustoelectric attenuation α/β (b) versus liquid conductivity σ and dielectric permittivity εlq calculated with formulas (2) and (3). Liquids: εlq/ε0 = 79.3 (1), 52.5 (2), 32.9 (3), and 24.2 (4). Plate: 128°YX-LiNbO3, h = 500 μm, εpl = 34ε0. Wave: λ = 300 μm, v = 14.093 m/s, k2 = 116 × 10−4, ω =109 s−1, displacements on the plate faces are {u1, u2, u3,} = {1, 0.36, 0.34}. β is the wave number.
Sensors 22 07231 g003
Figure 4. Insertion loss S12 (a) and phase ϕ (b) measured for AIR (without liquid) and conductive liquids (H2O + NaCl (5%)) with σ from 0 to 6 S/m deposited over the whole propagation path (30 mm). Plate: 128°YX-LiNbO3, h = 500 μm. Mode f = 63.9 MHz, λ = 200 μm. Test liquids: distilled water H2O (200 mg) and distilled water (200 mg) with drops of 0.9% NaCl (40 mg) in the water (arrows).
Figure 4. Insertion loss S12 (a) and phase ϕ (b) measured for AIR (without liquid) and conductive liquids (H2O + NaCl (5%)) with σ from 0 to 6 S/m deposited over the whole propagation path (30 mm). Plate: 128°YX-LiNbO3, h = 500 μm. Mode f = 63.9 MHz, λ = 200 μm. Test liquids: distilled water H2O (200 mg) and distilled water (200 mg) with drops of 0.9% NaCl (40 mg) in the water (arrows).
Sensors 22 07231 g004
Figure 5. Insertion loss S12 (a) and phase ϕ (b) measured for AIR (without liquid) and conductive liquids (H2O + NaCl (5%)) with σ from 0 to 6 S/m deposited only between transducers (18 mm). Plate: 128°YX-LiNbO3, h = 500 μm. Mode f = 63.9 MΓц, λ = 200 μm. Test liquid: distilled water H2O (200 mg) and distilled water (200 mg) with drops of 0.9% NaCl (40 mg) in the water (arrows).
Figure 5. Insertion loss S12 (a) and phase ϕ (b) measured for AIR (without liquid) and conductive liquids (H2O + NaCl (5%)) with σ from 0 to 6 S/m deposited only between transducers (18 mm). Plate: 128°YX-LiNbO3, h = 500 μm. Mode f = 63.9 MΓц, λ = 200 μm. Test liquid: distilled water H2O (200 mg) and distilled water (200 mg) with drops of 0.9% NaCl (40 mg) in the water (arrows).
Sensors 22 07231 g005
Figure 6. Insertion loss S12 of the delay lines with different propagation directions measured for various loadings on the path (L = 18 mm). Plate: 128°Y-LiNbO3, h/λ = 1.75 (h = 350 μm, λ = 200 μm). Propagation directions: Θ = 0° (a), 30° (b), 60° (c), and 90° (d) off X axis. Loadings: dashed lines correspond to case absence any liquid (AIR); black lines correspond to situation with presence of distilled water (σ = 0 S/m, 200 mg); red lines correspond to situation with presence of NaCl water solution (σ = 6 S/m, 200 mg). Bold arrows correspond to the modes with largest sensitivity towards liquid conductivity for each propagation direction.
Figure 6. Insertion loss S12 of the delay lines with different propagation directions measured for various loadings on the path (L = 18 mm). Plate: 128°Y-LiNbO3, h/λ = 1.75 (h = 350 μm, λ = 200 μm). Propagation directions: Θ = 0° (a), 30° (b), 60° (c), and 90° (d) off X axis. Loadings: dashed lines correspond to case absence any liquid (AIR); black lines correspond to situation with presence of distilled water (σ = 0 S/m, 200 mg); red lines correspond to situation with presence of NaCl water solution (σ = 6 S/m, 200 mg). Bold arrows correspond to the modes with largest sensitivity towards liquid conductivity for each propagation direction.
Sensors 22 07231 g006
Figure 7. Orientation dependence of the power flow angles Ψn for the first five modes propagating in 128°Y-LiNbO3 plate with normalized thickness h/λ = 1.75 (h = 350 μm, λ = 200 μm). The arrows indicate the angles Ψn for direction Θ = 30° off X axis.
Figure 7. Orientation dependence of the power flow angles Ψn for the first five modes propagating in 128°Y-LiNbO3 plate with normalized thickness h/λ = 1.75 (h = 350 μm, λ = 200 μm). The arrows indicate the angles Ψn for direction Θ = 30° off X axis.
Sensors 22 07231 g007
Figure 8. Schematic view of the fan-shaped power flow radiation from IDT on a plate surface for Lamb waves of different orders n (according to Table 2).
Figure 8. Schematic view of the fan-shaped power flow radiation from IDT on a plate surface for Lamb waves of different orders n (according to Table 2).
Sensors 22 07231 g008
Figure 9. Time variations in the insertion loss S12 (mode amplitude) (a) and the phase ϕ (mode velocity) (b) versus step-by-step increase in liquid conductivity σ (vertical arrows). Plate: 128°YX-LiNbO3, h = 350 μm, λ = 200 μm. Mode: f = 34.9 MHz, propagation direction Θ = 30° off X axis. Reference liquid: distilled water (200 mg, t < 0 s). Test liquid: distilled water (200 mg) with drops of 0.9% NaCl + H2O, 40 mg each (vertical arrows, 0 < t < 1200 s).
Figure 9. Time variations in the insertion loss S12 (mode amplitude) (a) and the phase ϕ (mode velocity) (b) versus step-by-step increase in liquid conductivity σ (vertical arrows). Plate: 128°YX-LiNbO3, h = 350 μm, λ = 200 μm. Mode: f = 34.9 MHz, propagation direction Θ = 30° off X axis. Reference liquid: distilled water (200 mg, t < 0 s). Test liquid: distilled water (200 mg) with drops of 0.9% NaCl + H2O, 40 mg each (vertical arrows, 0 < t < 1200 s).
Sensors 22 07231 g009
Figure 10. Calibration curves for one of the best modes towards liquid viscosity (a,b) and liquid conductivity (c,d) measured at 20 °C. Points are experimental data, solid lines (e,f) are numerical results calculated with formulas (2) and (3). Mode: f = 40.24 MHz, λ = 200 μm, v = 8.176 m/s, k2 = 110 × 10−4, surface displacements {u1, u2, u3,} = {1, 0.1, 0.5}. Plate: Y-LiNbO3, h/λ = 1.75, εpl = 34ε0. Propagation direction Θ = 60°. Liquids: water solutions of glycerin and NaCl, εlq0 = 79.3. Coated path L = 18 mm.
Figure 10. Calibration curves for one of the best modes towards liquid viscosity (a,b) and liquid conductivity (c,d) measured at 20 °C. Points are experimental data, solid lines (e,f) are numerical results calculated with formulas (2) and (3). Mode: f = 40.24 MHz, λ = 200 μm, v = 8.176 m/s, k2 = 110 × 10−4, surface displacements {u1, u2, u3,} = {1, 0.1, 0.5}. Plate: Y-LiNbO3, h/λ = 1.75, εpl = 34ε0. Propagation direction Θ = 60°. Liquids: water solutions of glycerin and NaCl, εlq0 = 79.3. Coated path L = 18 mm.
Sensors 22 07231 g010
Table 1. The characteristics of the delay lines shown in Figure 1.
Table 1. The characteristics of the delay lines shown in Figure 1.
IDT
Period,
μm
IDT
Aperture,
μm
Number of
Electrode Pairs
in IDT
Face-to-Face Distance between Input and Output IDTs, μmDistance L
Coated with Liquid, μm
Plate Thickness h/λ
20049002024,00018,0001.75 or 2.5
Table 2. Phase velocities vn, power flow angles Ψn, and own frequencies fn for the first five modes in 128°Y-LiNbO3 plate with normalized thickness h/λ = 1.75 (Θ = 30°).
Table 2. Phase velocities vn, power flow angles Ψn, and own frequencies fn for the first five modes in 128°Y-LiNbO3 plate with normalized thickness h/λ = 1.75 (Θ = 30°).
nvn, m/sΨn, degr.fn, MHz
03721.9−13.1°18.6
13804.95−8.9°19
24013.36−0.1°20
34272.112.7°21.4
44678.26−6.8°23.4
Table 3. Acoustic plate modes with enhanced sensitivity towards liquid conductivity measured for different plates and propagation directions.
Table 3. Acoustic plate modes with enhanced sensitivity towards liquid conductivity measured for different plates and propagation directions.
PlateΘ, degr.hf, MHzΔS1, dBΔS4/4, dBΔϕ1, degr.Δϕ4/4, degr.
128°Y-LNO90°2.536.55110.660°9.5°
128°Y-LNO1.7531.1870.864°8.5°
128°Y-LNO1.7543.460.2344°5.5°
128°Y-LNO30°1.7534.960.6559°6.5°
41°Y-LNO30°1.7535.925100.867°8.8°
41°Y-LNO90°1.7533.62151.5150°15°
Y-LNO30°1.7540.4290.840°21.3°
Y-LNO60°1.7540.24182102°32°
Y-LNO90°1.7538.724.60.637°6.1°
Y-LNO90°1.7552.48231.4282°14°
Y-LNO90°1.7558.46242.1300°10°
36°Y-LTaO2.562.21140.832°
36°Y-LTaO60°2.531.198.70.7582°7.5°
36°Y-LTaO90°2.52910176°13°
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Smirnov, A.; Anisimkin, V.; Voronova, N.; Shamsutdinova, E.; Li, P.; Ezzin, H.; Qian, Z.; Ma, T.; Kuznetsova, I. Multimode Design and Piezoelectric Substrate Anisotropy Use to Improve Performance of Acoustic Liquid Sensors. Sensors 2022, 22, 7231. https://doi.org/10.3390/s22197231

AMA Style

Smirnov A, Anisimkin V, Voronova N, Shamsutdinova E, Li P, Ezzin H, Qian Z, Ma T, Kuznetsova I. Multimode Design and Piezoelectric Substrate Anisotropy Use to Improve Performance of Acoustic Liquid Sensors. Sensors. 2022; 22(19):7231. https://doi.org/10.3390/s22197231

Chicago/Turabian Style

Smirnov, Andrey, Vladimir Anisimkin, Natalia Voronova, Elizaveta Shamsutdinova, Peng Li, Hamdi Ezzin, Zhenghua Qian, Tingfeng Ma, and Iren Kuznetsova. 2022. "Multimode Design and Piezoelectric Substrate Anisotropy Use to Improve Performance of Acoustic Liquid Sensors" Sensors 22, no. 19: 7231. https://doi.org/10.3390/s22197231

APA Style

Smirnov, A., Anisimkin, V., Voronova, N., Shamsutdinova, E., Li, P., Ezzin, H., Qian, Z., Ma, T., & Kuznetsova, I. (2022). Multimode Design and Piezoelectric Substrate Anisotropy Use to Improve Performance of Acoustic Liquid Sensors. Sensors, 22(19), 7231. https://doi.org/10.3390/s22197231

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop