• 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 low-dimensional nonadiabatic models to obtain optimal two- and three-pulse laser sequences for stimulating nuclear dynamics necessary for the CI identification. Our results provide insights on designing non-linear 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 frozen-width Gaussian wavepackets within the moving crude adiabatic representation with the on-the-fly 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 electron-nuclear 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 Born-Oppenheimer 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 two-dimensional two-state linear vibronic coupling model supporting BSCs. These modifications either add GP effects after the BOA or remove the GP within a two-state 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 N-dimensional linear vibronic coupling (ND-LVC) model. This model allows for the coordinate transformation encompassing all nonadiabatic effects within a two-dimensional (2D) subsystem while the other N-2 dimensions form a system of uncoupled harmonic oscillators identical for both electronic states and coupled bi-linearly with the subsystem coordinates. The 2D subsystem governs ultra-fast nonadiabatic dynamics through the CI and provides a convenient model for studying GP effects. Parameters of the original ND-LVC model define the Hamiltonian of the transformed 2D subsystem and thus influence GP effects directly. Our analysis reveals what values of ND-LVC 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 time-dependent 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 closely-spaced electronic states. We propose a mixed quantum-classical scheme that addresses this challenge by introducing collective electronic variables. These variables are defined through analytic block-diagonalization applied to the time-dependent Hamiltonian matrix governing the electronic dynamics. We compare our scheme with the Ehrenfest approach and with a full-memory electronic friction model on a one-dimensional "adatom + atomic chain" model. Our simulations demonstrate that collective-mode dynamics with only few (2-3) 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.
  • On-the-fly 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 non-trivial 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 time-dependent 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. frozen-width Gaussians). Unless a gauge transformation is used to enforce single-valued 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 semi-classical 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 valence-conduction 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 time-dependent density functional theory (TD-DFT) framework to partition molecular electronic energies of excited states in a rigorous manner. A molecular fragment is defined as a collection of atoms using Stratman-Scuseria-Frisch 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 two-electron integrals in the final expressions. The utility of this partitioning approach is demonstrated by examining several excited states of two bichromophoric compounds: 9-((1-naphthyl)-methyl)-anthracene and 4-((2-naphthyl)-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 wave-functions 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 two-dimensional 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 double-valued boundary conditions of the real electronic wave-functions. 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 ultra-fast excited state dynamics.
  • The role of electron-nuclear entanglement on the validity of the Born-Oppenheimer (BO) approximation is investigated. While nonadiabatic couplings generally lead to entanglement and to a failure of the BO approximation, surprisingly the degree of electron-nuclear 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 electron-nuclear 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 electron-nuclear eigenstates, and to be completely unentangled only in the limit in which the BO approximation becomes exact.
  • The diagonal Born-Oppenheimer 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 Born-Oppenheimer (BO) potential energy surfaces (PESs) is closing. Nevertheless, DBOC is typically neglected in mixed quantum-classical methods of simulating nonadiabatic dynamics (e.g., fewest-switch 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, phase-space surface-hopping (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 two-dimensional 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 off-diagonal nonadiabatic couplings while surface hopping approaches treat diagonal terms as PESs and off-diagonal 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 non-empirical 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 Franck-Condon 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 time-derivative non-adiabatic couplings (TDNACs). Evaluation of these determinants was the bottleneck in previous implementations of mixed quantum-classical methods using numerical differentiation of electronic wave functions in the Slater-determinant 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 finite-difference formula. Benchmark calculations prove the efficiency of the proposed scheme showing impressive several-order-of-magnitude speedups of the TDNAC calculation step for midsize molecules.
  • Methods of quantum nuclear wave-function dynamics have become very efficient in simulating large isolated systems using the time-dependent 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 non-stochastic open system Schrodinger equation (NOSSE) [L. Joubert-Doriol, 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 two-dimensional linear vibronic coupling model system interacting with a harmonic bath.
  • Adequate simulation of non-adiabatic dynamics through conical intersection requires account for a non-trivial geometric phase (GP) emerging in electronic and nuclear wave-functions in the adiabatic representation. Popular mixed quantum-classical (MQC) methods, surface hopping and Ehrenfest, do not carry a nuclear wave-function to be able to incorporate the GP into nuclear dynamics. Surprisingly, the MQC methods reproduce ultra-fast interstate crossing dynamics generated with the exact quantum propagation so well as if they contained information about the GP. Using two-dimensional 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 non-adiabatic 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 two-electron 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 two-electron integrals. We apply the developed algorithm to the 9-((1-naphthyl)-methyl)-anthracene (A1N) molecule by partitioning A1N into anthracenyl and CH2-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 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 wave-functions $\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 non-stochastic evolution equation preserves positive-definiteness of the system density and is applicable to both Markovian and non-Markovian system-bath treatments. Our formalism also resolves a long-standing problem of energy non-conservation in the time-dependent 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,6-bis(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 two-dimensional 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 so-called diagonal Born-Oppenheimer correction (DBOC). Second, the GP enhances the non-adiabatic transition probability for a wave-packet 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 N-dimensional systems we introduce a reduction of a general N-dimensional LVC model to an effective 2D LVC model using a mode transformation that preserves short-time dynamics of the original N-dimensional 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 short-time 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 quantum-classical 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 two-dimensional linear vibronic model where geometric phase (GP) effects arising from a conical intersection profoundly affect nuclear dynamics. We find that the Wigner-then-Adiabatic (WA) QCL approach captures GP effects, whereas the Adiabatic-then-Wigner (AW) QCL approach does not. Moreover, the Wigner transform in AW-QCL leads to an ill-defined Fourier transform of double-valued functions. The double-valued character of these functions stems from the nontrivial GP of adiabatic electronic states in the presence of a conical intersection. In contrast, WA-QCL avoids this issue by starting with the Wigner transform of single-valued quantities of the full problem. Since the WA-QCL approach uses solely the adiabatic potentials and non-adiabatic derivative couplings as an input, our results indicate that WA-QCL can capture GP effects in general two-state crossing problems using first-principles 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 wave-functions in the adiabatic representation. We study GP effects in nuclear dynamics of an N-dimensional linear vibronic coupling (LVC) model. The main impact of GP on low-energy 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 non-adiabatic effects within a two-dimensional subsystem interacting with an N-2 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 subsystem-environment 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 low-energy quantum nuclear dynamics at zero and low temperatures.
  • We report a method for the efficient evaluation of analytic infrared (IR) intensities within generalized Kohn-Sham 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 in-phase nuclear coordinate displacements. Benchmark calculations are presented for one-dimensional chains of water molecules and poly(paraphenylenevinylene).