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     

Spectroscopic and theoretical studies on solid 1,2-ethylenediamine dihydrochloride salt.

A structural study of [H3N(CH2)2NH3)]2+.2Cl-, the smallest element of the homologous series of the alpha,omega-diamine dihydrochlorides, was carried out by means of Raman and FTIR spectroscopy coupled to ab initio molecular orbital (MO) calculations. As a primary concern, an adequate molecular model for the representation of these solid amine salts was chosen. Thus, several models, varying in the number and position of the counterions as well as in the number of diamine units, were considered. It was found that the best molecular system (i.e., that yielding the best compromise between accuracy and computational requirements) consists of one ethylenediamine cation surrounded by six chloride ions in an arrangement based on the crystal structure reported in the literature for [H3N(CH2)2NH3)]2+.2Cl-. This conclusion will hopefully allow for a better understanding of the conformational preferences, in the solid state, of these biologically relevant linear polyamines.     

Structures and vibrational spectra of the sulfur-rich oxides SnO (n = 4-9): the importance of pi*-pi* interactions

The structures of a large number of isomers of the sulfur oxides S(n)O with n = 4-9 have been calculated at the G3X(MP2) level of theory. In most cases, homocyclic molecules with exocyclic oxygen atoms in an axial position are the global minimum structures. Perfect agreement is obtained with experimentally determined structures of S(7)O and S(8)O. The most stable S(4)O isomer as well as some less stable isomers of S(5)O and S(6)O are characterized by a strong pi*-pi* interaction between S==O and S==S groups, which results in relatively long S--S bonds with internuclear distances of 244-262 pm. Heterocyclic isomers are less stable than the global minimum structures, and this energy difference approximately increases with the ring size: 17 (S(4)O), 40 (S(5)O), 32 (S(6)O), 28 (S(7)O), 45 (S(8)O), and 54 kJ mol(-1) (S(9)O). Owing to a favorable pi*-pi* interaction, preference for an axial (or endo) conformation is calculated for the global energy minima of S(7)O, S(8)O, and S(9)O. Vapor-phase decomposition of S(n)O molecules to SO(2) and S(8) is strongly exothermic, whereas the formation of S(2)O and S(8) is exothermic if n<7, but slightly endothermic for S(7)O, S(8)O, and S(9)O. The calculated vibrational spectra of the most stable isomers of S(6)O, S(7)O, and S(8)O are in excellent agreement with the observed data.     

Hydrogen‐bonded complexes of nicotine with simple alcohols

Ab initio and density functional theory studies have been performed on the hydrogen‐bonded complexes of neutral and protonated nicotine with ethanol, methanol, and trifluromethanol to explore their relative stability in a systematic way. Among all the hydrogen‐bonded nicotine complexes considered here, protonated forms in nicotine–ethanol and nicotine–methanol, and neutral form in nicotine–trifluromethanol complexes have been found to be the most stable. In the former two complexes, the proton attached to the pyrrolidine nitrogen acts as a strong hydrogen bond donor, whereas the pyrrolidine nitrogen atom acts as a hydrogen bond acceptor in the latter case. Neutral complex of nicotine with trifluromethanol has been found to possess a very short hydrogen bond (1.57 Å) and basis set superposition error corrected hydrogen bond energy value of 19 kcal/mol. The nature of the various hydrogen bonds formed has been investigated through topological aspects using Bader's atoms in molecules theory. From the calculated topological results, excellent linear correlation is shown to exist among the hydrogen bond length, electron density, and its Laplacian at the bond critical points for all the complexes considered. The natural bond orbital analysis has been carried out to investigate the charge transfer in the nicotine alcohol complexes. In contrast to the blue shifting behavior that is generally exhibited by other C? H···O hydrogen bonds involving sp³ carbon atom, the C? H···O hydrogen bond in the protonated nicotine–ethanol and methanol complexes has been found to be proper with red shifting in nature. © 2011 Wiley Periodicals, Inc.     

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.     

Conformational and tautomeric preferences in 3‐aminoacrylaldehyde: A theoretical study

The conformational study of 3‐aminoacrylaldehyde were performed at various theoretical levels and the equilibrium conformations were determined. Furthermore, to have more reliable energies, the total energies of all forms recomputed at G2MP2 and CBS‐QB3. Theoretical calculations clearly show that the intramolecular hydrogen bond (IHB) is the origin of conformational preference and the resultant IHB order at HF/6‐311++G(d,p), MP2/6‐311++G(d,p), and CBS‐QB3 levels is different from the order which obtained from the B3LYP/6‐311++G(d,p), G2MP2, the geometrical parameters, AIM, and NBO analyses. Furthermore, our theoretical results reveal that the ketoamine (KA) tautomeric group is more stable than the enolimine (EI) and ketoimine (KI) ones. The IHB and tautomeric process could not rationalize the irregular stability of KA group with respect the others. But the population analyses of the possible conformations by NBO predict that the π‐electron delocalization, especially unusual π → π charge transfer, is the origin of tautomeric preference. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Ab initio HF-SCF study of naphthazarin: Geometries,isomerism, hydrogen bonding,and vibrational spectrum

The structural properties and intramolecular hydrogen bonding of a series of structures of naphthazarin molecule were investigated by ab initio HF-SCF methods. The geometries of theC _2v ,C _2h ,D _2h , andC _s symmetry structures were optimized using split-valence basis sets. MP2/6-31G^*// HF/6-31G single-point energy calculations indicate that theC _2v isomer (5,8-dihydroxy-1,4-naphthoquinone) is the lowest energy structure of the molecule and that theC _2h symmetry one (4,8-dihydroxy-1,5-naphthoquinone), lying 37 kJ/mol above theC _2v form, is the other stable isomer of naphthazarin. At the HF/6-31G level, the intramolecular proton exchange between two equivalentC _2v structures is a two-step process where each proton can be independently transferred through an unsymmetrical potential having a 1,5-quinone intermediate, theC _2h symmetry structure, and two equivalent transition states ofC _s symmetry, with a barrier height equal to 38 kJ/ mol (MP2/6-31G^*//HF/6-31G). The study of naphthazarin molecule is flanked by a theoretical investigation on theC _2v andC _2h isomers of the parent naphthoquinone and dihydroxynaphthalene molecules. The SCF vibrational spectrum of the ground state of naphthazarin, harmonic frequencies, and infrared and Raman band intensities were computed at the HF/6-31G level. The results of the calculations are compared with the matrix isolation FT-IR spectroscopy measurements and with the infrared and Raman spectra of the crystal molecule.     

The competition between the intramolecular hydrogen bond and π‐electron delocalization in trifluoroacetylacetone—A theoretical study

Conformational study of trifluoroacetylacetone was carried out using the HF, B3LYP, and MP2 methods with the 6‐31G(d, p) and 6‐311++G(d, p) basis sets. All of the results show that the chelated enol structures (E11 and E31) have extra stability with respect to the other forms and one of them (E11) is global minimum. The energy gap between the chelated forms is in the range 0.7–5.9 kJ mol^?1. Theoretical calculations show that this compound has an asymmetric double minimum potential energy surface which is in contrast with the electron diffraction result. Moreover, the computational results predict that due to the withdrawing effect of CF₃ group, hydrogen bond in trifluoroacetylacetone is weaker than the acetylacetone. Because of the more stability of E11, it is expected that the hydrogen bond energy in E11 is greater than the E31, but at all of the computational levels with most extended basis set the converse results were observed. These results clearly show that the hydrogen bond is not a superior parameter in conformational preference and the contribution of resonance is probably greater than the hydrogen bond. Finally, the analysis of this system by quantum theory of atoms in molecules and natural bond orbital methods fairly support the ab initio results. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Vibrational dynamics of hydrogen‐bonded HCN complexes with OH and NH acids: Computational DFT systematic study

Vibrational properties (band position, infrared [IR], and Raman intensities) of C?N stretching mode were studied in 65 gas phase hydrogen‐bonded 1:1 complexes of HCN with OH acids and NH acids using density functional theory (DFT) calculations at the B3LYP‐6‐311++G(d,p) level. Furthermore, general characteristics of the hydrogen bonds and vibrational changes in acids OH/NH stretching bands were also considered. Experimentally observed blue shift of the C?N stretching band promoted by hydrogen bonding, which shortens the triple bond length, is very well reproduced and quantitatively depends on the hydrogen bond length. Both IR and Raman ν(C?N) band intensities are enhanced, also in good agreement with the experimental results. IR intensity increase is a direct function of the hydrogen bond energy. However, the predicted Raman intensity raise is a more complex function, depending simultaneously on characteristics of both the hydrogen bond (C?N bond length) and the H‐donating acid (polarizability). With these two parameters, ν (C?N) Raman intensities of the complexes are explained with a mean error of ±2.4%. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Effect of hydrogen‐bonded interactions on the energetics and spectral properties of the astromolecule aminoacetonitrile

A detailed and systematic electronic structure calculation has been performed to analyze the hydrogen‐bonded interaction of aminoacetonitrile (H₂NCH₂CN) with hydrogen cyanide (HCN) and Glycine (H₂NCH₂COOH). Both HCN and aminoacetonitrile have already been detected in the interstellar medium (ISM) and their active role in the molecular mechanisms of glycine production has already been recognized. Four different density functional models have been used to study the effect of hydrogen bond formation on the energetic stability and vibrational spectra of the aminoacetonitrile‐HCN and aminoacetonitrile‐glycine complexes in gas phase. The aminoacetonitrile‐glycine dimer is energetically far more stable than all forms of aminoacetonitrile‐HCN dimers. Elastic and inelastic scattering of light off the hydrogen‐bonded clusters have been investigated in details via Rayleigh and Raman spectroscopic parameters. The dipole moments and depolarization ratios are found to be sensitive on the type of hydrogen‐bond network. The mean polarizabilty show appreciable dependence on the choice of the DFT‐model. In general, all the chemical groups (OH, CN, NH₂, and CH) that participate directly in the hydrogen bond formation suffer appreciable variation in the intensity of vibration.     

Spectra and structure of organophosphorus compounds: XLIX. Microwave,infrared, and Raman spectra,electric dipole moment,molecular structure,and ab initio calculations of dimethylphosphonothioic fluoride

The microwave spectra of (CH₃)₂PSF, (CH₃)(CD₃)PSF, (CD₃)₂PSF, and (CH₃)₂P³⁴SF have been investigated from 20.0 to 40.0 GHz. Botha-type R branch andc-type Q branch transitions have been measured in the ground states of each isotopic species. From a least-square adjustment to fit 12 rotational constants, the following structural parameters were obtained:r(P–F)=1.582 ± 0.003 Å;r(P=S)=1.902 ± 0.001 Å;r(P-C)=1.800 ± 0.001 Å;r(C-H)=1.088 ± 0.002 Å; HCP=109.28 ± 0.12°; SPF=114.50 ± 0.13°; and SPC=116.33 ± 0.06°. From Stark effect measurements, the dipole moment components have been determined to be ¦ _a ¦ =3.556 ± 0.005; ¦ _c ¦=2.026 ± 0.009; and ¦ _t ¦=4.093 ± 0.009 (D). The Raman spectra (3200 to 100 cm^–1) of each isotopic species have been measured for the solid, and liquid and qualitative depolarization values obtained. Additionally, the mid-infrared spectra (3200 to 500 cm^–1) of the solids have been recorded. Proposed assignments of the normal modes have been made on the basis of Raman depolarization values and group frequencies which are supported by normal coordinate analysis utilizing an ab initio force field. Optimized structural parameters have been obtained with both the 3-21G^* and 6-31G^* basis sets. These results are compared to the corresponding quantities for several similar molecules.For part XLVIII, seeJ. Raman Spectrosc.1922,23, 107.     

Molecular structure and vibrational spectra of phenobaraitone by density functional theory and ab initio hartree-Fock calculations

Quantum chemistry calculations have been performed using Gaussian03 program to compute optimized geometry, harmonic vibrational frequency along with intensities in IR and Raman spectra at RHF/6-31++G** and B3LYP/6-31++G** levels for phenobarbitone (C₁₂H₁₂N₂O₃) in the ground state. The scaled harmonic vibrational frequencies have been compared with experimental FT-IR and FT-Raman spectra. Theoretical vibrational spectra of the title compound were interpreted by means of potential energy distributions (PEDs) using MOLVIB program. A detailed interpretation of the infrared spectra of the title compound is reported. On the basis of the agreement between the calculated and observed results, the assignments of fundamental vibrational modes of phenobarbitone were examined and some assignments were proposed. The theoretical spectrograms for FT-IR and FT-Raman spectra of the title compound have been constructed.     

A theoretical study on three conformational structures of 2,6‐distyrylpyridine

Ab initio methods at the levels HF/cc‐pVDZ, HF/6‐31G(d,p), MP2/cc‐pVDZ, and MP2/6‐31G(d,p), as well as methods based on density functional theory (DFT) employing the hybrid functional B3LYP with the basis sets cc‐pVDZ and 6‐31G(d,p), have been applied to study the conformers of 2,6‐distyrylpyridine. Bond distances, bond angles, and dihedral angles have been calculated at the B3LYP level. The calculated values were in good agreement with those measured by X‐ray diffraction analysis of 2,6‐distyrylpyridine. The values calculated using the Hartree‐Fock method and second‐order perturbation theory (MP2) were inconsistent. The optimized lowest‐energy geometries were calculated from the reported X‐ray structural data by the B3LYP/cc‐pVDZ method. Three conformations, A, B, and C, were proposed for 2,6‐distyrylpyridine. Calculations at the three levels of theory indicated that conformation A was the most stable structure, with conformations C and B being higher in energy by 1.10 and 2.57 kcal/mol, respectively, using the same method and basis function. The same trend in the relative energies of the three possible conformations was observed at the two levels of theory and with the different basis sets employed. The reported X‐ray data were utilized to optimize total molecular energy of conformation A at the different calculation levels. The bond lengths, bond angles, and dihedral angles were then obtained from the optimized geometries by ab initio methods and by applying DFT using the two basis functions cc‐pVDZ and 6‐31G(d,p). The values were analyzed and compared. The calculated total energies, the relative energies of the molecular orbitals, the gap between them, and the dipole moment for each conformational structure proposed for 2,6‐distyrylpyridine are also reported. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Ab initio study of the structure,cooperativity, and vibrational properties in the mixed hydrogen‐bonded trimers of hydrogen isocyanide and water

The five trimers of H₂O···HNC···H₂O, H₂O···H₂O···HNC, HNC···H₂O···H₂O, H₂O···HNC···HNC, and HNC···HNC···H₂O have been studied with quantum chemical calculations. Their structures, harmonic vibrational frequencies and interaction energies have been calculated at the B3LYP and MP2 levels with the aug‐cc‐pVDZ and aug‐cc‐pVTZ basis sets. The cooperative effect on these properties has also been studied quantitatively. For HNC:(H₂O)₂ systems, the cyclic H₂O···H₂O···HNC trimer is most stable with an interaction energy of ?16.01 kcal/mol and a large cooperative energy of ?3.25 kcal/mol at the MP2/aug‐cc‐pVTZ level. For H₂O:(HNC)₂ systems, the interaction energy and cooperative energy in the H₂O···HNC···HNC trimer are larger than those in the HNC···HNC···H₂O trimer. The NH stretch frequency has a blue shift for the terminal HNC molecule in the HNC···H₂O···H₂O and HNC···HNC···H₂O trimers and a red shift in other cases. A many‐body analysis has also been performed to understand the interaction energies in these hydrogen‐bonded clusters. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Rapid Prediction of the Hydrogen Bond Cooperativity in N‐methylacetamide Chains

A method is proposed to rapidly predict the hydrogen bond cooperativity in N‐methylacetamide chains. The parameters needed are obtained from the fittings to the hydrogen bonding energies in the formamide chains containing 2 to 8 monomeric units. The scheme is then used to calculate the individual hydrogen bonding energies in N‐methylacetamide chains containing 2 to 7 monomeric units. The cooperativity predicted is in good agreement with those obtained from MP2/6‐31+G** calculations by including the BSSE correction. Our scheme is further employed to predict the individual hydrogen bonding energies in larger N‐methylacetamide chains containing up to 200 monomeric N‐methylacetamide units, to which the MP2 method cannot be applied. Based on our scheme, a cooperative effect of over 170 % of the dimer hydrogen bonding energy in long N‐methylacetamide chains is predicted. The method is also applied to heterogeneous chains containing formamide, acetamide, N‐methylformamide, and N‐methylacetamide. The individual hydrogen bonding energies in these heterogeneous chains are also in good agreement with those obtained from MP2 calculations with the BSSE correction, further demonstrating that our method is reasonable.     

Ab initio studies on vibrational spectra of XSO_2NCO (X=F,Cl):harmonic force fields and frequency assignments

ＴｈｅｆｌｕｏｒｏｓｕｌｆｏｎｙｌｉｓｏｃｙａｎａｔｅａｎｄｃｈｌｏｒｏｓｕｌｆｏｎｙｌｉｓｏｃｙａｎａｔｅｍｏｌｅｃｕｌｅｓａｒｅｔｈｅｍｏｓｔｒｅａｃｔｉｖｅｍｏｌｅｃｕｌｅｓｆｏｒｔｈｅｔｒａｎｓｆｅｒｅｎｃｅｏｆｔｈｅＮＣＯｇｒｏｕｐａｎｄａｒｅｖｅｒｙｉｍｐｏｒｔａｎｔｉｎｓｙｎｔｈｅｔｉｃｃｈｅｍｉｓｔｒｙ.Ｄｕｅｔｏｔｈｅｉｒｖｅｒｓａｔｉｌｉｔｙｉｎｃｈｅｍｉｃａｌｒｅａｃｔｉｏｎｓ,ｔｈｅｙｈａｖｅａｌｓｏｂｅｃｏｍｅｔｈｅｓｕｂｊｅｃｔｏｆｍａｎｙｓｔｒｕｃｔｕｒａｌｓｔｕｄｉ…     

Computational study of structures and proton transfer in hydrogen‐bonded ammonia complexes using semiempirical valence‐bond approach

In this study, the seGVB method was implemented for the N H bonding system, specifically for hydrogen‐bonded ammonia complexes, and the model well reproduces the MP2 geometries and energetics. A comparison between the ammonia dimer and water dimer is given from the viewpoint of valance‐bond structures in terms of the calculated bond energies and pair–pair interactions. The linear hydrogen bond is found to be stronger than the bent bonds in both cases, with the difference in energy between the linear and cyclic structures being comparable in both cases although the NH bonds are generally weaker. The energy decomposition clearly demonstrates that the changes in electronic energy are quite different in the two cases due to the presence of an additional lone pair on the water molecule, and it is this effect which leads to the net stabilization of the cyclic structure for the ammonia dimer. Proton‐transfer profiles for hydrogen‐bonded ammonia complexes [NH₂ H NH₂]⁻ and [NH₃ H NH₃]⁺ were calculated. The barrier for proton transfer in [NH₃ H NH₃]⁺ is larger than that in [NH₂ H NH₂]⁻, but smaller than that in the protonated water dimer. The different bonding structures substantially affect the barrier to proton transfer, even though they are isoelectronic systems. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 73: 357–367, 1999     

The vibrational spectra of the boron halides and their molecular complexes. Part 13. Ab initio studies of the complexes of boron trifluoride with formaldehyde and some of its analogs

The structure, the interaction energy, and the vibrational spectrum of the electron donor–acceptor complex formed between boron trifluoride, as a Lewis acid, and formaldehyde, as a Lewis base, have been determined by means of ab initio calculations at the second‐order level of Møller–Plesset perturbation theory, using a triple‐zeta basis set with polarization and diffuse functions on all atoms. The object was to examine the differences in the properties of the complexes formed between boron trifluoride and an oxygen base containing oxygen in the sp² hybrid state with those of some sp³ oxygen bases studied earlier. The investigation has then been extended to include the related bases thioformaldehyde and methanimine, to assess the effect of substituting the oxygen atom by a sulphur atom or an imino group. A further range of formaldehyde analogs, containing one and two methyl groups, one and two fluorine atoms, and one methyl group and one fluorine atom, has also been included. The preferred structure of each complex has been identified. The computed data have been compared, and the complex properties correlated with some relevant physical properties of the bases. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Comparative study on formamide–water complex

The hydrogen bonding of 1:1 complexes formed between formamide and water molecule have been investigated systematically using Hartree–Fock (HF), hybrid density functional theory (B3LYP), and post‐Hartree–Fock (MP2 and CCSD(T)) methods with range of basis sets 6‐31G(d), cc‐pVXZ (X = D, T, Q) and aug‐cc‐pVYZ (Y = D, T). Three stable structures are considered on the potential energy surface of formamide and water system. The optimized geometric parameters and interaction energies for various isomers at different levels are estimated. The IR frequencies, intensities, and frequency shifts are reported. This study shows that B3LYP/aug‐cc‐pVDZ method gives better performance for formamide‐water complexes. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010.     

A combined computational and experimental approach for investigating a hydrogen‐bonded supermolecular compound comprising benzimidazole and malonic acid

A hydrogen‐bonded tetramer supramolecular motif of 2:2 benzimidazole (BIZ) and malonic acid (MLA) has been synthesized and characterized by elemental analysis, infrared (IR), and X‐ray single crystal diffraction. Thermal stability analyses demonstrate that this supramolecular adduct is a new material and it is not the ordinary superposition of the original monomers. Density function theory (DFT) calculations for the models of dimers, trimers, and tetramer comprising BIZ and MLA have been carried out at B3LYP/6‐31G* and PBE1PBE/6‐31G* levels of theory, respectively. By comparing the calculated results with the experiments (single crystal structure, IR spectra, and thermal analysis) and based on the statistic thermodymnamic calculations, it is concluded that the dimers cannot exist at room temperature and the tetramer can simulate the title supramolecular complex better than the two trimers. Further studies on the model of tetramer indicate that the hydrogen bond of N···H? O is stronger than that of O···H? N. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008