5 research outputs found
Which Electronic Structure Method to Choose in Trajectory Surface Hopping Dynamics Simulations? Azomethane as a Case Study
Nonadiabatic
dynamics simulations have become a standard approach
to explore photochemical reactions. Such simulations require underlying
potential energy surfaces and couplings between them, calculated at
a chosen level of theory, yet this aspect is rarely assessed. Here,
in combination with the popular trajectory surface hopping dynamics
method, we use a high-accuracy XMS-CASPT2 electronic structure level
as a benchmark for assessing the performances of various post-Hartree–Fock
methods (namely, CIS, ADC(2), CC2, and CASSCF) and exchange–correlation
functionals (PBE, PBE0, and CAM-B3LYP) in a TD-DFT/TDA context, using
the isomerization around a double bond as test case. Different relaxation
pathways are identified, and the ability of the different methods
to reproduce their relative importance and time scale is discussed.
The results show that multireference electronic structure methods
should be preferred, when studying nonadiabatic decay between excited
and ground states. If not affordable, TD-DFT with TDA and hybrid functionals
and ADC(2) are efficient alternatives but overestimate the nonradiative
decay yield and thus may miss deexcitation pathways
Geometric Rotation of the Nuclear Gradient at a Conical Intersection: Extension to Complex Rotation of Diabatic States
Nonadiabatic dynamics
in the vicinity of conical intersections is
of essential importance in photochemistry. It is well known that if
the branching space is represented in polar coordinates, then for
a geometry represented by angle θ, the corresponding adiabatic
states are obtained from the diabatic states with the mixing angle
θ/2. In an equivalent way, one can study the relation between
the real rotation of diabatic states and the resulting nuclear gradient.
In this work, we extend the concept to allow a complex rotation of
diabatic states to form a nonstationary superposition of electronic
states. Our main result is that this leads to an elliptical transformation
of the effective potential energy surfaces; i.e., the magnitude of
the initial nuclear gradient changes as well as its direction. We
fully explore gradient changes that result from varying both θ
and ϕ (the complex rotation angle) as a way of electronically
controlling nuclear motion, through Ehrenfest dynamics simulations
for benzene cation
Electronic Control of Initial Nuclear Dynamics Adjacent to a Conical Intersection
Photoionization
can create a nonstationary electronic state and
therefore initiates coupled electron–nuclear dynamics in molecules.
Using a CASSCF implementation of the Ehrenfest method, we study the
nuclear dynamics following vertical ionization of toluene, starting
close to the conical intersection between ground and first excited
states of its cation. The results show how the initial nuclear dynamics
is controlled by the nonstationary electronic state character. In
particular, ionization of this system leading to an equal superposition
of the two lowest energy states can initiate nuclear dynamics in an
orthogonal direction in the branching space to dynamics on the ground
or first excited state potential energy surfaces alone
Charge migration engineered by localisation: electron-nuclear dynamics in polyenes and glycine
<p>We demonstrate that charge migration can be ‘engineered’ in arbitrary molecular systems if a single localised orbital – that diabatically follows nuclear displacements – is ionised. Specifically, we describe the use of natural bonding orbitals in Complete Active Space Configuration Interaction (CASCI) calculations to form cationic states with localised charge, providing consistently well-defined initial conditions across a zero point energy vibrational ensemble of molecular geometries. In Ehrenfest dynamics simulations following localised ionisation of -electrons in model polyenes (hexatriene and decapentaene) and -electrons in glycine, oscillatory charge migration can be observed for several femtoseconds before dephasing. Including nuclear motion leads to slower dephasing compared to fixed-geometry electron-only dynamics results. For future work, we discuss the possibility of designing laser pulses that would lead to charge migration that is experimentally observable, based on the proposed diabatic orbital approach.</p
OpenMolcas: From source code to insight
In this article we describe the OpenMolcas environment and invite the computational chemistry community to collaborate. The open-source project already
includes a large number of new developments realized during the transition from
the commercial MOLCAS product to the open-source platform. The paper initially
describes the technical details of the new software development platform. This is followed by brief presentations of many new methods, implementations, and features
of the OpenMolcas program suite. These developments include novel wave function methods such as stochastic complete active space self-consistent field, density
matrix renormalization group (DMRG) methods, and hybrid multiconfigurational wave function and density functional theory models. Some of these implementations
include an array of additional options and functionalities. The paper proceeds and
describes developments related to explorations of potential energy surfaces. Here
we present methods for the optimization of conical intersections, the simulation of
adiabatic and nonadiabatic molecular dynamics and interfaces to tools for semiclassical and quantum mechanical nuclear dynamics. Furthermore, the article describes
features unique to simulations of spectroscopic and magnetic phenomena such as
the exact semiclassical description of the interaction between light and matter, various X-ray processes, magnetic circular dichroism and properties. Finally, the paper
describes a number of built-in and add-on features to support the OpenMolcas platform with post calculation analysis and visualization, a multiscale simulation option
using frozen-density embedding theory and new electronic and muonic basis sets