Ground-state path integral Monte Carlo simulations of positive ions in (4)He clusters: bubbles or snowballs?

The local order around alkali (Li(+) and Na(+)) and alkaline-earth (Be(+), Mg(+), and Ca(+)) ions in (4)He clusters has been studied using ground-state path integral Monte Carlo calculations. The authors apply a criterion based on multipole dynamical correlations to discriminate between solidlike and liquidlike behaviors of the (4)He shells coating the ions. As it was earlier suggested by experimental measurements in bulk (4)He, their findings indicate that Be(+) produces a solidlike ("snowball") structure, similar to alkali ions and in contrast to the more liquidlike (4)He structure embedding heavier alkaline-earth ions.     

Ionization of doped helium nanodroplets: residual helium attached to diatomic cations and their clusters

Electron impact ionization of helium nanodroplets containing a dopant, M, can lead to the detection of both M(+) and helium-solvated cations of the type M(+)·He(n) in the gas phase. The observation of helium-doped ions, He(n)M(+), has the potential to provide information on the aftermath of the charge transfer process that leads to ion production from the helium droplet. Here we report on helium attachment to the ions from four common diatomic dopants, M = N(2), O(2), CO, and NO. For experiments carried out with droplets with an average size of 7500 helium atoms, the monomer cations show little tendency to attach and retain helium atoms on their journey out of the droplet. By way of contrast, the corresponding cluster cations, M(n)(+), where n ≥ 2, all show a clear affinity for helium and form He(m)M(n)(+) cluster ions. The stark difference between the monomer and cluster ions is attributed to more effective cooling of the latter in the aftermath of the ionization event.     

Hydration of alkali ions from first principles molecular dynamics revisited

Structural and dynamical properties of the hydration of Li(+), Na(+), and K(+) in liquid water at ambient conditions were studied by first principles molecular dynamics. Our simulations successfully captured the different hydration behavior shown by the three alkali ions as observed in experiments. The present analyses of the dependence of the self-diffusion coefficient and rotational correlation time of water on the ion concentration suggest that Li(+) (K(+)) is certainly categorized as a structure maker (breaker), whereas Na(+) acts as a weak structure breaker. An analysis of the relevant electronic structures, based on maximally localized Wannier functions, revealed that the dipole moment of H(2)O molecules in the first solvation shell of Na(+) and K(+) decreases by about 0.1 D compared to that in the bulk, due to a contraction of the oxygen lone pair orbital pointing toward the metal ion.     

Solvation of Na+, K+, and their dimers in helium

Helium atoms bind strongly to alkali cations which, when embedded in liquid helium, form so-called snowballs. Calculations suggest that helium atoms in the first solvation layer of these snowballs form rigid structures and that their number (n) is well defined, especially for the lighter alkalis. However, experiments have so far failed to accurately determine values of n. We present high-resolution mass spectra of Na(+)He(n), K(+)He(n), Na(2)(+)He(n) and K(2)(+)He(n), formed by electron ionization of doped helium droplets; the data allow for a critical comparison with several theoretical studies. For sodium and potassium monomers the spectra indicate that the value of n is slightly smaller than calculated. Na(2)(+)He(n) displays two distinct anomalies at n=2 and n=6, in agreement with theory; dissociation energies derived from experiment closely track theoretical values. K(2)(+)He(n) distributions are fairly featureless, which also agrees with predictions.     

Communication: the formation of helium cluster cations following the ionization of helium nanodroplets: influence of droplet size and dopant

The He(n)(+)/He(2)(+) (n ≥ 3) signal ratios in the mass spectra derived from electron impact ionization of pure helium nanodroplets are shown to increase with droplet size, reaching an asymptotic limit at an average droplet size of approximately 50,000 helium atoms. This is explained in terms of a charge hopping model, where on average the positive charge is able to penetrate more deeply into the liquid helium as the droplet size increases. The deeper the point where the charge localizes to form He(2)(+), the greater the likelihood of collisions with the surrounding helium as the ion begins to leave the droplet, thus increasing the probability that helium will be ejected in the form of He(n)(+) (n ≥ 3) cluster ions rather than He(2)(+). The addition of a dopant alters the He(n)(+)/He(2)(+) ratio for small helium droplets, an observation attributed to the potential energy gradient created by the cation-dopant interaction and its effect in drawing the positive charge towards the dopant in the interior of the droplet.     

Non-covalent interactions of alkali metal cations with singly charged tryptic peptides

The complexes formed by alkali metal cations (Cat(+) = Li(+), Na(+), K(+), Rb(+)) and singly charged tryptic peptides were investigated by combining results from the low-energy collision-induced dissociation (CID) and ion mobility experiments with molecular dynamics and density functional theory calculations. The structure and reactivity of [M + H + Cat](2+) tryptic peptides is greatly influenced by charge repulsion as well as the ability of the peptide to solvate charge points. Charge separation between fragment ions occurs upon dissociation, i.e. b ions tend to be alkali metal cationised while y ions are protonated, suggesting the location of the cation towards the peptide N-terminus. The low-energy dissociation channels were found to be strongly dependant on the cation size. Complexes containing smaller cations (Li(+) or Na(+)) dissociate predominantly by sequence-specific cleavages, whereas the main process for complexes containing larger cations (Rb(+)) is cation expulsion and formation of [M + H](+). The obtained structural data might suggest a relationship between the peptide primary structure and the nature of the cation coordination shell. Peptides with a significant number of side chain carbonyl oxygens provide good charge solvation without the need for involving peptide bond carbonyl groups and thus forming a tight globular structure. However, due to the lack of the conformational flexibility which would allow effective solvation of both charges (the cation and the proton) peptides with seven or less amino acids are unable to form sufficiently abundant [M + H + Cat](2+) ion. Finally, the fact that [M + H + Cat](2+) peptides dissociate similarly as [M + H](+) (via sequence-specific cleavages, however, with the additional formation of alkali metal cationised b ions) offers a way for generating the low-energy CID spectra of 'singly charged' tryptic peptides.     

Ionic dimers in He droplets: interaction potentials for Li2(+)-He,Na2(+)-He, and K2(+)-He and stability of the smaller clusters

We present post Hartree-Fock calculations of the potential energy surfaces (PESs) for the ground electronic states of the three alkali dimer ions Li(2) (+),Na(2) (+), and K(2) (+) interacting with neutral helium. The calculations were carried out for the frozen molecular equilibrium geometries and for an extensive range of the remaining two Jacobi coordinates, R and theta, for which a total of about 1000 points is generated for each surface. The corresponding raw data were then fitted numerically to produce analytic expressions for the three PESs, which were in turn employed to evaluate the bound states of the three trimers for their J=0 configurations: The final spatial features of such bound states are also discussed in detail. The possible behavior of additional systems with more helium atoms surrounding the ionic dopants is gleaned from further calculations on the structural stability of aggregates with up to six He atoms. The validity of a sum-of-potential approximation to yield realistic total energies of the smaller cluster is briefly discussed vis-a-vis the results from many-body calculations.     

Structure and energetics of poly(ethylene glycol) cationized by Li(+), Na(+), K(+) and Cs(+): a first-principles study

Density functional theoretical methods, including several basis sets and two functional, were used to collect information on the structure and energetic parameters of poly(ethylene glycol) (PEG), also referred to as poly(ethylene oxide) (PEO), coordinated by alkali metal ions. The oligomer chain is found to form a spiral around the alkali cation, which grows to roughly two helical turns when the oligomer size increases to about the decamer for each alkali ion. Above this size, the additional monomer units do not build the spiral further for Li(+) and Na(+); instead, they form less organized segments outside or next to the initial spiral. The distance of the first layer of co-ordinating O atoms from the alkali cation is 1.9-2.15 ? for Li(+), 2.3-2.5 ? for Na(+), 2.75-3.2 ? for K(+) and 3.5-3.8 ? for Cs(+) complexes. The number of O atoms in the innermost shell is five, six, seven and eleven for Li(+), Na(+), K(+) and Cs(+). The collision cross sections with He increase linearly with the oligomer to a very good approximation. No sign of leaning towards the 2/3 power dependence characterizing spherical particles is observed. The binding energy of the cation to the oligomer increases up to polymerization degree of about 10, where it levels off for each alkali-metal ion, indicating that this is approximately the limit of the oligomer size that can be influenced by the alkali cation. The binding energy-degree of polymerization curves are remarkably parallel for the four cations. The limiting binding energy at large polymerization degrees is about 544 kJ mol(-1), 460 kJ mol(-1), 356 kJ mol(-1) and 314 kJ mol(-1) for Li, Na, K and Cs, respectively. The geometrical features are compared with the X-ray and neutron diffraction data on crystalline and amorphous phases of conducting polymers formed by alkali-metal salts and PEG. The implications of the observations concerning collision cross sections and binding energies to ion mobility spectroscopy and mass spectrometry are discussed.     

Fragmentation of ionized doped helium nanodroplets: theoretical evidence for a dopant ejection mechanism

We report a theoretical study of the effect induced by a helium nanodroplet environment on the fragmentation dynamics of a dopant. The dopant is an ionized neon cluster Ne(n) (+) (n=4-6) surrounded by a helium nanodroplet composed of 100 atoms. A newly designed mixed quantum/classical approach is used to take into account both the large helium cluster zero-point energy due to the light mass of the helium atoms and all the nonadiabatic couplings between the Ne(n) (+) potential-energy surfaces. The results reveal that the intermediate ionic dopant can be ejected from the droplet, possibly with some helium atoms still attached, thereby reducing the cooling power of the droplet. Energy relaxation by helium atom evaporation and dissociation, the other mechanism which has been used in most interpretations of doped helium cluster dynamics, also exhibits new features. The kinetic energy distribution of the neutral monomer fragments can be fitted to the sum of two Boltzmann distributions, one with a low kinetic energy and the other with a higher kinetic energy. This indicates that cooling by helium atom evaporation is more efficient than was believed so far, as suggested by recent experiments. The results also reveal the predominance of Ne(2) (+) and He(q)Ne(2) (+) fragments and the absence of bare Ne(+) fragments, in agreement with available experimental data (obtained for larger helium nanodroplets). Moreover, the abundance in fragments with a trimeric neon core is found to increase with the increase in dopant size. Most of the fragmentation is achieved within 10 ps and the only subsequent dynamical process is the relaxation of hot intermediate He(q)Ne(2) (+) species to Ne(2) (+) by helium atom evaporation. The dependence of the ionic fragment distribution on the parent ion electronic state reached by ionization is also investigated. It reveals that He(q)Ne(+) fragments are produced only from the highest electronic state, whereas He(q)Ne(2) (+) fragments originate from all the electronic states. Surprisingly, the highest electronic states also lead to fragments that still contain the original ionic dopant species. A mechanism is conjectured to explain this fragmentation inhibition.     

Alkali metal ion binding to amino acids versus their methyl esters: affinity trends and structural changes in the gas phase

The relative alkali metal ion (M(+)) affinities (binding energies) between seventeen different amino acids (AA) and the corresponding methyl esters (AAOMe) were determined in the gas phase by the kinetic method based on the dissociation of AA-M(+)-AAOMe heterodimers (M=Li, Na, K, Cs). With the exception of proline, the Li(+), Na(+), and K(+) affinities of the other aliphatic amino acids increase in the order AAAAOMe is already observed for K(+). Proline binds more strongly than its methyl ester to all M(+) except Li(+). Ab initio calculations on the M(+) complexes of alanine, beta-aminoisobutyric acid, proline, glycine methyl ester, alanine methyl ester, and proline methyl ester show that their energetically most favorable complexes result from charge solvation, except for proline which forms salt bridges. The most stable mode of charge solvation depends on the ligand (AA or AAOMe) and, for AA, it gradually changes with metal ion size. Esters chelate all M(+) ions through the amine and carbonyl groups. Amino acids coordinate Li(+) and Na(+) ions through the amine and carbonyl groups as well, but K(+) and Cs(+) ions are coordinated by the O atoms of the carboxyl group. Upon consideration of these differences in favored binding geometries, the theoretically derived relative M(+) affinities between aliphatic AA and AAOMe are in good overall agreement with the above given experimental trends. The majority of side chain functionalized amino acids studied show experimentally the affinity order AAAAOMe. The latter ranking is attributed to salt bridge formation.     

Mg impurity in helium droplets

Within the diffusion Monte Carlo approach, we have determined the structure of isotopically pure and mixed helium droplets doped with one magnesium atom. For pure (4)He clusters, our results confirm those of Mella et al. [J. Chem. Phys. 123, 054328 (2005)] that the impurity experiences a transition from a surface to a bulk location as the number of helium atoms in the droplet increases. Contrarily, for pure (3)He clusters Mg resides in the bulk of the droplet due to the smaller surface tension of this isotope. Results for mixed droplets are presented. We have also obtained the absorption spectrum of Mg around the 3s3p?(1)P(1) ← 3s(2)?(1)S(0) transition.     

Alkali-ion microsolvation with benzene molecules

The target of this investigation is to characterize by a recently developed methodology, the main features of the first solvation shells of alkaline ions in nonpolar environments due to aromatic rings, which is of crucial relevance to understand the selectivity of several biochemical phenomena. We employ an evolutionary algorithm to obtain putative global minima of clusters formed with alkali-ions (M(+)) solvated with n benzene (Bz) molecules, i.e., M(+)-(Bz)(n). The global intermolecular interaction has been decomposed in Bz-Bz and in M(+)-Bz contributions, using a potential model based on different decompositions of the molecular polarizability of benzene. Specifically, we have studied the microsolvation of Na(+), K(+), and Cs(+) with benzene molecules. Microsolvation clusters up to n = 21 benzene molecules are involved in this work and the achieved global minimum structures are reported and discussed in detail. We observe that the number of benzene molecules allocated in the first solvation shell increases with the size of the cation, showing three molecules for Na(+) and four for both K(+) and Cs(+). The structure of this solvation shell keeps approximately unchanged as more benzene molecules are added to the cluster, which is independent of the ion. Particularly stable structures, so-called "magic numbers", arise for various nuclearities of the three alkali-ions. Strong "magic numbers" appear at n = 2, 3, and 4 for Na(+), K(+), and Cs(+), respectively. In addition, another set of weaker "magic numbers" (three per alkali-ion) are reported for larger nuclearities.     

UV and IR spectroscopic studies of cold alkali metal ion-crown ether complexes in the gas phase

We report UV photodissociation (UVPD) and IR-UV double-resonance spectra of dibenzo-18-crown-6 (DB18C6) complexes with alkali metal ions (Li(+), Na(+), K(+), Rb(+), and Cs(+)) in a cold, 22-pole ion trap. All the complexes show a number of vibronically resolved UV bands in the 36,000-38,000 cm(-1) region. The Li(+) and Na(+) complexes each exhibit two stable conformations in the cold ion trap (as verified by IR-UV double resonance), whereas the K(+), Rb(+), and Cs(+) complexes exist in a single conformation. We analyze the structure of the conformers with the aid of density functional theory (DFT) calculations. In the Li(+) and Na(+) complexes, DB18C6 distorts the ether ring to fit the cavity size to the small diameter of Li(+) and Na(+). In the complexes with K(+), Rb(+), and Cs(+), DB18C6 adopts a boat-type (C(2v)) open conformation. The K(+) ion is captured in the cavity of the open conformer thanks to the optimum matching between the cavity size and the ion diameter. The Rb(+) and Cs(+) ions sit on top of the ether ring because they are too large to enter the cavity of the open conformer. According to time-dependent DFT calculations, complexes that are highly distorted to hold metal ions open the ether ring upon S(1)-S(0) excitation, and this is confirmed by extensive low-frequency progressions in the UVPD spectra.     

Solvent effects on the conformation of DNA dodecamer segment: a simulation study

Different solvent temperatures with five kinds of counterions are used to investigate solvent effects on the DNA microscopic structure. The dodecamer d (CGCGAATTCGCG) DNA segment is merged into the solvents and its conformation transition is studied with the molecular dynamics simulations in detail. For the simple point charge model of water molecule with Na(+) counterions, as temperature increases from 200 K to 343 K, the duplex DNA changes from stiff B form to a state between A form and B form, which we define as mixed (A-B) structure, with a double helix unwinding. To study the counterions effects, other four alkali cations, Li(+), K(+), Rb(+), or Cs(+) ions, are substituted for Na(+) ions at 298 K and 343 K, respectively. For the cases of Li(+), Rb(+), and Cs(+) ions, the duplex DNA becomes more flexible with sugar configuration changing form C2'-endo to C1'-endo type and the width and depth of minor groove at CpG and GpC steps moving towards A values, as the mass of the counterions decreasing. For the case of K(+) ions, DNA-K(+) interaction widens the width of minor and major grooves at ApA steps and TpT steps, respectively. It seems that the light ions (Li(+) or Na(+)) prefer to interact with the free phosphate oxygen atoms while the heavier ions (Rb(+) and Cs(+)) strongly interact with the base pairs.     

The Influence of Alkali Metal Ions on Homogeneous Nucleation of Al(OH)(3) Crystals from Supersaturated Caustic Aluminate Solutions

Homogeneous nucleation of Al(OH)(3) crystals from synthetic, optically clear, caustic aluminate solutions and the influence of alkali metal ion (Na(+) versus K(+)) have been investigated under isothermal, batch crystallization conditions. The nucleation kinetics showed a seventh-order dependence upon Al(III) relative supersaturation and a strong temperature effect. Activation energy of 160 kJ mol(-1) and interfacial energy of 33 mJ m(-2) were estimated and found to be independent of alkali ion, as was the Al(OH)(3) equilibrium solubility. The nucleation rate, however, was faster in aging sodium than in potassium aluminate solutions. It appears that Na(+) ions provide greater stability for the formation and densification of Al(III)-containing, supramolecular clusters which grow more rapidly into Al(OH)(3) crystallites than do K(+) ions. The development of the Al-OH octahedral structure of Al(OH)(3) nuclei is an alkali metal ion-mediated, chemical reaction-controlled condensation process, displaying specific gibbsite (gamma-Al(OH)(3))-bayerite (alpha-Al(OH)(3)) dimorphism. Furthermore, significant differences in the level of alkali ion incorporation, reflecting in the purity and morphology of the crystalline product, were observed. Copyright 2000 Academic Press.     

Microsolvation of an ionic dopant in small (4)He clusters: OH(+)((3)sigma)((4)He)(N) via genetic algorithm optimizations.

The optimized spatial structures of the small clusters (with N up to 33) formed by an increasing number of (4)He atoms, which act as a microsolvent surrounding the OH(+) ionic molecular dopant, are obtained using a sum-of-potentials scheme corrected by three-body (3B) effects. The most stable structures are generated using the type of genetic algorithm described herein, and the sequential formation of regular shell structures is analyzed in detail. Possible quantum corrections for both the solvent distributions and the stable energetics are analyzed and discussed.     

On the influence of ion exchange on the local structure of the titanosilicate ETS-10

The effect of ion exchange with different monovalent cations (NH(4)(+), K(+), Na(+) and Cs(+)) on the local structure of the titanosilicate ETS-10 has been studied by (29)Si MAS NMR and Raman spectroscopy. Although X-ray diffraction shows no significant influence of ion exchange on the long range order, ammonium exchange is found to result in substantial damage to the local structure. Ion exchange experiments with alkali cations under significantly more acidic conditions clearly show that the structural damage brought about by ammonium exchange is not caused by the low pH of the exchange solution. The exchange with potassium and caesium ions also leads to significant changes in the (29)Si NMR and Raman spectra. However, these changes can largely be reversed by sodium back-exchange.     

Atomic energy levels and oscillator strengths calculated with a screened potential

A screened potential model (SPM ) is formulated with Slater-type functions to reproduce the density of the electron shells. The orbital exponent for the valence shell is optimized to calibrate the SPM against experimental valence and Rydberg s and p energy levels for the alkali metals (Li, Na, and K) rare gases (He, Ne, Ar, and Kr), second-row atoms, and Cu, Zn, Mg, S, and Ga. The resulting one-electron wave functions are used to calculate the dipole and velocity forms of the oscillator strengths and 〈1/r³〉 for spin-orbit coupling. The excellent agreement with the averaged experimental results suggests that the SPM atomic orbitals are a good starting point for the evaluation of one-electron properties.     

Investigation of homo- and heterodimer alkali metal cation complexes of resorc[4]arenes by electrospray ionization mass spectrometry

Resorc[4]arenes are compounds with interesting properties, mainly because of their ability to form host-guest complexes with the guest located inside the cavity. The size of the guest limits the complexation, as shown by a competition experiment with tetraalkylammonium ions of different size. By electroscopy ionization tandem mass spectrometric experiments on resorc[4]arene heterodimers bearing an alkali metal ion as guest, it was found that there must be two different binding mechanisms for alkali metal ions with high surface charge density (Li(+) and Na(+)) on the one hand compared with those with a lower surface charge density on the other hand (K(+), Rb(+), Cs(+)).