Ab initio direct dynamics study of cyclopropyl radical ring-opening

Quasiclassical direct dynamics simulations, at the CASSCF(3,3)/6-31G(d) level of theory, are used to study the stereochemistry of the electrocyclic ring-opening reaction of the cyclopropyl radical. The trajectories are initiated at the reaction's transition state (TS), with their initial conditions sampled from the TS's 174 degrees C Boltzmann distribution. Intrinsic reaction coordinate calculations predict the overall reaction to have disrotatory stereochemistry. Though this is the preferred initial reaction stereochemistry in the trajectories, 43% of the trajectories follow the conrotatory path. Four unique trajectory types are observed during 200 fs dynamics of the product allyl radical. Intramolecular vibrational energy redistribution and internal rotation are incomplete on this time scale, and a statistical distribution of the allyl isomers is not observed.     

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 and experimental studies of the large-amplitude motion in BF2OH

The ground state rotational spectrum of BF2OH was measured under high resolution by microwave Fourier transform spectroscopy (FTMW), and the small torsional splitting could be resolved for several lines. This splitting was analyzed using a phenomenological model previously developed for HNO3 [Coudert and Perrin, J. Mol. Spectrosc. 1995, 172, 352] and with the help of the geometries of the stationary points calculated ab initio. The torsional splitting was also calculated using the results of the calculations for the ground vibrational state, for the excited OH torsional states 91 and 92, and for the excited BOH bending state 41, and a satisfactory agreement with available experimental data was found.     

Ab initio quantum dynamics using coupled-cluster

The curse of dimensionality (COD) limits the current state-of-the-art ab initio propagation methods for non-relativistic quantum mechanics to relatively few particles. For stationary structure calculations, the coupled-cluster (CC) method overcomes the COD in the sense that the method scales polynomially with the number of particles while still being size-consistent and extensive. We generalize the CC method to the time domain while allowing the single-particle functions to vary in an adaptive fashion as well, thereby creating a highly flexible, polynomially scaling approximation to the time-dependent Schro?dinger equation. The method inherits size-consistency and extensivity from the CC method. The method is dubbed orbital-adaptive time-dependent coupled-cluster, and is a hierarchy of approximations to the now standard multi-configurational time-dependent Hartree method for fermions. A numerical experiment is also given.     

Ab initio photoionization dynamics of beta-alanine

Photoionization dynamics of beta-alanine is studied by the trajectory simulations using the ab initio potential energy surface. Vertical photoionization in the spirit of the Franck-Condon principle is assumed both for the adiabatic and thermostatic simulations. Both intramolecular proton transfer and fragmentation while only the proton transfer are found in the thermostatic and adiabatic simulations, respectively, for the conformer having the intramolecular hydrogen bond N...H-O. The theoretical predictions are in line with the experimental observations available in the literatures. It is reported for the first time that the thermostatic temperatures strongly affect the fragmentation processes induced by photoionization.     

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.     

Ab initio study of the cyclooctatetraenyl radical

Ab initio calculations have been carried out on the 1,3,5,7- and 1,2,4,7-tetraene configurations of the cyclooctatetraenyl radical at UHF, ROHF, MCSCF, ROCISD, QCISD, and CCSD(T) levels of theory with 6-311G(d,p) and cc-pVDZ basis sets. Although spin contamination is present, the ROCISD calculations support the energies obtained from less intensive, UHF-based coupled cluster calculations over the energies obtained from MCSCF analysis of the pi-electron orbitals. The 1,3,5,7-form is a local minimum at the coupled cluster levels, higher in energy than the resonance-stabilized 1,2,4,7-form by 10-13 kJ/mol, but bounded by a barrier of less than 0.5 kJ/mol. The isomerization surface connecting these two structures is described and results reported from integration of the vibrational Schr?dinger equation on that surface. Excited vibrational states at energies just above the isomerization barrier are dominated by the character of the 1,3,5,7-tetraenyl radical, which suggests that chemistry involving this intermediate at typical combustion temperatures may branch at this juncture.     

Ab initio dynamics study of poly-para-phenylene vinylene

We present an ab initio dynamics investigation within a density-functional perturbation theory framework of the properties of the conjugated polymer poly-para-phenylene vinylene (PPV) in both the isolated chain and crystalline states. The calculated results show that for an isolated chain, most of the vibrational modes correspond to experimentally observed modes in crystalline PPV. However, additional hitherto unidentified modes have been observed in experiment and our calculations on crystalline material have allowed us to assign these. We also present the results of calculations of the polarizability and permittivity tensors of the material, which are in reasonable agreement with the typical values for conjugated polymers. Dynamical Born effective charges [S. Baroni, S. de Gironcoli, A. Dal Corso, and P. Giannozzi, Rev. Mod. Phys. 73, 515 (2001)] are calculated and compared with atomic charges obtained from Mulliken population analysis [M. D. Segall, C. J. Pickard, R. Shah, and M. C. Payne, Mol. Phys. 89, 571 (1996)] and we conclude that effective charges are more appropriate for use in the study of the dynamics of the system. Notable differences are found in the infrared-absorption spectra obtained for the isolated chain and crystalline states, which can be attributed to the differences in the crystalline packing effects, which clearly play a key role in influencing the lattice dynamics of PPV.     

Ab initio simulation of the vibrationally resolved photoelectron spectrum of Si3-

Electron photodetachment spectra provide a wealth of information about the electronic and vibrational level structures of neutral molecules that form stable anions. Experiments carried out for the smallest polyatomic silicon cluster anion (Si3-+hupsilon-->Si3*+e-) show vibrational progressions in six observed electronic bands (X-E) of the neutral species. The authors have performed ab initio calculations using the MRCI+D/aug-cc-pVQZ level for the corresponding electronic states followed by variational calculations of the vibronic levels associated with these adiabatic potential energy surfaces. In contrast to previous approaches, the authors treat the nonadiabatic dynamics on the potential energy surfaces, which allows for a vastly improved reproduction of the experimental level structure and a corrected assignment for band A.     

Ab initio study of hydrazinyl radical: toward a DMBE potential energy surface

A series of stationary structures of the hydrazinyl radical have been characterized by optimization at the CCSD(T)/cc-pVTZ level of theory. CCSD(T)/aug-cc-pVXZ single-point calculations have also been carried out at the optimized geometries with basis sets of different cardinal numbers (X = T, Q), which were used to obtain accurate energies via extrapolation to the complete basis set limit. A discussion on the analytical modeling of the potential energy surface of hydrazinyl is also presented.     

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 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 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 dynamics with wave-packets and density matrices

Efficient methodologies to conduct simultaneous dynamics of electrons and nuclei are discussed. Particularly, attention is directed to a recent development that combines quantum dynamics with ab initio molecular dynamics. The two components of the methodology, namely, quantum dynamics and ab initio molecular dynamics, are harnessed together using a time-dependent self-consistent field-like coupling procedure. An approach to conduct quantum dynamics using an accurate banded, sparse and Toeplitz representation for the discrete free propagator is highlighted with suitable review of other related approaches. One notable feature of the method is that all important quantum dynamical effects including zero-point effects, tunneling as well as over-barrier reflections are accurately treated. Computational methodologies for improved efficiency of the quantum dynamics are also discussed. There exists a number of ways to carry out simultaneous ab initio molecular dynamics (such as Born–Oppenheimer dynamics and extended Lagrangian dynamics, Car–Parrinello dynamics being a prime example of the latter); our prime focus remains on atom-centered density-matrix propagation and Born–Oppenheimer dynamics. The electronic degrees of freedom are handled at accurate levels of density functional theory, using hybrid or gradient corrected approximations. Benchmark calculations are provided for a prototypical proton transfer system. Future generalizations and goals are discussed.     

Ab initio inter-radical potential for p-phenylenediamine radical cation dimer

Ab initio molecular orbital calculations were carried out to investigate the inter-radical interaction of the paired p-phenylenediamine radical cations in the singlet state. After initial optimization of the dimer in the parallel sandwich (D2h) and parallel displaced (Cs) configurations at the B3LYP/6-31G* theoretical level, the MP2/6-31G* and B3LYP/6-31G* single energies of the dimer were calculated as a function of the inter-radical distance R. The depths of the potential minima near R = 3.2 A were estimated to be in the order of the hydrogen bonding energy, assuming that the electrostatic contribution between the cations is canceled out by the attractive contributions due to the counter anions on the aspect of a simple electrostatic model. This can be related to the indications of the cation dimer formation in solution in the presence of counter anions at a low temperature reported previously in the literature by resonance Raman and electronic absorption spectra. The inter-radical (Raman active) frequencies of the dimer were calculated, one of which corresponds to the reported value at 161 cm(-1) observed in the resonance Raman spectrum in ethanol at 200 K by Yokoyama and Maeda (Chem. Phys. Lett. 48 (1977) 59).     

Ab initio study of water clustering in the presence of a methyl radical

The geometrical structure and binding energy of small clusters of methyl radical and water molecules (up to five water molecules) in gas phase and water media have been investigated at the MP2 level of theory using 6-311++G(2df,2p) basis set. The complexes characterized contain OH···O, CH···O, and OH···C attractive interactions with stabilization energies in the range 6–143 kJ mol^?1. The solvent has an enhancing influence on the stabilities of studied clusters. The atoms in molecules theory were also applied to explain the nature of the complexes. The interaction energies have been partitioned with the natural energy decomposition analysis showing that the most important attractive term corresponds to the charge transfer one.     

Ab initio molecular dynamics simulation of a room temperature ionic liquid

Ab initio molecular dynamics simulations have been performed for the first time on the room-temperature organic ionic liquid dimethyl imidazolium chloride [DMIM][Cl] using density functional theory. The aim is to compare the local liquid structure with both that obtained from two different classical force fields and from neutron scattering experiments. The local structure around the cation shows significant differences compared to both the classical calculations and the neutron results. In particular, and unlike in the gas-phase ion pair, chloride ions tend to be located near a ring C-H proton in a position suggesting hydrogen bonding. The results are used to suggest ways in which the classical potentials may be improved.