Infrared spectroscopy of ozone-water complex in a neon matrix

Matrix isolation infrared spectroscopy has been applied to study an ozone-water complex of atmospheric interest. The complex was identified in the spectral region of three normal modes of ozone and water. Ab initio calculation at MP4(SDQ), QCISD, and CCSD(T) levels indicates the existence of only one stable conformer, which accords with the present experimental result. This conformer belongs to the Cs symmetry group where two molecular planes of ozone and water are perpendicular to the Cs symmetry plane. The binding energy was calculated to be 1.89 kcal/mol at the CCSD(T)/6-311++G(3df,3pd)//CCSD(T)/6-311++G(d,p) level of theory. The formation constant and atmospheric abundance of the ozone-water complex are estimated using the thermodynamic and spectroscopic data obtained.     

Valence orbital response to pseudorotation of tetrahydrofuran: a snapshot using dual space analysis

The pseudorotation of tetrahydrofuran (THF) (C(4)H(8)O) has been studied using density functional theory, with respect to the valence orbital responses to the ionization potentials and to orbital electron and momentum distributions. Three conformations of THF, the global minimum structure C(s), local minimum structure C(2), and a transition state structure C(1), which are characteristic configurations on the potential energy surface, are examined using the SAOP/et-pVQZ//B3LYP/6-311++G** models with the aforementioned dual space analysis. It is noted in the ionization energy spectra that the minimum structures C(s) and C(2) are not directly connected by pseudorotation, but through the transition state structure C(1). As a result, some orbitals of the C(s) conformer are able to "correlate" to orbitals of the C(2) conformer without a strict symmetry constraint, i.e., orbital 7a' of the C(s) conformer is correlated to orbital 5b of the C(2) conformer. It is also noted that although the valence orbital ionization potentials are not significantly altered by the pseudorotation of THF, their spectra (mainly due to excitation) are quite different indeed. Detailed orbital analysis based on dual space analysis is given. The valence orbital behavior of the conformations is orbital dependent. It can be approximately divided into three groups: the "signature group" is associated with orbitals experiencing significant changes. The frontier orbitals are in this group. The "nearly identical group" includes orbitals without apparent changes across the conformations. Most of the orbitals showing a certain degree of distortion during the pseudorotation process belong to the third group. The present study demonstrates that a comprehensive understanding of the pseudorotation of THF and its dynamics requires multidimensional information and that the information gained from momentum space is complementary to that from the more familiar coordinate space.     

Conformational analysis of the methyl ester of alanine by gas electron diffraction and Ab initio geometry optimization

The gas electron diffraction data of the methyl ester of -alanine were recorded at an approximate nozzle temperature of 80°C and analyzed with geometrical constraints taken from 4-21G ab initio gradient geometry refinements. The molecular intensities are consistent with a model for the most populated conformational state in which the N---C---C=O torsion (τ) is syn-peri-planar and a bifurcated hydrogen bond exists between NH₂ and C=O. Details of the intensity data make it impossible to rule out the presence of significant concentrations of a second conformer, (τ + 150°), in contrast to the ab initio calculations which predict greater stability of the system at τ −150° than +150°.     

Experimental and theoretical electron momentum spectroscopic study of the valence electronic structure of tetrahydrofuran under pseudorotation

The most populated structure of tetrahydrofuran (THF) has been investigated in our previous study using electron momentum spectroscopy (EMS). Because of the relatively low impact energy (600 eV) and low energy resolution (DeltaE = 1.20 eV) in the previous experiment, only the highest occupied molecular orbital (HOMO) of THF was investigated. The present study reports the most recent high-resolution EMS of THF in the valence space for the first time. The binding energy spectra of THF are measured at 1200 and 2400 eV plus the binding energies, respectively, for a series of azimuthal angles. The experimentally obtained binding energy spectra and orbital momentum distributions (MDs) are employed to study the orbital responses of the pseudorotation motion of THF. The outer valence Greens function (OVGF), the OVGF/6-311++G** model, and density function theory (DFT)-based SAOP/et-pVQZ model are employed to simulate the binding energy spectra. The orbital momentum distributions (MDs) are produced using the DFT-based B3LYP/aug-cc-pVTZ model, incorporating thermodynamic population analysis. Good agreement between theory and experiment is achieved. Orbital MDs of valence orbitals exhibit only slight differences with respect to the impact energies at 1200 and 2400 eV, indicating validation of the plane wave impulse approximation (PWIA). The present study has further discovered that the orbital MDs of the HOMO in the low-momentum region (p < 0.70 a.u) change significantly with the pseudorotation angle, phi, giving a v-shaped cross section, whereas the innermost valence orbital of THF does not vary with pseudorotation, revealing a very different bonding mechanism from the HOMO. The present study explores an innovative approach to study pseudorotation of sugar puckering, which sheds a light to study other biological systems with low energy barriers among ring-puckering conformations.     

The conformational composition of 3-buten-1-ol, the importance of intramolecular hydrogen bonding

The conformations of 3-buten-1-ol (1), and its model compounds cis-6-methyl-3-cyclohexen-1-ol (6), 3-cyclopenten-1-ol (7) and epicholesterol (9) have been investigated by FT-IR and ¹H NMR spectroscopy. The energies and geometries of 1, 6 and 7 were also investigated by molecular mechanics, semiempirical molecular orbital and ab initio calculations, while 9 was investigated by molecular mechanics only. The objective of the work was to study the conformational composition and importance of intramolecular OH…π hydrogen bonding for this composition in 1. Only two conformers of 1 have a geometrical possibility for intramolecular hydrogen bonding: Conformers 12 and 13 (Fig. 1). Compounds 6 and 7 were used as models for Conformer 12, while 9 was used as a model for Conformer 13. The investigations showed that Conformer 13 is the only hydrogen-bonded conformer, and that Conformer 12 is not intramolecularly hydrogen bonded. Conformer 13 was the most populated conformer, while Conformer 12 was hardly populated. The combination of experimental and theoretical data, and the use of model compounds was found necessary to obtain this conclusion.     

Unsaturated didehydrodeoxycytidine drugs. 1. Impact of C=C positions in the sugar ring

The chemical names of a pair of recently synthesized antitumor drugs are given in the present study as 1',2'-didehydro-3',4'-deoxycytidine and 3',4'-didehydro-2',4'-deoxycytidine. The order of stabilities, geometries, and ionization potentials of the unsaturated sugar-modified cytidine derivatives is investigated quantum mechanically. Our density functional theory calculations based on the B3LYP/6-311++G** model reveal that 3',4'-didehydro-2',4'-deoxycytidine (SD-C2) is slightly more stable than its isomer, 1',2'-didehydro-3',4'-deoxycytidine, by an energy of 5.28 kJ x mol(-1) in isolation. The isomers structurally differ by only the C=C location in the sugar ring. However, the compounds exhibit an unusual orientation with a less puckered sugar ring; that is, 3',4'-didehydro-2',4'-deoxycytidine is determined to be a beta-nucleoside, which is a C1'-endo, north conformer with an anticlinal sugar ring, whereas 1',2'-didehydro-3',4'-deoxycytidine is neither an alpha-nucleoside nor a beta-nucleoside but is a C4'-endo, south conformer with an antiperiplanar sugar ring. The present study further indicates that the C=C double bond location imposes significant effects on their ionization potentials (IPs) and other important molecular properties such as molecular electrostatic potential (MEP). In addition, inner shell binding energy spectral variations with respect to the C=C bond exhibit more site dependence. The valence shell binding energy spectral changes are, on the other hand, significant and delocalized. The latter indicates that such changes in valence space are not isolated effects but are within the entire nucleoside. Finally, the present study suggests that the nearly 0.6 eV difference in the first ionization potentials (highest occupied molecular orbital) of the isomers is sufficiently large to identify them by further spectroscopic measures.     

Infrared and Raman spectra, conformational stability, ab initio calculations of structure, and vibrational assignment of ethynylmethyl cyclobutane.

The infrared spectra (3200-300 cm(-1)) of the gas and solid and the Raman spectra (3200-30 cm(-1) of the liquid with quantitative depolarization values and the solid have been recorded of ethynylmethyl cyclobutane (cyclobutylmethyl acetylene), c-C4H7CH2C[ triple bond]CH. Both the equatorial and the axial conformers have been identified in the fluid phases and both the gauche and trans conformations of the methyl acetylenic group have been identified for each ring conformer. Variable temperature (-105 to -150 degrees C) studies of the infrared spectra of the sample dissolved in liquid krypton have been carried out. From these data the enthalpy differences have been determined to the 112 +/- 11 cm(-1) (1.34 +/- 0.13 kJ mol) between the most stable equatorial-trans (Et) conformer and the equatorial-gauche (Eg) conformer which is the second most stable conformer and 327 +/- 35 cm(-1) (3.91 + 0.42 kJ/mol) with the axial-gauche (Ag) conformer which is the least stable conformer. The enthalpy difference between the axial-trans (At) and the equatorial-gauche(Eg) is 56 +/- 6 cm(-1). At ambient temperature there is approximately 33% of the Et conformer, 38% of the Eg form, 15% of the At and 14% of the Ag conformer. For the polycrystalline solid the Eg conformer is the only form present which probably results form the packing in the crystal. A complete vibrational assignment is proposed for both equatorial conformers and additionally a few of the fundamentals of the At and Ag conformers have been assigned. The vibrational assignments are supported by normal coordinate calculations utilizing the force constants from ab initio MP2/6-31G(d) calculations. Complete equilibrium geometries have been determined for all four possible rotamers by ab initio calculations employing the 6-31G(d) and 6-311 + + G(d,p) basis sets at levels of restricted Hartree-Fock (RHF) and /or Moller-Plesset (MP2) with full electron correlation by the perturbation method to second order. The results are discussed and compared to those obtained for some similar molecules.     

Density functional theory investigations on the conformational stability and vibrational spectra of the model compound of vitamin K

Becke 3-Lee–Yang–Parr density functional theory (DFT/B3LYP) using 6-31G(d) and 6-311G(d) basis sets and the scaled quantum mechanics (SQM) force field method are used to study molecular conformations and vibrational spectra of a model compound of vitamin K (VK). In this molecule, there are six conformers on the torsional potential energy map of the dihedral angles C₈C₁₄C₁₅C₁₆ () and C₇C₈C₁₄C₁₅ (β). It is shown that the VK_1 conformer ( = 237.7° and β = 274.2°) is the most stable form. For this lowest energy conformer, the harmonic force fields calculated by B3LYP/6-311G(d) and B3LYP/6-31G(d) methods are scaled with one scale factor of 0.9623 and a set of different scale factors transferred from the previous studies for the other similar molecules, respectively. The vibrational frequencies, IR intensities, and Raman activities are obtained for the lowest energy conformer. On the basis of B3LYP calculations, the normal mode analysis is performed to assign the vibrational fundamental frequencies according to the potential energy distributions. The computational frequencies are in good agreement with the observed results.     

2,4,6-三甲氧基苯-1-O-D葡萄糖苷的核磁共振谱的理论研究

The 1H and 13C-NMR of 2,4,6-trimethoxyphenol-1-O-D-glucopyranoside(Compound 1) isolated from Celastrus angulatus (Celastraceae) was calculated theoretically at the both levels HF/6-311+G(2d,p)//B3LYP/6-31G(d) and HF/6-311+G(2d,p)//B3LYP/6-31G(d,p) using the GIAO (gauge-independent atomic orbital) method. Statistical error analysis for theoretically predicted δH and δC values versus those experimentally observed for compound 1 was discussed. The results show that the theoretically predicted δH and δC values of β conformer of compound 1 are more close to the experimentally observed values than α conformer, and the β conformer of compound 1 is more stable than α conformer according to molecular energy theoretically calculated. So compound 1 is assigned to be 2,4,6-trimethoxyphenol-1-O-β-D-glucopyranoside, which is in good consistence with the conclusion deduced by the anomeric proton signal (δH=4.80, J=7.3 Hz) experimentally observed.     

A conformational analysis of the six isomers of thiobispyridine

The conformational behaviour of the six isomers of thiobispyridine has been investigated using ab initio STO-3G*//rigid-roto, STO-3G*//STO-3G* and 6–31G**//STO-3G* molecular orbital models. The analysis reveals both the importance of optimising critical structure parameters and the basis set dependence of calculated rotational barrier heights. The most reliable model (6–31G**//STO-3G*) clearly indicates that the minimum energy conformers are not planar and that energy barriers between 30–100 kJ mol^?1 restrict inter-conversion to planar structures, thereby preventing conjugation between the p-electrons of the sulfur atom and the π system of both pyridine rings. From the calculated barrier heights, two mechanisms can be employed to explain conformer interconversion about the C? S bond: a disrotatory one-ring flip or a conrotatory two-ring flip mechanism. Where comparisons can be made (eg. 2,2′-thiobispyridine), dipole moment calculations are shown to be in good agreement with experiment. Finally, of the six isomers, appropriately substituted 2,2′, 2,3′- and 2,4′-thiobispyridines are most prone to a Smiles rearrangment.     

Probing furanose ring conformation by gas-phase computational methods: energy profile and structural parameters in methyl beta-D-arabinofuranoside as a function of ring conformation

The potential energy surface of methyl beta-D-arabinofuranoside (3) has been studied by ab initio molecular orbital (HF/6-31G) and density functional theory (B3LYP/6-31G) calculations via minimization of the 10 possible envelope conformers. The partial potential energy surface identified that the global minimum and lowest energy northern conformer was E(2). In the HF calculations, (2)E was the most stable southern conformer, while the density functional theory methods identified (4)E as the local minimum in this hemisphere. Additional calculations at higher levels of theory showed that the B3LYP-derived energies of many of the envelope conformers of 3 are dependent upon the basis set used. It has also been demonstrated that B3LYP/6-31+G//B3LYP/6-31G single point energies are essentially the same as those obtained from full geometry optimizations at the B3LYP/6-31+G level. The northern and southern minima of the B3LYP/6-31+G surface are, respectively, the E(2) and (2)E conformers. The B3LYP/6-31G geometries were used to study the relationship between ring conformation and various structural parameters including bond angles, dihedral angles, bond lengths, and interatomic distances.     

Conformations of alpha-aminobutyric acid in the gas phase.

A laser-ablation molecular-beam Fourier transform microwave (LA-MB-FTMW) spectrometer has been successfully applied to the structural study of alpha-aminobutyric acid. Three neutral conformers have been identified in the gas phase by comparing their experimental rotational and 14N nuclear quadrupole coupling parameters with those predicted by ab initio calculations at the MP2/6-311++G(d,p) level. The most stable conformer is stabilized by a bifurcated amine-to-carbonyl hydrogen bond (N--HO=C) and a cis-COOH group, and the side-chain adopts a configuration with a torsion angle tau(C(gamma)-C(beta)-C(alpha)-C') of about 180 degrees. The second most stable conformer exhibits the same configuration for the amino acid skeleton but adopts a different orientation for the side chain with tau(C(gamma)-C(beta)-C(alpha)-C') approximately -60 degrees. In the third conformer an intramolecular hydrogen bond is established between the hydroxyl group and the nitrogen atom (NH--O), with a side-chain orientation similar to that of the most stable conformer.     

Gas-phase structure of N,N-dimethylglycine.

Three conformers of the neutral amino acid N,N-dimethylglycine [(CH3)2NCH2COOH] were detected in a supersonic expansion by a combination of laser ablation (LA) and molecular-beam Fourier transform microwave (MB-FTMW) spectroscopy. A bifurcated methyl-to-carbonyl (C--HO==C) weak intramolecular hydrogen bond might stabilise the most stable conformer of C(s) symmetry. The second most stable conformer of C1 symmetry has a hydrogen bond between the hydroxyl group and the lone pair at the nitrogen atom (NH--O). The r(s) and r0 structures were derived for this conformer from the rotational data for the parent and six minor 13C, 15N and OD isotopomers. A third conformer exhibits a cis-carboxyl functional group and C1 symmetry. Ab initio MP2/6-311++G(d,p) predictions of the spectroscopic parameters were useful in analysing the spectra. In particular, the agreement of the nuclear quadrupole coupling constants with those calculated was conclusive in identifying the different conformers.     

Conformational properties of penicillins: quantum chemical calculations and molecular dynamics simulations of benzylpenicillin

Herein, we present theoretical results on the conformational properties of benzylpenicillin, which are characterized by means of quantum chemical calculations (MP2/6-31G* and B3LYP/6-31G*) and classical molecular dynamics simulations (5 ns) both in the gas phase and in aqueous solution. In the gas phase, the benzylpenicillin conformer in which the thiazolidine ring has the carboxylate group oriented axially is the most favored one. Both intramolecular CH. O and dispersion interactions contribute to stabilize the axial conformer with respect to the equatorial one. In aqueous solution, a molecular dynamics simulation predicts a relative population of the axial:equatorial conformers of 0.70:0.30 in consonance with NMR experimental data. Overall, the quantum chemical calculations as well as the simulations give insight into substituent effects, the conformational dynamics of benzylpenicillin, the frequency of ring-puckering motions, and the correlation of side chain and ring-puckering motions.     

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.     

The role of chalcogen-chalcogen interactions in the intrinsic basicity acidity of beta-chalcogenovinyl(thio)aldehydes HC([double bond]X)[bond]CH[double bond]CH[bond]CYH

The intrinsic acidity and basicity of a series of beta-chalcogenovinyl(thio)aldehydes HC([double bond]X)[bond]CH[double bond]CH[bond]CYH (X=O, S; Y=Se, Te) were investigated by B3LYP/6-311+G(3df,2p) density functional and G2(MP2) calculations on geometries optimized at the B3LYP/6-31G(d) level for neutral molecules and at the B3LYP/6-31+G(d) level for anions. The results showed that selenovinylaldehyde and selenovinylthioaldehyde should behave as Se bases in the gas phase, because the most stable neutral conformer is stabilized by an X[bond]H...Se (X=O, S) intramolecular hydrogen bond (IHB). In contrast the Te-containing analogues behave as oxygen or sulfur bases, because the most stable conformer is stabilized by typical X...Y[bond]H chalcogen-chalcogen interactions. These compounds have a lower basicity than expected because either chalcogen-chalcogen interactions or IHBs become weaker upon protonation. Similarly, they are also weaker acids than expected because deprotonation results in a significantly destabilized anion. Loss of the proton from the X[bond]H or Y[bond]H groups is a much more favorable than from the C[bond]H groups. Therefore, for Se compounds the deprotonation process results in loss of the X[bond]H...Se (X=O, S) IHBs present in the most stable neutral conformer, while for Te-containing compounds the stabilizing X...Y[bond]H chalcogen-chalcogen interaction present in the most stable neutral conformer becomes repulsive in the corresponding anion.     

Theoretical conformational analysis for neurotransmitters in the gas phase and in aqueous solution. Norepinephrine

The natural neurotransmitter (R)-norepinephrine takes the monocationic form in 93% abundance at the physiological tissue pH of 7.4. Ab initio and DFT/B3LYP calculations were performed for 12 protonated conformers of (R)-norepinephrine in the gas phase with geometry optimizations up to the MP2/6-311++G level, and with single-point calculations up to the QCISD(T) level at the HF/6-31G-optimized geometries. Four monohydrates were studied at the MP2/6-31G//HF/6-31G level. In the gas phase, the G1 conformer is the most stable with phenyl.NH(3)(+) gauche and HO(alc).NH(3)(+) gauche arrangements. A strained intramolecular hydrogen bond was found for conformers (G1 and T) with close NH(3)(+) and OH groups. Upon rotation of the NH(3)(+) group as a whole unit about the C(beta)-C(alpha) axis, a 3-fold potential was calculated with free energies for barriers of 3-12 kcal/mol at the HF/6-31G level. Only small deviations were found in MP2/6-311++G single-point calculations. A 2-fold potential was calculated for the phenyl rotation with free energies of 11-13 kcal/mol for the barriers at T = 310 K and p = 1 atm. A molecular mechanics docking study of (R)-norepinephrine in a model binding pocket of the beta-adrenergic receptor shows that the ligand takes a conformation close to the T(3) arrangement. The effect of aqueous solvation was considered by the free energy perturbation method implemented in Monte Carlo simulations. There are 4-5 strongly bound water molecules in hydrogen bonds to the conformers. Although hydration stabilizes mostly the G2 form with gauche phenyl.NH(3)(+) arrangement and a water-exposed NH(3)(+) group, the conformer population becomes T > G1 > G2, in agreement with the PMR spectroscopy measurements by Solmajer et al. (Z. Naturforsch. 1983, 38c, 758). Solvent effects reduce the free energies for barriers to 3-6 and 9-12 kcal/mol for rotations about the C(beta)-C(alpha) and the C(1)(ring)-C(beta) axes, respectively.     

Ab initio study of the substituent effects on the relative stability of the E and Z conformers of phenyl esters. Stereoelectronic effects on the reactivity of the carbonyl group

Equilibria between the Z (tau1= 0 degrees) and E (tau1= 180 degrees) conformers of p-substituted phenyl acetates 4 and trifluoroacetates 5 (X = OMe, Me, H, Cl, CN, NO2) were studied by ab initio calculations at the HF/6-31G* and MP2/6-31G* levels of theory. The preference for the Z conformer, DeltaE(HF), was calculated to be 5.36 kcal mol(-1) and 7.50 kcal mol(-1) for phenyl acetate and phenyl trifluoroacetate (i.e., with X = H), respectively. The increasing electron-withdrawing ability of the phenyl substituent X increases the preference of the Z conformer. An excellent correlation with a negative slope was observed for both series between DeltaE of the E-Z equilibrium and the Hammett sigma constant. By using an appropriate isodesmic reaction, it was shown that electron-withdrawing substituents decrease the stability of both conformers, but the effect is higher with the E conformer. Electron-withdrawing phenyl substituents decrease the delocalization of the lone pair of the ether oxygen to the C=O antibonding orbital (nO--> pi*C=O) in both the E and Z forms and in both series studied; this effect is higher in the E conformer than in the Z conformer. The nO --> pi*C=O electron donation has a minimum value with tau1= 90 degrees and a maximum value with tau1= 0 degrees (the Z conformer), the value with tau1= 180 degrees (the E conformer) being between these two values, obviously due to steric hindrance. The effects of the phenyl substituents on the reactivity of the esters studied are discussed in terms of molecular orbital interactions. ED/EW substituents adjust the availability of the pi*C=O antibonding orbital to interact with the lone pair orbital of the attacking nucleophile and therefore affect the reactivity: EW substituents increase and ED substituents decrease it. Excellent correlations were observed between the rate coefficients of nucleophilic acyl substitutions and pi*C=O occupancies of the ester series 4 and 5.     

Topological study of intramolecular hydrogen bonding in beta-hydroxyethylperoxy radical and beta-hydroxyethoxy radical along its dissociation pathway

The presence of intramolecular hydrogen bond (IHB) in beta-hydroxyethylperoxy and beta-hydroxyethoxy radicals was investigated using the QTAIM topological study on UB3LYP/6-311++G(2d,2p) charge densities. Only one of the two conformers of beta-hydroxyethylperoxy radical which were previously considered to present IHB displays a bond critical point (BCP) associated to an IHB. Furthermore, the atomic energies and electron populations indicate no evidence of IHB in the second conformer. Nevertheless, very small differences in molecular energies were obtained using several one-step and multi-step methods (G3, G3B3) between both conformers. No BCP is found between the hydroxyl hydrogen and the oxygen in the most stable conformer or in the transition state for the dissociation path of beta-hydroxyethoxy radical. However, a BCP is formed in the last steps of this path, thereby yielding H-bonded products.     

吡啶-BH~3相互作用复合物的理论研究

对吡啶-BH~3复合物分别用MP2/6-31+G^*和B3LYP/6-31+G^*进行理论计算以预测该复合物的构型及解离能，得到四种构型，在MP2优化构型基础上作CCSD/6-31+G^*单点能量计算以验证MP2与B3LYP结果的可靠性，然后用B3LYP作振动频率分析，计算了各构型的垂直电离势，最后用更大基组作单点能量计算和自然键轨道(NBO)分析。结果表明，N-B直接相连的构型最稳定，其解离能为141.50kJ/mol,MP2和B3LYP对N-H接近的构型结果相关较大，另外两种构型稳定性介于二者之间，解离能分别为15.18kJ/mol,14.06kJ/mol(MP2/6-31+G^*)。