Dipole effects on cation-pi interactions: absolute bond dissociation energies of complexes of alkali metal cations to N-methylaniline and N,N-dimethylaniline

Threshold collision-induced dissociation of M (+)( nMA) x with Xe is studied using guided ion beam mass spectrometry, where nMA = N-methylaniline and N, N-dimethylaniline and x = 1 and 2. M (+) includes the following alkali metal cations: Li (+), Na (+), K (+), Rb (+), and Cs (+). In all cases, the primary dissociation pathway corresponds to the endothermic loss of an intact nMA ligand. The primary cross section thresholds are interpreted to yield 0 and 298 K bond dissociation energies (BDEs) for ( nMA) x-1 M (+)-( nMA) after accounting for the effects of multiple ion-neutral collisions, the internal and kinetic energy distributions of the reactants, and the dissociation lifetimes. Density functional theory calculations at the B3LYP/6-31G* level of theory are used to determine the structures of these complexes, which are also used in single-point calculations at the MP2(full)/6-311+G(2d,2p) level to determine theoretical BDEs. The results of these studies are compared to previous studies of the analogous M (+)(aniline) x complexes to examine the effects of methylation of the amino group on the binding interactions. Comparisons are also made to a wide variety of cation-pi complexes previously studied to elucidate the contributions that ion-dipole, ion-induced-dipole, and ion-quadrupole interactions make to the overall binding.     

Cation-pi interactions with a model for the side chain of tryptophan: structures and absolute binding energies of alkali metal cation-indole complexes

Threshold collision-induced dissociation techniques are employed to determine bond dissociation energies (BDEs) of mono- and bis-complexes of alkali metal cations, Li+, Na+, K+, Rb+, and Cs+, with indole, C8H7N. The primary and lowest energy dissociation pathway in all cases is endothermic loss of an intact indole ligand. Sequential loss of a second indole ligand is observed at elevated energies for the bis-complexes. Density functional theory calculations at the B3LYP/6-31G level of theory are used to determine the structures, vibrational frequencies, and rotational constants of these complexes. Theoretical BDEs are determined from single point energy calculations at the MP2(full)/6-311+G(2d,2p) level using the B3LYP/6-31G* geometries. The agreement between theory and experiment is very good for all complexes except Li+ (C8H7N), where theory underestimates the strength of the binding. The trends in the BDEs of these alkali metal cation-indole complexes are compared with the analogous benzene and naphthalene complexes to examine the influence of the extended pi network and heteroatom on the strength of the cation-pi interaction. The Na+ and K+ binding affinities of benzene, phenol, and indole are also compared to those of the aromatic amino acids, phenylalanine, tyrosine, and tryptophan to elucidate the factors that contribute to the binding in complexes to the aromatic amino acids. The nature of the binding and trends in the BDEs of cation-pi complexes between alkali metal cations and benzene, phenol, and indole are examined to help understand nature's preference for engaging tryptophan over phenylalanine and tyrosine in cation-pi interactions in biological systems.     

Experimental and theoretical investigation of alkali metal cation interactions with hydroxyl side-chain amino acids

Absolute bond dissociation energies of serine (Ser) and threonine (Thr) to alkali metal cations are determined experimentally by threshold collision-induced dissociation of M+AA complexes, where M+=Li+, Na+, and K+ and AA=Ser and Thr, with xenon in a guided ion beam tandem mass spectrometer. Experimental results show that the binding energies of both amino acids to the alkali metal cations are very similar to one another and follow the order of Li+>Na+>K+. Quantum chemical calculations at three different levels, B3LYP, B3P86, and MP2(full), using the 6-311+G(2d,2p) basis set with geometries and zero-point energies calculated at the B3LYP/6-311+G(d,p) level show good agreement with the experimental bond energies. Theoretical calculations show that all M+AA complexes have charge-solvated structures (nonzwitterionic) with [CO, N, O] tridentate coordination.     

A simple model for metal cation-phosphate interactions in nucleic acids in the gas phase: Alkali metal cations and trimethyl phosphate

Threshold collision-induced dissociation techniques are employed to determine the bond dissociation energies (BDEs) of complexes of alkali metal cations to trimethyl phosphate, TMP. Endothermic loss of the intact TMP ligand is the only dissociation pathway observed for all complexes. Theoretical calculations at the B3LYP/6-31G* level of theory are used to determine the structures, vibrational frequencies, and rotational constants of neutral TMP and the M+(TMP) complexes. Theoretical BDEs are determined from single point energy calculations at the B3LYP/6-311+G(2d,2p) level using the B3LYP/6-31G* optimized geometries. The agreement between theory and experiment is reasonably good for all complexes except Li+(TMP). The absolute M+-(TMP) BDEs are found to decrease monotonically as the size of the alkali metal cation increases. No activated dissociation was observed for alkali metal cation binding to TMP. The binding of alkali metal cations to TMP is compared with that to acetone and methanol.     

Noncovalent interactions of Cu+ with N-donor ligands (pyridine, 4,4-dipyridyl, 2,2-dipyridyl, and 1,10-phenanthroline): collision-induced dissociation and theoretical studies

Collision-induced dissociation of complexes of Cu+ bound to a variety of N-donor ligands (N-L) with Xe is studied using guided ion beam tandem mass spectrometry. The N-L ligands examined include pyridine, 4,4-dipyridyl, 2,2-dipyridyl, and 1,10-phenanthroline. In all cases, the primary and lowest-energy dissociation channel observed corresponds to the endothermic loss of a single intact N-L ligand. Sequential dissociation of additional N-L ligands is observed at elevated energies for the pyridine and 4,4-dipyridyl complexes containing more than one ligand. Ligand exchange processes to produce Cu+Xe are also observed as minor reaction pathways in several systems. The primary cross section thresholds are interpreted to yield 0 and 298 K bond dissociation energies (BDEs) after accounting for the effects of multiple ion-neutral collisions, the kinetic and internal energy distributions of the reactants, and dissociation lifetimes. Density functional theory calculations at the B3LYP/6-31G* level are performed to obtain model structures, vibrational frequencies, and rotational constants for the neutral N-L ligands and the Cu+(N-L)x complexes. The relative stabilities of the various conformations of these N-L ligands and Cu+(N-L)x complexes as well as theoretical BDEs are determined from single point energy calculations at the B3LYP/6-311+G(2d,2p) level of theory using B3LYP/6-31G* optimized geometries. Excellent agreement between theory and experiment is observed for all complexes containing one or two N-L ligands, while theory systematically underestimates the strength of binding for complexes containing more than two N-L ligands. The ground-state structures of these complexes and the trends in the sequential BDEs are explained in terms of stabilization gained from sd-hybridization and repulsive ligand-ligand interactions. The nature of the binding interactions in the Cu+(N-L)x complexes are examined via natural bond orbital analyses.     

Absolute thermodynamic measurements of alkali metal cation interactions with a simple dipeptide and tripeptide

Absolute bond dissociation energies (BDEs) of glycylglycine (GG) and glycylglycylglycine (GGG) to sodium and potassium cations and sequential bond energies of glycine (G) in Na+G2 were determined experimentally by threshold collision-induced dissociation (TCID) in a guided ion beam tandem mass spectrometer. Experimental results showed that the binding energies follow the order of Na+ > K+ and M+GGG > M+GG > M+G. Theoretical calculations at the B3LYP/6-311+G(d) level show that all complexes had charge-solvated structures (nonzwitterionic) with either [CO,CO] bidentate or [N,CO,CO] tridentate coordination for M+GG complexes, [CO,CO,CO] tridentate or [N,CO,CO,CO] tetradentate coordination for M+GGG complexes, and [N,CO,N,CO] tetradentate coordination for Na+G2. Ab initio calculations at three different levels of theory (B3LYP, B3P86, and MP2(full) using the 6-311+G(2d,2p) basis set with geometries and zero-point energies calculated at the B3LYP/6-311+G(d) level) show good agreement with the experimental bond energies. This study demonstrates for the first time that TCID measurements of absolute BDEs can be successfully extended to biological molecules as complex as a tripeptide.     

Solvation of copper ions by acetone. structures and sequential binding energies of Cu+(acetone)x,x = 1-4 from collision-induced dissociation and theoretical studies

Collision-induced dissociation of Cu+(acetone)(x), x = 1-4, with Xe is studied as a function of kinetic energy using guided ion beam mass spectrometry. In all cases, the primary and lowest energy dissociation channel observed is endothermic loss of one acetone molecule. The primary cross section thresholds are interpreted to yield 0 and 298 K bond energies after accounting for the effects of multiple ion-neutral collisions, internal energy of the complexes, and dissociation lifetimes. Density functional calculations at the B3LYP/6-31G* level of theory are used to determine the structures of these complexes and provide molecular constants necessary for the thermodynamic analysis of the experimental data. Theoretical bond dissociation energies are determined from single point calculations at the B3LYP/6-311+G(2d,2p) and MP2(full)/6-311+G(2d,2p) levels, using the B3LYP/6-31G* optimized geometries. The experimental bond energies determined here are in good agreement with previous experimental measurements made in a high-pressure mass spectrometer for the sum of the first and second bond energy (i.e., Cu+(acetone)2 --> Cu+ + 2 acetone) when these results are properly anchored. The agreement between theory and experiment is reasonable in all cases, but varies both with the size of the cluster and the level of theory employed. B3LYP does an excellent job for the x = 1 and 3 clusters, but is systematically low for the x = 2 and 4 clusters such that the overall trends in sequential binding energies are not parallel. In contrast, all MP2 values are somewhat low, but the overall trends parallel the measured values for all clusters. The trends in the measured Cu+(acetone), binding energies are explained in terms of 4s-3d sigma hybridization effects and ligand-ligand repulsion in the clusters.     

Cation-pi interactions: structures and energetics of complexation of Na+ and K+ with the aromatic amino acids, phenylalanine, tyrosine, and tryptophan

Threshold collision-induced dissociation of M(+)(AAA) with Xe is studied using guided ion beam tandem mass spectrometry. M(+) include the alkali metal ions Na(+) and K(+). The three aromatic amino acids are examined, AAA = phenylalanine, tyrosine, or tryptophan. In all cases, endothermic loss of the intact aromatic amino acid is the dominant reaction pathway. The threshold regions of the cross sections are interpreted to extract 0 and 298 K bond dissociation energies for the M(+)-AAA complexes after accounting for the effects of multiple ion-neutral collisions, internal energy of the reactant ions, and dissociation lifetimes. Density functional theory calculations at the B3LYP/6-31G level of theory are used to determine the structures of the neutral aromatic amino acids and their complexes to Na(+) and K(+) and to provide molecular constants required for the thermochemical analysis of the experimental data. Theoretical bond dissociation energies are determined from single-point energy calculations at the B3LYP/6-311++G(3df,3pd) level using the B3LYP/6-31G geometries. Good agreement between theory and experiment is found for all systems. The present results are compared to earlier studies of these systems performed via kinetic and equilibrium methods. The present results are also compared to the analogous Na(+) and K(+) complexes to glycine, benzene, phenol, and indole to elucidate the relative contributions that each of the functional components of these aromatic amino acids make to the overall binding in these complexes.     

An evaluation of various computational methods for the treatment of organoselenium compounds

A reliable computational method for the prediction of organoselenium geometries and bond dissociation energies (BDEs) has been determined on the basis of the performance of density functional theory (DFT: B3LYP and B3PW91) and ab initio molecular orbital procedures (Hartree-Fock (HF)) in conjunction with various Pople basis sets including (but not limited to) the 6-31G(d), 6-31G(d,p), 6-311G(d), 6-311G(d,p), 6-311G(2df,p), and 6-311G(3df,3pd) sets. Predicted geometries and BDEs are compared with available experimental data and quadratic configuration interaction including single and double substitutions (QCISD) results. The B3PW91/6-311G(2df,p) level of theory is recommended for the prediction of the geometries and energetics of organoselenium compounds.     

Modeling metal cation-phosphate interactions in nucleic acids in the gas phase via alkali metal cation-triethyl phosphate complexes

Threshold collision-induced dissociation techniques are employed to determine the bond dissociation energies (BDEs) of complexes of alkali metal cations, Na+, K+, Rb+, and Cs+, to triethyl phosphate (TEP). The primary and lowest energy dissociation pathway in all cases is the endothermic loss of the neutral TEP ligand. Theoretical electronic structure calculations at the B3LYP/6-311+G(2d,2p)//B3LYP/6-31G* level of theory are used to determine the structures, molecular parameters, and theoretical estimates for the BDEs of these complexes. For the complexes to Rb+ and Cs+, theoretical calculations were performed using hybrid basis sets in which the effective core potentials and valence basis sets of Hay and Wadt were used to describe the alkali metal cation, while the standard basis sets were used for all other atoms. The agreement between theory and experiment is excellent for the complexes to Na+ and K+ and is somewhat less satisfactory for the complexes to the heavier alkali metal cations, Rb+ and Cs+, where effective core potentials were used to describe the cation. The trends in the binding energies are examined. The binding of alkali metal cations to triethyl phosphate is compared with that to trimethylphosphate.     

The special five-membered ring of proline: An experimental and theoretical investigation of alkali metal cation interactions with proline and its four- and six-membered ring analogues

The interaction of the alkali metal cations, Li+, Na+, and K+, with the amino acid proline (Pro) and its four- and six-membered ring analogues, azetidine-2-carboxylic acid (Aze) and pipecolic acid (Pip), are examined in detail. Experimentally, threshold collision-induced dissociation of the M+(L) complexes, where M = Li, Na, and K and L = Pro, Aze, and Pip, with Xe are studied using a guided ion beam tandem mass spectrometer. From analysis of the kinetic energy dependent cross sections, M(+)-L bond dissociation energies are measured. These analyses account for unimolecular decay rates, internal energy of reactant ions, and multiple ion-molecule collisions. Ab initio calculations for a number of geometric conformations of the M+(L) complexes were determined at the B3LYP/6-311G(d,p) level with single-point energies calculated at MP2(full), B3LYP, and B3P86 levels using a 6-311+G(2d,2p) basis set. Theoretical bond energies show good agreement with the experimental bond energies, which establishes that the zwitterionic form of the alkali metal cation/amino acid, the lowest energy conformation, is formed in all cases. Despite the increased conformational mobility in the Pip systems, the Li+, Na+, and K+ complexes of Pro show higher binding energies. A meticulous examination of the zwitterionic structures of these complexes provides an explanation for the stability of the five-membered ring complexes.     

Bond dissociation energies and equilibrium structures of Cu+(MeOH)x, x = 1-6, in the gas phase: competition between solvation of the metal ion and hydrogen-bonding interactions

The solvation of Cu+ by methanol (MeOH) was studied via examination of the kinetic energy dependence of the collision-induced dissociation of Cu+(MeOH)x complexes, where x = 1-6, with Xe in a guided ion beam tandem mass spectrometer. In all cases, the primary and lowest-energy dissociation channel observed is the endothermic loss of a single MeOH molecule. The primary cross section thresholds are interpreted to yield 0 and 298 K bond dissociation energies (BDEs) after accounting for the effects of multiple ion-neutral collisions, kinetic and internal energy distributions of the reactants, and lifetimes for dissociation. Density functional theory calculations at the B3LYP/6-31G* level are performed to obtain model structures, vibrational frequencies, and rotational constants for the Cu+(MeOH)x complexes and their dissociation products. The relative stabilities of various conformations and theoretical BDEs are determined from single-point energy calculations at the B3LYP/6-311+G(2d,2p) level of theory using B3LYP/6-31G*-optimized geometries. The relative stabilities of the various conformations of the Cu+(MeOH)x complexes and the trends in the sequential BDEs are explained in terms of stabilization gained from sd hybridization, hydrogen-bonding interactions, electron donor-acceptor natural bond orbital stabilizing interactions, and destabilization arising from ligand-ligand repulsion.     

DFT study of the C-Cl bond dissociation enthalpies and electronic structure of substituted chlorobenzene compounds

Quantum chemical calculations were used to estimate the bond dissociation energies (BDEs) for 13 substituted chlorobenzene compounds. These compounds were studied by employing the hybrid density functional theory methods (B3LYP, B3PW1, B3P86) with 6-31G** and 6-311G** basis sets. It was demonstrated that B3P86/6-311G** method is the best method for computing the reliable BDEs for substituted chlorobenzene compounds which contain the C-Cl bond. It was found that the C-Cl BDE depends strongly on a computational method and basis set used. Substitution effect on the C-Cl BDE of substituted chlorobenzene compounds is further discussed. It is shown that the effects of substitution on the C-Cl BDE of substituted chlorobenzene compounds are very insignificant. Frontier orbital energy gap of studied compounds was also investigated. From the data on frontier orbital energies gap, we estimated the relative thermal stability of substituted chlorobenzene compounds.     

Experimental and theoretical studies of potassium cation interactions with the acidic amino acids and their amide derivatives

The binding of K(+) to aspartic acid (Asp), glutamic acid (Glu), asparagine (Asn), and glutamine (Gln) is examined in detail by studying the collision-induced dissociation (CID) of the four potassium cation-bound amino acid complexes with Xe using a guided ion beam tandem mass spectrometer (GIBMS). Formed by electrospray ionization, these complexes have energy-dependent CID cross sections that are analyzed to provide 0 K bond energies after accounting for unimolecular decay rates, internal energy of reactant ions, and multiple ion-molecule collisions. Quantum chemical calculations for a number of geometric conformations of each K(+)(L) complex are determined at the B3LYP/6-311+G(d,p) level with single-point energies calculated at B3LYP, B3P86, and MP2(full) levels using a 6-311+G(2d,2p) basis set. Theoretical bond dissociation energies are in good agreement with the experimental values. This coordinated examination of both experimental work and quantum chemical calculations allows for a comprehensive understanding of the molecular interactions of K(+) with the Asx and Glx amino acids. K(+) binding affinities for the amide complexes are systematically stronger than those for the acid complexes by 9+/-1 kJ/mol, which is attributed to an inductive effect of the OH group in the carboxylic acid side chain. Additionally, the K(+) binding affinity for the longer-chain amino acids (Glx) is enhanced by 5+/-1 kJ/mol compared to the shorter-chain Asx because steric effects are reduced. Further, a detailed comparison between experimental and theoretical results reveals interesting differences in the binding of K(+) and Na(+) to these amino acids.     

Experimental and theoretical investigation of the decomposition of lithiated hydroxyl side-chain amino acids

Lithium cation complexes with serine (Ser) and threonine (Thr) are collisionally activated with xenon in a guided ion beam tandem mass spectrometer and are observed to exhibit a variety of decomposition pathways in addition to a loss of the intact ligand. Prominent pathways include a loss of H2O, CO2, and aldehydes (XCHO where X=H for Ser and CH3 for Thr). Quantum chemical calculations at the B3LYP/6-311+G(d,p) level are used to explore the reaction mechanisms for these processes in detail. Complete potential energy surfaces for all three processes are elucidated, including all intermediates and transition states. Theoretical molecular parameters for the rate-limiting transition states are then used to analyze the threshold energies in the experimental data, providing experimental measurements of the energies of these transition states. These experimental energies are compared with single-point energies calculated at three different levels, B3LYP, B3P86, and MP2(full), using the 6-311+G(2d,2p) basis set with geometries and zero-point energies calculated at the B3LYP/6-311+G(d,p) level. Good agreement between experiment and theory (especially MP2(full)) suggests that the reaction mechanisms have been reasonably elucidated.     

Structural and Energetic Effects in the Molecular Recognition of Acetylated Amino Acids by 18-Crown-6

Absolute 18-crown-6 (18C6) binding affinities of four protonated acetylated amino acids (AcAAs) are determined using guided ion beam tandem mass spectrometry techniques. The AcAAs examined in this work include: N-terminal acetylated lysine (N_???CAcLys), histidine (N_???CAcHis), and arginine (N_???CAcArg) as well as side chain acetylated lysine (N_???CAcLys). The kinetic-energy-dependent cross sections for collision-induced dissociation (CID) of the (AcAA)H⁺(18C6) complexes are analyzed using an empirical threshold law to extract absolute 0 and 298?K (AcAA)H⁺?18C6 bond dissociation energies (BDEs) after accounting for the effects of multiple collisions, kinetic and internal energy distributions of the reactants, and unimolecular dissociation lifetimes. Theoretical electronic structure calculations are performed to determine stable geometries and energetics for neutral and protonated 18C6 and the AcAAs as well as the proton bound complexes of these species, (AcAA)H⁺(18C6), at the B3LYP/6-311+G(2d,2p)//B3LYP/6-31?G* and M06/6-311+G(2d,2p)//B3LYP/6-31G* levels of theory. For all four (AcAA)H⁺(18C6) complexes, loss of neutral 18C6 corresponds to the most favorable dissociation pathway. At elevated energies, products arising from sequential dissociation of the primary CID product, H⁺(AcAA), are also observed. Protonated N_???CAcLys exhibits a greater 18C6 binding affinity than other protonated N_???CAcAAs, suggesting that the side chains of Lys residues are the preferred binding sites for 18C6 complexation to peptides and proteins. N_???CAcLys exhibits a greater 18C6 binding affinity than N_???CAcLys, suggesting that binding of 18C6 to the side chain of Lys residues is more favorable than to the N-terminal amino group of Lys.     

Theoretical study of the N-H…O red-shifted and blue-shifted hydrogen bonds

Theoretical calculations are performed to study the nature of the hydrogen bonds in complexes HCHO…HNO, HCOOH…HNO, HCHO…NH3, HCOOH…NH3, HCHO…NH2F and HCOOH…NH2F. The geometric structures and vibrational frequencies of these six complexes at the MP2/6-31+G(d,p), MP2/6-311++G(d,p), B3LYP/6-31+G(d,p) and B3LYP/6-311++G(d,p) levels are calculated by standard and counterpoise-corrected methods, respectively. The results indicate that in complexes HCHO…HNO and HCOOH…HNO the N-H bond is strongly contracted and N-H…O blue-shifted hydrogen bonds are observed. While in complexes HCHO…NH3, HCOOH…NH3, HCHO…NH2F and HCOOH…NH2F, the N-H bond is elongated and N-H…O red-shifted hydrogen bonds are found. From the natural bond orbital analysis it can be seen that the X-H bond length in the X-H…Y hydrogen bond is controlled by a balance of four main factors in the opposite directions: hyperconjugation, electron density redistribution, rehybridization and structural reorganization. Among them hyperconjugation has the effect of elongating the X-H bond, and the other three factors belong to the bond shortening effects. In complexes HCHO…HNO and HCOOH…HNO, the shortening effects dominate which lead to the blue shift of the N-H stretching frequencies. In complexes HCHO…NH3, HCOOH…NH3, HCHO…NH2F and HCOOH…NH2F where elongating effects are dominant, the N-H…O hydrogen bonds are red-shifted.     

气相中O3与HSO自由基间的氢键复合物

气相中O3与HSO自由基之间的相互作用及其反应在大气化学中非常重要.在DFT-B3LYP/6-311++G**和MP2/6-311++G**水平上求得O3+HSO复合物势能面上的稳定构型,B3LYP方法得到了三种构型(复合物Ⅰ,Ⅱ和Ⅲ),而MP2方法只能得到一种构犁(复合物Ⅱ).在复合物Ⅰ和Ⅲ中,HSO单元中的1H原子作为质子供体.与O3分子中的端基O原子作为质子受体相互作用,形成红移氢键复合物;而在复合物Ⅱ中,虽与复合物Ⅰ和Ⅲ中具有相间的质子供体和质子受体,却形成了蓝移氢键复合物.B3LYP/6-311++G**水平上计算的单体间相互作用能的计算考虑了基组重甍误差(BSSE)和零点振动能(ZPVE)校正,其值在-3.37到-4.55 kJ·mol-1之间.采用自然键轨道理论(NBO)对单体间相互作用的本质进行了考查,并通过分子中原子理论(AIM)分析了三种复合物中氢键的电子密度拓扑性质.     

CH and Chalogen bond dissociation energies for fluorinated and chlorinated methane evaluated with hybrid B3LYP density functional theory methods and their comparison with experimental data and the CBS-Q ab initio computational approach

The CH and Chalogen bond dissociation energies (BDEs) were computed with the hybrid B3LYP/6-311 + G(2d,2p) theory model for chlorinated and fluorinated methane. All computed values were substantially lower (5–10 kcal mol⁻¹) than the experimental values. To obtain better agreement, a correlation factor was introduced. When this factor was applied, excellent agreement between the B3LYP/6-311 + G(2d,2p) computed energies and the experimental BDEs was observed. On the other hand, the CBS-Q ab initio computational approach generated BDEs which are in good agreement with experimental values without a correction factor.     

气相中O3与HSO自由基间的氢键复合物

气相中O3与HSO自由基之间的相互作用及其反应在大气化学中非常重要. 在DFT-B3LYP/6-311++G**和MP2/6-311++G**水平上求得O3+HSO复合物势能面上的稳定构型, B3LYP方法得到了三种构型(复合物I, II和III), 而MP2方法只能得到一种构型(复合物II). 在复合物I和III中, HSO单元中的1H原子作为质子供体, 与O3分子中的端基O原子作为质子受体相互作用, 形成红移氢键复合物; 而在复合物II中, 虽与复合物I和III中具有相同的质子供体和质子受体, 却形成了蓝移氢键复合物. B3LYP/6-311++G**水平上计算的单体间相互作用能的计算考虑了基组重叠误差(BSSE)和零点振动能(ZPVE)校正, 其值在-3.37到-4.55 kJ·mol-1之间. 采用自然键轨道理论(NBO)对单体间相互作用的本质进行了考查, 并通过分子中原子理论(AIM)分析了三种复合物中氢键的电子密度拓扑性质.