DNA base trimers: empirical and quantum chemical ab initio calculations versus experiment in vacuo

A complete scan of the potential‐energy surfaces for selected DNA base trimers has been performed by a molecular dynamics/quenching technique using the force field of Cornell et al. implemented in the AMBER7 program. The resulting most stable/populated structures were then reoptimized at a correlated ab initio level by employing resolution of the identity, Møller–Plesset second‐order perturbation theory (RI‐MP2). A systematic study of these trimers at such a complete level of electronic structure theory is presented for the first time. We show that prior experimental and theoretical interpretations were incorrect in assuming that the most stable structures of the methylated trimers corresponded to planar systems characterized by cyclic intermolecular hydrogen bonding. We found that stacked structures of two bases with the third base in a T‐shape arrangement are the global minima in all of the methylated systems: they are more stable than the cyclic planar structures by about 10 kcal mol^?1. The different behaviors of nonmethylated and methylated trimers is also discussed. The high‐level geometries and interaction energies computed for the trimers serve also as a reference for the testing of recently developed density functional theory (DFT) functionals with respect to their ability to correctly describe the balance between the electrostatic and dispersion contributions that bind these trimers together. The recently reported M052X functional with a polarized triple‐zeta basis set predicts 11 uracil trimer interaction energies with a root‐mean‐square error of 2.3 kcal mol^?1 relative to highly correlated ab initio theoretical calculations.     

Pure trinuclear 4f single-molecule magnets: synthesis, structures, magnetism and ab initio investigation

A family of linear Dy₃ and Tb₃ clusters have been facilely synthesized from the reactions of DyCl₃, the polydentate 3‐methyloxysalicylaldoxime (MeOsaloxH₂) ligand with auxiliary monoanionic ligands, such as trichloroacetate, NO₃^?, OH^?, and Cl^?. Complexes 1 – 5 contain a nearly linear Ln₃ core, with similar Ln???Ln distances (3.6901(4)–3.7304(3) Å for the Dy₃ species, and 3.7273(3)–3.7485(5) Å for the Tb₃ species) and Ln???Ln???Ln angles of 157.036(8)–159.026(15)° for the Dy₃ species and 157.156(8)–160.926(15)° for the Tb₃ species. All three Ln centers are bridged by the two doubly‐deprotonated [MeOsalox]^2? ligands and two of the four [MeOsaloxH]^? ligands through the N,O‐η²‐oximato groups and the phenoxo oxygen atoms (Dy‐O‐Dy angles=102.28(16)–106.85(13)°; Tb‐O‐Tb angles=102.00(11)–106.62(11)°). The remaining two [MeOsaloxH]^? ligands each chelate an outer Ln^III center through their phenoxo oxygen and oxime nitrogen atoms. Magnetic studies reveal that both Dy₃ and Tb₃ clusters exhibit significant ferromagnetic interactions and that the Dy₃ species behave as single‐molecule magnets, expanding upon the recent reports of the pure 4f type SMMs.     

Complexes of ammonia with propane and cyclopropane: electrostatic guidelines for ab initio treatment

A model based on the molecular electrostatic potential (MESP) is employed for the investigation of structures and energies of complexes of ammonia with propane and cyclopropane. The electrostatic model geometries are employed as starting points for an ab initio investigation at the self-consistent field and second-order M?ller-Plesset (MP2) levels. The most stable structures of C₃H₆..NH₃ and C₃H₈..NH₃ complexes have the interaction energies of 10.07 kJ/mol and 8.15 kJ/mol, respectively, at the MP2/6-31G(d,p) level. The energy rank order of the structures is not altered with the use of the 6-31++G(d,p) basis set, and the basis␣set superposition error has little effect. The interaction energy decomposition analysis shows that the electrostatic component is dominant over the other ones. MESP topography thus seems to offer valuable hints for predicting the structures of weakly bonded complexes. Received: 8 July 1998 / Accepted: 4 August 1998 / Published online: 2 November 1998     

Electron delocalization in mixed-valence Keggin polyoxometalates. Ab initio calculation of the local effective transfer integrals and its consequences on the spin coupling

We present a quantitative evaluation of the influence of the electron transfer on the magnetic properties of mixed-valence polyoxometalates reduced by two electrons. For that purpose, we extract from valence-spectroscopy ab initio calculations on embedded fragments the value of the transfer integrals between W nearest-neighbor atoms in a mixed-valence alphaPW(12)O(40) polyoxowolframate Keggin anion. In contradiction with what is usually assumed, we show that the electron transfer between edge-sharing and corner-sharing WO(6) octahedra have very close values. Considering fragments of various ranges, we analyze the accuracy of calculations on fragments based on only two WO(5) pyramids which should allow a low cost general study of transfer parameters in polyoxometalates. Finally, these parameters are introduced in an extended Hubbard Hamiltonian that models the whole anion. It permits to prove that electron transfers induce a large energy gap between the singlet ground state and the lowest triplet states providing a clear explanation of the diamagnetic properties of the mixed-valence Keggin ions reduced by two electrons.     

Interaction energies for the purine inhibitor roscovitine with cyclin-dependent kinase 2: correlated ab initio quantum-chemical, DFT and empirical calculations

The interaction between roscovitine and cyclin-dependent kinase 2 (cdk2) was investigated by performing correlated ab initio quantum-chemical calculations. The whole protein was fragmented into smaller systems consisting of one or a few amino acids, and the interaction energies of these fragments with roscovitine were determined by using the MP2 method with the extended aug-cc-pVDZ basis set. For selected complexes, the complete basis set limit MP2 interaction energies, as well as the coupled-cluster corrections with inclusion of single, double and noninteractive triples contributions [CCSD(T)], were also evaluated. The energies of interaction between roscovitine and small fragments and between roscovitine and substantial sections of protein (722 atoms) were also computed by using density-functional tight-binding methods covering dispersion energy (DFTB-D) and the Cornell empirical potential. Total stabilisation energy originates predominantly from dispersion energy and methods that do not account for the dispersion energy cannot, therefore, be recommended for the study of protein-inhibitor interactions. The Cornell empirical potential describes reasonably well the interaction between roscovitine and protein; therefore, this method can be applied in future thermodynamic calculations. A limited number of amino acid residues contribute significantly to the binding of roscovitine and cdk2, whereas a rather large number of amino acids make a negligible contribution.     

DFT and ab initio calculations on two reactions between hydrogen atoms and the fire suppressants 2-H heptafluoropropane and CF3Br

Reaction enthalpies and barrier heights of the reactions CF3Br+H-->CF3+HBr {reaction (1)} and CF3CHFCF3+H-->CF3CFCF3+H2 {reaction (2)} have been calculated at the near state-of-the-art ab initio level, and also by employing the B3LYP, BH&HLYP, BB1K, MPW1K, MPWB1K and TPSS1KCIS functionals. In addition, the integrated molecular orbital+molecular orbital (IMOMO) method has been used to study reaction (2). The ab initio benchmark values of the reaction enthalpy (298 K) and barrier height (0 K) of reaction (2) are reported for the first time {-(0.7+/-0.7) and 13.3+/-0.5 kcal/mole respectively}. When density functional theory (DFT) results are compared with ab initio benchmarks for both reactions (1) and (2), the MPWB1K functional is found to have the best performance of the six functionals used. The IMOMO method with the RCCSD/aug-cc-pVTZ and/or RCCSD(T)/aug-cc-pVTZ levels, as the high levels of calculation on the model system, gives reaction enthalpies and barrier heights of reaction (2), which agree with ab initio benchmark values to within 1 kcal/mole. Computed key geometrical parameters and imaginary vibrational frequencies of the transition state structures of reactions (1) and (2) obtained at different levels of calculation are compared. The magnitudes of the computed imaginary vibrational frequencies of the transition states of both reactions considered are found to be very sensitive to the levels of calculation used to obtain them. The heat of formation (298 K) of CF3CFCF3 calculated at the near state-of-the-art level has a value of -(318+/-3) kcal/mole.     

The mechanism of formamide hydrolysis in water from ab initio calculations and simulations

The neutral hydrolysis of formamide in water is a suitable reference to quantify the efficiency of proteolytic enzymes. However, experimental data for this reaction has only very recently been obtained and the kinetic constant determined experimentally is significantly higher than that predicted by previous theoretical estimations. In this work, we have investigated in detail the possible mechanisms of this reaction. Several solvent models have been considered that represent a considerable improvement on those used in previous studies. Density functional and ab initio calculations have been carried out on a system which explicitly includes the first solvation shell of the formamide molecule. Its interaction with the bulk has been treated with the aid of a dielectric continuum model. Molecular dynamics simulations at the combined density functional/molecular mechanics level have been carried out in parallel to better understand the structure of the reaction intermediates in aqueous solution. Overall, the most favored mechanism predicted by our study involves two reaction steps. In the first step, the carbonyl group of the formamide molecule is hydrated to form a diol intermediate. The corresponding transition structure involves two water molecules. From this intermediate, a water-assisted proton transfer occurs from one of the hydroxy groups to the amino group. This reaction step may lead either to the formation of a new reaction intermediate with a marked zwitterionic character or to dissociation of the system into ammonia and formic acid. The zwitterionic intermediate dissociates quite easily but its lifetime is not negligible and it could play a role in the hydrolysis of substituted amides or peptides. The predicted pseudo-first-order kinetic constant for the rate-limiting step (the first step) of the hydrolysis reaction at 25 degrees C (3.9x10(-10) s(-1)) is in excellent agreement with experimental data (1.1x10(-10) s(-1)).     

FA(I):A(+) and FA(II):Cu(+) laser activity and photographic sensitization at the low coordinated surfaces of AgBr ab initio calculations

The twofold potentials of F(A)(I):Au(+) and F(A)(II)Cu(+) color centers at the low coordinated surfaces of AgBr thin films in providing tunable laser activity and photographic sensitization were investigated using ab initio methods of molecular electronic structure calculations. Clusters of variable size were embedded in simulated Coulomb fields that closely approximated the Madelung fields of the host surfaces, and the nearest neighbor ions to the F(A) defect site were allowed to relax to equilibrium in each case. Based on the calculated Stokes shifted optical transition bands and horizontal shifts along the configuration coordinate diagrams, both F(A)(I):Au(+) and F(A)(II):Cu(+) color centers were found to be laser active. The laser activity faded quickly as the bromide ion coordination decreased from 5 (flat) to 4 (edge) to 3 (corner) and as the size of the impurity cation increased from Cu(+) to Au(+). The latter relation was explainable in terms of the axial perturbation of the impurity cation. The smallest calculated Stokes-shift at the corner surface suggested that emission had the same oscillator strength as absorption. All relaxed excited states RESs of the defect containing surfaces were deep below the lower edges of the conduction bands of the defect free ground state surfaces, indicating that F(A)(I):Au(+) and F(A)(II):Cu(+) are suitable laser defects. The probability of orientational destruction of the two centers attributed to the assumed RES saddle point ion configurations along the <110> axis was found to be directly proportional to the size of the impurity cation, with activation energy barriers of about 0.655-3.294 eV for Cu(+), and about 1.887-3.404 eV for Au(+). The possibility of exciton (energy) transfer from the sites of higher coordination to those of lower coordination is demonstrated. The more laser active F(A)(II):Cu(+) center was more easily formed than the less laser active F(A)(I):Au(+) center. The Glasner-Tompkins empirical relation was generalized to include F(A) centers at the low coordinated surfaces of silver bromide thin film. As far as color photographic sensitization is concerned, the lowest unoccupied molecular orbitals of the selected dye molecules in the excited states were high enough for electron injection. F(A) defect formation and rotational diffusion of silver clusters reduced the energy gaps between the excited dye molecules and the lower edges of the conduction bands and allowed for hole injection. About 54-60% of the reduction of silver ions at the flat surface of AgBr was attributed to the host anions and F(A) defect formation, leaving about 40-46% for the reduction of photoelectrons as well as the electrons of the developer or dye molecules. The unrelaxed rotational diffusions of the central Ag(4) by 90 degrees decreased the latter percentage, but were severely hindered by activation energy barriers.     

Electron-spin magnetic moments of the 2Σ+ ions Li2+, Li2−, and Be2+: An ab initio ROHF study

For the ²Σ⁺ ground states of the ions Li₂⁺, Li₂⁻, and Be₂⁺, the dependence of the magnetic moment (parametrized by g-shifts) on the bond length R was studied at the ROHF level. The Δ g-values were calculated via a perturbative approach (complete to second order in Breit-Pauli interactions) using quadruple-zeta AO basis sets augmented by semidiffuse and polarization functions. All Δ g-values in these systems are negative. The parallel component Δ g_∥ generally changes little with R, remaining close to the g-shift of the corresponding ²S atomic dissociation product. For Li₂⁺ and Be₂⁺, the perpendicular component Δ g _⟂ is more sensitive to geometry than is Δ g_∥, mainly because of the second-order magnetic coupling with excited ²Π states. For Li₂⁻, Δ g _⟂ and Δ g_∥ are similar due to the large size of the 2σ_u, SOMO, resulting in g-values close to that of a free electron. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 63: 511–521, 1997     

Spectroscopic properties in the liquid phase: combining high-level ab initio calculations and classical molecular dynamics.

We present an integrated computational tool, rooted in density functional theory, the polarizable continuum model, and classical molecular dynamics employing spherical boundary conditions, to study the spectroscopic observables of molecules in solution. As a test case, a modified OPLS-AA force field has been developed and used to compute the UV and NMR spectra of acetone in aqueous solution. The results show that provided the classical force fields are carefully reparameterized and validated, the proposed approach is robust and effective, and can also be used by nonspecialists to provide a general and powerful complement to experimental techniques.     

Single‐ and double‐electron reductions of CO2 by using superalkalis: An ab initio study

CO₂, a major contributor to global warming, can be balanced by converting it into fuels. The reduction of CO₂ has been difficult due to its extremely high stability. Recently, single‐electron reduction of CO₂ by superalkalis has been proposed using quantum chemical methods. Herein, we report a systematic study on the single‐reduction of CO₂ by using typical superalkalis. Superalkalis are hypervalent species possessing lower ionization energies than alkali atoms. We have studied the interaction of CO₂ with FLi₂, OLi₃, and NLi₄ superalkalis using ab initio MP2 calculations. We notice that this interaction leads to stable superalkali‐CO₂ complexes in which the structure of CO₂ is bent due to electron transfer from superalkalis. This clearly reveals that the CO₂ can successfully be reduced to the anion. It has been also noticed that the size of superalkalis plays a crucial in the single‐electron reduction of CO₂. For instance, the binding energy of superalkali‐CO₂ complex and charge transfer to CO₂ decreases monotonically with the increase in the size of superalkali. We have also proposed that CO₂ can be further reduced to in case of the anionic complex such as (FLi₂ CO₂)‾. Thus, FLi₂ superalkali is also capable of double‐electron reduction of CO₂. These findings should provide new insights into CO₂‐activation as well as motivate further research in this direction.     

The singlet-triplet separation in CF2: state-of-the-art ab initio calculations and Franck-Condon simulations including anharmonicity.

Geometrical parameters, vibrational frequencies and relative electronic energies of the X2B1 state of CF2- and the X1A1 and ?3B1 states of CF2 have been calculated. Core-electron effects on the computed minimum-energy geometries and relative electronic energies have been investigated, and relativistic contributions to the computed relative electronic energies calculated. Potential energy functions of the X2B1 state of CF2- and the X1A1 and ?3B1 states of CF2 have been determined, and anharmonic vibrational wavefunctions of these states calculated variationally. Franck-Condon factors including anharmonicity and Duschinsky rotation have been computed and used to simulate the ?-X emission spectrum of CF2 determined by S. Koda [Chem. Phys. Lett. 1978, 55, 353] and the 364 nm laser photodetachment spectrum of CF2- obtained by R. L. Schwartz et al. [J. Phys. Chem. A 1999, 103, 8213]. Comparison between theory and experiment shows that the theoretical approach benchmarked in the present study is able to give highly reliable positions for the CF2(X1A1) + e <-- CF2-(X2B1) and CF2(?3B1) + e <-- CF2-(X2B1) bands in the photoelectron spectrum of CF2- and a reliable singlet-triplet gap for CF2. It is therefore concluded that the same theoretical approach should give reliable simulated CCl2(X1A1) + e <-- CCl2-(X2B1) and CCl2(?3B1) + e <-- CCl2-(X2B1) bands in the photodetachment spectrum of CCl2- and a reliable singlet-triplet gap for CCl2.     

Rationalization of the behavior of M2(CH3CS2)4I (M = Ni,Pt) chains at room temperature from periodic density functional theory and ab initio cluster calculations

The electrical conductivities and plausible charge‐ordering states in the room temperature (r.t.) phase for MMX chains [Ni₂(dta)₄I]_∞ and [Pt₂(dta)₄I]_∞ (dta = CH₃CS) have been analyzed with periodic density functional theory (DFT) and correlated ab initio calculations combined with the effective Hamiltonian theory. Periodic DFT calculations show a more delocalized nature of the ground state in [Pt₂(dta)₄I]_∞ compared to [Ni₂(dta)₄I]_∞, which features a rather large energy gap between the occupied and empty bands, and charge polarized dimer units. A larger electrical conductivity for the Pt chain can be expected, especially because the Fermi level lies within a band with contributions from Pt and I orbitals. Electronic structure parameters extracted from ab initio cluster calculations show that the large difference between the observed conductivities at 300 K for Ni and Pt compounds, of 3 orders of magnitude, cannot be explained from the parameters extracted from an embedded M₂(dta)₄I₂ dimer fragment alone. When tetramer fragments are considered, we observe that the interdimer transfer integral (t) between neighboring M₂ units connected by an iodine atom at correlated level is comparable in both chains. On the other hand, the energy to transfer an electron from a dimer to the neighboring one (Coulomb repulsion U) is three times larger in the Ni compound with respect to the Pt chain, in line with the poor conductivity of the former. The electronic structure of the M₄(dta)₈I₃ fragment points to an alternate charge‐polarization state for Ni and an average valence state for Pt when the r.t. X‐ray structure is considered. © 2012 Wiley Periodicals, Inc.