Ab Initio calculation of 19F NMR chemical shielding for alkaline‐earth‐metal fluorides

Gauge‐independent atomic orbital (GIAO) method at Hartree‐Fock (HF) and density functional theory (DFT) levels, respectively, was employed to calculate ¹⁹F NMR chemical shieldings of solid state alkaline‐earth‐metal fluorides MF₂ (M = Mg, Ca, Sr, Ba). The results show that, although the calculated ¹⁹F chemical shieldings tend to be larger than the experimental values, they have a fairly good linear relationship with the observed ones. The calculated results based on different combinations of basis sets show that the B3LYP (a hybrid of DFT with HF) predictions are greatly superior to the HF predictions. When a basis set of metal atom with effective core potential (ECP) has well representation of valence wavefunction, especially wavefunction of d component, and proper definition of core electron number, it can be applied to obtain ¹⁹F chemical shielding which is close to that of all‐electron calculation. The variation of ¹⁹F chemical shielding of alkaline‐earth‐metal fluorides correlates well with the lattice factor A/R².     

Ab initio investigation of the Young's modulus of polyamide‐6

The Young's modulus of the α form of polyamide‐6 has been calculated using the supermolecule model. The crystalline polymer was represented by a single‐chain molecule, divided into a head, body, and tail part. The body of the model contains an even number of polyamide‐6 units (4–16 units) and is representative for a polyamide‐6 chain. The periodicity of the system was not explicitly taken into account, but in a few tests the effect of a linear constraint has been evaluated. An n‐butyl and n‐pentyl group have been used as head and tail, respectively. In a number of additional calculations the length of the head and tail has been varied. For all supermolecule models the equilibrium and elongated structures have been optimized using ab initio Hartree–Fock calculations with a 6‐31G** basis set. From the energy values of the optimized structure a Young's modulus of 334 GPa has been extrapolated for both the unconstrained and linearly constrained models. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

Ab initio quantum chemical calculations of a cluster C8H12

Using unrestricted Hartree–Fock–Roothaan approximation and density functional theory on 6‐31G* basis, ab initio calculation cluster C₈H₁₂ has been carried out. It was shown as a result of calculations that a ground state for the conformation given is a septet state. This cluster can be used for investigation of the kinematic mechanism of magnetic exchange and magnetic ordering in polyradicals. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002     

Ab initio fragment orbital‐based theory

A new formulation of ab initio theory is presented that treats a large molecule in terms of wave functions of its constituent molecular subunits (to be called fragments). The method aims to achieve near conventional ab initio accuracy but using a truncated set of fragment orbitals with a consequent drastic reduction of computing time and storage requirement. Illustrative calculations are presented for the molecule amino‐nitro‐stilbene. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

Ab initio prediction of the vibrational spectra of the hydrogen‐bonded complexes between CO and HONO2

The vibrational characteristics (vibrational frequencies and infrared intensities) for free and complexed CO and HONO₂ have been predicted using ab initio calculations at SCF and MP2 levels with different basis sets and B3LYP/6?31G(d,p) calculations. The ab initio calculations show that the complexation between HONO₂ and CO leads to two stable structures: CO … HONO₂ (1A) and OC … HONO₂ (1B). The changes in the vibrational characteristics from free monomers to complexes have been estimated. It was established that the most sensitive to the complexation is the stretching O? H vibration. In agreement with the experiment, its vibrational frequency in the complexes is shifted to lower frequency (Δν = ?123 cm^?1). The magnitude of the wave number shift is indicative of relatively strong hydrogen‐bonded interaction. The ab initio calculations at different levels predict an increase of the infrared intensity of the stretching O? H vibration for structure 1A more than five times and for structure 1B more than nine times. The most consistent agreement between the computed values of the frequency shifts for structure 1B and those experimentally observed suggests that this structure is preferred. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

Ab initio Studies on Li4+xTi5O12 Compounds as Anode Materials for Lithium‐Ion Batteries

The structural and electronic properties of Li_4+xTi₅O₁₂ compounds (with 0≤x≤6)—to be used as anode materials for lithium‐ion batteries—are studied by means of first principles calculations. The results suggest that Li₄Ti₅O₁₂ can be lithiated to the state Li_8.5Ti₅O₁₂, which provides a theoretical capacity that is about 1.5 times higher than that of the compound lithiated to Li₇Ti₅O₁₂. Further insertion of lithium species into the Li_8.5Ti₅O₁₂ lattice results in a clear structural distortion. The small lattice expansion observed upon lithium insertion (about 0.4 % for the lithiated material Li_8.5Ti₅O₁₂) and the retained [Li₁Ti₅]^16dO₁₂ framework indicate that the insertion/extraction process is reversible. Furthermore, the predicted intercalation potentials are 1.48 and 0.05 V (vs Li/Li⁺) for the Li₄Ti₅O₁₂/Li₇Ti₅O₁₂ and Li₇Ti₅O₁₂/Li_8.5Ti₅O₁₂ composition ranges, respectively. Electronic‐structure analysis shows that the lithiated states Li_4+xTi₅O₁₂ are metallic, which is indicative of good electronic‐conduction properties.     

Ab initio calculation of the interaction energy in the P2 binding pocket of HIV‐1 protease

Interactions at the P2 binding pocket of human immunodeficiency virus type 1 (HIV‐1) protease have been studied using calculated interaction energies for model systems that mimic this binding pocket. Models were built for the P2 pocket of HIV‐1 protease in complex with TMC114, nelfinavir, and amprenavir. A two‐step procedure was applied. In the first step, the size of the model system was confined to ～40 atoms, and the interaction energy was calculated at different computational levels. In the second step, the size of the system was increased to 138 atoms, and the calculations were only performed at the HF/6‐31G** level. The interaction energy of the HIV‐1 protease/TMC114 complex was found to be more favorable than the interaction energies of the other complexes because of the additional hydrogen bond interaction this inhibitor is able to make with the HIV‐1 protease backbone. The results of the calculations are supported by stockholder charges and electrostatic potential maps. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Ab initio study of C4H3 potential energy surface and reaction of ground‐state carbon atom with propargyl radical

The potential energy surface for the reaction of the ground‐state carbon atom [C(³P_j)] with the propargyl radical [HCCCH₂(X²B₁)] is investigated using the G2M(RCC,MP2) method. Numerous local minima and transition states for various isomerization and dissociation pathways of doublet C₄H₃ are studied. The results show that C(³P_j) attacks the π system of the propargyl radical at the acetylenic carbon atom and yields the n‐C₄H₃(²A′) isomer i3 after an 1,2‐H atom shift. This intermediate either splits a hydrogen atom and produces singlet diacetylene, [HCCCCH ( p1 )+H] or undergoes (to a minor amount) a 1,2‐H migration to i‐C₄H₃(²A′) i5 , which in turn dissociates to p1 plus an H atom. Alternatively, atomic carbon adds to the triple C?C bond of the propargyl radical to form a three‐member ring C₄H₃ isomer i1 , which ring opens to i3 . Diacetylene is concluded to be a nearly exclusive product of the C(³P_j)+HCCCH₂ reaction. At the internal energy of 10.0 kcal/mol above the reactant level, Rice–Ramsperger–Kassel–Marcus calculations show about 91.7% of HCCCCH comes from fragmentation of i3 and 8.3% from i5 . The other possible minor channels are identified as HCCCC+H₂ and C₂H+HCCH. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1522–1535, 2001     

Ab initio multireference investigation of disjoint diradicals: singlet versus triplet ground states

We present improved virtual orbital (IVO) complete active space (CAS) configuration interaction (IVO‐CASCI) and IVO‐CASCI‐based multireference Møller–Plesset perturbation theory (MRMPPT) calculations with an aim to elucidate the electronic structure of tetramethyleneethane (TME) in its lowest singlet and triplet state and to quantify their order and extent of splitting. The potential surfaces of singlet and triplet states for the twisting of TME are also studied. We found that the triplet state is higher in energy than the singlet one in the whole range of twisting angles with the energy gap minimum at a twisting angle of about 45°. Harmonic vibrational frequencies of TME have also been calculated for both the states. We also report the ground to first excited triplet state transition energies. Our results are analyzed with respect to the results available in the literature to illustrate the efficacy of our methods employed. We also demonstrate that the spin character of the ground state of disjoint, TME‐like diradicals can be manipulated by using appropriate selection of annulenic spacer to separate the allyl groups of TME.     

Ab initio and DFT study of the electronic structures and spectroscopic properties of pyrene ligands and their cyclometalated complexes

Two ligands 1‐diphenylphosphinopyrene (1‐PyP) ( L ₁ ), 1,6‐bis(diphenylphosphino)‐pyrene (1,6‐PyP) ( L ₂ ) and their cyclometalated complexes [Pt(dppm)(1‐PyP‐H)]⁺ ( 1 ), [Pt₂(dppm)₂(1,6‐PyP‐H₂)]²⁺ (dppm = bis(diphenylphosphino)methane ( 2 ), and [Pd(dppe)(1‐PyP‐H)⁺ (dppe = bis(diphenylphosphino)ethane) ( 3 ) are investigated theoretically to explore their electronic structures and spectroscopic properties. The ground‐ and excited‐state structures are optimized by the density functional theory (DFT) and single‐excitation configuration interaction method, respectively. At the time‐dependent DFT (TDDFT) and B3LYP level, the absorption and emission spectra in solution are obtained. As revealed from the calculations, the lowest‐energy absorptions of 1 and 3 are attributed to the mixing ligand‐to‐metal charge transfer (CT)/intraligand (IL)/ligand‐to‐ligand CT transitions, while that of 2 is attributed to the IL transition. The lowest‐energy phosphorescent emissions of the cyclometalated complexes are attributed to coming from the ³ILCT transitions. With the increase of the spin‐orbit coupling effect, the phosphorescence intensities and the emissions wavelength are correspondingly increased. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Ab initio study of the transition-metal carbene cations

The geometries and bonding characteristics of the first-row transition-metal carbene cations MCH_2~ were investigated by ab initio molecular orbital theory (HF/LANL2DZ). All of MCH_2~ are coplanar. In the closed shell structures the C bonds to M with double bonds; while in the open shell structures the partial double bonds are formed, because one of the σ and π orbitals is singly occupied. It is mainly the π-type overlap between the 2p_x orbital of C and 4p_x, 3d_(xz), orbitals of M~ that forms the π orbitals. The dissociation energies of C—M bond appear in periodic trend from Sc to Cu. Most of the calculated bond dissociation energies are close to the experimental ones.     

Carbon‐Centered Radical Addition and β‐Scission Reactions: Modeling of Activation Energies and Pre‐exponential Factors

A consistent set of group additive values ΔGAV° for 46 groups is derived, allowing the calculation of rate coefficients for hydrocarbon radical additions and β-scission reactions. A database of 51 rate coefficients based on CBS-QB3 calculations with corrections for hindered internal rotation was used as training set. The results of this computational method agree well with experimentally observed rate coefficients with a mean factor of deviation of 3, as benchmarked on a set of nine reactions. The temperature dependence on the resulting ΔGAV°s in the broad range of 300–1300 K is limited to ±4.5 kJ mol⁻¹ on activation energies and to ±0.4 on logA (A: pre-exponential factor) for 90 % of the groups. Validation of the ΔGAV°s was performed for a test set of 13 reactions. In the absence of severe steric hindrance and resonance effects in the transition state, the rate coefficients predicted by group additivity are within a factor of 3 of the CBS-QB3 ab initio rate coefficients for more than 90 % of the reactions in the test set. It can thus be expected that in most cases the GA method performs even better than standard DFT calculations for which a deviation factor of 10 is generally considered to be acceptable.     

Ab initio and DFT investigation of the structures and properties of chloromethyl and chlorofluoromethyl peroxyl radicals

Ab initio and Density Functional Theory (DFT) calculations were performed to determine the equilibrium geometries, charge distributions, spin density distributions, dipole moments, electron affinities (EAs), and C–O bond dissociation energies (BDEs) of CH₂ClO₂^? CHCl₂O₂^?, CCl₃O₂^?, CF₂ClO₂^?, CFCl₂O₂^?, and CHFClO₂^? peroxyl radicals. The C–H BDEs of the parent methanes were calculated using the same levels of theories. Both MP2(full) and B3LYP methods, using the 6‐31G(d,p) basis set, were found to be capable of accurately predicting the geometries of peroxyl radicals. The B3LYP/6‐31G(d,p) method was found to be comparable to high ab initio levels in predicting C–O BDEs of studied peroxyl radicals and C–H BDEs of the parent alkanes. The progressive chlorine substitution of hydrogen atoms in methyl peroxyl radicals results in an increase (decrease) of the spin density on the terminal (inner) oxygen, a decrease in dipole moments, and an increase in electron affinities. The substitution of fluorine by chlorine in the series CF₃O₂^? – CCl₃O₂^? was found to lengthen (destabilize) the C–O bonds. Both C–O BDEs and EAs of peroxyl radicals (RO₂^?) were found to correlate well with Taft σ* substituent constants for the R groups. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Ab initio study of reactions of hydroxyl radicals with chloro‐ and fluoro‐substituted methanes

Ab initio calculations at the unrestricted Hartree–Fock (UHF) level have been performed to investigate the hydrogen abstraction reactions of ^? OH radicals with methane and nine halogen‐substituted methanes (F, Cl). Geometry optimization and vibrational frequency calculations have been performed on all reactants, adducts, products, and transition states at the UHF/6‐31G* level. Single‐point energy calculations at the MP2/6‐31++G* level using the UHF/6‐31G* optimized geometries have also been carried out on all species. Pre‐ and postreaction adducts have been detected on the UHF/6‐31G* potential energy surfaces of the studied reactions. Energy barriers, ΔE^?, reaction energies, ΔE_r, reaction enthalpies, ΔH_r, and activation energies, E_a, have been determined for all reactions and corrected for zero‐point energy effects. Both E_a and ΔH_r come into reasonable agreement with the experiment when correlation energy is taken into account and when more polarized and diffuse basis sets are used. The E_a values, estimated at the PMP2/6‐31++G* level, are found to be in good agreement with the experimental ones and correctly reproduce the experimentally observed trends in fluorine and chlorine substitution effects. A linear correlation between E_a and ΔH_r is obtained, suggesting the presence of an Evans–Polanyi type of relationship. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem 84: 426–440, 2001     

Conformational properties of methyl vinyl sulfone: Ab initio geometry optimization and vibrational analysis

An ab initio conformational analysis of methyl vinyl sulfone (CH₂CHSO₂CH₃) has been carried out. Molecular geometry optimizations have been performed at the HF and MP2 levels of the theory. Relative energies of the stationary points have been determined by using different approaches, including electron correlation corrections up to the third order. The IR and Raman spectra of the liquid have been measured and a vibrational assignment is proposed. The height of the barrier to the methyl group internal rotation has been estimated. Theoretical calculations and vibrational spectra have shown that the predominant conformation of methyl vinyl sulfone has the C=C bond eclipsed with one of the S=O bonds. Similar eclipsed forms have been found in vinyl fluoro sulfone, vinyl chloro sulfone and divinyl sulfone by ab initio HF calculations.     

Wave‐function frozen‐density embedding: Approximate analytical nuclear ground‐state gradients

We report the derivation of approximate analytical nuclear ground‐state uncoupled frozen density embedding (FDEu) gradients for the resolution of identity (RI) variant of the second‐order approximate coupled cluster singles and doubles (RICC2) as well as density functional theory (DFT), and an efficient implementation thereof in the KOALA program. In order to guarantee a computationally efficient treatment, those gradient terms are neglected which would require the exchange of orbital information. This approach allows for geometry optimizations of single molecules surrounded by numerous molecules with fixed nuclei at RICC2‐in‐RICC2, RICC2‐in‐DFT, and DFT‐in‐DFT FDE level of theory using a dispersion correction, required due to the DFT‐based treatment of the interaction in FDE theory. Accuracy and applicability are assessed by the example of two case studies: (a) the Watson‐Crick pair adenine‐thymine, for which the optimized structures exhibit a maximum error of about 0.08 Å for our best scheme compared to supermolecular reference calculations, (b) carbon monoxide on a magnesium oxide surface model, for which the error amount up to 0.1 Å for our best scheme. Efficiency is demonstrated by successively including environment molecules and comparing to an optimized conventional supermolecular implementation, showing that the method is able to outperform conventional RICC2 schemes already with a rather small number of environment molecules, gaining significant speed up in computation time. © 2016 Wiley Periodicals, Inc.     

Peptide models XXXI. Conformational properties of hydrophobic residues shaping the core of proteins. An ab initio study of N‐formyl‐L‐valinamide and N‐formyl‐L‐phenylalaninamide

Employing introductory (3‐21G RHF) and medium‐size (6‐311++G** B3LYP) ab initio calculations, complete conformational libraries, containing as many as 27 conformers, have been determined for diamide model systems incorporating the amino acids valine (Val) and phenylalanine (Phe). Conformational and energetic properties of these libraries were analyzed. For example, significant correlation was found between relative energies from 6‐311++G** B3LYP and single‐point B3LYP/6‐311++G**//RHF/3‐21G calculations. Comparison of populations of molecular conformations of hydrophobic aromatic and nonaromatic residues, based on their ab initiorelative energies, with their natural abundance indicates that, at least for the hydrophobic core of proteins, the conformations of Val (Ile, Leu) and Phe (Tyr, Trp) are controlled by the local energetic preferences of the respective amino acids. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 732–751, 2001     

Lithium hydroxide phase transition under high pressure: an ab initio molecular dynamics study.

The high-pressure phase transition in the deuterated lithium hydroxide crystalline state has been studied by Car-Parrinello molecular dynamics simulations, in the constant-pressure, constant-temperature ensemble. The recently developed metadynamics approach has been applied to encourage the system to transform into different phases in an affordable simulation time. A previously not completely characterized high-pressure phase has been obtained. The structural and spectroscopic properties have been studied and compared with the neutron scattering, infrared and Raman measurements. It has been found that the calculated structure differs slightly from the experimental hypothesis, and that the presence of strong hydrogen bonds is the source of the red shift and of the characteristic features of the OD-stretching bands in both IR and Raman spectra.