View source for Physics:Surface hopping

{{Short description|Molecular dynamics simulation technique}}
'''Surface hopping''' is a [[Chemistry:Mixed quantum-classical dynamics|mixed quantum-classical technique]] that incorporates [[Physics:Quantum mechanics|quantum mechanical]] effects into [[Physics:Molecular dynamics|molecular dynamics]] simulations.<ref name=herman_1984 /><ref name=tully_1990 /><ref name=":0">{{Cite book|title=Quantum simulations of complex many-body systems: from theory to algorithms : winter school, 25 February - 1 March 2002, Rolduc Conference Centre, Kerkrade, the Netherlands ; lecture notes|date=2002|publisher=NIC-Secretariat|others=Grotendorst, Johannes., Winter school (2002.02.25-03.01 : Kerkrade)|isbn=3000090576|location=Jülich|oclc=248502198}}</ref><ref name=barbatti_2011 /> Traditional molecular dynamics assume the Born-Oppenheimer approximation, where the lighter electrons adjust instantaneously to the motion of the nuclei. Though the Born-Oppenheimer approximation is applicable to a wide range of problems, there are several applications, such as photoexcited dynamics, [[Physics:Electron transfer|electron transfer]], and surface chemistry where this approximation falls apart. Surface hopping partially incorporates the non-adiabatic effects by including excited adiabatic surfaces in the calculations, and allowing for 'hops' between these surfaces, subject to certain criteria.

== Motivation ==

Molecular dynamics simulations numerically solve the classical [[Physics:Equations of motion|equations of motion]]. These simulations, though, assume that the forces on the electrons are derived solely by the ground adiabatic surface. Solving the time-dependent [[Physics:Schrödinger equation|Schrödinger equation]] numerically incorporates all these effects, but is computationally unfeasible when the system has many degrees of freedom. 
To tackle this issue, one approach is the mean field or Ehrenfest method, where the molecular dynamics is run on the average potential energy surface given by a linear combination of the adiabatic states. This was applied successfully for some applications, but has some important limitations. When the difference between the adiabatic states is large, then the dynamics must be primarily driven by only one surface, and not an average potential. In addition, this method also violates the principle of microscopic reversibility.<ref name=":0" />

Surface hopping accounts for these limitations by propagating an ensemble of trajectories, each one of them on a single adiabatic surface at any given time. The trajectories are allowed to 'hop' between various adiabatic states at certain times such that the [[Physics:Probability amplitude|quantum amplitudes]] for the adiabatic states follow the time dependent Schrödinger equation. The probability of these hops are dependent on the coupling between the states, and is generally significant only in the regions where the difference between adiabatic energies is small.

== Theory behind the method ==

The formulation described here is in the adiabatic representation for simplicity.<ref name=sharon_1994 /> It can easily be generalized to a different representation.
The coordinates of the system are divided into two categories: quantum (<math>\mathbf{q}</math>) and classical (<math>\mathbf{R}</math>). The [[Physics:Hamiltonian (quantum mechanics)|Hamiltonian]] of the quantum [[Physics:Degrees of freedom (mechanics)|degrees of freedom]] with [[Physics:Mass|mass]] <math>m_n</math> is defined as:
:<math> H = \sum_n -\frac{\hbar^2}{2m_n}\nabla_{q_n}^2 + V(\mathbf{q},\mathbf{R}) </math>,
where <math>V</math> describes the [[Scalar potential|potential]] for the whole system. The eigenvalues of <math>H</math> as a function of <math>\mathbf{R}</math> are called the adiabatic surfaces :<math>\phi_n(\mathbf{q};\mathbf{R})</math>. Typically, <math>\mathbf{q}</math> corresponds to the electronic degree of freedom, light atoms such as [[Chemistry:Hydrogen|hydrogen]], or high frequency vibrations such as O-H stretch. The forces in the molecular dynamics simulations are derived only from one adiabatic surface, and are given by:
:<math>\begin{align}
\mathbf{F}_{\mathbf{R}} &= -\nabla_{\mathbf{R}}\langle\phi_n|H|\phi_n\rangle \\
&= -\langle\phi_n|\nabla_{R}H|\phi_n\rangle,
\end{align}</math>
where <math>n</math> represents the chosen adiabatic surface. The last equation is derived using the Hellmann-Feynman theorem. The brackets show that the [[Integral|integral]] is done only over the quantum degrees of freedom. Choosing only one adiabatic surface is an excellent approximation if the difference between the adiabatic surfaces is large for energetically accessible regions of <math>\mathbf{R}</math>. When this is not the case, the effect of the other states become important. This effect is incorporated in the surface hopping algorithm by considering the wavefunction of the quantum degrees of freedom at time t as an expansion in the adiabatic basis:
:<math>\psi(\mathbf{q};\mathbf{R},t)=\sum_n c_n(t)\phi_n(\mathbf{q};\mathbf{R})</math>,
where <math>c_n(t)</math> are the expansion coefficients. Substituting the above equation into the time dependent Schrödinger equation gives
:<math> i\hbar\dot{c_j}=\sum_n c_n\left(V_{jn}-i\hbar\dot{\mathbf{R}}.\mathbf{d}_{jn} \right)  </math>,
where <math>V_{jn}</math> and the [[Physics:Vibronic coupling|nonadiabatic coupling vector]] <math>\mathbf{d}_{jn}</math> are given by
:<math>\begin{align}
V_{jn}&=\langle\phi_j|H|\phi_n\rangle=\langle\phi_j|H|\phi_j\rangle \delta_{jn}\\
\mathbf{d}_{jn}&=\langle\phi_j|\nabla_{\mathbf{R}}\phi_n\rangle
\end{align}</math>
The adiabatic surface can switch at any given time t based on how the quantum probabilities <math>|c_j(t)|^2</math> are changing with time. The rate of change of <math>|c_j(t)|^2</math> is given by:
:<math> \dot{|c_j(t)|^2} = \sum_n \frac{2}{\hbar} Im(a_{nj}V_{jn}) - 2Re(a_{nj}\dot{\mathbf{R}}.\mathbf{d}_{jn}) </math>,
where <math>a_{nj}=c_nc_j^*</math>. For a small time interval dt, the fractional change in <math>|c_j(t)|^2</math> is given by
:<math> \frac{|c_j(t+dt)|^2-|c_j(t)|^2}{|c_j(t)|^2} \approx \frac{dt}{a_{jj}} \sum_n \frac{2}{\hbar} Im(a_{nj}V_{jn}) - 2Re(a_{nj}\dot{\mathbf{R}}.\mathbf{d}_{jn}) </math>.
This gives the net change in flux of population from state <math>j</math>. Based on this, the probability of hopping from state j to n is proposed to be
:<math> P_{j\to n} = \frac{dt}{a_{jj}} \left(\frac{2}{\hbar} Im(a_{nj}V_{jn}) - 2Re(a_{nj}\dot{\mathbf{R}}.\mathbf{d}_{jn}) \right)</math>.
This criterion is known as the "fewest switching" algorithm, as it minimizes the number of hops required to maintain the population in various adiabatic states.

Whenever a hop takes place, the velocity is adjusted to maintain [[Physics:Conservation of energy|conservation of energy]]. To compute the direction of the change in velocity, the nuclear forces in the transition is
:<math> \begin{align}
\langle\phi_j|\nabla_{\mathbf{R}}H|\phi_n\rangle &= \nabla_{\mathbf{R}}\langle\phi_j|H|\phi_n\rangle - \langle\nabla_{\mathbf{R}}\phi_j|H|\phi_n\rangle - \langle\phi_j|H|\nabla_{\mathbf{R}}\phi_n\rangle\\
&= \nabla_{\mathbf{R}} E_j \delta_{jn} + (E_j-E_n)\mathbf{d}_{jn},
\end{align} </math>
where <math>E_j=\langle\phi_j|H|\phi_j\rangle</math> is the eigen value. For the last equality, <math>d_{jn}=-d_{nj}</math> is used. This shows that the nuclear forces acting during the hop are in the direction of the nonadiabatic coupling vector <math>\mathbf{d}_{jn}</math>. Hence <math>\mathbf{d}_{jn}</math> is a reasonable choice for the direction along which velocity should be changed.

=== Frustrated hops ===

If the velocity reduction required to conserve energy while making a hop is greater than the component of the velocity to be adjusted, then the hop is known as frustrated. In other words, a hop is frustrated if the system does not have enough energy to make the hop. Several approaches have been suggested to deal with these frustrated hops. The simplest of these is to ignore these hops.<ref name=tully_1990 /> Another suggestion is not to change the adiabatic state, but reverse the direction of the component of the velocity along the nonadiabatic coupling vector.<ref name=sharon_1994 />  Yet another approach is to allow the hop to happen if an allowed hopping point is reachable within [[Physics:Quantum fluctuation|uncertainty time]] <math>\delta t=\hbar/2\Delta E </math>, where <math>\Delta E</math> is the extra energy that the system needed to make the hop possible.<ref name=jasper_2002 /> Ignoring forbidden hops without any form of velocity reversal does not recover the correct scaling for [[Physics:Marcus theory|Marcus theory]] in the nonadiabatic limit, but a velocity reversal can usually correct the errors <ref name = Jain_2015/>

=== Decoherence time ===

Surface hopping can develop nonphysical coherences between the quantum coefficients over large time which can degrade the quality of the calculations, at times leading the incorrect scaling for [[Physics:Marcus theory|Marcus theory]].<ref name=landry_2011 /> To eliminate these errors, the quantum coefficients for the inactive state can be damped or set to zero after a predefined time has elapsed after the trajectory crosses the region where hopping has high probabilities.<ref name=sharon_1994 />

== Outline of the algorithm ==

The state of the system at any time <math>t</math> is given by the [[Physics:Phase space|phase space]] of all the classical particles, the quantum amplitudes, and the adiabatic state. The simulation broadly consists of the following steps:

Step 1. Initialize the state of the system. The classical positions and velocities are chosen based on the [[Physics:Statistical ensemble (mathematical physics)|ensemble]] required.

Step 2. Compute forces using Hellmann-Feynman theorem, and integrate the [[Physics:Equations of motion|equations of motion]] by time step <math>\Delta t</math> to obtain the classical phase space at time <math>t+\Delta t</math>.

Step 3. Integrate the Schrödinger equation to evolve quantum amplitudes from time <math>t</math> to <math>t+\Delta t</math> in increments of <math>\delta t</math>. This time step <math>\delta t</math> is typically much smaller than <math>\Delta t</math>.

Step 4. Compute probability of hopping from current state to all other states. Generate a random number, and determine whether a switch should take place. If a switch does occur, change velocities to conserve energy. Go back to step 2, till trajectories have been evolved for the desired time.

== Applications ==

The method has been applied successfully to understand dynamics of systems that include tunneling, conical intersections and [[Physics:Electron excitation|electronic excitation]].<ref name=tapavicza_2007 /><ref name=jiang_2012 /><ref name=muller_1997 /><ref name=Martens_2016 />

== Limitations and foundations ==
In practice, surface hopping  is computationally feasible only for a limited number of quantum degrees of freedom. In addition, the trajectories must have enough energy to be able to reach the regions where probability of hopping is large.

Most of the formal critique of the surface hopping method comes from the unnatural separation of classical and quantum degrees of freedom. Recent work has shown, however, that the surface hopping algorithm can be partially justified
by comparison with the Quantum Classical Liouville Equation.<ref name = Subotnik_2013/> It has further been demonstrated that spectroscopic observables can be calculated in close agreement with the formally exact hierarchical equations of motion.<ref name = Tempelaar_2013/>

== See also ==

*[[Chemistry:Computational chemistry|Computational chemistry]]
*[[Physics:Molecular dynamics|Molecular dynamics]]
*Path integral molecular dynamics
*[[Chemistry:Quantum chemistry|Quantum chemistry]]

== References ==

{{reflist|
refs=<ref name=sharon_1994>{{cite journal|last=Hammes-Schiffer|first=Sharon|author2=Tully, John C. |title=Proton transfer in solution: Molecular dynamics with quantum transitions|journal=The Journal of Chemical Physics|date=1994|volume=101|issue=6|pages=4657|doi=10.1063/1.467455|bibcode = 1994JChPh.101.4657H }}</ref>

<ref name=tully_1990>{{cite journal|last=Tully|first=John C.|title=Molecular dynamics with electronic transitions|journal=The Journal of Chemical Physics|date=1990|volume=93|issue=2|pages=1061–1071|doi=10.1063/1.459170|bibcode = 1990JChPh..93.1061T |s2cid=15191625 }}</ref>

<ref name=tapavicza_2007>{{cite journal|last=Tapavicza|first=Enrico|author2=Tavernelli, Ivano|author3=Rothlisberger, Ursula|title=Trajectory surface hopping within linear response time-dependent density-functional theory|journal=Physical Review Letters|date=2007|volume=98|issue=2|pages=023001|doi=10.1103/PhysRevLett.98.023001|pmid=17358601|bibcode=2007PhRvL..98b3001T}}</ref>

<ref name=jiang_2012>{{cite journal|last=Jiang|first=Ruomu|author2=Sibert, Edwin L. |title=Surface hopping simulation of vibrational predissociation of methanol dimer|journal=The Journal of Chemical Physics|date=2012|volume=136|issue=22|pages=224104|doi=10.1063/1.4724219|pmid=22713033|bibcode = 2012JChPh.136v4104J }}</ref>

<ref name=muller_1997>{{cite journal|last=Müller|first=Uwe|author2=Stock, Gerhard |title=Surface-hopping modeling of photoinduced relaxation dynamics on coupled potential-energy surfaces|journal=The Journal of Chemical Physics|date=22 October 1997|volume=107|issue=16|pages=6230–6245|doi=10.1063/1.474288|bibcode = 1997JChPh.107.6230M }}</ref>

<ref name=herman_1984>{{cite journal|last=Herman|first=Michael F.|title=Nonadiabatic semiclassical scattering. I. Analysis of generalized surface hopping procedures|journal=The Journal of Chemical Physics|date=1984|volume=81|issue=2|pages=754–763|doi=10.1063/1.447708|bibcode = 1984JChPh..81..754H }}</ref>

<ref name=jasper_2002>{{cite journal|last=Jasper|first=Ahren W.|author2=Stechmann, Samuel N. |author3=Truhlar, Donald G. |title=Fewest-switches with time uncertainty: A modified trajectory surface-hopping algorithm with better accuracy for classically forbidden electronic transitions|journal=The Journal of Chemical Physics|date=2002|volume=116|issue=13|pages=5424|doi=10.1063/1.1453404|bibcode = 2002JChPh.116.5424J }}</ref>

<!--unused<ref name=book_2002>{{cite book|first=ed. by Johannes Grotendorst|title=Quantum simulations of complex many-body systems: from theory to algorithms : winter school, 25 February - 1 March 2002, Rolduc Conference Centre, Kerkrade, the Netherlands ; lecture notes|url=https://archive.org/details/quantumsimulatio00jgro|url-access=limited|date=2002|publisher=NIC-Secretariat|location=Jülich|isbn=978-3-00-009057-8|pages=[https://archive.org/details/quantumsimulatio00jgro/page/n386 377]}}</ref>-->

<ref name=barbatti_2011>{{cite journal|last=Barbatti|first=Mario|title=Nonadiabatic dynamics with trajectory surface hopping method|journal=Wiley Interdisciplinary Reviews: Computational Molecular Science|volume=1|issue=4|pages=620–633|doi=10.1002/wcms.64|year=2011|s2cid=123626773 }}</ref>

<ref name=Jain_2015>{{cite journal|last=Jain|first=Amber|author2=Subotnik, Joseph  |title=Surface hopping, transition state theory, and decoherence. II. Thermal rate constants and detailed balance|journal=The Journal of Chemical Physics|date=2015|volume=143|issue=13|pages=134107|doi=10.1063/1.4930549|pmid=26450292|bibcode = 2015JChPh.143m4107J |s2cid=205207864 }}</ref>

<ref name=landry_2011>{{cite journal|last=Landry|first=Brian R.|author2=Subotnik, Joseph  |title=Standard surface hopping predicts incorrect scaling for Marcus' golden-rule rate: The decoherence problem cannot be ignored|journal=The Journal of Chemical Physics|date=2015|volume=135|issue=19|pages=191101|doi=10.1063/1.3663870|pmid=22112058|bibcode = 2011JChPh.135s1101L |doi-access=free}}</ref>

<ref name=Subotnik_2013>{{cite journal|last=Subotnik|first=Joseph E.|author2=Wenjun Ouyang | author3 = Brian R. Landry  |title=Can we derive Tully's surface-hopping algorithm from the semiclassical quantum Liouville equation? Almost, but only with decoherence|journal=The Journal of Chemical Physics|date=2013|volume=139|issue=21|pages=214107|doi=10.1063/1.4829856|pmid=24320364|bibcode = 2013JChPh.139u4107S }}</ref>

<ref name=Tempelaar_2013>{{cite journal|last=Tempelaar|first=Roel|author2=van de Vegte, Cornelis | author3 = Knoester, Jasper | author4 = Jansen, Thomas L. C.  |title=Surface hopping modeling of two-dimensional spectra|journal=The Journal of Chemical Physics|date=2013|volume=138|issue=16|pages=164106 | doi=10.1063/1.4801519 |pmid=23635110|bibcode=2013JChPh.138p4106T|url=https://pure.rug.nl/ws/files/15846607/2013JChemPhysTempelaar.pdf}}</ref>

<ref name=Martens_2016>{{cite journal|last=Martens|first=Craig C.|date=2016-07-07|title=Surface Hopping by Consensus|journal=The Journal of Physical Chemistry Letters|volume=7|issue=13|pages=2610–2615|doi=10.1021/acs.jpclett.6b01186|pmid=27345103|issn=1948-7185}}</ref>
}}

== External links ==
* [http://www.newtonx.org Newton-X: A package for Newtonian dynamics close to the crossing seam.]
* [https://www.youtube.com/user/mbarbatti Movie examples of surface hopping.]

[[Category:Quantum mechanics]]
[[Category:Molecular dynamics]]

{{Sourceattribution|Surface hopping}}