The solvation of Li+ and Na+ in acetonitrile from ab initio-derived many-body ion–solvent potentials

Several Li+- and Na+-acetonitrile models were derived from ab initio calculations at the counterpoise-corrected MP2/TZV++(d,p) level for distorted ion-(MeCN)n clusters with n=1, 4 and 6. Two different many-body ion-acetonitrile models were constructed: an effective three-body potential for use with the six-site effective pair model of Bohm et al., and an effective polarizable many-body model. The polarizable acetonitrile model used in the latter model is a new empirical model which was also derived in the present paper. Mainly for comparative purposes, two ion-acetonitrile pair potentials were also constructed from the ab initio cluster calculations: one pure pair potential and one effective pair potential. Using all these potential models, MD simulations in the NPT ensemble were performed for the pure acetonitrile liquid and for Li+(MeCN) and Na+(MeCN) solutions with 1 ion in 512 solvent molecules and with a simulation time of at least 120 ps per system. Thermodynamic properties, solvation-shell structure and the self-diffusion coefficient of the ions and of the solvent molecules were calculated and compared between the different models and with experimental data, where available. The Li+ ion is found to be four-coordinated when the new many-body potentials are used, in contrast to the six-coordinated structure obtained for the pure pair and effective pair potentials. The coordination number of Na+ is close to six for all the models derived here, although the coordination number becomes slightly smaller with the many-body potentials. For both ions, the solvent molecules in the first shell point their nitrogen ends towards the cation, while in the second shell the opposite orientation is the most common.

[1]  P. Kusalik,et al.  Structure in liquid water: A study of spatial distribution functions , 1993 .

[2]  M. Månsson,et al.  Enthalpies of combustion and formation of acetonitrile , 1983 .

[3]  P. Jedlovszky,et al.  Reverse Monte Carlo simulation of a heteronuclear molecular liquid: Structural study of acetonitrile , 1994 .

[4]  Ian R. McDonald,et al.  An effective pair potential for liquid acetonitrile , 1983 .

[5]  Mark S. Gordon,et al.  Self‐consistent molecular orbital methods. XXIII. A polarization‐type basis set for second‐row elements , 1982 .

[6]  Eric Barth,et al.  Algorithms for constrained molecular dynamics , 1995, J. Comput. Chem..

[7]  E. Hawlicka Acetonitrile-Water Solutions of Sodium Halides: Viscosity and Self-Diffusion of CH3CN and H2O , 1988 .

[8]  E. Guàrdia,et al.  On the solvation shell of Na+ and Cl- ions in acetonitrile: a computer simulation study , 2000 .

[9]  R. Impey,et al.  Ionic solvation in nonaqueous solvents: the structure of lithium ion and chloride in methanol, ammonia, and methylamine , 1987 .

[10]  L. Degrève,et al.  From ionic aqueous solvation shell to bulk fluid: A structural‐energetic stability problem , 1994 .

[11]  P. C. Hariharan,et al.  The influence of polarization functions on molecular orbital hydrogenation energies , 1973 .

[12]  L. Pejov Interaction of acetonitrile with alkaline metal cations: A density functional, coupled‐cluster, and quadratic configuration interaction study , 2002 .

[13]  A. D. McLean,et al.  Contracted Gaussian basis sets for molecular calculations. I. Second row atoms, Z=11–18 , 1980 .

[14]  Michael J. Frisch,et al.  Self‐consistent molecular orbital methods 25. Supplementary functions for Gaussian basis sets , 1984 .

[15]  Harold Basch,et al.  Compact effective potentials and efficient shared‐exponent basis sets for the first‐ and second‐row atoms , 1984 .

[16]  Kurt Kremer,et al.  Molecular dynamics simulation of a polymer chain in solution , 1993 .

[17]  M. Dixon,et al.  Interionic potentials in alkali halides and their use in simulations of the molten salts , 1976 .

[18]  J. R. Carl,et al.  Atom dipole interaction model for molecular polarizability. Application to polyatomic molecules and determination of atom polarizabilities , 1972 .

[19]  S. Nosé A molecular dynamics method for simulations in the canonical ensemble , 1984 .

[20]  Steven J. Stuart,et al.  Potentials and Algorithms for Incorporating Polarizability in Computer Simulations , 2003 .

[21]  E. Cabaleiro-Lago,et al.  Ab initio study of M(CH3CN)n clusters (M=Li+, Na+, Mg2+) in the gas phase , 2000 .

[22]  W. Li,et al.  Determination of the transfer thermodynamic functions for some monovalent ions from water to N,N-dimethylthioformamide, and for some anions from water to methanol, dimethyl sulfoxide, acetonitrile and pyridine, and standard electrode potentials of some M+/M(s) couples in N,N-dimethylthioformamide , 1994 .

[23]  K. Hermansson,et al.  Effective three-body potentials for Li+(aq) and Mg2+(aq) , 2003 .

[24]  E. Hawlicka,et al.  MD Simulation Studies of Selective Solvation in Methanol−Water Mixtures: An Effect of the Charge Density of a Solute , 2002 .

[25]  P. Kollman,et al.  Structure and Properties of Neat Liquids Using Nonadditive Molecular Dynamics: Water, Methanol, and N-Methylacetamide , 1995 .

[26]  Mark S. Gordon,et al.  General atomic and molecular electronic structure system , 1993, J. Comput. Chem..

[27]  L. Troxler,et al.  Conformation and Dynamics of 18-Crown-6, Cryptand 222, and Their Cation Complexes in Acetonitrile Studied by Molecular Dynamics Simulations , 1994 .

[28]  Gerhard Hummer,et al.  Free Energy of Ionic Hydration , 1996 .

[29]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[30]  J. Pople,et al.  Self‐Consistent Molecular‐Orbital Methods. IX. An Extended Gaussian‐Type Basis for Molecular‐Orbital Studies of Organic Molecules , 1971 .

[31]  W. C. Swope,et al.  A computer simulation method for the calculation of equilibrium constants for the formation of physi , 1981 .

[32]  A. Sum,et al.  Prediction of the phase behavior of acetonitrile and methanol with ab initio pair potentials. I. Pure components , 2002 .

[33]  Mark A. Fryling,et al.  Investigation of lithium−water interactions in acetonitrile solutions using proton nuclear magnetic resonance, Raman, and infrared spectroscopies and extended Hückel molecular orbital calculations , 1989 .

[34]  V. P. Sazonov,et al.  IUPAC-NIST Solubility Data Series. 78. Acetonitrile Binary Systems , 2002 .

[35]  R. Fuoss,et al.  Electrolyte-solvent interaction. XX. Picric acid in mixtures of acetonitrile and basic solvents , 1969 .

[36]  D Fincham,et al.  Shell model simulations by adiabatic dynamics , 1993 .

[37]  Michel Dupuis,et al.  The general atomic and molecular electronic structure system hondo: Version 7.0 , 1989 .

[38]  J. Barthel,et al.  Conductance of 1,1-electrolytes in acetonitrile solutions from −40° to 35°C , 1990 .

[39]  John Burgess,et al.  Metal Ions in Solution , 1978 .

[40]  K. Szalewicz,et al.  Comment on “On the importance of the fragment relaxation energy terms in the estimation of the basis set superposition error correction to the intermolecular interaction energy” [J. Chem. Phys. 104, 8821 (1996)] , 1998 .

[41]  William L. Jorgensen,et al.  Monte Carlo simulations of liquid acetonitrile with a three-site model , 1988 .

[42]  E. Kamieńska-Piotrowicz,et al.  Determination of the transfer thermodynamic functions for the zinc(II), cadmium(II), mercury(II) and mercury(I) ions from water to methanol, dimethyl sulfoxide, acetonitrile, pyridine and N,N-dimethylthioformamide and of standard electrode potentials of the M2+/M(s) couples in these solvents , 1994 .

[43]  P. Kebarle,et al.  Ionic solvation by aprotic solvents. Gas phase solvation of the alkali ions by acetonitrile , 1976 .

[44]  J. Barthel,et al.  FTIR study of ion solvation and ion-pair formation in alkaline and alkaline earth metal salt solutions in acetonitrile , 1994 .

[45]  J. Reimers,et al.  The Solvation of Acetonitrile , 1999 .

[46]  H. C. Andersen Rattle: A “velocity” version of the shake algorithm for molecular dynamics calculations , 1983 .

[47]  E. Guàrdia,et al.  Molecular dynamics study of Na+ and Cl− in methanol , 1996 .

[48]  F. Hörzenberger,et al.  Gibbs energies, enthalpies and entropies of transfer of cations from acetonitrile and N,N-dimethylformamide into water , 1992 .

[49]  Y. Marcus Thermodynamic functions of transfer of single ions from water to nonaqueous and mixed solvents: Part I - Gibbs free energies of transfer to nonaqueous solvents , 1983 .

[50]  W. Gordy,et al.  Precision measurement of dipole moments and other spectral constants of normal and deuterated methyl fluoride and methyl cyanide , 1966 .

[51]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[52]  Edward L Cussler,et al.  Diffusion: Mass Transfer in Fluid Systems , 1984 .

[53]  A. W. Overhauser,et al.  Theory of the Dielectric Constants of Alkali Halide Crystals , 1958 .

[54]  T. Radnai,et al.  A mass spectrometric study of solvated clusters of ions and ion pairs generated from lithium halide solutions in polar solvents: Acetonitrile compared to methanol , 2003 .

[55]  J. Beauchamp,et al.  Intrinsic acid-base properties of molecules. Binding energies of lithium(1+) ion to .pi.- and n-donor bases , 1975 .

[56]  Surjit Singh,et al.  Infrared spectral studies on preferential solvation of lithium ions in binary mixtures of acetonitrile and dimethyl formamide , 1996 .

[57]  W. L. Jorgensen,et al.  Quantum and statistical mechanical studies of liquids. 21. The nature of dilute solutions of sodium and methoxide ions in methanol , 1982 .

[58]  H. Krienke,et al.  The solvation of ions in acetonitrile and acetone. II. Monte Carlo simulations using polarizable solvent models , 2002 .

[59]  F. Millero Molal volumes of electrolytes , 1971 .

[60]  Masaaki Tabata,et al.  Liquid Structure of Acetonitrile−Water Mixtures by X-ray Diffraction and Infrared Spectroscopy , 1998 .

[61]  W. Kunz,et al.  Lithium bromide in acetonitrile and water: a neutron scattering study , 1991 .

[62]  E. Hawlicka Self-Diffusion of Sodium, Chloride and Iodide Ions in Acetonitrile-Water Mixtures , 1987 .

[63]  J. G. Snijders,et al.  Solving the Dirac equation, using the large component only, in a Dirac-type Slater orbital basis set , 1995 .

[64]  Ingemar Wadsö,et al.  Enthalpies of Vaporization of Organic Compounds. IV. Alkyl Nitriles. , 1970 .

[65]  H. Krienke,et al.  The solvation of ions in acetonitrile and acetone: A molecular Ornstein–Zernike study , 1998 .

[66]  W. Fawcett,et al.  Attenuated total reflection fourier-transform infrared spectroscopic study of ion–solvent and ion–ion interactions in alkali-metal perchlorate–acetonitrile solutions , 1993 .

[67]  W. L. Jorgensen,et al.  The nature of dilute solutions of sodium ion in water, methanol, and tetrahydrofuran , 1982 .