A hybrid Schrödinger/Gaussian beam solver for quantum barriers and
surface hopping

Shi Jin; Peng Qi

doi:10.3934/krm.2011.4.1097

Previous Article
Analysis of a diffusive effective mass model for nanowires
KRM Home
This Issue
Next Article
Averaged kinetic models for flows on unstructured networks

December 2011, 4(4): 1097-1120. doi: 10.3934/krm.2011.4.1097

A hybrid Schrödinger/Gaussian beam solver for quantum barriers and surface hopping

Shi Jin ^1, and Peng Qi ^2,

1.	Department of Mathematics, University of Wisconsin-Madison, 480 Lincoln Drive, Madison, WI 53706, United States
2.	Department of Mathematics, University of Wisconsin, 480 Lincoln Drive, Madison, WI 53706, United States

Received May 2011 Revised October 2011 Published November 2011

Abstract
Related Papers

In this paper, we propose a hybrid method coupling a Schrödinger solver and the Gaussian beam method for the numerical simulation of quantum tunneling through potential barriers or surface hopping across electronic potential energy surfaces. The idea is to use a Schrödinger solver near potential barriers or zones where potential energy surfaces cross, and the Gaussian beam method--which is much more efficient than a direct Schrödinger solver--elsewhere. Buffer zones are used to convert data between the Schrödinger solver and the Gaussian beam solver. Numerical examples show that this method indeed captures quantum tunneling and surface hopping accurately, with a computational cost much lower than a direct quantum solver in the entire domain.

Keywords: Gaussian beam, Surface hopping, Time-splitting spectral method, quantum tunneling..

Mathematics Subject Classification: Primary: 58F15, 58F17; Secondary: 53C3.

Citation: Shi Jin, Peng Qi. A hybrid Schrödinger/Gaussian beam solver for quantum barriers and surface hopping. Kinetic & Related Models, 2011, 4 (4) : 1097-1120. doi: 10.3934/krm.2011.4.1097

References:

[1]	M. Baer, "Beyond Born-Oppenheimer: Electronic Nonadiabatic Coupling Terms and Conical Intersections,'', Wiley, (2006). doi: 10.1002/0471780081. Google Scholar
[2]	W. Bao, S. Jin and P. Markowich, On the time-splitting spectral approximations for the Schrödinger equation in the semiclassical regime,, J. Comput. Phys., 175 (2002), 487. doi: 10.1006/jcph.2001.6956. Google Scholar
[3]	M. Born and R. Oppenheimer, Zur Quantentheorie der Molekeln,, Ann. Phys., 84 (1927), 457. doi: 10.1002/andp.19273892002. Google Scholar
[4]	N. Ben Abdallah, P. Degond and I. M. Gamba, Coupling one-dimensional time-dependent classical and quantum transport models,, J. Math. Phys., 43 (2002), 1. doi: 10.1063/1.1421635. Google Scholar
[5]	M. V. Berry, Quantal phase factors accompanying adiabatic changes,, Proc. R. Soc. Lond Ser. A, 392 (1984), 45. doi: 10.1098/rspa.1984.0023. Google Scholar
[6]	V. Cerveny, M. M. Popov and I. Psencik, Computation of wave fields in inhomogeneous media-Gaussian beam approach,, Geophys. J. R. Astr. Soc., 70 (1982), 109. Google Scholar
[7]	K. Drukker, Basics of surface hopping in mixed quantum/classical simulations,, J. Comp. Phys., 153 (1999), 225. doi: 10.1006/jcph.1999.6287. Google Scholar
[8]	P. Gérard, P. A. Markowich, N. J. Mauser and F. Poupaud, Homogenization limits and Wigner transforms,, Communications on Pure and Applied Mathematics, 50 (1997), 323. doi: 10.1002/(SICI)1097-0312(199704)50:4<323::AID-CPA4>3.0.CO;2-C. Google Scholar
[9]	E. J. Heller, Cellular dynamics: A new semiclassical approach to time-dependent quantum mechanics,, J. Chem. Phys., 94 (1991), 2723. doi: 10.1063/1.459848. Google Scholar
[10]	N. R. Hill, Gaussian beam migration,, Geophys., 55 (1990), 1416. doi: 10.1190/1.1442788. Google Scholar
[11]	I. Horenko, C. Salzmann, B. Schmidt and C. Schutte, Quantum-classical Liouville approach to molecular dynamics: Surface hopping Gaussian phase-space packets,, J. Chem. Phys., 117 (2002), 11075. doi: 10.1063/1.1522712. Google Scholar
[12]	S. Jin, P. Markowich and C. Sparber, Mathematical and computational methods for semiclassical Schrödinger equations,, Acta Numerica, 20 (2011), 211. Google Scholar
[13]	S. Jin and K. Novak, A semiclassical transport model for thin quantum barriers,, Multiscale Modeling and Simulation, 5 (2006), 1063. doi: 10.1137/060653214. Google Scholar
[14]	S. Jin and K. Novak, A semiclassical transport model for two-dimensional thin quantum barriers,, J. Comp. Phys., 226 (2007), 1623. doi: 10.1016/j.jcp.2007.06.006. Google Scholar
[15]	S. Jin and K. Novak, A coherent semiclassical transport model for pure-state quantum scattering,, Comm. Math. Sci., 8 (2010), 253. Google Scholar
[16]	S. Jin, P. Qi and Z. Zhang, An Eulerian surface hopping method for the Schrödinger equation with conical crossings,, SIAM Multiscale Modeling & Simulation, 9 (2011), 258. Google Scholar
[17]	S. Jin, H. Wu and X. Yang, Gaussian beam methods for the Schrödinger equation in the semi-classical regime: Lagrangian and Eulerian formulations,, Comm. Math. Sci., 6 (2008), 995. Google Scholar
[18]	S. Jin and D. Yin, Computation of high frequency wave diffraction by a half plane via the Liouville equation and geometric theory of diffraction,, Communications in Computational Physics, 4 (2008), 1106. Google Scholar
[19]	L. Landau, Zur Theorie der Energiebertragung II,, Physics of the Soviet Union, 2 (1932), 46. Google Scholar
[20]	C. F. Kammerer and C. Lasser, Wigner measures and codimension two crossings,, Jour. Math. Phys., 44 (2003), 507. doi: 10.1063/1.1527221. Google Scholar
[21]	C. F. Kammerer and C. Lasser, Single switch surface hopping for molecular dynamics with transitions,, J. Chem. Phys., 128 (2008), 5. Google Scholar
[22]	C. Lasser, T. Swart and S. Teufel, Construction and validation of a rigorous surface hopping algorithm for conical crossings,, Comm. Math. Sci., 5 (2007), 789. Google Scholar
[23]	S. Leung, J. Qian and R. Burridge, Eulerian Gaussian beams for high frequency wave propagation,, Geophysics, 72 (2007), 61. doi: 10.1190/1.2752136. Google Scholar
[24]	J. Lu and X. Yang, Frozen Gaussian approximation for high frequency wave propagation,, Commun. Math. Sci., 9 (2011), 663. Google Scholar
[25]	M. Motamed and O. Runborg, Taylor expansion and discretization errors in Gaussian beam superposition,, Wave Motion, 47 (2010), 421. doi: 10.1016/j.wavemoti.2010.02.001. Google Scholar
[26]	W. H. Miller and Thomas F. George, Semiclassical theory of electronic transitions in low energy atomic and molecular collisions involving several nuclear degrees of freedom,, Journal of Chemical Physics, 56 (1972), 5637. doi: 10.1063/1.1677083. Google Scholar
[27]	B. N. Parlett, "The Symmetric Eigenvalue Problem,'', Classics ed., (1998). Google Scholar
[28]	M. M. Popov, A new method of computation of wave fields using Gaussian beams,, Wave Motion, 4 (1982), 85. doi: 10.1016/0165-2125(82)90016-6. Google Scholar
[29]	J. Qian and L. Ying, Fast Gaussian wavepacket transforms and Gaussian beams for the Schrödinger equation,, J. Comput. Phys., 229 (2010), 7848. doi: 10.1016/j.jcp.2010.06.043. Google Scholar
[30]	J. Ralston, Gaussian beams and the propagation of singularities,, in, 23 (1982), 206. Google Scholar
[31]	D. Sholla and J. Tully, A generalized surface hopping method,, J. Chem. Phys., 109 (1998). doi: 10.1063/1.477416. Google Scholar
[32]	J. Tully and R. Preston, Trajectory surface hopping approach to nonadiabatic molecular collisions: The reaction of $H^+$ with $D_2$,, J. Chem. Phys., 55 (1971), 562. doi: 10.1063/1.1675788. Google Scholar
[33]	J. Tully, Molecular dynamics with electronic transitions,, J. Chem. Phys., 93 (1990). doi: 10.1063/1.459170. Google Scholar
[34]	N. Tanushev, Superpositions and higher order Gaussian beams,, Commun. Math. Sci., 6 (2008), 449. Google Scholar
[35]	N. Tanushev, J. Qian and J. Ralston, Mountain waves and Gaussian beams,, SIAM J. Multiscale Model. Simul., 6 (2007), 688. doi: 10.1137/060673667. Google Scholar
[36]	N. M. Tanushev, R. Tsai and B. Engquist, Coupled finite difference-Gaussian beam method for high frequency wave propagation,, UCLA CAM Reports, (2010), 10. Google Scholar
[37]	N. M. Tanushev, B. Engquist and R. Tsai, Gaussian beam decomposition of high frequency wave fields,, Journal of Computational Physics, 228 (2009), 8856. doi: 10.1016/j.jcp.2009.08.028. Google Scholar
[38]	A. Voronin, J. Marques and A. Varandas, Trajectory surface hopping study of the $Li$ + $Li_2$ ($X^1\Sigma_g^+$) Dissociation reaction,, J. Phys. Chem. A, 102 (1998), 6057. doi: 10.1021/jp9805860. Google Scholar
[39]	D. Wei and X. Yang, Gaussian beam method for high frequency wave propagation in heterogeneous media,, preprint 2011., (2011). Google Scholar
[40]	W. R. E. Weiss and G. A. Hagedorn, Reflection and transmission of high freuency pulses at an interface,, Transport Theory and Statistical Physics, 14 (1985), 539. doi: 10.1080/00411458508211692. Google Scholar
[41]	D. Yin and C. Zheng, Gaussian beam formulations and interface conditions for the one-dimensional linear Schrödinger equation,, Wave Motion, 48 (2011), 310. doi: 10.1016/j.wavemoti.2010.11.006. Google Scholar
[42]	D. Yin and C. Zheng, Composite Gaussian beam approximation method for multi-phased wave functions,, preprint, (2011). Google Scholar
[43]	C. Zener, Non-adiabatic crossing of energy levels,, Proceedings of the Royal Society of London, 137 (1932), 692. Google Scholar

show all references

References:

[1]	M. Baer, "Beyond Born-Oppenheimer: Electronic Nonadiabatic Coupling Terms and Conical Intersections,'', Wiley, (2006). doi: 10.1002/0471780081. Google Scholar
[2]	W. Bao, S. Jin and P. Markowich, On the time-splitting spectral approximations for the Schrödinger equation in the semiclassical regime,, J. Comput. Phys., 175 (2002), 487. doi: 10.1006/jcph.2001.6956. Google Scholar
[3]	M. Born and R. Oppenheimer, Zur Quantentheorie der Molekeln,, Ann. Phys., 84 (1927), 457. doi: 10.1002/andp.19273892002. Google Scholar
[4]	N. Ben Abdallah, P. Degond and I. M. Gamba, Coupling one-dimensional time-dependent classical and quantum transport models,, J. Math. Phys., 43 (2002), 1. doi: 10.1063/1.1421635. Google Scholar
[5]	M. V. Berry, Quantal phase factors accompanying adiabatic changes,, Proc. R. Soc. Lond Ser. A, 392 (1984), 45. doi: 10.1098/rspa.1984.0023. Google Scholar
[6]	V. Cerveny, M. M. Popov and I. Psencik, Computation of wave fields in inhomogeneous media-Gaussian beam approach,, Geophys. J. R. Astr. Soc., 70 (1982), 109. Google Scholar
[7]	K. Drukker, Basics of surface hopping in mixed quantum/classical simulations,, J. Comp. Phys., 153 (1999), 225. doi: 10.1006/jcph.1999.6287. Google Scholar
[8]	P. Gérard, P. A. Markowich, N. J. Mauser and F. Poupaud, Homogenization limits and Wigner transforms,, Communications on Pure and Applied Mathematics, 50 (1997), 323. doi: 10.1002/(SICI)1097-0312(199704)50:4<323::AID-CPA4>3.0.CO;2-C. Google Scholar
[9]	E. J. Heller, Cellular dynamics: A new semiclassical approach to time-dependent quantum mechanics,, J. Chem. Phys., 94 (1991), 2723. doi: 10.1063/1.459848. Google Scholar
[10]	N. R. Hill, Gaussian beam migration,, Geophys., 55 (1990), 1416. doi: 10.1190/1.1442788. Google Scholar
[11]	I. Horenko, C. Salzmann, B. Schmidt and C. Schutte, Quantum-classical Liouville approach to molecular dynamics: Surface hopping Gaussian phase-space packets,, J. Chem. Phys., 117 (2002), 11075. doi: 10.1063/1.1522712. Google Scholar
[12]	S. Jin, P. Markowich and C. Sparber, Mathematical and computational methods for semiclassical Schrödinger equations,, Acta Numerica, 20 (2011), 211. Google Scholar
[13]	S. Jin and K. Novak, A semiclassical transport model for thin quantum barriers,, Multiscale Modeling and Simulation, 5 (2006), 1063. doi: 10.1137/060653214. Google Scholar
[14]	S. Jin and K. Novak, A semiclassical transport model for two-dimensional thin quantum barriers,, J. Comp. Phys., 226 (2007), 1623. doi: 10.1016/j.jcp.2007.06.006. Google Scholar
[15]	S. Jin and K. Novak, A coherent semiclassical transport model for pure-state quantum scattering,, Comm. Math. Sci., 8 (2010), 253. Google Scholar
[16]	S. Jin, P. Qi and Z. Zhang, An Eulerian surface hopping method for the Schrödinger equation with conical crossings,, SIAM Multiscale Modeling & Simulation, 9 (2011), 258. Google Scholar
[17]	S. Jin, H. Wu and X. Yang, Gaussian beam methods for the Schrödinger equation in the semi-classical regime: Lagrangian and Eulerian formulations,, Comm. Math. Sci., 6 (2008), 995. Google Scholar
[18]	S. Jin and D. Yin, Computation of high frequency wave diffraction by a half plane via the Liouville equation and geometric theory of diffraction,, Communications in Computational Physics, 4 (2008), 1106. Google Scholar
[19]	L. Landau, Zur Theorie der Energiebertragung II,, Physics of the Soviet Union, 2 (1932), 46. Google Scholar
[20]	C. F. Kammerer and C. Lasser, Wigner measures and codimension two crossings,, Jour. Math. Phys., 44 (2003), 507. doi: 10.1063/1.1527221. Google Scholar
[21]	C. F. Kammerer and C. Lasser, Single switch surface hopping for molecular dynamics with transitions,, J. Chem. Phys., 128 (2008), 5. Google Scholar
[22]	C. Lasser, T. Swart and S. Teufel, Construction and validation of a rigorous surface hopping algorithm for conical crossings,, Comm. Math. Sci., 5 (2007), 789. Google Scholar
[23]	S. Leung, J. Qian and R. Burridge, Eulerian Gaussian beams for high frequency wave propagation,, Geophysics, 72 (2007), 61. doi: 10.1190/1.2752136. Google Scholar
[24]	J. Lu and X. Yang, Frozen Gaussian approximation for high frequency wave propagation,, Commun. Math. Sci., 9 (2011), 663. Google Scholar
[25]	M. Motamed and O. Runborg, Taylor expansion and discretization errors in Gaussian beam superposition,, Wave Motion, 47 (2010), 421. doi: 10.1016/j.wavemoti.2010.02.001. Google Scholar
[26]	W. H. Miller and Thomas F. George, Semiclassical theory of electronic transitions in low energy atomic and molecular collisions involving several nuclear degrees of freedom,, Journal of Chemical Physics, 56 (1972), 5637. doi: 10.1063/1.1677083. Google Scholar
[27]	B. N. Parlett, "The Symmetric Eigenvalue Problem,'', Classics ed., (1998). Google Scholar
[28]	M. M. Popov, A new method of computation of wave fields using Gaussian beams,, Wave Motion, 4 (1982), 85. doi: 10.1016/0165-2125(82)90016-6. Google Scholar
[29]	J. Qian and L. Ying, Fast Gaussian wavepacket transforms and Gaussian beams for the Schrödinger equation,, J. Comput. Phys., 229 (2010), 7848. doi: 10.1016/j.jcp.2010.06.043. Google Scholar
[30]	J. Ralston, Gaussian beams and the propagation of singularities,, in, 23 (1982), 206. Google Scholar
[31]	D. Sholla and J. Tully, A generalized surface hopping method,, J. Chem. Phys., 109 (1998). doi: 10.1063/1.477416. Google Scholar
[32]	J. Tully and R. Preston, Trajectory surface hopping approach to nonadiabatic molecular collisions: The reaction of $H^+$ with $D_2$,, J. Chem. Phys., 55 (1971), 562. doi: 10.1063/1.1675788. Google Scholar
[33]	J. Tully, Molecular dynamics with electronic transitions,, J. Chem. Phys., 93 (1990). doi: 10.1063/1.459170. Google Scholar
[34]	N. Tanushev, Superpositions and higher order Gaussian beams,, Commun. Math. Sci., 6 (2008), 449. Google Scholar
[35]	N. Tanushev, J. Qian and J. Ralston, Mountain waves and Gaussian beams,, SIAM J. Multiscale Model. Simul., 6 (2007), 688. doi: 10.1137/060673667. Google Scholar
[36]	N. M. Tanushev, R. Tsai and B. Engquist, Coupled finite difference-Gaussian beam method for high frequency wave propagation,, UCLA CAM Reports, (2010), 10. Google Scholar
[37]	N. M. Tanushev, B. Engquist and R. Tsai, Gaussian beam decomposition of high frequency wave fields,, Journal of Computational Physics, 228 (2009), 8856. doi: 10.1016/j.jcp.2009.08.028. Google Scholar
[38]	A. Voronin, J. Marques and A. Varandas, Trajectory surface hopping study of the $Li$ + $Li_2$ ($X^1\Sigma_g^+$) Dissociation reaction,, J. Phys. Chem. A, 102 (1998), 6057. doi: 10.1021/jp9805860. Google Scholar
[39]	D. Wei and X. Yang, Gaussian beam method for high frequency wave propagation in heterogeneous media,, preprint 2011., (2011). Google Scholar
[40]	W. R. E. Weiss and G. A. Hagedorn, Reflection and transmission of high freuency pulses at an interface,, Transport Theory and Statistical Physics, 14 (1985), 539. doi: 10.1080/00411458508211692. Google Scholar
[41]	D. Yin and C. Zheng, Gaussian beam formulations and interface conditions for the one-dimensional linear Schrödinger equation,, Wave Motion, 48 (2011), 310. doi: 10.1016/j.wavemoti.2010.11.006. Google Scholar
[42]	D. Yin and C. Zheng, Composite Gaussian beam approximation method for multi-phased wave functions,, preprint, (2011). Google Scholar
[43]	C. Zener, Non-adiabatic crossing of energy levels,, Proceedings of the Royal Society of London, 137 (1932), 692. Google Scholar

[1]	Dongsheng Yin, Min Tang, Shi Jin. The Gaussian beam method for the wigner equation with discontinuous potentials. Inverse Problems & Imaging, 2013, 7 (3) : 1051-1074. doi: 10.3934/ipi.2013.7.1051
[2]	Zhong Tan, Huaqiao Wang, Yucong Wang. Time-splitting methods to solve the Hall-MHD systems with Lévy noises. Kinetic & Related Models, 2019, 12 (1) : 243-267. doi: 10.3934/krm.2019011
[3]	Masaru Ikehata. On finding the surface admittance of an obstacle via the time domain enclosure method. Inverse Problems & Imaging, 2019, 13 (2) : 263-284. doi: 10.3934/ipi.2019014
[4]	Chaoxu Pei, Mark Sussman, M. Yousuff Hussaini. A space-time discontinuous Galerkin spectral element method for the Stefan problem. Discrete & Continuous Dynamical Systems - B, 2018, 23 (9) : 3595-3622. doi: 10.3934/dcdsb.2017216
[5]	Zhonghua Qiao, Xuguang Yang. A multiple-relaxation-time lattice Boltzmann method with Beam-Warming scheme for a coupled chemotaxis-fluid model. Electronic Research Archive, 2020, 28 (3) : 1207-1225. doi: 10.3934/era.2020066
[6]	Hongjun Gao, Fei Liang. On the stochastic beam equation driven by a Non-Gaussian Lévy process. Discrete & Continuous Dynamical Systems - B, 2014, 19 (4) : 1027-1045. doi: 10.3934/dcdsb.2014.19.1027
[7]	Jean-Luc Akian, Radjesvarane Alexandre, Salma Bougacha. A Gaussian beam approach for computing Wigner measures in convex domains. Kinetic & Related Models, 2011, 4 (3) : 589-631. doi: 10.3934/krm.2011.4.589
[8]	Corinna Burkard, Roland Potthast. A time-domain probe method for three-dimensional rough surface reconstructions. Inverse Problems & Imaging, 2009, 3 (2) : 259-274. doi: 10.3934/ipi.2009.3.259
[9]	Jingwei Hu, Shi Jin. On kinetic flux vector splitting schemes for quantum Euler equations. Kinetic & Related Models, 2011, 4 (2) : 517-530. doi: 10.3934/krm.2011.4.517
[10]	Eskil Hansen, Alexander Ostermann. Dimension splitting for time dependent operators. Conference Publications, 2009, 2009 (Special) : 322-332. doi: 10.3934/proc.2009.2009.322
[11]	Moulay Rchid Sidi Ammi, Ismail Jamiai. Finite difference and Legendre spectral method for a time-fractional diffusion-convection equation for image restoration. Discrete & Continuous Dynamical Systems - S, 2018, 11 (1) : 103-117. doi: 10.3934/dcdss.2018007
[12]	Jie Tang, Ziqing Xie, Zhimin Zhang. The long time behavior of a spectral collocation method for delay differential equations of pantograph type. Discrete & Continuous Dynamical Systems - B, 2013, 18 (3) : 797-819. doi: 10.3934/dcdsb.2013.18.797
[13]	Yanzhao Cao, Li Yin. Spectral Galerkin method for stochastic wave equations driven by space-time white noise. Communications on Pure & Applied Analysis, 2007, 6 (3) : 607-617. doi: 10.3934/cpaa.2007.6.607
[14]	Liejune Shiau, Roland Glowinski. Operator splitting method for friction constrained dynamical systems. Conference Publications, 2005, 2005 (Special) : 806-815. doi: 10.3934/proc.2005.2005.806
[15]	Petra Csomós, Hermann Mena. Fourier-splitting method for solving hyperbolic LQR problems. Numerical Algebra, Control & Optimization, 2018, 8 (1) : 17-46. doi: 10.3934/naco.2018002
[16]	Zhili Ge, Gang Qian, Deren Han. Global convergence of an inexact operator splitting method for monotone variational inequalities. Journal of Industrial & Management Optimization, 2011, 7 (4) : 1013-1026. doi: 10.3934/jimo.2011.7.1013
[17]	Lijian Jiang, Craig C. Douglas. Analysis of an operator splitting method in 4D-Var. Conference Publications, 2009, 2009 (Special) : 394-403. doi: 10.3934/proc.2009.2009.394
[18]	Kai Wang, Lingling Xu, Deren Han. A new parallel splitting descent method for structured variational inequalities. Journal of Industrial & Management Optimization, 2014, 10 (2) : 461-476. doi: 10.3934/jimo.2014.10.461
[19]	Hao Chen, Kaitai Li, Yuchuan Chu, Zhiqiang Chen, Yiren Yang. A dimension splitting and characteristic projection method for three-dimensional incompressible flow. Discrete & Continuous Dynamical Systems - B, 2019, 24 (1) : 127-147. doi: 10.3934/dcdsb.2018111
[20]	Xiao Ding, Deren Han. A modification of the forward-backward splitting method for maximal monotone mappings. Numerical Algebra, Control & Optimization, 2013, 3 (2) : 295-307. doi: 10.3934/naco.2013.3.295

2019 Impact Factor: 1.311

American Institute of Mathematical Sciences

A hybrid Schrödinger/Gaussian beam solver for quantum barriers and surface hopping

References:

References:

Tools

Metrics

Other articles
by authors

Tools

Metrics

Other articles
by authors

American Institute of Mathematical Sciences

A hybrid Schrödinger/Gaussian beam solver for quantum barriers and surface hopping

References:

References:

Tools

Metrics

Other articlesby authors

Tools

Metrics

Other articlesby authors

Other articles
by authors

Other articles
by authors