Molecular response properties from explicitly time-dependent configuration interaction methods

In this paper we report the calculation of molecular electric response properties with the help of explicitly time-dependent configuration interaction (TD-CI) methods. These methods have the advantage of being applicable (within the limitations of the time-dependent Schrodinger equation) to time-dependent perturbations of arbitrary shape and strength. Three variants are used to solve the time-dependent electronic Schrodinger equation, namely, the TD-CIS (inclusion of single excitations only), TD-CISD (inclusion of single and double excitations), and TD-CIS(D) (single excitations and perturbative treatment of double excitations) methods and applied for illustration to small molecules, H(2) and H(2)O. In the calculation, slowly varying off-resonant electric fields are applied to the molecules and linear (polarizabilities) and nonlinear (hyperpolarizabilities, harmonic generation) response properties are determined from the time-dependent dipole moments.     

Vibronic energies and spectra of molecular dimers

We consider three distinct methods of calculating the vibronic levels and absorption spectra of molecular dimers coupled by dipole-dipole interactions. The first method is direct diagonalization of the vibronic Hamiltonian in a basis of monomer eigenstates. The second method is to use creation and annihilation operators leading in harmonic approximation to the Jaynes-Cummings Hamiltonian. The adiabatic approximation to this problem provides insight into spectral behavior in the weak and strong coupling limits. The third method, which serves as a check on the accuracy of the previous methods, is a numerically exact solution of the time-dependent Schrodinger equation. Using these methods, dimer spectra are calculated for three separate dye molecules and show good agreement with measured spectra.     

Time-dependent approach to electronically excited states of molecules with the multiconfiguration time-dependent Hartree-Fock method

In this paper the authors show how the multiconfiguration time-dependent Hartree-Fock (MCTDHF) method can be used for the calculation of electronic properties of molecules associated with the population of excited states. In contrast to other methods for correlated electron dynamics, such as configuration interaction, MCTDHF does not rely on a solution of the electronic Schrodinger equation prior to the propagation. The authors apply this approach to the calculation of vertical excitation energies, transition dipole moments, and oscillator strengths for two test molecules, lithium hydride and methane.     

Multiscale quantum propagation using compact-support wavelets in space and time

Orthogonal compact-support Daubechies wavelets are employed as bases for both space and time variables in the solution of the time-dependent Schrodinger equation. Initial value conditions are enforced using special early-time wavelets analogous to edge wavelets used in boundary-value problems. It is shown that the quantum equations may be solved directly and accurately in the discrete wavelet representation, an important finding for the eventual goal of highly adaptive multiresolution Schrodinger equation solvers. While the temporal part of the basis is not sharp in either time or frequency, the Chebyshev method used for pure time-domain propagations is adapted to use in the mixed domain and is able to take advantage of Hamiltonian matrix sparseness. The orthogonal separation into different time scales is determined theoretically to persist throughout the evolution and is demonstrated numerically in a partially adaptive treatment of scattering from an asymmetric Eckart barrier.     

Multidimensional reactive scattering with quantum trajectories: dynamics with Morse vibrational modes

The reactive scattering of a wave packet is studied by the quantum trajectory method for a model system with up to 25 Morse vibrational modes. The equations of motion are formulated in curvilinear reaction path coordinates with the restriction to a planar reaction path. Spatial derivatives are evaluated by the least squares method using contracted basis sets. Dynamical results, including trajectory evolution and time-dependent reaction probabilities, are presented and analyzed. For the case of one Morse vibrational mode, the results are in good agreement with those derived through direct numerical integration of the time-dependent Schrodinger equation.     

Spectral difference Lanczos method for efficient time propagation in quantum control theory

Spectral difference methods represent the real-space Hamiltonian of a quantum system as a banded matrix which possesses the accuracy of the discrete variable representation (DVR) and the efficiency of finite differences. When applied to time-dependent quantum mechanics, spectral differences enhance the efficiency of propagation methods for evolving the Schrodinger equation. We develop a spectral difference Lanczos method which is computationally more economical than the sinc-DVR Lanczos method, the split-operator technique, and even the fast-Fourier-Transform Lanczos method. Application of fast propagation is made to quantum control theory where chirped laser pulses are designed to dissociate both diatomic and polyatomic molecules. The specificity of the chirped laser fields is also tested as a possible method for molecular identification and discrimination.     

Double wells

The new concept of the stability of the quantum mechanical state is introduced for the study of the time-dependent process of the quasistationary state. A stable solution of the time-dependent Schrodinger equation predicts the presence of the new sharply localized long-living state at the most repulsive part of the effective potential.     

Nonvariational time-dependent multiconfiguration self-consistent field equations for electronic dynamics in laser-driven molecules

A time-dependent multiconfiguration self-consistent field (TDMCSCF) scheme is developed to describe the time-resolved electron dynamics of a laser-driven many-electron atomic or molecular system, starting directly from the time-dependent Schrodinger equation for the system. This nonvariational formulation aims at the full exploitations of concepts, tools, and facilities of existing, well-developed quantum chemical MCSCF codes. The theory uses, in particular, a unitary representation of time-dependent configuration mixings and orbital transformations. Within a short-time, or adiabatic approximation, the TDMCSCF scheme amounts to a second-order split-operator algorithm involving generically the two noncommuting one-electron and two-electron parts of the time-dependent electronic Hamiltonian. We implement the scheme to calculate the laser-induced dynamics of the two-electron H2 molecule described within a minimal basis, and show how electron correlation is affected by the interaction of the molecule with a strong laser field.     

Quantum and electromagnetic propagation with the conjugate symmetric Lanczos method

The conjugate symmetric Lanczos (CSL) method is introduced for the solution of the time-dependent Schrodinger equation. This remarkably simple and efficient time-domain algorithm is a low-order polynomial expansion of the quantum propagator for time-independent Hamiltonians and derives from the time-reversal symmetry of the Schrodinger equation. The CSL algorithm gives forward solutions by simply complex conjugating backward polynomial expansion coefficients. Interestingly, the expansion coefficients are the same for each uniform time step, a fact that is only spoiled by basis incompleteness and finite precision. This is true for the Krylov basis and, with further investigation, is also found to be true for the Lanczos basis, important for efficient orthogonal projection-based algorithms. The CSL method errors roughly track those of the short iterative Lanczos method while requiring fewer matrix-vector products than the Chebyshev method. With the CSL method, only a few vectors need to be stored at a time, there is no need to estimate the Hamiltonian spectral range, and only matrix-vector and vector-vector products are required. Applications using localized wavelet bases are made to harmonic oscillator and anharmonic Morse oscillator systems as well as electrodynamic pulse propagation using the Hamiltonian form of Maxwell's equations. For gold with a Drude dielectric function, the latter is non-Hermitian, requiring consideration of corrections to the CSL algorithm.     

Charge and energy dynamics in photo-excited poly(para-phenylenevinylene) systems

We report results from simulations of charge and energy dynamics in poly(para-phenylenevinylene) (PPV) and PPV interacting with C60. The simulations were performed by solving the time-dependent Schrodinger equation and the lattice equation of motion simultaneously and nonadiabatically. The electronic system and the coupling of the electrons to the lattice were described by an extended three-dimensional version of the Su-Schrieffer-Heeger model, which also included an external electric field. Electron and lattice dynamics following electronic excitations at different energies have been simulated. The effect of additional lattice energy was also included in the simulations. Our results show that both exciton diffusion and transitions from high to lower lying excitations are stimulated by increasing the lattice energy. Also field induced charge separation occurs faster if the lattice energy is increased. This separation process is highly nonadiabatic and involves a significant rearrangement of the electron distribution. In the case of PPV coupled to C60, we observe a spontaneous charge separation. The separation time is in this case limited by the local concentration of C60 molecules close to the PPV chain.     

Two-electron reduced density matrices from the anti-Hermitian contracted Schrodinger equation: enhanced energies and properties with larger basis sets

Two-electron reduced density matrices (2-RDMs) have recently been directly determined from the solution of the anti-Hermitian contracted Schrodinger 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 Schrodinger 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.     

Gradient symplectic algorithms for solving the radial Schrodinger equation

The radial Schrodinger equation for a spherically symmetric potential can be regarded as a one-dimensional classical harmonic oscillator with a time-dependent spring constant. For solving classical dynamics problems, symplectic integrators are well known for their excellent conservation properties. The class of gradient symplectic algorithms is particularly suited for solving harmonic-oscillator dynamics. By use of Suzuki's rule [Proc. Jpn. Acad., Ser. B: Phys. Biol. Sci. 69, 161 (1993)] for decomposing time-ordered operators, these algorithms can be easily applied to the Schrodinger equation. We demonstrate the power of this class of gradient algorithms by solving the spectrum of highly singular radial potentials using Killingbeck's method [J. Phys. A 18, 245 (1985)] of backward Newton-Ralphson iterations.     

Symplectic splitting operator methods for the time-dependent Schrodinger equation

We present a family of symplectic splitting methods especially tailored to solve numerically the time-dependent Schrodinger equation. When discretized in time, this equation can be recast in the form of a classical Hamiltonian system with a Hamiltonian function corresponding to a generalized high-dimensional separable harmonic oscillator. The structure of the system allows us to build highly efficient symplectic integrators at any order. The new methods are accurate, easy to implement, and very stable in comparison with other standard symplectic integrators.     

Solving for the time-dependent multi-configuration hartree fock equations based on sine-discrete variable representation

A parallel quantum electrons wave packet computer code has been developed to study laser-atom interaction in the nonperturbative regime with attosecond resolution. The motion equations of the multi-configuration time-dependent hartree fock (MCTDHF) based on a sine discrete variable representation were solved by using an adaptive stepsize Runge-Kutta integrator of eight orders. Some efficient algorithms and strategies to accelerate the calculation velocity are introduced and discussed in details. Some illustrated imaginary time propagation and real time propagation have been respectively done in the paper. Single ionization probabilities calculated by using this one dimension MCTDHF model underestimate the accurate results calculated by solving time-dependent Schrodinger equation directly.     

Time-dependent density-functional theory beyond the adiabatic approximation: insights from a two-electron model system

Most applications of time-dependent density-functional theory (TDDFT) use the adiabatic local-density approximation (ALDA) for the dynamical exchange-correlation potential V(xc)(r,t). An exact (i.e., nonadiabatic) extension of the ground-state LDA into the dynamical regime leads to a V(xc)(r,t) with a memory, which causes the electron dynamics to become dissipative. To illustrate and explain this nonadiabatic behavior, this paper studies the dynamics of two interacting electrons on a two-dimensional quantum strip of finite size, comparing TDDFT within and beyond the ALDA with numerical solutions of the two-electron time-dependent Schrodinger equation. It is shown explicitly how dissipation arises through multiple particle-hole excitations, and how the nonadiabatic extension of the ALDA fails for finite systems but becomes correct in the thermodynamic limit.     

Spectral convergence of the quadrature discretization method in the solution of the Schrodinger and Fokker-Planck equations: comparison with sinc methods

Spectral methods based on nonclassical polynomials and Fourier basis functions or sinc interpolation techniques are compared for several eigenvalue problems for the Fokker-Planck and Schrodinger equations. A very rapid spectral convergence of the eigenvalues versus the number of quadrature points is obtained with the quadrature discretization method (QDM) and the appropriate choice of the weight function. The QDM is a pseudospectral method and the rate of convergence is compared with the sinc method reported by Wei [J. Chem. Phys., 110, 8930 (1999)]. In general, sinc methods based on Fourier basis functions with a uniform grid provide a much slower convergence. The paper considers Fokker-Planck equations (and analogous Schrodinger equations) for the thermalization of electrons in atomic moderators and for a quartic potential employed to model chemical reactions. The solution of the Schrodinger equation for the vibrational states of I2 with a Morse potential is also considered.     

用双激光脉冲操纵N2分子取向

利用求解含时薛定谔方程的方法, 研究了双原子N2分子在激光强度为1.5×10¹⁴ W·cm^-2, 脉冲宽度为50 fs激光脉冲作用下的取向行为. 研究结果表明, 保持总激光强度不变, 将一束激光脉冲分成具有同样脉宽, 强度比为0.3636的两束激光脉冲, 当第一束脉冲产生的分子取向即将达到最大时, 加上有一定延迟时间的第二束脉冲能够使N2分子达到最大的取向程度.     

Theory of three-dimensional alignment by intense laser pulses

We introduce a theoretical framework for study of three-dimensional alignment by moderately intense laser pulses and discuss it at an elementary level. Several features of formal interest are noted and clarified. Our approach is nonperturbative, treating the laser field within classical and the material system within quantum mechanics. The theory is implemented numerically using a basis set of rotational eigenstates, transforming the time-dependent Schrodinger equation to a set of coupled differential equations where all matrix elements are analytically soluble. The approach was applied over the past few years to explore different adiabatic and nonadiabatic three-dimensional alignment approaches in conjunction with experiments, but its formal details and numerical implementation were not reported in previous studies. Although we provide simple numerical examples to illustrate the content of the equations, our main goal is to complement previous reports through an introductory discussion of the underlying theory.     

A computationally efficient exact pseudopotential method. I. Analytic reformulation of the Phillips-Kleinman theory

Even with modern computers, it is still not possible to solve the Schrodinger equation exactly for systems with more than a handful of electrons. For many systems, the deeply bound core electrons serve merely as placeholders and only a few valence electrons participate in the chemical process of interest. Pseudopotential theory takes advantage of this fact to reduce the dimensionality of a multielectron chemical problem: the Schrodinger equation is solved only for the valence electrons, and the effects of the core electrons are included implicitly via an extra term in the Hamiltonian known as the pseudopotential. Phillips and Kleinman (PK) [Phys. Rev. 116, 287 (1959)]. demonstrated that it is possible to derive a pseudopotential that guarantees that the valence electron wave function is orthogonal to the (implicitly included) core electron wave functions. The PK theory, however, is expensive to implement since the pseudopotential is nonlocal and its computation involves iterative evaluation of the full Hamiltonian. In this paper, we present an analytically exact reformulation of the PK pseudopotential theory. Our reformulation has the advantage that it greatly simplifies the expressions that need to be evaluated during the iterative determination of the pseudopotential, greatly increasing the computational efficiency. We demonstrate our new formalism by calculating the pseudopotential for the 3s valence electron of the Na atom, and in the subsequent paper, we show that pseudopotentials for molecules as complex as tetrahydrofuran can be calculated with our formalism in only a few seconds. Our reformulation also provides a clear geometric interpretation of how the constraint equations in the PK theory, which are required to obtain a unique solution, are themselves sufficient to calculate the pseudopotential.     

三电子原子的超球坐标分析和基态能态计算

用超球坐标表示多电子原子的薛定谔方程,从符合Pauli原理的全反称波函数出发,推导出相应的超球径耦合微分方程,其形式与双电子原子的相似。以九维超球谐为投影函数构造S₃群不可约表示的对称基,将超球谐-广义Laguerre展开法(HHGL)推广到Li原子,首次对其Schrdinger方程的精确求解进行了研究,基态能量的计算值与实验值接近。最后,还从波函数出发对提高计算的精度进行了讨论。