Research Article | Open Access
B. Nagarjun, A. V. Sarma, G. V. Rama Rao, C. Rambabu, "Thermodynamic and Acoustic Study on Molecular Interactions in Certain Binary Liquid Systems Involving Ethyl Benzoate", Journal of Thermodynamics, vol. 2013, Article ID 285796, 9 pages, 2013. https://doi.org/10.1155/2013/285796
Thermodynamic and Acoustic Study on Molecular Interactions in Certain Binary Liquid Systems Involving Ethyl Benzoate
Speeds of sound and density for binary mixtures of ethyl benzoate (EB) with N,N-dimethylformamide (NNDMF), N,N-dimethyl acetamide (NNDMAc), and N,N-dimethylaniline (NNDMA) were measured as a function of mole fraction at temperatures 303.15, 308.15 K, 313.15 K, and 318.15 K and atmospheric pressure. From the experimental data, adiabatic compressibility (), intermolecular free length (), and molar volume () have been computed. The excess values of the above parameters were also evaluated and discussed in light of molecular interactions. Deviation in adiabatic compressibilities and excess intermolecular free length () are found to be negative over the molefraction of ethyl benzoate indicating the presence of strong interactions between the molecules. The negative excess molar volume values are attributed to strong dipole-dipole interactions between unlike molecules in the mixtures. The binary data of , , and were correlated as a function of molefraction by using the Redlich-Kister equation.
The thermophysical study of esters is of increasing interest due to their wide range in flavouring, perfumery, artificial essences, and cosmetics. Esters are also important solvents in pharmaceutical, paint, and plastic industries . Recently, interest in the study of liquid mixtures containing esters as one of the components [2–5] has increased. These studies are of great significance because one can get information regarding structural changes that occur in pure ester because of mixing.
Among the different types of esters, ethyl benzoate (aromatic ester) is a polar solvent () and is not a strongly associated liquid. On the other hand, NNDMF, to some extent, is associated by means of dipole-dipole interactions. Significant structural effects are absent due to the lack of hydrogen bonds. Therefore, it acts as an aprotic protophilic solvent . NNDMAc is a dipolar aprotic solvent and is moderately structured . The N,N-dialkylamides have no significant intermolecular hydrogen bonding capability but are highly polar with high percentage of ionic character, making oxygen of C=O group strongly negative . NNDMA is a unassociated liquid . The interactions of ethyl benzoate with N,N-dialkylamides and N,N-dimethylaniline are important to enable analysis of their dissimilar geometric structures on mixture properties.
Rao  calculated the excess isentropic compressibilities for the binary mixtures of N,N-dimethylformamide and N,N-dimethylacetamide and various normal and branched esters. The negative excess values are attributed to the presence of n-n interactions between unlike molecules. Aminabhavi et al.  studied the physicochemical behavior of binary mixtures of ethyl benzoate with diethylene glycol and dimethyl ether. Rathnam et al.  reported the densities, viscosities, and refractive indices of ethyl benzoate with o-xylene, m-xylene, p-xylene, and ethyl benzene. Lien et al.  have measured the excess molar enthalpies of the binary mixture of ethyl benzoate with 1-Octanol at 298.15 K. Joshi et al.  measured the densities and viscosities of the binary mixture of bromoform with ethyl benzoate.
We report herein the experimental values of densities and ultrasonic velocities of mixtures of ethyl benzoate with NNDMF, NNDMAc, and NNDMA at temperatures 303.15, 308.15, 313.15, and 318.15 K covering the whole miscibility range to enable analysis of the effect of their dissimilar geometric structures on the mixture properties. The non-rectilinear behavior of ultrasonic velocity, compressibility, and other thermodynamic parameters of liquid mixtures reveal the strength of interactions. To get additional information about the nature and strength of molecular interactions, other parameters such as free length, adiabatic compressibility, and their excess parameters have been calculated in the liquid mixtures. These results were correlated with the Redlich-Kister polynomial equation  to derive the coefficients and standard deviation.
2. Materials and Experimental Details
Ethyl benzoate (SD Fine Chemicals, purity > 99%) was used without any further treatment. NNDMF, NNDMAc, and NNDMA were purchased from SD Fine or Merck and were purified by the recommended methods [16, 17]. Further, the purities were ascertained from their ultrasonic speed and density values at 303.15 K which agreed with the literature values (Table 1).
|Data taken from . bData taken from . cData taken from . dData taken from .|
Mixtures were prepared by mixing appropriate volumes of liquids in airtight bottles and weighed in a single-pan Mettler balance to an accuracy of ±0.001 mg. Preferential evaporation losses of solvents from mixtures were kept to a minimum as evidenced by repeated measurement of the physical properties over an interval of 2-3 days, during which no changes in physical properties were observed. The possible error in mole fractions is estimated to be around ±0.0001.
The densities of liquids and their mixtures were measured by a 25 mL specific gravity bottle, calibrated with redistilled water. The average uncertainty in measurement in the measured density is ±0.05 kgm−3. With the fluctuation of ±0.05 K, temperature was controlled by a water thermostat. Ultrasonic velocity was measured using the ultrasonic interferometer (Model M-83) provided by Mittal Enterprises, New Delhi. The values agree closely with the values given in the literature.
The experimental values of ultrasonic speeds (), densities (), adiabatic compressibilities (), molar volumes (), and free length () of pure ethyl benzoate, NNDMF, NNDMAc, NNDMA, and those of their binary mixtures over the entire composition range and at 303.15 K, expressed by mole fraction of ethyl benzoate, are listed in Table 2.
3. Results and Discussion
The experimental density () values of binary mixtures were used to calculate the excess molar volumes as where is the molar mass; subscripts 1 and 2 stand for pure components, ethyl benzoate and NNDMF, NNDMAc, and NNDMA, respectively. The uncertainty in is estimated to be within .
Assuming that ultrasonic absorption is negligible, adiabatic compressibility () can be obtained from the density and velocity of ultrasonic sound () using the relation Deviation in adiabatic compressibility was calculated using the relation where is the compressibility of the mixture.
Intermolecular free length () has been calculated using the relation where is Jacobson’s temperature-dependant constant and is equal to .
The excess intermolecular free length () at a given mole fraction is the difference between mean free length and the sum of the fractional contributions of the two liquids given by where and are the individual intermolecular free length values of pure liquids in the binary mixtures.
The experimental values of ultrasonic velocity () and values of density (), at the four temperatures, namely, K, 308.15 K, 313.15 K, and 318.15 K, along with the derived values of adiabatic compressibility (), intermolecular free length (), and molar volume () and their excess parameters are given in Tables 2, 3, and 4.
The excess properties were fitted to a Redlich, Kister-type  polynomial equation The optimum number of coefficients, , was ascertained from an examination of the variation of the standard deviation . The values of coefficients, , were evaluated by using the method of least squares, with all points weighted equally. The coefficients , , and along with standard deviations, of fit for all the mixtures are listed in Table 5.
In this investigation, the values of ultrasonic velocity decrease with increase in the concentration of ethyl benzoate and decrease with increase in temperature at any particular concentration for all the three systems (Figure 1). The ultrasonic velocity values decrease with increase of temperature due to the breaking of hetero-and homomolecular clusters at high temperatures . Lagemann and Duban  were the first to point out the ultrasonic velocity approach for qualitative estimation of the interaction in liquids.
The intermolecular free length is the distance between the surfaces of the neighboring molecules. The variation of ultrasonic velocity in a solution depends upon the increase or decrease of intermolecular free length after mixing the components. The interdependence of intermolecular free length and ultrasonic velocity was evolved from a model for sound propagation proposed by Kincaid and Eyring . The ultrasonic velocity should decrease if the intermolecular free length increases as a result of mixing of components. This fact is observed in the present investigation for ethyl benzoate + NNDMF, ethyl benzoate + NNDMAc, and ethyl benzoate + NNDMA systems. Figure 2 represents the variation of excess intermolecular free length with mole fraction of ethyl benzoate. The excess intermolecular free length values are negative, and the curves appear to reach negative peak value at about 0.45 mole fraction of ethyl benzoate. According to Ramamurthy and Sastry , the negative values indicate that sound wave has to travel a longer distance. This may be attributed to the dominant nature of interactions between unlike molecules.
The compressibility behavior of solutes, which is the second derivative of the Gibbs energy, is a very sensitive indicator of molecular interactions and can provide useful information about this phenomenon [26–28]. The structural change of molecules takes place due to the existence of electrostatic field between interacting molecules. The change in adiabatic compressibility value in liquids and liquid mixtures may be ascribed to the strength of intermolecular attraction. The effect of depolymerization increases the compressibility of the system . Electrostatic attraction and association decrease the compressibility. The compressibility of the mixtures is the result of these two effects depending upon the predominance . According to Jacobson [31, 32], the adiabatic compressibility can be studied more through intermolecular free length.
Figure 3 represents the deviation in adiabatic compressibility with the mole fraction of ethyl benzoate. In the present investigation, deviation in adiabatic compressibilities is found to be negative over the mole fraction of ethyl benzoate indicating the presence of strong interactions between the molecules. The negative deviation in adiabatic compressibility reaches a peak at about 0.45 mole fraction of ethyl benzoate in all the three systems chosen.
The deviation in adiabatic compressibility can be explained by taking into consideration the following factors.(a)Loss of dipolar association and difference in size and shape of component molecules which lead to decrease in velocity and increase in compressibility.(b)Dipole-dipole interaction or hydrogen-bonded complex formation between unlike molecules which lead to increase in sound velocity and decrease of compressibility.
The actual deviation depends on the resultant effect. The strength of the interaction between the components increases when excess values tend to become increasingly negative. This may be qualitatively interpreted in terms of closer approach of unlike molecules leading to reductions in compressibility and volume [33, 34]. This type of interactions for the binary mixtures has been already reported previously . The deviations in adiabatic compressibilities are found to increase with increasing temperature which is in agreement with the previously reported results . The and minima occur at the same concentrations further strengthen the occurrence of molecular associations .
The sign of of a system depends upon the relative magnitude of expansion and contraction of the two liquids due to mixing .
The negative arises due to dominance of the following factors.(a)Chemical interaction between constituent molecules such as heteromolecular association through the formation of H-bond, often termed as strong specific interaction.(b)Association through weaker physical forces such as dipolar force or any other forces of this kind.(c)Accommodation of molecules of one component into the interstitial positions of the structural network of molecules of the other component.(d)Geometry of the molecular structure that favors fitting of the component molecules with each other.
The values are negative over the entire mole fraction range and at all temperatures investigated for all binary systems under study (Figure 4). The observed negative values of for the three mixtures indicate the presence of specific interactions between ethyl benzoate and amide molecules. The negative values are attributed to strong dipole-dipole interactions between unlike molecules in the mixtures. The values are more negative for ethyl benzoate + DMAc than those for other two systems.
The electron density at oxygen atom of the carbonyl atom of DMAc is greater than that of DMF due to the presence of methyl group at carbon atom of carbonyl group in DMF resulting in stronger interaction in thesystem [39, 40]. The values of ethyl benzoate + DMA are in between these two. The strength of interactions between ethyl benzoate + DMA is stronger than ethyl benzoate + DMF due to the fact that negative charge on nitrogen in DMA is more than that of nitrogen in DMF due to conjugation with C=O group (this is also evident from the fact that the DMA is more basic than DMF).
Rao and Reddy  reported increased negative values of excess molar volumes with the increase in carbon chain length of the ester, in the binary mixtures of DMF, and aliphatic esters. The observed higher negative values of () for aromatic ester with DMF in the present investigation, when compared to that of ethyl acetate and DMF () at the same temperature (303.15 K), indicate much stronger interactions between the unlike molecules of components in this binary mixture due to the formation of not only dipole-dipole interactions but also of induced polar interaction like dipole-induced dipole and dipole interactions between aromatic ester and DMF. Similar interactions are reported in the mixtures of DMF and polycyclic aromatic hydrocarbons by Nikam and Kharat [42, 43] and by Ramadevi and Prabhakara Rao . Also, formation of induced polar interactions like dipole-induced dipole interactions between polycyclic aromatic hydrocarbons and NMP (acyclic amide) is reported by Sugiura and Ogawa . This further supports that amides react with aromatic compounds more strongly than the corresponding aliphatic counter parts.
The values decrease (become more negative) with increase in temperature for all the 3 systems. This is attributed due to the dissociation of self-associated amide molecules resulting in the more favorable fitting of the smaller molecules into large voids of bigger ethyl benzoate molecules leading to contraction in volume and hence resulting in more negative values with rise in temperature.
The adiabatic compressibility () and intermolecular free length () both have an inverse relationship with ultrasonic velocity (). Occurrence of maxima, , and minima at the same concentrations indicates the strong interaction through dipole-dipole interactions between the components. The negative values are attributed to strong dipole-dipole interactions between unlike molecules in the mixtures.
The authors sincerely thank the University Grants Commission, India, for funding the current research work under UGC Scholarship Assistance Program (SAP) in the Department of Physics, Andhra University.
- Y.-W. Sheu and C. H. Tu, “Densities and viscosities of binary mixtures of ethyl acetoacetate, ethyl isovalerate, methyl benzoate, benzyl acetate, ethyl salicylate, and benzyl propionate with ethanol at T = (288.15, 298.15, 308.15, and 318.15) K,” Journal of Chemical and Engineering Data, vol. 51, no. 2, pp. 545–553, 2006.
- J. N. Nayak, M. I. Aralaguppi, and T. M. Aminabhavi, “Density, viscosity, refractive index, and speed of sound in the binary mixtures of 1,4-dioxane + ethyl acetoacetate, + diethyl oxalate, + diethyl phthalate, or + dioctyl phthalate at 298.15, 303.15, and 308.15 K,” Journal of Chemical and Engineering Data, vol. 48, no. 6, pp. 1489–1494, 2003.
- W. V. Steele, R. D. Chirico, A. B. Cowell, S. E. Knipmeyer, and A. Nguyen, “Thermodynamic properties and ideal-gas enthalpies of formation for methyl benzoate, ethyl benzoate, (R)-(+)-limonene, tert-amyl methyl ether, trans-crotonaldehyde, and diethylene glycol,” Journal of Chemical and Engineering Data, vol. 47, no. 4, pp. 667–688, 2002.
- N. Indrswati, Mudjiati, F. Wicaksana, H. Hindarso, and S. IsmadjiS, “Measurements of density and viscosity of binary mixtures of several flavor compounds with 1-butanol and 1-pentanol at 293.15 K, 303.15 K, 313.15 K, and 323.15 K,” Journal of Chemical and Engineering Data, vol. 46, no. 3, pp. 696–702, 2001.
- M. Hasan, A. P. Hiray, U. B. Kadam, D. F. Shirude, K. J. Kurhe, and A. B. Sawant, “Densities, sound speed, and IR studies of (methanol + 1-acetoxybutane) and (methanol + 1,1-dimethylethyl ester) at (298.15, 303.15, 308.15, and 313.15) K,” Journal of Chemical and Engineering Data, vol. 55, no. 1, pp. 535–538, 2010.
- M. M. H. Bhuiyan and M. H. Uddin, “Excess molar volumes and excess viscosities for mixtures of N,N-dimethylformamide with methanol, ethanol and 2-propanol at different temperatures,” Journal of Molecular Liquids, vol. 138, no. 1–3, pp. 139–146, 2008.
- P. J. Victor, D. Das, and D. K. Hazra, “Excess molar volumes, viscosity deviations and isentropic compressibility changes in binary mixtures of N,N-dimethylacetamide with 2-methoxyethanol and water in the temperature range 298.15 to 318.15 K,” Journal of the Indian Chemical Society, vol. 81, no. 12, pp. 1045–1050, 2004.
- R. K. Bhardwaj and A. Pal, “Spectroscopic studies of binary liquid mixtures of alkoxyethanols with substituted and cyclic amides at 298.15 K,” Journal of Molecular Liquids, vol. 118, no. 1–3, pp. 37–39, 2005.
- R. Palepu and J. H. MacNeil, “Viscosities and densities of binary-mixtures of 2,2,2-trichloroethanol with substituted anilines,” Australian Journal of Chemistry, vol. 41, no. 5, pp. 791–797, 1988.
- K. P. C. Rao, “Excess isentropic compressibilities of N,N-dimethylformamide and N,N-dimethylacetamide with aliphatic esters at 303.15 K,” Ultrasonics, vol. 28, no. 2, pp. 120–124, 1990.
- T. M. Aminabhavi, H. T. S. Phayde, R. S. Khinnavar, B. Gopalakrishna, and K. C. Hansen, “Densities, refractive indices, speeds of sound, and shear viscosities of diethylene glycol dimethyl ether with ethyl acetate, methyl benzoate, ethyl benzoate, and diethyl succinate in the temperature range from 298.15 to 318.15 K,” Journal of Chemical and Engineering Data, vol. 39, no. 2, pp. 251–260, 1994.
- M. V. Rathnam, S. Mohite, and M. S. S. Kumar, “Viscosity, density, and refractive index of some (ester + hydrocarbon) binary mixtures at 303.15 K and 313.15 K,” Journal of Chemical and Engineering Data, vol. 50, no. 2, pp. 325–329, 2005.
- P.-J. Lien, S. T. Wu, M. J. Lee, and H. M. Lin, “Excess molar enthalpies for dimethyl carbonate with o-xylene, m-xylene, p-xylene, ethylbenzene, or ethyl benzoate at 298.15 K and 10.2 MPa,” Journal of Chemical and Engineering Data, vol. 48, no. 3, pp. 632–636, 2003.
- S. S. Joshi, T. M. Aminabhavi, and S. S. Shukla, “Densities and viscosities of binary mixtures of bromoform with anisole, acetophenone, ethyl benzoate, 1,2- dichloroethane and 1,1,2,2-tetrachloroethane from 298.15 to 313.15 K,” Indian Journal of Technology, vol. 29, pp. 319–326, 1991.
- O. Redlich and A. T. Kister, “Algebraic representation of thermodynamic properties and the classification of solutions,” Industrial and Engineering Chemistry, vol. 40, no. 2, pp. 345–348, 1948.
- A. I. Vogel, A Text Book of Practical Organic Chemistry, John Willey, New York, NY, USA, 5th edition, 1989.
- J. A. Riddick, W. B. Bunger, and T. K. Sokano, Techniques in Chemistry, Vol 2, Organic Solvents, John Willey, NewYork, NY, USA, 4th edition, 1986.
- G. V. Rama Rao, A. V. Sarma, and C. Rambabu, “Evaluation of excess thermodynamic properties in some binary mixtures of o-chlorophenol,” Indian Journal of Chemistry A, vol. 43, no. 12, pp. 2518–2528, 2004.
- D. Papamatthaiakis, F. Aroni, and V. Havredaki, “Isentropic compressibilities of (amide + water) mixtures: a comparative study,” Journal of Chemical Thermodynamics, vol. 40, no. 1, pp. 107–118, 2008.
- R. Kumar, Ultrasonic and spectroscopic investigations of molecular interactions in liquid mixtures [Ph.D. thesis], University of Madras, Chennai, India, 2009.
- M. Katz, P. W. Lobo, A. S. Minano, and H. Solimo, “Viscosities, densities, and refractive indices of binary liquid mixtures,” Canadian Journal of Chemistry, vol. 49, no. 15, pp. 2605–2609, 1971.
- P. S. Nikam, B. S. Jagdale, A. B. Sawant, and M. Hasan, Journal of Pure and Applied Ultrasonics, vol. 22, p. 115, 2000.
- R. T. Lagemann and W. S. Duban, Journal of Chemical Physics, vol. 49, p. 423, 1945.
- J. F. Kincaid and H. Eyring, “Free volumes and free angle ratios of molecules in liquids,” The Journal of Chemical Physics, vol. 6, no. 10, pp. 620–629, 1938.
- M. Ramamurthy and O. S. Sastry, Indian Journal of Pure and Applied Physics, vol. 21, p. 579, 1983.
- I. Mozo, J. A. Gonzalez, I. G. de la Fuente, J. C. Cobos, N. Riesco, and N. Riesco, “Thermodynamics of mixtures containing alkoxyethanols. Part XXV. Densities, excess molar volumes and speeds of sound at 293.15, 298.15 and 303.15 K, and isothermal compressibilities at 298.15 K for 2-alkoxyethanol + 1-butanol systems,” Journal of Molecular Liquids, vol. 140, no. 1–3, pp. 87–100, 2008.
- V. Ulagendran, R. Kumar, S. Jayakumar, and V. Kannappan, “Ultrasonic and spectroscopic investigations of charge-transfer complexes in ternary liquid mixtures,” Journal of Molecular Liquids, vol. 148, no. 2-3, pp. 67–72, 2009.
- P. J. Victor, P. K. Muhuri, B. Das, and D. K. Hazra, “Thermodynamics of ion association and solvation in 2-methoxyethanol: behavior of tetraphenylarsonium, picrate, and tetraphenylborate ions from conductivity and ultrasonic data,” Journal of Physical Chemistry B, vol. 103, no. 50, pp. 11227–11232, 1999.
- M. V. Kaulgud and K. J. Patil, Acustica, vol. 28, p. 130, 1973.
- K. Vijayalakshmi, V. Lalitha, and R. Gandhimathi, Journal of Pure and Applied Ultrasonics, vol. 28, pp. 154–156, 2006.
- B. Jacobson, “Intermolecular free lengths in liquids in relation to compressibility, surface tension and viscosity,” Acta Chemica Scandanavia, vol. 5, pp. 1214–1216, 1951.
- B. Jacobson, Acta Chemica Scandanavia, vol. 34, p. 121, 1975.
- S. S. Yadava and A. Yadav, “Ultrasonic study on binary liquid mixtures between some bromoalkanes and hydrocarbons,” Ultrasonics, vol. 43, no. 9, pp. 732–735, 2005.
- R. J. Fort and W. R. Moore, “Adiabatic compressibilities of binary liquid mixtures,” Transactions of the Faraday Society, vol. 61, pp. 2102–2111, 1965.
- M. Yasmin, K. P. Singh, S. Parveen, M. Gupta, and J. P. Shukla, “Thermoacoustical excess properties of binary liquid mixtures-a comparative experimental and theoretical study,” Acta Physica Polonica A, vol. 115, no. 5, pp. 890–900, 2009.
- K. Subbarangaiah, N. M. Murthy, and S. V. Subrahmanyam, Acustica, vol. 55, p. 105, 1985.
- A. M. E. Raj, L. B. Resmi, V. B. Jothy, M. Jayachandran, and C. Sanjeeviraja, “Ultrasonic study on binary mixture containing dimethylformamide and methanol over the entire miscibility range () at temperatures 303–323 K,” Fluid Phase Equilibria, vol. 281, no. 1, pp. 78–86, 2009.
- M. M. H. Bhuiyan and M. H. Uddin, “Excess molar volumes and excess viscosities for mixtures of N,N-dimethylformamide with methanol, ethanol and 2-propanol at different temperatures,” Journal of Molecular Liquids, vol. 138, no. 1–3, pp. 139–146, 2008.
- A. Ali and A. K. Nain, “Ultrasonic and volumetric study of binary mixtures of benzyl alcohol with amides,” Bulletin of the Chemical Society of Japan, vol. 75, no. 4, pp. 681–687, 2002.
- N. K. Kim, H. J. Lee, K. H. Choi et al., “Substituent effect of N,N-dialkylamides on the intermolecular hydrogen bonding with thioacetamide,” Journal of Physical Chemistry A, vol. 104, no. 23, pp. 5572–5578, 2000.
- K. P. C. Rao and K. S. Reddy, “Excess volumes of N,N-dimethyl formamide and N,N-dimethyl acetamide with aliphatic esters at room temperature,” Physics and Chemistry of Liquids, vol. 18, no. 1, pp. 75–79, 1988.
- P. S. Nikam and S. J. Kharat, “Densities, viscosities, and thermodynamic properties of (N,N-dimethylformamide + benzene + chlorobenzene) ternary mixtures at (298.15, 303.15, 308.15, and 313.15) K,” Journal of Chemical and Engineering Data, vol. 48, no. 5, pp. 1202–1207, 2003.
- P. S. Nikam and S. J. Kharat, “Density and viscosity studies of binary mixtures of N,N-dimethylformamide with toluene and methyl benzoate at (298.15, 303.15, 308.15, and 313.15) K,” Journal of Chemical and Engineering Data, vol. 50, no. 2, pp. 455–459, 2005.
- R. S. Ramadevi and M. V. Prabhakara Rao, “Excess volumes of substituted benzenes with N,N-dimethylformamide,” Journal of Chemical and Engineering Data, vol. 40, no. 1, pp. 65–67, 1995.
- T. Sugiura and H. Ogawa, “Enthalpies of solution, partial molar volumes and isentropic compressibilities of polycyclic aromatic hydrocarbons in 1-methyl-2-pyrrolidone at 298.15 K,” Fluid Phase Equilibria, vol. 277, no. 1, pp. 29–34, 2009.
Copyright © 2013 B. Nagarjun et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.