Ab initio calculations of electronically excited states of cyano-substituted polyacetylene cations

The transition energies for the lowest energy pi --> pi* electronic excitations are calculated with the complete active space self-consistent field method (CASSCF) and with the complete active space second-order perturbation theory method (CASPT2) for the linear cyano-substituted polyacetylene cations, H-Cn-CN+, n = 4-11, and NC-Cn-CN+, n = 2-10. These systems are models for an important class of interstellar species. We demonstrate the utility of the theoretical calculations in assigning the experimental spectra.     

Following the excited state relaxation dynamics of indole and 5-hydroxyindole using time-resolved photoelectron spectroscopy

Time-resolved photoelectron spectroscopy was used to obtain new information about the dynamics of electronic relaxation in gas-phase indole and 5-hydroxyindole following UV excitation with femtosecond laser pulses centred at 249 nm and 273 nm. Our analysis of the data was supported by ab initio calculations at the coupled cluster and complete-active-space self-consistent-field levels. The optically bright (1)L(a) and (1)L(b) electronic states of (1)ππ? character and spectroscopically dark and dissociative (1)πσ? states were all found to play a role in the overall relaxation process. In both molecules we conclude that the initially excited (1)L(a) state decays non-adiabatically on a sub 100 fs timescale via two competing pathways, populating either the subsequently long-lived (1)L(b) state or the (1)πσ? state localised along the N-H coordinate, which exhibits a lifetime on the order of 1 ps. In the case of 5-hydroxyindole, we conclude that the (1)πσ? state localised along the O-H coordinate plays little or no role in the relaxation dynamics at the two excitation wavelengths studied.     

Ab initio mass tensor molecular dynamics

Mass tensor molecular dynamics method was first introduced by Bennett [J. Comput. Phys. 19, 267 (1975)] for efficient sampling of phase space through the use of generalized atomic masses. Here, we show how to apply this method to ab initio molecular dynamics simulations with minimal computational overhead. Test calculations on liquid water show a threefold reduction in computational effort without making the fixed geometry approximation. We also present a simple recipe for estimating the optimal atomic masses using only the first derivatives of the potential energy.     

On the electronically excited states of uracil

Vertical excitation energies in uracil in the gas phase and in water solution are investigated by the equation-of-motion coupled-cluster and multireference configuration interaction methods. Basis set effects are found to be important for converged results. The analysis of electronic wave functions reveals that the lowest singlet states are predominantly of a singly excited character and are therefore well described by single-reference equation-of-motion methods augmented by a perturbative triples correction to account for dynamical correlation.Our best estimates for the vertical excitation energies for the lowest singlet n --> pi* and pi --> pi* are 5.0 +/- 0.1 eV and 5.3 +/- 0.1 eV, respectively. The solvent effects for these states are estimated to be +0.5 eV and +/- 0.1 eV, respectively. We attribute the difference between the computed vertical excitations and the maximum of the experimental absorption to strong vibronic interaction between the lowest A" and A' states leading to intensity borrowing by the forbidden transition.     

Ultrafast nonradiative dynamics in electronically excited hexafluorobenzene by femtosecond time-resolved mass spectrometry

The fast nonradiative decay dynamics of the lowest two excited pipi(*) electronic states (S(2) and S(3)) of hexafluorobenzene have been investigated by using femtosecond time-resolved time-of-flight mass spectrometry. The molecules were excited at wavelengths between 265 nm > or = lambda(pump) > or = 217 nm and probed by four- and three-photon ionization at lambda(probe)=775 nm. The observed temporal profiles exhibit two exponential decay times (tau(1)=0.54-0.1 ps and tau(2)=493-4.67 ps, depending on the excitation wavelength) and a superimposed coherent oscillation with vibrational frequency nu(osc)=97 cm(-1) and damping time tau(D) that is two to three times longer than the respective tau(1). The first decay component (tau(1)) is assigned to rapid radiationless transfer from the excited optically bright pipi(*) electronic state (S(2) or S(3), respectively) through a conical intersection (CI) to the lower-lying optically dark pisigma(*) state (S(1)) of the molecule; the second component (tau(2)) is attributed to the subsequent slower relaxation from the S(1) state back to the electronic ground state (S(0)). tau(2) dramatically decreases with increasing vibronic excitation energy up to the CI connecting the pisigma(*) with the S(0) state. The coherent oscillation is identified as nuclear motion along the out-of-plane vibration nu(16a) (notation as for benzene), which has e(2u) symmetry and acts as coupling mode between the pipi(*) and pisigma(*) states.     

Ab initio molecular dynamics simulation of ionic liquids

Ab initio Car-Parinnello molecular dynamics is used to simulate the structure and the dynamics of 1-butyl-3-methylimidazolium iodide ([bmim]I) ionic liquid at 300 K. Site-site pair correlation functions reveal that the anion has a strong interaction with any three C-H's of the imidazolium ring. The ring bends over and wraps around the anion such that the two nitrogen atoms take a distance to the anion. Electron donating butyl group contributes the electronic polarization in addition to geometrical (out-of-plane) polarization of the ring due to the liquid environment. This facilitates bending of the ring along the axis passing through nitrogen atoms. The average bending angle depends largely on the alkyl chain length and slightly on the anion type. Redistribution of electron density over the ring caused by the electron donating alkyl group provides additional independent evidence to the instability of lattice structure, hence the low melting point of the ionic liquid. Simulated viscosity and diffusion coefficients of [bmim]I are in quite agreement with the experiments.     

Ab initio molecular dynamics study of liquid methanol

We present a density-functional theory based molecular-dynamics study of the structural, dynamical, and electronic properties of liquid methanol under ambient conditions. The calculated radial distribution functions involving the oxygen and hydroxyl hydrogen show a pronounced hydrogen bonding and compare well with recent neutron diffraction data, except for an underestimate of the oxygen–oxygen correlation. We observe that, in line with infrared spectroscopic data, the hydroxyl stretching mode is significantly red-shifted in the liquid. A substantial enhancement of the dipole moment is accompanied by significant fluctuations due to thermal motion. Our results provide valuable data for improvement of empirical potentials.     

A gas phase ab initio excited state geometry optimization study of thymine, cytosine and uracil

Molecular geometries of the nucleic acid bases thymine, cytosine and uracil in the ground and the lowest two singlet excited states were optimized using the ab initio approach employing the 4-31G basis set for all the atoms except the amino group of cytosine for which the 6-311+G^* basis set was used. The excited state calculations were performed employing configuration interaction involving singly excited configurations (CIS). Vibrational frequencies were computed in order to examine the nature of the stationary points on the potential energy surfaces obtained by geometry optimization. While the ground state geometries of uracil and thymine (except the methyl group hydrogens) are planar, the corresponding excited state geometries were found to be significantly nonplanar. In the case of cytosine, the amino group is pyramidal and the rest of the molecule is only slightly nonplanar in the ground state, but the excited state geometries are appreciably nonplanar. In particular, consequent to the S₂(n–π^*) excitation of cytosine, the amino group plane is strongly rotated. While thymine is stable in the S₂(π–π^*) excited state, uracil appears to be dissociative in the corresponding excited state.     

Ab initio study of molecular structures and excited states in anthocyanidins

The structural and electronic characters of four types of hydroxyl group-substituted anthocyanidins (pelargonidin, cyanidin, delphinidin, and aurantinidin) were examined using quantum chemical calculations. For these cationic molecules, both the planar and non-planar structures in the electronic ground state were determined at the B3LYP/D95 level of theory. We revealed that the planar structure is slightly more stable than the non-planar structure for each molecule. For the optimized planar structures, single excitation-configuration interaction (SE-CI) based on the restricted Hartree-Fock (RHF) wave function was evaluated and the electronic character in the low-excited states was discussed in terms of the MO theory. Symmetry adapted cluster (SAC)/SAC-CI calculations were also carried out to estimate the excitation energies precisely. The results showed that hydroxylation of the phenyl group causes a change in the excitation energies without taking the solvent effects into account. The results are in agreement with spectral experiments and previous MO calculations.     

Ab initio molecular dynamics using hybrid density functionals

Ab initio molecular dynamics simulations with hybrid density functionals have so far found little application due to their computational cost. In this work, an implementation of the Hartree-Fock exchange is presented that is specifically targeted at ab initio molecular dynamics simulations of medium sized systems. We demonstrate that our implementation, which is available as part of the CP2K/Quickstep program, is robust and efficient. Several prescreening techniques lead to a linear scaling cost for integral evaluation and storage. Integral compression techniques allow for in-core calculations on systems containing several thousand basis functions. The massively parallel implementation respects integral symmetry and scales up to hundreds of CPUs using a dynamic load balancing scheme. A time-reversible multiple time step scheme, exploiting the difference in computational efficiency between hybrid and local functionals, brings further time savings. With extensive simulations of liquid water, we demonstrate the ability to perform, for several tens of picoseconds, ab initio molecular dynamics based on hybrid functionals of systems in the condensed phase containing a few thousand Gaussian basis functions.     

Ab initio molecular dynamics calculations to study catalysis

The modern versions of the density functional theory (DFT), especially those using the generalized gradient approximation (GGA), have reached (almost) chemical accuracy and thus can be applied to study problems of real chemical interest such as catalysis. The important equations for the DFT, the local density approximation (LDA), and GGA are given. The full-potential linearized augmented plane wave method (LAPW) is used to check the accuracy of GGA in solids. The basic concepts of the ab initio molecular dynamics (MD) method by Car and Parrinello and its recent implementation using the projector augmented Wave (PAW) method which use a similar augmentation as LAPW are described. PAW applications to ferrocene and beryllocene are summarized, which illustrate that vibrational or fluxional behavior are well described. Sodalite, an aluminosilicate, is discussed as a generic zeolite in comparison with gmelinite. A study of the dynamics of such a system allows the determination of, e.g., the proton stretch vibrations which can be related to infrared spectra. This is illustrated for the OH stretch vibration of the acid site in silicon-rich sodalite. With this methodology, we are able to study the interaction of methanol trapped inside the cage structure of silicon-rich sodalite and to gain new insight into crucial steps of catalytic reactions, namely, the hydrogen-bonding and the possible protonation in this system, or a proton-exchange reaction. The strategies for parallelizing the PAW code are outlined. © 1997 John Wiley & Sons, Inc.     

Ab initio molecular dynamics study of formate ion hydration

We perform ab initio molecular dynamics simulations of the aqueous formate ion. The mean number of water molecules in the first solvation shell, or the hydration number, of each formate oxygen is found to be consistent with recent experiments. Our ab initio pair correlation functions, however, differ significantly from many classical force field results and hybrid quantum mechanics/molecular mechanics predictions. They yield roughly one less hydrogen bond between each formate oxygen and water than force field or hybrid methods predict. Both the BLYP and PW91 exchange correlation functionals give qualitatively similar results. The time dependence of the hydration numbers are examined, and Wannier function techniques are used to analyze electronic configurations along the molecular dynamics trajectory.     

Ab initio molecular dynamics of liquid 1,3-dimethylimidazolium chloride

Density-functional-based Car-Parrinello molecular dynamics (CPMD) simulations have been performed for the ionic liquid 1,3-dimethylimidazolium chloride, [dmim]Cl, at 438 K. The local structure of the liquid is described in terms of various partial radial distribution functions and anisotropic spatial distributions, which reveal a significant extent of hydrogen bonding. The cation-anion distribution simulated with the BP86 functional is in qualitative agreement with the structural model derived from neutron diffraction data for the liquid, whereas the theoretical cation-cation distribution shows less satisfactory accord. Population analyses indicate noticeable charge transfer from anions to cations, and specific CH...Cl hydrogen bonds are characterized in terms of donor-acceptor interactions between lone pairs on Cl and antibonding sigma(CH) orbitals.     

Ab initio molecular dynamics of heme in cytochrome c

Ab initio molecular dynamics (AIMD) calculations, based on the Car-Parrinello method, have been carried out for three models of heme c that is present in cytochrome c. Both the reduced (Fe(II)) and oxidized (Fe(III)) forms have been analyzed. The simplest models (1R and 1O, respectively) consist of a unsubstituted porphyrin (with no side chains) and two axially coordinated imidazole and ethylmethylthioether ligands. Density functional theory optimizations of these models confirm the basic electronic features and are the starting point for building more complex derivatives. AIMD simulations were performed after reaching the thermal stability at T = 300 K. The evolution of the Fe-L(ax) bond strengths is examined together with the relative rotations of the imidazole and methionine about the axial vector, which appear rather independent from each other. The next models (2R and 2O) contain side chains at the heme to better simulate the actual active site. It is observed that two adjacent propionate groups induce some important effects. The axial Fe-Sdelta bond is only weakened in 2R but is definitely cleaved in the oxidized species 2O. Also the mobility of the Im ligand seems to be reduced by the formation of a strong hydrogen bond that involves the Im Ndelta1-Hdelta1 bond and one carboxylate group. In 2O the interaction becomes so strong that a proton transfer occurs and the propionic acid is formed. Finally, the models 3 include a free N-methyl-acetamide molecule to mimic a portion of the protein backbone. This influences the orientation of carboxylate groups and limits the amount of their hydrogen bonding with the Im ligand. Residual electrostatic interactions are maintained, which are still able to modulate the dissociation of the methionine from the heme.     

Ab initio molecular dynamics with dual basis set methods

On-the-fly, ab initio classical molecular dynamics are demonstrated with an underlying dual basis set potential energy surface. Dual-basis self-consistent field (Hartree-Fock and density functional theory) and resolution-of-the-identity second-order M?ller-Plesset perturbation theory (RI-MP2) dynamics are tested for small systems, including the water dimer. The resulting dynamics are shown to be faithful representations of their single-basis analogues for individual trajectories, as well as vibrational spectra. Computational cost savings of 58% are demonstrated for SCF methods, even relative to Fock-extrapolated dynamics, and savings are further increased to 71% with RI-MP2. Notably, these timings outperform an idealized estimate of extended-Lagrangian molecular dynamics. The method is subsequently demonstrated on the vibrational absorption spectrum of two NO(+)(H?O)? isomers and is shown to recover the significant width of the shared-proton bands observed experimentally.     

Ab initio Calculations on cytosine,thymine and adenine

An ab initio LCAO SCF calculation has been performed in a small contracted GTO basis set on three bases of ADN. A check of the representativity of the basis set is reported for formamide and pyrrole.

---

Zusammenfassung Es wurde eine ab initio LCAO-SCF-Rechnung mit einem kleinen, kontrahierten GTO-Basissatz für drei ADN-Basen durchgeführt. Eine Prüfung des Einflusses dieses Basissatzes wird für Formamid und Pyrrol beschrieben.

---

Résumé Un calcul ab initio dans une petite base de gaussiennes contractée a été fait pour trois bases de l'ADN. Un calcul-test sur la representativité de la base atomique est donné sur les exemples de la formamide et du pyrrole.

---

---

This work was supported by grant 67-00-532 of the Délégation Générale à la Recherche Scientifique et Technique (Comité de Biologie Moléculaire).     

Ab initio Ehrenfest dynamics

We present an ab initio direct Ehrenfest dynamics scheme using a three time-step integrator. The three different time steps are implemented with nuclear velocity Verlet, nuclear-position-coupled midpoint Fock integrator, and time-dependent Hartree-Fock with a modified midpoint and unitary transformation algorithm. The computational cost of the ab initio direct Ehrenfest dynamics presented here is found to be only a factor of 2-4 larger than that of Born-Oppenheimer (BO) dynamics. As an example, we compute the vibration of the NaCl molecule and the intramolecular torsional motion of H2C=NH2+ by Ehrenfest dynamics compared with BO dynamics. For the vibration of NaCl with an initial kinetic energy of 1.16 eV, Ehrenfest dynamics converges to BO dynamics with the same vibrational frequency. The intramolecular rotation of H2C=NH2+ produces significant electronic excitation in the Ehrenfest trajectory. The amount of nonadiabaticity, suggested by the amplitude of the coherent progression of the excited and ground electronic states, is observed to be directly related to the strength of the electron-nuclear coupling. Such nonadiabaticity is seen to have a significant effect on the dynamics compared with the adiabatic approximation.     

Ab initio calculations of triplet excited states and potential-energy surfaces of vinyl chloride: insights into spectroscopy and photodissociation dynamics

The equilibrium geometries and harmonic vibrational frequencies of three low-lying triplet excited states of vinyl chloride have been calculated using the state-averaged complete active space self-consistent field (CASSCF) method with the 6-311++G(d,p) basis set and an active space of four electrons distributed in 13 orbitals. Both adiabatic and vertical excitation energies have been obtained using the state-averaged CASSCF and the multireference configuration-interaction methods. The potential-energy surfaces of six low-lying singlet states have also been calculated. While the 3(pi, pi*) state has a nonplanar equilibrium structure, the 3(pi, 3s) and 3(pi, sigma*) states are planar. The calculated vertical excitation energy of the 3(pi, pi*) state is in agreement with the experiment. The singlet excited states are found to be multiconfigurational, in particular, the first excited state is of (pi, 3s) character at the planar equilibrium structure, of (pi, sigma*) as the C-Cl bond elongates, and of (pi, pi*) for highly twisted geometries. Avoided crossings are observed between the potential-energy surfaces of the first three singlet excited states. The absorption spectra of vinyl chloride at 5.5-6.5 eV can be unambiguously assigned to the transitions from the ground state to the first singlet excited state. The dissociation of Cl atoms following 193-nm excitation is concluded to take place via two pathways: one is through (pi, sigma*) at planar or nearly planar structures leading to fast Cl atoms and the other through (pi, pi*) at twisted geometries from which internal conversion to the ground state and subsequent dissociation produces slow Cl atoms.