
Conical intersections (CIs) have been widely studied using spectroscopic
techniques. However, CIs have mainly been identified by rapid internal
conversion transitions that take place after the photoexcitation. Such
identifications cannot distinguish various types of intersections as well as to
separate the actual intersection from an avoided crossing. In this paper, we
investigate how ultrafast IR laser pulses can be utilized to stimulate nuclear
dynamics revealing geometric phase features associated with CIs. We consider
two lowdimensional nonadiabatic models to obtain optimal two and threepulse
laser sequences for stimulating nuclear dynamics necessary for the CI
identification. Our results provide insights on designing nonlinear
spectroscopic schemes for subsequent probes of the nuclear wavepackets by
ultrafast electron diffraction techniques to unambiguously detect CIs in
molecules.

A new methodology of simulating nonadiabatic dynamics using frozenwidth
Gaussian wavepackets within the moving crude adiabatic representation with the
onthefly evaluation of electronic structure is presented. The main feature of
the new approach is elimination of any global or local model representation of
electronic potential energy surfaces, instead, the electronnuclear interaction
is treated explicitly using the Gaussian integration. As a result, the new
scheme does not introduce any uncontrolled approximations. The employed
variational principle ensures the energy conservation and leaves the number of
electronic and nuclear basis functions as the only parameter determining the
accuracy. To assess performance of the approach, a model with two electronic
and two nuclear spacial degrees of freedom containing conical intersections
between potential energy surfaces has been considered. Dynamical features
associated with nonadiabatic transitions and nontrivial geometric (or Berry)
phases were successfully reproduced within a limited basis expansion.

Bound states in the continuum (BSCs) were reported in a linear vibronic
coupling model with a conical intersection (CI) [Cederbaum et al. Phys. Rev.
Lett. 90, 013001 (2003)]. It was also found that these states are destroyed
within the BornOppenheimer approximation (BOA). We investigate whether a
nontrivial topological or geometric phase (GP) associated with the CI is
responsible for BSCs. To address this question we explore modifications of the
original twodimensional twostate linear vibronic coupling model supporting
BSCs. These modifications either add GP effects after the BOA or remove the GP
within a twostate problem. Using the stabilization graph technique we shown
that the GP is crucial for emergence of BSCs.

We investigate geometric phase (GP) effects in nonadiabatic transitions
through a conical intersection (CI) in an Ndimensional linear vibronic
coupling (NDLVC) model. This model allows for the coordinate transformation
encompassing all nonadiabatic effects within a twodimensional (2D) subsystem
while the other N2 dimensions form a system of uncoupled harmonic oscillators
identical for both electronic states and coupled bilinearly with the subsystem
coordinates. The 2D subsystem governs ultrafast nonadiabatic dynamics through
the CI and provides a convenient model for studying GP effects. Parameters of
the original NDLVC model define the Hamiltonian of the transformed 2D
subsystem and thus influence GP effects directly. Our analysis reveals what
values of NDLVC parameters can introduce symmetry breaking in the 2D subsystem
that diminishes GP effects.

We propose a systematic approach to the basis set extension for nonadiabatic
dynamics of entangled combination of nuclear coherent states (CSs) evolving
according to the timedependent variational principle (TDVP). TDVP provides a
rigorous framework for fully quantum nonadiabatic dynamics of closed systems,
however, quality of results strongly depends on available basis functions.
Starting with a single nuclear CS replicated vertically on all electronic
states, our approach clones this function when replicas of the CS on different
electronic states experience increasingly different forces. Created clones move
away from each other (decohere) extending the basis set. To determine a moment
for cloning we introduce generalized forces based on derivatives that maximally
contribute to a variation of the total quantum action and thus account for
entanglement of all basis functions.

An accurate description of nonadiabatic dynamics of molecular species on
metallic surfaces poses a serious computational challenge associated with a
multitude of closelyspaced electronic states. We propose a mixed
quantumclassical scheme that addresses this challenge by introducing
collective electronic variables. These variables are defined through analytic
blockdiagonalization applied to the timedependent Hamiltonian matrix
governing the electronic dynamics. We compare our scheme with the Ehrenfest
approach and with a fullmemory electronic friction model on a onedimensional
"adatom + atomic chain" model. Our simulations demonstrate that collectivemode
dynamics with only few (23) electronic variables is robust and can describe a
variety of situations: from a chemisorbed atom on an insulator to an atom on a
metallic surface. Our molecular model also reveals that the friction approach
is prone to unpredictable and catastrophic failures.

Onthefly quantum nonadiabatic dynamics for large systems greatly benefits
from the adiabatic representation readily available from the electronic
structure programs. However, frequently occurring in this representation
conical intersections introduce nontrivial geometric or Berry phases which
require a special treatment for adequate modelling of the nuclear dynamics. We
analyze two approaches for nonadiabatic dynamics using the timedependent
variational principle and the adiabatic representation. The first approach
employs adiabatic electronic functions with global parametric dependence on the
nuclear coordinates. The second approach uses adiabatic electronic functions
obtained only at the centres of moving localized nuclear basis functions (e.g.
frozenwidth Gaussians). Unless a gauge transformation is used to enforce
singlevalued boundary conditions, the first approach fails to capture the
geometric phase. In contrast, the second approach accounts for the geometric
phase naturally because of the absence of the global nuclear coordinate
dependence in the electronic functions.

Most of the textbooks explaining electric conductivity in the context of
quantum mechanics provide either incomplete or semiclassical explanations that
are not connected with the elementary concepts of quantum mechanics. We
illustrate the conduction phenomena using the simplest model system in quantum
dynamics, a particle in a box (PIB). To induce the particle dynamics, a linear
potential tilting the bottom of the box is introduced, which is equivalent to
imposing a constant electric field for a charged particle. Although the PIB
model represents a closed system that cannot have a flow of electrons through
the system, we consider the oscillatory dynamics of the particle probability
density as the analogue of the electric current. Relating the amplitude and
other parameters of the particle oscillatory dynamics with the gap between the
ground and excited states of the PIB model allows us to demonstrate one of the
most basic dependencies of electric conductivity on the valenceconduction band
gap of the material.

We extend the localized operator partitioning method (LOPM) [J. Nagesh, A.F.
Izmaylov, and P. Brumer, J. Chem. Phys. 142, 084114 (2015)] to the
timedependent density functional theory (TDDFT) framework to partition
molecular electronic energies of excited states in a rigorous manner. A
molecular fragment is defined as a collection of atoms using
StratmanScuseriaFrisch atomic partitioning. A numerically efficient scheme
for evaluating the fragment excitation energy is derived employing a resolution
of the identity to preserve standard one and twoelectron integrals in the
final expressions. The utility of this partitioning approach is demonstrated by
examining several excited states of two bichromophoric compounds:
9((1naphthyl)methyl)anthracene and 4((2naphthyl)methyl)benzaldehyde.
The LOPM is found to provide nontrivial insights into the nature of electronic
energy localization that are not accessible using simple density difference
analysis.

Electronic wavefunctions in the adiabatic representation acquire nontrivial
geometric phases (GPs) when corresponding potential energy surfaces undergo
conical intersection (CI). These GPs have profound effects on the nuclear
quantum dynamics and cannot be eliminated in the adiabatic representation
without changing the physics of the system. To define dynamical effects arising
from the GP presence the nuclear quantum dynamics of the CI containing system
is compared with that of the system with artificially removed GP. We explore a
new construction of the system with removed GP via a modification of the
diabatic representation for the original CI containing system. Using an
absolute value function of diabatic couplings we remove the GP while preserving
adiabatic potential energy surfaces and CI. We assess GP effects in dynamics of
a twodimensional linear vibronic coupling model both for ground and excited
state dynamics. Results are compared with those obtained with a conventional
removal of the GP by ignoring doublevalued boundary conditions of the real
electronic wavefunctions. Interestingly, GP effects appear similar in two
approaches only for the low energy dynamics. In contrast with the conventional
approach, a new approach does not have substantial GP effects in the ultrafast
excited state dynamics.

The role of electronnuclear entanglement on the validity of the
BornOppenheimer (BO) approximation is investigated. While nonadiabatic
couplings generally lead to entanglement and to a failure of the BO
approximation, surprisingly the degree of electronnuclear entanglement is
found to be uncorrelated with the degree of validity of the BO approximation.
This is because while the degree of entanglement of BO states is determined by
their deviation from the corresponding states in the crude BO approximation,
the accuracy of the BO approximation is dictated, instead, by the deviation of
the BO states from the exact electronnuclear states. In fact, in the context
of a minimal avoided crossing model, extreme cases are identified where an
adequate BO state is seen to be maximally entangled, and where the BO
approximation fails but the associated BO state remains approximately
unentangled. Further, the BO states are found to not preserve the entanglement
properties of the exact electronnuclear eigenstates, and to be completely
unentangled only in the limit in which the BO approximation becomes exact.

The diagonal BornOppenheimer correction (DBOC) stems from the diagonal
second derivative coupling term in the adiabatic representation, and it can
have an arbitrary large magnitude when a gap between neighbouring
BornOppenheimer (BO) potential energy surfaces (PESs) is closing.
Nevertheless, DBOC is typically neglected in mixed quantumclassical methods of
simulating nonadiabatic dynamics (e.g., fewestswitch surface hopping (FSSH)
method). A straightforward addition of DBOC to BO PESs in the FSSH method,
FSSH+D, has been shown to lead to numerically much inferior results for models
containing conical intersections. More sophisticated variation of the DBOC
inclusion, phasespace surfacehopping (PSSH) was more successful than FSSH+D
but on model problems without conical intersections. This work comprehensively
assesses the role of DBOC in nonadiabatic dynamics of two electronic state
problems and the performance of FSSH, FSSH+D, and PSSH methods in variety of
one and twodimensional models. Our results show that the inclusion of DBOC
can enhance the accuracy of surface hopping simulations when two conditions are
simultaneously satisfied: 1) nuclei have kinetic energy lower than DBOC and 2)
PESs are not strongly nonadiabatically coupled. The inclusion of DBOC is
detrimental in situations where its energy scale becomes very high or even
diverges, because in these regions PESs are also very strongly coupled. In this
case, the true quantum formalism heavily relies on an interplay between
diagonal and offdiagonal nonadiabatic couplings while surface hopping
approaches treat diagonal terms as PESs and offdiagonal ones stochastically.

Photolyases are proteins capable of harvesting the sunlight to repair DNA
damages caused by UV light. In this work we focus on the first step in the
repair process of the cyclobutane pyrimidine dimer photoproduct (CPD) lesion,
which is an electron transfer (ET) from a flavine cofactor to CPD, and study
the role of various nuclear degrees of freedom (DOF) in this step. The ET step
has been experimentally studied using transient spectroscopy and the
corresponding data provide excellent basis for testing the quality of quantum
dynamical models. Based on previous theoretical studies of electronic structure
and conformations of the protein active site, we present a procedure to build a
diabatic Hamiltonian for simulating the ET reaction in a molecular complex
mimicking the enzyme's active site. We generate a reduced nuclear dimensional
model that provides a first nonempirical quantum dynamical description of the
structural features influencing the ET rate. By varying the nuclear DOF
parametrization in the model to assess the role of different nuclear motions,
we demonstrate that the low frequency flavin butterfly bending mode slows ET by
reducing FranckCondon overlaps between donor and acceptor states and also
induces decoherence.

We have developed a numerical differentiation scheme which eliminates
evaluation of overlap determinants in calculating the timederivative
nonadiabatic couplings (TDNACs). Evaluation of these determinants was the
bottleneck in previous implementations of mixed quantumclassical methods using
numerical differentiation of electronic wave functions in the
Slaterdeterminant representation. The central idea of our approach is, first,
to reduce the analytic time derivatives of Slater determinants to time
derivatives of molecular orbitals, and then to apply a finitedifference
formula. Benchmark calculations prove the efficiency of the proposed scheme
showing impressive severalorderofmagnitude speedups of the TDNAC calculation
step for midsize molecules.

Methods of quantum nuclear wavefunction dynamics have become very efficient
in simulating large isolated systems using the timedependent variational
principle (TDVP). However, a straightforward extension of the TDVP to the
density matrix framework gives rise to methods that do not conserve the energy
in the isolated system limit and the total system population for open systems
where only energy exchange with the environment is allowed. These problems
arise when the system density is in a mixed state and is simulated using an
incomplete basis. Thus, the basis set incompleteness, which is inevitable in
practical calculations, creates artificial channels for energy and population
dissipation. To overcome this unphysical behavior, we have introduced a
constrained Lagrangian formulation of TDVP applied to the nonstochastic open
system Schrodinger equation (NOSSE) [L. JoubertDoriol, I. G. Ryabinkin, and A.
F. Izmaylov, J. Chem. Phys. 141, 234112 (2014)]. While our formulation can be
applied to any variational ansatz for the system density matrix, derivation of
working equations and numerical assessment are done within the variational
multiconfiguration Gaussian approach for a twodimensional linear vibronic
coupling model system interacting with a harmonic bath.

Adequate simulation of nonadiabatic dynamics through conical intersection
requires account for a nontrivial geometric phase (GP) emerging in electronic
and nuclear wavefunctions in the adiabatic representation. Popular mixed
quantumclassical (MQC) methods, surface hopping and Ehrenfest, do not carry a
nuclear wavefunction to be able to incorporate the GP into nuclear dynamics.
Surprisingly, the MQC methods reproduce ultrafast interstate crossing dynamics
generated with the exact quantum propagation so well as if they contained
information about the GP. Using twodimensional linear vibronic coupling models
we unravel how the MQC methods can effectively mimic the most significant
dynamical GP effects: 1) compensation for repulsive diagonal second order
nonadiabatic couplings and 2) transfer enhancement for a fully cylindrically
symmetric component of a nuclear distribution.

The localized operator partitioning method [Y. Khan and P. Brumer, J. Chem.
Phys. 137, 194112 (2012)] rigorously defines the electronic energy on any
subsystem within a molecule and gives a precise meaning to the subsystem ground
and excited electronic energies, which is crucial for investigating electronic
energy transfer from first principles. However, an efficient implementation of
this approach has been hindered by complicated one and twoelectron integrals
arising in its formulation. Using a resolution of the identity in the
definition of partitioning we reformulate the method in a computationally
efficient manner that involves standard one and twoelectron integrals. We
apply the developed algorithm to the 9((1naphthyl)methyl)anthracene (A1N)
molecule by partitioning A1N into anthracenyl and CH2naphthyl groups as
subsystems, and examine their electronic energies and populations for several
excited states using Configuration Interaction Singles method. The implemented
approach shows a wide variety of different behaviors amongst the excited
electronic states.

We propose an extension of the Schr\"odinger equation for a quantum system
interacting with environment. This equation describes dynamics of auxiliary
wavefunctions $\mathbf{m}$, from which the system density matrix can be
reconstructed as $\hat{\rho} = \mathbf{m} \mathbf{m}^\dagger$. We formulate a
compatibility condition, which ensures that the reconstructed density satisfies
a given quantum master equation for the system density. The resulting
nonstochastic evolution equation preserves positivedefiniteness of the system
density and is applicable to both Markovian and nonMarkovian systembath
treatments. Our formalism also resolves a longstanding problem of energy
nonconservation in the timedependent variational principle applied to mixed
states of closed systems.

We consider a fully quadratic vibronic model Hamiltonian for studying
photoinduced electronic transitions through conical intersections. Using a
second order perturbative approximation for diabatic couplings we derive an
analytical expression for the time evolution of electronic populations at a
given temperature. This formalism extends upon a previously developed
perturbative technique for a linear vibronic coupling Hamiltonian. The
advantage of the quadratic model Hamiltonian is that it allows one to use
separate quadratic representations for potential energy surfaces of different
electronic states and a more flexible representation of interstate couplings.
We explore features introduced by the quadratic Hamiltonian in a series of 2D
models, and then apply our formalism to the 2,6bis(methylene) adamantyl
cation, and its dimethyl derivative. The Hamiltonian parameters for the
molecular systems have been obtained from electronic structure calculations
followed by a diabatization procedure. The evolution of electronic populations
in the molecular systems using the perturbative formalism shows a good
agreement with that from variational quantum dynamics.

We investigate the role of the geometric phase (GP) in an internal conversion
process when the system changes its electronic state by passing through a
conical intersection (CI). Local analysis of a twodimensional linear vibronic
coupling (LVC) model Hamiltonian near the CI shows that the role of the GP is
twofold. First, it compensates for a repulsion created by the socalled
diagonal BornOppenheimer correction (DBOC). Second, the GP enhances the
nonadiabatic transition probability for a wavepacket part that experiences a
central collision with the CI. To assess the significance of both GP
contributions we propose two indicators that can be computed from parameters of
electronic surfaces and initial conditions. To generalize our analysis to
Ndimensional systems we introduce a reduction of a general Ndimensional LVC
model to an effective 2D LVC model using a mode transformation that preserves
shorttime dynamics of the original Ndimensional model. Using examples of the
bis(methylene) adamantyl and butatriene cations, and the pyrazine molecule we
have demonstrated that their effective 2D models reproduce the shorttime
dynamics of the corresponding full dimensional models, and the introduced
indicators are very reliable in assessing GP effects.

An efficient computational algorithm to implement a local operator approach
to partitioning electronic energy in general molecular systems is presented.
This approach, which rigorously defines the electronic energy on any subsystem
within a molecule, gives a precise meaning to the subsystem ground and excited
electronic energies, which is crucial for investigating electronic energy
transfer from first principles. We apply the technique to the
$9$(($1$naphthyl)$$methyl)anthracene (A1N) molecule by partitioning A1N
into anthracenyl and CH$_2$naphthyl groups as subsystems, and examine their
electronic energies and populations for several excited states using
Configuration Interaction Singles method. The implemented approach shows a wide
variety of different behaviors amongst these excited electronic states.

We analyze two approaches to the quantumclassical Liouville (QCL) formalism
that differ in the order of two operations: Wigner transformation and
projection onto adiabatic electronic states. The analysis is carried out on a
twodimensional linear vibronic model where geometric phase (GP) effects
arising from a conical intersection profoundly affect nuclear dynamics. We find
that the WignerthenAdiabatic (WA) QCL approach captures GP effects, whereas
the AdiabaticthenWigner (AW) QCL approach does not. Moreover, the Wigner
transform in AWQCL leads to an illdefined Fourier transform of doublevalued
functions. The doublevalued character of these functions stems from the
nontrivial GP of adiabatic electronic states in the presence of a conical
intersection. In contrast, WAQCL avoids this issue by starting with the Wigner
transform of singlevalued quantities of the full problem. Since the WAQCL
approach uses solely the adiabatic potentials and nonadiabatic derivative
couplings as an input, our results indicate that WAQCL can capture GP effects
in general twostate crossing problems using firstprinciples electronic
structure calculations without prior diabatization or introduction of explicit
phase factors.

In molecular systems containing conical intersections (CIs), a nontrivial
geometric phase (GP) appears in the nuclear and electronic wavefunctions in
the adiabatic representation. We study GP effects in nuclear dynamics of an
Ndimensional linear vibronic coupling (LVC) model. The main impact of GP on
lowenergy nuclear dynamics is reduction of population transfer between the
local minima of the LVC lower energy surface. For the LVC model, we proposed an
isometric coordinate transformation that confines nonadiabatic effects within
a twodimensional subsystem interacting with an N2 dimensional environment.
Since environmental modes do not couple electronic states, all GP effects
originate from nuclear dynamics within the subsystem. We explored when the GP
affects nuclear dynamics of the isolated subsystem, and how the
subsystemenvironment interaction can interfere with GP effects. Comparing
quantum dynamics with and without GP allowed us to devise simple rules to
determine significance of the GP for nuclear dynamics in this model.

We show that finite systems with conical intersections can exhibit
spontaneous symmetry breaking which manifests itself in spatial localization of
eigenstates. This localization has a geometric phase origin and is robust
against variation of model parameters. The transition between localized and
delocalized eigenstate regimes resembles a continuous phase transition. The
localization slows down the lowenergy quantum nuclear dynamics at zero and low
temperatures.

We report a method for the efficient evaluation of analytic infrared (IR)
intensities within generalized KohnSham density functional theory using
Gaussian orbitals and periodic boundary conditions. A discretized form of the
Berry phase is used to evaluate a periodic dipole moment and its derivatives
with respect to inphase nuclear coordinate displacements. Benchmark
calculations are presented for onedimensional chains of water molecules and
poly(paraphenylenevinylene).