An analytical approach for computing Franck‐Condon integrals of harmonic oscillators with arbitrary dimensions

We have developed an analytical approach for computing Franck‐Condon integrals (FCIs) of harmonic oscillators (HOs) with arbitrary dimensions in which the mode‐mixing Duschinsky effect is taken into account. A general formula of FCIs of HOs was obtained and was applied to study the photoelectron spectroscopy of vinyl alcohol and ovalene (C₃₂H₁₄). The equilibrium geometries, harmonic vibrational frequencies and normal modes of vinyl alcohol, ovalene, and their cations were computed at the B3LYP/aug‐cc‐pVTZ or the B3LYP/6‐31G(d) level, from which Franck‐Condon factors were calculated and photoelectron spectra were simulated. The adiabatic ionization energies of vinyl alcohol were also computed by extrapolating the CCSD(T) energies to the complete basis set limit with aug‐cc‐pVXZ (X = D, T, Q, 5). The simulated photoelectron spectra of both vinyl alcohol and ovalene are in agreement with the experiments. The computed adiabatic ionization energies of syn‐ and anti‐vinyl alcohol are in consistent with the experiment within 0.008 eV and 0.014 eV, respectively. We show, for the first time, that the analytical approach of computing FCIs is also efficient and promising for the studies of vibronic spectra of macrosystems. © 2012 Wiley Periodicals, Inc.     

Connection between the regular approximation and the normalized elimination of the small component in relativistic quantum theory

The regular approximation to the normalized elimination of the small component (NESC) in the modified Dirac equation has been developed and presented in matrix form. The matrix form of the infinite-order regular approximation (IORA) expressions, obtained in [Filatov and Cremer, J. Chem. Phys. 118, 6741 (2003)] using the resolution of the identity, is the exact matrix representation and corresponds to the zeroth-order regular approximation to NESC (NESC-ZORA). Because IORA (=NESC-ZORA) is a variationally stable method, it was used as a suitable starting point for the development of the second-order regular approximation to NESC (NESC-SORA). As shown for hydrogenlike ions, NESC-SORA energies are closer to the exact Dirac energies than the energies from the fifth-order Douglas-Kroll approximation, which is much more computationally demanding than NESC-SORA. For the application of IORA (=NESC-ZORA) and NESC-SORA to many-electron systems, the number of the two-electron integrals that need to be evaluated (identical to the number of the two-electron integrals of a full Dirac-Hartree-Fock calculation) was drastically reduced by using the resolution of the identity technique. An approximation was derived, which requires only the two-electron integrals of a nonrelativistic calculation. The accuracy of this approach was demonstrated for heliumlike ions. The total energy based on the approximate integrals deviates from the energy calculated with the exact integrals by less than 5 x 10(-9) hartree units. NESC-ZORA and NESC-SORA can easily be implemented in any nonrelativistic quantum chemical program. Their application is comparable in cost with that of nonrelativistic methods. The methods can be run with density functional theory and any wave function method. NESC-SORA has the advantage that it does not imply a picture change.     

Vector processing algorithm for electron repulsion integrals in Ab Initio HF calculation based upon the PK supermatrix

An efficient vector processing algorithm generating PK supermatrices has been developed, in particular aiming at calculations on large molecules. The algorithm utilizes the recurrence relations for electron repulsion integrals. The PK supermatrices are listed in a nearly canonical order so that the Fock matrix generation is efficiently vectorized, no temporary ERI and PK files being used. This is effected by partition of the basis set (atomic orbitals) into subsets of certain appropriate sizes, and the partition approach is named as the three-dimensional partial space method. A high-speed Hartree–Fock calculation including integrals and SCF procedures is achieved. © 1992 by John Wiley & Sons, Inc.     

Density functional calculations of potential energy surface and charge transfer integrals in molecular triphenylene derivative HAT6

We investigate the effect of structural fluctuations on charge transfer integrals, overlap integrals, and site energies in a system of two stacked molecular 2,3,6,7,10,11-hexakishexyloxytriphenylene (HAT₆), which is a model system for conducting devices in organic photocell applications. A density functional based computational study is reported. Accurate potential energy surface calculations are carried out using an improved meta-hybrid density functional to determine the most stable configuration of the two weakly bound HAT₆ molecules. The equilibrium parameters in terms of the twist angle and co-facial separation are calculated. Adopting the fragment approach within the Kohn–Sham density functional framework, these parameters are combined to a lateral slide, to mimic structural/conformational fluctuations and variations in the columnar phase. The charge transfer and spatial overlap integrals, and site energies, which form the matrix element of the Kohn–Sham Hamiltonian are derived. It is found that these quantities are strongly affected by the conformational variations. The spatial overlap between stacked molecules is found to be of considerable importance since charge transfer integrals obtained using the fragment approach differ significantly from those using the dimer approach.     

Auxiliary-field quantum Monte Carlo calculations of molecular systems with a Gaussian basis

We extend the recently introduced phaseless auxiliary-field quantum Monte Carlo (QMC) approach to any single-particle basis and apply it to molecular systems with Gaussian basis sets. QMC methods in general scale favorably with the system size as a low power. A QMC approach with auxiliary fields, in principle, allows an exact solution of the Schrodinger equation in the chosen basis. However, the well-known sign/phase problem causes the statistical noise to increase exponentially. The phaseless method controls this problem by constraining the paths in the auxiliary-field path integrals with an approximate phase condition that depends on a trial wave function. In the present calculations, the trial wave function is a single Slater determinant from a Hartree-Fock calculation. The calculated all-electron total energies show typical systematic errors of no more than a few millihartrees compared to exact results. At equilibrium geometries in the molecules we studied, this accuracy is roughly comparable to that of coupled cluster with single and double excitations and with noniterative triples [CCSD(T)]. For stretched bonds in H(2)O, our method exhibits a better overall accuracy and a more uniform behavior than CCSD(T).     

Ab-initio calculations on large molecules and solids by desirable computational procedures

Desirable computational procedures developed here recently for ab-initio calculations on large molecules are outlined. Effective core model potentials (MODPOT) permit calculations of valence electrons only explicitly, yet accurately; a charge-conserving integral prescreening evaluation to decide whether a block of integrals will be larger than a preset threshold and thus be calculated explicitly is effective for spatially extended systems; an efficient MERGE technique to save and reuse common invariant skeletal integrals is useful for geometry variations and for adding basis fcuntions, substituent groups and molecules; and an effective configuration interaction (CI) Hamiltonian into which are folded the effects of the occupied molecular orbitals from which no excitations are allowed is useful for molecular decompositions and intermolecular reactions. These techniques have been extended for CI calculations on breaking a chemical bond in a molecule in a crystal or solid; atom-class/atomic-class potential functions and dispersion calculations have been added. In a new program, POLY-CRYST, all the integral strategies for large molecules are meshed.     

General harmonic oscillator integrals by the operator method

The operator technique with a minimum of commutator algebra is employed to calculate matrix elements of any number of operators between distorted, displaced harmonic oscillator wavefunctions. The results are valid for multidimensional integrals, and regardless of the extent of the Duschinsky effect. General recursion relations useful in machine calculations are given. The formalism is illustrated for the well-known one-dimensional Franck–Condon integrals.     

Electrostatic energies and forces computed without explicit interparticle interactions: a linear time complexity formulation

A rapid method for the calculation of the electrostatic energy of a system without a cutoff is described in which the computational time grows linearly with the number of particles or charges. The inverse of the distance is approximated as a polynomial, which is then transformed into a function whose terms involve individual particles, instead of particle pairs, by a partitioning of the double sum. In this way, the electrostatic energy that is determined by the interparticle interactions is obtained without explicit calculation of these interactions. For systems of positive charges positioned on a face-centered cubic lattice, the calculation of the energy by the new method is shown to be faster than the calculation of the exact energy, in many cases by an order of magnitude, and to be accurate to within 1-2%. The application of this method to increase the accuracy of conventional truncation-based calculations in condensed-phase systems is also demonstrated by combining the approximated long-range electrostatic interactions with the exact short-range interactions in a "hybrid" calculation. For a 20-A sphere of water molecules, the forces are shown to be six times as accurate using this hybrid method as those calculated with conventional truncation of the electrostatic energy function at 12 A. This is accomplished with a slight increase in speed, and with a sevenfold increase in speed relative to the exact all-pair calculation. Structures minimized with the hybrid function are shown to be closer to structures minimized with an exact all-pair electrostatic energy function than are those minimized with a conventional 13-A cutoff-based electrostatic energy function. Comparison of the energies and forces calculated with the exact method illustrate that the absolute errors obtained with standard truncation can be very large. The extension of the current method to other pairwise functions as well as to multibody functions, is described.     

Ab initio molecular orbital calculations on large lattice cluster models: Use of translational symmetry

Summary Translational symmetry has been shown to be useful in the calculation of electronic structures of large lattice models. The number of unique integrals has been derived for cases of different dimensionality. For the unique integrals zero screening and approximation methods are described. The method has been applied to arrays of hydrogen atoms and to a zincblende surface model. When the size of the system is increased the translationally unique integrals are shown to become either zero or they can be calculated by simple coulombic approximations.     

Half-numerical evaluation of pseudopotential integrals

A half-numeric algorithm for the evaluation of effective core potential integrals over Cartesian Gaussian functions is described. Local and semilocal integrals are separated into two-dimensional angular and one-dimensional radial integrals. The angular integrals are evaluated analytically using a general approach that has no limitation for the l-quantum number. The radial integrals are calculated by an adaptive one-dimensional numerical quadrature. For the semilocal radial part a pretabulation scheme is used. This pretabulation simplifies the handling of radial integrals, makes their calculation much faster, and allows their easy reuse for different integrals within a given shell combination. The implementation of this new algorithm is described and its performance is analyzed.     

Damped response theory description of two-photon absorption

Damped response theory is applied to the calculation of two-photon absorption (TPA) spectra, which are determined directly, at each frequency, from a modified damped cubic response function. The TPA spectrum may therefore be evaluated for selected frequency ranges, making the damped TPA approach attractive for calculations on large molecules with a high density of states, where the calculation of TPA using standard theory is more problematic. Damped response theory can also be applied to the case of intermediate state resonances, where the standard TPA expression is divergent. Both exact damped response theory and its application within density functional theory are discussed. The latter is implemented using an atomic-orbital based density matrix formulation, which makes the approach especially suitable for studies on large systems. A test preliminary study is presented for the TPA spectrum of R-(+)-1,1'-bi(2-naphtol).     

Promoting-mode free formalism for excited state radiationless decay process with Duschinsky rotation effect

In the present work, through the path integral of Gaussian type correlation function, a new formalism based on Fermi-Golden Rule for calculating the rate constant of nonradiative decay process with Duschinsky rotation effect in polyatomic molecules is developed. The advantage of the present path-integral formalism is promoting-mode free. In order to get the rate constant, a "transition rate matrix" needs to be calculated. The rate constant calculated previously is only an approximation of diagonal elements of our "transition rate matrix " . The total rate should be the summation over all the matrix elements.     

An unconventional scf method for calculations on large molecules

An unconventional SCF method for calculations on large molecules with more than 100 basis functions is described. Storage problems which arise in conventional SCF schemes when storing more than 10⁷ integrals are avoided by repeated calculation of integrals. The resulting increase in computational times is kept at a reasonable level by (a) improving the initial guess, (b) accelerating convergence, (c) employing a recursive construction of the Fock matrix, and (d) eliminating insignificant integrals from the calculation by a density-weighted cutoff criterion. Sample calculations show that, compared with conventional SCF calculations, computational times increase by 25%–75% depending on the basis set and the shape of the molecule.     

Møller-Plesset (MP2) perturbation theory for large molecules

Summary A novel formulation of MP2 theory is presented which starts from the Laplace transform MP2 ansatz, and subsequently moves from a molecular orbital (MO) representation to an atomic orbital (AO) representation. Consequently, the new formulation is denoted AO-MP2. As in traditional MP2 approaches electron repulsion integrals still need to be transformed. Strict bounds on the individual MP2 energy contribution of each intermediate four-index quantity allow to screen off numerically insignificant integrals with a single threshold parameter. Implicit in our formulation is a bound to two-particle density matrix elements. For small molecules the computational cost for AO-MP2 calculations is about a factor of 100 higher than for traditional MO-based approaches, but due to screening the computational effort in larger systems will only grow with the fourth power of the size of the system (or less) as is demonstrated both in theory and in application. MP2 calculations on (non-metallic) crystalline systems seem to be a feasible extension of the Laplace transform approach. In large molecules the AO-MP2 ansatz allows massively parallel MP2 calculations without input/output of four-index quantities provided that each processor has in-core memory for a limited number of two-index quantities. Energy gradient formulas for the AO-MP2 approach are derived.Dedicated to Prof. W. Kutzelnigg whose books on theoretical chemistry aroused my interest in this field     

Quality of approximate electron densities and internal consistency of molecular alignment algorithms in molecular quantum similarity

The calculation of molecular quantum similarity measures using the molecular electron density requires the electron density and molecular alignment between two molecules. To obtain meaningful quantum similarity matrices, the electron density should be calculated efficiently and accurately and the alignment should be internally consistent. The internal consistency of the alignment for a series of molecules is investigated through distance geometry concepts. The calculation of the quantum similarity matrix requires the calculation of a quadratic number of similarity integrals, and a scheme to obtain these efficiently is developed. Both the alignment procedure and the ASA method for approximate molecular electron densities are tested for a set of steroid molecules.     

Density fitting of two-electron integrals in extended systems with translational periodicity: the Coulomb problem

Density fitting approach to Coulomb integrals for infinite systems with translational periodicity is reformulated in direct space. Despite of the Coulomb infinite decay of some integrals, direct-space calculation is shown to be feasible. Moreover, we show that the direct-space ansatz is completely equivalent to our previous formulation in reciprocal space. Computational demands scale linearly with the number of unit cells. In addition, direct-space treatment has some practical advantages over the reciprocal-space formulation. The efficiency of our scheme is demonstrated on systems with translational periodicity in one dimension. Computation time takes only a small fraction of the conventional calculation with exact integrals. We show that for infinite systems auxiliary basis sets of equally good quality as for molecules can be constructed in a systematic way.     

Unified analytical treatment for calculation of the two‐dimensional Franck–Condon factors using the duschinsky transformation

Including binomial expansion theorems, we present an analytical formula for calculating Franck–Condon (FC) factors of two‐dimensional (2D) harmonic oscillators including the Duschinsky effect. The FC principle has various practical applications in quantum modeling of electronic spectra of polyatomic molecules. The 2D FC factors are expressed through the binomial coefficients. Use of the memory of the computer for the calculation of binomial coefficients may extend the limits to large arguments for users and result in speeder calculation, should such limits be required in practice. Accurate numerical results are provided to validate the proposed algorithm. © 2012 Wiley Periodicals, Inc.     

Perturbative calculation of Franck-Condon integrals: new hints for a rational implementation

Perturbative approaches to the mode mixing effects in the calculation of Franck-Condon integrals are analyzed and discussed. The zero order multidimensional Franck-Condon integrals are factorized into products of one-dimensional ones, so that recurrence relations can be used without need of storing a huge number of data. Calculations on model systems show that at the second order of perturbation, the method gives results in very good agreement with the exact ones, even in the case of significantly large Duschinsky effect. The accuracy of the results can be substantially improved by grouping together all those modes which are strongly mixed with each other, usually a few ones, for which Franck-Condon integrals can be computed exactly, and using the perturbative approach for treating the smaller mixing between all the other modes.     

Excited state radiationless decay process with Duschinsky rotation effect: formalism and implementation

Duschinsky rotation effect is a simple and effective way to characterize the difference between the ground state and excited state potential energy surfaces. For complex molecules, harmonic oscillator model is still the practical way to describe the dynamics of excited states. Based on the first-order perturbation theory a la Fermi golden rule, the authors have applied the path integral of Gaussian type for the correlation function to derive an analytic formalism to calculate the internal conversion rate process with Duschinsky rotation effect being taken into account. The validity of their formalism is verified through comparison with previous work, both analytically for the case of neglecting Duschinsky rotation and numerically for the ethylene molecules with two-mode mixing. Their expression is derived for multimode mixing.