丁烯和丁烷自由基阳离子的分子结构和超精细结构的 理论研究

用密度函数理论B3LYP方法和6-31G(d,p),6-311G(d,p)及6-311+G(d,p)基组，分别对1-C4H^+~8,2-C4H^+~8和C4H^+~10进行了构型优化和频率分析计算，预言1-C4H^+~8具有非平面构型，与以往报道的从头算和密度函数理论计算结果不同。在各自由基阳离子的B3LYP构型上，进行了B3LYP、MP2及MRSDCI方法的超精细偶合常数计算，得到了比以往更好的结果，特别是MP2/B3LYP计算值是至今与实验值符合得最好的理论计算结果。     

槲皮素与桑色素ESI-MS裂解行为的比较分析

运用量子化学方法辅助解析并比较槲皮素与桑色素在电喷雾离子阱质谱(ESI-MS)负离子模式下的裂解行为。依据密度泛函理论(Density functional theory,DFT),在B3LYP/6-31G(d)水平,对槲皮素与桑色素的分子空间构型进行优化,确定稳定的几何构型与去质子化位点,在RB3LYP/6-31+G(2d,2p)水平,计算相对碎裂电压下的二级质谱中碎片离子处于稳定状态时的能量,通过比较准分子离子稳定构型并结合基组重叠误差(Basis set superposition error,BSSE)校正后的键解离能(Bond dissociation energy,BDE),推导了质谱碎裂过程。结果显示:槲皮素的稳定构型为A,B,C环处于同一平面,桑色素上的2'-OH使得B环与AC环之间翻转一定角度,二面角D(1,2,1',6')为-134.662 4°。槲皮素与桑色素的质谱裂解过程主要通过C环跨环裂解产生,且具有多种开裂方式,开裂先后顺序为:1,2开裂、0,2开裂、1,3开裂、1,4开裂与0,4开裂,分别生成碎片离子1,2A-,0,2A-,1,3A-,1,4A-与0,4A-,并逐步进行后续裂解,而2'-OH的存在促进了桑色素的裂解。该研究为进一步揭示黄酮醇类化合物的质谱裂解规律提供了理论依据。     

山奈酚的电喷雾质谱裂解途径

采用电喷雾质谱技术对山奈酚的质谱过程进行了表征,并用量子化学方法对山奈酚的质谱裂解途径进行了理论研究。 依据密度泛函理论,在B3LYP/6-31G(d)水平,对山奈酚的各质谱碎片离子进行了几何结构优化,确定了各碎片的稳定结构;然后,在ROB3LYP/ 6-311++G(2d,2p)水平计算了形成各碎片离子所需的键断裂能,进而推导出了山奈酚的质谱碎裂途径。 计算结果显示,山奈酚主要通过C环碎裂而发生裂解,出现碎片m/z 284.7、256.7、228.7、210.7、184.8、168.7和150.7,其中m/z 210.7的碎片离子键断裂能最小,m/z 150.7的碎片离子键断裂能最大,说明前者很容易由母离子碎片形成,后者较难由母离子碎片形成。     

Collision induced dissociations of non‐derivatized and trimethylsilyl‐derivatized estradiols: similarities in fragmentation patterns

Fragmentation mechanisms of estradiol and trimethylsilyl (TMS)‐derivatized estradiol were studied by triple quadrupole tandem mass spectrometry (MSMS) and density functional theory (DFT) at B3LYP/6‐311G(d,p) level. Collision induced dissociations (CID) of estradiol give product ions that are associated with the cleavage of B, C and D rings. Characteristic fragments from the cleavage of the aromatic ring A were not identified, and this was confirmed with both labeled estradiol and trimethylsilyl (TMS)‐derivatized estradiol. The mechanisms are based on charge‐site directed, radical‐directed and charge remote fragmentations that are consistent with previous studies of steroids. CID spectra show ion pairs at m/z: 145/146, 157/158, 185/186, 211/213 and 225/226 with significant intensities, suggesting that these pairs are not from isotopic contributions. The mechanisms show similarities with some minor differences in the fragmentation patterns between the non‐derivatized and the TMS‐derivatized estradiol. Copyright © 2015 John Wiley & Sons, Ltd.     

Mechanistic insights into triterpene synthesis from quantum mechanical calculations. Detection of systematic errors in B3LYP cyclization energies

Most quantum mechanical studies of triterpene synthesis have been done on small models. We calculated mPW1PW91/6-311+G(2d,p)//B3LYP/6-31G* energies for many C30H51O+ intermediates to establish the first comprehensive energy profiles for the cationic cyclization of oxidosqualene to lanosterol, lupeol, and hopen-3beta-ol. Differences among these 3 profiles were attributed to ring strain, steric effects, and proton affinity. Modest activation energy barriers and the ample exothermicity of most annulations indicated that the cationic intermediates rarely need enzymatic stabilization. The course of reaction is guided by hyperconjugation of the carbocationic 2p orbital with parallel C-C and C-H bonds. Hyperconjugation for cations with a horizontal 2p orbital (in the plane of the ABCD ring system) leads to annulation and ring expansion. If the 2p orbital becomes vertical, hyperconjugation fosters 1,2-methyl and hydride shifts. Transition states leading to rings D and E were bridged cyclopropane/carbonium ions, which allow ring expansion/annulation to bypass formation of undesirable anti-Markovnikov cations. Similar bridged species are also involved in many cation rearrangements. Our calculations revealed systematic errors in DFT cyclization energies. A spectacular example was the B3LYP/6-311+G(2d,p)//B3LYP/6-31G* prediction of endothermicity for the strongly exothermic cyclization of squalene to hopene. DFT cyclization energies for the 6-311+G(2d,p) basis set ranged from reasonable accuracy (mPW1PW91, TPSSh with 25% HF exchange) to underestimation (B3LYP, HCTH, TPSS, O3LYP) or overestimation (MP2, MPW1K, PBE1PBE). Despite minor inaccuracies, B3LYP/6-31G* geometries usually gave credible mPW1PW91 single-point energies. Nevertheless, DFT energies should be used cautiously until broadly reliable methods are established.     

C alpha-C backbone fragmentation dominates in electron detachment dissociation of gas-phase polypeptide polyanions

Fragmentation of peptide polyanions by electron detachment dissociation (EDD) has been induced by electron irradiation of deprotonated polypeptides [M-nH](n-) with >10 eV electrons. EDD has been found to lead preferentially to a* and x fragment ions (C(alpha)-C backbone cleavage) arising from the dissociation of oxidized radical anions [M-nH]((n-1)-*. We demonstrate that C(alpha)-C cleavages, which are otherwise rarely observed in tandem mass spectrometry, can account for most of the backbone fragmentation, with even-electron x fragments dominating over radical a* ions. Ab initio calculations at the B3 LYP level of theory with the 6-311+G(2 p,2 d)//6-31+G(d,p) basis set suggested a unidirectional mechanism for EDD (cleavage always N-terminal to the radical site), with a*, x formation being favored over a, x* fragmentation by 74.2 kJ mol(-1). Thus, backbone C(alpha)-C bonds N-terminal to proline residues should be immune to EDD, in agreement with the observations. EDD may find application in mass spectrometry for such tasks as peptide sequencing and localization of labile post-translational modifications, for example, those introduced by sulfation and phosphorylation. EDD can now be performed not only in Fourier transform mass spectrometry, but also in far more widely used quadrupole (Paul) ion traps.     

Structural elucidation of nitro-substituted five-membered aromatic heterocycles utilizing GIAO DFT calculations

The GIAO (Gauge Including Atomic Orbitals) DFT (Density Functional Theory) method is applied at the B3LYP/6-31+G(d,p)//B3LYP/6-31+G(d), B3LYP/6-311++G(d,p)//B3LYP/6-31+G(d), B3LYP/6-311+G (2d,p)//B3LYP/6-31+G(d) and B3LYP/6-311++G(d,p)//B3LYP/6-311++G(d,p) levels of theory for the calculation of proton and carbon chemicals shifts and coupling constants for 25 nitro-substituted five-membered heterocycles. Difference (1D NOE) spectra in combination with long-range gHMBC experiments were used as tools for the structural elucidation of nitro-substituted five-membered heterocycles. The assigned NMR data (chemical shifts and coupling constants) for all compounds were found to be in good agreement with theoretical calculations using the GIAO DFT method. The magnitudes of one-bond (1JCH) and long-range (nJCH, n>1) coupling constants were utilized for unambiguous differentiation between regioisomers of nitro-substituted five-membered heterocycles.     

Inductive effects in radicals calculated from DFT energies; substituted bicyclo[2.2.2]octan-1-yloxy radicals

Energies of a series of 4-substituted 1-oxybicyclo[2.2.2]octan-1-yloxy radicals with 18 various substituents were calculated within the framework of the DFT theory at the levels UB3LYP/6-311+G(d,p)//UB3LYP/6-311+G(d,p) and UB3LYP/6-311++G(2df,p)//UB3LYP/6-311+G(d,p) and compared with similar series of the parent alcohols, their deprotonated and protonated forms calculated at the levels B3LYP/6-311+G(d,p)//B3LYP/6-311+G(d,p) and B3LYP/6-311++G(2df,p)//B3LYP/6-311+G(d,p). The two levels are of the same performance and both are sufficient for molecules of this type according to comparison with scarce experimental gas-phase acidities and basicities. The substituent effects were analyzed in terms of isodesmic equations. In addition to strong dependence on the substituent inductive effect, a slight dependence on the electronegativity of the first atom of the substituent was proven in certain cases. In all aspects, there is no qualitative difference between the effects on radicals and on similar closed shell species. Radicals behave as slightly electron deficient; the substituent effect is weaker than that on the ions but stronger than on neutral molecules.     

Theoretical study of the conformational energy hypersurface of cyclotrisarcosyl

The multidimensional Conformational Potential Energy Hypersurface (PEHS) of cyclotrisarcosyl was comprehensively investigated at the DFT (B3LYP/6-31G(d), B3LYP/6-31G(d,p) and B3LYP/6-311++G(d,p)), levels of theory. The equilibrium structures, their relative stability, and the Transition State (TS) structures involved in the conformational interconversion pathways were analyzed. Aug-cc-pVTZ//B3LYP/6-311++G(d,p) and MP2/6-31G(d)//B3LYP/6-311++G(d,p) single point calculations predict a symmetric cis-cis-cis crown conformation as the energetically preferred form for this compound, which is in agreement with the experimental data. The conformational interconversion between the global minimum and the twist form requires 20.88 kcal mol-1 at the MP2/6-31G(d)//B3LYP/6-311++G(d,p) level of theory. Our results allow us to form a concise idea about the internal intricacies of the PEHSs of this cyclic tripeptide, describing the conformations as well as the conformational interconversion processes in this hypersurface. In addition, a comparative analysis between the conformational behaviors of cyclotrisarcosyl with that previously reported for cyclotriglycine was carried out     

Theoretical investigation on the GaCl3-catalyzed ring-closing metathesis reaction of N-2,3-butadienyl-2-propynyl-1-amine: three-membered ring versus four-membered ring mechanism

The gallium chloride (GaCl(3))-catalyzed ring-closing metathesis reaction mechanism of N-2,3-butadienyl-2-propynyl-1-amine has been studied at the Becke three-parameter hybrid functional combined with Lee-Yang-Parr correlation functional (B3LYP)/6-31G(d), B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p)//B3LYP/ 6-31G(d) and the second-order M?ller-Plesset perturbation (MP2)/6-311++G(d,p)//B3LYP/6-31+G(d,p) levels. It was found that the final metathesis product can be yielded via a three-membered or four-membered ring mechanism. The three-membered ring pathway is favorable due to its low energy barrier at the rate determining step. The whole reaction is stepwise and strongly exothermic.     

A novel salt bridge mechanism highlights the need for nonmobile proton conditions to promote disulfide bond cleavage in protonated peptides under low-energy collisional activation

The gas-phase fragmentation mechanisms of small models for peptides containing intermolecular disulfide links have been studied using a combination of tandem mass spectrometry experiments, isotopic labeling, structural labeling, accurate mass measurements of product ions, and theoretical calculations (at the MP2/6-311 + G(2d,p)//B3LYP/3-21G(d) level of theory). Cystine and its C-terminal derivatives were observed to fragment via a range of pathways, including loss of neutral molecules, amide bond cleavage, and S-S and C-S bond cleavages. Various mechanisms were considered to rationalize S-S and C-S bond cleavage processes, including charge directed neighboring group processes and nonmobile proton salt bridge mechanism. Three low-energy fragmentation pathways were identified from theoretical calculations on cystine N-methyl amide: (1) S-S bond cleavage dominated by a neighboring group process involving the C-terminal amide N to form either a protonated cysteine derivative or protonated sulfenyl amide product ion (44.3 kcal mol(-1)); (2) C-S bond cleavage via a salt bridge mechanism, involving abstraction of the alpha-hydrogen by the N-terminal amino group to form a protonated thiocysteine derivative (35.0 kcal mol(-1)); and (3) C-S bond cleavage via a Grob-like fragmentation process in which the nucleophilic N-terminal amino group forms a protonated dithiazolidine (57.9 kcal mol(-1)). Interestingly, C-S bond cleavage by neighboring group processes have high activation barriers (63.1 kcal mol(-1)) and are thus not expected to be accessible during low-energy CID experiments. In comparison to the energetics of simple amide bond cleavage, these S-S and C-S bond cleavage reactions are higher in energy, which helps rationalize why bond cleavage processes involving the disulfide bond are rarely observed for low-energy CID of peptides with mobile proton(s) containing intermolecular disulfide bonds. On the other hand, the absence of a mobile proton appears to "switch on" disulfide bond cleavage reactions, which can be rationalized by the salt bridge mechanism. This potentially has important ramifications in explaining the prevalence of disulfide bond cleavage in singly protonated peptides under MALDI conditions.     

Mass spectrometric characterization of 4-oxopentanoic acid and gas-phase ion fragmentation mechanisms studied using a triple quadrupole and time-of-flight analyzer hybrid system and density functional theory

4-Oxopentanoic acid was characterized experimentally by electrospray ionization using a triple quadrupole and time-of-flight analyzer hybrid system. This compound was chosen as a model substance for small organic compounds bearing an acetyl and a carboxyl group. Collision-induced dissociation experiments at different activation energies were performed to elucidate possible fragmentation pathways. These pathways were also studied on the theoretical level using density functional theory (DFT) B3LYP/6-311++G(3df,3pd)//B3LYP/6-31+G(d)+ZPVE calculations. CO2 ejection from the [M-H](-) anion of 4-oxopentanoic acid was observed and the fragmentation pathway studied by DFT reveals a new concerted mechanism for CO2 elimination accompanied by an intramolecular proton transfer within a pentagonal transition state structure. Successive elimination of water and CO from the [M-H](-) anion of 4-oxopentanoic acid was also observed. A rearrangement in the primary deprotonated ketene anion produced after water elimination was found on the theoretical level and leads to CO elimination from the primary product anion [M-H-H2O](-). Energy diagrams along the reaction coordinates of the fragmentation pathways are presented and discussed in detail. Mulliken charge distributions of some important structures are presented.     

NMR spectra, GIAO and charge density calculations of five-membered aromatic heterocycles

The B3LYP/6-31+G(d) molecular geometry optimized structures of 17 five-membered heterocycles were employed together with the gauge including atomic orbitals (GIAO) density functional theory (DFT) method at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p) and B3LYP/6-311+G(2d,p) levels of theory for the calculation of proton and carbon chemicals shifts and coupling constants. The method of geometry optimization for pyrrole (1), N-methylpyrrole (2) and thiophene (7) using the larger 6-311++G(d,p) basis sets at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,p) and B3LYP/cc-pVTZ levels of theory gave little difference between calculated and experimental values of coupling constants. In general, the (1)H and 13C chemical shifts for all compounds are in good agreement with theoretical calculations using the smaller 6-31 basis set. The values of nJHH(n=3, 4, 5) and rmnJ(CH)(n=1, 2, 3, 4) were predicted well using the larger 6-31+G(d,p) and 6-311++G(d,p) basis sets and at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,2p) levels of theory. The computed atomic charges [Mülliken; Natural Bond Orbital Analysis (NBO); Merz-Kollman (MK); CHELP and CHELPG] for the B3LYP/6-311++G(d,p) geometry optimized structures of 1-17 were used to explore correlations with the experimental proton and carbon chemical shifts.     

Theoretical Study on the Structures and Isomerization of Germylenoid GeH_2LiCl

Using geometrical optimization and DFT method at the B3LYP/6-311++G (3df,3pd) level, four equilibrium geometries and one transition state of GeH2LiCl were identified, and the structures at the MP2/6-311++G(3df,3pd) level were calculated simultaneously. We also studied the solvent effects on the structures of Germylenoid GeH2LiCl at the B3LYP/6-311++G (3df,3pd) level. The two more stable forms are suggested to be the p-complex and three-membered ring. The vibrational frequencies and infrared intensities were computed at the B3LYP/6-311++G (3df,3pd) level.     

戊烯自由基阳离子的密度泛函理论研究

使用密度泛函理论B3LYP方法和6-31G(d,p)、6-31＋G(d,p)、6-311G(d,p)及6-311＋G(d,p)基组，分别对2-C5H10＋和1-C5H10＋的各种构象进行了几何构型优化，并用B3LYP/6-311G(d,p)进行了频率分析计算.计算预言1-C5H10＋具有非平面构型，与以往报导的从头算计算结论相反.在两个自由基阳离子的各种构象的B3LYP几何构型上，进行了B3LYP和UMP2(full)方法的超精细偶合常数计算，得到了比以往更好的结果.     

Mass spectrometric characterization of beta-caryophyllene ozonolysis products in the aerosol studied using an electrospray triple quadrupole and time-of-flight analyzer hybrid system and density functional theory

The components of the organic aerosol formed due to gas-phase beta-caryophyllene ozonolysis were characterized by the use of a triple quadrupole and time-of-flight analyzer hybrid system coupled to an electrospray ionization source operated in the negative ion mode. A reversed-phase high-performance liquid chromatography (HPLC) column was used to achieve chromatographic separations at neutral pH which has been proved to induce ionization of organic compounds bearing aldehyde moieties. In addition to the detected oxo- and dicarboxylic acids, isomeric oxidation products, which bear multi-functional groups such as aldehyde, carbonyl and hydroxyl groups, could be differentiated by examining their corresponding collision-induced dissociation (CID) fragmentation pathways. Proposed fragmentation mechanisms were drawn for the experimentally observed fragmentation pathways in all the CID experiments. Cyclic oxidation products could also be discerned and their fragmentation behaviour under low energy collisional conditions was studied in detail. Gas-phase deprotonation potentials were calculated by the use of DFT B3LYP/6-311+G(2d,p)//B3LYP/6-31+G(d) + ZPVE to estimate the most thermodynamically favourable deprotonation site for efficient negative ion formation in the ion source. The optimized gas-phase geometries for the most prominent oxidation products reveal a strong intramolecular interaction between the upper and lower C4 carbon chains, which are formed after the decomposition of the primary ozonide generated by ozone attack of the reactive endocyclic C==C bond.     

CH2F自由基与HNCO反应机理的理论研究

The reaction mechanism of CH2F radical with HNCO was investigated by density functional theory (DFT)at the B3LYP/6-311++G(d,p) level. The geometries of the reactants, the intermediates, the transition states and the products were optimized. The transition states were verified through the vibration analysis.The relative energies were calculated at the QCISD(T)/6-311++G**//B3LYP/6-311++G(d,p) level. Seven feasible reaction pathways of the reaction were studied. The results indicate that the pathway (5) is the most favorable to occur, so it is the main pathway of the reaction.     

Theoretical and experimental study of the absolute configuration of helical structure of (2R,3S)-Rubiginone A2 analog

A novel analogue of (2R,3S)-Rubiginone A₂ was synthesized as a chiral helical model compound via an eight-step procedure (2.7% overall yield). Quantum methods, such as density functional theory (DFT) at different basis sets of 6-311+(d), 6-311++G(2d,p), were used to compute its optical rotation and electronic circular dichroism at the B3LYP/6-311++G(2d,p) level in the gas phase and in solution using PCM model, respectively. UV corrections were performed in electronic circular dichroism (ECD) simulations to match the experimental ECD well. The suitable computational methods, e.g., B3LYP/6-311++G(2d,p)//B3LYP/6-311++G(2d,p) in the gas phase using zero-point energy in Boltzmann statistics, were found and suggested for optical rotation and circular dichroism computations that can be used for absolute configuration determination of chiral helical compounds.     

Reaction of ketenyl radical with acetylene: a promising route for cyclopropenyl radical

The reaction of the ketenyl radical (HCCO) with acetylene (C(2)H(2)) is very relevant to the oxygen-acetylene flames and fuel-rich combustion process for nitrogen-containing compounds. Unfortunately, except for several rate constant measurements, the mechanism is completely unknown for this reaction. In this paper, detailed theoretical investigations are performed for the HCCO + C(2)H(2) reaction at the G3B3 level using the B3LYP/6-31G(d), B3LYP/6-311++G(d,p), and QCISD/6-31G(d) geometries. The exclusive fragmentation channel is the formation of the cyclopropenyl radical (c-C(3)H(3)) and carbon monoxide (CO) via the chainlike OCCHCHCH and three-membered ring OC-cCHCHCH intermediates. Thus, the mass spectroscopic peak of C(3)H(3)(+) in a previous experiment can be explained. The calculated overall reaction barrier is 4.4, 4.4, and 5.3 kcal/mol at the G3B3//B3LYP/6-31G(d), G3B3//B3LYP/6-311++G(d,p), and G3B3//QCISD/6-31G(d) levels, respectively. The title reaction may provide an effective route for generating the long-sought cyclopropenyl radical in the laboratory, which has been the long-standing subject of numerous theoretical studies as the simplest cyclic conjugate radical, and its bulky derivatives were already known. Future experimental investigations for the HCCO + C(2)H(2) reaction are greatly desired to test the predicted fragmentation channel. The implication of the present study in combustion and interstellar processes is discussed.     

The Dowd-Beckwith ring expansion: a theoretical study

UB3LYP/6-311++G(d,p) and ROMP2/6-311++G(d,p)//UB3LYP/6-311++G(d,p) calculations including the effect of benzene solvent through the PCM-UAHF method render a concerted mechanism without fragmentation as the most favourable one for the Dowd-Beckwith radical ring expansion of the bromomethyl adduct of methyl cyclopentanone-2-carboxylate to yield methyl cyclohexanone-3-carboxylate. The corresponding concerted TS is a bicyclic alcoxy radical.