Molecular dynamics simulation with an ab initio potential energy function and finite element interpolation: The photoisomerization of cis-stilbene in solution

Approximating Periodic Potential Energy Surfaces with Sparse Trigonometric Interpolation

The potential energy surface (PES) describes the energy of a chemical system as a function of its geometry and is a fundamental concept in computational chemistry. A PES provides much useful information about the system, including the structures and energies of various stationary points, such as local minima, maxima, and transition states. Construction of full-dimensional PESs for molecules with more than 10 atoms is computationally expensive and often not feasible. Previous work in our group used sparse interpolation with polynomial basis functions to construct a surrogate reduced-dimensional PESs along chemically significant reaction coordinates, such as bond lengths, bond angles, and torsion angles. However, polynomial interpolation does not preserve the periodicity of the PES gradient with respect to angular components of geometry, such as torsion angles, which can lead to nonphysical phenomena. In this work, we construct a surrogate PES using trigonometric basis functions, for a system where the selected reaction coordinates all correspond to the torsion angles, resulting in a periodically repeating PES. We find that a trigonometric interpolation basis not only guarantees periodicity of the gradient but also results in slightly lower approximation error than polynomial interpolation.

Molecular Dynamics Simulations on Relaxed Reduced-Dimensional Potential Energy Surfaces

Molecular dynamics (MD) simulations with full-dimensional potential energy surfaces (PESs) obtained from high-level ab initio calculations are frequently used to model reaction dynamics of small molecules (i.e., molecules with up to 10 atoms). Construction of full-dimensional PESs for larger molecules is, however, not feasible since the number of ab initio calculations required grows rapidly with the increase of dimension. Only a small number of coordinates are often essential for describing the reactivity of even very large systems, and reduced-dimensional PESs with these coordinates can be built for reaction dynamics studies. While analytical methods based on transition-state theory framework are well established for analyzing the reduced-dimensional PESs, MD simulation algorithms capable of generating trajectories on such surfaces are more rare. In this work, we present a new MD implementation that utilizes the relaxed reduced-dimensional PES for standard microcanonical (NVE) and canonical (NVT) MD simulations. The method is applied to the pyramidal inversion of a NH₃ molecule. The results from the MD simulations on a reduced, three-dimensional PES are validated against the ab initio MD simulations, as well as MD simulations on full-dimensional PES and experimental data.

Landscapes of four-enantiomer conical intersections for photoisomerization of stilbene: CASSCF calculation

The photoisomerization of cis- and trans-stilbene through conical intersections (CI) is mainly governed by four dihedral angles around central C═C double bonds. The two of them are C-C═C-C and H-C═C-H dihedral angles that are found to form a mirror rotation coordinate, and the mirror plane appears at the two dihedral angles equal to zeros with which the middle state is defined through partial optimization. There exist the first-type of hula-twist-CI enantiomers, the second-type of hula-twist-CI enantiomers, the first-type of one-bond-flip-CI enantiomers, and the second type of one-bond-flip-CI enantiomers as well as cis-enantiomers and trans-enantiomers with respect to this mirror plane. The complete active space self-consistent field method is employed to calculate minimum potential energy profile along the mirror rotation coordinate for each enantiomers, and it is found that the left-hand manifold and the right-hand manifold of potential energy surfaces can be energetically transferred via photoisomerization. Furthermore, two-dimensional potential energy surfaces in terms of the branching plane g-h coordinates are constructed at vicinity of each conical intersection, and the landscapes of conical intersections show distinct feature, and in excited-state four potential wells separated in different section of g-h plane related to different conical intersections which indicate different photoisomerization pathways.

Photoisomerization of Stilbene: The Detailed XMCQDPT2 Treatment

We report the detailed XMCQDPT2/cc-pVTZ study of trans-cis photoisomerization in one of the core systems of both experimental and computational photochemistry-the stilbene molecule. For the first time, the potential energy surface (PES) of the S1 state has been directly optimized and scanned using a multistate multiconfiguration second-order perturbation theory. We characterize the trans-stilbene, pyramidalized (phantom), and DHP-cis-stilbene geometric domains of the S1 state and describe their stationary points including the transition states between them, as well as S1/S0 intersections. Also reported are the minima and the activation barriers in the ground state. Our calculations correctly predict the kinetic isotope effect due to H/D exchange at ethylenic hydrogens, the dynamic behavior of excited cis-stilbene, and trans-cis branching ratio after relaxation to S0 through a rather unsymmetric conical intersection. In general, the XMCQDPT2 results confirm the qualitative adequacy of the TDDFT (especially SF-TDDFT) picture of the excited stilbene but also reveal quantitative discrepancies that deserve further exploration.

Photoisomerization of stilbene: a spin-flip density functional theory approach

The photoisomerization process of 1,2-diphenylethylene (stilbene) is investigated using the spin-flip density functional theory (SFDFT), which has recently been shown to be a promising approach for locating conical intersection (CI) points (Minezawa, N.; Gordon, M. S. J. Phys. Chem. A2009, 113, 12749). The SFDFT method gives valuable insight into twisted stilbene to which the linear response time-dependent DFT approach cannot be applied. In contrast to the previous SFDFT study of ethylene, a distinct twisted minimum is found for stilbene. The optimized structure has a sizable pyramidalization angle and strong ionic character, indicating that a purely twisted geometry is not a true minimum. In addition, the SFDFT approach can successfully locate two CI points: the twisted-pyramidalized CI that is similar to the ethylene counterpart and another CI that possibly lies on the cyclization pathway of cis-stilbene. The mechanisms of the cis--trans isomerization reaction are discussed on the basis of the two-dimensional potential energy surface along the twisting and pyramidalization angles.

Photoinduced conformational dynamics of a photoswitchable peptide: a nonequilibrium molecular dynamics simulation study

Employing nonequilibrium molecular dynamics simulations, a comprehensive computational study of the photoinduced conformational dynamics of a photoswitchable bicyclic azobenzene octapeptide is presented. The calculation of time-dependent probability distributions along various global and local reaction coordinates reveals that the conformational rearrangement of the peptide is rather complex and occurs on at least four timescales: 1) After photoexcitation, the azobenzene unit of the molecule undergoes nonadiabatic photoisomerization within 0.2 ps. 2) On the picosecond timescale, the cooling (13 ps) and the stretching (14 ps) of the photoexcited peptide is observed. 3) Most reaction coordinates exhibit a 50-100 ps component reflecting a fast conformational rearrangement. 4) The 500-1000 ps component observed in the simulation accounts for the slow diffusion-controlled conformational equilibration of the system. The simulation of the photoinduced molecular processes is in remarkable agreement with time-resolved optical and infrared experiments, although the calculated cooling as well as the initial conformational rearrangements of the peptide appear to be somewhat too slow. Based on an ab initio parameterized vibrational Hamiltonian, the time-dependent amide I frequency shift is calculated. Both intramolecular and solvent-induced contributions to the frequency shift were found to change by < or = 2 cm(-1), in reasonable agreement with experiment. The potential of transient infrared spectra to characterize the conformational dynamics of peptides is discussed in some detail.