Sherman, M. C. | Engineering Papers

DOE OSTI2016

Nonadiabatic transition path sampling

Fewest-switches surface hopping (FSSH) is combined with transition path sampling (TPS) to produce a new method called nonadiabatic path sampling (NAPS). The NAPS method is validated on a model electron transfer system coupled to a Langevin bath. Numerically exact rate constants are computed using the reactive flux (RF) method over a broad range of solvent frictions that span from the energy diffusion (low friction) regime to the spatial diffusion (high friction) regime. The NAPS method is shown to quantitatively reproduce the RF benchmark rate constants over the full range of solvent friction. Integrating FSSH within the TPS framework expands the applicability of both approaches and creates a new method that will be helpful in determining detailed mechanisms for nonadiabatic reactions in the condensed-phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

Thermal equilibrium properties of surface hopping with an implicit Langevin bath

The ability of fewest switches surface hopping (FSSH) approach, where the classical degrees of freedom are coupled to an implicit Langevin bath, to establish and maintain an appropriate thermal equilibrium was evaluated in the context of a three site model for electron transfer. The electron transfer model consisted of three coupled diabatic states that each depends harmonically on the collective bath coordinate. This results in three states with increasing energy in the adiabatic representation. The adiabatic populations and distributions of the collective solvent coordinate were monitored during the course of 250 ns FSSH-Langevin (FSSH-L) simulations performed at a broad range of temperatures and for three different nonadiabatic coupling strengths. The agreement between the FSSH-L simulations and numerically exact results for the adiabatic population ratios and solvent coordinate distributions was generally favorable. The FSSH-L method produces a correct Boltzmann distribution of the solvent coordinate on each of the adiabats, but the integrated populations are slightly incorrect because FSSH does not rigorously obey detailed balance. The overall agreement is better at high temperatures and for high nonadiabatic coupling, which agrees with a previously reported analytical and simulation analysis [J. R. Schmidt, P. V. Parandekar, and J. C. Tully, J. Chem. Phys. 129, 044104 (2008)] on a two-level system coupled to a classical bath.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗