def2-TZVPPD Basis Sets Research Articles

Analysis of two overlapping fragmentation approaches in density matrix construction: GMBE-DM vs. ADMA.

In this study, we conduct a comparative analysis of two density matrix construction methods: the generalized many-body expansion for building density matrices (GMBE-DM) based on the set-theoretical principle of inclusion/exclusion and the adjustable density matrix assembler (ADMA) based on the Mulliken-Mezey ansatz. We apply these methods to various noncovalent clusters, including water clusters, ion-water clusters, and ion-pair clusters, using both small 6-31G(d) and large def2-TZVPPD basis sets. Our findings reveal that the GMBE-DM method, particularly when combined with the purification scheme and truncation at the one-body level [GMBE(1)-DM-P], exhibits superior performance across all test systems and basis sets. In contrast, all ADMA set of methods show reasonable results only with small and compact basis sets. For example, GMBE(1)-DM-P outperforms the best ADMA method by at least 4 and 16 times with small and large basis sets, respectively, in the case of (H2O)N=6-55. This highlights the significance of the basis set choice for ADMA, which is even more critical than the fragmentation scheme, such as the size of subsystems, while GMBE-DM consistently produces accurate results irrespective of the chosen basis set. Consequently, the efficient and robust GMBE(1)-DM-P approach is recommended as a fragmentation method for generating accurate absolute and relative energies across different binding patterns and basis sets for noncovalent clusters.

Understanding the molecular mechanism of the chlorine atom transfer between ammonia and hypochlorous acid with electron localisation function (ELF)

The nature of the chlorine atom transferred between ammonia and hypochlorous acid was investigated using the bonding evolution theory (BET). The application of topological analysis of electron localisation function (ELF) and Thom's catastrophe theory enabled a detailed description of the evolution of chemical bonds and non-bonding electron density. The BET analysis was performed at the DFT(MN15)/def2-TZVPPD computational level. The Hausdorff distance was used as a qualitative view on the evolution of the ELF field. The molecular mechanism of the chlorination reaction consisted of 11 structural stability domains (SSDs) and 10 ELF catastrophes. Four SSDs were observed before the transition state (TS). The TS was found in SSD V, and six SSDs were distinguished after the TS. The energy barriers were analysed with six different DFT functionals (B3LYP, CAM-B3LYP, LC-ωPBE, M11, MN15, ωB97XD) and def2-TZVPPD basis set. The energy barrier between transition state and reagents was smaller than transition state and products, thus the chlorination of ammonia is thermodynamically favourable. Alongside BET analysis, NPA (natural population analysis), AIM (atoms in molecules) and APT (atomic polar tensor) charge analysis was also applied to the transferred chlorine atom.

Thermochemistry of the Ir+ + SO2 reaction using guided ion beam tandem mass spectrometry and theory.

Kinetic energy dependences of the reactions of Ir+ (5F5) with SO2 were studied using a guided ion beam tandem mass spectrometer and theory. The observed cationic products are IrO+, IrS+, and IrSO+, formed in endothermic reactions. Bond dissociation energies (BDEs) of the products are determined by modeling the kinetic energy dependent product cross sections: D0(Ir+-O) = 4.27 ± 0.11 eV, D0(Ir+-S) = 4.03 ± 0.06 eV, and D0(Ir+-SO) ≥ 2.95 ± 0.06 eV. The oxide BDE agrees well with literature values, whereas the two latter results are novel measurements. Quantum mechanical calculations are performed at the B3LYP level of theory using the def2-TZVPPD basis set for all product BDEs with additional calculations for IrS+, IrO2 +, and IrSO+ at the coupled cluster with single, double, and perturbative triple excitation levels with def2-QZVPPD and aug-cc-pVXZ (X = T and Q and for IrS+, also X = 5) basis sets and complete basis set extrapolations. These theoretical BDEs agree reasonably well with the experimental values. 1A1 (IrO2 +), 5Δ4 (IrS+), and 3A″/1A' (IrSO+) are found to be the ground states after including empirical spin-orbit corrections. The potential energy surfaces including intermediates and transition states for each reaction are also calculated at the B3LYP/def2-TZVPPD level. The formation of MO+ (M = Re, Os, and Ir) from M+ + SO2 reactions is compared with those from the M+ + O2 and M+ + CO reactions, where interesting trends in cross sections are observed. Overall, these studies suggest that the M+ + O2 reactions had restrictions associated with reactions along A' and A″ surfaces.

Thermochemical studies of reactions of Re+ with SO2 using guided ion beam experiments and theory.

The kinetic energy dependent reactions of Re+ with SO2 were studied with guided ion beam tandem mass spectrometry. ReO+, ReO2+, and OReS+ species were observed as products, all in endothermic reactions. Modeling of the kinetic energy dependent cross sections yields 0 K bond dissociation energies (BDEs, in eV) of 4.78 ± 0.06 (Re+-O), 5.75 ± 0.02 (Re+-O2), and 4.35 ± 0.14 (Re+-SO). The latter two values can be combined with other information to derive the additional values 6.05 ± 0.05 (ORe+-O) and 4.89 ± 0.19 (ORe+-S). BDEs of ReO+ and ReO2+ agree with literature values whereas the values for OReS+ are the first measurements. The former result is obtained even though formation of ground state ReO+ is spin-forbidden. Quantum mechanical calculations at the B3LYP level of theory with a def2-TZVPPD basis set yield results that agree reasonably well with experimental values. Additional calculations at the BP86 and CCSD(T) levels of theory using def2-QZVPPD and aug-cc-pVxZ (x = T, Q, and 5) basis sets were performed to compare thermochemistry with experiment to determine that ReO2+ rather than the isobaric ReS+ is formed. Product ground states are 3Δ3(ReO+), 3B1(OReO+), 5Π-1(ReS+), and 3A''(OReS+) after including empirical spin-orbit corrections, which means that formation of ground state products is spin-forbidden for all three product channels. The potential energy surfaces for the ReSO2+ system were also explored at the B3LYP/def2-TZVPPD level and exhibited no barriers in excess of the endothermicities for all products. BDEs for rhenium oxide and sulfide diatomics and triatomics are compared and discussed. The present result for formation of ReO+ is compared to that for formation of ReO+ in the reactions of Re+ + O2 and CO, where the former system exhibited interesting dual cross section features. Results are consistent with the hypothesis that the distinction of in-plane and out-of-plane CS symmetry in the triatomic systems might be the explanation for the two endothermic features observed in the Re+ + O2 reaction.

Hypervalent halogen hydrides HalHn (Hal = Cl, Br, I; n = 3, 5, 7): DFT and ab initio stability prediction

Stability—both kinetic and thermodynamic—and structure of hypervalent halogen hydrides and respective ionic species were investigated using B3LYP, B3PW91 and MPW1PW91 DFT functionals and second-order Moller-Plesset perturbation theory (MP2) with Def2-TZVPPD basis set (for I atoms) and aug-cc-pVTZ basis set (for Br, Cl and H atoms). Various dissociation pathways were considered: respective reaction enthalpies and H2 elimination activation enthalpies were calculated. At these levels of theory, all of the HalHn species (Hal = Cl, Br, I; n = 3, 5, 7) have covalently bound local minima (except for ClH7) and have a relatively significant kinetic stability with respect to H2 elimination, which indicates the possibility of them being, at the very least, spectroscopically observable, if not isolable. The first candidates for isolation/observation appear to be IH3 and IH5. To aid future identification of hypervalent halogen hydrides, harmonic frequencies for these compounds were calculated. Also, NBO analyses were performed and natural electron configurations, Wiberg bond indices and HOMO–LUMO gaps for these species were determined.

Bond dissociation energy of Au2 +: A guided ion beam and theoretical investigation.

Guided ion beam tandem mass spectrometry was employed to measure the kinetic energy-dependent product ion cross sections for the collision induced dissociation of Au2 + with Xe. Gold dimer cations were formed in a glow discharge flow tube source that should create ions in their ground electronic state with thermal internal energies. Analysis of the endothermic kinetic energy dependent cross section accounts for multiple collisions, lifetime effects, and the internal energy of the reactant ion. The value obtained for the bond dissociation energy (BDE) of Au2 + is the first direct measurement and is reported here as 2.20 ± 0.21 eV. For comparison with experimental results, theoretical calculations were also completed at the B3LYP, M06-2X, and coupled cluster singles, doubles, and perturbative triples [CCSD(T,full)] levels of theory using the def2-TZVPPD basis set and at the CCSD(T)-F12/correlation-consistent polarized valence triple zeta basis with pseudopotential level. These results predict a 2Σg + electronic ground state for Au2 + with BDEs calculated at the B3LYP and both CCSD(T) levels of theory in agreement with the experiment within the uncertainty. Several electronically excited states are also evaluated theoretically.

Thermodynamic properties of 2–methylindole: Experimental and computational results for gas-phase entropy and enthalpy of formation

Abstract Measurements leading to the calculation of thermodynamic properties in the ideal-gas state for 2-methylindole (Chemical Abstracts registry number [95–20–5]) are reported. Experimental methods were adiabatic heat-capacity calorimetry, differential scanning calorimetry (d.s.c.), comparative ebulliometry, inclined-piston manometry, and oxygen bomb calorimetry. The critical temperature of 2-methylindole was determined experimentally with d.s.c. Molar thermodynamic functions for the condensed and ideal-gas states were derived from the experimental results. Statistical calculations were performed based on molecular geometry optimization and vibrational frequencies using B3LYP hybrid density functional theory with the def2-TZVPPD basis set. Excellent accord between computed and experimentally-derived ideal-gas entropies is shown. The enthalpy of formation for 2-methylindole in the gas phase was computed with an atomization-based protocol described recently, and excellent agreement with the experimental values is seen. The experimental literature for enthalpies of formation in the gas phase for 1- and 2-ring pyrrollic compounds is reviewed, and comparisons with computed values further support the findings here. All experimental results are compared with property values reported in the literature, where possible.

Noble gas supported boron-pentagonal clusters B5Ngn3+: exploring the structures and bonding.

A novel type of trivalent BNg five-membered cational species B5Ngn3+(Ng = He~Rn, n = 1~5) has been found and investigated theoretically using the B3LYP and MP2 methods with the def2-QZVPPD and def2-TZVPPD basis sets. The geometry, harmonic vibrational frequencies, bond energies, charge distribution, bond nature, aromaticity, and energy decomposition analysis of these structures were reported. The calculated B-Ng bond energy is quite large (the averaged bond energy is in the range of 209.2~585.76 kJ mol-1) for heavy rare gases and increases with the Ng atomic number. The analyses of the molecular wavefunction show that in the BNg compounds of heavy Ng atoms Ar~Rn, the B-Ng bonds are of typical covalent character. Nuclear independent chemical shifts display that both B53+ and B5Ngn3+(n=1~5) have obvious aromaticity. Energy decomposition analysis shows that these BNg compounds are mainly stabilized by the σ-donation from the Ng valence p orbital to the B53+ LUMO. These findings offer valuable clues toward the design and synthesis of new stable Ng-containing compounds.

Thermodynamic properties of pyrrole, 1-methylpyrrole, 2,4-dimethylpyrrole, and 2,5-dimethylpyrrole: Experimental and computational results

New measurements of critical temperature, saturated liquid heat capacity from temperature T≈315K to T≈550K, and saturated liquid density from T=323K to T≈475K are reported for pyrrole (Chemical Abstracts registry number [109-97-7]), 1-methylpyrrole [96-54-8], 2,4-dimethylpyrrole [625-82-1], and 2,5-dimethylpyrrole [625-84-3]. New measurements of vapor pressure are reported for 2,4-dimethylpyrrole {338<(T/K)<477} and 2,5-dimethylpyrrole {341<(T/K)<479}, as well as the enthalpy of combustion determined with oxygen-bomb calorimetry for 2,4-dimethylpyrrole. These new measurements are combined with literature values to calculate thermodynamic properties in the ideal-gas state for extended ranges of temperature for all compounds: pyrrole {298<(T/K)<550}, 1-methylpyrrole {298<(T/K)<530}, 2,4-dimethylpyrrole {298<(T/K)<600}, and 2,5-dimethylpyrrole {298<(T/K)<560}. Molar thermodynamic functions (enthalpies, entropies, and Gibbs energies) for the liquid and ideal-gas states were derived from the experimental studies at selected temperatures. Statistical calculations were performed with optimized geometries, scaled vibrational frequencies, and methyl rotational potentials performed using B3LYP hybrid density functional theory with the def2-TZVPPD basis set. Methyl torsional barriers were evaluated at the DLPNO-CCSD(T)/def2-QZVP level of theory using the B3LYP/def2-TZVPPD geometries. Computed ideal-gas properties are shown to be in excellent accord with ideal-gas entropies derived from thermophysical property measurements, as well as with reported experimental heat capacities for the ideal-gas state. All experimental results are compared with property values available in the literature.

DSD-PBEP86-NL and DOD-PBEP86-NL functionals for noncovalent interactions: Basis set effects and tentative applications to large noncovalent systems

Basis set effects on the DSD-PBEP86-NL and DOD-PBEP86-NL functionals for noncovalent interactions have been extensively studied in this work. The cc-pVXZ (X = D, T, Q, 5, 6) and augmented aug-cc-pVXZ (X = D, T, Q) basis sets are systematically tested without counterpoise (CP) corrections against the well-known S66 database. Additionally, the basis sets of def2-TZVPP and def2-TZVPPD are also examined. Based on our computations, the performances of the aug-cc-pVQZ, cc-pV5Z, and cc-pV6Z basis sets are very approximate to those obtained with the def2-QZVP basis set for both the DSD-PBEP86-NL and DOD-PBEP86-NL functionals. Note that the short-range attenuation parameters for these two functionals were directly optimized using the def2-QZVP basis set without CP corrections against the S66 database. Generally speaking, the cc-pVXZ (X = D, T, Q), aug-cc-pVXZ (X = D, T, Q), def2-TZVPP, and def2-TZVPPD basis sets favor half CP correction for these two functionals. Nevertheless, the aug-cc-pVQZ basis set already performs well without any CP correction, especially for the DOD-PBEP86-NL functional. With respect to accuracy and computational cost, the cc-pVTZ and def2-TZVPP basis sets with half CP corrections are recommended for these two functionals to evaluate interaction energies of large noncovalent complexes.

Matrix-isolation and computational study of H2CCCl and H2CCBr radicals.

We report on two new radicals, H2CCCl and H2CCBr, prepared in low-temperature noble-gas matrices and characterized using infrared spectroscopy. These radicals are made by UV photolysis of HCCCl and HCCBr and subsequent thermal annealing to mobilize hydrogen atoms in the matrices and promote their reaction with the residual precursor molecules. Three characteristic infrared bands are observed for each radical. The assignments are supported by quantum chemical calculations at the B3LYP and CCSD(T) levels of theory with the def2-TZVPPD basis set.

The electronic structure of the xenon insertion compounds XXe–MX2 (X = F, Cl, Br, I; M = B, Al, Ga)

Nature of chemical bonding in the XXe–MX2 (X=F, Cl, Br, I; M=B, Al, Ga) molecules have been investigated using topological analysis of Electron Localisation Function (ELF) and electron density. The wave function has been approximated by the DFT(B3LYP, M062X) and the CCSD methods using all-electron basis sets (TZVPPall, TZVall+) and the ecp-28 approximation for Xe, I with Def2-TZVPPD basis set. Electronic structure of the XXeMX2 consists of the [Xe–MX2] fragment, bound non-covalently to the halogen atom, Xδ−⋯[Xe–MX2]δ+ (δ=0.4–0.7e). The xenon–icosagen bond, Xe–M displays covalent nature in all studied molecules, confirmed by bonding disynaptic attractor V(Xe,M) for the ELF field with basin population ranging from 1.1 to 1.8e. The basin population of the Xe–M bond decreases together with the r(Xe–M) bond length increase and the increase of the polarisability of the respective halogen. The Xe–X bond always shows a non-covalent nature, confirmed by an absence of bonding attractor in the region between the core basins C(Xe) and C(X). M–X bonds can have either covalent or ionic nature depending on the type of icosagen and halogen involved in bonding. The Al–F and Ga–X (X=F, Cl, Br) bonds have ionic nature, except for the Ga–I bond, showing slightly covalent but highly delocalised character, confirmed by the presence of the bonding basin V(Ge,X) with the population of 0.64e. In addition to bonding population calculations, the analysis of the bond polarity, based on topological analysis of electron density combined with ELF fields, and delocalisation of electron density in bonding basins has been performed and discussed.

Matrix-isolation and ab initio study of HKrCCCl and HXeCCCl.

We report on two new noble-gas molecules, HKrCCCl and HXeCCCl, prepared in low-temperature Kr and Xe matrices. These molecules are made by UV photolysis of HCCCl in the matrices and subsequent thermal annealing. The HCCCl precursor is produced by microwave discharge of a mixture of a matrix gas with trichloroethylene (HClC=CCl2). The assignments of the new noble-gas molecules are supported by deuteration experiments and quantum chemical calculations at the MP2(full) and CCSD(T) levels of theory with the def2-TZVPPD basis set. No evidence of ClXeCCH, which is computationally reliably stable, is found in the experiments. ClKrCCH as well as the Ar compounds HArCCCl and ClArCCH are not observed either, which is in agreement with the calculations.

Alkali metal cation binding affinities of cytosine in the gas phase: revisited.

Binding of metal cations to the nucleobases can influence base pairing, base stacking and nucleobase tautomerism. Gas-phase condensation of dc discharge generated alkali metal cations and thermally vaporized cytosine (DC/FT) has been found to produce kinetically trapped excited tautomeric conformations of the M(+)(cytosine) complexes, which influences the threshold collision-induced dissociation (TCID) behavior. In order to elucidate the effects of the size of alkali metal cation on the strength of binding to the canonical form of cytosine, the binding affinities of Na(+) and K(+) to cytosine are re-examined here, and studies are extended to include Rb(+) and Cs(+) again using TCID techniques. The M(+)(cytosine) complexes are generated in an electrospray ionization source, which has been shown to produce ground-state tautomeric conformations of M(+)(cytosine). The energy-dependent cross sections are interpreted to yield bond dissociation energies (BDEs) using an analysis that includes consideration of unimolecular decay rates, the kinetic and internal energy distributions of the reactants, and multiple M(+)(cytosine)-Xe collisions. Revised BDEs for the Na(+)(cytosine) and K(+)(cytosine) complexes exceed those previously measured by 31.9 and 25.5 kJ mol(-1), respectively, consistent with the hypothesis proposed by Yang and Rodgers that excited tautomeric conformations are accessed when the complexes are generated by DC/FT ionization. Experimentally measured BDEs are compared to theoretical values calculated at the B3LYP and MP2(full) levels of theory using the 6-311+G(2d,2p)_HW* and def2-TZVPPD basis sets. The B3LYP/def2-TZVPPD level of theory is found to provide the best agreement with the measured BDEs, suggesting that this level of theory can be employed to provide reliable energetics for similar metal-ligand systems.

Metal Cation Dependence of Interactions with Amino Acids: Bond Dissociation Energies of Rb+ and Cs+ to the Acidic Amino Acids and Their Amide Derivatives

Metal cation-amino acid interactions are key components controlling the secondary structure and biological function of proteins, enzymes, and macromolecular complexes comprising these species. Determination of pairwise interactions of alkali metal cations with amino acids provides a thermodynamic vocabulary that begins to quantify these fundamental processes. In the present work, we expand a systematic study of such interactions by examining rubidium and cesium cations binding with the acidic amino acids (AA), aspartic acid (Asp) and glutamic acid (Glu), and their amide derivatives, asparagine (Asn) and glutamine (Gln). These eight complexes are formed using electrospray ionization and their bond dissociation energies (BDEs) are determined experimentally using threshold collision-induced dissociation with xenon in a guided ion beam tandem mass spectrometer. Analyses of the energy-dependent cross sections include consideration of unimolecular decay rates, internal energy of the reactant ions, and multiple ion-neutral collisions. Quantum chemical calculations are conducted at the B3LYP, MP2(full), and M06 levels of theory using def2-TZVPPD basis sets, with results showing reasonable agreement with experiment. At 0 and 298 K, most levels of theory predict that the ground-state conformers for M(+)(Asp) and M(+)(Asn) involve tridentate binding of the metal cation to the backbone carbonyl, amino, and side-chain carbonyl groups, although tridentate binding to the carboxylic acid group and side-chain carbonyl is competitive for M(+)(Asn). For the two longer side-chain amino acids, Glu and Gln, multiple structures are competitive. A comparison of these results to those for the smaller alkali cations, Na(+) and K(+), provides insight into the trends in binding energies associated with the molecular polarizability and dipole moment of the side chain. For all four metal cations, the BDEs are inversely correlated with the size of the metal cation and follow the order Asp < Glu < Asn < Gln.

Metal Cation Dependence of Interactions with Amino Acids: Bond Energies of Rb+ and Cs+ to Met, Phe, Tyr, and Trp

The interactions of rubidium and cesium cations with four amino acids (AA) including methionine (Met), phenylalanine (Phe), tyrosine (Tyr), and tryptophan (Trp) are examined in detail. Experimentally, the bond dissociation energies (BDEs) are determined using threshold collision-induced dissociation of the Rb(+)(AA) and Cs(+)(AA) complexes with xenon in a guided ion beam tandem mass spectrometer. Analyses of the energy dependent cross sections include consideration of unimolecular decay rates, internal energy of the reactant ions, and multiple ion-neutral collisions. 0 K BDEs of 121.0 ± 7.0 (102.8 ± 6.6), 123.8 ± 7.2 (112.9 ± 5.5), 125.8 ± 7.4 (115.6 ± 6.9), and 138.1 ± 7.5 (125.0 ± 6.8) kJ/mol are determined for complexes of Rb(+) (Cs(+)) with Met, Phe, Tyr, and Trp, respectively. Quantum chemical calculations are conducted at the B3LYP, MP2(full), and M06 levels of theory with geometries and zero point energies calculated at the B3LYP level using def2-TZVPPD basis sets. Results obtained using all three levels show good agreement with experiment, with B3LYP values being systematically low and MP2(full) and M06 values being systematically high. At 0 and 298 K, theory predicts the ground-state conformers for M(+)(Met) either have tridentate binding of the metal cation to the carbonyl, amino, and sulfur groups (MP2 and M06) or to both oxygens of a zwitterionic conformation (B3LYP). At 298 K, binding to the carboxylic acid group and the sulfur also becomes competitive. For the aromatic amino acids at 0 K, most levels of theory favor tridentate binding of the metal ions to the backbone carbonyl and amino groups along with the π-cloud of the ring, whereas for Rb(+)(Trp) and Cs(+)(AA), B3LYP theory favors binding to only the carbonyl and ring groups. At 298 K, B3LYP favors the latter binding mode for all three Rb(+)(aromatic AA) complexes. Comparison of these results to those for the smaller alkali cations provides insight into the trends in binding affinities and structures associated with metal cation variations.

Thermochemistry of Alkali Metal Cation Interactions with Histidine: Influence of the Side Chain

The interactions of alkali metal cations (M(+) = Na(+), K(+), Rb(+), Cs(+)) with the amino acid histidine (His) are examined in detail. Experimentally, bond energies are determined using threshold collision-induced dissociation of the M(+)(His) complexes with xenon in a guided ion beam tandem mass spectrometer. Analyses of the energy dependent cross sections provide 0 K bond energies of 2.31 ± 0.11, 1.70 ± 0.08, 1.42 ± 0.06, and 1.22 ± 0.06 eV for complexes of His with Na(+), K(+), Rb(+), and Cs(+), respectively. All bond dissociation energy (BDE) determinations include consideration of unimolecular decay rates, internal energy of reactant ions, and multiple ion-neutral collisions. These experimental results are compared to values obtained from quantum chemical calculations conducted previously at the MP2(full)/6-311+G(2d,2p), B3LYP/6-311+G(2d,2p), and B3P86/6-311+G(2d,2p) levels with geometries and zero point energies calculated at the B3LYP/6-311+G(d,p) level where Rb and Cs use the Hay-Wadt effective core potential and basis set augmented with additional polarization functions (HW*). Additional calculations using the def2-TZVPPD basis set with B3LYP geometries were conducted here at all three levels of theory. Either basis set yields similar results for Na(+)(His) and K(+)(His), which are in reasonable agreement with the experimental BDEs. For Rb(+)(His) and Cs(+)(His), the HW* basis set and ECP underestimate the experimental BDEs, whereas the def2-TZVPPD basis set yields results in good agreement. The effect of the imidazole side chain on the BDEs is examined by comparing the present results with previous thermochemistry for other amino acids. Both polarizability and the local dipole moment of the side chain are influential in the energetics.

Alkali metal cation interactions with 12-crown-4 in the gas phase: Revisited

Quantitative interactions of alkali metal cations with the cyclic 12-crown-4 polyether ligand (12C4) are studied. Experimentally, Rb+(12C4) and Cs+(12C4) complexes are formed using electrospray ionization and their bond dissociation energies (BDEs) determined using threshold collision-induced dissociation of these complexes with xenon in a guided ion beam tandem mass spectrometer. The energy-dependent cross sections thus obtained are interpreted using an analysis that includes consideration of unimolecular decay rates, internal energy of the reactant ions, and multiple ion-neutral collisions. 0K BDEs of 151.5±9.7 and 137.0±8.7kJ/mol, respectively, are determined and exceed those previously measured by 60 and 54kJ/mol, respectively, consistent with the hypothesis proposed there that excited conformers had been studied. In order to provide comparable thermochemical results for the Na+(12C4) and K+(12C4) systems, the published data for these systems are reinterpreted using the same analysis techniques, which have advanced since the original data were acquired. Revised BDEs for these systems are obtained as 243.9±12.6 and 182.0±17.3kJ/mol, respectively, which are within experimental uncertainty of the previously reported values. In addition, quantum chemical calculations are conducted at the B3LYP and MP2(full) levels of theory with geometries and zero point energies calculated at the B3LYP level using both HW*/6-311+G(2d,2p) and def2-TZVPPD basis sets. The theoretical results are in reasonable agreement with experiment, with B3LYP/def2-TZVPPD values being in particularly good agreement. Computations also allow the potential energy surfaces for dissociation of the M+(12C4) complexes to be elucidated. These are used to help explain why the previous studies formed excited conformers of Rb+(12C4) and Cs+(12C4) but apparently not of Na+(12C4) and K+(12C4).

Metal Cation Dependence of Interactions with Amino Acids: Bond Energies of Cs+ to Gly, Pro, Ser, Thr, and Cys

The interactions of cesium cations with five amino acids (AA) including glycine (Gly), proline (Pro), serine (Ser), threonine (Thr), and cysteine (Cys) are examined in detail. Experimentally, the bond dissociation energies (BDEs) are determined using threshold collision-induced dissociation of the Cs(+)(AA) complexes with xenon in a guided ion beam tandem mass spectrometer. Analyses of the energy-dependent cross sections include consideration of unimolecular decay rates, internal energy of the reactant ions, and multiple ion-neutral collisions. Bond dissociation energies (0 K) of 93.3 ± 2.5, 107.9 ± 4.6, 102.3 ± 4.1, 105.4 ± 4.3, and 96.8 ± 4.2 kJ/mol are determined for complexes of Cs(+) with Gly, Pro, Ser, Thr, and Cys, respectively. Quantum chemical calculations are conducted at the B3LYP, B3P86, MP2(full), and M06 levels of theory with geometries and zero-point energies calculated at the B3LYP level using both HW*/6-311+G(2d,2p) and def2-TZVPPD basis sets. Results obtained using the former basis sets are systematically low compared to the experimental bond energies, whereas the latter basis sets show good agreement. For Cs(+)(Gly), theory predicts the ground-state conformer has the cesium cation binding to the carbonyl group of the carboxylic acid. For Cs(+)(Pro), the secondary nitrogen accepts the carboxylic acid hydrogen to form the zwitterionic structure, and the metal cation binds to both oxygens. Cs(+)(Ser), Cs(+)(Thr), and Cs(+)(Cys) are found to have tridentate binding at the MP2(full) level, whereas the density functional approaches slightly prefer bidentate binding of Cs(+) at the carboxylic acid moiety. Comparison of these results to those for the smaller alkali cations provides insight into the trends in binding affinities and structures associated with metal cation variations.