Search

Quantum dynamics simulation and computational spectroscopy serve as indispensable tools for the theoretical understanding of various fundamental physical and chemical processes, ranging from charge transfer to photochemical reactions. When simulating realistic systems, the primary challenge stems from the overwhelming number of degrees of freedom and the pronounced many-body correlations. Here, we present Kylin-V, an innovative quantum dynamics package designed for accurate and efficient simulations of dynamics and spectroscopic properties of vibronic Hamiltonians for molecular systems and their aggregates. Kylin-V supports various quantum dynamics and computational spectroscopy methods, such as time-dependent density matrix renormalization group and our recently proposed single-site and hierarchical mapping approaches, as well as vibrational heat-bath configuration interaction. In this paper, we introduce the methodologies implemented in Kylin-V and illustrate their performances through a diverse collection of numerical examples. [ABSTRACT FROM AUTHOR]

Fewest-switches surface hopping (FSSH) has emerged as one of the leading methods for modeling the quantum dynamics of molecular systems. While its original formulation was limited to adiabatic populations, the growing interest in the application of FSSH to coherent phenomena prompts the question of how one should construct a complete density matrix based on FSSH trajectories. A straightforward solution is to define adiabatic coherences based on wavefunction coefficients. In this paper, we demonstrate that inconsistencies introduced in the density matrix through such treatment may lead to a violation of positivity. We furthermore show that a recently proposed coherent generalization of FSSH results in density matrices that satisfy positivity while yielding improved accuracy throughout much (but not all) of parameter space. [ABSTRACT FROM AUTHOR]

The generalized quantum master equation (GQME) provides a general and formally exact framework for simulating the reduced dynamics of open quantum systems. The recently introduced modified approach to the GQME (M-GQME) corresponds to a specific implementation of the GQME that is geared toward simulating the dynamics of the electronic reduced density matrix in systems governed by an excitonic Hamiltonian. Such a Hamiltonian, which is often used for describing energy and charge transfer dynamics in complex molecular systems, is given in terms of diabatic electronic states that are coupled to each other and correspond to different nuclear Hamiltonians. Within the M-GQME approach, the effect of the nuclear degrees of freedom on the time evolution of the electronic density matrix is fully captured by a memory kernel superoperator, which can be obtained from short-lived (compared to the time scale of energy/charge transfer) projection-free inputs. In this paper, we test the ability of the M-GQME to predict the energy transfer dynamics within a seven-state benchmark model of the Fenna–Matthews–Olson (FMO) complex, with the short-lived projection-free inputs obtained via the Ehrenfest method. The M-GQME with Ehrenfest-based inputs is shown to yield accurate results across a wide parameter range. It is also found to dramatically outperform the direct application of the Ehrenfest method and to provide better-behaved convergence with respect to memory time in comparison to an alternative implementation of the GQME approach previously applied to the same FMO model. [ABSTRACT FROM AUTHOR]

In this paper, we investigate the time evolution of quantum coherence—the off-diagonal elements of the density matrix of a multistate quantum system—from the perspective of the Wigner–Moyal formalism. This approach provides an exact phase space representation of quantum mechanics. We consider the coherent evolution of nuclear wavepackets in a molecule with two electronic states. For harmonic potentials, the problem is analytically soluble for both a fully quantum mechanical description and a semiclassical description. We highlight the serious deficiencies of the semiclassical treatment of coherence for general systems and illustrate how even qualitative accuracy requires higher order terms in the Moyal expansion to be included. The model provides an experimentally relevant example of a molecular Schrödinger's cat state. The alive and dead cats of the exact two-state quantum evolution collapse into a "zombie" cat in the semiclassical limit—an averaged behavior, neither alive nor dead, leading to significant errors. The inclusion of the Moyal correction restores a faithful simultaneously alive and dead representation of the cat that is experimentally observable. [ABSTRACT FROM AUTHOR]

The time dependent density matrix renormalization group (TD-DMRG) has become one of the cutting edge methods of quantum dynamics for complex systems. In this paper, we comparatively study the accuracy of three time evolution schemes in the TD-DMRG, the global propagation and compression method with the Runge-Kutta algorithm (P&C-RK), the time dependent variational principle based methods with the matrix unfolding algorithm (TDVP-MU), and with the projector-splitting algorithm (TDVP-PS), by performing benchmarks on the exciton dynamics of the Fenna-Matthews-Olson complex. We show that TDVP-MU and TDVP-PS yield the same result when the time step size is converged and they are more accurate than P&C-RK4, while TDVP-PS tolerates a larger time step size than TDVP-MU. We further adopt the graphical processing units to accelerate the heavy tensor contractions in the TD-DMRG, and it is able to speed up the TDVP-MU and TDVP-PS schemes by up to 73 times. [ABSTRACT FROM AUTHOR]

The molecular dynamics following the electronic B Π u 3 0 + ⟵ X Σ g + 1 photoexcitation of the iodine molecule embedded in solid krypton are studied quantum mechanically using the Gaussian variant of the multiconfigurational time-dependent Hartree method (G-MCTDH). The accuracy of the Gaussian wave packet approximation is validated against numerically exact MCTDH simulations for a fully anharmonic seven-dimensional model of the I2Kr18 cluster in a crystal Kr cage. The linear absorption spectrum, time-evolving vibrational probability densities, and I2 energy expectation value are accurately reproduced by the numerically efficient G-MCTDH approach. The reduced density matrix of the chromophore is analyzed in the coordinate, Wigner and energy representations, so as to obtain a multifaceted dynamical view of the guest-host interactions. Vibrational coherences extending over the bond distance range 2.7 Å < RI–I < 4.0 Å are found to survive for several vibrational periods, despite extensive dissipation. The present results prepare the ground for the simulation of time-resolved coherent Raman spectroscopy of the I2-krypton system addressed in Paper II. [ABSTRACT FROM AUTHOR]

This paper presents a new approach to propagating the density matrix based on a time stepping procedure arising from a Trotter factorization and combining the forward and backward incremental propagators. The sums over intermediate states of the discrete quantum subsystem are implemented by a Monte Carlo surface hopping-like procedure, while the integrals over the continuous variables are performed using a linearization in the difference between the forward and backward paths of these variables leading to classical-like equations of motion with forces determined by the quantum subsystem states. The approach is tested on several models and numerical convergence is explored. [ABSTRACT FROM AUTHOR]

In a recent Letter [Chem. Phys. Lett. 221, 482 (1994)], we demonstrated that the dynamics of reduced density matrices for systems in contact with dissipative harmonic environments can be obtained in an iterative fashion by multiplication of a propagator tensor. The feasibility of iterative procedures in reduced dimension spaces arises from intrinsic features of the dissipative influence functional in Feynman’s path integral formulation of quantum dynamics. Specifically, the continuum of frequencies characteristic of broad condensed phase spectra disrupts phase coherence to a large extent, such that the dynamics of an augmented reduced density tensor becomes Markovian. In a preceding article [J. Chem. Phys. 102, 4600 (1995)] we examined in detail the formal properties of the tensor propagator. In the present paper we show that the tensor propagator can be further decomposed into a product of small rank tensors, resulting in an extremely simple and efficient numerical scheme that scales almost linearly with the dimension of the augmented reduced density tensor. Numerical application to a model electron transfer reaction is presented. © 1995 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The switching of single cyclooctadiene molecules chemisorbed on a Si(100) surface between two stable conformations, can be achieved with a scanning tunneling microscope [Nacci et al., Phys. Rev. B 77, 121405(R) (2008)]. Recently, it was shown by quantum chemical and quantum dynamical simulations that major experimental facts can be explained by a single-mode model with switching enforced by inelastic electron tunneling (IET) excitations and perturbed by vibrational relaxation [Nacci et al., Nano Lett. 9, 2997 (2009)]. In the present paper, we extend the previous theoretical work in several respects: (1) The model is generalized to a two-mode description in which two C2H4 units of COD can move independently; (2) contributions of dipole and, in addition, (cation and anion) resonance-IET rates are considered; (3) the harmonic-linear vibrational relaxation model used previously is generalized to anharmonic vibrations. While the present models highlight generic aspects of IET-switching between two potential minima, they also rationalize specific experimental findings for COD/Si(100): (1) A single-electron excitation mechanism with a linear dependence of the switching rate on tunneling current I, (2) the capability to switch both at negative and positive sample biases, and (3) a crossover temperature around ∼60 K from an IET-driven, T-independent atom tunneling regime, to classical over-the-barrier isomerization with exponential T-dependence at higher temperatures for a bias voltage of +1.5 V and an average tunneling current of 0.73 nA. [ABSTRACT FROM AUTHOR]

The reformulation of the standard golden rule approach considered in this paper for treating reactive tunneling reduces the computation of the reaction rate to a derivation of band shapes for energy levels of reactant and product states. This treatment is based on the assumption that the medium environment is actively involved as a partner in the energy exchange with the reactive subsystem but its reorganization effect is negligible. Starting from the quantum relaxation equation for the density matrix, the required band shapes are represented in terms of the spectral density function, exhibiting the continuum spectrum inherent to the interaction between the reactants and the medium in the total reactive system. The simplest Lorentzian spectral bands, obtained under Redfield approximation, proved to be unsatisfactory because they produced a divergent rate expression at low temperature. The problem is resolved by invoking a refined spectral band shape, which behaves as Lorentzian one at the band center but decays exponentially at its tails. The corresponding closed non-Markovian rate expression is derived and investigated taking as an example the photochemical H-transfer reaction between fluorene and acridine proceeding in the fluorene molecular crystal. The kinetics in this reactive system was thoroughly studied experimentally in a wide temperature range [B. Prass et al., Ber. Bunsenges. Phys. Chem. 102, 498 (1998)]. [ABSTRACT FROM AUTHOR]

Errors in the control of quantum systems may be classified as unitary, decoherent, and incoherent. Unitary errors are systematic, and result in a density matrix that differs from the desired one by a unitary operation. Decoherent errors correspond to general completely positive superoperators, and can only be corrected using methods such as quantum error correction. Incoherent errors can also be described, on average, by completely positive superoperators, but can nevertheless be corrected by the application of a locally unitary operation that “refocuses” them. They are due to reproducible spatial or temporal variations in the system’s Hamiltonian, so that information on the variations is encoded in the system’s spatiotemporal state and can be used to correct them. In this paper liquid-state nuclear magnetic resonance is used to demonstrate that such refocusing effects can be built directly into the control fields, where the incoherence arises from spatial inhomogeneities in the quantizing static magnetic field as well as the radio-frequency control fields themselves. Using perturbation theory, it is further shown that the eigenvalue spectrum of the completely positive superoperator exhibits a characteristic spread that contains information on the Hamiltonians’ underlying distribution. © 2003 American Institute of Physics. [ABSTRACT FROM AUTHOR]

We develop an accurate and numerically efficient non-adiabatic path-integral approach to simulate the non-linear spectroscopy of exciton–polariton systems. This approach is based on the partial linearized density matrix approach to model the exciton dynamics with explicit propagation of the phonon bath environment, combined with a stochastic Lindblad dynamics approach to model the cavity loss dynamics. Through simulating both linear and polariton two-dimensional electronic spectra, we systematically investigate how light–matter coupling strength and cavity loss rate influence the optical response signal. Our results confirm the polaron decoupling effect, which is the reduced exciton–phonon coupling among polariton states due to the strong light–matter interactions. We further demonstrate that the polariton coherence time can be significantly prolonged compared to the electronic coherence outside the cavity. [ABSTRACT FROM AUTHOR]

For common condensed phase problems described by a low-dimensional system coupled to a harmonic bath, Feynman’s path integral formulation of time-dependent quantum mechanics leads to expressions for the reduced density matrix of the system where the effects of the harmonic environment enter through an influence functional that is nonlocal in time. In a recent Letter [Chem. Phys. Lett. 221, 482 (1994)], we demonstrated that the range of the nonlocal interactions is finite even at zero temperature, such that the nonlocal kernel extends over only a few time steps if the path integral is expressed in terms of accurate quasiadiabatic propagators. This feature arises from disruption of phase coherence in macroscopic environments and leads to Markovian dynamics for an augmented reduced density tensor, permitting iterative time evolution schemes. In the present paper we analyze the structure and properties of the relevant tensor propagator. Specifically, we show that the tensor multiplication scheme rigorously conserves the trace of the reduced density matrix, and that in cases of short-range nonlocality it leads to Redfield-type equations which are correct to all orders in perturbation theory and which take into account memory effects. We also argue that a simple eigenvector analysis reveals (without actual iteration) the nature of the dynamics and of the equilibrium state, and directly yields quantum reaction or relaxation rates. © 1995 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Kohn–Sham density functional theory calculations using conventional diagonalization based methods become increasingly expensive as temperature increases due to the need to compute increasing numbers of partially occupied states. We present a density matrix based method for Kohn–Sham calculations at high temperatures that eliminates the need for diagonalization entirely, thus reducing the cost of such calculations significantly. Specifically, we develop real-space expressions for the electron density, electronic free energy, Hellmann–Feynman forces, and Hellmann–Feynman stress tensor in terms of an orthonormal auxiliary orbital basis and its density kernel transform, the density kernel being the matrix representation of the density operator in the auxiliary basis. Using Chebyshev filtering to generate the auxiliary basis, we next develop an approach akin to Clenshaw–Curtis spectral quadrature to calculate the individual columns of the density kernel based on the Fermi operator expansion in Chebyshev polynomials and employ a similar approach to evaluate band structure and entropic energy components. We implement the proposed formulation in the SPARC electronic structure code, using which we show systematic convergence of the aforementioned quantities to exact diagonalization results, and obtain significant speedups relative to conventional diagonalization based methods. Finally, we employ the new method to compute the self-diffusion coefficient and viscosity of aluminum at 116 045 K from Kohn–Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods. [ABSTRACT FROM AUTHOR]

Tensor network decompositions of path integrals for simulating open quantum systems have recently been proven to be useful. However, these methods scale exponentially with the system size. This makes it challenging to simulate the non-equilibrium dynamics of extended quantum systems coupled with local dissipative environments. In this work, we extend the tensor network path integral (TNPI) framework to efficiently simulate such extended systems. The Feynman–Vernon influence functional is a popular approach used to account for the effect of environments on the dynamics of the system. In order to facilitate the incorporation of the influence functional into a multisite framework (MS-TNPI), we combine a matrix product state (MPS) decomposition of the reduced density tensor of the system along the sites with a corresponding tensor network representation of the time axis to construct an efficient 2D tensor network. The 2D MS-TNPI network, when contracted, yields the time-dependent reduced density tensor of the extended system as an MPS. The algorithm presented is independent of the system Hamiltonian. We outline an iteration scheme to take the simulation beyond the non-Markovian memory introduced by solvents. Applications to spin chains coupled to local harmonic baths are presented; we consider the Ising, XXZ, and Heisenberg models, demonstrating that the presence of local environments can often dissipate the entanglement between the sites. We discuss three factors causing the system to transition from a coherent oscillatory dynamics to a fully incoherent dynamics. The MS-TNPI method is useful for studying a variety of extended quantum systems coupled with solvents. [ABSTRACT FROM AUTHOR]

Pulay’s original expression for the basis-set dependent adjustment term to the Hellmann–Feynman force in electronic structure theory, which occurs for nonorthogonal local basis-set representations, is based on the idempotency condition of a pure ensemble. At finite electronic temperatures with a fractional occupation of the states, the conventional expression of the Pulay force is therefore no longer valid. Here we derive a simple and computationally efficient expression for a generalized Pulay force, which is suitable for large-scale ab initio simulations at finite electronic temperatures using local nonorthogonal basis-set representations. The generalized Pulay force expression is given in terms of the temperature-dependent density matrix. For the construction of the density matrix, we propose a recursive Fermi operator expansion algorithm that automatically converges to the correct chemical potential. [ABSTRACT FROM AUTHOR]

We introduce a practical method for compacting the time evolution of the quantum state of a closed physical system. The density matrix is specified as a function of a few time-independent observables where their coefficients are time-dependent. The key mathematical step is the vectorization of the surprisal, the logarithm of the density matrix, at each time point of interest. The time span used depends on the required spectral resolution. The entire course of the system evolution is represented as a matrix where each column is the vectorized surprisal at the given time point. Using the singular value decomposition (SVD) of this matrix, we generate realistic approximations for the time-independent observables and their respective time-dependent coefficients. This allows for a simplification of the algebraic procedure for determining the dominant constraints (the time-independent observables) in the sense of the maximal entropy approach. A non-stationary coherent initial state of a Morse oscillator is used to introduce the approach. We derive the analytical exact expression for the surprisal as a function of time, and this offers a benchmark for comparison with the accurate but approximate SVD results. We discuss two examples of a Morse potential of different anharmonicities, H2 and I2 molecules. We further demonstrate the approach for a two-coupled electronic state problem, the well-studied non-radiative decay of pyrazine from its bright state. Five constraints are found to be enough to capture the ultrafast electronic population exchange and to recover the dynamics of the wave packet in both electronic states. [ABSTRACT FROM AUTHOR]

This work explores the possibility of simulating an electron transfer process between a donor and an acceptor in real time using time-dependent density functional theory electron dynamics. To achieve this objective, a central issue to resolve is the definition of the initial state. This must be a non-equilibrium electronic state able to trigger the charge transfer dynamics; here, two schemes are proposed to prepare such states. One is based on the combination of the density matrices of the donor and acceptor converged separately with appropriate charges (for example, −1 for the donor and +1 for the acceptor). The second approach relied on constrained DFT to localize the charge on each fragment. With these schemes, electron transfer processes are simulated in different model systems of increasing complexity: an atomic hydrogen dimer, a polyacetylene chain, and the active site of the T. cruzi hybrid type A heme peroxidase, for which two possible electron transfer paths have been postulated. For the latter system, the present methodology applied in a hybrid Quantum Mechanics - Molecular Mechanics framework allows us to establish the relative probabilities of each path and provides insight into the inhibition of the electron transfer provoked by the substitution of tryptophan by phenylalanine in the W233F mutant. [ABSTRACT FROM AUTHOR]

Marcus theory has been successfully applied to molecular design for organic semiconductors with the aid of quantum chemistry calculations for the molecular parameters: the intermolecular electronic coupling V and the intramolecular charge reorganization energy λ. The assumption behind this is the localized nature of the electronic state for representing the charge carriers, being holes or electrons. As far as the quantitative description of carrier mobility is concerned, the direct application of Marcus semiclassical theory usually led to underestimation of the experimental data. A number of effects going beyond such a semiclassical description will be introduced here, including the quantum nuclear effect, dynamic disorder, and delocalization effects. The recently developed quantum dynamics simulation at the time-dependent density matrix renormalization group theory is briefly discussed. The latter was shown to be a quickly emerging efficient quantum dynamics method for the complex system. [ABSTRACT FROM AUTHOR]

The multiconfiguration time-dependent Hartree (MCTDH) method is a powerful method for solving the time-dependent Schrödinger equation in quantum molecular dynamics. It is, however, hampered by the so-called curse of dimensionality which results in exponential scaling with respect to the number of degrees of freedom in the system and, thus, limits its applicability to small- and medium-sized molecules. To avoid this scaling, we derive equations of motion for a series of truncated MCTDH methods using a many-mode second-quantization formulation where the configuration space is restricted based on mode-combination levels as also done in the vibrational configuration interaction and vibrational coupled cluster methods for solving the time-independent Schrödinger equation. The full MCTDH wave function is invariant with respect to the choice of constraint (or gauge) operators, but restricting the configuration space removes this invariance. We, thus, analyze the remaining redundancies and derive equations for variationally optimizing the non-redundant matrix elements of the constraint operators. As an alternative, we also present a constraint that keeps the density matrices block diagonal during the propagation and the two choices are compared. Example calculations are performed on formyl fluoride and a series of high-dimensional Henon–Heiles potentials. The results show that the MCTDH[n] methods can be applied to large systems and that an optimal choice of constraint operators is key to obtaining the correct physical behavior of the wave function. [ABSTRACT FROM AUTHOR]

We extend the coupled-cluster method to correlated quantum dynamics of both closed and open systems at finite temperatures using the thermofield formalism. The approach expresses the time-dependent density matrix in an exponential ansatz and describes time-evolution along the Keldysh path contour. A distinct advantage of the approach is exact trace-preservation as a function of time, ensuring conservation of probability and particle number. Furthermore, the method avoids the computation of correlated bra-states, simplifying the computational implementation. We develop the method in a thermal quasiparticle representation, which allows seamless connection to the projection method and diagrammatic techniques of the traditional coupled-cluster formalism. For comparison, we also apply the thermofield framework to the density-matrix renormalization-group method to obtain reference results for closed and open systems at finite temperature. We test the singles and doubles approximation to the density-matrix coupled-cluster method on the correlated electronic dynamics of the single-impurity Anderson model, demonstrating that the new method successfully captures the correlated dynamics of both closed systems at finite temperature and driven-dissipative open systems. This encouraging performance motivates future applications to nonequilibrium quantum many-body dynamics in realistic systems. [ABSTRACT FROM AUTHOR]

In this paper, we investigate coherent quantum dynamics in a nonequilibrium environment. We focus on a two-state quantum system strongly coupled to a single classical environmental oscillator, and explore the effect of nonstationary statistical properties of the oscillator on the quantum evolution. A simple nonequilibrium model, consisting of an oscillator with a well-defined initial phase which undergoes subsequent diffusion, is introduced and studied. Approximate but accurate analytic expressions for the evolution of the off-diagonal density matrix element of the quantum system are derived in the second-order cumulant approximation. The effect of the initial phase choice on the subsequent quantum evolution is quantified. It is observed that the initial phase can have a significant effect on the preservation of coherence on short time scales, suggesting this variable as a control parameter for optimizing coherence in many-body quantum systems. [ABSTRACT FROM AUTHOR]

Density matrix embedding theory (DMET) is a quantum embedding theory for strongly correlated systems. From a computational perspective, one bottleneck in DMET is the optimization of the correlation potential to achieve self-consistency, especially for heterogeneous systems of large size. We propose a new method, called projected DMET (p-DMET), which achieves self-consistency without needing to optimize the correlation potential. We demonstrate the performance of p-DMET on the two-dimensional Hubbard model. [ABSTRACT FROM AUTHOR]

Copyright of Journal of Chemical Physics is the property of American Institute of Physics and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

The iterative decomposition of the blip-summed path integral [N. Makri, J. Chem. Phys. 141, 134117 (2014)] is described. The starting point is the expression of the reduced density matrix for a quantum system interacting with a harmonic dissipative bath in the form of a forward-backward path sum, where the effects of the bath enter through the Feynman-Vernon influence functional. The path sum is evaluated iteratively in time by propagating an array that stores blip configurations within the memory interval. Convergence with respect to the number of blips and the memory length yields numerically exact results which are free of statistical error. In situations of strongly dissipative, sluggish baths, the algorithm leads to a dramatic reduction of computational effort in comparison with iterative path integral methods that do not implement the blip decomposition. This gain in efficiency arises from (i) the rapid convergence of the blip series and (ii) circumventing the explicit enumeration of between-blip path segments, whose number grows exponentially with the memory length. Application to an asymmetric dissipative two-level system illustrates the rapid convergence of the algorithm even when the bath memory is extremely long. [ABSTRACT FROM AUTHOR]

We present a formalism that explicitly unifies the commonly used Nakajima-Zwanzig approach for reduced density matrix dynamics with the more versatile Mori theory in the context of nonequilibrium dynamics. Employing a Dyson-type expansion to circumvent the difficulty of projected dynamics, we obtain a self-consistent equation for the memory kernel which requires only knowledge of normally evolved auxiliary kernels. To illustrate the properties of the current approach, we focus on the spin-boson model and limit our attention to the use of a simple and inexpensive quasi-classical dynamics, given by the Ehrenfest method, for the calculation of the auxiliary kernels. For the first time, we provide a detailed analysis of the dependence of the properties of the memory kernels obtained via different projection operators, namely, the thermal (Redfield-type) and population based (NIBA-type) projection operators. We further elucidate the conditions that lead to short-lived memory kernels and the regions of parameter space to which this program is best suited. Via a thorough analysis of the different closures available for the auxiliary kernels and the convergence properties of the self-consistently extracted memory kernel, we identify the mechanisms whereby the current approach leads to a significant improvement over the direct usage of standard semi- and quasi-classical dynamics. [ABSTRACT FROM AUTHOR]

We propose a way of obtaining effective low energy Hubbard-like model Hamiltonians from ab initio quantum Monte Carlo calculations for molecular and extended systems. The Hamiltonian parameters are fit to best match the ab initio two-body density matrices and energies of the ground and excited states, and thus we refer to the method as ab initio density matrix based downfolding. For benzene (a finite system), we find good agreement with experimentally available energy gaps without using any experimental inputs. For graphene, a two dimensional solid (extended system) with periodic boundary conditions, we find the effective on-site HubbardU./t to be 1.3 ¡¾ 0.2, comparable to a recent estimate based on the constrained random phase approximation. For molecules, such parameterizations enable calculation of excited states that are usually not accessible within ground state approaches. For solids, the effective Hamiltonian enables large-scale calculations using techniques designed for lattice models. [ABSTRACT FROM AUTHOR]

The exact conditions for density functionals and density matrix functionals in terms of fractional charges and fractional spins are known, and their violation in commonly used functionals has been shown to be the root of many major failures in practical applications. However, approximate functionals are designed for physical systems with integer charges and spins, not in terms of the fractional variables. Here we develop a general framework for extending approximate density functionals and many-electron theory to fractional-charge and fractional-spin systems. Our development allows for the fractional extension of any approximate theory that is a functional of G0, the one-electron Green's function of the non-interacting reference system. The extension to fractional charge and fractional spin systems is based on the ensemble average of the basic variable, G0. We demonstrate the fractional extension for the following theories: (1) any explicit functional of the one-electron density, such as the local density approximation and generalized gradient approximations; (2) any explicit functional of the one-electron density matrix of the non-interacting reference system, such as the exact exchange functional (or Hartree-Fock theory) and hybrid functionals; (3) many-body perturbation theory; and (4) random-phase approximations. A general rule for such an extension has also been derived through scaling the orbitals and should be useful for functionals where the link to the Green's function is not obvious. The development thus enables the examination of approximate theories against known exact conditions on the fractional variables and the analysis of their failures in chemical and physical applications in terms of violations of exact conditions of the energy functionals. The present work should facilitate the calculation of chemical potentials and fundamental bandgaps with approximate functionals and many-electron theories through the energy derivatives with respect to the fractional charge. It should play an important role in developing accurate approximate density functionals and many-body theory. [ABSTRACT FROM AUTHOR]

Quantum-mechanical evolution of systems with periodic time-modulated Hamiltonians is often described by effective interactions. Such average Hamiltonians, calculated as few terms of an expansion in powers of the interaction, are sometimes difficult to relate to experimental observations. We propose a frequency-domain approach to this problem, which offers certain advantages and produces an approximate solution for the density matrix, better linked to measurable quantities. The formalism is suitable for calculating the intensities of narrowed spectral peaks. Fast magic-angle-spinning NMR spectra of solids are used to experimentally illustrate the method. [ABSTRACT FROM AUTHOR]

In the study of open quantum systems, the polaron transformation has recently attracted a renewed interest as it offers the possibility to explore the strong system-bath coupling regime. Despite this interest, a clear and unambiguous analysis of the regimes of validity of the polaron transformation is still lacking. Here we provide such a benchmark, comparing second order perturbation theory results in the original untransformed frame, the polaron frame, and the variational extension with numerically exact path integral calculations of the equilibrium reduced density matrix. Equilibrium quantities allow a direct comparison of the three methods without invoking any further approximations as is usually required in deriving master equations. It is found that the second order results in the original frame are accurate for weak system-bath coupling; the results deteriorate when the bath cut-off frequency decreases. The full polaron results are accurate for the entire range of coupling for a fast bath but only in the strong coupling regime for a slow bath. The variational method is capable of interpolating between these two methods and is valid over a much broader range of parameters. [ABSTRACT FROM AUTHOR]

A novel scheme for the steady state solution of the standard Redfield quantum master equation is developed which yields agreement with the exact result for the corresponding reduced density matrix up to second order in the system-bath coupling strength. We achieve this objective by use of an analytic continuation of the off-diagonal matrix elements of the Redfield solution towards its diagonal limit. Notably, our scheme does not require the provision of yet higher order relaxation tensors. Testing this modified method for a heat bath consisting of a collection of harmonic oscillators we assess that the system relaxes towards its correct coupling-dependent, generalized quantum Gibbs state in second order. We numerically compare our formulation for a damped quantum harmonic system with the nonequilibrium Green's function formalism: we find good agreement at low temperatures for coupling strengths that are even larger than expected from the very regime of validity of the second-order Redfield quantum master equation. Yet another advantage of our method is that it markedly reduces the numerical complexity of the problem; thus, allowing to study efficiently large-sized system Hilbert spaces. [ABSTRACT FROM AUTHOR]

Electronic energy transfer in the condensed phase, such as that occurring in photosynthetic complexes, frequently occurs in regimes where the energy scales of the system and environment are similar. This situation provides a challenge to theoretical investigation since most approaches are accurate only when a certain energetic parameter is small compared to others in the problem. Here we show that in these difficult regimes, the Ehrenfest approach provides a good starting point for a dynamical description of the energy transfer process due to its ability to accurately treat coupling to slow environmental modes. To further improve on the accuracy of the Ehrenfest approach, we use our reduced density matrix hybrid framework to treat the faster environmental modes quantum mechanically, at the level of a perturbative master equation. This combined approach is shown to provide an efficient and quantitative description of electronic energy transfer in a model dimer and the Fenna-Matthews-Olson complex and is used to investigate the effect of environmental preparation on the resulting dynamics. [ABSTRACT FROM AUTHOR]

Several methods for density matrix propagation in parallel computing environments are proposed and evaluated. It is demonstrated that the large communication overhead associated with each propagation step (two-sided multiplication of the density matrix by an exponential propagator and its conjugate) may be avoided and the simulation recast in a form that requires virtually no inter-thread communication. Good scaling is demonstrated on a 128-core (16 nodes, 8 cores each) cluster. [ABSTRACT FROM AUTHOR]

We present a new approach to calculate real-time quantum dynamics in complex systems. The formalism is based on the partitioning of a system's environment into 'core' and 'reservoir' modes with the former to be treated quantum mechanically and the latter classically. The presented method only requires the calculation of the system's reduced density matrix averaged over the quantum core degrees of freedom which is then coupled to a classically evolved reservoir to treat the remaining modes. We demonstrate our approach by applying it to the spin-boson problem using the noninteracting blip approximation to treat the system and core, and Ehrenfest dynamics to treat the reservoir. The resulting hybrid methodology is accurate for both fast and slow baths, since it naturally reduces to its composite methods in their respective regimes of validity. In addition, our combined method is shown to yield good results in intermediate regimes where neither approximation alone is accurate and to perform equally well for both strong and weak system-bath coupling. Our approach therefore provides an accurate and efficient methodology for calculating quantum dynamics in complex systems. [ABSTRACT FROM AUTHOR]

Starting from the equation of motion in the density matrix formulation, we reformulate the analytical gradient of the excited-state energy at the time-dependent density functional theory level in the nonorthogonal Gaussian atom-centered orbital (AO) basis. Analogous to the analytical first derivative in molecular-orbital (MO) basis, a Z-vector equation has been derived with respect to the reduced one-electronic density matrix in AO basis, which provides a potential possibility to exploit quantum locality of the density matrix and avoids the matrix transformation between the AO and the MO basis. Numerical tests are finished for the excited-state geometry optimization and adiabatic excitation energy calculation of a series of small molecules. The results demonstrate the computational efficiency and accuracy of the current AO-based energy gradient expression in comparison with the MO-based scheme. [ABSTRACT FROM AUTHOR]

Despite the fact that idempotent one-particle reduced density matrices are pervasive in quantum chemistry, the understanding of a general energy functional of such idempotent density matrices for the ground state energy has been lacking. By a constrained search, we show the structure of the general functional, illuminating the contributions from various terms. For the examples of the 'best idempotent density matrix' and Kohn-Sham idempotent density matrices, we contrast the functional forms and suggest how the best idempotent density matrix approach may be a good starting point for further development. [ABSTRACT FROM AUTHOR]

A double-index atomic partitioning of the molecular first-order density matrix is proposed. Contributions diagonal in the atomic indices correspond to atomic density matrices, whereas off-diagonal contributions carry information about the bonds. The resulting matrices have good localization properties, in contrast to single-index atomic partitioning schemes of the molecular density matrix. It is shown that the electron density assigned to individual atoms, when derived from the density matrix partitioning, can be made consistent with well-known partitions of the electron density over atom in the molecule basins, either with sharp or with fuzzy boundaries. The method is applied to a test set of about 50 molecules, representative for various types of chemical binding. A close correlation is observed between the trace of the bond matrices and the shared electron density index. [ABSTRACT FROM AUTHOR]

We propose a spin-free approach to the cumulant decomposition of reduced density matrices of singlet and spin-rotation or SU(2) invariant ensemble of nonsinglet states as in [W. Kutzelnigg and D. Mukherjee, J. Chem. Phys. 110, 2800 (1999); 116, 4787 (2002)]. We provide a simple recursive procedure to obtain expressions which relate different spin components of spin-orbital reduced density matrices and cumulants of such states to the spin-free counterparts. These results are used to set up a spin-summation procedure to arrive at a definition of spin-free cumulants of any order. Alternatively, an analytic formula for the spin-free form resulting from a spin summation involving product of two spin-orbital cumulants is derived and its utility in spin-free cumulant decomposition of reduced density matrices is demonstrated. This leads to suitable definitions of spin-free analog of multireference normal ordering and the associated Wicks theorem. The results of this formulation are expected to be useful in investigations of spin-free multireference internally contracted coupled-cluster methods where cumulant approximations to the active reduced density matrices are employed. [ABSTRACT FROM AUTHOR]

We develop monotonically convergent free-time and fixed end-point optimal control theory (OCT) in the density-matrix representation to deal with quantum systems showing dissipation. Our theory is more general and flexible for tailoring optimal laser pulses in order to control quantum dynamics with dissipation than the conventional fixed-time and fixed end-point OCT in that the optimal temporal duration of laser pulses can also be optimized exactly. To show the usefulness of our theory, it is applied to the generation and maintenance of the vibrational entanglement of carbon monoxide adsorbed on the copper (100) surface, CO/Cu(100). We demonstrate the numerical results and clarify how to combat vibrational decoherence as much as possible by the tailored shapes of the optimal laser pulses. It is expected that our theory will be general enough to be applied to a variety of dissipative quantum dynamics systems because the decoherence is one of the quantum phenomena sensitive to the temporal duration of the quantum dynamics. [ABSTRACT FROM AUTHOR]

This article presents an efficient and parallelized implementation of the density matrix renormalization group (DMRG) algorithm for quantum chemistry calculations. The DMRG method as a large-scale multireference electronic structure model is by nature particularly efficient for one-dimensionally correlated systems, while the present development is oriented toward applications for polynuclear transition metal compounds, in which the macroscopic one-dimensional structure of electron correlation is absent. A straightforward extension of the DMRG algorithm is proposed with further improvements and aggressive optimizations to allow its application with large multireference active space, which is often demanded for metal compound calculations. Special efficiency is achieved by making better use of sparsity and symmetry in the operator and wave function representations. By accomplishing computationally intensive DMRG calculations, the authors have found that a large number of renormalized basis states are required to represent high entanglement of the electron correlation for metal compound applications, and it is crucial to adopt auxiliary perturbative correction to the projected density matrix during the DMRG sweep optimization in order to attain proper convergence to the solution. Potential energy curve calculations for the Cr2 molecule near the known equilibrium precisely predicted the full configuration interaction energies with a correlation space of 24 electrons in 30 orbitals [denoted by (24e,30o)]. The energies are demonstrated to be accurate to 0.6mEh (the error from the extrapolated best value) when as many as 10 000 renormalized basis states are employed for the left and right DMRG block representations. The relative energy curves for [Cu2O2]2+ along the isomerization coordinate were obtained from DMRG and other correlated calculations, for which a fairly large orbital space (32e,62o) is modeled as a full correlation space. The DMRG prediction nearly overlaps with the energy curve from the coupled cluster with singles, doubles, and perturbative triple [CCSD(T)] calculations, while the multireference complete active space self-consistent field (CASSCF) calculations with the small reference configuration (8e,8o) are found to overestimate the biradical character of the electronic state of [Cu2O2]2+ according to the one-electron density matrix analysis. [ABSTRACT FROM AUTHOR]

We have investigated the radial electron pair probability distributions (REPPDs) of the helium dimer within the Piris natural orbital functional (PNOF) theory. The analytical formulas to evaluate intracule densities, Fermi, Coulomb, and total correlation holes using our reconstruction functional PNOF-2 [J. Chem. Phys. 126, 214103 (2007)] are derived. The Löwdin’s Coulomb holes from PNOF-2 and full configuration interaction calculations are analyzed showing a very similar behavior. New definitions of the Coulomb and Fermi holes based on the cumulant expansion of the two-particle reduced density matrix are presented. The holes are defined in terms of the exact one-particle reduced density matrix and the two-particle cumulant without any reference to the Hartree–Fock state. Through these definitions, we analyze separately the contribution of each component to the total REPPD at several values of the internuclear distance. A straight connection between the Coulomb hole and dispersion interactions is observed. [ABSTRACT FROM AUTHOR]

Linear scaling density matrix perturbation theory [A. M. N. Niklasson and M. Challacombe, Phys. Rev. Lett. 92, 193001 (2004)] is extended to basis-set-dependent quantum response calculations for a nonorthogonal basis set representation. The generalization is achieved by a perturbation-dependent congruence transform, derived from the factorization of the inverse overlap matrix, which transforms the generalized eigenvalue problem to an orthogonal, standard form. With this orthogonalization transform the basis-set-dependent perturbation in the overlap matrix is included in the orthogonalized Hamiltonian, which is expanded in orders of the perturbation. In this way density matrix perturbation theory developed for an orthogonal representation can be applied also to basis-set-dependent response calculations. The method offers an alternative to the previous solution of the basis-set-dependent response problem, based on a nonorthogonal generalization of the density matrix perturbation theory, where the calculations are performed within a purely nonorthogonal setting [A. M. N. Niklasson et al., J. Chem. Phys. 123, 44107 (2005)]. [ABSTRACT FROM AUTHOR]

The authors have investigated the description of the dispersion interaction within the Piris natural orbital functional (PNOF) theory. The PNOF arises from an explicit antisymmetric approach for the two-particle cumulant in terms of two symmetric matrices, Δ and Λ. The functional forms of these matrices are obtained from the generalization of the two-particle system expressions, except for the off-diagonal elements of Δ. The mean value theorem and the partial sum rule obtained for the off-diagonal elements of Δ provide a prescription for deriving practical functionals. In particular, the previous employed approximation {Jpp/2} for the mean values {Jp*} affords several molecular properties but it is incapable to account for dispersion effects. In this work, the authors analyze a new approach for Jp* obtained by factorization of the matrix Δ within the bounds on its off-diagonal elements imposed by the positivity conditions of the two-particle reduced density matrix. Additional terms for the matrix elements of Λ proportional to the square root of the holes are again introduced to describe properly the occupation numbers of the lowest occupied levels. The authors have found that the cross products between weakly occupied orbitals must be removed from the functional form of Λ to obtain a correct long-range asymptotic behavior. The PNOF is used to predict the binding energy as well as the equilibrium distance of the helium dimer. The results are compared with the full configuration-interaction calculations and the corresponding experimental data. [ABSTRACT FROM AUTHOR]

Two-electron reduced density matrices (2-RDMs) have recently been directly determined from the solution of the anti-Hermitian contracted Schrödinger equation (ACSE) to obtain 95%–100% of the ground-state correlation energy of atoms and molecules, which significantly improves upon the accuracy of the contracted Schrödinger equation (CSE) [D. A. Mazziotti, Phys. Rev. Lett. 97, 143002 (2006)]. Two subsets of the CSE, the ACSE and the contraction of the CSE onto the one-particle space, known as the 1,3-CSE, have two important properties: (i) dependence upon only the 3-RDM and (ii) inclusion of all second-order terms when the 3-RDM is reconstructed as only a first-order functional of the 2-RDM. The error in the 1,3-CSE has an important role as a stopping criterion in solving the ACSE for the 2-RDM. Using a computationally more efficient implementation of the ACSE, the author treats a variety of molecules, including H2O, NH3, HCN, and HO3-, in larger basis sets such as correlation-consistent polarized double- and triple-zeta. The ground-state energy of neon is also calculated in a polarized quadruple-zeta basis set with extrapolation to the complete basis-set limit, and the equilibrium bond length and harmonic frequency of N2 are computed with comparison to experimental values. The author observes that increasing the basis set enhances the ability of the ACSE to capture correlation effects in ground-state energies and properties. In the triple-zeta basis set, for example, the ACSE yields energies and properties that are closer in accuracy to coupled cluster with single, double, and triple excitations than to coupled cluster with single and double excitations. In all basis sets, the computed 2-RDMs very closely satisfy known N-representability conditions. [ABSTRACT FROM AUTHOR]

Several techniques to solve a hierarchical set of equations of motion for propagating a reduced density matrix coupled to a thermal bath have been developed in recent years. This is either done using the path integral technique as in the original proposal by Tanimura and Kubo [J. Phys. Soc. Jpn. 58, 101 (1998)] or by the use of stochastic fields as done by Yan et al. [Chem. Phys. Lett. 395, 216 (2004)]. Based on the latter ansatz a compact derivation of the hierarchy using a decomposition of the spectral density function is given in the present contribution. The method is applied to calculate the time evolution of the reduced density matrix describing the motion in a harmonic, an anharmonic, and two coupled oscillators where each system is coupled to a thermal bath. Calculations to several orders in the system-bath coupling with two different truncations of the hierarchy are performed. The respective density matrices are used to calculate the time evolution of various system properties and the results are compared and discussed with a special focus on the convergence with respect to the truncation scheme applied. [ABSTRACT FROM AUTHOR]

One-electron density matrices resulting from the explicitly connected commutator expansion of the expectation value were implemented at the singles and doubles coupled cluster (CCSD) level. In the proposed approach the one-electron density matrix is obtained at a little extra cost in comparison to the calculation of the CCSD correlation energy. Therefore, in terms of the computational time the new method is significantly less demanding than the conventional linear-response CCSD theory which requires additionally an expensive calculation of the left-hand solution of the CCSD equations. The quality of the new density matrices was investigated by computing a set of one-electron properties for a series of molecules of varying sizes and comparing the results with data obtained using the full configuration interaction method or higher level coupled cluster theory. It has been found that the results obtained using the new approach are of the same quality as those predicted by the linear-response CCSD method. The novel one-electron density matrices have also been applied to study the energy of the electrostatic interaction for a number of van der Waals complexes, including the benzene and azulene dimers. [ABSTRACT FROM AUTHOR]

Starting point is the energy expectation value as a functional of the one-particle density matrix γ and the two-particle density cumulant λ2. We decompose γ into a best idempotent approximation κ and a correction τ, that is entirely expressible in terms of λ2. So we get the energy E as a functional of κ and λ2, which can be varied independently. Approximate n-representability conditions, derived by perturbation theory are imposed on the variation of λ2. A nonlinear system of equations satisfied by λ2 is derived, the linearized version of which turns out to be equivalent to the coupled electron-pair approximation, variant zero. The start for κ is Hartree-Fock, but κ is then updated to become the best idempotent approximation of γ. Relations to density matrix functional theory and Kohn-Sham type density functional theory are discussed. [ABSTRACT FROM AUTHOR]

The application of the quantum-chemical density-matrix renormalization group (DMRG) algorithm is cumbersome for complex electronic structures with many active orbitals. The high computational cost is mainly due to the poor convergence of standard DMRG calculations. A factor which affects the convergence behavior of the calculations is the choice of the start-up procedure. In this start-up step matrix representations of operators have to be calculated in a guessed many-electron basis of the DMRG environment block. Different possibilities for the construction of these basis states exist, and we first compare four procedures to approximate the environment states using Slater determinants explicitly. These start-up procedures are applied to DMRG calculations on a sophisticated test system: the chromium dimer. It is found that the converged energies and the rate of convergence depend significantly on the choice of the start-up procedure. However, since already the most simple start-up procedure, which uses only the Hartree-Fock determinant, is comparatively good, Slater determinants, in general, appear not to be a good choice as approximate environment basis states for convergence acceleration. Based on extensive test calculations it is demonstrated that the computational cost can be significantly reduced if the number of total states m is successively increased. This is done in such a way that the environment states are built up stepwise from system states of previous truncated DMRG sweeps for slowly increasing m values. [ABSTRACT FROM AUTHOR]

The ab initio atom-centered density-matrix propagation approach and the multistate empirical valence bond method have been employed to study the structure, dynamics, and rovibrational spectrum of a hydrated proton in the “magic” 21 water cluster. In addition to the conclusion that the hydrated proton tends to reside on the surface of the cluster, with the lone pair on the protonated oxygen pointing “outwards,” it is also found that dynamical effects play an important role in determining the vibrational properties of such clusters. This result is used to analyze and complement recent experimental and theoretical studies. [ABSTRACT FROM AUTHOR]

The exchange-correlation energy is studied using the density scaling proposed by Chan and Handy [G. K.-L. Chan and N. C. Handy, Phys. Rev. A 59, 2670 (1999)]. It is shown that there exists a value of the scaling factor for which the correlation energy disappears. The optimized potential method and the Krieger-Li-Iafrate approach are generalized to incorporate correlation. [ABSTRACT FROM AUTHOR]