Hydration energies and structures of alkaline earth metal ions, M2+(H2O)n, n = 5-7, M = Mg, Ca, Sr, and Ba

The evaporation of water from hydrated alkaline earth metal ions, produced by electrospray ionization, was studied in a Fourier transform mass spectrometer. Zero-pressure-limit dissociation rate constants for loss of a single water molecule from the hydrated divalent metal ions, M(2+)(H(2)O)(n) (M = Mg, Ca, and Sr for n = 5-7, and M = Ba for n = 4-7), are measured as a function of temperature using blackbody infrared radiative dissociation. From these values, zero-pressure-limit Arrhenius parameters are obtained. By modeling the dissociation kinetics using a master equation formalism, threshold dissociation energies (E(o)) are determined. These reactions should have a negligible reverse activation barrier; therefore, E(o) values should be approximately equal to the binding energy or hydration enthalpy at 0 K. For the hepta- and hexahydrated ions at low temperature, binding energies follow the trend expected on the basis of ionic radii: Mg > Ca > Sr > Ba. For the hexahydrated ions at high temperature, binding energies follow the order Ca > Mg > Sr > Ba. The same order is observed for the pentahydrated ions. Collisional dissociation experiments on the tetrahydrated species result in relative dissociation rates that directly correlate with the size of the metals. These results indicate the presence of two isomers for hexahydrated magnesium ions: a low-temperature isomer in which the six water molecules are located in the first solvation shell, and a high-temperature isomer with the most likely structure corresponding to four water molecules in the inner shell and two water molecules in the second shell. These results also indicate that the pentahydrated magnesium ions have a structure with four water molecules in the first solvation shell and one in the outer shell. The dissociation kinetics for the hexa- and pentahydrated clusters of Ca(2+), Sr(2+), and Ba(2+) are consistent with structures in which all the water molecules are located in the first solvation shell.     

Structures and energetics of the cation-pi interactions of Li+, Na+, and K+ with cup-shaped molecules: effect of ring addition to benzene and cavity selectivity

The interactions of alkali metal cations (Li (+), Na (+), and K (+)) with the cup-shaped molecules, tris(bicyclo[2.2.1]hepteno)benzene and tris(7-azabicyclo[2.2.1]hepteno)benzene have been investigated using MP2(FULL)/6-311+G(d,p)//MP2/6-31G(d) level of theory. The geometries and interaction energies obtained for the metal ion complexation with the cup-shaped systems trindene and benzotripyrrole are compared with the results for benzene-metal ion complexes to examine the effect of ring addition to the benzene on structural and binding affinities. The cup-shaped molecules exhibit two faces or cavities (top and bottom). Except for one of the conformers of tris(7-azabicyclo[2.2.1]hepteno)benzene), the metal ions prefer to bind with the top face over bottom face of the cup-shaped molecules. The selectivity of the top face is due to strong interaction of the cation with the pi cloud not only from the central six-membered ring but also from the pi electrons of rim C=C bonds. In contrast, the metal ions under study exhibit preference to bind with the bottom face rather than top face of tris(7-azabicyclo[2.2.1]hepteno)benzene) when the lone pair of electrons of three nitrogen atoms participates in binding with metal ions. This bottom face selectivity could be ascribed to the combined effect of the cation-pi and strong cation-lone pair interactions. As evidenced from the values of pyramidalization angles, the host molecule becomes deeper bowl when the lone pair of electrons of nitrogen atoms participates in binding with cation. Molecular electrostatic potential surfaces nicely explain the cavity selectivity in the cup-shaped systems and the variation of interaction energies for different ligands. Vibrational frequency analysis is useful in characterizing different metal ion complexes and to distinguish top and bottom face complexes of metal ions with the cup-shaped molecules.     

Dual-scale modeling of benzene adsorption onto Ni(111) and Au(111) surfaces in explicit water.

We present a multiscale modeling approach for studying interactions of organic molecules with metal surfaces in explicit water. The approach is based on combining adsorption energies of isolated molecules on transition metal surfaces calculated by ab initio density functional methods and classical molecular dynamics simulations using atomistically detailed force fields. The interaction of benzene with Ni(111) and Au(111) surfaces was studied. It is shown that a strong affinity of water for the hydrophilic surfaces makes benzene adsorption on Au thermodynamically unfavorable, while on Ni there is no preference. The work presented here serves as a first step in modeling the interactions of larger organic molecules with metal surfaces.     

Hydrogen bonds in gas-phase clusters between halide ions and olefins

Gas-phase clustering reactions of halide ions (X- = F-, Cl-, Br-, and I-) with ethylene (C2H4) and propylene (C3H6) were studied with a pulsed electron beam mass spectrometer. Bonding energies of all cluster ions were found to be less than 10 kcal/mol, i.e., no anion-initiated polymerization of C2H4 and C3H6 took place. For the cluster F-(C2H4)n, a small gap in the binding energy is observed between n = 4 and 5 suggesting that the first shell is completed with n = 4. For larger halide ions, the bond energies for the clusters X-(C2H4)n were found to be nearly n independent. For Cl-(C3H6)n a steep decrease in binding energies was observed between n = 2 and 3 and n = 3 and 4. The structure of the cluster ions was investigated by ab initio calculations. X-(C2H4)n complexes were calculated to have hydrogen-bond geometries regardless of the identity of the halide ions, and bidentate (chelate) type geometries of X-(C3H6)1 were found.     

Probing the structural and magnetic properties of transition metal-benzene anion complexes

Two types of transition metal-benzene anion complexes, (titanium)(n)(benzene)(m)? and (cobalt)(n)(benzene)(m)? (n ≤ 2, m ≤ 3) have been determined using density functional theory. The photoelectron spectra of Ti(n)Bz(m)? and Co(n)Bz(m)? (n ≤ 2, m ≤ 3) were discussed from the perspective of quantum chemical calculations of the vertical detachment energies (VDEs) of several low-energy isomers obtained by the structural optimization procedure. The binding of Ti and Co atoms to benzene molecules is accounted by 3d-π bonds, as revealed by the molecular orbitals. The topology of the electronic density has been analyzed, suggesting that the C-C bonds were weakened in the transition metal-benzene complexes in comparison to those in free benzene. Spin density distribution results show the spin densities for Ti(n)Bz(m)? and Co(n)Bz(m)? (n ≤ 2, m ≤ 3) reside mainly on the metal Ti and Co centers (70%-90%). A shift to lower magnetic moment with respect to the pure titanium/cobalt cluster anions indicates the solvent benzene molecule acts to demagnetize the bare titanium/cobalt cluster anions.     

Gas phase hydration of organic ions

In this work, we study the hydration phenomenon on a molecular level in the gas phase where a selected number of water molecules can interact with the organic ion of interest. The stepwise binding energies (DeltaH degrees (n-1,n)) of 1-7 water molecules to the phenyl acetylene cation are determined by equilibrium measurements using an ion mobility drift cell. The stepwise hydration energies DeltaH degrees (n-1,n) are nearly constant at 39.7 +/- 6.3 kJ mol(-1) from n = 1 to 7. The entropy change is larger in the n = 7 step, suggesting cyclic or cage-like water structures. No water addition is observed on the ionized phenyl acetylene trimer consistent with cyclization of the trimer ion to form triphenyl benzene cations C(24)H(18) (+) which are expected to interact weakly with the water molecules due to steric interactions and the delocalization of the charge on the large organic ion. The work demonstrates that hydration studies of organic ions can provide structural information on the organic ions.     

Reactivity and infrared spectroscopy of gaseous hydrated trivalent metal ions

Hydrated trivalent rare earth metal ions containing yttrium and all naturally abundant lanthanide metals are formed using electrospray ionization, and the structures and reactivities of these ions containing 17-21 water molecules are probed using blackbody infrared radiative dissociation (BIRD) and infrared action spectroscopy. With the low-energy activation conditions of BIRD, there is an abrupt transition in the dissociation pathway from the exclusive loss of a single neutral water molecule to the exclusive loss of a small protonated water cluster via a charge-separation process. This transition occurs over a narrow range of cluster sizes that differs by only a few water molecules for each metal ion. The effective turnover size at which these two dissociation rates become equal depends on metal ion identity and is poorly correlated with the third ionization energies of the isolated metals but is well correlated with the hydrolysis constants of the trivalent metal ions in bulk aqueous solution. Infrared action spectra of these ions at cluster sizes near the turnover size are largely independent of the specific identity of the trivalent metal ion, suggesting that any differences in the structures of the ions present in our experiment are subtle.     

New relativistic ANO basis sets for transition metal atoms

New basis sets of the atomic natural orbital (ANO) type have been developed for the first, second, and third row transition metal atoms. The ANOs have been obtained from the average density matrix of the ground and lowest excited states of the atom, the positive and negative ions, and the atom in an electric field. Scalar relativistic effects are included through the use of a Douglas-Kroll-Hess Hamiltonian. Multiconfigurational wave functions have been used with dynamic correlation included using second order perturbation theory (CASSCF/CASPT2). The basis sets are applied in calculations of ionization energies, electron affinities, and excitation energies for all atoms and polarizabilities for spherically symmetric atoms. These calculations include spin-orbit coupling using a variation-perturbation approach. Computed ionization energies have an accuracy better than 0.2 eV in most cases. The accuracy of computed electron affinities is the same except in cases where the experimental values are smaller than 0.5 eV. Accurate results are obtained for the polarizabilities of atoms with spherical symmetry. Multiplet levels are presented for some of the third row transition metals.     

Influence of d orbital occupation on the binding of metal ions to adenine

Threshold collision-induced dissociation of M(+)(adenine) with xenon is studied using guided ion beam mass spectrometry. M(+) includes all 10 first-row transition metal ions: Sc(+), Ti(+), V(+), Cr(+), Mn(+), Fe(+), Co(+), Ni(+), Cu(+), and Zn(+). For the systems involving the late metal ions, Cr(+) through Cu(+), the primary product corresponds to endothermic loss of the intact adenine molecule, whereas for Zn(+), this process occurs but to form Zn + adenine(+). For the complexes to the early metal ions, Sc(+), Ti(+), and V(+), intact ligand loss competes with endothermic elimination of purine and of HCN to form MNH(+) and M(+)(C(4)H(4)N(4)), respectively, as the primary ionic products. For Sc(+), loss of ammonia is also a prominent process at low energies. Several minor channels corresponding to formation of M(+)(C(x)H(x)N(x)), x = 1-3, are also observed for these three systems at elevated energies. The energy-dependent collision-induced dissociation cross sections for M(+)(adenine), where M(+) = V(+) through Zn(+), are modeled to yield thresholds that are directly related to 0 and 298 K bond dissociation energies for M(+)-adenine after accounting for the effects of multiple ion-molecule collisions, kinetic and internal energy distributions of the reactants, and dissociation lifetimes. The measured bond energies are compared to those previously studied for simple nitrogen donor ligands, NH(3) and pyrimidine, and to results for alkali metal cations bound to adenine. Trends in these results and theoretical calculations on Cu(+)(adenine) suggest distinct differences in the binding site propensities of adenine to the alkali vs transition metal ions, a consequence of s-dsigma hybridization on the latter.     

Ab initio electron correlated studies on the intracluster reaction of NO+ (H2O)(n) → H3O+ (H2O)(n-2) (HONO) (n = 4 and 5)

Hydrated nitrosonium ion clusters NO(+)(H(2)O)(n) (n = 4 and 5) were investigated by using MP2/aug-cc-pVTZ level of theory to clarify isomeric reaction pathways for formation of HONO and fully hydrated hydride ions. We found some new isomers and transition state structures in each hydration number, whose lowest activation energies of the intracluster reactions were found to be 4.1 and 3.4 kcal mol(-1) for n = 4 and n = 5, respectively. These thermodynamic properties and full quantum mechanical molecular dynamics simulation suggest that product isomers with HONO and fully hydrated hydride ions can be obtained at n = 4 and n = 5 in terms of excess hydration binding energies which can overcome these activation barriers.     

Ions in size-selected aqueous nanodrops: sequential water molecule binding energies and effects of water on ion fluorescence

The effects of water on ion fluorescence were investigated, and average sequential water molecule binding energies to hydrated ions, M(z)(H(2)O)(n), at large cluster size were measured using ion nanocalorimetry. Upon 248-nm excitation, nanodrops with ~25 or more water molecules that contain either rhodamine 590(+), rhodamine 640(+), or Ce(3+) emit a photon with average energies of approximately 548, 590, and 348 nm, respectively. These values are very close to the emission maxima of the corresponding ions in solution, indicating that the photophysical properties of these ions in the nanodrops approach those of the fully hydrated ions at relatively small cluster size. As occurs in solution, these ions in nanodrops with 8 or more water molecules fluoresce with a quantum yield of ~1. Ce(3+) containing nanodrops that also contain OH(-) fluoresce, whereas those with NO(3)(-) do not. This indirect fluorescence detection method has the advantages of high sensitivity, and both the size of the nanodrops as well as their constituents can be carefully controlled. For ions that do not fluoresce in solution, such as protonated tryptophan, full internal conversion of the absorbed 248-nm photon occurs, and the average sequential water molecule binding energies to the hydrated ions can be accurately obtained at large cluster sizes. The average sequential water molecule binding energies for TrpH(+)(H(2)O)(n) and a doubly protonated tripeptide, [KYK + 2H](2+)(H(2)O)(n), approach asymptotic values of ~9.3 (n ≥ 11) and ~10.0 kcal/mol (n ≥ 25), respectively, consistent with a liquidlike structure of water in these nanodrops.     

Modulation of Stacking Interactions by Transition‐Metal Coordination: Ab Initio Benchmark Studies

A series of ab initio calculations are used to determine the C? H???π and π???π‐stacking interactions of aromatic rings coordinated to transition‐metal centres. Two model complexes have been employed, namely, ferrocene and chromium benzene tricarbonyl. Benchmark data obtained from extrapolation of MP2 energies to the basis set limit, coupled with CCSD(T) correction, indicate that coordinated aromatic rings are slightly weaker hydrogen‐bond acceptors but are significantly stronger hydrogen‐bond donors than uncomplexed rings. It is found that π???π stacking to a second benzene is stronger than in the free benzene dimer, especially in the chromium case. This is assigned, by use of energy partitioning in the local correlation method, to dispersion interactions between metal d and benzene π orbitals. The benchmark data is also used to test the performance of more efficient theoretical methods, indicating that spin‐component scaling of MP2 energies performs well in all cases, whereas various density functionals describe some complexes well, but others with errors of more than 1 kcal mol^?1.     

Characterization of ginseng saponins using electrospray mass spectrometry and collision-induced dissociation experiments of metal-attachment ions

Electrospray mass spectrometry (ESMS) and collision-induced dissociation (CID) methodologies have been developed for the structural characterization of ginseng saponins (ginsenosides). Ginsenosides are terpene glycosides containing a triterpene core to which one to four sugars may be attached. They are neutral molecules which readily form molecular metal-attachment ions in positive ion ESMS experiments. In the presence of ammonium hydroxide intense deprotonated ions are generated. Both positive and negative ion ESMS experiments were found to be useful for molecular mass and structure determination of ten ginsenoside standards. Negative ion experiments made possible the determination of the molecular mass of each ginsenoside standard, the mass of the triterpene core and the masses and sequences of the sugar residues. Positive ion ESMS experiments with the alkali metal cations Li+ or Na+ and the transition metal cations Co2+, Ni2+ and Zn2+ were also useful in determining molecular masses. These alkali and transition metal cations form strongly bonded attachment ions with the ginsenosides. As a result, the CID mass spectra of the metal attachment ions show a variety of (structure characteristic) fragmentations. These experiments can be used to determine the identity of the triterpene core, the types and attachment points of sugars to the core and the nature of the O-glycosidic linkages in the appended disaccharides. Combining the results from the negative and positive ion experiments provides a promising approach to the structure analysis of this class of natural products.     

Binding in transition metal complexes: Reduced multireference coupled-cluster study of the MCH2+ (M=Sc to Cu) compounds

The recently developed reduced multireference coupled-cluster method with singles and doubles (RMR CCSD), which is perturbatively corrected for triples [RMR CCSD(T)], is employed to compute binding energies of nine transition metal ions with CH2. Unlike analogous compounds involving main-group elements, the MCH2+ (M=Sc to Cu) transition metal complexes often exhibit a non-negligible multireference character. The authors thus employ the RMR CCSD(T) method, which represents an extension of the standard single-reference (SR) CCSD(T) method and can account for multireference effects, while employing only small reference spaces. In this way the role of quasidegeneracy effects on the binding energies of these complexes can be assessed at a higher SD(T) level than is possible with the widely used ab initio methods, namely, with the standard SR CCSD(T) approach, and provide a new benchmark for these quantities. The difference between the RMR and the standard CCSD(T) methods becomes particularly evident when considering nonequilibrium geometries.     

Structures and stabilities of 3d-transition metal-benzene organometallic clusters

In the gas phase, we have successfully synthesized organometallic clusters, M_n(benzene)_m (M=3d transition metal atoms), by using a laser vaporization method. The measurements of mass spectra and ionization energies (E_i) have revealed that the organometallic clusters can take two types of structures; layered sandwich structures (m = n + 1) and metal clusters saturatedly covered with benzenes. For early transition metals of Sc, Ti, and V, only the multiple decker sandwich structure clusters were preferentially produced, in which benzene and metal atoms are alternately piled up. For late transition metals of Co and Ni, the metal clusters saturatedly surrounded by benzenes were also produced as well as the sandwich clusters. Furthermore, the E_is of M₁(benzene)₂ (M = Sc-Ni) were systematically measured and their electronic properties will be discussed.     

Investigation of the metal binding site in methionine aminopeptidase by density functional theory

All methionine aminopeptidases exhibit the same conserved metal binding site. The structure of this site with either Co²⁺ions or Zn²⁺ions was investigated using density functional theory. The calculations showed that the structure of the site was not influenced by the identity of the metal ions. This was the case for both of the systems studied; one based on the X-ray structure of the human methionine aminopeptidase type 2 (hMetAP-2) and the other based on the X-ray structure of the E. colimethionine aminopeptidase type 1 (eMetAP-1). Another important structural issue is the identity of the bridging oxygen, which is part of either a water molecule or a hydroxide ion. Within the site of hMetAP-2 the results strongly indicate that a hydroxide ion bridges the metal ions. By contrast, the nature of the oxygen bridging the metal ions within the metal binding site of eMetAP-1 cannot be determined based on the results here, due to the similar structural results obtained with a bridging water molecule and a bridging hydroxide ion.     

Interference between the hydrogen bonds to the two rings of nicotine

The biochemical transport and binding of nicotine depends on the hydrogen bonding between water and binding site residues to the pyridine ring and the protonated pyrrolidinium ring. To test the independence of these two moderately separated hydrogen-bonding sites, we have calculated the structures of clusters of protonated nicotine with water and a bicarbonate anion, benzene, indole, or a second water molecule. Unprotonated nicotine-water clusters have also been studied for contrast. The potential energy surfaces are first explored with an intermolecular anisotropic atom-atom model potential. Full geometry optimizations are then carried out using density functional theory to include nonadditive terms in the interaction energies. The presence of the charge on the pyrrolidine nitrogen removes the conventional hydrogen-bonding site on the pyridine ring. The hydrogen-bond ability of this site is nearly recovered when the protonated pyrrolidinium ring is bound to a bicarbonate anion, whereas its interaction with benzene shows a much smaller effect. Indole appears to partially restore the hydrogen-bond ability of the pyridine nitrogen, although indole and benzene both pi-bond to the pyrrolidinium ring. A second hydrogen-bonding water produces a significant conformational distortion of the nicotine. This demonstrates the limitations of the conventional qualitative predictions of hydrogen bonding based on the independence of molecular fragments. It also provides benchmarks for the development of atomistic modeling of biochemical systems.     

Spectroscopy of benzene complexes with perylene and other aromatic species

The fluorescence excitation spectra are reported of complexes of benzene with perylene and anthracene. Some correlations are established between the results of simple potential energy calculations and the experimental data, in regard to the structures and binding energies of a number of different complexes. In particular, while the anthracenebenzene complex is somewhat similar to one postulated form of the benzene dimer where parallel, but displaced rings are found, benzene appears to occupy a center-of-mass position on perylene. The spectrum of perylenebenzene also shows evidence of strong coupling between internal modes of the two molecules near 1600 cm⁻¹. The other major perturbation of the spectrum involves damping of the out-of-plane “butterfly” motion of perylene by the adsorbed benzene molecule. The principal low-frequency mode, known to be a v = 0 → v = 2 transition in the bare molecule, at 95 cm⁻¹, is replaced in the benzene 1:1 and 2:1 complexes by a transition, at 68 cm⁻¹. Furthermore, unlike the bare molecule, the ground state frequency of the perturbed out-of-plane mode is very similar to that of the excited state. Indications from these data support the idea that the equilibrium out-of-plane distortion of perylene in a complex with benzene is rather different from that observed in the bare perylene molecule.     

Ultraviolet photodepletion spectroscopy of dibenzo-18-crown-6-ether complexes with alkaline Earth metal divalent cations

Ultraviolet photodepletion spectra of dibenzo-18-crown-6-ether complexes with alkaline earth metal divalent cations (A(2+)-DB18C6, A = Ba, Sr, Ca, and Mg) were obtained in the gas phase using electrospray ionization quadrupole ion-trap reflectron time-of-flight mass spectrometry. Each spectrum exhibits the lowest energy absorption band in the wavenumber region of 35?400-37?800 cm(-1), which is tentatively assigned as the origin of the S(0)-S(1) transition of A(2+)-DB18C6. This origin band shows a red shift as the size of the metal dication increases from Mg(2+) to Ba(2+). The binding energies of the metal dications to DB18C6 at the S(0) state were calculated at the lowest energy structures optimized by the density functional theory and employed with the experimental energies of the origin bands to estimate the binding energies at the S(1) state. We suggest that the red shifts of the origin bands arise from the decrease in the binding energies of the metal dications at the S(1) state by nearly constant ratios with respect to the binding energies at the S(0) state, which decrease with increasing size of the metal dication. This unique relationship of the binding energies between the S(0) and S(1) states gives rise to a linear correlation between the relative shift of the origin band of A(2+)-DB18C6 and the binding energy of the metal dication at the S(0) state. The size effects of the metal cations on the properties of metal-DB18C6 complex ions are also manifested in the linear plot of the relative shift of the origin band as a function of the size to charge ratio of the metal cations, where the shifts of the origin bands for all DB18C6 complexes with alkali and alkaline earth metal cations are fit to the same line.     

Bonding preferences of C2X4-bridged bimetallic transition metal complexes of Ti, Cu, and Ag

The bonding preference of transition metal species of general formula [(PH3)2M]2(mu-C2X4), where M = Cu or Ag and X = O, S, Se, or Te, and (Cp2Ti)2(mu-C2X4), where X = S or Se, are explored using density functional theory. The relative energies of metal binding to the bridging ligand in a dithiolene-like vs dithiocarbamate-like manner are evaluated. In all cases, the most stable structure corresponds to dithiolene-like (or side-side) bonding, consistent with the vast majority of these compounds which have been experimentally characterized. However, for M = Ag and X = S, Se, or Te, the two isomers are nearly degenerate.