Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

The Journal of Chemical Physics

Reaction coordinate method for quantum tunneling

P.G. de Buda, Morton D. Kostin

Published Feb 15, 1980• 5 citations
Publisher

Abstract

A method based on the concept of reaction coordinates is presented for calculating the difference between the first two energy eigenvalues of the Schrödinger equation for a symmetric double-well potential. The method consists of using reaction coordinates to reduce the multidimensional Schrödinger equation to a simpler set of equations, which can then be solved numerically. The method is applied to a two-dimensional problem and the results are compared with those obtained for a special case where an exact solution is available. The agreement is very good.

DOI
10.1063/1.439478

Related papers

Exact Solutions of the Quantum Double-Square-Well Potential

Enrique Peacock-López

For a symmetrical quantum double square well potential, we find analytical expressions satisfied by the quantized energies. Graphical or numerical solutions of the former relations allow us to calculate normalization constants and construct the first eight solutions of the Schrödinger equation. With these exact solutions, we analyze quantum tunneling across a potential barrier and compare our results with the experimental data for ammonia.

The Chemical Educator•Oct 24, 2006•1 citation

Low-lying states of two-dimensional double-well potentials

Amlan K. Roy, Ajit J. Thakkar, Bimalendu Deb

Wavefunctions, energies and selected expectation values of the low-lying stationary states of two-dimensional double-well potentials are obtained from the long-time solutions of the corresponding time-dependent Schrödinger equation. The numerical method consists of transformation to a diffusion-like equation which is then solved by an alternating-direction, implicit, finite-difference method. The method is tested for cases in which the energies have been obtained by other methods. Then the dependence of the energies and other properties on the potential parameters is discussed on the basis of results for four states of 35 different sets of potential parameters including some that lead to pseudo-degeneracies.

Journal of Physics A Mathematical and General•Feb 24, 2005•13 citations

Oscillations through a barrier tunneling potential for a particle in a central potential (double potential well)

Hasan Hüseyin Erbil, Ahmet Saim Selvi

In some molecules such as ammonia, double vibrations are observed at small energies. These small vibrations are thought to arise from double well potentials. To explain these small energy gaps, double potential wells are investigated, and complex calculations are made. Until recently, since the exact solution of the Schrödinger equation was not known, approximate solutions were always obtained. The exact solution of the Schrödinger equation is now known. In previous studies, a simple procedure for the general solution of the radial Schrödinger Equation has been found for spherical symmetric potentials without making any approximation. In this article, the Schrodinger equation was solved with this new solution method. It has been applied to the oscillations through a barrier tunneling potential for a particle in a central quadratic deformed potential well (double potential well). It is seen that there is no need to look for complex potentials to explain the double wells observed in small vibrations in ammonia and similar molecules, they were compared with experimentally measured energies. The results were found to be fully compatible. It has been seen that there can always be double well potentials at small energies for all molecules, atomic nuclei, and similar particles with angular momentum other than zero.

Open Access Journal of Science•Jan 23, 2024•0 citations

Quasilinearization method applied to multidimensional quantum tunneling

Mohsen Razavy, Vincent J. Cote

We apply the quasilinearization method of Bellman and Kalaba [Quasilinearization and Nonlinear Boundary-Value Problems (Elsevier, New York, 1965)] to find approximate solutions for the multidimensional quantum tunneling for separable as well as nonseparable wave equations. By introducing the idea of the complex ``semiclassical trajectory'' which is valid for the motion over and under the barrier, and which, in the proper limit, reduces to the real classical trajectory in the allowed region, we obtain an eigenvalue equation for the characteristic wave numbers. This eigenvalue equation is similar to the corresponding equation obtained from the WKB approximation and yields complex eigenvalues with negative imaginary parts. When the barrier changes very rapidly as a function of the radial distance, we can replace the concept of the semiclassical trajectory, which may not be applicable in this case, by the concept of a complex ``quantum trajectory.'' The trajectory defined either way depends on a constant of integration, and by minimizing the action with respect to this constant we can obtain the minimum escape path. The case of two-dimensional tunneling is discussed as an example of this method.

Physical Review A•Apr 1, 1994•4 citations

Tunneling in Chemical Reactions

Hiroki Nakamura, Gennady Mil'nikov

In the case of tunneling in chemical reactions [see figure1.1(c)], the total energy of reagents cannot be generally considered to be small; treating it as the first-order term of h¯ as was successfully done in Chapters 6 and 7 does not lead to any sensible semiclassical theory. The semiclassical mechanics requires (1) construction of Lagrange manifolds, (2) determination of amultidimensional caustic surface, and (3) connection of wave functions between classically allowed and forbidden regions [59]. Unfortunately, these problems have not yet been solved mathematically. As was discussed in Section 3.1, in the framework of the ordinary WKB approach we can construct Lagrange manifolds and caustics in the case of two dimensions and the connection problem can be treated relatively well by assuming the local separability. Extension of this method to higher dimensions, however, cannot be practical. The complexvalued (complex coordinate and/or momentum) WKB theory has been devised, as briefly explained there, but cannot be simply extended to multidimensions (see, for instance, [15]). This formidable task can hardly be accomplished mainly because of the complexity of the Lagrange manifold and the fact that analytical properties of any realistic potential energy surfaces extended into the complex-valued coordinate space cannot be guaranteed. The use of a complex trajectory was nicely materialized by the double-ended complex trajectory in the classical S-matrix theory [206], in which the quantized boundary conditions are satisfied at both (initial and final) ends. Unfortunately, this method turned out not to be practically useful, because finding such trajectories is also a formidable task in a multidimensional complex coordinate (and/or momentum) space. We can, however, make a break through the difficulty of detecting caustics in multidimensional space. An efficient way of doing that has been devised along classical trajectory [47]. This is not a trivial task at all, since it is not possible to find caustics in high-dimensional space by just shooting classical trajectories. Since the caustics can be found along each classical trajectory, chemical dynamics with tunneling effects incorporated can now be treated by the on-the-fly method. In this subsection, this method is explained.

CRC Press eBooks•Jul 18, 2013•2 citations

Numerical solution of the Schrödinger equation for two-dimensional double-well oscillators

Amlan K. Roy, Ajit J. Thakkar, Bimalendu Deb

Wave functions, energies and selected expectation values of the low-lying stationary states of two-dimensional double well potentials are obtained from the long-time solutions of the corresponding time-dependent Schrödinger equation. The latter is transformed to a diffusion-like equation which is then solved by an alternating-direction, implicit, finite-difference method.

CRC Press eBooks•Apr 29, 2019•0 citations
A method based on the concept of reaction coordinates is presented for calculating the difference between the first two energy eigenvalues of the Schrödinger equation for a symmetric double-well potential. The method consists of using reaction coordinates to reduce the multidimensional Schrödinger equation to a simpler set of equations, which can then be solved numerically. The method is applied to a two-dimensional problem and the results are compared with those obtained for a special case where an exact solution is available. The agreement is very good.
Wavefunctions, energies and selected expectation values of the low-lying stationary states of two-dimensional double-well potentials are obtained from the long-time solutions of the corresponding time-dependent Schrödinger equation. The numerical method consists of transformation to a diffusion-like equation which is then solved by an alternating-direction, implicit, finite-difference method. The method is tested for cases in which the energies have been obtained by other methods. Then the dependence of the energies and other properties on the potential parameters is discussed on the basis of results for four states of 35 different sets of potential parameters including some that lead to pseudo-degeneracies.
In some molecules such as ammonia, double vibrations are observed at small energies. These small vibrations are thought to arise from double well potentials. To explain these small energy gaps, double potential wells are investigated, and complex calculations are made. Until recently, since the exact solution of the Schrödinger equation was not known, approximate solutions were always obtained. The exact solution of the Schrödinger equation is now known. In previous studies, a simple procedure for the general solution of the radial Schrödinger Equation has been found for spherical symmetric potentials without making any approximation. In this article, the Schrodinger equation was solved with this new solution method. It has been applied to the oscillations through a barrier tunneling potential for a particle in a central quadratic deformed potential well (double potential well). It is seen that there is no need to look for complex potentials to explain the double wells observed in small vibrations in ammonia and similar molecules, they were compared with experimentally measured energies. The results were found to be fully compatible. It has been seen that there can always be double well potentials at small energies for all molecules, atomic nuclei, and similar particles with angular momentum other than zero.
We apply the quasilinearization method of Bellman and Kalaba [Quasilinearization and Nonlinear Boundary-Value Problems (Elsevier, New York, 1965)] to find approximate solutions for the multidimensional quantum tunneling for separable as well as nonseparable wave equations. By introducing the idea of the complex ``semiclassical trajectory'' which is valid for the motion over and under the barrier, and which, in the proper limit, reduces to the real classical trajectory in the allowed region, we obtain an eigenvalue equation for the characteristic wave numbers. This eigenvalue equation is similar to the corresponding equation obtained from the WKB approximation and yields complex eigenvalues with negative imaginary parts. When the barrier changes very rapidly as a function of the radial distance, we can replace the concept of the semiclassical trajectory, which may not be applicable in this case, by the concept of a complex ``quantum trajectory.'' The trajectory defined either way depends on a constant of integration, and by minimizing the action with respect to this constant we can obtain the minimum escape path. The case of two-dimensional tunneling is discussed as an example of this method.
In the case of tunneling in chemical reactions [see figure1.1(c)], the total energy of reagents cannot be generally considered to be small; treating it as the first-order term of h¯ as was successfully done in Chapters 6 and 7 does not lead to any sensible semiclassical theory. The semiclassical mechanics requires (1) construction of Lagrange manifolds, (2) determination of amultidimensional caustic surface, and (3) connection of wave functions between classically allowed and forbidden regions [59]. Unfortunately, these problems have not yet been solved mathematically. As was discussed in Section 3.1, in the framework of the ordinary WKB approach we can construct Lagrange manifolds and caustics in the case of two dimensions and the connection problem can be treated relatively well by assuming the local separability. Extension of this method to higher dimensions, however, cannot be practical. The complexvalued (complex coordinate and/or momentum) WKB theory has been devised, as briefly explained there, but cannot be simply extended to multidimensions (see, for instance, [15]). This formidable task can hardly be accomplished mainly because of the complexity of the Lagrange manifold and the fact that analytical properties of any realistic potential energy surfaces extended into the complex-valued coordinate space cannot be guaranteed. The use of a complex trajectory was nicely materialized by the double-ended complex trajectory in the classical S-matrix theory [206], in which the quantized boundary conditions are satisfied at both (initial and final) ends. Unfortunately, this method turned out not to be practically useful, because finding such trajectories is also a formidable task in a multidimensional complex coordinate (and/or momentum) space. We can, however, make a break through the difficulty of detecting caustics in multidimensional space. An efficient way of doing that has been devised along classical trajectory [47]. This is not a trivial task at all, since it is not possible to find caustics in high-dimensional space by just shooting classical trajectories. Since the caustics can be found along each classical trajectory, chemical dynamics with tunneling effects incorporated can now be treated by the on-the-fly method. In this subsection, this method is explained.