Conformational behavior of simple furanosides studied by optical rotation

Experimental and theoretical specific optical rotations (OR) of anhydro, epithio, and epiminoderivatives of methyl tetrofuranosides in chloroform solutions have been compared and used as a tool for exploring their conformational behavior. The potential energy surfaces of these saccharides with reduced flexibility were examined with the density functional theory and the MP2 and CCSD(T) wavefunctions methods. Theoretical ORs were obtained by Boltzmann averaging of values calculated for local minima. Resultant rotations could be used to assess the quality of the DFT and MP2 relative conformer energies. OR values calculated for equilibrium geometries in vacuum were significantly improved when the solvent was accounted for by a polarizable continuum model and first and diagonal second OR derivatives were used for an anharmonic vibrational averaging. The DFT used as a default method reproduced the experimental data fairly well. A modified B3LYP functional containing 70% of HF exchange further improved the results. Because of the strong dependence of OR on the conformation, not only the absolute configuration could be determined, but also the conformational populations were estimated. Likewise, the predicted dependence of OR on the light wavelength well agreed with experiment. The increasing precision of the contemporary computational methods thus makes it possible to relate the specific rotation to more detailed features in molecular structure. © 2010 Wiley Periodicals, Inc. J Comput Chem 2010     

Structure and conformational stability of CH2CHCH2X (X=F, Cl and Br) molecules—post Hartree–Fock and density functional theory methods

The molecular structure and conformational stability of CH₂CHCH₂X (X=F, Cl and Br) molecules were studied using ab initio and density functional theory (DFT) methods. The molecular geometries of 3-fluoropropene were optimized employing BLYP and B3LYP levels of theory of DFT method implementing 6-311+G(d,p) basis set. The MP2/6-31G^*, BLYP and B3LYP levels of theory of ab initio and DFT methods were used to optimize the 3-chloropropene and 3-bromopropene molecules. The structural and physical parameters of the molecules are discussed with the available experimental values. The rotational potential energy surface of the above molecules were obtained at MP2/6-31G^* and B3LYP/6-311+G(d,p) levels of theory. The Fourier decomposition of the rotational potentials were analyzed. The HF/6-31G^* and MP2/6-31G^* levels of theory have predicted the cis conformer as the minimum energy structure for 3-fluoropropene, which is in agreement with the experimental values, whereas the BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory reverses the order of conformation. The ΔE values calculated for 3-chloropropene at MP2/6-31G^*, BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory show that the gauche form is more stable than the cis form, which is in agreement with the experimental value. The same levels of theory have also predicted that the gauche form is stable than cis for 3-bromopropene molecule. The maximum hardness principle has been able to predict the stable conformer of 3-fluoropropene at HF/6-31G^* level of theory, but the same level of theory reverses the conformational stability of 3-chloropropene and 3-bromopropene molecules and MP2/6-31G^* level of theory predicted the stable conformer correctly.     

Electronic structure studies of tetrazolium-based ionic liquids

New energetic ionic liquids are investigated as potential high energy density materials. Ionic liquids are composed of large, charge-diffuse cations, coupled with various (usually oxygen containing) anions. In this work, calculations have been performed on the tetrazolium cation with a variety of substituents. Density functional theory (DFT) with the B3LYP functional, using the 6-311G(d,p) basis set was used to optimize geometries. Improved treatment of dynamic electron correlation was obtained using second-order perturbation theory (MP2). Heats of formation of the cation with different substituent groups were calculated using isodesmic reactions and Gaussian-2 calculations on the reactants. The cation was paired with oxygen rich anions ClO4-, NO3-, or N(NO2)2- and those structures were optimized using both DFT and MP2. The reaction pathway for proton transfer from the cation to the anion was investigated.     

Gas phase study of the trans and gauche rotamers of 1,2-dicyanoethane, novel 1,2-dicyanodisilane and cyano(cyanomethyl)silane by ab initio and density functional methods

The gauche and trans rotamers of 1,2-dicyanoethane, novel 1,2-dicyanodisilane and cyano(cyanomethyl)silane have been studied theoretically in the gas phase. The methods used are second order M?ller-Plesset theory (MP2) and density functional theory (DFT). The basis set used is 6-311++G(d,p) for all atoms. B3LYP is the functional used for the DFT method and G2/MP2 calculation has also carried out using the MP2 optimised structure. All calculations have been done using Gaussian 03W. All structures have been fully optimised and the optimised geometries, dipole moments, moment of inertia and energies are reported. Energies of the optimised structures have been used to obtain the energy difference (DeltaE) between the trans and gauche rotamers. The optimised structures have been used for calculations of vibrational frequencies and these frequencies are reported with appropriate assignments. The computed parameters for 1,2-dicyanoethane compare satisfactorily with experimental literature values. However, the literature for 1,2-dicyanodisilane and cyano(cyanomethyl)silane, in terms of conformational studies, is limited and therefore the data of this work should also be appropriate for them. The results indicate that in general, the energy difference for these molecules is in the order 1,2-dicyanoethane>cyano(cyanomethyl)silane>1,2-dicyanodisilane.     

High level of ab initio and density functional theory computational study of structural and energetic properties of tetrafluorhydrazine rotamers

The HF, MP2, MP3, MP4, and QCISD ab initio methods were compared with local, hybrid, and gradient-corrected density functional theory (DFT) methods for computing structures and energies of N2F4 rotamers. In all DFT calculations 6-311 + G(2d) basis set was used. The generated structures energies of trans- and gauche-N₂F₄ rotamers, and their dissociation energies to nitrogen difluoride were compared with experimental data. Suitable hybrid and gradient-corrected DFT methods for determining structures and energies for these and similar molecular systems were discussed.     

Coupled cluster calculations of optical rotatory dispersion of (S)-methyloxirane

Coupled cluster (CC) and density-functional theory (DFT) calculations of optical rotation, [alpha](lambda), have been carried out for the difficult case of (S)-methyloxirane for comparison to recently published gas-phase cavity ringdown polarimetry data. Both theoretical methods are exquisitely sensitive to the choice of one-electron basis set, and diffuse functions have a particularly large impact on the computed values of [alpha](lambda). Furthermore, both methods show a surprising sensitivity to the choice of optimized geometry, with [alpha](355) values varying by as much as 15 deg dm(-1) (g/mL)(-1) among molecular structures that differ only negligibly. Although at first glance the DFT/B3LYP values of [alpha](355) appear to be superior to those from CC theory, the success of DFT in this case appears to stem from a significant underestimation of the lowest (Rydberg) excitation energy in methyloxirane, resulting in a shift of the first-order pole in [alpha](lambda) (the Cotton effect) towards the experimentally chosen incident radiation lines. This leads to a fortuitous positive shift in the value of [alpha](355) towards the experimental result. The coupled cluster singles and doubles model, on the other hand, correctly predicts the position of the absorption pole (to within 0.05 eV of the experimental result), but fails to describe correctly the shape/curvature of the ORD region lambda=355, resulting in an incorrect prediction of both the magnitude and the sign of the optical rotation.     

Stabilization energies of the hydrogen-bonded and stacked structures of nucleic acid base pairs in the crystal geometries of CG, AT, and AC DNA steps and in the NMR geometry of the 5'-d(GCGAAGC)-3' hairpin: Complete basis set calculations at the MP2 and CCSD(T) levels

Stabilization energies of the H-bonded and stacked structures of a DNA base pair were studied in the crystal structures of adenine-thymine, cytosine-guanine, and adenine-cytosine steps as well as in the 5'-d(GCGAAGC)-3' hairpin (utilizing the NMR geometry). Stabilization energies were determined as the sum of the complete basis set (CBS) limit of MP2 stabilization energies and the Delta E(CCSD(T)) - Delta E(MP2) correction term evaluated with the 6-31G*(0.25) basis set. The CBS limit was determined by a two-point extrapolation using the aug-cc-pVXZ basis sets for X = D and T. While the H-bonding energies are comparable to those of base pairs in a crystal and a vacuum, the stacking energies are considerably smaller in a crystal. Despite this, the stacking is still important and accounts for a significant part of the overall stabilization. It contributes equally to the stability of DNA as does H-bonding for AT-rich DNAs, while in the case of GC-rich DNAs it forms about one-third of the total stabilization. Interstrand stacking reaches surprisingly large values, well comparable to the intrastrand ones, and thus contributes significantly to the overall stabilization. The hairpin structure is characterized by significant stacking, and both guanine...cytosine pairs possess stacking energies larger than 11.5 kcal/mol. A high portion of stabilization in the studied hairpin comes from stacking (similar to that found for AT-rich DNAs) despite the fact that it contains two GC Watson-Crick pairs having very large H-bonding stabilization. The DFT/B3LYP/6-31G** method yields satisfactory values of interaction energies for H-bonded structures, while it fails completely for stacking.     

Solvent effects on the S(N)2 reaction: Application of the density functional theory-based effective fragment potential method

The performance of the density functional theory (DFT)-based effective fragment potential (EFP) method is assessed using the S(N)2 reaction: Cl- + nH2O + CH3Br = CH3Cl + Br- + nH2O. The effect of the systematic addition of water molecules on the structures and relative energies of all species involved in the reaction has been studied. The EFP1 method is compared with second-order perturbation theory (MP2) and DFT results for n = 1, 2, and 3, and EFP1 results are also presented for four water molecules. The incremental hydration effects on the barrier height are the same for all methods. However, only full MP2 or MP2 with EFP1 solvent molecules are able to provide an accurate treatment of the transition state (TS) and hence the central barriers. Full DFT and DFT with EFP1 solvent molecules both predict central barriers that are too small. The results illustrate that the EFP1-based DFT method gives reliable results when combined with an accurate quantum mechanical (QM) method, so it may be used as an efficient alternative to fully QM methods in the treatment of larger microsolvated systems.     

Structure, stability, thermodynamic properties, and infrared spectra of the protonated water octamer H(+)(H2O)8

We investigated various two-dimensional (2D) and three-dimensional (3D) structures of H (+)(H 2O) 8, using density functional theory (DFT), Moller-Plesset second-order perturbation theory (MP2), and coupled cluster theory with single, double, and perturbative triple excitations (CCSD(T)). The 3D structure is more stable than the 2D structure at all levels of theory on the Born-Oppenheimer surface. With the zero-point energy (ZPE) correction, the predicted structure varies depending on the level of theory. The DFT employing Becke's three parameters with Lee-Yang-Parr functionals (B3LYP) favors the 2D structure. At the complete basis set (CBS) limit, the MP2 calculation favors the 3D structure by 0.29 kcal/mol, and the CCSD(T) calculation favors the 3D structure by 0.27 kcal/mol. It is thus expected that both 2D and 3D structures are nearly isoenergetic near 0 K. At 100 K, all the calculations show that the 2D structure is much more stable in free binding energy than the 3D structure. The DFT and MP2 vibrational spectra of the 2D structure are consistent with the experimental spectra. First-principles Car-Parrinello molecular dynamics (CPMD) simulations show that the 2D Zundel-type vibrational spectra are in good agreement with the experiment.     

The microwave spectrum of a two-top peptide mimetic: the N-acetyl alanine methyl ester molecule

The rotational spectrum of N-acetyl alanine methyl ester, a derivative of the biomimetic, N-acetyl alanine N'-methyl amide or alanine dipeptide, has been measured using a mini Fourier transform spectrometer between 9 and 25 GHz as part of a project undertaken to determine the conformational structures of various peptide mimetics from the torsion-rotation parameters of low-barrier methyl tops. Torsion-rotation splittings from two of the three methyl tops capping the acetyl end of the -NH-C(=O)- and the methoxy end of -C(=O)-O- groups account for most of the observed lines. In addition to the AA state, two E states have been assigned and include an AE state having a torsional barrier of 396.45(7) cm(-1) (methoxy rotor) and an EA state having a barrier of 64.96(4) cm(-1) (acetyl rotor). The observed torsional barriers and rotational constants of alanine dipeptide and its methyl ester are compared with predictions from M?ller-Plesset second-order perturbation theory (MP2) and density functional theory (DFT) in an effort to explore systematic errors at the two levels of theory. After accounting for zero-point energy differences, the torsional barriers at the MP2/cc-pVTZ level are in excellent agreement with experiment for the acetyl and methoxy groups while DFT predictions range from 8% to 80% too high or low. DFT is found to consistently overestimate the overall molecular size while MP2 methods give structures that are undersized. Structural discrepancies of similar magnitude are evident in previous DFT results of crystalline peptides.     

Treating dispersion effects in extended systems by hybrid MP2:DFT calculations--protonation of isobutene in zeolite ferrierite

We propose use of a hybrid method to study problems that involve both bond rearrangements and van-der-Waals interactions. The method combines second-order M?ller-Plesset perturbation theory (MP2) calculations for the reaction site with density functional theory (DFT) calculations for a large system under periodic boundary conditions. Hybrid MP2:DFT structure optimisation for a cluster embedded in the periodic model is the first of three steps in a multi-level approach. The second step is extrapolation of the MP2 energy to the complete basis set limit. The third step is extrapolating the high-level (MP2) correction to the limiting case of the full periodic structure. This is done by calculating the MP2 correction for a series of cluster models of increasing size, fitting an analytic expression to these energy corrections, and applying the fitted expression to the full periodic structure. We assume that, up to a constant, the high-level correction is described by a damped dispersion expression. Combining the results of all three steps yields an estimate of the MP2 reaction energy for the full periodic system at the complete basis set level. The method is designed for a reaction between a small or medium sized substrate molecule and a very large chemical system. For adsorption of isobutene in zeolite H-ferrierite, the energies obtained for the formation of different structures, the pi-complex, the isobutoxide, the tert-butoxide, and the tert-butyl carbenium ion, are -78, -73, -48, and -21 kJ mol(-1), respectively. This corresponds to corrections of the pure DFT (PBE functional) results by -62, -70, -67, and -29 kJ mol(-1), respectively. Hence, the MP2 corrections are substantial and, perhaps more importantly, not the same for the different hydrocarbon species in the zeolite. Coupled-cluster (CCSD(T)) calculations change the MP2 energies by -4 kJ mol(-1) (tert-butyl cation) or less (below +/-1 kJ mol(-1) for the other species).     

Structures and properties of stable Al4, Al4+, and Al4- comparatively studied by ab initio theories

Previous high-level theoretical calculations of aluminum clusters mostly relied on density functional theory (DFT) or theories less sophisticated than it. Here, we point out that a second-order M?ller-Plesset perturbation (MP2) method is more appropriate and is the minimum level of theory in the predictions of property such as geometries, electronic structures, and IR and Raman spectra of Al4, Al4+, and Al4- clusters. The theoretical electron affinities and ionization potentials predicted with MP2 geometries are closer to experimental ones than those by DFT. The p-electron characters and single valence state of aluminum atoms and IR and Raman spectra of the aluminum clusters were also reliably predicted by MP2 and could be based on for further experimental and theoretical studies.     

Performance of theoretical methods and basis sets on the molecular structure, atomisation and ionisation energies, electron affinity, and vibrational spectrum of silylene

This work compares the performance of theoretical methods and basis sets on the molecular structure, atomisation and ionisation energies, electron affinity, and vibrational spectrum of silylene. Silylene, its cation and anion have been studied in ¹ A ₁, ² A ₁ and ² B ₁ states, respectively, in the gas phase and C_2v symmetry. The methods considered are second-order Møller-Plesset perturbation theory (MP2), the density functional theory (DFT), Gaussian-2 (G2) and complete basis set methods (CBS-4M and CBS-Q). The basis sets used are 6-31G(d,p), 6-311G(d,p), 6-31++G(d,p) and 6-311++G(d,p). The functional used for the DFT method is B3LYP. Silylene and its cation and anion have been optimised using the MP2 and DFT methods and the named basis sets. Single-point energy calculations (G2, CBS-4M and CBS-Q) were performed using MP2/6-311++G(d,p) structures and these energies have been used to calculate atomisation energy, ionisation energy and adiabatic electron affinity. Frequency calculations were also done and the raw vibrational frequencies were assigned. It is interesting to note the close similarity between the predicted parameters and some of the available literature values. The results obtained are consistent and converge with different basis sets with improved size and quality. However, the parameters obtained are very much method dependent.     

Highly Accurate CCSD(T) and DFT–SAPT Stabilization Energies of H‐Bonded and Stacked Structures of the Uracil Dimer

The CCSD(T) interaction energies for the H‐bonded and stacked structures of the uracil dimer are determined at the aug‐cc‐pVDZ and aug‐cc‐pVTZ levels. On the basis of these calculations we can construct the CCSD(T) interaction energies at the complete basis set (CBS) limit. The most accurate energies, based either on direct extrapolation of the CCSD(T) correlation energies obtained with the aug‐cc‐pVDZ and aug‐cc‐pVTZ basis sets or on the sum of extrapolated MP2 interaction energies (from aug‐cc‐pVTZ and aug‐cc‐pVQZ basis sets) and extrapolated ΔCCSD(T) correction terms [difference between CCSD(T) and MP2 interaction energies] differ only slightly, which demonstrates the reliability and robustness of both techniques. The latter values, which represent new standards for the H‐bonding and stacking structures of the uracil dimer, differ from the previously published data for the S22 set by a small amount. This suggests that interaction energies of the S22 set are generated with chemical accuracy. The most accurate CCSD(T)/CBS interaction energies are compared with interaction energies obtained from various computational procedures, namely the SCS–MP2 (SCS: spin‐component‐scaled), SCS(MI)–MP2 (MI: molecular interaction), MP3, dispersion‐augmented DFT (DFT–D), M06–2X, and DFT–SAPT (SAPT: symmetry‐adapted perturbation theory) methods. Among these techniques, the best results are obtained with the SCS(MI)–MP2 method. Remarkably good binding energies are also obtained with the DFT–SAPT method. Both DFT techniques tested yield similarly good interaction energies. The large magnitude of the stacking energy for the uracil dimer, compared to that of the benzene dimer, is explained by attractive electrostatic interactions present in the stacked uracil dimer. These interactions force both subsystems to approach each other and the dispersion energy benefits from a shorter intersystem separation.     

Conformer Pair Contributions to Optical Rotations in a Series of Chiral Linear Aliphatic Alcohols

The chain length effect of four chiral aliphatic alcohols,(S)-2-butanol,(S)-2-pentanol,(S)-2-hexanol and (S)-2-heptanol,on their specific optical rotations(OR)was studied experimentally and theoretically via quantum theory.Many conformations of each chiral alcohol exist as conformer pairs in solution.The OR sum from these pairs of conformers has much smaller contributions to OR values than that contributed by the most stable conformation. These four alcohols'OR values were also investigated using the matrix model,which employs each substituent's comprehensive mass,radii,electronegativity and symmetry number as the elements in the matrix.These are all particle properties.This matrix determinant is proportional to its OR values within a closely related structural series of chiral compounds.The experimental OR values and the matrix determinants of these four alcohols were compared with the predicted OR values obtained from quantum theory wave functions.The ORs predicted by the matrix method, which is based on particle function statistics,agreed with the results from quantum theory.The agreement between OR predictions by the matrix method and DFT calculations illustrates the wave-particle duality of polarized light that is operating in these predictions.     

Torsional potential of 4,4'-bipyridine: ab initio analysis of dispersion and vibrational effects

Ab initio calculations using restricted Hartree-Fock, second-order M?ller-Plesset perturbation theory (MP2), density-functional theory (DFT), and coupled-cluster methods have been done to obtain the torsional potential-energy profile of the aza-aromatic molecule 4,4'-bipyridine. The torsional potential is evaluated adiabatically by fixing the normalized sum of the dihedral angles through the C-C inter-ring bond at several values along the torsional path and relaxing the remaining degrees of freedom. Previous discrepancies between MP2 and DFT internal rotation barrier heights are removed, and seen to be mostly due to the underestimation of the dispersion energy in the coplanar conformer by MP2 when using relatively small basis sets. The calculations indicate that the barrier height between the twisted global minimum and the 0 degrees conformer is around 1.5-1.8 kcal mol-1 while that corresponding to the 90 degrees one is about 2.0-2.2 kcal mol-1. This same relative energy ordering of the coplanar and perpendicular conformers was experimentally derived from nuclear magnetic resonance (NMR) measurements of 1H dipolar couplings on 4,4'-bipyridine solutions in a nematic liquid crystal, although the barrier heights are much lower than those estimated from NMR experiments in the gas phase. The DFT infrared spectrum and zero-point vibrational energy corrections to the torsional energy profile have also been calculated, the latter having a small influence on the torsional potential-energy profiles.     

A theoretical study of nonlinearity in heterocyclic hydrogen-bonded complexes

Ab initio calculations are performed at the MP2/6-311++G(d,p) and DFT/B3LYP/6-311++G(d,p) theoretical levels to obtain geometries, H-bond energies and harmonic infrared vibrational properties for the Cs symmetry structures of heterocyclic hydrogen-bonded complexes, CnHmY-HX. The H-bond lengths in DFT/B3LYP calculation level are in better agreement with the experimental values than the MP2 results. The geometry optimization are interpreted in terms of hydrogen bond nonlinearity represented by theta; and phi angles, once the hydrogen bond is formed among n-electrons pairs of the heteroatom in heterocyclic and the hydrogen atom in HX. The hydrogen bond energy after of the zero-point vibrational energy (ZPE) and basis set superposition error (BSSE) corrections are overestimated at DFT/B3LYP, whereas the MP2 BSSE corrections are very large than corresponding DFT/B3LYP. For example, the BSSE corrections for the C2H4S-HNC complex are 7.60 and 0.09 kJ mol(-1) in MP2 and DFT/B3LYP calculations levels, respectively. The new vibrational modes in infrared harmonic spectrum arising from complexation show several interesting features, especially the intermolecular stretching mode.     

How Important Is Molecular Rigidity for the Complex Stability of Artificial Host–Guest Systems? A Theoretical Study on Self‐Assembly of Gas‐Phase Arginine

Arginine forms much less stable dimers than 2-(guanidiniocarbonyl)-1H-pyrrole-5-carboxylate although the principal binding interactions are very similar. The reasons for this difference are addressed in this work by state-of-the-art ab initio computations. The investigation shows that the extraordinary high stability of the 2-(guanidiniocarbonyl)-1H-pyrrole-5-carboxylate dimer results to about 50 % from the rigidity of its monomer. Within this study monomer and dimer conformers of arginine were calculated leading to new low lying structures which have not been reported before as well as new global minima are predicted. In these structures stacking interactions with the guanidinium moiety are especially important. For the monomer we predict the energy minimum to be the canonical form with the lowest lying zwitterionic structure being only 9 kJ mol(-1) less stable. During the course of these calculations we found that DFT did not predict the structures and their relative energy correctly in comparison to perturbation theory (MP2) and some potential reasons for the failure of DFT in these cases are discussed. Vibrational frequencies of the various structures are presented and a suitable wavenumber region for an experimental determination of the global minimum of the arginine monomer is identified. The effect of molecular rigidity on the self-assembly is probed using a local minimum of the arginine monomer which does not possess any intramolecular stabilizing effects. Our results suggest that the deliberate control of the conformational flexibility is a powerful instrument to steer the complex affinity of artificial hosts.     

Energetics and structure of glycine and alanine based model peptides: Approximate SCC-DFTB,AM1 and PM3 methods in comparison with DFT,HF and MP2 calculations

We calculate relative energies and geometries of important secondary structural elements for small glycine and alanine based polypeptides containing up to eight residues. We compare the performance of the approximate methods AM1, PM3 and self-consistent charge, density-functional tight-binding (SCC-DFTB) to density-functional theory (DFT), Hartree–Fock (HF) and MP2. The SCC-DFTB is able to reproduce structures and relative energies of various peptide models reliably compared to DFT results. The AM1 and PM3 methods show deficiencies in describing important secondary structure elements like extended, helical or turn structures. The discrepancies between different ab initio (HF, MP2) and DFT (B3LYP) methods for medium sized basis sets (6-31G*) also show the need for higher level calculations, since systematic errors found for small molecules may add up when investigating longer polypeptides.     

The determination of the equilibrium structures of oxygen,ozone, and hydrogen peroxide using the ab initio and density functional theory methods

The geometries and energies of small oxygen containing molecules are studied by both the ab initio and density functional theory (DFT) methods. The RHF, MP2, and QCISD(T) ab initio methods, BHandH, BHandHLYP, BeckeSLYP, Becke3P86 DFT hybrid methods, BLYP, and the BP86 non-local DFT methods with the 3-21G¹, 6-31G(d,p), 6-311 + G(2d,2p) and 6-311 + + G(3df,3pd) basis sets were used for the computational study. The obtained results from the different methods were compared to the experimental values. The suitability of the DFT methods for reproducing experimental data were discussed.