A study of the conformational stability and the vibrational spectra of 2,3-dichloro-1-propanol

The conformational stability and the three rotor internal rotations in 2,3-dichloro-1-propanol were investigated at DFT-B3LYP/6-311 + G**, MP2/6-311 + G** and MP4(SDQ) levels of theory. From the calculated potential energy surface, ten distinct minima were located all of which were predicted to have real frequencies at the B3LYP level of theory. The calculated lowest energy minima in the potential curves of the molecule were predicted to correspond to the Ggg and Gtg1 structures. The observed broad and very intense infrared band centered at about 3370 cm^?1 supports the existence of the strong intermolecular H-bonding in 2,3-dichloro-1-propanol. The equilibrium constants for the conformational interconversion in the molecule were estimated from the calculated Gibb's energies at the B3LYP/6-311 + G** level of calculation and found to correspond to an equilibrium mixture of about 49% Ggg, 27% Gtg1, 5% Ggt and 5% Tgg conformations at 298.15 K.     

Reaction of CH3CHO with Y+: A density functional theoretical study

The reaction mechanism of the Y⁺ cation with CH₃CHO has been investigated with a DFT approach. All the stationary points are determined at the UB3LYP/ECP/6-311++G** level of the theory. Both ground and excited state potential energy surfaces are investigated in detail. The present results show that the title reaction start with the formation of a CH₃CHO-metal complex followed by C-C, aldehyde C-H, methyl C-H and C-O activation. These reactions can lead to four different products (Y⁺CH₄ + CO, Y⁺CO + CH₄, Y⁺COCH₂ + H₂ and Y⁺O + C₂H₄). The minimum energy reaction path is found to involve the spin inversion in the different reaction steps, this potential energy curve-crossing dramatically affects reaction exothermic. The present results may be helpful in understanding the mechanism of the title reaction and further experimental investigation of the reaction.     

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.     

Potential Energy Surface of Cooperative Proton Tunneling in the C2H+ 3 Cation

The nine-dimensional potential energy surface for proton tunneling in the nonrigid C₂H⁺ ₃ cation was constructed from quantum-chemical data [MP4SDQ(T)/6-311++G(3df,3pd)] on the equilibrium geometry, energy, frequencies, and eigenvectors of the normal vibrations at the stationary points and transition states using the theory of isodynamic symmetry groups along the tunneling path.     

Conformers and intramolecular hydrogen bonding of the oxalic acid monomer and its anions

Density functional theory, B3LYP/6‐31G** and B3LYP/6‐311+G(2d,p), and ab initio MP2/6‐31G** calculations have been carried out to investigate the conformers, transition states, and energy barriers of the conformational processes of oxalic acid and its anions. QCISD/6‐31G** geometrical optimization is also performed in the stable forms. Its calculated energy differences between the two most stable conformers are very near to the related observed value at 7.0 kJ/mol. It is found that the structures and relative energies of oxalic acid conformers predicted by these methods show similar results, and that the conformer L1 (C_2h) with the double‐interfunctional‐groups hydrogen bonds is the most stable conformer. The magnitude of hydrogen bond energies depends on the energy differences of various optimized structures. The hydrogen bond energies will be about 32 kJ/mol for interfunctional groups, 17 kJ/mol for weak interfunctional groups, 24 kJ/mol for intra‐COOH in (COOH)₂, and 60 kJ/mol for interfunctional groups in (COOH)COO⁻¹ ion if calculated using the B3LYP/6‐311+G(2d,p) method. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 76: 541–551, 2000     

Is tetrahedral H42+ a minimum? Anomalous behavior of popular basis sets with the standard p exponents on hydrogen

The nature of the tetrahedral H₄²⁺ stationary point (minimum or triply degenerate saddle) depends remarkably upon the theoretical level employed. Harmonic vibrational analyses with, e.g., the 6-31G** (and 6-31 + +G**) and Dunning's [4s2p1d;2s1p] [D95(d,p)] basis sets using the standard p exponent suggest (erroneously) that the T_d geometry is a minimum at both the HF and MP2 levels. This is not the case at definitive higher levels. The C₃H₄²⁺ structure with an apical H is another example of the failure of the calculations with the 6-31G**, 6-311G**, and D95(d,p) basis sets. Even at MP2/6-31G** and MP2/ cc-pVDZ levels, the C_3v structure has no negative eigenvalues of the Hessian. Actually, this form is a second-order saddle point as shown by the MP2/6-31G** calculation with the optimized exponent. The D_4h methane dication structure is also an example of the misleading performance of the 6-31G** basis set. In all these cases, energy-optimized hydrogen p exponents give the correct results, i.e., those found with more extended treatments. Optimized values of the hydrogen polarization function exponents eliminate these defects in 6-31G** calculations. Species with higher coordinate hydrogens may also be calculated reliably by using more than one set of p functions on hydrogen [e.g., the 6-31G(d,2p) basis set]. Not all cases are critical. A survey of examples, also including some boron compounds, provides calibration. © 1993 John Wiley & Sons, Inc.     

Theoretical study of the relative stabilities of H+(X)2 and H+(X)3 conformers and their clustering energies: X  CO and N2

Third-order Møller–Plesset perturbation theory (MP 3) with a 6-31G** basis set was applied to study the relative stabilities of H⁺(X)₂ conformations (X ? CO and N₂) and their clustering energies. The effect of both basis set extensions and electron correlation is not negligible on the relative stabilities of the H⁺(CO)₂ clusters. The most stable conformation of H⁺(CO)₂ is found to be a C_∞v structure in which a carbon atom of CO bonds to the proton of H⁺(CO), whereas that of H⁺(N₂)₂ is a symmetry D_∞h structure. The second lowest energy conformations of H⁺(CO)₂ and H⁺(N₂)₂ lie within 2 kcal/mol above the energies of the most stable structures. Clustering energies computed using MP 3 method with the 6-31G** basis set are in good agreement with the experimental findings of Hiraoka, Saluja, and Kebarle. The low-lying singlet conformations of H⁺(X)₃ (X ? CO and N₂) have been studied by the use of the Hartree–Fock MO method with the 6-31G** basis set and second-order Møller–Plesset perturbation theory with a 4-31G basis set. The most stable structure is a T-shaped structure in which a carbon atom of CO (or a nitrogen atom of N₂) attacks the proton of the most stable conformation of H⁺(X)₂ clusters.     

Theoretical studies on the reaction of pentafulvenone with water in gas phase and aqueous solvent

In gas phase, the hydrations of pentafulvenone to generate three types of cyclopentadienyl carboxylic acids are studied theoretically at the MP2/6-311+G**//B3LYP/6-311+G** level. A water molecule attacking the C=O double bond of pentafulvenone can yield cyclopentadienyl carboxylic acids via the formation of fulvenediols, and attacking the C=C double bond of pentafulvenone can directly yield cyclopentadienyl carboxylic acid. The barriers of rate-determining transition states are 42.2 and 30.4 kcal mol⁻¹, respectively. The barriers of rate-determining transition states for two water molecules system are 20.2 and 19.6 kcal mol⁻¹, respectively. The products can isomerize to each other. In aqueous solvent, the hydrations of pentafulvenone are investigated using PCM-UAHF model at the MP2 (PCM)/6-311+G**// B3LYP (PCM)/6-311+G** and MP2 (PCM)/6-311+G**// B3LYP/6-311+G** levels. The barriers of all rate-determining transition states are decreased. The added water molecule acts as catalyst in both gas phase and aqueous solvent. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Ab initio study of the formation of C3H3+ from the reaction of CH3+ with acetylene

Ab initio molecular orbital theory has been used to study the mechanism of the formation of C₃H₃⁺ from the reaction of CH₃⁺ with acetylene. The highest level geometry optimizations and frequencies were computed at MP2-FC/6-31G**; single point energies of all the critical structures were computed to the MP4-FC/6-31G**//MP2-FC/6-31G** theory level. One of the three alternative transition structures leading to the formation of C₃H₃⁺ gives the cyclopropenyl cation and the other two the propargyl cation. The proportions of C₃H₂D⁺ and C₃HD₂⁺ obtained when CD₃⁺ reacts with acetylene, and the composite nature of the metastable peak observed for the [C₃H₅]⁺→[C₃H₃]⁺ + H₂ fragmentation are explained by assuming a different degree of deuterium scrambling depending on the energy of the system. © 1996 by John Wiley & Sons, Inc.     

The interaction model between metal cation and tropylium: a quantum chemistry predication

The aromatic cation tropylium, C₇H₇⁺, predicted at the MP2/6-31G** level, is capable of binding with metal cations Be²⁺ or Mg²⁺, forming M²⁺–C₇H₇⁺ complexes. The obstacle for their binding is almost electrostatic repulsion, and the binding is from polarization and charge transfer. The orbital interaction between the M²⁺ and C₇H₇⁺ is mainly the s–π and p–π interactions. Interestingly, Be²⁺ is possible to pass through the ring of C₇H₇⁺, while Mg²⁺ is not. The intrinsic IR band of the M²⁺–C₇H₇⁺ complex is below 600 cm⁻¹, which results from the vibration of the M²⁺ along the normal axis of C₇H₇⁺.     

On the additivity of basis set effects in some simple fluorine containing systems

The reactions F + H₂ → HF + H, HF → H + F, F → F⁺ + e^? and F + e^? → F^? were used as simple test cases to assess the additivity of basis set effects on reaction energetics computed at the MP4 level. The 6-31G and 6-311G basis sets were augmented with 1, 2, and 3 sets of polarization functions, higher angular momentum polarization functions, and diffuse functions (27 basis sets from 6-31Gd, p) to 6-31 ++ G(3df, 3pd) and likewise for the 6-311G series). For both series substantial nonadditivity was found between diffuse functions on the heavy atom and multiple polarization functions (e.g., 6-31 + G(3d, 3p) vs. 6-31 + G(d, p) and 6-31G(3d, 3p)). For the 6-311G series there is an extra nonadditivity between d functions on hydrogen and multiple polarization functions. Provided that these interactions are taken into account, the remaining basis set effects are additive to within ±0.5 kcal/mol for the reactions considered. Large basis set MP4 calculations can also be estimated to within ±0.5 kcal/mol using MP2 calculations, est. E_MP4(6-31 ++ G(3df, 3pd)) ≈ E_MP4(6-31G(d, p)) + E_MP2(6-31 ++ G(3df, 3pd)) – E_MP2(6-31G(d, p)) or E_MP4(6-31 + G(d, p) + E_MP2(6-31 ++ G(3df, 3pd)) – E_MP2(6-31 + G(d, p)) and likewise for the 6-311G series.     

Counterion effects on the isotopic properties of the Cl<Superscript>−</Superscript> and Li<Superscript>+</Superscript> ions in aqueous solutions

The reduced partition function ratios between isotopic forms (β-factors) were calculated by the ab initio RHF/6-311++G**(3df, 3p) and MP2/6-311++G**(3df, 3p) quantum-chemical methods for hydrated chloride ion and ion pair hydrates NaCl·nH₂O and LiOH · nH₂O. The influence of the Na⁺ cation on the β-factor value and the chlorine isotope separation factor in the precipitation of NaCl from concentrated aqueous solutions was found to be substantial. At the same time, the presence of OH^? counterions had no noticeable effect on the β-factor of the hydrated Li⁺ cation.     

Study of the temperature‐dependent conformational averaging of 1H NMR resonances in vinylcyclopropane through the use of ab initio methodology and Boltzmann statistics

The temperature dependence of the ¹ H NMR resonance of the C‐4 olefinic proton in vinylcyclopropane was investigated through a combination of ab initio calculations and Boltzmann statistics. A torsional energy profile as a function of the 〈?〉 dihedral angle was obtained using HF methodology with a 6–311G** basis set, while the corresponding ¹ H chemical shift profiles for the C‐4 proton were computed using the GIAO approach and either HF, DFT (B3LYP) or MP2 methods at the 6–311G** level of theory. Chemical shifts at different temperatures calculated as canonical ensemble averages in which the different ab initio ¹ H chemical shift profiles and a Boltzmann factor defined by the HF/ 6–311G** energy function are employed reproduce remarkably well the temperature dependence observed experimentally. Attempts to perform a similar study using only the GIAO‐MP2 ¹ H chemical shift profile and 〈?〉 dihedral angle trajectories obtained from molecular dynamics simulations at different temperatures failed to reproduce the experimental trends. This shortcoming was attributed to the inability of the force fields employed, Tripos 6.0 and MMFF94, to reproduce properly the three‐well torsional potential of vinylcyclopropane. The application of both methodologies to the calculation of population‐dependent chemical shifts in other systems is discussed. Copyright © 2002 John Wiley & Sons, Ltd.     

FTIR and quantum chemical study of LiBr solvation in acetonitrile solutions

FTIR spectroscopy and quantum chemical calculations at the RTF + MP2/6-311G** level of theory with solvation model density (SMD) corrections were used to study ion solvation and association in LiBr/acetonitrile solutions. The aim of this study was to establish the composition and geometry of the predominant ionic species solvated by acetonitrile molecules and to analyse their spectroscopic signatures. The results obtained make it possible to propose an equilibrium between Li⁺Br⁻(CH₃CN)₃, Li⁺(CH₃CN)₄, and anionic Br⁻(CH₃CN)_n complexes with an undetermined n value and bent coordination of the solvent molecules. The calculated wavenumbers and the geometric parameters of the solvated ionic species were found to be in excellent agreement with the experimental data.     

Ab initio investigation on stability and properties of XYCO... HZ complexes. II: Post hartree-fock studies on H2CO... HF

Ab initio MP2 level of theory in conjunction with three basis sets of a triple-zeta quality was applied to study the molecular geometry and stability of the H₂CO... HF complex. An interaction energy predicted for this system at the highest, MP4(SDTQ)/6-311 + +G(2df, 2pd)//MP2/6-311 + +G(2df, 2pd), level corrected for the BSSE and ZPE contributions amounts to -4.85 kcal/ mol. BSSE contributes significantly to the interaction energies at all applied levels. Reliable MP2/ 6-311 + +G(2df, 2pd) level harmonic vibrational frequencies, IR intensities, and the predicted isotopic shifts upon deuteration and¹⁸O substitution are presented in order to facilitate experimental studies on the IR spectrum of the title complex.     

Conformational study of the structure of 12-crown-4-alkali metal cation complexes

A conformational search was performed for the 12-crown-4 (12c4)-alkali metal cation complexes using two different methods, one of them is the CONFLEX method, whereby eight conformations were predicted. Computations were performed for the eight predicted conformations at the HF/6-31+G*, MP2/6-31+G*//HF/6-31+G*, B3LYP/6-31+G*, MP2/6-31+G*//B3LYP/6-31+G*, and MP2/6-31+G* levels. The calculated energies predict a C4 conformation for the 12c4-Na+, -K+, -Rb+, and -Cs+ complexes and a C(s) conformation for the 12c4-Li+ complex to be the lowest energy conformations. For most of the conformations considered, the relative energies, with respect to the C4 conformation, at the MP2/6-31+G*//B3LYP/6-31+G* are overestimated, compared to those at the MP2/6-31+G* level, the highest level of theory considerd in this report, by 0.2 kcal/mol. Larger relative energy differences are attributed to larger differences between the B3LYP and MP2 optimized geomtries. Binding enthalpies (BEs) were calculated at the above-mentioned levels for the eight conformations. The agreement between the calculated and experimental BEs is discussed.     

Does CH prefer a C2v rather than a Cs structure?

The closely related C_s ( 1 ) and C_2v ( 3 ) structures of CH have been reinvestigated at many ab initio levels using MP2/6-31G** and MP2/6-311 + + G(2df, 2pd) geometries. The largest basis sets employed were 6-311G(3df, 2p), 6-311 + + G(3df, 3pd), and the Dunning “correlation consistent” polarized triple-split valence basis set (cc-pVTZ). Electron correlation was probed at the MP4 level, but the QCISD method was also used with the largest basis sets. While electron correlation favors 3 over 1 by about 2 kcal/mol, the correlated relative energies with all basis sets employed range from 0.36–1.03 kcal/mol in favor of 1 . The best estimate of this difference, 0.86 kcal/mol, is essentially identical with the (scaled) zero-point energy difference, 0.84 kcal/mol, favoring 3 over 1 . These results indicate that 1 and 3 have almost exactly the same energy at 0 K. Our best value for the dissociation energy of CH is 42.0 kcal/mol [QCISD(T)/6-311 + + G(3df, 3pd)//MP2(fu)/6-311 + + G(2df, 2pd), corrected to 298 K], which agrees very well with the experimental value. © 1992 by John Wiley & Sons, Inc.     

Synthesis,characterization, and quantum chemical calculational studies on 3-<Emphasis Type="Italic">p</Emphasis>-methylphenyl-4-amino- 1, 2, 4-triazole-5-thione

The title compound of 3-p-methylphenyl-4-amino-1, 2, 4-triazole-5-thione was synthesized and characterized by elemental analysis, IR, electronic spectra, and X-ray single crystal diffraction. Quantum chemical calculations of the structure, natural bond orbital, and thermodynamic functions of the title compound were performed by using B3LYP/6-311G** and HF-6-311G** methods. Both the methods can well simulate the molecular structure. Vibrational frequencies were predicted, assigned and compared with the experimental values, and B3LYP/6-311G** method is superior to HF/6-311G** method to predict the vibrational frequencies. Electronic absorption spectra calculated by B3LYP/6-311G** method have some red shifts compared with the experimental ones and natural bond orbitals analyses indicate that the two absorption bands are mainly derived from the contribution of n → π^* and π → π^* transitions. On the basis of vibrational analyses, the thermodynamic properties of the title compound at different temperatures have been calculated, revealing the correlations between C ⁰ _p,m , S ⁰ _m , H ⁰ _m , and temperatures.     

Quantum-Chemical Study of Alkyl Carbenium Ions in 100% Sulfuric Acid

A quantum-chemical study of neutral and protonated monoalkyl sulfates RHSO₄and [RH₂SO₄]⁺(where R = CH₃, C₂H₅, iso-C₃H₇, and tert-C₄H₉) is carried out. Calculations are performed using the Hartree–Fock method in the 6-31G** and 6-31++G** basis sets taking into account electron correlation according to the Müller–Plesset perturbation theory MP2/6-31+G*//6-31+G*. Protonated tert-butyl sulfate was also calculated by the DFT B3LYP/6-31++G** method. It was found that monoalkyl sulfates are covalent compounds, and the complete abstraction of alkyl carbenium ions from them has substantial energy cost: 196.4, 161.7, 150.8 and 136.0 kcal/mol, respectively. Protonated methyl and ethyl sulfates are also covalent compounds according to the calculation. They have lower but still high energies of heterolytic dissociation (65.0 and 33.5 kcal/mol, respectively). The energy of R⁺abstraction from protonated isopropyl sulfate is much lower: 23.6 kcal/mol. The main covalent state and the ion–molecular pair, which is a carbenium ion [C(CH₃)₂H]⁺solvated by the H₂SO₄molecule, have about the same energy. The ground state of protonated tert-butyl sulfate corresponds to the ion–molecular complex [C(CH₃)⁺ ₃H₂SO₄] with still lower energy of carbenium ion [C(CH₃)₃]⁺abstraction, which is equal to 10.0 kcal/mol. Calculation according to the DFT B3LYP/6-31++G** method shows the absence of a minimum for the protonated tert-butyl sulfate with a covalent structure on the potential energy surface.     

HCN exchange on [Cu(HCN)4]+: a quantum chemical investigation

Density functional (B3LYP, B3PW91, X3LYP, BP86, PBEPBE, PW91PW91, and M06) and ab initio (MP2, MP4sdq, CCSD, and CCSD(T)) calculations with extended basis sets (6-311+G**, TZVP, LANL2DZ+p, and SDD+p, the latter including extra polarization and diffuse functions) indicate that HCN exchange on [Cu(HCN)₄]⁺ proceeds via an associative interchange (I_a) mechanism and a D_3h transition structure {[Cu(HCN)₅]⁺}^?. The activation barrier, relative to the model complex [Cu(HCN)₄]⁺·HCN, varies modestly, depending on the computational level. Typical values are 8.0?kcal?M^?1 (B3LYP/6-311+G**), 6.0?kcal?M^?1 (M06/6-311+G**), and 4.8?kcal?M^?1 (CCSD(T)/6-311+G**//MP2(full)/6-311+G**). Inclusion of an implicit solvent model (B3LYP(CPCM)/6-311+G**) leads to an activation barrier of 5.8?kcal?mol^?1. Comparison of the HCN exchange mechanisms on [Li(HCN)₄]⁺ (limiting associative, A) and [Cu(HCN)₄]⁺ (associative interchange, I_a) reveals that π back donation in the equatorial Cu–N bonds in the transition state determines the mechanism.