Changes in the composition of the solvation shell of lithium ions in γ-butyrolactone upon addition of 15-crown-5 as studied by quantum chemistry

Quantum chemical modeling of Li⁺ ion transfer from the solvation shell of γ-butyrolactone (GBL) as the solvent to the cavity of 15-crown-5 (15C5) macrocyclic ligand was carried out. Calculations were performed using the PBE nonempirical density functional and an extended basis for the SBK pseudopotential. The solvation energy was included in the framework of the polarizable continuum model. The calculated geometric parameters of GBL and 15C5 molecules are in good agreement with experimental X-ray data. The energies and structures of the Li(GBL) _n ⁺ (n = 1–5) complexes and Li(GBL) _m (15C5)⁺ (m = 0–3) mixed complexes were calculated. The binding energy of the fifth GBL molecule is low; therefore, the Li⁺ ion is mainly surrounded by four GBL molecules. The formation of mixed complexes by consecutive displacement of GBL molecules from the solvation shell of the lithium ion leads to structures with the coordination number 5. The equilibrium constants of these processes were used to determine the dependence of the composition of the solvation complexes on the concentration of 15C5 in the system. The concentrations of the Li(15C5)⁺ and Li(GBL)(15C5)⁺ complexes appeared to be comparable. The revealed structural features of the Li⁺ solvation complexes in the GBL-15C5 system were used to analyze the operating efficiency of lithium power sources.     

A detailed investigation on the global minimum structures of mixed rare‐gas clusters: Geometry,energetics, and site occupancy

We performed a global minimum search of mixed rare‐gas clusters by applying an evolutionary algorithm (EA), which was recently proposed for binary atomic systems (Marques and Pereira, Chem. Phys. Lett. 2010, 485, 211). Before being applied to the potentials used in this work, the EA was further tested against results previously reported for the Ar_NXe_38?N clusters and several new putative global minima were discovered. We employed either simple Lennard‐Jones (LJ) potentials or more realistic functions to describe pair interactions in Ar_NKr_38?N, Ar_NXe_38?N, and Kr_NXe_38?N clusters. The long‐range tail of the pair‐potentials shows some influence on the energetic features and shape of the structure of clusters. In turn, core–shell type structures are mostly observed for global minima of the binary rare‐gas clusters, for both accurate and LJ potentials. However, the long‐range tail of the potential may have influence on the type of atoms that segregate on the surface or form the core of the cluster. While relevant differences for the preferential site occupancy occur between the two potentials for Ar_NKr_38?N (for N > 21), the type of atoms that segregate on the surface for Ar_NXe_38?N and Kr_NXe_38?N clusters is unaffected by the accuracy of the long‐range part of the interaction in almost all cases. Moreover, the global minimum search for model‐potentials in binary systems reveals that the surface‐site occupancy is mainly determined by the combination of two parameters: the size ratio of the two types of particles forming the cluster and the minimum‐energy ratio corresponding to the pair‐interactions between unlike atoms. © 2012 Wiley Periodicals, Inc.     

Theoretical prediction for the structures of gas phase lithium oxide clusters: (Li2O)n (n = 1–8)

We apply genetic algorithm combining directly with density functional method to search the potential energy surface of lithium‐oxide clusters (Li₂O)_n up to n = 8. In (Li₂O)_n (n = 1–8) clusters, the planar structures are found to be global minimum up to n = 2, and the global minimum structures are all three‐dimensional at n ≥ 3. At n ≥ 4, the tetrahedral unit (TU) is found in most of the stable structures. In the TU, the central Li is bonded with four O atoms in sp³ interactions, which leads to unusual charge transformation, and the probability of the central Li participating in the bonding is higher by adaptive natural density partitioning analysis, so the central Li is in particularly low positive charge. At large cluster size, distortion of structures is viewed, which breaks the symmetry and may make energy higher. The global minimum structures of (Li₂O)₂, (Li₂O)₆, and (Li₂O)₇ clusters are the most stable magic numbers, where the first one is planar and the later both have stable structural units of tetrahedral and C_4v. © 2012 Wiley Periodicals, Inc.     

A diatomics-in-molecules model for singly ionized neon clusters

A minimum-basis diatomics-in-molecules (DIM) model previously developed for singly-ionized argon clusters is applied to neon clusters, Ne _n ⁺ , forn=3, 4,...,22. A search for the global minimum energy of each cluster yields structures with the positive charge localised on a dimer-ion. This appears to be due largely to the linear unsymmetrical configuration which the model finds for Ne ₃ ⁺ . For this reason, the structures of the clusters at their minimum energy are different from those for Ar _n ⁺ computed with the same model. On the other hand, the behaviour of the charge distribution as a function of the geometrical configuration is similar to that for Ar _n ⁺ , as are the overall shapes of the potential energy surfaces. The results are discussed in terms of the charge distributions and the ratios of equilibrium properties of the dimers and dimer-ions which constitute the input to the model.     

Ion‐Pair SN2 Reaction of OH− and CH3Cl: Activation Strain Analyses of Counterion and Solvent Effects

We have theoretically studied the non‐identity S_N2 reactions of M_nOH^(n?1)+CH₃Cl (M⁺=Li⁺, Na⁺, K⁺, and MgCl⁺; n=0, 1) in the gas phase and in THF solution at the OLYP/6‐31++G(d,p) level using polarizable continuum model (PCM) implicit solvation. We want to explore and understand the effect of the metal counterion M⁺ and solvation on the reaction profile and the stereoselectivity of these processes. To this end, we have explored the potential energy surfaces of the backside (S_N2‐b) and frontside (S_N2‐f) pathways. To explain the computed trends, we have carried out analyses with an extended activation strain model (ASM) of chemical reactivity that includes the treatment of solvation effects.     

Electrostatic,sequential bond energies and structures of Li+·(N2)n complexes: computational study

The MP2 and CCSD calculations of the geometries and binding energies of the Li⁺·(N₂)_n (n?=?1–4) complexes are obtained. The potential energy surface showed that these complexes exhibit one minimum state and one transition state. The mono- and di-ligated complexes exhibit linear configurations with a binding energy of 11.1 and 21.2 kcal mol^?1, respectively. Trigonal planar and tetrahedral configurations are obtained for tri- and tetra-ligated complexes, respectively. The computed sequential bond dissociation energies (BDEs) of Li⁺·(N₂)_n (n?=?1–4) complexes are also calculated in which the mono-ligated complex has the largest BDE value. The obtained trend is mainly dependent on the variation in the ion-quadrupole interaction of these ion complexes. These calculations predict that these complexes are of purely electrostatic nature.

     

Comparison of line search minimization algorithms for exploring topography of multidimensional potential energy surfaces: Mg+Arn case

The most robust numerical algorithms for unconstrained optimization that involve a line search are tested in the problem of locating stable structures and transition states of atomic microclusters. Specifically, the popular quenching technique is compared with conjugate gradient and variable metric algorithms in the Mg⁺Ar_n clusters. It is found that the variable metric method BFGS combined with an approximate line minimization routine is the most efficient, and it shows global convergence properties. This technique is applied to find a few hundred stationary points of Mg⁺Ar₁₂ and to locate isomerization paths between the two most stable icosahedral structures found for Mg⁺Ar₁₂. The latter correspond to a solvated and a nonsolvated ion, respectively. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 :1011–1022, 1997     

Two-photon ionization of Li2: isotopic separation and determination of IP(Li2) and De(Li+2)

Complete isotope separation is achieved by two-photon ionization of Li₂ by a single mode Ar⁺ laser. With the use of two Ar⁺ lasers, the ionization potential of Li₂ is found to be 5.174 ± 0.013 eV, and the dissociation energy D_e(Li⁺₂) to be 1.274 ± 0.019 eV.     

FTIR and quantum chemical study of LiBr solvation in acetonitrile solutions

FTIR spectroscopy and quantum chemical calculations at the RTF + MP2/6-311G** level of theory with solvation model density (SMD) corrections were used to study ion solvation and association in LiBr/acetonitrile solutions. The aim of this study was to establish the composition and geometry of the predominant ionic species solvated by acetonitrile molecules and to analyse their spectroscopic signatures. The results obtained make it possible to propose an equilibrium between Li⁺Br⁻(CH₃CN)₃, Li⁺(CH₃CN)₄, and anionic Br⁻(CH₃CN)_n complexes with an undetermined n value and bent coordination of the solvent molecules. The calculated wavenumbers and the geometric parameters of the solvated ionic species were found to be in excellent agreement with the experimental data.     

Non-additivity in water-ion-water interactions

The three-body system Li⁺(H₂O)₂ was analyzed to study that non-additive part of the interaction potential which can be obtained by the Hartree-Fock approximation.For long and intermediate distances the three-body correction was found to be well represented by the induction energy, where bond dipoles are induced on each water molecule by point charges located on the (unpolarizable) lithium ion and on the other molecule respectively: for shorter distances this approximation was corrected by means of an exponential repulsive term. Such a potential model for non-additive interactions was extended to the more general situation Li⁺(H₂O)_n, and Monte-Carlo calculations were carried out on clusters containing up to six water molecules; comparison with other simulation results and with available data showed a significantly improved agreement with experiment. Tentative values for H are presented for n =7, 8,..., 20, where experimental data are not available.     

Ab initio potential energy surface for Ar 3 +

The potential energy surface (PES) of linear Ar ₃ ⁺ is calculated at the MP4/6-31G* level including all single, double, triple and quadruple excitations. The results show that the PES of the linear Ar ₃ ⁺ has a very flat valley along the asymmetric stretching vibration normal mode, ν₃. A higher level quadratic configuration interaction calculation including single, double and triple substitutions QCISD (T) along this flat valley suggests that an asymmetric geometry energy minimum reported earlier based on MP2 [1] is due to symmetry breaking in UHF. The global minimum of the PES is found to be for the symmetric geometry atR _ab =R _bc =2.66±0.01 Å, which is in good agreement with the MRD-CI calculation [2] and expectations from our earlier photodissociation experiments [3]. The calculational results are compared with other theoretical calculations, and are discussed in the context of the photodissociation and dynamics of dissociation experiments conducted on Ar ₃ ⁺ .     

Vibrational spectroscopic and density functional theory studies on ion solvation and ion association of lithium tetrafluoroborate in N,N-dimethylcarbamoyl chloride-based solvents

Solvation and association interactions in solutions of LiBF₄/DMCC (DMCC for N,N-dimethylcarbamoyl chloride) and LiBF₄/DMCC–DME (DME for 1,2-dimethoxyethane) have been studied as a function of concentration of lithium tetrafluoroborate by infrared and Raman spectroscopy. Strong interactions between Li⁺ and solvent molecules or BF₄⁻ anions are observed. The apparent solvation numbers of Li⁺ in LiBF₄/DMCC solutions were deduced. Band-fitting to the B–F stretching band of BF₄⁻ anion permits detailed assess of the ion pairing. Based on the calculations of density function theory, optimal structures of Li⁺(DMCC)_n (n = 1–3) were suggested. It is found that the lithium ion was preferentially solvated by DME in DMCC–DME binary solvents. This finding is supported by quantum chemistry calculations.     

Calculation of solvation free energies of Li+ and O2 ? ions and neutral lithium–oxygen compounds in acetonitrile using mixed cluster/continuum models

Solvation effects play a major role in determining the cycling characteristics of the non-aqueous rechargeable Li-air battery. We use a mixed cluster/continuum solvent model with varying number of explicit solvent molecules (n?=?4–10) to calculate the solvation free energies ( $ \Updelta G_{\text{solv}}^{*} $ ) of Li⁺ and O₂ ^? ions and neutral LiO₂, Li₂O₂, LiO, and Li₂O species in acetonitrile solvent. Calculations for complexes with the full first solvation shell around Li⁺ (n?=?4) and O₂ ^? (n?=?8) show excellent agreement with the solvation free energies obtained using the cluster pair approximation (the error is below 2.0?kcal/mol). The use of the pure continuum model fitted to reproduce the experimental values of $ \Updelta G_{\text{solv}}^{*} $ (Li⁺) and $ \Updelta G_{\text{solv}}^{*} $ (O₂ ^?) gives the solvation free energies of various lithium–oxygen species (Li _x O _y ; x, y?=?1, 2) that are in excellent agreement with the results obtained using mixed cluster/continuum models (n?≥?8). This provides a theoretical framework for including solvent effects in the theoretical models of oxygen reduction and evolution reactions in the aprotic Li-air battery.     

Multiphoton ionization and photodissociation of (NO)mArn clusters

Picosecond multiphoton ionization of (NO)_mAr_n clusters produced in a supersonic expansion of NO/Ar gas mixtures has been studied using time-of-flight mass spectrometry. Two-photon ionization with 266 nm photons show that dilute gas mixtures (1% NO/Ar) yield clusters limited to m≤7, but with as many as 37 argon atoms. Magic numbers are observed for NO⁺Ar₁₂, NO⁺Ar₁₈, (NO) ₂ ⁺ Ar₁₇, NO⁺Ar₂₂, and (NO) ₂ ⁺ Ar₂₁ and are understood in terms of solvation of the NO⁺ and (NO) ₂ ⁺ by argon in icosahedral arrangements. Four-photon ionization with 532 nm light produces dissociation of the clusters to yield only NO⁺Ar_n with n up to 54. This distribution exhibits an additional magic number at n=54, consistent with the completion of a second solvation sphere about the NO⁺. The known wavelength dependence for photodissociation of (NO) ₂ ⁺ and (NO) ₃ ⁺ and comparison of MPI spectra obtained with 266, 355, and 532 nm light indicate that the dissociation is occurring in the cluster ions.     

Formation of positive cluster ions LinBr (n = 2–7) and ionization energies studied by thermal ionization mass spectrometry

Clusters of the type Li_nX (X = halides) can be considered as potential building blocks of cluster‐assembly materials. In this work, Li_nBr (n = 2–7) clusters were obtained by a thermal ionization source of modified design and selected by a magnetic sector mass spectrometer. Positive ions of the Li_nBr (n = 4–7) cluster were detected for the first time. The order of ion intensities was Li₂Br⁺ > Li₄Br⁺ > Li₅Br⁺ > Li₆Br⁺ > Li₃Br⁺. The ionization energies (IEs) were measured and found to be 3.95 ± 0.20 eV for Li₂Br, 3.92 ± 0.20 eV for Li₃Br, 3.93 ± 0.20 eV for Li₄Br, 4.08 ± 0.20 eV for Li₅Br, 4.14 ± 0.20 eV for Li₆Br and 4.19 ± 0.20 eV for Li₇Br. All of these clusters have a much lower ionization potential than that of the lithium atom, so they belong to the superalkali class. The IEs of Li_nBr (n = 2–4) are slightly lower than those in the corresponding small Li_n or Li_nH clusters, whereas the IEs of Li_nBr are very similar to those of Li_n or Li_nH for n = 5 and 6. The thermal ionization source of modified design is an important means for simultaneously obtaining and measuring the IEs of Li_nBr (n = 2–7) clusters (because their ions are thermodynamically stable with respect to the loss of lithium atoms in the gas phase) and increasingly contributes toward the development of clusters for practical applications. Copyright © 2012 John Wiley & Sons, Ltd.     

Correlating the Solvating Power of Solvents with the Strength of Ion-Dipole Interaction in Electrolytes of Lithium-ion Batteries

The solvation structure of Li⁺ plays a significant role in determining the physicochemical properties of electrolytes. However, to date, there is still no clear definition of the solvating power of different electrolyte solvents, and even the solvents that preferentially participate in the solvation structure remain controversial. In this study, we comprehensively discuss the solvating power and solvation process of Li⁺ ions using both experimental characterizations and theoretical calculations. Our findings reveal that the solvating power is dependent on the strength of the Li⁺-solvent (ion-dipole) interaction. Additionally, we uncover that the anions tend to enter the solvation sheath in most electrolyte systems through Li⁺-anion (ion-ion) interaction, which is weakened by the shielding effect of solvents. The competition between the Li⁺-solvent and Li⁺-anion interactions ultimately determines the final solvation structures. This insight into the fundamental understanding of the solvation structure of Li⁺ provides inspiration for the design of multifunctional mixed-solvent electrolytes for advanced batteries.     

Ab initio pair potentials for the ionic lithium-formate system

The potential energy surfaces for the interatomic interaction in the Li⁺HCOO^– system have been investigated byab initio methods within the rigid-molecule approximation. Analytical potential expressions were fitted to 133 calculated SCF energies for the Li⁺-HCOO^– interaction, 42 SCF energies for the Li⁺-Li⁺ interaction, and 332 SCF energies for the HCOO^–-HCOO^– interaction. The global minimum on the Li⁺-HCOO^– surface is –170 kcal/mol and corresponds to the lithium ion lying on the C₂ axis of the formate ion at 2.2 Å from the carbon atom on the oxygen side. The cation-cation and anion-anion interactions are repulsive everywhere, although the potential surface is markedly anisotropic for the HCOO^–-HCOO^– interaction.     

Ab initio copper–water interaction potential for the simulation of aqueous solutions

A new ab initio effective two-body potential that aims at mimicking the average copper–water interaction energy of the first solvation shell was developed. This new potential, together with the MCY water–water potential and a three-body ion–water–water induction potential, is tested in simulations of gas-phase clusters [Cu²⁺? (H₂O)₂₀] and diluted solutions [Cu²⁺? (H₂O)₂₀₀] at T = 298 K. The results of simulations with conventional ab initio pair potentials, with and without three-body induction corrections, are also presented. The different types of copper–water interaction potentials are evaluated comparatively and the efficiency of the newly proposed effective pair potential is discussed. © 1993 John Wiley & Sons, Inc.     

Effect of Cation Size on Solvation and Association with Superoxide Anion in Aprotic Solvents

Influence of cation size on solvation strength, diffusion, and kinetics of the association reaction with anions O₂⁻ in aprotic solvents, such as acetonitrile and dimethyl sulfoxide, has been investigated by means of molecular dynamics simulations. The work is motivated by the need to understand the molecular nature of the solvent-induced changes in capacity of Li-air batteries. We have shown that the dependence of the solvation shell stability on the cation size has a maximum at a particular ion radius that corresponds to a solvent coordination number of 4. The shell stability maximum coincides with the diffusion coefficient minimum. The variation of the cation shell stability has a crucial impact on the kinetics of the cation-O₂⁻ association. We have demonstrated that profound inhibition of the association reaction for Li⁺ in dimethyl sulfoxide is a result of the lock-and-key effect that cannot be described in the framework of Hard Soft Acid Base theory.     

One-Electron Pseudo-Potential Determination of Stable Isomers of the Li*Ar n Excited Clusters: Absorption Spectra

We present pseudo-potential calculations of geometrical structures of stable isomers of LiAr _n clusters with both an electronic ground state and excited states of the lithium atom. The Li atom is perturbed by argon atoms in LiAr _n clusters. Its electronic structure obtained as the eigenfunctions of a single-electron operator describing the electron in the field of a Li⁺Ar _n core, the Li⁺ and Ar atoms are replaced by pseudo-potentials. These pseudo-potentials include core-polarization operators to account for the polarization and correlation of the inert core with the valence Lithium electron [J Chem Phys 116, 1839 1]. The geometry optimization of the ground and excited states of LiAr _n (n = 1–12) clusters is carried out via the Basin-Hopping method of Wales et al. [J Phys Chem 101, 5111 2; J Chem Phys 285, 1368 3]. The geometries of the ground and ionic states of LiAr _n clusters were used to determine the energy of the high excited states of the neutral LiAr _n clusters. The variation of the excited state energies of LiAr _n clusters as a function of the number of argon atoms shows an approximate Rydberg character, corresponding to the picture of an excited electron surrounding an ionic cluster core, is already reached for the 3s state. The result of optical transitions calculations shows that the absorption spectral features are sensitive to isomer structure. It is clearly the case for transitions close to the 2p levels of Li which are distorted by the cluster environment.