Ar solvation shells in K(+)-HFBz: from cluster rearrangement to solvation dynamics

The effect of some leading intermolecular interaction components on specific features of weakly bound clusters involving an aromatic molecule, a closed shell ion, and Ar atoms is analyzed by performing molecular dynamics simulations on potential energy surfaces properly formulated in a consistent way. In particular, our investigation focuses on the three-dimensional Ar distributions around the K(+)-hexafluorobenzene (K(+)-HFBz) dimer, in K(+)-HFBz-Ar(n) aggregates (n ≤ 15), and on the gradual evolution from cluster rearrangement to solvation dynamics when ensembles of 50, 100, 200, and 500 Ar atoms are taken into account. Results indicate that the Ar atoms compete to be placed in such a way to favor an attractive interaction with both K(+) and HFBz, occupying positions above and below the aromatic plane but close to the cation. When these positions are already occupied, the Ar atoms tend to be placed behind the cation, at larger distances from the center of mass of HFBz. Accordingly, three different groups of Ar atoms are observed when increasing n, with two of them surrounding K(+), thus, disrupting the K(+)-HFBz equilibrium geometry and favoring the dissociation of the solvated cation when the temperature increases. The selective role of the leading intermolecular interaction components directly depending on the ion size repulsion is discussed in detail by analyzing similarities and differences on the behavior of the Ar-solvated K(+)-HFBz and Cl(-)-Bz aggregates.     

Isomerization dynamics and thermodynamics of ionic argon clusters

The dynamics and thermodynamics of small Ar(n) (+) clusters, n=3, 6, and 9, are investigated using molecular dynamics (MD) and exchange Monte Carlo (MC) simulations. A diatomic-in-molecule Hamiltonian provides an accurate model for the electronic ground state potential energy surface. The microcanonical caloric curves calculated from MD and MC methods are shown to agree with each other, provided that the rigorous conservation of angular momentum is accounted for in the phase space density of the MC simulations. The previously proposed projective partition of the kinetic energy is used to assist MD simulations in interpreting the cluster dynamics in terms of inertial, internal, and external modes. The thermal behavior is correlated with the nature of the charged core in the cluster by computing a dedicated charge localization order parameter. We also perform systematic quenches to establish a connection with the various isomers. We find that the Ar(3) (+) cluster is very stable in its linear ground state geometry up to about 300 K, and then isomerizes to a T-shaped isomer in which a quasineutral atom lies around a charged dimer. In Ar(6) (+) and Ar(9) (+), the covalent trimer core is solvated by neutral atoms, and the weakly bound solvent shell melts at much lower energies, occasionally leading to a tetramer or pentamer core with weakly charged extremities. At high energies the core itself becomes metastable and the cluster transforms into Ar(2) (+) solvated by a fluid of neutral argon atoms.     

Solvation effects on the intracluster elimination channels in M+(L)n, where M+= Mg+ and Ca+, L = CH3OH, and NH3, and n = 2-6

The methanol and ammonia solvated Ca (+) or Mg (+) clusters are known to go through intracluster H or CH 3 eliminations which are typically switched on just below n = 6. By first principles calculations at the B3LYP/6-311+G** level, we have identified the transition structures, activation barriers, and energy changes in these reactions for clusters with 2-6 solvent molecules. The activation barrier is crucial to explain the previously reported experimental results. While increasing number of solvent molecules stabilizes a transition structure, the increasing presence of solvent molecules in the first solvation shell makes it difficult for the metal ion to assist the bond breaking through its interaction with the departing H atom or CH 3 group. The balance of these two factors determines whether a particular elimination channel could be switched on.     

Size-specific interaction of alkali metal ions in the solvation of M+-benzene clusters by Ar atoms

The size-specific influence of the M+ alkali ion (M = Li, Na, K, Rb, and Cs) in the solvation process of the M+-benzene clusters by Ar atoms is investigated by means of molecular dynamic simulations. To fully understand the behavior observed in M+-bz-Ar(n) clusters, solvation is also studied in clusters containing either M+ or benzene only. The potential energy surfaces employed are based on a semiempirical bond-atom decomposition, which has been developed previously by some of the authors. The outcome of the dynamics is analyzed by employing radial distribution functions, studying the evolution of the distances between the Ar atoms and the alkali ion M+ or the benzene molecule for all M+-bz-Ar(n) clusters. For all members, in the M+-bz series, the benzene molecule (bz) is found to remain strongly bound to M+ even in the presence of solvent atoms. The radial distribution functions for the heavier clusters (K+-bz, Rb+-bz, and Cs+-bz), are found to be different than for the lighter (Na+-bz and Li+-bz) ones.     

Structure and photoabsorption properties of cationic alkali dimers solvated in neon clusters

We present a theoretical investigation of the structure and optical absorption of M(2)(+) alkali dimers (M=Li,Na,K) solvated in Ne(n) clusters for n=1 to a few tens Ne atoms. For all these alkali, the lowest-energy isomers are obtained by aggregation of the first Ne atoms at the extremity of the alkali molecule. This particular geometry, common to other M(2)(+)-rare gas clusters, is intimately related to the shape of the electronic density of the X (2)Σ(g)(+) ground state of the bare M(2)(+) molecules. The structure of the first solvation shell presents equilateral Ne(3) and capped pentagonal Ne(6) motifs, which are characteristic of pure rare gas clusters. The size and geometry of the complete solvation shell depend on the alkali and were obtained at n=22 with a D(4h) symmetry for Li and at n=27 with a D(5h) symmetry for Na. For K, our study suggests that the closure of the first solvation shell occurs well beyond n=36. We show that the atomic arrangement of these clusters has a profound influence on their optical absorption spectrum. In particular, the XΣ transition from the X (2)Σ(g)(+) ground state to the first excited (2)Σ(u)(+) state is strongly blueshifted in the Frank-Condon area.     

Disparate effects of Cu and V on structures of exohedral transition metal-doped silicon clusters: a combined far-infrared spectroscopic and computational study

The growth mechanisms of small cationic silicon clusters containing up to 11 Si atoms, exohedrally doped by V and Cu atoms, are described. We find that as dopants, V and Cu follow two different paths: while V prefers substitution of a silicon atom in a highly coordinated position of the cationic bare silicon clusters, Cu favors adsorption to the neutral or cationic bare clusters in a lower coordination site. The different behavior of the two transition metals becomes evident in the structures of Si(n)M(+) (n = 4-11 for M = V, and n = 6-11 for M = Cu), which are investigated by density functional theory and, for several sizes, confirmed by comparison with their experimental vibrational spectra. The spectra are measured on the corresponding Si(n)M(+)·Ar complexes, which can be formed for the exohedrally doped silicon clusters. The comparison between experimental and calculated spectra indicates that the BP86 functional is suitable to predict far-infrared spectra of these clusters. In most cases, the calculated infrared spectrum of the lowest-lying isomer fits well with the experiment, even when various isomers and different electronic states are close in energy. However, in a few cases, namely Si(9)Cu(+), Si(11)Cu(+), and Si(10)V(+), the experimentally verified isomers are not the lowest in energy according to the density functional theory calculations, but their structures still follow the described growth mechanism. The different growth patterns of the two series of doped Si clusters reflect the role of the transition metal's 3d orbitals in the binding of the dopant atoms.     

Formation and photodissociation of M(+)-C(6)H(6) (M(+) = V(+) and Ta(+)) and Ta(+)-C(6)H(4) complexes in a time-of-flight mass spectrometer

A series of cyclic hydrocarbons were introduced to react with V(+) and Ta(+) using a pulsed beam expansion source in a time-of-flight mass spectrometer. The third-row metal Ta(+) displayed high reactivity in dehydrogenation to form benzyne complexes, whereas benzene complexes were the terminal products for V(+). M(+)-C(6)H(6) (M(+) = V(+) and Ta(+)) and Ta(+)-C(6)H(4) were selected to perform the photodissociation experiments. In contrast to the V(+) fragment formation via simple cleavage of the V(+)-C(6)H(6) bond, a photoinduced loss of C(2)H(2) occurred in both the Ta(+)-C(6)H(6) and Ta(+)-C(6)H(4) complexes. Plausible explanations involved in the formation of Ta(+)-C(6)H(6) and Ta(+)-C(6)H(4) complexes are given for observing such photo-induced dissociation. The observed photodissociation in Ta(+)-C(6)H(6) is analogous to the dissociative process previously investigated in metal ion-molecule reactions. The photodissociation spectrum of Ta(+)-C(6)H(4) was obtained by recording the appearance of Ta(+)-C(4)H(2) as a function of wavelength and yielded a dissociation energy of 91 +/- 1 kcal mol(-1).     

Theoretical studies on the electronic structure, charge distribution and vibrational spectra of diglyme-M(+)-AsF(6)(-) (M=Li, Na, K)

Electronic structure and the vibrational spectra of CH(3)(OCH(2)CH(2))(2)OCH(3)-M(+)-AsF(6)(-) (M=Li, Na, K) have been obtained using the density functional theory. Lithium ion exhibits a pentavalent coordination via 3 oxygens from diglyme and two fluorines of AsF(6)(-) whereas Na(+) and K(+) exhibit coordinate number 6 with 3 fluorines of the anion binding to alkali metal in these complexes. Analysis of calculated spectra reveal that the CH(2) wag (840-1120 cm(-1)) vibrations in the complex are sensitive to metal ion coordination. A frequency downshift relative to the free anion has been predicted for the vibrations of AsF(6)(-) anion when the fluorines are directly bonded (denoted by F) to metal ion. Consequent reorganization of electron density in the complex engenders a frequency shift in the opposite direction for As-F vibrations wherein the fluorine atoms are not coordinating to the alkali metal ion. An approach based on the molecular electron density topography coupled with the difference electron density map explains the direction of the frequency shifts of C-O-C and the As-F stretchings compared to those of free diglyme or AsF(6) anion. A new method, which includes the color-mapping function for the difference molecular electron density (MED), superimposed on the bond critical points in MED topography has been suggested to explain the direction of the frequency shifts in a single attempt.     

Negative ions of transition metal-halogen clusters

A systematic density functional theory based study of the structure and spectroscopic properties of neutral and negatively charged MX(n) clusters formed by a transition metal atom M (M=Sc,Ti,V) and up to seven halogen atoms X (X=F,Cl,Br) has revealed a number of interesting features: (1) Halogen atoms are bound chemically to Sc, Ti, and V for n≤n(max), where the maximal valence n(max) equals to 3, 4, and 5 for Sc, Ti, and V, respectively. For n>n(max), two halogen atoms became dimerized in the neutral species, while dimerization begins at n=5, 6, and 7 for negatively charged clusters containing Sc, Ti, and V. (2) Magnetic moments of the transition metal atoms depend strongly on the number of halogen atoms in a cluster and the cluster charge. (3) The number of halogen atoms that can be attached to a metal atom exceeds the maximal formal valence of the metal atom. (4) The electron affinities of the neutral clusters abruptly rise at n=n(max), reaching values as high as 7 eV. The corresponding anions could be used in the synthesis of new salts, once appropriate counterions are identified.     

Size resolved infrared spectroscopy of Na(CH3OH)n (n = 4-7) clusters in the OH stretching region: unravelling the interaction of methanol clusters with a sodium atom and the emergence of the solvated electron

Size resolved IR action spectra of neutral sodium doped methanol clusters have been measured using IR excitation modulated photoionisation mass spectroscopy. The Na(CH(3)OH)(n) clusters were generated in a supersonic He seeded expansion of methanol by subsequent Na doping in a pick-up cell. A combined analysis of IR action spectra, IP evolutions and harmonic predictions of IR spectra (using density functional theory) of the most stable structures revealed that for n = 4, 5 structures with an exterior Na atom showing high ionisation potentials (IPs) of ~4 eV dominate, while for n = 6, 7 clusters with lower IPs (~3.2 eV) featuring fully solvated Na atoms and solvated electrons emerge and dominate the IR action spectra. For n = 4 simulations of photoionisation spectra using an ab initio MD approach confirm the dominance of exterior structures and explain the previously reported appearance IP of 3.48 eV by small fractions of clusters with partly solvated Na atoms. Only for this cluster size a shift in the isomer composition with cluster temperature has been observed, which may be related to kinetic stabilisation of less Na solvated clusters at low temperatures. Features of slow fragmentation dynamics of cationic Na(+)(CH(3)OH)(6) clusters have been observed for the photoionisation near the adiabatic limit. This finding points to the relevance of previously proposed non-vertical photoionisation dynamics of this system.     

Static and dynamic properties of anionic intermolecular aggregates: the I−–benzene–Ar n case

Molecular dynamics simulations on the I^?–benzene–Ar _n clusters have been carried out using an atom(ion)-bond model to describe the nonelectrostatic contribution to the total interaction. Results for I^?–benzene–Ar and I^?–benzene–Ar _n (n = 3, 18 and 25) are presented and some predicted properties are compared with those of the alkali cation–benzene clusters solvated by Ar atoms.     

Unimolecular reactions of dihydrated alkaline earth metal dications M2+(H2O)2, M = Be, Mg, Ca, Sr, and Ba: salt-bridge mechanism in the proton-transfer reaction M2+(H2O)2 --> MOH+ + H3O

The unimolecular reactivity of M(2+)(H(2)O)(2), M = Be, Mg, Ca, Sr, and Ba, is investigated by density functional theory. Dissociation of the complex occurs either by proton transfer to form singly charged metal hydroxide, MOH(+), and protonated water, H(3)O(+), or by loss of water to form M(2+)(H(2)O) and H(2)O. Charge transfer from water to the metal forming H(2)O(+) and M(+)(H(2)O) is not favorable for any of the metal complexes. The relative energetics of these processes are dominated by the metal dication size. Formation of MOH(+) proceeds first by one water ligand moving to the second solvation shell followed by proton transfer to this second-shell water molecule and subsequent Coulomb explosion. These hydroxide formation reactions are exothermic with activation energies that are comparable to the water binding energy for the larger metals. This results in a competition between proton transfer and loss of a water molecule. The arrangement with one water ligand in the second solvation shell is a local minimum on the potential energy surface for all metals except Be. The two transition states separating this intermediate from the reactant and the products are identified. The second transition state determines the height of the activation barrier and corresponds to a M(2+)-OH(-)-H(3)O(+) "salt-bridge" structure. The computed B3LYP energy of this structure can be quantitatively reproduced by a simple ionic model in which Lewis charges are localized on individual atoms. This salt-bridge arrangement lowers the activation energy of the proton-transfer reaction by providing a loophole on the potential energy surface for the escape of H(3)O(+). Similar salt-bridge mechanisms may be involved in a number of proton-transfer reactions in small solvated metal ion complexes, as well as in other ionic reactions.     

Solvation of Na+ in argon clusters

The structures and stabilities of Ar(n)Na+ clusters (n < or = 54) are investigated using atomistic potentials fitted to reproduce ab initio calculations performed at the coupled-cluster level on the smaller clusters. Polarization effects are described using either the interaction between dipoles induced by the sodium ion, or a small charge transfer in the framework of a fluctuating charges model. In both models, extra three-body contributions of the Axilrod-Teller type are also included between the sodium ion and all pairs of argon atoms. The two models predict essentially similar growth patterns, and a transition in the metal ion coordination from 8 (square antiprism) to 12 (icosahedron) is seen to occur near n = 50, in response to the intrasolvent constraints.     

Analogy of the coordination chemistry of alkaline earth metal and lanthanide Ln(2+) ions: the isostructural zoo of mixed metal cages [IM(OtBu)4{Li(thf)}4(OH)] (M=Ca, Sr, Ba, Eu), [MM'6(OPh)8(thf)6] (M=Ca, Sr, Ba, Sm, Eu, M'=Li, Na), and their derivatives with 1,2-dimethoxyethane

As previously shown, alkali and alkaline earth metal iodides in nonaqueous, aprotic solvents behave like transition metal halides, forming cis- and trans-dihalides with various neutral O-donor ligands. These compounds can be used as precursors for the synthesis of new mixed alkali/alkaline earth metal aggregates. We show here that Ln2+ ions form isostructural cluster compounds. Thus, with LiOtBu, 50% of the initial iodide can be replaced in MI2, M=Ca, Sr, Ba, Eu, to generate the mixed-metal alkoxide aggregates [IM(OtBu)4{Li(thf)}4(OH)], for which the M--OH contacts were investigated by theoretical methods. With M'OPh (M'=Li, Na), a new mixed-metal aryloxide cluster type [MM'6(OPh)8(thf)6] is obtained for M=Ca, Sr, Ba, Sm, Eu. Their stability versus DME (DME=1,2-dimethoxyethane) as bidentate ligand is studied.     

Estimating the hydration enthalpies of neutral alkali metal atoms

Using existing data on the ionization energies of alkali metal atoms in small clusters of water, a thermodynamic cycle is proposed from which the hydration enthalpies of the neutral metal atoms can be estimated. Where comparisons are possible, the results are in reasonable agreement with those obtained using both experimental and ab initio methods. Application of the thermodynamic cycle to neutral alkali metal atoms solvated in ammonia yields solvation enthalpies that are significantly lower than those obtained for water.     

Isotope effects, dynamics, and the mechanism of solvolysis of aryldiazonium cations in water

The mechanism of the heterolytic solvolysis of p-tolyldiazonium cation in water was studied by a combination of kinetic isotope effects, theoretical calculations, and dynamics trajectories. Significant (13)C kinetic isotope effects were observed at the ipso (k(12)C/k(13)C = 1.024), ortho (1.017), and meta (1.013) carbons, indicative of substantial weakening of the C(2)-C(3) and C(5)-C(6) bonds at the transition state. This is qualitatively consistent with a transition state forming an aryl cation, but on a quantitative basis, simple S(N)1 heterolysis does not account best for the isotope effects. Theoretical S(N)2Ar transition structures for concerted displacement of N(2) by a single water molecule lead to poor predictions of the experimental isotope effects. The best predictions of the (13)C isotope effects arose from transition structures for the heterolytic process solvated by clusters of water molecules. These structures, formally saddle points for concerted displacements on the potential energy surface, may be described as transition structures for solvent reorganization around the aryl cation. Quasiclassical dynamics trajectories starting from these transition structures afforded products very slowly, compared to a similar S(N)2 displacement, and the trajectories often afforded long-lived aryl cation intermediates. Critical prior evidence for aryl cation intermediates is reconsidered with the aid of DFT calculations. Overall, the nucleophilic displacement process for aryldiazonium ions in water is at the boundary between S(N)2Ar and S(N)1 mechanisms, and an accurate view of the reaction mechanism requires consideration of dynamic effects.     

Relativistic-consistent electron densities of the coinage metal clusters M2, M4, M4(2-), and M4Na2 (M = Cu, Ag, Au): a QTAIM study

We employ second-order M?ller-Plesset perturbation theory level in combination with recently developed pseudopotential-based correlation consistent basis sets to obtain accurate relativistic-consistent electron densities for small coinage metal clusters. Using calculated electron densities, we employ Bader's quantum theory of atoms in molecules (QTAIM) to gain insights into the nature of metal-metal bonding in the clusters M(2), M(4), M(4)(2-), and M(4)Na(2) (M = Cu, Ag, Au). For the simplest case of the metal dimer, M(2), we correlate the strength of the metal-metal bond with the value of the electron density at the bond critical point, the total energy density at the bond critical point, the sharing (delocalization) index, and the values of the two principle negative curvatures. We then consider changes to the metal-metal bonding and charge density distribution upon the addition of two metal atoms to form the metal tetramer, M(4), and then followed by the addition of an electron pair to form M(4)(2-) and finally followed by the addition of two alkali metal (sodium) ions to form M(4)Na(2). Using topological properties of the electron density, we present evidence for the existence of σ-aromaticity in Au(4)(2-). We also report the existence of two non-nuclear attractors in the molecular graph of Cu(4)(2-) and large negative charge accumulation in the nonbonded Cu basins of this cluster.     

Study of ion specific interactions of alkali cations with dicarboxylate dianions

Alkali metal cations often show pronounced ion-specific interactions and selectivity with macromolecules in biological processes, colloids, and interfacial sciences, but a fundamental understanding about the underlying microscopic mechanism is still very limited. Here we report a direct probe of interactions between alkali metal cations (M(+)) and dicarboxylate dianions, (-)O(2)C(CH(2))(n)CO(2)(-) (D(n)(2-)) in the gas phase by combined photoelectron spectroscopy (PES) and ab initio electronic structure calculations on nine M(+)-D(n)(2-) complexes (M = Li, Na, K; n = 2, 4, 6). PES spectra show that the electron binding energy (EBE) decreases from Li(+) to Na(+) to K(+) for complexes of M(+)-D(2)(2-), whereas the order is Li(+) < Na(+) ≈ K(+) when M(+) interacts with a more flexible D(6)(2-) dianion. Theoretical modeling suggests that M(+) prefers to interact with both ends of the carboxylate -COO(-) groups by bending the flexible aliphatic backbone, and the local binding environments are found to depend upon backbone length n, carboxylate orientation, and the specific cation M(+). The observed variance of EBEs reflects how well each specific dicarboxylate dianion accommodates each M(+). This work demonstrates the delicate interplay among several factors (electrostatic interaction, size matching, and strain energy) that play critical roles in determining the structures and energetics of gaseous clusters as well as ion specificity and selectivity in solutions and biological systems.