Critical appraisal of the fewest switches algorithm for surface hopping

In this paper the authors address the problem of internal consistency in trajectory surface hopping methods, i.e., the requirement that the fraction of trajectories running on each electronic state equals the probabilities computed by the electronic time-dependent Schrodinger equation, after averaging over all trajectories. They derive a formula for the hopping probability in Tully's "fewest switches" spirit that would yield a rigorously consistent treatment. They show the relationship of Tully's widely used surface hopping algorithm with the "exact" prescription that cannot be applied when running each trajectory independently. They also bring out the connection of the consistency problem with the coherent propagation of the electronic wave function and the artifacts caused by coherent Rabi-type oscillations of the state probabilities in weak coupling regimes. A real molecule (azobenzene) and two ad hoc models serve as examples to illustrate the above theoretical arguments. Following a proposal by Truhlar's group [Zhu et al., J. Chem. Phys. 121, 7658 (2004) Zhu et al., J. Chem. Theory Comput. 1, 527 (2005)], they apply a decoherence correction to the state probabilities, in conjunction with Tully's algorithm, and they obtain satisfactory results in terms of internal consistency and of agreement with the outcomes of quantum wave packet calculations.

The Photoisomerization Mechanism of Azobenzene: A Semiclassical Simulation of Nonadiabatic Dynamics

We have simulated the photoisomerization dynamics of azobenzene, taking into account internal conversion and geometrical relaxation processes, by means of a semiclassical surface hopping approach. Both n-->pi* and pi-->pi* excitations and both cis-->trans and trans-->cis conversions have been considered. We show that in all cases the torsion around the N==N double bond is the preferred mechanism. The quantum yields measured are correctly reproduced and the observed differences are explained as a result of the competition between the inertia of the torsional motion and the premature deactivation of the excited state. Recent time-resolved spectroscopic experiments are interpreted in the light of the simulated dynamics.

Manipulating azobenzene photoisomerization through strong light-molecule coupling

The formation of hybrid light–molecule states (polaritons) offers a new strategy to manipulate the photochemistry of molecules. To fully exploit its potential, one needs to build a toolbox of polaritonic phenomenologies that supplement those of standard photochemistry. By means of a state-of-the-art computational photochemistry approach extended to the strong-coupling regime, here we disclose various mechanisms peculiar of polaritonic chemistry: coherent population oscillations between polaritons, quenching by trapping in dead-end polaritonic states and the alteration of the photochemical reaction pathway and quantum yields. We focus on azobenzene photoisomerization, that encompasses the essential features of complex photochemical reactions such as the presence of conical intersections and reaction coordinates involving multiple internal modes. In the strong coupling regime, a polaritonic conical intersection arises and we characterize its role in the photochemical process. Our chemically detailed simulations provide a framework to rationalize how the strong coupling impacts the photochemistry of realistic molecules.

Manipulation of the photochemistry of molecules is traditionally achieved through synthetic chemical modifications. Here the authors use computational photochemistry to show how to control azobenzene photoisomerization through hybrid light–molecule states (polaritons).

Semiempirical Hamiltonian for simulation of azobenzene photochemistry

We present a semiempirical Hamiltonian that provides an accurate description of the first singlet and triplet potential energy surfaces of azobenzene for use in direct simulations of the excited-state dynamics. The parameterization made use of spectroscopic and thermochemical data and the best ab initio results available to date. Two-dimensional potential energy surfaces based on constrained geometry optimizations are presented for the states that are most relevant for the photochemistry of azobenzene, namely, S(0), S(1), and S(2). In order to run simulations of the photodynamics of azobenzene in hydrocarbons or hydroxylic solvents, we determined the interactions of methane and methanol with the azo group by ab initio calculations and fitted the interactions with a QM/MM interaction Hamiltonian.

Are azobenzenophanes rotation-restricted?

We simulated the photoisomerization dynamics of an azobenzenophane with a semiclassical surface hopping approach and a semiempirical reparametrized quantum mechanics/molecular mechanics Hamiltonian. Only one of the two azobenzene chromophores in the molecule is taken into account quantum mechanically: the other one is treated by molecular mechanics. Both n-->pi* and pi-->pi* excitations are considered. Our results show that the photoisomerization reaction mainly involves the rotation around the N=N double bond. The excited state relaxation features are in qualitative agreement with experimental time-resolved fluorescence results.

Newton-X Platform: New Software Developments for Surface Hopping and Nuclear Ensembles

Newton-X is an open-source computational platform to perform nonadiabatic molecular dynamics based on surface hopping and spectrum simulations using the nuclear ensemble approach. Both are among the most common methodologies in computational chemistry for photophysical and photochemical investigations. This paper describes the main features of these methods and how they are implemented in Newton-X. It emphasizes the newest developments, including zero-point-energy leakage correction, dynamics on complex-valued potential energy surfaces, dynamics induced by incoherent light, dynamics based on machine-learning potentials, exciton dynamics of multiple chromophores, and supervised and unsupervised machine learning techniques. Newton-X is interfaced with several third-party quantum-chemistry programs, spanning a broad spectrum of electronic structure methods.

Dynamics of Azobenzene Dimer Photoisomerization: Electronic and Steric Effects

While azobenzenes readily photoswitch in solution, their photoisomerization in densely packed self-assembled monolayers (SAMs) can be suppressed. Reasons for this can be steric hindrance and/or electronic quenching, e.g., by exciton coupling. We address these possibilities by means of nonadiabatic molecular dynamics with trajectory surface hopping calculations, investigating the trans → cis isomerization of azobenzene after excitation into the ππ* absorption band. We consider a free monomer, an isolated dimer and a dimer embedded in a SAM-like environment of additional azobenzene molecules, imitating in this way the gradual transition from an unconstrained over an electronically coupled to an electronically coupled and sterically hindered, molecular switch. Our simulations reveal that in comparison to the single molecule the quantum yield of the trans → cis photoisomerization is similar for the isolated dimer, but greatly reduced in the sterically constrained situation. Other implications of dimerization and steric constraints are also discussed.