Next Article in Journal
Smart Homes as Enablers for Depression Pre-Diagnosis Using PHQ-9 on HMI through Fuzzy Logic Decision System
Next Article in Special Issue
Comparing of Frequency Shift and Impedance Analysis Method Based on QCM Sensor for Measuring the Blood Viscosity
Previous Article in Journal
A Novel Optimized V-VLC Receiver Sensor Design Using μGA in Automotive Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

FEM Analysis of Various Multilayer Structures for CMOS Compatible Wearable Acousto-Optic Devices

1
Department of Electrical and Electronics Engineering, Universiti Teknologi PETRONAS, Seri Iskandar 32610, Perak, Malaysia
2
Faculty of Technical Sciences, University of Novi Sad, Trg Dositeja Obradovica 6, 21000 Novi Sad, Serbia
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(23), 7863; https://doi.org/10.3390/s21237863
Submission received: 16 October 2021 / Revised: 20 November 2021 / Accepted: 22 November 2021 / Published: 26 November 2021

Abstract

:
Lately, wearable applications featuring photonic on-chip sensors are on the rise. Among many ways of controlling and/or modulating, the acousto-optic technique is seen to be a popular technique. This paper undertakes the study of different multilayer structures that can be fabricated for realizing an acousto-optic device, the objective being to obtain a high acousto-optic figure of merit (AOFM). By varying the thicknesses of the layers of these materials, several properties are discussed. The study shows that the multilayer thin film structure-based devices can give a high value of electromechanical coupling coefficient (k2) and a high AOFM as compared to the bulk piezoelectric/optical materials. The study is conducted to find the optimal normalised thickness of the multilayer structures with a material possessing the best optical and piezoelectric properties for fabricating acousto-optic devices. Based on simulations and studies of SAW propagation characteristics such as the electromechanical coupling coefficient (k2) and phase velocity (v), the acousto-optic figure of merit is calculated. The maximum value of the acousto-optic figure of merit achieved is higher than the AOFM of all the individual materials used in these layer structures. The suggested SAW device has potential application in wearable and small footprint acousto-optic devices and gives better results than those made with bulk piezoelectric materials.

1. Introduction

In the last few decades, photonic on-chip sensors have gained significant attention amongst the various attractive photonic integrated circuits (PICs) applications [1,2,3,4]. Moreover, due to the advancement in optic structures (gratings, waveguides, etc.), photonics-based sensors have also gained considerable interest in several other fields like telecommunication and bio-medical instrumentation [5]. In the area of integrated optics, surface acousto-optic (AO) devices such as tunable optical filters, modulators, and optical switches are extensively used [6,7,8]. As the demand for AO devices increases, the interest in using AO material with high sensitivity towards acoustic or optic stimulants also increases.
The application of acousto-optic tunable filters (AOTFs) in spectroscopy has been well recognised in the past several years. Spectrometers based on AOTF technology have several appealing design characteristics. Among these advantages, the most prominent ones are high spectral resolution, less rf driving power, high durability, compact size, maintenance-free functioning, lightweight and quick data gathering [9]. Another important feature of AOTFs is, that they can simultaneously and individually filter a wide range of optical wavelengths by operating the acoustic transducer at various corresponding frequencies. They are thus suited for use in spectrometer instruments as the tuning element. A schematic of a typical miniaturized AOTF spectrometer is shown in Figure 1. In the past few years, efforts have been made to incorporate spectrometers in wearable sensors. AOTF-based spectrometers possess the potential to be incorporated into such sensors [10]. However, to design such a device the features and design parameters (acousto-optic) of an AOTF must be taken into consideration. As they play a vital role in the miniaturization of the device [11].
Previously, acousto-optic parameters were primarily investigated in the domain of crystallography; however, the focus has now switched to device optimisation [12]. Most of the photonics-based sensors work on the principle of acousto-optic effect in which there is a surface acoustic wave (SAW) device that acts as piezoelectric and optic devices. The most attractive features of SAW devices are that they consume less power and have a fast response time. So far, typical piezoelectric materials utilised in acousto-optic devices are langasite (La3Ga5SiO14), lead molybdate (PbMoO4), quartz (SiO2), tellurium dioxide (TeO2), lithium tantalite (LiTaO3),) and lithium niobate (LiNbO3) [4]. These bulk piezoelectric materials are used in most acoustic devices, but they are generally expensive and difficult to integrate with electronics, and some are even toxic in nature. Moreover, they are usually delicate and fragile if polished into extremely thin structures for high-frequency applications. In addition to this, using these large bulk materials is impractical in terms of cost, size, acousto-optic figure of merit, and integration [13]. The optimum thickness of the substrate required for SAW applications is normally around 10 λ or more. This thickness is chosen so that all the SAW travel on the top surface (piezoelectric material) and eventually dies out when it travels deeper into the substrate. However, for acousto-optic purposes where SAW interacts with light over an optical waveguide, the waveguide thickness can be of the order of 1 or more wavelengths of the infrared maximum range—nearly 5 microns or less. Hence, the multilayer substrate thickness is expected to be in excess of 10 s of the wavelength while the optical waveguide thickness can be about 1 wavelength or so. Attempts to replace these bulk materials with thin films have produced promising results [14]. Here, we emphasise the need to evaluate and optimise these thin film structures for achieving a better acousto-optic figure of merit.

2. Literature Review

The use of thin film-based optical devices [15,16,17,18] has increased recently due to the advancement of thin-film deposition technology that can produce thin films with high uniformity, homogeneity, etc. These devices typically use single-layer thin films and can achieve remarkable performance. However, the performance of these devices can be further improved by considering alternative thin films that have been made available recently [19]. Using these thin films in a multilayer structure can significantly enhance the device’s performance. For instance, in designing optical filters or SAW devices, multilayer structures improve phase velocity, temperature stability, and mechanical strength [20,21,22]. In addition to this, the use of multilayer thin films is also preferable for broad tunability compared to the doped semiconductors, as these multilayer thin films give the desired response for various targeted applications. A typical example of it is an acousto-optic tunable filter (AOTF) [23,24,25,26]. The performance of devices utilising multilayered structures is dependent on various properties of a material medium such as acoustic, optic, and acousto-optic properties using both the sound and light wave interaction. For such devices to significantly improve the productivity of acousto-optic interaction, strong SAW interaction and optical waveguiding with low loss are essential. Thus, the determination of a reasonable layered structure is necessary. The layered structure utilised in most of the acousto-optic devices ought to have the characteristics such as high electromechanical coupling coefficient (k2) and phase velocity (v) [27]. These parameters are then formulated into the acousto-optic figure of merit that can measure the suitability of such materials for AO devices.
The acousto-optic figure of merit measures a material’s capacity to modulate diffraction intensity. For the effective realisation of such devices, it is better to make the blend of layer structure by utilising the materials having said properties [7,19]. Zinc oxide (ZnO) thin films exhibit good piezoelectric characteristics with high k2, sensitivity, reliability and attractive optical properties. [28,29]. On account of promising optical, electro-optic, piezoelectric, photoelastic, photorefractive and elastic properties [30], ferroelectric lithium niobate (LiNbO3) is also broadly utilised in guided-wave optics. Notwithstanding that, another component tends to be developed on various substrates, incorporating silicon (Si). Additionally, aluminium nitride (AlN)-based layered SAW has attracted significant interest [25,26,27,28]. It is mainly popular among other piezoelectric materials because of its high acoustic velocity (5760 ms−1). Another important feature of AlN thin film is its ease of integration into complementary metal-oxide-semiconductor (CMOS) platforms when contrasted with other piezoelectric materials like LiNbO3, ZnO, LiTaO3 and quartz, as these materials create contamination in integrated circuit (IC) processes [31,32]. Within the last several years, there has been a tremendous rate of increase in interest in CMOS compatible surface acoustic wave (SAW) devices. The main reason to prefer CMOS compatible SAW devices over other devices is to eliminate the parasitic capacitance effect (mainly occurs when discrete SAW devices are integrated circuit connections are bonded). Furthermore, another reason is to minimise the footprints by merging the SAW device and the integrated circuit on a single chip.
In this manuscript, we investigated three multilayered structures A (ZnO/SiO2/Si), B (LiNbO3/SiO2/Si), and C (AlN/SiO2/Si). There are three piezoelectric materials used in this study (i.e., AlN, ZnO, and LiNbO3). AlN and ZnO are in c-axis orientation with texture 11 20 ̲ and for LiNbO3, 128° Y-cut Z-propagating is used. The top layers, i.e., LiNbO3, AlN, ZnO, are commonly used as bulk materials in AO devices with good piezoelectric and optic properties [20,21,22,23]. By comparing the propagation properties of SAW within said structures, the acousto-optic figure of merit of these three multilayer structures is then calculated.

3. Materials and Methods

This work uses the COMSOL Multiphysics software to perform 2D finite element method (FEM) simulations (modal and harmonic analysis) on the layered SAW device. Researchers used, validated, and reported on the results of the FEM study on different piezoelectric materials [33,34,35,36]. The complexity of modelling the SAW modes with longitudinal and shear vertical particle displacements can be significantly reduced using 2D structures. Furthermore, the use of an infinite set of IDTs helps in simulating a standing wave pattern that results when a resonator is employed. So, a unit cell with periodic boundary conditions on either side of the cell becomes a useful tool to simulate a SAW resonator with an infinite set of IDTs [32,37,38,39,40].
To determine the acousto-optic figure of merit, various SAW propagation characteristics, such as phase velocity and k2, are calculated and used. The electromechanical coupling coefficient (k2) is a numerical measure of the conversion efficiency between electrical and acoustic energy in piezoelectric materials. In this study, all three structures work on the principle of the acousto-optic effect in which there is a surface acoustic wave (SAW) device that acts as a piezoelectric and optic device. In designing an efficient SAW-based acousto-optic device, one of the concerns is the driving power. The input driving power is directly related to the electromechanical coupling coefficient. Higher the k2, lesser is the input power needed to convert the electrical energy to acoustic energy, which, as a result, enhances the device efficiency. Hence, the electromechanical coupling coefficient is the parameter that enables us to know the efficiency of the acousto-optic device. Additionally, the electromechanical coupling coefficient, which is a parameter that also depends on the piezoelectric film thickness, is also related to acoustic bandwidth. Higher bandwidth allows for a wider range of modulation for RF carrier, which, in turn, means a wider range of tunability of acousto-optic filters. Thus, this study analyses the FEM 2D models A, B, and C (Figure 2). Table 1 summarises the dimensions of the SAW device used in the simulation, while Figure 2 details the geometry of the 2D single-cell multilayered SAW device.
Hence, for simulating a SAW device, a single cell measuring one λ is taken. A total of 10 λ is taken as the depth of a unit cell because SAW is generated on the top surface and dies out in the lower boundary as it travels down into the substrate. The material used for IDTs in this device is aluminium (Al). The properties of Al are obtained from the inbuilt COMSOL library, which is the density of 2700 kilogrammes per cubic metre, 70 GPa Young’s modulus, and a Poisson’s ratio of 0.33. The ratio of metallisation is taken as 50% (2a/λ), while 2.5% is the electrode’s relative thickness (b/λ). Here, the width and height of electrodes are a (λ/4) and b (100 nm), respectively. To improve the accuracy of the results, a very thin triangular mesh is used for the simulations. Both the electrodes, from the top and sides, are free in terms of mechanical boundaries. The boundaries of the left electrode are grounded, and for the right one, it is set to a floating potential with no charge accumulation on the surface. The boundary conditions utilised in the simulations are shown in Table 2, and the material constants used in the simulations are listed in Table 3. All the material constants mentioned in Table 3 are taken from reference [41,42,43,44].
In general, the equation for determining the acoustic wave velocity is as follows:
v = f λ
The central frequency is denoted by f, and the wavelength of the acoustic wave is represented by λ. To find the velocity in eigenmode, the following equation is used:
v = f r e s + f a n t i p
Here, both f r e s and f a n t i frequencies denote the resonant and anti-resonant frequencies, respectively, and p represents the pitch (electrode spacing + electrode width). In our study, the pitch value is 2 µm (λ/2), and the resonant and anti-resonant frequencies are eigenfrequencies we obtainthrough simulation analysis. The equation used to find k is as follows:
k 2 = 2   v v m v
Where vm is the phase velocity with the metal short, and v is the free surface velocity without it.
ZnO, LiNbO3, and AlN thin films act as optic and piezoelectric layers in the proposed layer structure as shown in Figure 2a–c, respectively, with IDTs on top. The interaction of light and acoustic waves occurs on the top layer in all the structures. In all of these structures, a SiO2 layer is deposited on the silicon (Si) substrate. It is incorporatedbecause the index of refraction of Si is greater than the indexes of refraction of ZnO, LiNbO3 and AlN, whereas the index of refraction of SiO2 is less than the other materials deposited on top of it. Therefore, the waveguide action cannot be achieved without it [45]. Variation of the electric field in the IDTs, travelling in the top layer of all the structures, produces SAW and penetrates SiO2 from the top layer [46]. The light waves fall on the top layer (acting as the optic layer). It creates the acousto-optic effect after interaction with SAW, which is the basic concept of acousto-optic devices [47]. ZnO, LiNbO3, and AlN are used in these structures because of their strong piezoelectric and optical characteristics, as well as the ease with which they are deposited [48,49,50,51,52,53,54,55,56]. A material’s k2 is dependent on characteristics such as dielectric, piezoelectric, and elastic properties. Therefore, the materials are selected based on these characteristics.
Furthermore, while choosing a material for the acousto-optic devices, it should have a suitable acousto-optic figure of merit, which is determined by the photoelastic constant (p), density (ρ), refractive index (n), and acoustic phase velocity (v). The acousto-optic figure of merit measures the suitability of a material to modulate the diffraction intensity. The acousto-optic figure of merit (M) can be expressed mathematically as:
M = n 6 p 2 ρ v 3  
All the parameters in Equation (4) should be considered when selecting the top layer material. The refractive index and acoustic velocity are the most essential. The slower the acoustic and optical waves in the material, the more interaction is possible [57]. Some materials offer good optical properties like less losses and high transmittance, but the acousto-optic figure of merit becomes low because of having high acoustic velocity. There are other figures of merit related to acousto-optic devices; however, the acousto-optic figure of merit referred to in Equation (4) is used primarily for gauging the power efficiency of acousto-optic materials. It is preferable to use a combination of materials to achieve a high acousto-optic figure of merit. Therefore, to decrease the overall structure’s velocity and increase the acousto-optic figure of merit, it is better to incorporate the structure with a material possessing low velocity. The velocity of the layer structure is the combined velocity of all the materials. With different combinations of these materials, the acousto-optic figure of merit is calculated.

4. Results and Discussion

To validate our FEM simulation method used in the structures under study, we first began by simulating the AlN/Diamond (slow–fast) structure, which was previously described [58]. In that study, the propagation characteristics of SAW were theoretically computed using PC acoustic wave software developed by McGill University. SAW propagation in multilayer structures is calculated with this software using the transfer matrix approach, which is implemented in this software. We obtained their numerical data by employing a semi-automated open-source plot digitiser application, i.e., WebPlotDigitizer v. 3.12, which has been utilised in several published papers. To conduct our investigation, we used the unit cell shown in Figure 2 (but with two layers, namely AlN/Diamond) and the boundary conditions shown in Table 2 to simulate the AlN/Diamond structure. The same material constants as those given were utilised [58]. Equations (1)–(3) are used to compute the SAW propagation characteristics. Figure 3 depicts the simulated results, which are somewhat similar to those presented in [58]. The slight inaccuracy can be attributed to the fact that the procedure employed in this case is different from the work presented in [58]. The approach utilised in [58] was based on a transfer matrix method that made assumptions about numerous thin layers’ behaviour to approximate the behaviour of interfacial layers. In contrast, the method employed in this study is based on a finite element method (FEM).
A mesh convergence investigation was also carried out to justify the appropriateness of the chosen mesh density. The mesh density was increased until the SAW velocity values became constant. In this case, the most triangular mesh elements possible were selected (i.e., 120,503). However, increasing the number of mesh elements also increases the processing time. Therefore, to achieve accurate results, it has to be performed at the expense of processing time.
In this work, we calculated phase velocity and electromechanical coupling coefficient of structures A (ZnO/SiO2/Si), B (LiNbO3/SiO2/Si), and C (AlN/SiO2/Si) on different thicknesses of SiO2 (0.4 µm, 1 µm, 2 µm, and 3 µm) with changing thicknesses of ZnO, LiNbO3 and AlN (0.4 µm, 0.8 µm, 1.6 µm, 2.4 µm, 3.2 µm and 4 µm) using COMSOL. The values of resonant and anti-resonant frequencies are obtained using COMSOL. Afterward, using Equations (2) and (3) the SAW phase velocity and electromechanical coupling coefficient are calculated, respectively. The polarised mode studied in this work is the Rayleigh mode which occurs in the sagittal plane. The displacement profile for structures A (ZnO/SiO2/Si), B (LiNbO3/SiO2/Si), and C (AlN/SiO2/Si) is shown in Figure 4, Figure 5 and Figure 6, respectively.
Figure 4 depicts the displacement profile of structure A (ZnO/SiO2/Si). Figure 4a–c shows the displacement profile with the normalised thickness of ZnO as 0.01 λ, 0.1 λ and 2 λ, respectively. It can be vividly seen that by increasing the thickness of ZnO that the surface acoustic waves start to confine in the ZnO layer. It can be seen from Figure 4b that the whole wave is confined in ZnO and SiO2 layer and none of it is going in Silicon substrate.
Figure 5 shows the displacement profile of structure B (LiNbO3/SiO2/Si). The shapes of the mode in the displacement profile help to identify the Rayleigh wave modes. Figure 5a–c shows the results of Eigen frequency analysis for anti-resonance modes at t LiNbO3/λ = 0.01, 0.1, and 2 respectively. Moreover, it can be seen clearly from Figure 5a–c that by increasing the thickness of the LiNbO3 layer, surface acoustic wave no longer travels in Silicon substrate. Additionally, if it keeps on increasing, it gets entirely confined only in the LiNbO3 layer.
Similarly, the displacement profile of structure C (AlN/SiO2/Si) is summarised in Figure 6. The displacement profile’s mode shapes help in identifying Rayleigh wave modes. The results in Figure 6a–c were obtained using eigenfrequency analysis in anti-resonance modes for tAlN/λ = 0.01, 0.1 and 2. It can be seen that increasing the thickness of AlN thin film value of eigenfrequency increases and the surface acoustic wave no longer passes in Silicon substrate and only travels in AlN and SiO2 thin film, as shown in Figure 6b. It continues to do so as the thickness of the AlN layer increases until the entire acoustic wave is contained within the AlN substrate (as seen in Figure 6c). On this point, it reaches a velocity of 5507 m/s at tAlN/λ = 2, which is close to AlN’s theoretical acoustic velocity (5600 m/s) [44].
The simulation results for SAW propagation characteristics of the multilayer structures A, B and C are displayed in Figure 7 and Figure 8. In this study, the normalised thickness chosen for SiO2 is 0.1 λ, 0.25 λ, 0.5 λ, 0.75 λ and 1 λ. Whereas, 0.1 λ, 0.2 λ, 0.4 λ, 0.6 λ, 0.8 λ and 1.0 λ normalised thickness have been chosen in the case of AlN, ZnO and LiNbO3.
It is clear from Figure 7 that the overall phase velocity decreases as the thickness of SiO2 increases. All three piezoelectric materials used in structures A, B and C possess high acoustic velocities as compared to SiO2. The decrease in the overall behaviour of phase velocity in the structures A, B and C, due to the increase in the normalised thickness of SiO2 can be explained using the fact that there is a significant difference between the velocities of SiO2 (3750 m/s) with respect to the velocities of AlN (5760 m/s), ZnO (6133 m/s) and LiNbO3 (6487 m/s). However, if the change in the normalised thickness of AlN, ZnO and LiNbO3 is observed, the case is different. In the case of AlN, when the thickness of AlN increases, the phase velocity also increases compared to the case when the thickness of SiO2 was increased. It happens because when the normalised thicknesses of AlN increases, the entire surface acoustic wave travels only in AlN layers, and none of them reaches the substrate. All surface acoustic waves become confined in the piezoelectric layer.
On the contrary, when the thickness of ZnO or LiNbO3 increases in structures A and B, as shown in Figure 7a,b, a declined trend is observed. However, by increasing the thickness by more than 0.2 λ, there is no significant change in the behaviour of the phase velocity observed. It happens because the surface acoustic wave becomes confined in the top piezoelectric layer only.
Figure 8 shows the behaviour of the electromechanical coupling coefficient with respect to normalised thicknesses of SiO2 and the top layers in the structures A, B and C (i.e., ZnO, LiNbO3 and AlN, respectively). The overall electromechanical coupling coefficient behaviour showed a decreasing trend with increased normalised thicknesses of SiO2. Figure 8a shows the change in the electromechanical coupling coefficient as ZnO’s thickness increases. It is observed that by increasing the normalised thickness of ZnO, the electromechanical coupling coefficient decreased sharply. However, after the normalised value of 0.4, no significant change in its value is observed. At the smaller values of the normalised thickness of ZnO (tZnO/λ = 0.1 and 0.2), surface acoustic waves are also travelling in SiO2. Whereas, when the normalised thickness is increased (greater than tZnO/λ = 0.2), the surface acoustic waves become confined in ZnO, due to which the electromechanical coupling coefficient becomes almost constant. Moreover, in all three structures, the highest value of k2 is achieved at 0.4 µm thickness of SiO2. Additionally, among all three structures, the highest value of electromechanical coupling coefficient is attained in structure AlN/SiO2/Si (at 0.4 µm thickness of SiO2 and AlN), which is 7.15.
Table 4 illustrates the acousto-optic figure of merit, phase velocity, and electromechanical coupling coefficient of structures A, B, and C. It is evident from Table 4 that among all the structures, structure C (AlN/SiO2/Si) exhibits the highest value of electromechanical coupling coefficient (7.15 at the normalised thickness of 0.1 λ of SiO2 and AlN layer) with the acousto-optic figure of merit 7.04 × 10−14 s3/kg. Moreover, it is to be noted that the acousto-optic figure of merit achieved in structure AlN/SiO2/Si with the thicknesses as mentioned earlier is higher than the acousto-optic figure of merit of all the individual materials used in these layer structures (i.e., ZnO, LiNbO3, AlN and SiO2) [29,49].
In light of the results discussed above, it can be concluded that to fabricate an acousto-optic device, structure C (AlN/SiO2/Si) would give promising results. However, there is a trade-off between normalised thickness, the acousto-optic figure of merit, acoustic velocity, and electromechanical coupling coefficient. Changing the value of normalised thicknesses, the above said properties could be tailored according to the requirements. Furthermore, AlN is well-suited for use in CMOS technologies since it can be readily deposited by physical vapour deposition at low temperatures, making it a cost-effective option [59,60].

5. Conclusions

The study undertaken was targeted towards studying various CMOS compatible multilayer structures that would produce better acousto-optic figures of merit for an AOTF application. SAW propagation characteristics were investigated in three multilayer structures: ZnO/SiO2/Si, LiNbO3/SiO2/Si, and AlN/SiO2/Si. It is found that the AlN/SiO2/Si multilayer structure gave the highest electromechanical coupling coefficient, i.e., 7.15 at the normalised thickness of 0.1 λ (for SiO2 and AlN layer). It also has the highest acousto-optic figure of merit, i.e., 7.04 × 10−14 s3/kg (which is also higher than the AOFM of individual bulk materials found in literature and mentioned in Table 3). Hence, it is concluded that by using an AlN/SiO2/Si multilayer structure, a high acousto-optic figure of merit can be achieved, which is not possible by using any choice of single bulk material. Moreover, it makes it possible to realise lab-on-chip applications using structure C (AlN/SiO2/Si), which is a CMOS compatible SAW structure.

Author Contributions

The work given here was primarily conducted by M.H. under the guidance of V.J., S.Q. and M.R.A.; M.H. contributed to the device design and multiphysics simulation. S.Q. and M.Z.A. gave helpful recommendations to the research.; conceptualisation, V.J. and M.H.; methodology, M.H., S.Q. and V.J.; software, M.H., S.Q. and M.Z.A.; validation, S.Q., V.J., M.R.A. and G.S.; formal analysis, V.J. and M.H.; investigation, M.H.; writing—original draft preparation, M.H.; writing—review and editing, M.R.A., V.J. and G.S.; supervision, M.R.A. and V.J.; All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the European Union’s Horizon 2020 research and innovation programme under grant agreement No. 854194.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This study was financially supported by the Postgraduate Studies programme of the Universiti Teknologi PETRONAS, Malaysia, and also the funding received from the European Union’s Horizon 2020 research and innovation programme under grant agreement No. 854194. The authors would also like to further acknowledge Annas-ul-Hassan Qureshi, Muhammad Wahyu Nugraha, Muhammad Hamza Azam and Danish Mahmood Khan for their reviewing contribution.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Washburn, A.L.; Bailey, R.C. Photonics-on-a-chip: Recent advances in integrated waveguides as enabling detection elements for real-world, lab-on-a-chip biosensing applications. Analyst 2011, 136, 227–236. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Hänsel, A.; Heck, M.J. Opportunities for photonic integrated circuits in optical gas sensors. J. Phys. Photonics 2020, 2, 012002. [Google Scholar] [CrossRef]
  3. Zinoviev, K.; Carrascosa, L.G.; Sánchez del Río, J.; Sepulveda, B.; Domínguez, C.; Lechuga, L.M. Silicon photonic biosensors for lab-on-a-chip applications. Adv. Opt. Technol. 2008, 2008, 1–6. [Google Scholar] [CrossRef]
  4. Zhu, H.; Goi, K.; Ishikura, N.; Omichi, K. Silicon Photonics Based System-On-Chip Gas Sensor. In Proceedings of the Optical Fiber Sensors, Zurich, Switzerland, 2–5 July 2018; Optical Society of America: Washington, DC, USA, 2018; p. ThE46. [Google Scholar]
  5. Micó, G.; Gargallo, B.; Pastor, D.; Muñoz, P. Integrated optic sensing spectrometer: Concept and design. Sensors 2019, 19, 1018. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Munk, D.; Katzman, M.; Hen, M.; Priel, M.; Feldberg, M.; Sharabani, T.; Levy, S.; Bergman, A.; Zadok, A. Surface acoustic wave photonic devices in silicon on insulator. Nat. Commun. 2019, 10, 4214. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Poveda, A.C.; Bühler, D.D.; Sáez, A.C.; Santos, P.V.; Lima, M.M.d., Jr. Optical waveguide devices modulated by surface acoustic waves. arXiv 2018, arXiv:1811.03051. [Google Scholar] [CrossRef] [Green Version]
  8. Van Der Slot, P.J.; Porcel, M.A.; Boller, K.-J. Surface acoustic waves for acousto-optic modulation in buried silicon nitride waveguides. Opt. Express 2019, 27, 1433–1452. [Google Scholar] [CrossRef] [PubMed]
  9. Li, J.; Gui, Y.; Xu, R.; Zhang, Z.; Liu, W.; Lv, G.; Wang, M.; Li, C.; He, Z. Applications of AOTF Spectrometers in In Situ Lunar Measurements. Materials 2021, 14, 3454. [Google Scholar] [CrossRef]
  10. Yang, Z.; Albrow-Owen, T.; Cai, W.; Hasan, T. Miniaturization of optical spectrometers. Science 2021, 371. [Google Scholar] [CrossRef] [PubMed]
  11. Govindaraj, R.; Ganesh, G. Development and integration of an AOTF based NIR spectrophotometer. J. Opt. 2018, 47, 132–142. [Google Scholar] [CrossRef]
  12. Fernandes, M.E. Acousto-Optic Effect and Its Use in Signal Processing. Available online: https://www.semanticscholar.org/paper/Acousto-Optic-Effect-and-Its-use-in-Signal-Fernandes/9ce8b3947de1e06d1f7e83066244098cf10145b6 (accessed on 21 November 2021).
  13. Honardoost, A. Thin-Film Lithium Niobate Integrated Photonics on Silicon for Electro-and Nonlinear-Optic Applications. Ph.D. Thesis, University of Central Florida, Orlando, FL, USA, May 2020. [Google Scholar]
  14. Muralt, P.; Polcawich, R.G.; Trolier-McKinstry, S. Piezoelectric thin films for sensors, actuators, and energy harvesting. MRS Bull. 2009, 34, 658–664. [Google Scholar] [CrossRef] [Green Version]
  15. Meillaud, F.; Boccard, M.; Bugnon, G.; Despeisse, M.; Hänni, S.; Haug, F.-J.; Persoz, J.; Schüttauf, J.-W.; Stuckelberger, M.; Ballif, C. Recent advances and remaining challenges in thin-film silicon photovoltaic technology. Mater. Today 2015, 18, 378–384. [Google Scholar] [CrossRef]
  16. Dharmadasa, I. Advances in Thin-Film Solar Cells; Jenny Stanford Publishing: Singapore, 2018. [Google Scholar]
  17. Melville, O.A.; Lessard, B.H.; Bender, T.P. Phthalocyanine-based organic thin-film transistors: A review of recent advances. ACS Appl. Mater. Interfaces 2015, 7, 13105–13118. [Google Scholar] [CrossRef]
  18. Fortunato, E.; Barquinha, P.; Martins, R. Oxide semiconductor thin-film transistors: A review of recent advances. Adv. Mater. 2012, 24, 2945–2986. [Google Scholar] [CrossRef]
  19. Fu, Y.Q.; Luo, J.; Nguyen, N.-T.; Walton, A.; Flewitt, A.J.; Zu, X.-T.; Li, Y.; McHale, G.; Matthews, A.; Iborra, E. Advances in piezoelectric thin films for acoustic biosensors, acoustofluidics and lab-on-chip applications. Prog. Mater. Sci. 2017, 89, 31–91. [Google Scholar] [CrossRef] [Green Version]
  20. Juneau, J.-M. The Simulation, Design, and Fabrication of Optical Filters; Rose Hulman Institute of Technology: Terre Haute, IN, USA, 2017. [Google Scholar]
  21. Keshmiri, S.; Mirsalehi, M. Multilayer Thin-Film Optical Filters: Design, Fabrication, and Applications. In Physics and Technology of Thin Films: IWTF 2003; World Scientific: Singapore, 2004; pp. 306–317. [Google Scholar]
  22. Kalantar-Zadeh, K.; Powell, D.A.; Wlodarski, W.; Ippolito, S.; Galatsis, K. Comparison of layered based SAW sensors. Sens. Actuators B Chem. 2003, 91, 303–308. [Google Scholar] [CrossRef]
  23. Chang, I. Tunable acousto-optic filters: An overview. In Proceedings of the Acousto-Optics: Device Development/Instrumentation/Applications, San Diego, CA, USA, 26–27 August 1976; International Society for Optics and Photonics: Bellingham, DC, USA, 1976; pp. 12–22. [Google Scholar]
  24. Mantsevich, S.; Korablev, O.; Kalinnikov, Y.K.; Ivanov, A.Y.; Kiselev, A. Wide-aperture TeO2 AOTF at low temperatures: Operation and survival. Ultrasonics 2015, 59, 50–58. [Google Scholar] [CrossRef] [PubMed]
  25. Georgiev, G.; Glenar, D.A.; Hillman, J.J. Spectral characterization of acousto-optic filters used in imaging spectroscopy. Appl. Opt. 2002, 41, 209–217. [Google Scholar] [CrossRef]
  26. Mukhopadhyay, T.; Mahata, A.; Adhikari, S.; Zaeem, M.A. Effective mechanical properties of multilayer nano-heterostructures. Sci. Rep. 2017, 7, 15818. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Pfeiffer, J.B.; Wagner, K.H. Acousto-optic figure of merit search. Phys. Procedia 2015, 70, 762–765. [Google Scholar] [CrossRef] [Green Version]
  28. Yu, C.; Yu, X.-X.; Zheng, D.-S.; Yin, H. Piezoelectric potential enhanced photocatalytic performance based on ZnO with different nanostructures. Nanotechnology 2021, 32, 135703. [Google Scholar] [CrossRef] [PubMed]
  29. Yang, A.; Qiu, Y.; Yang, D.; Lin, K.; Guo, S. Piezoelectric property comparison of two-dimensional ZnO nanostructures for energy harvesting devices. RSC Adv. 2021, 11, 3363–3370. [Google Scholar] [CrossRef]
  30. Weis, R.; Gaylord, T. Lithium niobate: Summary of physical properties and crystal structure. Appl. Phys. A 1985, 37, 191–203. [Google Scholar] [CrossRef]
  31. Kaletta, U.C.; Wipf, C.; Fraschke, M.; Wolansky, D.; Schubert, M.A.; Schroeder, T.; Wenger, C. AlN/SiO2/Si3N4/Si(100)-based CMOS compatible surface acoustic wave filter with− 12.8-dB minimum insertion loss. IEEE Trans. Electron Devices 2015, 62, 764–768. [Google Scholar] [CrossRef]
  32. Kaletta, U.C.; Wenger, C. FEM simulation of Rayleigh waves for CMOS compatible SAW devices based on AlN/SiO2/Si (1 0 0). Ultrasonics 2014, 54, 291–295. [Google Scholar] [CrossRef] [PubMed]
  33. Aslam, M.Z.; Jeoti, V.; Karuppanan, S.; Malik, A.F.; Iqbal, A. FEM analysis of sezawa mode SAW sensor for VOC based on CMOS compatible AlN/SiO2/Si multilayer structure. Sensors 2018, 18, 1687. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Zhang, Z.; Wen, Z.; Wang, C. Investigation of surface acoustic waves propagating in ZnO–SiO2–Si multilayer structure. Ultrasonics 2013, 53, 363–368. [Google Scholar] [CrossRef]
  35. Shu, L.; Peng, B.; Li, C.; Gong, D.; Yang, Z.; Liu, X.; Zhang, W. The characterization of surface acoustic wave devices based on AlN-metal structures. Sensors 2016, 16, 526. [Google Scholar] [CrossRef] [PubMed]
  36. Mohanan, A.A.; Islam, M.S.; Ali, S.H.M.; Parthiban, R.; Ramakrishnan, N. Investigation into mass loading sensitivity of sezawa wave mode-based surface acoustic wave sensors. Sensors 2013, 13, 2164–2175. [Google Scholar] [CrossRef]
  37. Maouhoub, S.; Aoura, Y.; Mir, A. FEM simulation of AlN thin layers on diamond substrates for high frequency SAW devices. Diam. Relat. Mater. 2016, 62, 7–13. [Google Scholar] [CrossRef]
  38. Hofer, M.; Finger, N.; Kovacs, G.; Schoberl, J.; Langer, U.; Lerch, R. Finite element simulation of bulk-and surface acoustic wave (SAW) interaction in SAW devices. In Proceedings of the 2002 IEEE Ultrasonics Symposium, Munich, Germany, 8–11 October 2002; pp. 53–56. [Google Scholar]
  39. Buchner, M.; Ruile, W.; Dietz, A.; Dill, R. FEM analysis of the reflection coefficient of SAWs in an infinite periodic array. In Proceedings of the IEEE 1991 Ultrasonics Symposium, Orlando, FL, USA, 8–11 December 1991; pp. 371–375. [Google Scholar]
  40. Yong, Y.-K. Analysis of periodic structures for BAW and SAW resonators. In Proceedings of the 2001 IEEE Ultrasonics Symposium. Proceedings. An International Symposium (Cat. No. 01CH37263), Atlanta, GA, USA, 7–10 October 2001; pp. 781–790. [Google Scholar]
  41. McIntosh, R.; Bhalla, A.S.; Guo, R. Finite element modeling of acousto-optic effect and optimization of the figure of merit. In Proceedings of the Photonic Fiber and Crystal Devices: Advances in Materials and Innovations in Device Applications VI, San Diego, CA, USA, 12–13 August 2021; p. 849703. [Google Scholar] [CrossRef]
  42. Davydov, S.Y. Evaluation of physical parameters for the group III nitrates: BN, AlN, GaN, and InN. Semiconductors 2002, 36, 41–44. [Google Scholar] [CrossRef]
  43. Aslam, M.Z.; Jeoti, V.; Karuppanan, S.; Malik, A.F. FEM Simulation Analysis of AlN/SiO 2/Si Multilayer Structure and Effect of IDT Configuration on SAW Propagation Modes and Characteristics. In Proceedings of the 2018 International Conference on Intelligent and Advanced System (ICIAS), Kuala Lumpur, Malaysia, 13–14 August 2018; pp. 1–4. [Google Scholar]
  44. Malik, A.F.; Jeoti, V.; Fawzy, M.; Iqbal, A.; Aslam, Z.; Pandian, M.S.; Marigo, E. Estimation of SAW velocity and coupling coefficient in multilayered piezo-substrates AlN/SiO2/Si. In Proceedings of the 2016 6th International Conference on Intelligent and Advanced Systems (ICIAS), Kuala Lumpur, Malaysia, 15–17 August 2016; pp. 1–5. [Google Scholar]
  45. Jain, S.; Mansingh, A. Thin film layered structure for acousto-optic devices. J. Phys. D Appl. Phys. 1992, 25, 1116. [Google Scholar] [CrossRef]
  46. Nayak, R.; Nayak, A.; Gupta, V.; Sreenivas, K. Optical interactions in ZnO-TeO/sub 2/bi-layer for AO device applications. In Proceedings of the IEEE Symposium on Ultrasonics 2003, Honolulu, HI, USA, 5–8 October 2003; pp. 1129–1132. [Google Scholar]
  47. Stedham, C.; Draper, M.; Ward, J.; Wachman, E.; Pannell, C. A novel acousto-optic tunable filter for use in hyperspectral imaging systems. In Proceedings of the Physics and Simulation of Optoelectronic Devices XVI, San Jose, CA, USA, 22 February 2008; p. 68891. [Google Scholar]
  48. Kitui, M.; Mwamburi, M.M.; Gaitho, F.; Maghanga, C.M. Optical Properties of TiO2 Based Multilayer Thin Films: Application to Optical Filters. Int. J. Thin Film Sci. Technol. 2015, 4, 1. [Google Scholar]
  49. Maouhoub, S.; Aoura, Y.; Mir, A. FEM simulation of Rayleigh waves for SAW devices based on ZnO/AlN/Si. Microelectron. Eng. 2015, 136, 22–25. [Google Scholar] [CrossRef]
  50. Mahmood, T.; Cao, C.; Khan, W.S.; Usman, Z.; Butt, F.K.; Hussain, S. Electronic, elastic, optical properties of rutile TiO2 under pressure: A DFT study. Phys. B Condens. Matter 2012, 407, 958–965. [Google Scholar] [CrossRef]
  51. Rathore, B.P.S.; Prakash, R.; Kaur, D. Effect of AlN layer on the resistive switching properties of TiO2 based ReRAM memory devices. Curr. Appl. Phys. 2018, 18, 102–106. [Google Scholar] [CrossRef]
  52. Yoshida, F.; Nagashima, K.; Tsubouchi, M.; Maruyama, M.; Ochi, Y. THz pulse generation using a contact grating device composed of TiO2/SiO2 thin films on LiNbO3 crystal. J. Appl. Phys. 2016, 120, 183103. [Google Scholar] [CrossRef]
  53. Hao, L.; Li, Y.; Zhu, J.; Wu, Z.; Long, F.; Liu, X.; Zhang, W. Microstructure and memory characteristics of ferroelectric LiNbO3/ZnO composite thin films on Pt/TiO2/SiO2/Si substrates. J. Alloy. Compd. 2014, 590, 205–209. [Google Scholar] [CrossRef]
  54. Shirazi, M.; Hosseinnejad, M.; Zendehnam, A.; Ghorannevis, Z.; Ghoranneviss, M. Deposition of ZnO multilayer on LiNbO3 single crystals by DC-magnetron sputtering. Appl. Surf. Sci. 2011, 257, 10233–10238. [Google Scholar] [CrossRef]
  55. Xiong, S.; Liu, X.; Zhou, J.; Liu, Y.; Shen, Y.; Yin, X.; Wu, J.; Tao, R.; Fu, Y.; Duan, H. Stability studies of ZnO and AlN thin film acoustic wave devices in acid and alkali harsh environments. RSC Adv. 2020, 10, 19178–19184. [Google Scholar] [CrossRef]
  56. Kischkat, J.; Peters, S.; Gruska, B.; Semtsiv, M.; Chashnikova, M.; Klinkmüller, M.; Fedosenko, O.; Machulik, S.; Aleksandrova, A.; Monastyrskyi, G. Mid-infrared optical properties of thin films of aluminum oxide, titanium dioxide, silicon dioxide, aluminum nitride, and silicon nitride. Appl. Opt. 2012, 51, 6789–6798. [Google Scholar] [CrossRef] [PubMed]
  57. Newnham, R.E. Properties of Materials: Anisotropy, Symmetry, Structure; Oxford University Press on Demand: Oxford, UK, 2005. [Google Scholar]
  58. Benetti, M.; Cannata, D.; Di Pictrantonio, F.; Verona, E. Growth of AlN piezoelectric film on diamond for high-frequency surface acoustic wave devices. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2005, 52, 1806–1811. [Google Scholar] [CrossRef] [PubMed]
  59. Takagaki, Y.; Santos, P.; Wiebicke, E.; Brandt, O.; Schönherr, H.-P.; Ploog, K. Superhigh-frequency surface-acoustic-wave transducers using AlN layers grown on SiC substrates. Appl. Phys. Lett. 2002, 81, 2538–2540. [Google Scholar] [CrossRef] [Green Version]
  60. Bu, G.; Ciplys, D.; Shur, M.; Schowalter, L.J.; Schujman, S.; Gaska, R. Surface acoustic wave velocity in single-crystal AlN substrates. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2006, 53, 251–254. [Google Scholar] [CrossRef] [PubMed]
Figure 1. A schematic of the experimental setup of an AOTF-based spectrometer.
Figure 1. A schematic of the experimental setup of an AOTF-based spectrometer.
Sensors 21 07863 g001
Figure 2. A two-dimensional unit cell geometry used in FEM simulation for the multilayer structure (a) (ZnO/SiO2/Si), (b) (LiNbO3/SiO2/Si) and (c) (AlN/SiO2/Si).
Figure 2. A two-dimensional unit cell geometry used in FEM simulation for the multilayer structure (a) (ZnO/SiO2/Si), (b) (LiNbO3/SiO2/Si) and (c) (AlN/SiO2/Si).
Sensors 21 07863 g002
Figure 3. Simulated results of AlN/Diamond structure (a) phase velocity (Rayleigh mode) versus normalised thickness of AlN (b) k2 (Rayleigh mode) versus normalised thickness of AlN.
Figure 3. Simulated results of AlN/Diamond structure (a) phase velocity (Rayleigh mode) versus normalised thickness of AlN (b) k2 (Rayleigh mode) versus normalised thickness of AlN.
Sensors 21 07863 g003
Figure 4. Zoom-in view of wave confinement and displacement profiles of structure A (ZnO/SiO2/Si) in Rayleigh mode with the normalised thickness of ZnO (a) 0.01, (b) 0.1 and (c) 2.
Figure 4. Zoom-in view of wave confinement and displacement profiles of structure A (ZnO/SiO2/Si) in Rayleigh mode with the normalised thickness of ZnO (a) 0.01, (b) 0.1 and (c) 2.
Sensors 21 07863 g004
Figure 5. Zoom-in view of wave confinement and displacement profiles of structure B (LiNbO3/SiO2/Si) in Rayleigh mode with normalised thickness of LiNbO3 (a) 0.01, (b) 0.1 and (c) 2.
Figure 5. Zoom-in view of wave confinement and displacement profiles of structure B (LiNbO3/SiO2/Si) in Rayleigh mode with normalised thickness of LiNbO3 (a) 0.01, (b) 0.1 and (c) 2.
Sensors 21 07863 g005
Figure 6. Zoom-in view of wave confinement and displacement profiles of structure C (AlN/SiO2/Si) in Rayleigh mode with normalized thickness of AlN (a) 0.01, (b) 0.1 and (c) 2.
Figure 6. Zoom-in view of wave confinement and displacement profiles of structure C (AlN/SiO2/Si) in Rayleigh mode with normalized thickness of AlN (a) 0.01, (b) 0.1 and (c) 2.
Sensors 21 07863 g006
Figure 7. Phase velocity versus normalised thickness of SiO2 on different normalised thicknesses of (a) ZnO (b) LiNbO3 (c) AlN.
Figure 7. Phase velocity versus normalised thickness of SiO2 on different normalised thicknesses of (a) ZnO (b) LiNbO3 (c) AlN.
Sensors 21 07863 g007
Figure 8. Electromechanical coupling coefficient (k2) versus normalised thickness of SiO2 on different normalised thicknesses of (a) ZnO (b) LiNbO3 (c) AlN.
Figure 8. Electromechanical coupling coefficient (k2) versus normalised thickness of SiO2 on different normalised thicknesses of (a) ZnO (b) LiNbO3 (c) AlN.
Sensors 21 07863 g008
Table 1. Dimensions of the device.
Table 1. Dimensions of the device.
DimensionsValue (µm)
Wavelength4 (λ)
Pitch of electrode2 (λ/2)
Width of IDT1 (λ/4)
Substrate thickness40 (10 λ)
Table 2. Boundary conditions (BC).
Table 2. Boundary conditions (BC).
Mechanical BCElectrical BC
Γ1FreeZero Charge/Symmetry
Γ2ContinuityContinuity
Γ3ContinuityContinuity
ΓFixed ConstraintGround
ΓL, ΓRPeriodic BC-
Table 3. Constants of materials used (AlN, LiNbO3, ZnO, SiO2 and Si).
Table 3. Constants of materials used (AlN, LiNbO3, ZnO, SiO2 and Si).
ParameterSymbolAlNLiNbO3ZnOSiO2Si
Density (kg/m3)ρ32604644560622002330
Elastic constants (GPa)c113452.03209.778.5166
c121250.5312.116.164
c131200.75105.416.164
c333952.45211.278.5166
c441180.642.431.280
c66110--31.280
Piezoelectric constants (C/m2)e15−0.484.1607−0.45--
e31−0.450.8661−0.51
e33−1.553.71.22
Dielectric constants (10−11 F/m)ε11943.68.553.3210.62
ε331129.1610.23.3210.62
Refractive indexN2.12.2032.0151.53.88
Acoustic velocity (m/s)v57606487.66133.537505000
Acousto-optic figure of merit 10−15 [s3/kg]M2-10.4182.5720.59-
Table 4. AOFM, acoustic phase velocity and k2 of structures A, B and C (with thickness 0.4 µm for SiO2, ZnO, LiNbO3, and AlN layer).
Table 4. AOFM, acoustic phase velocity and k2 of structures A, B and C (with thickness 0.4 µm for SiO2, ZnO, LiNbO3, and AlN layer).
StructureFigure of Merit (s3/kg)Phase Velocity (m/s)Electromechanical Coupling Coefficient
AZnO/SiO2/Si1.23 × 10−143.63 × 1034.20
BLiNbO3/SiO2/Si5.09 × 10−143.83 × 1035.03
CAlN/SiO2/Si7.04 × 10−144.42 × 1037.15
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hanif, M.; Jeoti, V.; Ahmad, M.R.; Aslam, M.Z.; Qureshi, S.; Stojanovic, G. FEM Analysis of Various Multilayer Structures for CMOS Compatible Wearable Acousto-Optic Devices. Sensors 2021, 21, 7863. https://doi.org/10.3390/s21237863

AMA Style

Hanif M, Jeoti V, Ahmad MR, Aslam MZ, Qureshi S, Stojanovic G. FEM Analysis of Various Multilayer Structures for CMOS Compatible Wearable Acousto-Optic Devices. Sensors. 2021; 21(23):7863. https://doi.org/10.3390/s21237863

Chicago/Turabian Style

Hanif, Mehwish, Varun Jeoti, Mohamad Radzi Ahmad, Muhammad Zubair Aslam, Saima Qureshi, and Goran Stojanovic. 2021. "FEM Analysis of Various Multilayer Structures for CMOS Compatible Wearable Acousto-Optic Devices" Sensors 21, no. 23: 7863. https://doi.org/10.3390/s21237863

APA Style

Hanif, M., Jeoti, V., Ahmad, M. R., Aslam, M. Z., Qureshi, S., & Stojanovic, G. (2021). FEM Analysis of Various Multilayer Structures for CMOS Compatible Wearable Acousto-Optic Devices. Sensors, 21(23), 7863. https://doi.org/10.3390/s21237863

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