A comparative study of the energetics, structures, and mechanisms of the HCN H HNC and LiCN <-> LiNC isomerizations

The potential energy surfaces of the HCN<->HNC and LiCN<->LiNC isomerization processes were determined by ab initio theory using fully optimized triple-zeta double polarization types of basis sets. Both the MP2 corrections and the QCISD level of calculations were performed to correct for the electron correlation. Results show that electron correlation has a considerable influence on the energetics and structures. Analysis of the intramolecular bond rearrangement processes reveals that, in both cases, H (or Li+) migrates in an almost elliptic path in the plane of the molecule. In HCN<->HNC, the migrating hydrogen interacts with the in-plane pi,pi* orbitals of CN, leading to a decrease in the C-N bond order. In LiCN<->LiNC, Li+ does not interact with the corresponding pi,pi* orbitals of CN.

The hydrogen bond: a molecular beam microwave spectroscopist's view with a universal appeal

In this manuscript, we propose a criterion for a weakly bound complex formed in a supersonic beam to be characterized as a `hydrogen bonded complex'. For a `hydrogen bonded complex', the zero point energy along any large amplitude vibrational coordinate that destroys the orientational preference for the hydrogen bond should be significantly below the barrier along that coordinate so that there is at least one bound level. These are vibrational modes that do not lead to the breakdown of the complex as a whole. If the zero point level is higher than the barrier, the `hydrogen bond' would not be able to stabilize the orientation which favors it and it is no longer sensible to characterize a complex as hydrogen bonded. Four complexes, Ar-2-H2O, Ar-2-H2S, C2H4-H2O and C2H4-H2S, were chosen for investigations. Zero point energies and barriers for large amplitude motions were calculated at a reasonable level of calculation, MP2(full)/aug-cc-pVTZ, for all these complexes. Atoms in molecules (AIM) theoretical analyses of these complexes were carried out as well. All these complexes would be considered hydrogen bonded according to the AIM theoretical criteria suggested by Koch and Popelier for C-H center dot center dot center dot O hydrogen bonds (U. Koch and P. L. A. Popelier, J. Phys. Chem., 1995, 99, 9747), which has been widely and, at times, incorrectly used for all types of contacts involving H. It is shown that, according to the criterion proposed here, the Ar-2-H2O/H2S complexes are not hydrogen bonded even at zero kelvin and C2H4-H2O/H2S complexes are. This analysis can naturally be extended to all temperatures. It can explain the recent experimental observations on crystal structures of H2S at various conditions and the crossed beam scattering studies on rare gases with H2O and H2S.

Design, Synthesis, and DNA Binding Properties of Photoisomerizable Azobenzene−Distamycin Conjugates: An Experimental and Computational Study

Here, we present the synthesis, photochemical, and DNA binding properties of three photoisomerizable azobenzene−distamycin conjugates in which two distamycin units were linked via electron-rich alkoxy or electron-withdrawing carboxamido moieties with the azobenzene core. Like parent distamycin A, these molecules also demonstrated AT-specific DNA binding. Duplex DNA binding abilities of these conjugates were found to depend upon the nature and length of the spacer, the location of protonatable residues, and the isomeric state of the conjugate. The changes in the duplex DNA binding efficiency of the individual conjugates in the dark and with their respective photoirradiated forms were examined by circular dichroism, thermal denaturation of DNA, and Hoechst displacement assay with poly[d(A-T).d(T-A)] DNA in 150 mM NaCl buffer. Computational structural analyses of the uncomplexed ligands using ab initio HF and MP2 theory and molecular docking studies involving the conjugates with duplex d[(GC(AT)10CG)]2 DNA were performed to rationalize the nature of binding of these conjugates.

Ab initio studies of the structure and vibrational spectra of methyl isocyanide-borane and acetonitrile-borane

The SCF/DZP and MP2/DZP methods of ab initio quantum chemistry have been utilized to study the structure, vibrational spectra, binding energy, and barrier to internal rotation of methyl isocyanide-borane and acetonitrile-borane adducts. The eclipsed conformation of the complexes was predicted to be a minimum, and the staggered form is a transition state with a barrier height of about 10 cal/mol. The vibrational analyses of CH3NC-BH3 and CH3CN-BH3 and several of their isotopomers have been carried out by the GF matrix method. Computations have also been carried out for free CH3NC and CH3CN in order to investigate the changes in CH3NC and CH3CN as a result of their complex formation with BH3. To obtain an acceptable set of force constants, a recently proposed procedure ''RECOVES'' has been utilized. The increase in the N=C/C=N stretching force constant of CH3NC/CH3CN on adduct formation is interpreted with the help of Parr and Borkman's model. The binding energies for the two adducts have been determined taking basis set superposition error (BSSE) into consideration. The effect of the BSSE on structure, dipole moment, and vibrational frequencies of CH3CN and CH3NC is also evaluated. The predicted infrared band intensities for the two complexes are in good agreement with the experimentally observed features, and they have been utilized in the assignment of vibrational frequencies.

Electron donor-acceptor complexes of I-2 with diethyl ether and diethyl sulphide. An ab initio MO study

The charge-transfer complexes of I-2 with the n-donors diethyl ether and diethyl sulfide were studied at the Hartree-Fock and MP2 levels. The structures were fully optimized using the 3-21G((*)) basis set as well as with effective core potentials. The calculations consistently yield a C-2v structure for the ether-I-2 complex, but an unsymmetrical form for the sulfide-I-2 complex. A natural bond orbital analysis and the BSSE-corrected complexation energies reveal stronger interactions in the sulfide complex. The computed orbital energies of the monomers and complexes reproduce the trends in experimentally observed vertical ionization potentials.

Electron donor-acceptor complex of ICl with diethyl ether - He I photoelectron spectroscopy and ab initio molecular orbital study

The He I photoelectron spectrum of the diethyl ether-ICl complex has been obtained. The oxygen orbitals are shifted to higher binding energies and that of ICl to lower binding energies owing to complex formation. Ab initio molecular orbital (MO) calculations of the complex molecule showed that the bonding is between the sigma-type lone pair of oxygen and the I atom and that the complex has C-2v symmetry. The binding energy of the complex is computed to be 8.06 kcal mol(-1) at the MP2/3-21G* level. The orbital energies obtained from the photoelectron spectra of the complex are compared and assigned with orbital energies obtained by MO calculations. Natural bond orbital analysis (NBO) shows that charge transfer is from the sigma-type oxygen lone pair to the iodine atom and the magnitude of charge transfer is 0.0744 e.

SCF and electron correlation studies on structure, force constants and vibrational spectra of borane monoammoniate complex

Effects of basis set and electron correlation on the equilibrium geometry, force constants and vibrational spectra of BH3NH3 have been studied. A series of basis sets ranging from double zeta to triple zeta including polarization and diffuse functions have been utilized. All the SCF based calculations overestimate the dative B-N bond distance and considerable improvement occurs when the treatment for electron correlation is introduced. Detailed vibrational analysis for BH3NH3 has been carried out. The mean absolute percentage deviation of the ab initio predicted vibration frequencies of (BH3NH3)-B-11 from the experiment is about 10% for the SCF based calculations and the MP2 method shows better agreement, the overall deviation being 5-6%. The ground state effective force constants of BH3NH3 were obtained using RECOVES procedure. The RECOVES sets of force constants are found to be highly satisfactory for the prediction of the vibrational frequencies of different isotopomers of BH3NH3. The mean absolute percentage deviation of the calculated frequencies of different isotopomers from the experiment is much less than 1%. The RECOVES-MP2/augDZP set of force constants was found to be the best set among the different sets for this molecule. Theoretical infrared intensities are in fair agreement with the observed spectral features.

Structure and stability of Li2F

The equilibrium geometries and fundamental vibration frequencies of the Li2F system were calculated by ab initio methods at the MP2 = full/6-311(+ +)G** and CCSD(T) levels. Two isomers were observed and are best described as salts of the Li-2(+) cation with F-. A linear isomer with an arrangement of atoms such as Li-Li-F and a bent C-2v structure are predicted. The stability of these structures are discussed in terms of charge resonance between Li and Li+. (C) 1999 Elsevier Science B.V. All rights reserved.

Unimolecular HCl elimination from 1,2-dichloroethane: A single pulse shock tube and ab initio study

Thermal decomposition of 1,2-dichloroethane (1,2-DCE) has been studied in the temperature range of 10501175 K behind reflected shock waves in a single pulse shock tube. The unimolecular elimination of HCl is found to be the major channel through which 1,2-DCE decomposes under these conditions. The rate constant for the unimolecular elimination of HCl from 1,2-dichloroethane is found to be 10(13.98+/-0.80) exp(-57.8+/-2.0/RT) s(-1), where the activation energy is given in kcal mol(-1) and is very close to that value for CH3CH2Cl (EC). Ab initio (HF and MP2) and DFT calculations have been carried out to find the activation barrier and the structure of the transition state for this reaction channel from both EC and 1,2-DCE. The preexponential factors calculated at various levels of theory (BF/6-311++G**, MP2/6-311++G**, and B3LYP/6-311++G**) are (approximate to10(15) s(-1)) significantly larger than the experimental results. If the torsional mode in the ground state is treated as free internal rotation the preexponential factors reduce significantly, giving excellent agreement with experimental values. The DFT results are in excellent (fortuitous?) agreement with the experimental value for activation energy for 1,2-DCE while the MP2 and HF results seem to overestimate the barrier. However, DFT results for EC is 4.5 kcal mol(-1) less than the previously reported experimental values. At all levels, theory predicts an increase in HCI elimination barrier on beta-Cl substitution on EC.

Ab initio, DFT and transition state theory calculations on 1,2-HF, HCI and CIF elimination reactions from CH2F-CH2Cl

This paper reports ab intio, DFT and transition state theory (TST) calculations on HF, HCI and CIF elimination reactions from CH2Cl-CH2F molecule. Both the ground state and the transition state for HX elimination reactions have been optimized at HF, MP2 and DFT calculations with 6-31G*, 6-31G** and 6-311++G** basis sets. In addition, CCSD(T) single point calculations were carried out with MP2/6-311++G** optimized geometry for more accurate determination of the energies of the minima and transition state, compared to the other methods employed here. Classical barriers are converted to Arrhenius activation energy by TST calculations for comparisons with experimental results. The pre-exponential factors, A, calculated at all levels of theory are significantly larger than the experimental values. For activation energy, E-a DFT gives good results for HF elimination, within 4-8 W mol(-1) from experimental values. None of the methods employed, including CCSD(T), give comparable results for HCI elimination reactions. However, rate constants calculated by CCSD(T) method are in very good agreement with experiment for HCI elimination and they are in reasonable agreement for HF elimination reactions. Due to the strong correlation between A and E., the rate constants could be fit to a lower A and E-a (as given by experimental fitting, corresponding to a tight TS) or to larger A and E-a (as given by high level ab initio calculations, corresponding to a loose TS). The barrier for CIF elimination is determined to be 607 U mol(-1) at HF level and it is unlikely to be important for CH2FCH2Cl. Results for other CH2X-CH2Y (X,Y = F/Cl) are included for comparison.

Microwave spectrum and structure of C(6)H(5)CCH center dot center dot center dot H(2)S complex

Rotational spectra of C(6)H(5)CCH center dot center dot center dot H(2)S, C(6)H(5)CCH center dot center dot center dot H(2)(34)S, C(6)H(5)CCH center dot center dot center dot HDS, C(6)H(5)CCH center dot center dot center dot D(2)S and C(6) H(5)CCD center dot center dot center dot H(2)S complexes have been observed using a pulsed nozzle Fourier transform microwave spectrometer. The observed spectrum is consistent with a structure in which hydrogen sulfide is located over the phenyl ring pi cloud and the distance between the centers of masses of the two monomers is 3.74 +/- 0.01 angstrom. In the complex, the H(2)S unit is shifted from the phenyl ring center towards the acetylene group. The vibrationally averaged structure has an effective Cs symmetry. Ab initio calculations were performed at MP2/aug-cc-pVDZ level of theory to locate the possible geometries of the complex. The calculations reveal the experimentally observed structure to be more stable than a coplanar arrangement of the monomers, which was observed for the C(6)H(5)CCH center dot center dot center dot H(2)O complex. Atoms in molecule theoretical analysis shows the presence of S-H center dot center dot center dot pi hydrogen bond. For the parent isotopologue, each transition frequency was found to split into two resulting from an interchange of the equivalent hydrogens of H(2)S unit in the complex. (C) 2011 Elsevier Inc. All rights reserved.

Structure of the Pyridine-Chloranil Complex in Solution: A Surprise from Depolarized Hyper-Rayleigh Scattering Measurements

In this article, we report the structure of a 1:1 charge transfer complex between pyridine (PYR) and chloranil (CHL) in solution (CHCl(3)) from the measurement of hyperpolarizability (beta(HRS)) and linear and circular depolarization ratios, D and D', respectively, by the hyper-Rayleigh scattering technique and state-of-the-art quantum chemical calculations. Using linearly (electric field vector along X) and circularly polarized incident light, respectively, we have measured two macroscopic depolarization ratios D = I(X,X)(2 omega)/I(X,Z)(2 omega) and D' = I(X,C)(2 omega)/I(Z,C)(2 omega) in the laboratory fixed XYZ frame by detecting the second harmonic (SH) scattered light in a polarization resolved fashion. The stabilization energy and the optical gap calculated through the MP2/cc-pVDZ method using Gaussian09 were not significantly different to distinguish between the cofacial and T-shape structures. Only when the experimentally obtained beta(HRS) and the depolarization ratios, D and D', were matched with the theoretically computed values from single and double configuration interaction (SDCI) calculations performed using the ZINDO-SCRF technique, we concluded that the room temperature equilibrium structure of the complex is cofacial. This is in sharp contrast to an earlier theoretical prediction of the T-shape structure of the complex.

Conformation-Changing Aggregation in Hydroxyacetone: A Combined Low-Temperature FTIR, Jet, and Crystallographic Study

Aggregation in hydroxyacetone (HA) is studied using low-temperature FTIR, supersonic jet expansion, and X-ray crystallographic (in situ cryocrystallization) techniques. Along with quantum chemical methods (MP2 and DFT), the experiments unravel the conformational preferences of HA upon aggregation to dinners and oligomers. The O-H center dot center dot center dot O=C intramolecular hydrogen bond present in the gas-phase monomer partially opens upon aggregation in supersonic expansions, giving rise to intermolecular cooperatively enhanced O-H center dot center dot center dot O-H hydrogen bonds in competition with isolated O-H center dot center dot center dot O=C hydrogen bonds. On the other hand, low-temperature IR studies on the neat solid and X-ray crystallographic data reveal that HA undergoes profound conformational changes upon crystallization, with the HOCC dihedral angle changing from similar to 0 degrees in the gas phase to similar to 180 degrees in the crystalline phase, hence giving rise to a completely new conformation. These conclusions are supported by theoretical calculations performed on the geometry derived from the crystalline phase.

Thermal decomposition of haloethanols: single pulse shock tube and ab initio studies

This paper reports single pulse shock tube and ab initio studies on thermal decomposition of 2-fluoro and 2-chloroethanol at T=1000–1200 K. Both molecules have HX (X = F/Cl) and H2O molecular elimination channels. The CH3CHO formed by HX elimination is chemically active and undergoes secondary decomposition resulting in the formation of CH4, C2H6, and C2H4. A detailed kinetic simulation indicates that the formation of C2H4 could not be quantitatively explained as arising exclusively from secondary CH3CHO decomposition. Contributions from primary radical processes need to be considered to explain C2H4 quantitatively. Ab initio calculations on HX and H2O elimination reactions from the haloethanols at HF, MP2, and DFT levels with various basis sets up to 6/311++G**are reported. It is pointed out that due to strong correlations between A and Eα, comparison of these two parameters between experimental and theoretical results could be misleading.

Adsorption on Edge-Functionalized Bilayer Graphene Nanoribbons: Assessing the Role of Functional Groups in Methane Uptake

A combination of ab initio and classical Monte Carlo simulations is used to investigate the effects of functional groups on methane binding. Using Moller-Plesset (MP2) calculations, we obtain the binding energies for benzene functionalized with NH2, OH, CH3, COOH, and H2PO3 and identify the methane binding sites. In all cases, the preferred binding sites are located above the benzene plane in the vicinity of the benzene carbon atom attached to the functional group. Functional groups enhance methane binding relative to benzene (-6.39 kJ/mol), with the largest enhancement observed for H2PO3 (-8.37 kJ/mol) followed by COOH and CH3 (-7.77 kJ/mol). Adsorption isotherms are obtained for edge-functionalized bilayer graphene nanoribbons using grand canonical Monte Carlo simulations with a five-site methane model. Adsorbed excess and heats of adsorption for pressures up to 40 bar and 298 K are obtained with functional group concentrations ranging from 3.125 to 6.25 mol 96 for graphene edges functionalized with OH, NH2, and COOH. The functional groups are found to act as preferred adsorption sites, and in the case of COOH the local methane density in the vicinity of the functional group is found to exceed that of bare graphene. The largest enhancement of 44.5% in the methane excess adsorbed is observed for COOH-functionalized nanoribbons when compared to H terminated ribbons. The corresponding enhancements for OH- and NH2-functionalized ribbons are 10.5% and 3.7%, respectively. The excess adsorption across functional groups reflects the trends observed in the binding energies from MP2 calculations. Our study reveals that specific site functionalization can have a significant effect on the local adsorption characteristics and can be used as a design strategy to tailor materials with enhanced methane storage capacity.