The structure of the hydrated electron. Part 1. Magnetic resonance of internally trapping water anions: a density functional theory study

Density functional theory is used to rationalize magnetic parameters of hydrated electron trapped in alkaline glasses as observed using electron paramagnetic resonance (EPR) and electron spin echo envelope modulation (ESEEM) spectroscopies. To this end, model water cluster anions (n=4-8 and n=20, 24) that localize the electron internally are examined. It is shown that hyperfine coupling tensors of H/D nuclei in the water molecules are defined mainly by the cavity size and the coordination number of the electron; the water molecules in the second solvation shell play a relatively minor role. An idealized model of the hydrated electron (that is usually attributed to L. Kevan) in which six hydroxyl groups arranged in an octahedral pattern point toward the common center is shown to provide the closest match to the experimental parameters, such as isotropic and anisotropic hyperfine coupling constants for the protons (estimated from ESEEM), the second moment of the EPR spectra, and the radius of gyration. The salient feature is the significant transfer (10-20%) of spin density into the frontal O 2p orbitals of water molecules. Spin bond polarization involving these oxygen orbitals accounts for small, negative hyperfine coupling constants for protons in hydroxyl groups that form the electron-trapping cavity. In Part 2, these results are generalized for more realistic geometries of core anions obtained using a dynamic one-electron mixed quantum/classical molecular dynamics model.     

Reaction dynamics following electron capture of chlorofluorocarbon adsorbed on water cluster: a direct density functional theory molecular dynamics study

The electron capture dynamics of halocarbon and its water complex have been investigated by means of the full dimensional direct density functional theory molecular dynamics method in order to shed light on the mechanism of electron capture of a halocarbon adsorbed on the ice surface. The CF(2)Cl(2) molecule and a cyclic water trimer (H(2)O)(3) were used as halocarbon and water cluster, respectively. The dynamics calculation of CF(2)Cl(2) showed that both C-Cl bonds are largely elongated after the electron capture, while one of the Cl atoms is dissociated from CF(2)Cl(2) (-) as a Cl(-) ion. Almost all total available energy was transferred into the internal modes of the parent CF(2)Cl radical on the product state, while the relative translational energy of Cl(-) was significantly low due to the elongation of two C-Cl bonds. In the case of a halocarbon-water cluster system, the geometry optimization of neutral complex CF(2)Cl(2)(H(2)O)(3) showed that one of the Cl atoms interacts with n orbital of water molecules of trimer and the other Cl atom existed as a dangling Cl atom. After the electron capture, only one C-Cl bond (dangling Cl atom) was rapidly elongated, whereas the other C-Cl bond is silent during the reaction. The dangling Cl atom was directly dissociated from CF(2)Cl(2) (-)(H(2)O)(3) as Cl(-). The fast Cl(-) ion was generated from CF(2)Cl(2) (-)(H(2)O)(3) on the water cluster. The mechanism of the electron capture of halocarbon on water ice was discussed on the basis of the theoretical results.     

Dynamic behavior of chemical reactivity indices in density functional theory: A Bohn-Oppenheimer quantum molecular dynamics study

Dynamic behaviors of chemical concepts in density functional theory such as frontier orbitals (HOMO/LUMO), chemical potential, hardness, and electrophilicity index have been investigated in this work in the context of Bohn-Oppenheimer quantum molecular dynamics in association with molecular conformation changes. Exemplary molecular systems like CH ₅ ⁺ , Cl⁻ (H₂O)₃₀ and Ca²⁺ (H₂O)1₅ are studied at 300 K in the gas phase, demonstrating that HOMO is more dynamic than LUMO, chemical potential and hardness often fluctuate concurrently. It is argued that DFT concepts and indices may serve as a good framework to understand molecular conformation changes as well as other dynamic phenomena.     

An open‐source framework for analyzing N‐electron dynamics. II. Hybrid density functional theory/configuration interaction methodology

In this contribution, we extend our framework for analyzing and visualizing correlated many‐electron dynamics to non‐variational, highly scalable electronic structure method. Specifically, an explicitly time‐dependent electronic wave packet is written as a linear combination of N‐electron wave functions at the configuration interaction singles (CIS) level, which are obtained from a reference time‐dependent density functional theory (TDDFT) calculation. The procedure is implemented in the open‐source Python program det CI@ORBKIT, which extends the capabilities of our recently published post‐processing toolbox (Hermann et al., J. Comput. Chem. 2016, 37, 1511). From the output of standard quantum chemistry packages using atom‐centered Gaussian‐type basis functions, the framework exploits the multideterminental structure of the hybrid TDDFT/CIS wave packet to compute fundamental one‐electron quantities such as difference electronic densities, transient electronic flux densities, and transition dipole moments. The hybrid scheme is benchmarked against wave function data for the laser‐driven state selective excitation in LiH. It is shown that all features of the electron dynamics are in good quantitative agreement with the higher‐level method provided a judicious choice of functional is made. Broadband excitation of a medium‐sized organic chromophore further demonstrates the scalability of the method. In addition, the time‐dependent flux densities unravel the mechanistic details of the simulated charge migration process at a glance. © 2017 Wiley Periodicals, Inc.     

Exploring the role of decoherence in condensed-phase nonadiabatic dynamics: a comparison of different mixed quantum/classical simulation algorithms for the excited hydrated electron

Mixed quantum/classical (MQC) molecular dynamics simulation has become the method of choice for simulating the dynamics of quantum mechanical objects that interact with condensed-phase systems. There are many MQC algorithms available, however, and in cases where nonadiabatic coupling is important, different algorithms may lead to different results. Thus, it has been difficult to reach definitive conclusions about relaxation dynamics using nonadiabatic MQC methods because one is never certain whether any given algorithm includes enough of the necessary physics. In this paper, we explore the physics underlying different nonadiabatic MQC algorithms by comparing and contrasting the excited-state relaxation dynamics of the prototypical condensed-phase MQC system, the hydrated electron, calculated using different algorithms, including: fewest-switches surface hopping, stationary-phase surface hopping, and mean-field dynamics with surface hopping. We also describe in detail how a new nonadiabatic algorithm, mean-field dynamics with stochastic decoherence (MF-SD), is to be implemented for condensed-phase problems, and we apply MF-SD to the excited-state relaxation of the hydrated electron. Our discussion emphasizes the different ways quantum decoherence is treated in each algorithm and the resulting implications for hydrated-electron relaxation dynamics. We find that for three MQC methods that use Tully's fewest-switches criterion to determine surface hopping probabilities, the excited-state lifetime of the electron is the same. Moreover, the nonequilibrium solvent response function of the excited hydrated electron is the same with all of the nonadiabatic MQC algorithms discussed here, so that all of the algorithms would produce similar agreement with experiment. Despite the identical solvent response predicted by each MQC algorithm, we find that MF-SD allows much more mixing of multiple basis states into the quantum wave function than do other methods. This leads to an excited-state lifetime that is longer with MF-SD than with any method that incorporates nonadiabatic effects with the fewest-switches surface hopping criterion.     

Temperature dependence of structure and dynamics of the hydrated Ca2+ ion according to ab initio quantum mechanical charge field and classical molecular dynamics

Simulations using ab initio quantum mechanical charge field molecular dynamics (QMCF MD) and classical molecular dynamics using two‐body and three‐body potentials were performed to investigate the hydration of the Ca²⁺ ion at different temperatures. Results from the simulations demonstrate significant effects of temperature on solution dynamics and the corresponding composition and structure of hydrated Ca²⁺. Substantial increase in ligand exchange events was observed in going from 273.15 K to 368.15 K, resulting in a redistribution of coordination numbers to lower values. The effect of temperature is also visible in a red‐shift of the ion‐oxygen stretching frequencies, reflecting weakened ligand binding. Even the moderate increase from ambient to body temperature leads to significant changes in the properties of Ca²⁺ in aqueous environment. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

A theoretical study on the cyclopropane adsorption onto the copper surfaces by density functional theory and quantum chemical molecular dynamics methods

We report a theoretical study on the cyclopropane adsorption onto Cu(1 1 1) surfaces by density functional theory (DFT) and quantum chemical molecular dynamics methods. The equilibrium geometry of the physisorbed species was obtained using both periodic and cluster models by DFT methods that employ Cambridge serial total energy package (CASTEP), DMol ab initio quantum chemistry software of Accelrys’ materials studio (DMol), and Amsterdam density functional (ADF) program. It was found that the adsorbate molecule was tilted towards the metal surface with one C---C bond (upwards) parallel to the surface and that the physisorption occurred via a third carbon atom pointing (downwards) towards the surface. The electronic distribution and geometrical structure of physisorbed cyclopropane were slightly deviated from its gas phase molecule. The calculated vibrational frequencies and adsorption energies are close to experimental data, confirming the reliability of our DFT results. The adsorption process was simulated using our novel tight-binding quantum chemical molecular dynamics program, ‘Colors’. The calculation results indicated that both the adsorption and desorption processes of cyclopropane took place molecularly. The electron transfer and structural properties of equilibrium position obtained by ‘Colors’ are consistent with those by the first principles DFT methods.     

Optimal laser control of ultrafast photodissociation of I2- in water: mixed quantum/classical molecular dynamics simulation

A linearized optimal control method in combination with mixed quantum/classical molecular dynamics simulation is used for numerically investigating the possibility of controlling photodissociation wave packets of I(2)(-) in water. Optimal pulses are designed using an ensemble of photodissociation samples, aiming at the creation of localized dissociation wave packets. Numerical results clearly show the effectiveness of the control although the control achievement is reduced with an increase in the internuclear distance associated with a target region. We introduce effective optimal pulses that are designed using a statistically averaged effective dissociation potential, and show that they semiquantitatively reproduce the control achievements calculated by using optimal pulses. The control mechanisms are interpreted from the time- and frequency-resolved spectra of the effective optimal pulses.     

The onset of calcium carbonate nucleation: a density functional theory molecular dynamics and hybrid microsolvation/continuum study

Density functional theory (Perdew-Burke-Ernzerhof) based methods have been used to study the structure and hydration environment of the building blocks of CaCO 3 in aqueous solutions of calcium bicarbonate and calcium carbonate. Car-Parrinello molecular dynamics simulations of Ca(2+)/CO3(2-) and Ca (2+)/HCO3(-) in explicit water were performed to investigate the formation of CaCO3 and the hydration shell of the solvated hetero-ion pair. Our simulations show that the formation of the monomer of CaCO3 occurs with an associative mechanism and that the dominant building block of calcium (bi)carbonate in aqueous solution is Ca[eta(1)-(H)CO3](H2O)5, i.e., the preferred hydration number is five, while the (bi)carbonate is coordinated to the calcium in a monodentate mode. This result agrees with static calculations, where a hybrid approach using a combination of explicit solvent molecules and a polarizable continuum model has been applied to compute the solvation free energies of calcium bicarbonate species. Furthermore, the discrete-continuum calculations predict that the Ca(HCO3)2 and Ca(HCO3)3(-) species are stable in an aqueous environment preferentially as Ca(HCO3)2(H2O)4 and Ca(HCO3)3(H2O)2(-), respectively.     

A density functional theory based study of the microscopic structure and dynamics of aqueous HCl solutions

The aqueous solvation of hydrochloric acid is studied using density functional theory based molecular dynamics simulations at two concentrations. The large simulation boxes that we use allow us to investigate larger-scale structures such as the water-bridged chloride ion network. We find a strong concentration dependence for almost all structural and dynamical properties. Excess protons are mostly present both as Eigen and Zundel structures, either as a direct hydronium-chloride contact-ion pair or a solvent-separated ion pair. Increasing the concentration has a detrimental effect on the natural hydrogen bonded network of water molecules. This effect is visible in our studies as a decrease in the persistence time of the solvation shells around the chloride ions. Also the number of proton hops, determined by a new and well defined identification procedure, suffers from the breakdown of the natural hydrogen bond network.     

The cluster Ir4 and its interaction with a hydrogen impurity. A density functional study

To contribute to the understanding of how iridium particles act as catalysts for hydrogenation and dehydrogenation of hydrocarbons, we have determined structures and binding energies of various isomers of Ir(4) as well as HIr(4) on the basis of relativistic density functional theory. The most stable isomer of Ir(4) showed a square planar structure with eight unpaired electrons. The tetrahedral structure, experimentally suggested for supported species, was calculated 49 kJ mol(-1) less stable. Hydrogen coordinates preferentially to a single Ir center of the planar cluster with a binding energy of up to 88 kJ mol(-1) with respect to the atom in the H(2) molecule. Terminal interaction of hydrogen with an Ir(4) tetrahedron causes the cluster to open to a butterfly structure. We calculated terminal binding of hydrogen at different Ir(4) isomers to be more stable than bridge coordination, at variance with earlier studies.     

Hydrogen-bonding interaction of urea with DNA bases: A density functional theory study

This work deals with the interaction between urea and DNA bases (adenine, thymine, guanine, and cytosine). The optimized geometries, binding energies, and harmonic vibrational frequencies are calculated using the DFT/B3LYP functional combined with the 6–31+G(d,p) basis set. Their interactions are studied aiming to understand more about the nature of the intercalation binding forces between urea and DNA. Fourteen stable complexes are found on the potential energy surface. The structures are cyclic; they are stabilized by NH...O/N and CH...O interactions. The binding energies range from −19.9 kJ·mol⁻¹ to −74.0 kJ·mol⁻¹. The obtained formation energies indicate that Urea:G and Urea:C are more favorable than Urea:T and Urea:A. In addition, the Atoms in Molecules theory is performed to study the hydrogen bonds in the complexes.     

The structure of 1-chlorosilatrane: An ab initio molecular orbital and a density functional theory study

The molecular geometries of the 1-chloro-, 1-fluoro-, 1-methyl-, and 1-hydrogenosilatranes were fully optimized by the restricted Hartree-Fock (HF) method supplemented with 3-21G, 3-21G(d), 6-31G(d), and CEP-31G(d) basis sets; by MP2 calculations using 6-31G(d) and CEP-31G(d) basis sets; and by GGA-DFT calculations using 6-31G(d5) basis set with the aim of locating the positions of the local minima on the energy hypersurface. The HF/6-31G(d) calculations predict long (>254 pm) and the MP2/CEP calculations predicted short (∼225 pm) equilibrium Si(SINGLE BOND)N distances. The present GGA-DFT calculations reproduce the available gas phase experimental Si(SINGLE BOND)N distances correctly. The solid phase experimental results predict that the Si(SINGLE BOND)N distance is shorter in 1-chlorosilatrane than in 1-fluorosilatrane. In this respect the HF results show a strong basis set dependence, the MP2/CEP results contradict the experiment, and the GGA-DFT results in electrolytic medium agree with the experiment. The latter calculations predict that 1-chlorosilatrane is more polarizable than 1-fluorosilatrane and also support a general Si(SINGLE BOND)N distance shortening trend for silatranes during the transition from gas phase to polar liquid or solid phase. The calculations predict that the ethoxy links of the silatrane skeleton are flexible. Consequently, it is difficult to measure experimentally the related bond lengths and bond and torsion angles. This is the probable origin of the surprisingly large differences for the experimental structural parameters. On the basis of experimental analogies, ab initio calculations, and density functional theory (DFT) calculations, a gas phase equilibrium (r_e) geometry is predicted for 1-chlorosilatrane. The semiempirical methods predict a so-called exo minimum (at above 310 pm Si(SINGLE BOND)N distance); however, the ab initio and GGA-DFT calculations suggest that this form is nonexistent. The GGA-DFT geometry optima were characterized by frequency analysis. © 1996 by John Wiley & Sons, Inc.     

Interaction of acetonitrile with alkaline metal cations: A density functional,coupled‐cluster,and quadratic configuration interaction study

A systematic quantum chemical study of CH₃CN and its CH₃CN?M⁺ 1:1 model adducts (M⁺∈{Li⁺, Na⁺}) is presented, with respect to binding energetics, structural and vibrational force field changes. Several gradient‐corrected density functional levels of theory were employed (of both “pure” and “hybrid” character), together with the coupled cluster including double substitutions from the Hartree–Fock determinant (CCD) and quadratic configuration interaction including single and double substitutions (QCISD) methods [with the rather large 6‐311G(d,p) basis set], and their performances compared. The binding energy decompositions according to the Kitaura–Morokuma scheme and the reduced variational space self‐consistent field (RVS‐SCF) method have shown that the electrostatic plus polarization interaction terms are primarily responsible for overall stabilization, while the charge‐transfer term is negligibly small and virtually identical for both adducts. The computed harmonic vibrational frequencies for acetonitrile correlate excellently with the experimental ones (r²>0.9998 for almost all cases, while for the BLYP level, r²=1). It is shown for the first time that the experimentally observed blue shifts of the ν_CN mode are caused even by formation of 1:1 adducts, contrary to the previously accepted opinions. The CCD and QCISD, as well as the BPW91 and BP86 levels of theory predict almost excellently the ν_CN mode blue shift upon adduct formation, while the BLYP and B3LYP levels perform significantly poorer. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Self-interaction corrected density functional theory for the study of intramolecular electron transfer dynamics in radical carbocations

We derived an orbital dependent Kohn-Sham based scheme for the correction of the self-interaction error in DFT, which is particularly suited for the study of open shell molecular systems. Our approach is based on a weighted form of the Perdew and Zunger (PZ) self-interaction correction scheme, in which an empirical coefficient is introduced in front of the SIC term to remove overcorrections. The method is used to investigate the first principle molecular dynamics of the intramolecular electron transfer (IET) in bis(methylene) adamanthyl radical cation and allows the analysis of the free energy surface that governs the IET process at room temperature (300 K). Indeed, the thermal activation of all molecular degrees of freedom provides important additional information about the mechanisms involved in the IET process. Our study confirms and extends previous results obtained with CASSCF and shows that there is no predominant degeneracy-lifting mode even at room temperature. However, we also identified regions in the phase space for which there is a significant probability for a productive IET event. In addition, we performed thermodynamic integrations along selected reaction coordinates to determine an estimate of the activation free energy barrier for the IET process.     

Time-dependent density functional theory excited state nonadiabatic dynamics combined with quantum mechanical/molecular mechanical approach: photodynamics of indole in water

We present a combination of time-dependent density functional theory with the quantum mechanical/molecular mechanical approach which can be applied to study nonadiabatic dynamical processes in molecular systems interacting with the environment. Our method is illustrated on the example of ultrafast excited state dynamics of indole in water. We compare the mechanisms of nonradiative relaxation and the electronic state lifetimes for isolated indole, indole in a sphere of classical water, and indole + 3H(2)O embedded in a classical water sphere. In the case of isolated indole, the initial excitation to the S(2) electronic state is followed by an ultrafast internal conversion to the S(1) state with a time constant of 17 fs. The S(1) state is long living (>30 ps) and deactivates to the ground state along the N-H stretching coordinate. This deactivation mechanism remains unchanged for indole in a classical water sphere. However, the lifetimes of the S(2) and S(1) electronic states are extended. The inclusion of three explicit water molecules opens a new relaxation channel which involves the electron transfer to the solvent, leading eventually to the formation of a solvated electron. The relaxation to the ground state takes place on a time scale of 60 fs and contributes to the lowering of the fluorescence quantum yield. Our simulations demonstrate the importance of including explicit water molecules in the theoretical treatment of solvated systems.     

Binding of pertechnetate to uranyl(VI) in aqueous solution. A density functional theory molecular dynamics study

According to constrained Car-Parrinello molecular dynamics simulations and thermodynamic integration, the free binding energy between uranyl hydrate and pertechnetate in aqueous solution is significantly lower than that between uranyl and nitrate, namely, by 1.7 kcal mol(-1). This is the first study of the differential binding of these two ligands to uranyl, which can have implications for the separability of uranium and technetium during the reprocessing of nuclear waste.     

A study of electron affinities using the initiator approach to full configuration interaction quantum Monte Carlo

For the atoms with Z ≤ 11, energies obtained using the "initiator" extension to full configuration interaction quantum Monte Carlo (i-FCIQMC) come to within statistical errors of the FCIQMC results. As these FCIQMC values have been shown to converge onto FCI results, the i-FCIQMC method allows similar accuracy to be achieved while significantly reducing the scaling with the size of the Slater determinant space. The i-FCIQMC electron affinities of the Z ≤ 11 atoms in the aug-cc-pVXZ basis sets are presented here. In every case, values are obtained to well within chemical accuracy [the mean absolute deviation (MAD) from the relativistically corrected experimental values is 0.41 mE(h)], and significantly improve on coupled cluster with singles, doubles and perturbative triples [CCSD(T)] results. Since the only remaining source of error is basis set incompleteness, we have investigated using CCSD(T)-F12 contributions to correct the i-FCIQMC results. By doing so, much faster convergence with respect to basis set size may be achieved for both the electron affinities and the FCIQMC ionization potentials presented in a previous paper. With this F12 correction, the MAD can be further reduced to 0.13 mE(h) for the electron affinities and 0.31 mE(h) for the ionization potentials.     

Stereoselective aziridination of imines via ammonium ylides: A molecular electron density theory study

The highly trans-stereoselective reaction of ammonium salt ( AS 2 ) toward (E)-N-benzylidene-4-methylbenzenesulfonamide ( IM 4 ) in the presence of Na₂CO₃ leading to an aziridine derivative, trans -Az 6 , was theoretically studied using Molecular Electron Density Theory (MEDT) at the B3LYP/6-31G(d) computational level to probe energy transformation, selectivities, and molecular mechanism. The reaction starts by a nucleophilic substitution reaction between 1,4-diazabicyclo [2.2.2]octane ( DABCO ) and phenacyl bromide ( PB ) to form AS 2 which undergoes a proton abstraction by Na₂CO₃ to produce Ammonium Ylide ( AY 3) . Subsequently, nucleophilic addition of AY 3 to the double bond in IM 4 leads to form a betaine-like intermediate, namely, IN-Ta which named IN-Ta . Finally, trans -Az 6 is produced as a result of the nucleophilic attack of the negatively charged nitrogen atom on the carbon atom bearing DABCO in IN-Ta . Analysis of the relative Gibbs free energies shows that the ring closure step is the rate-determining step (RDS). By an investigation of the conceptual density functional theory, CDFT reactivity indices AY 3 and IM 4 are classified as a strong nucleophilic molecular system and as a strong electrophilic molecular system, respectively, which implies that the addition step of AY 3 to IM 4 has a high polar character. An analysis of the calculated electrophilic and nucleophilic Parr functions at the reactive sites of reagents clarifies the regioselectivity observed experimentally within the C₁–C₂ bond creation process.