Vibrational frequencies and structural determinations of 1,4-thioxane.

We present a detailed analysis of the structure and infrared spectra of 1,4-thioxane. The vibrational frequencies of the 1,4-thioxane molecule were analyzed using standard quantum chemical techniques. Frequencies were calculated at the MP2 and DFT levels of theory using the standard 6-31G* basis set. The structural transformation of the chemical agent bis (2-chloroethyl) sulfide (HD, mustard gas) and the related symmetry to a previously study compound(4) makes the symmetry of the 1,4-thioxane molecule an interesting candidate for study. The molecule exists normally in a Cs configuration similar to the chair form of cyclohexane. High-energy forms of 1,4-thioxane with C1 and C2 symmetry also exist.     

Kinetics of ion-molecule reactions with 2-chloroethyl ethyl sulfide at 298 K: a search for CIMS schemes for mustard gas

The rate constants and product ion branching ratios have been measured in a selected ion flow tube (SIFT) at 298 K for a variety of positive and negative ions reacting with 2-chloroethyl ethyl sulfide (2-CEES), a surrogate for mustard gas (HD). This series of experiments is designed to elucidate ion-molecule reactions that have large rate constants and produce unique product ions to guide the development of chemical ionization mass spectrometry (CIMS) detection methods for the chemical weapon agent using the surrogate instead. The negative ions typically used in CIMS instruments are essentially unreactive with 2-CEES, that is, SF 6 (-), SF 4 (-), CF 3O (-), and CO 3 (-). A few negative ions such as NO 2 (-) and NO 3 (-) undergo three-body association to give a unique product ion, but the bimolecular rate constants are small in the SIFT. Positive ions typically react at the collisional limit, primarily by charge and proton transfer, some of which is dissociative. For ions with high proton binding energies, association with 2-CEES has also been observed. Many of these reactions produced ions with the 2-CEES intact, including the parent cation, the protonated cation, and clusters. G3(MP2) calculations of the thermochemical properties for 2-CEES and mustard have been performed, along with calculations of the structures for the observed product cations. Reacting a series of protonated neutral molecules with 2-CEES brackets the proton affinity (PA) to between 812 ((CH 3) 2CO) and 854 (NH 3) kJ mol (-1). G3(MP2) calculations give a PA for 2-CEES of 823 kJ mol (-1) and a PA for mustard of 796 kJ mol (-1), indicating that the present results for 2-CEES should be directly transferable to mustard to design a CIMS detection scheme.     

Theoretical calculation of the OH vibrational overtone spectra of 1,5-pentanediol and 1,6-hexanediol

It is well-known that intramolecular hydrogen bonding affects the relative energetics of conformers, as well as the OH stretching peak positions, intensities, and width. In this study we simulated the Δv(OH) = 3, 4 overtone spectra of 1,5-pentanediol (PeD) and 1,6-hexanediol (HD) using the peak positions, intensities, and width calculated from the B3LYP/6-31+G(d,p) method. Furthermore, room temperature free energy calculations were performed using B3LYP/6-31+G(d,p) MP2/6-31+G(d,p), and MP2/6-311++G(3df,3pd) to obtain the relative population of the conformers. From the calculation of 109 and 381 distinct conformers for PeD and HD, respectively, we find that for these long chain diols the intramolecular hydrogen bonded conformers are not the most dominant conformation at room temperature. This is in stark contrast with shorter chain diols such as ethylene glycol for which the hydrogen bonded conformer dominates the population at room temperature. On the other hand, we found that the correlation between the hydrogen bonded OH red shift versus the homogeneous width, Γ = 0.0155(Δω)(1.36), which was derived for shorter chain diols, is valid even for these longer chain diols. We also showed that the intramolecular hydrogen bonded OH initially decays through the CCOH torsion and COH bending mode no matter how long the alkanediol chain length is for 1,n-alkanediols for n up to 6.     

Hydride affinity scale of various substituted arylcarbeniums in acetonitrile

Combined with the integral equation formalism polarized continuum model (IEFPCM), the hydride affinities of 96 various acylcarbenium ions in the gas phase and CH(3)CN were estimated by using the B3LYP/6-31+G(d)//B3LYP/6-31+G(d), B3LYP/6-311++G(2df,2p)//B3LYP/6-31+G(d), and BLYP/6-311++G(2df,2p)//B3LYP/6-31+G(d) methods for the first time. The results show that the combination of the BLYP/6-311++G(2df,2p)//B3LYP/6-31+G(d) method and IEFPCM could successfully predict the hydride affinities of arylcarbeniums in MeCN with a precision of about 3 kcal/mol. On the basis of the calculated results from the BLYP method, it can be found that the hydride affinity scale of the 96 arylcarbeniums in MeCN ranges from -130.76 kcal/mol for NO(2)-PhCH(+)-CN to -63.02 kcal/mol for p-(Me)(2)N-PhCH(+)-N(Me)(2), suggesting most of the arylcarbeniums are good hydride acceptors. Examination of the effect of the number of phenyl rings attached to the carbeniums on the hydride affinities shows that the increase of the hydride affinities takes place linearly with increasing number of benzene rings in the arylcarbeniums. Analyzing the effect of the substituents on the hydride affinities of arylcarbeniums indicates that electron-donating groups decrease the hydride affinities and electron-withdrawing groups show the opposite effect. The hydride affinities of arylcarbeniums are linearly dependent on the sum of the Hammett substituent parameters σ(p)(+). Inspection of the correlation of the solution-phase hydride affinities with gas-phase hydride affinities and aqueous-phase pK(R)(+) values reveals a remarkably good correspondence of ΔG(H(-)A)(R(+)) with both the gas-phase relative hydride affinities only if the α substituents X have no large electron-donating or -withdrawing properties and the pK(R)(+) values even though the media are dramatically different. The solution-phase hydride affinities also have a linear relationship with the electrophilicity parameter E, and this dependence can certainly serve as one of the most effective ways to estimate the new E values from ΔG(H(-)A)(R(+)) or vice versa. Combining the hydride affinities and the reduction potentials of the arylcarbeniums, we obtained the bond homolytic dissociation Gibbs free energy changes of the C-H bonds in the corresponding hydride adducts in acetonitrile, ΔG(HD)(RH), and found that the effects of the substituent on ΔG(HD)(RH) are very small. Simple thermodynamic analytic platforms for the three C-H cleavage modes were constructed. It is evident that the present work would be helpful in understanding the nature of the stabilities of the carbeniums and mechanisms of the hydride transfers between carbeniums and other hydride donors.     

Theoretical study of light-metal tetrahydroaluminates

Ab initio calculations of the structural, energetic, vibrational, and magnetic characteristics of the lowest-lying structures for isolated molecules and ions of light-metal tetrahydroaluminates (Li, Na, Be, Mg, and Al) have been performed by the perturbation theory (MP2), quadratic configuration interaction (QCISD(T)), coupled cluster (CCSD(T)), and density functional theory (B3LYP) methods using the 6-31G*, 6-31G**, 6-311+G**, and 6-311++G** basis sets. The trends in the behavior of the molecular characteristics have been analyzed in various related series of these compounds. The results are compared with the data on analogous light-metal tetrahydroborates calculated at the same levels of theory. The differences in structure and stability between analogous hydroborate and alanate complexes are examined. The economical approximation B3LYP/6-311++G**//B3LYP/6-311+G** has been shown to adequately reproduce the results obtained at the higher level of theory CCSD(T)/6-311++G**//MP2/6-31G* even though it requires considerably shorter CPU times and smaller amounts of memory.     

二氟硅杂环丙烯异构化反应的理论研究

近十几年来，硅烯（：SN小'一到作为活性中间体的研究引起化学界的广泛兴趣，形成了内容丰富的硅烯化学．1975年Ilass＊等研究了硅烯和乙炔的加成反应，指出该反应的中间体为硅杂环丙烯，硅杂环丙烯异构化为硅甲基乙炔．Boatz问等利用：j－ZIG（d）基组对金属杂环丙烯小［＊xZC     

Selective and sensitive trace analysis of sulfur mustard with thermal desorption and two-dimensional gas chromatography-mass spectrometry

An improved method is presented for the trace analysis of sulfur mustard (HD) in biological samples, such as blood and tissue from laboratory animals. Using the internal standard method and liquid-liquid extraction with ethyl acetate, up to 400 microL of the extract was injected by thermal desorption from Tenax and analyzed by two-dimensional GC-MS/EI in SIM mode. The analysis was compared with a direct GC injection. Reversed thermal desorption was used as a tool for handling heavily contaminated (fat) samples, thus preventing contamination of the injection system and pre-column. A successful analytical configuration has been set up for the bioanalysis of HD at the low, toxicologically relevant pM level. A detection limit of 10 pg mL(-1) blood or pg g(-1) tissue of sulfur mustard (S/N=3) was established by using this configuration.     

2-丁基-四氢噻吩亚砜13C-NMR的理论研究

在RHF/6-31G和B3LYP/6-31G水平上对顺式(Cis-)与反式(Trans-)2-丁基-四氢噻吩亚砜(BTHTO)进行几何优化,应用规范不变原子轨道法(GIAO)在6-31G、6-31+G、6-31++G和6-31+G(2d,p)水平上计算了Cis-和Trans-BTHTO的13C-NMR,对13C-NMR谱进行了归属。结果表明,BTHTO噻吩五元环的稳定构象呈半椅式,Cis-和Trans-BTHTO中与硫原子直接碳原子13C-NMR的显著差异主要是由于空间构型不同引起分子的静电势场对相应碳原子的屏蔽作用不同所致。     

A density functional theory study of the structure and vibrational spectra of beta-carotene, capsanthin, and capsorubin

A theoretical study of the structure and the vibrational spectra of the beta-carotene molecule and its derivatives capsanthin and capsorubin is carried out. We first investigate systematically the theoretical method which provides the best results for beta-carotene by performing ab initio calculations at the HF/6-31G(d), SVWN/6-31G(d), PBE0/6-31G(d), BLYP/6-31G(d), B3LYP/6-31G(d), B3LYP/6-31G(d,p), B3LYP/6-311G(d), and B3LYP/6-311G(d,p) levels and by using previous theoretical results available in the literature obtained at the AM1 and BPW91/6-31G(d) levels. The influence of both the level of calculation and the size of the basis set used in the geometry optimization and in the determination of the IR and Raman spectra of this molecule is thus analyzed. It is confirmed that the hybrid functional B3LYP with the basis 6-31G(d) is the method that gives the best results as a whole. By use of this level of calculation, we next optimize the molecular geometries of related molecules of capsanthin and capsorubin, which to the best of our knowledge have only been studied at the semiempirical AM1 level. In addition we calculate the IR and Raman spectra of these molecules at the B3LYP/6-31G(d) level of theory. The results obtained for capsanthin show on the one hand that the double bond of the beta-ionone ring is outside the polyene chain plane, due to the repulsion between the hydrogen atoms of the ring methyl groups and the hydrogen atoms of the polyene chain, and on the other hand that the carbonyl double bond in the other headgroup is very close to planarity with the polyene chain, since in this case such a repulsion does not exist. For the molecule of capsorubin the two carbonyl groups also take the same coplanar orientation relative to the polyene chain. The IR and Raman spectra theoretically computed for these two molecules are finally compared with their experimental spectra and the vibrational normal modes of the main signals are interpreted.     

An ab initio study of intermolecular interactions of nitromethane dimer and nitromethane trimer

Different geometries of nitromethane dimer and nitromethane trimer have been fully optimized employing the density functional theory B3LYP method and the 6-31++G** basis set. Three-body interaction energy has been obtained with the ab initio supermolecular approach at the levels of MP2/6-31++G**//B3LYP/6-31++G** and MP2/aug-cc-pVDZ//B3LYP/6-31++G**. The internal rotation of methyl group induced by intermolecular interaction has been observed theoretically. For the optimized structures of nitromethane dimer, the strength of C--H...O--N H-bond ranges from -9.0 to -12.4 kJ mol(-1) at the MP2/aug-cc-pVDZ//B3LYP/6-31++G** level, and the B3LYP method underestimates the interaction strength compared with the MP2 method, while MP2/6-31++G**//B3LYP/6-31++G** calculated DeltaE(C) is within 2.5 kJ mol(-1) of the corresponding value at the MP4(SDTQ)/6-31G**//B3LYP/6-31++G** level. The analytic atom-atom intermolecular potential has been successfully regressed by using the MP2/6-31++G**//B3LYP/6-31++G** calculated interaction energies of nitromethane dimer. For the optimized structures of nitromethane trimer the three-body interaction energies occupy small percentage of corresponding total binding energies, but become important for the compressed nitromethane explosive. In addition, it has been discovered that the three-body interaction energy in the cyclic nitromethane trimer is more and more negative as intermolecular distances decrease from 2.2 to 1.7 A.     

双自由基CF2与O3的反应机理

主要用作致冷剂和发泡剂的氯氟烃（CFCs）是破坏臭氧层的主要物质之一．对氯氟烃类化合物及其降解产物（包括光解、光氧化、化学反应产物等）在大气中行为问题的研究是大气化学研究的重要内容．前人[1-3]从理论和实验两方面研究了自由基与臭氧的反应机制，但是氯氟烃光解过程中     

硝酸甲与硝酸乙酯的结构、振动频率和热力学性质的密度函数 理论(DFT)研究

用密度函数理论(DFT)的BLYP和B3LYP方法,取6-31G,6-31G^*,6-31G^*^*,6-311G,6-311G^*和6-311G^*^*六种基组,对硝酸甲酯和硝酸乙酯的几何构型和红外振动频率进行了计算研究.结果表明,B3LYP方法在采用极化基组(6-31G^*,6-31G^*^*,6-311G^*和6-311G^*^*)时计算得到的结果均较好,适用于硝酸酯类化合物的研究.而BLYP方法无论采用何种基组均不适用;运用校正后的B3LYP/6-31G^*频率(校正因子0.975)计算得到的热力学性质(C^o~p,H^o和S^o)与实验结果较吻合。     

Structure and electronic properties of tris(4-hydroxy-1,5-naphthyridinato) aluminum (AlND3) and its methyl derivatives: a theoretical study

Molecular structures of the ground state (S₀) of tris(4-hydroxy-1,5-naphthyridinato) aluminum (AlND3) and its methyl derivatives have been optimized using B3LYP/6-31G(D) and the first singlet excited state (S₁) using CIS/6-31G(D) method, respectively. The frontier molecular orbitals characteristics of these Al-complexes have been analyzed systematically in order to understand the electronic transitions. It is seen from these results that in all these complexes, like in earlier reported mer-Alq3, the highest occupied molecular orbital (HOMO) is localized on the pyridine-4-ol ring of A-ligand while lowest unoccupied molecular orbital (LUMO) is on the pyridyl ring of B-ligand for S₀ states irrespective of the methyl substitution present on the ligands. The absorption and emission wavelengths have been evaluated at the TD-PBE0/6-31G(D) level and found to be comparable with the experiment. The charge transfer integrals have been calculated for AlND3, and results reveals that electron transport is larger than hole transport. The reorganization energies have been calculated at B3LYP/6-31G(D) level for these complexes, and the results show that the charge mobilities of these complexes are comparable with mer-Alq3.     

Ab initio benchmark calculations on Ca(II) complexes and assessment of density functional theory methodologies

A set of benchmark results for the geometries, binding energies, and protonation affinities of 24 complexes of small organic ligands with Ca(II) is provided. The chosen level of theory is CCSD(T)/CBS obtained by means of a composite procedure. The performance of four density functionals, namely, PW91, PBE, B3LYP, and TPSS and several Pople-type basis sets, namely, 6-31G(d), 6-31+G(d), 6-31+G(2d,p) and 6-311+G(d) have been assessed. Additionally, the nature of the metal ligand bonding has been analyzed by means of the Symmetry Adapted Perturbation Theory (SAPT). We have found that the B3LYP hybrid functional, in conjunction with either the polarized double-ζ 6-31+G(2d,p) or the triple-ζ 6-311+G(d) basis sets, yields the closest results compared to the benchmark data. The SAPT analysis stresses the importance of induction effects in the binding of these complexes and suggests that consideration of classical electrostatic contributions alone may not be reliable enough for the prediction of relative binding energies for Ca(II) complexes.     

Interaction of anions with perfluoro aromatic compounds

The complexes formed by a variety of anions with perfluoro derivatives of benzene, naphthalene, pyridine, thiophene, and furan have been calculated using DFT (B3LYP/6-31++G**) and MP2 (MP2/6-31++G** and MP2/6-311++G**) ab initio methods. The minimum structures show the anion interacting with the pi-cloud of the aromatic compounds. The interaction energies obtained range between -8 and -19 kcal mol(-1). The results obtained at the MP2/6-31++G** and MP2/6-311++G** levels are similar. However, the B3LYP/6-31++G** results provide longer interaction distances and smaller interaction energies than do the MP2 results. The interaction energies have been partitioned using an electrostatic, polarization, and van der Waals scheme. The AIM analysis of the electron density shows a variety of topologies depending on the aromatic system considered.     

Gas-phase structure, rotational barrier, and vibrational properties of methyl methanethiosulfonate, CH3SO2SCH3: an experimental and computational study

The molecular structure of methyl methanethiosulfonate, CH3SO2SCH3, has been determined in the gas phase from electron-diffraction data supplemented by ab initio (HF, MP2) and density functional theory (DFT) calculations using 6-31G(d), 6-311++G(d,p), and 6-311G(3df,3pd) basis sets. Both experimental and theoretical data indicate that although both anti and gauche conformers are possible by rotating about the S-S bond, the preferred conformation is gauche. The barrier to internal rotation in the CSSC skeleton has been calculated using the RHF/6-31G(d), MP2/6-31G(d), and B3LYP/6-31G(d) methods as well as MP2 with a 6-31G(3df) basis set on sulfur and 6-31G(d) on C, H, and O. A 6-fold decomposition of the rotational barrier has been performed in terms of a Fourier-type expansion, enabling us to analyze the nature of the potential function, showing that the coefficients V1 and V2 are the dominant terms; V1 is associated with nonbonding interactions, and V2 is associated with hyperconjugative interactions. A natural bond orbital analysis showed that the lone pair --> sigma* hyperconjugative interactions favor the gauche conformation. Furthermore, the infrared spectra for the liquid and solid phases and the Raman spectrum for the liquid have been recorded, and the observed bands have been assigned to the vibrational normal modes. The experimental vibrational data, along with calculated theoretical force constants, were used to define a scaled quantum mechanical force field for the target system that enabled us to estimate the measured frequencies with a final root-mean-square deviation of 6 cm-1.     

2,3,7,8-四氯苯并二英分子结构与热力学性质的理论研究

用密度泛函理论(DFT)和从头算(ab initio)方法,在B3LYP/6-31G、 B3LYP/6-31G*、 B3LYP/6-311G*和MP2/6-31G*水平上全优化计算了2,3,7,8-四氯苯并二英(2,3,7,8-TCDD)的几何构型、电子结构和振动频率,并用校正后的频率计算了298～1500 K的标准热力学函数,同时用半经验的PM3 SCF-MO进行了同样的计算,计算结果与实验值及文献值较好地吻合.     

Electronic, optical, and charge transfer properties of donor–bridge–acceptor hydrazone sensitizers

The ground state geometries have been computed by using density functional theory (DFT) at B3LYP/6-31G*, B3LYP/6-31G**, and PCM-B3LYP/6-31G* level of theories. The highest occupied molecular orbitals (HOMOs) are delocalized on whole of the molecule and the lowest unoccupied molecular orbitals (LUMOs) are localized on the tricarbonitrile. The lowest HOMO and LUMO energies have been observed for Dye1 while highest for Dye4. The LUMO energies of Dye1–Dye4 are above the conduction band of TiO₂ and HOMOs are below the redox couple. The absorption spectra have been computed in solvent (methanol) and without solvent by using time-dependant DFT at TD-B3LYP/6-31G*, TD-B3LYP/6-31G**, and PCM-TD-B3LYP/6-31G* level of theories. The calculated maximum absorption wavelengths of the spectra in methanol are in good agreement with experimental evidences. The maximum absorption wavelengths of new designed sensitizers are red shifted compared to parent molecule. The electronic coupling constant and electron injection have been computed by first principle investigations. The improved electronic coupling constant and electron injection revealed that new modeled systems would be efficient sensitizers.     

The shape of gaseous n-butylbenzene: assessment of computational methods and comparison with experiments

Conformational landscape of neutral and ionized n-butylbenzene has been examined. Geometries have been optimized at the B3LYP/6-31G(d), B3LYP/6-31+G(d,p), B3LYP-D/6-31+G(d,p), B2PLYP/6-31+G(d,p), B2PLYP-D/6-31+G(d,p), B97-D/6-31+G(d,p), and M06-2X/6-31+G(d,p) levels. This study is complemented by energy computations using 6-311++G(3df,2p) basis set and CBS-QB3 and G3MP2B3 composite methods to obtain accurate relative enthalpies. Five distinguishable conformers have been identified for both the neutral and ionized systems. Comparison with experimentally determined rotational constants shows that the best geometrical parameters are provided by B3LYP-D and M06-2X functionals, which include an explicit treatment of dispersion effects. Composite G3MP2B3 and CBS-QB3 methods, and B2PLYP-D, B3LYP-D, B97-D, and M06-2X functionals, provide comparable relative energies for the two sets of neutral and ionized conformers of butyl benzene. An exception is noted however for conformer V(+) the stability of which being overestimated by the B3LYP-D and B97-D functionals. The better stability of neutral conformers I, III, and IV, and of cation I(+) , demonstrated by our computations, is in perfect agreement with conclusions based on micro wave, fluorescence, and multiphoton ionization experiments.     

FT-IR, FT-Raman spectra, density functional computations of the vibrational spectra and molecular geometry of butylated hydroxy toluene

The FT-IR spectrum of 2,6-di-tert-butyl-4-methylphenol [butylated hydroxy toluene] was recorded in the region 4000-400 cm(-1). The FT-Raman spectrum of butylated hydroxy toluene was also recorded in the region 3500-50 cm(-1). The molecular structure and vibrational frequencies of butylated hydroxy toluene (BHT) have been investigated with combined experimental and theoretical study. Two stable conformers of the title compound were obtained from the result of geometry optimizations of these possible conformers. The conformer 1 is (approximately 2.6 kcal/mol) more stable than conformer 2. Geometry optimizations and vibrational frequency calculations were performed by BLYP and B3LYP methods using 6-31G(d), 6-31G(d,p) and 6-31+G(d,p) as basis sets. The scaled frequencies were compared with experimental spectrum and on the basis of this comparison; assignments of fundamental vibrational modes were examined. Comparison of the experimental spectra with harmonic vibrational wavenumbers indicates that B3LYP/6-31G(d) results are more accurate. Predicted electronic absorption spectra of BHT from TD-DFT calculation have been analyzed and compared with the experimental UV-vis spectrum. The calculated HOMO and LUMO energies show that the charge transfer occurs within the molecule.