Ab initio correlation corrections to the Hartree-Fock quasi band-structure of periodic systems employing Wannier-type orbitals

A size-consistent ab initio formalism to calculate correlation corrections to ionization potentials as well as electron affinities of periodic systems is presented. Our approach is based on a Hartree-Fock scheme which directly yields local orbitals without any a posteriori localization step. The use of local orbitals implies non-zero off-diagonal matrix elements of the Fock operator, which are treated as an additional perturbation and give rise to localization diagrams. Based on the obtained local orbitals, an effective Bloch Hamiltonian is constructed to second order of perturbation theory with all third-order localization diagrams included. In addition, the summation of certain classes of diagrams up to infinite order in the off-diagonal Fock elements as well as the Epstein-Nesbet partitioning of the full Hamiltonian are discussed. The problem of intruder states, frequently encountered in many-body perturbation theory, is dealt with by employing the theory of intermediate Hamiltonians. As model systems we have chosen cyclic periodic structures up to an oligoethylene ring in double-zeta basis; however, the theory presented here straightforwardly carries over to infinite periodic systems. Received: 30 April 1998 / Accepted: 27 July 1998 /  Published online: 7 October 1998     

Potential energy surfaces and dynamics of Ni2+ ion aqueous solution: molecular dynamics simulation of the electronic absorption spectrum

We develop a model effective Hamiltonian for describing the electronic structures of first-row transition metals in aqueous solutions using a quasidegenerate perturbation theory. All the states consisting of 3d(n) electronic configurations are determined by diagonalizing a small effective Hamiltonian matrix, where various intermolecular interaction terms such as the electrostatic, polarization, exchange, charge transfer, and three-body interactions are effectively incorporated. This model Hamiltonian is applied to constructing the ground and triplet excited states potential energy functions of Ni(2+) in aqueous solution, based on the ab initio multiconfiguration quasidegenerate perturbation theory calculations. We perform molecular dynamics simulation calculations for the ground state of Ni(2+) aqueous solution to calculate the electronic absorption spectral shape as well as the ground state properties. Agreement between the simulation and experimental spectra is satisfactory, indicating that the present model can well describe the Ni(2+) excited state potential surfaces in aqueous solution.     

Theory of radiationless multiphonon transitions in the electron shells of molecules

The probability of a multiphonon transition is calculated in the second approximation of the nonadiabatic perturbation theory. Allowance is made for resonance of the adiabatic potentials and for the change in the frequencies of the normal modes resulting from the transition.Read at the Symposium on Quantum Chemistry, Palanga, June 1965.     

Different forms of perturbation theory for the calculation of the correlation energy

Different forms of perturbation theory for the calculation of correlation energy in both closed-and open-shell systems are discussed. For closed-shell systems, Epstein–Nesbet perturbation theory is compared with Møller–Plesset (MP ) perturbation theory based on canonical Hartree–Fock orbitals and with MP theory based on internally consistent SCF orbitals. The traditional MP theory gives superior results despite its use of an inferior zeroth-order Hamiltonian. This behavior is rationalized in terms of the larger denominators present in the traditional MP theory. These conclusions are used to support the restricted open-shell perturbation methods proposed recently by Murray and Davidson, and these new methods are compared with spin-restricted Epstein–Nesbet theory and the unrestricted MP (UMP ) approach. © 1992 John Wiley & Sons, Inc.     

Calculation of vibrational transition frequencies and intensities in water dimer: comparison of different vibrational approaches

We have calculated frequencies and intensities of fundamental and overtone vibrational transitions in water and water dimer with use of different vibrational methods. We have compared results obtained with correlation-corrected vibrational self-consistent-field theory and vibrational second-order perturbation theory both using normal modes and finally with a harmonically coupled anharmonic oscillator local mode model including OH-stretching and HOH-bending local modes. The coupled cluster with singles, doubles, and perturbative triples ab initio method with augmented correlation-consistent triple-zeta Dunning and atomic natural orbital basis sets has been used to obtain the necessary potential energy and dipole moment surfaces. We identify the strengths and weaknesses of these different vibrational approaches and compare our results to the available experimental results.     

Improved perturbative treatment of electronic energies from a minimal-norm approach to many-body perturbation theory

We propose a zeroth-order Hamiltonian for many-body perturbation theory based on the unitary decomposition of the two-particle reduced Hamiltonian. For the zeroth-order Hamiltonian constrained to be diagonal in the Hartree-Fock basis set, the two-particle reduced perturbation matrix is chosen to have a minimal Frobenius norm. When compared with the M?ller-Plesset partitioning, the method yields more accurate second-order energies.     

Femtosecond primary events in bacteriorhodopsin and its retinal modified analogs: revision of commonly accepted interpretation of electronic spectra of transient intermediates in the bacteriorhodopsin photocycle

Femtosecond primary events in bacteriorhodopsin (BR) and its retinal modified analogs are discussed. Ultrafast time resolved electronic spectra of the primary intermediates induced in the BR photocycle are discussed along with spectral and kinetic inconsistencies of the previous models proposed in the literature. The theoretical model proposed in this paper based on vibrational coupling between the electronic transition of the chromophore and intramolecular vibrational modes allows us to calculate the equilibrium electronic absorption band shape and the hole burning profiles. The model is able to rationalize the complex pattern of behavior for the primary events in BR and explain the origin of the apparent inconsistencies between the experiment and the previous theoretical models. The model presented in the paper is based on the anharmonic coupling assumption in the adiabatic approximation using the canonical transformation method for diagonalization of the vibrational Hamiltonian instead of the commonly used perturbation theory. The electronic transition occurs between the Born-Oppenheimer potential energy surfaces with the electron involved in the transition being coupled to the intramolecular vibrational modes of the molecule (chromophore). The relaxation of the excited state occurs by indirect damping (dephasing) mechanisms. The indirect dephasing is governed by the time evolution of the anharmonic coupling constant driven by the resonance energy exchange between the intramolecular vibrational mode and the bath. The coupling with the intramolecular vibrational modes results in the Franck-Condon progression of bands that are broadened due to the vibrational dephasing mechanisms. The electronic absorption line shape has been calculated based on the linear response theory whereas the third order nonlinear response functions have been used to analyze the hole burning profiles obtained from the pump-probe time-resolved measurements. The theoretical treatment proposed in this paper provides a basis for a substantial revision of the commonly accepted interpretation of the primary events in the BR photocycle that exists in the literature.     

Tunneling and quantum localization in chaos‐driven symmetric triple well potential: An approach using quantum theory of motion

For a symmetric triple well potential, driven by the forces associated with the bifurcation diagram of a logistic map, the tunneling and quantum localization are studied using quantum theory of motion and time‐dependent Fourier grid Hamiltonian methods. Detailed analysis reveals that application of only asymmetric or symmetric perturbation results into either quantum localization or over‐barrier transition and no tunneling while application of mixed symmetry perturbation gives either tunneling or over‐barrier transition, depending on temporal nature and initial position of the particle. For bifurcative and chaotic symmetric‐asymmetric perturbation, with truncation of mixed symmetry perturbation, a sudden jump in energy causes a transition from the tunneling phenomenon to the over‐barrier transition. With particle located initially near to either of the minima of the unperturbed well, quantum localization, or over‐barrier transition is observed, depending on types of perturbation used.     

Effective Floquet Hamiltonian for spin I = 1 in magic angle spinning NMR using contact transformation

Contact transformation is an operator transformation method in time-independent perturbation theory which is used successfully in molecular spectroscopy to obtain an effective Hamiltonian. Floquet theory is used to transform the periodic time-dependent Hamiltonian, to a time-independent Floquet Hamiltonian. In this article contact transformation method has been used to get the analytical representation of Floquet Hamiltonian for quadrupolar nuclei with spin I = 1 in the presence of an RF field and first order quadrupolar interaction in magic angle spinning NMR experiments. The eigenvalues of contact transformed Hamiltonian as well as Floquet Hamiltonian have been calculated and a comparison is made between the eigenvalues obtained using the two Hamiltonians.     

Calculation of nonadiabatic couplings with restricted open-shell Kohn-Sham density-functional theory

This paper generalizes the recently proposed approaches for calculating the derivative couplings between adiabatic states in density-functional theory (DFT) based on a Slater transition-state density to transitions such as singlet-singlet excitations, where a single-determinant ansatz is insufficient. The proposed approach is based on restricted open-shell Frank et al. [J. Chem. Phys. 108, 4060 (1998)] theory used to describe a spin-adapted Slater transition state. To treat the dependence of electron-electron interactions on the nuclear positions, variational linear-response density-functional perturbation theory is generalized to reference states with an orbital-dependent Kohn-Sham Hamiltonian and nontrivial occupation patterns. The methods proposed in this paper are not limited to the calculation of derivative coupling vectors, but can also be used for the calculation of other transition matrix elements. Moreover, they can be used to calculate the linear response of open-shell systems to arbitrary external perturbations in DFT.     

Simulation of inversion motion and N-H stretching overtone spectra of aniline

A curvilinear internal coordinate Hamiltonian is used to simulate the N-H stretching overtone spectra and the associated inversion splittings in aniline. A simple local mode type model is applied to the N-H stretching and H-N-H bending modes. Geometric algebra is employed to derive the kinetic energy operator for the large amplitude inversion motion. Electronic structure calculations at the Moller-Plesset second order perturbation theory and correlation consistent aug-cc-pVTZ basis set level are used to obtain model parameters, some of which have been optimized with the least-squares method using experimental vibrational term values as data. The observed N-H stretching overtone vibrational levels and the inversional tunneling splittings are well reproduced with our approach.     

The OH-stretching and OOH-bending overtone spectrum of HOONO

We have simulated the HOONO vibrational overtone spectrum with use of a local mode Hamiltonian that includes the OH-stretching, OOH-bending, and NOOH-torsional modes and coupling between all three modes. The local mode parameters and the dipole moment function are calculated with coupled-cluster ab initio theory and an augmented Dunning-type triple-zeta basis set. We investigate the accuracy of the local mode parameters obtained from two different potential-energy fitting routines, as well as the sensitivity of these parameters to the basis set employed. We compare our simulated spectra to previously published action spectra in the first and second OH-stretching overtone regions. In addition we have recorded the spectrum in the OH-stretch and OOH-bend combination region around 7700 cm-1 and we also compare to this. Our simulated spectrum is in qualitative agreement with experiment in the first and second OH-stretching overtone and in the stretch-bend regions.     

Bending dynamics of acetylene: new modes born in bifurcations of normal modes

Semiclassical techniques are used to analyze highly excited pure bending vibrational dynamics from spectra of C2H2. An analytic bifurcation approach is developed, based on critical points of a classical version of the quantum fitting Hamiltonian. At high energy four new types of anharmonic modes are born in bifurcations of the normal modes: local, orthogonal, precessional, and counter-rotator. Visual insight into their nature is obtained with the help of computer-generated three-dimensional animations. The connection between the local mode and the acetylene-vinylidene isomerization "reaction mode" is considered.     

The calculation of vibrational energy levels of polyatomic molecules including anharmonic effect using contact transformation perturbation method

Using contact transformation perturbation method based on the Taylor expansion of the potential energy function in terms of dimensionless normal coordinates up to sixth‐order, the vibrational energy levels in terms of force constants are derived. The contact transformation theory has been applied to simplify the calculation of perturbation effects. To calculate the second‐order vibrational energy correction, the third and fourth‐order terms of potential function have been placed in the first‐order perturbation Hamiltonian and the second‐order Hamiltonian contains hexatic ones. We present expressions which give relations between the fourth‐ and sixth‐order terms in dimensionless normal coordinates of the potential and the anharmonicity coefficients. For illustration, a set of vibrational energies levels of SO₂, and H₂O molecules including anharmonic effects has been calculated. © 2013 Wiley Periodicals, Inc.     

Numerical-analytic implementation of the higher-order canonical Van Vleck perturbation theory for the interpretation of medium-sized molecule vibrational spectra

Anharmonic vibrational states of semirigid polyatomic molecules are often studied using the second-order vibrational perturbation theory (VPT2). For efficient higher-order analysis, an approach based on the canonical Van Vleck perturbation theory (CVPT), the Watson Hamiltonian and operators of creation and annihilation of vibrational quanta is employed. This method allows analysis of the convergence of perturbation theory and solves a number of theoretical problems of VPT2, e.g., yields anharmonic constants y(ijk), z(ijkl), and allows the reliable evaluation of vibrational IR and Raman anharmonic intensities in the presence of resonances. Darling-Dennison and higher-order resonance coupling coefficients can be reliably evaluated as well. The method is illustrated on classic molecules: water and formaldehyde. A number of theoretical conclusions results, including the necessity of using sextic force field in the fourth order (CVPT4) and the nearly vanishing CVPT4 contributions for bending and wagging modes. The coefficients of perturbative Dunham-type Hamiltonians in high-orders of CVPT are found to conform to the rules of equality at different orders as earlier proven analytically for diatomic molecules. The method can serve as a good substitution of the more traditional VPT2.     

A nondiagonal quasidegenerate fourth-order perturbation theory

A canonical quasidegenerate Rayleigh-Schrödinger perturbation theory, correct through fourth order in the energy, is explored for a block-diagonal unperturbed Hamiltonian. The theory is developed completely within a Lie Algebra in Hilbert space. Explicit equations forn-particle transition elements in terms of solutions of simultaneous linear equations are presented. A two-dimensional anisotropic anharmonic oscillator is used to provide numerical results. The perturbation theory is shown to be stable under small separation of model and complement spaces. An iterative variant of the fourth order perturbation theory is introduced; the iterative variant is related to the non-iterative one in much the same way as nondegenerate coupled cluster theories are related to nondegenerate perturbation theory. The quasidegenerate coupled cluster theory appears to be stable in the presence of multiple intruder states.     

Including anharmonicity in the calculation of rate constants. 1. The HCN/HNC isomerization reaction

A method for calculating anharmonic vibrational energy levels in asymmetric top and linear systems that is based on second-order perturbation theory in curvilinear coordinates is extended to the bound generalized normal modes at nonstationary points along a reaction path. Explicit formulas for the anharmonicity coefficients, x(ij), and the constant term, E0, are presented, and the necessary modifications for resonance cases are considered. The method is combined with variational transition state theory with semiclassical multidimensional tunneling approximations to calculate thermal rate constants for the HCN/HNC isomerization reaction. Although the results for this system are not very sensitive to the choice of coordinates, we find that the inclusion of anharmonicity leads to a substantial improvement in the vibrational energy levels. We also present detailed comparisons of rate constants computed with and without anharmonicity, with various approximations for incorporating tunneling along the reaction path, and with a more practical approach to calculating the vibrational partition functions needed for larger systems.     

Quantum simulation of solution phase intramolecular electron transfer rates in betaine-30

Mixed quantum-classical atomistic simulations have been carried out to investigate the mechanistic details of excited state intramolecular electron transfer in a betaine-30 molecule in acetonitrile. The key electronic degrees of freedom of the solute molecule are treated quantum mechanically using the semiempirical Pariser-Parr-Pople Hamiltonian, including the solvent influence on electronic structure. The intramolecular vibrational modes are also treated explicitly at a quantum level, with the remaining elements treated classically using empirical potentials. The electron-transfer rate, corresponding to S1 --> S0 relaxation, is evaluated via time-dependent perturbation theory with the explicit inclusion of the dynamics of solvation and intramolecular conformation. The calculations reveal that, while solvation dynamics is critical to the rate, the intramolecular torsional dynamics also plays an important role. The importance of the use of multiple high-frequency quantum modes is also discussed.     

Optimal diabatic bases via thermodynamic bounds

Describing kinetic processes within a perturbation theory approach such as Fermi's golden rule requires an understanding of the initial and final states of the system. A number of different methods have been proposed for obtaining these diabatic-like states, but a robust criterion for evaluating their accuracy has not been established. Here, we approach the problem of determining the most appropriate set of diabatic states for use in incoherent rate expressions. We develop a method that rotates an initial set of diabats into an optimized set beginning with a zeroth-order diabatic Hamiltonian and choosing the rotation that minimizes the effect of non-diabatic terms on the thermodynamic free energy. The Gibbs-Bogoliubov (GB) bound on the Helmholtz free energy is thus used as the diabatic criterion. We first derive the GB free energy for a two site system and then find an expression general for any electronic system Hamiltonian. Efficient numerical methods are used to perform the minimization subject to orthogonality constraints, and we examine the resulting diabats for system Hamiltonians in various parameter regimes. The transition from localized to delocalized states is clearly seen in these calculations, and some interesting features are discussed.     

Multi‐state multi‐reference Møller–Plesset second‐order perturbation theory for molecular calculations

This work presents multi‐state multi‐reference Møller–Plesset second‐order perturbation theory as a variant of multi‐reference perturbation theory to treat electron correlation in molecules. An effective Hamiltonian is constructed from the first‐order wave operator to treat several strongly interacting electronic states simultaneously. The wave operator is obtained by solving the generalized Bloch equation within the first‐order interaction space using a multi‐partitioning of the Hamiltonian based on multi‐reference Møller–Plesset second‐order perturbation theory. The corresponding zeroth‐order Hamiltonians are nondiagonal. To reduce the computational effort that arises from the nondiagonal generalized Fock operator, a selection procedure is used that divides the configurations of the first‐order interaction space into two sets based on the strength of the interaction with the reference space. In the weaker interacting set, only the projected diagonal part of the zeroth‐order Hamiltonian is taken into account. The justification of the approach is demonstrated in two examples: the mixing of valence Rydberg states in ethylene, and the avoided crossing of neutral and ionic potential curves in LiF. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006