TDMP2 calculation of dynamic multipole polarizabilities and dispersion coefficients for the halogen anions F−, Cl−, Br− and I−

TURBOMOLE: Today and Tomorrow

TURBOMOLE is a highly optimized software suite for large-scale quantum-chemical and materials science simulations of molecules, clusters, extended systems, and periodic solids. TURBOMOLE uses Gaussian basis sets and has been designed with robust and fast quantum-chemical applications in mind, ranging from homogeneous and heterogeneous catalysis to inorganic and organic chemistry and various types of spectroscopy, light-matter interactions, and biochemistry. This Perspective briefly surveys TURBOMOLE's functionality and highlights recent developments that have taken place between 2020 and 2023, comprising new electronic structure methods for molecules and solids, previously unavailable molecular properties, embedding, and molecular dynamics approaches. Select features under development are reviewed to illustrate the continuous growth of the program suite, including nuclear electronic orbital methods, Hartree-Fock-based adiabatic connection models, simplified time-dependent density functional theory, relativistic effects and magnetic properties, and multiscale modeling of optical properties.

Pressure in Molecular Simulations with Scaled Charges. 1. Ionic Systems

Charge scaling, rationalized as MDEC (molecular dynamics in electronic continuum) or ECC (electronic continuum correction), has become a widely used simple approach to how to avoid self-consistent induced dipoles yet approximately take into account the effects of electronic polarizability. It has been assumed that the continuum permittivity does not depend on density; in turn, pressure is calculated by standard formulas. In this work, we elaborate a complementary approximation of density-independent molecular polarizability and derive formulas for pressure corrections within the MDEC framework; real behavior lies between these two extremes. The pressure corrections for test ionic systems are huge and negative, leading to sizable densities in constant-pressure MDEC simulations. A comparison of MDEC results with equivalent polarizable systems gives a good pressure match for a crystal but very low MDEC pressures for ionic liquids. These results witness about the importance of a correct density dependence not only of continuum permittivity in MDEC simulations but also of polarizability in polarizable simulations.

The first example of fluorescence of the solid individual compounds of Eu2+ ion: EuCl2 , EuI2 , EuBr2

Fluorescence (FL) from Eu²⁺ ion solid individual compounds was observed for the first time using the examples EuCl₂ , EuI₂ and EuBr₂ . Fluorescence was detected at excitation by ultraviolet (UV) light from powdered samples of Eu (II) halides at room temperature (RT) in an argon atmosphere. In air, FL of all Eu²⁺ compounds studied was stable, and intensity persisted for weeks. Depending on the nature of the halide anion and due to a nephelauxetic effect, the position of the maxima in the FL spectra underwent a red shift in the series: EuCl₂ (409 nm) < EuBr₂ (428 mm) < EuI₂ (432 nm). The lifetimes of the excited states of the Eu²⁺ * ions for EuI₂ , EuCl₂ and EuBr₂ were 355, 76 and 54 ns, respectively.

Developing multisite empirical force field models for Pt(II) and cisplatin

We have developed empirical force field parameters for Pt(II) and cisplatin. Two force field frameworks were used-modified OPLS-AA and our second-order polarizable POSSIM. A seven-site model was used for the Pt(II) ion. The goal was to create transferable parameter sets compatible with the force field models for proteins and general organic compounds. A number of properties of the Pt(II) ion and its coordination compounds have been considered, including geometries and energies of the complexes, hydration free energy, and radial distribution functions in water. Comparison has been made with experimental and quantum mechanical results. We have demonstrated that both versions are generally capable of reproducing key properties of the system, but the second-order polarizable POSSIM formalism permits more accurate quantitative results to be obtained. For example, the energy of formation of cisplatin as calculated with the modified OPLS-AA exhibited an error of 9.9%, while the POSSIM error for the same quantity was 6.2%. The produced parameter sets are transferable and suitable to be used in protein-metal binding simulations in which position or even coordination of the ion does not have to be constrained using preexisting knowledge. © 2016 Wiley Periodicals, Inc.

Structural, dynamical, and transport properties of the hydrated halides: How do At(-) bulk properties compare with those of the other halides, from F(-) to I(-)?

The properties of halides from the lightest, fluoride (F(-)), to the heaviest, astatide (At(-)), have been studied in water using a polarizable force-field approach based on molecular dynamics (MD) simulations at the 10 ns scale. The selected force-field explicitly treats the cooperativity within the halide-water hydrogen bond networks. The force-field parameters have been adjusted to ab initio data on anion/water clusters computed at the relativistic Möller-Plesset second-order perturbation theory level of theory. The anion static polarizabilities of the two heaviest halides, I(-) and At(-), were computed in the gas phase using large and diffuse atomic basis sets, and taking into account both electron correlation and spin-orbit coupling within a four-component framework. Our MD simulation results show the solvation properties of I(-) and At(-) in aqueous phase to be very close. For instance, their first hydration shells are structured and encompass 9.2 and 9.1 water molecules at about 3.70 ± 0.05 Å, respectively. These values have to be compared to the F(-), Cl(-), and Br(-) ones, i.e., 6.3, 8.4, and 9.0 water molecules at 2.74, 3.38, and 3.55 Å, respectively. Moreover our computations predict the solvation free energy of At(-) in liquid water at ambient conditions to be 68 kcal mol(-1), a value also close the I(-) one, about 70 kcal mol(-1). In all, our simulation results for I(-) are in excellent agreement with the latest neutron- and X-ray diffraction studies. Those for the At(-) ion are predictive, as no theoretical or experimental data are available to date.

Nuclear quantum effects in water exchange around lithium and fluoride ions

We employ classical and ring polymer molecular dynamics simulations to study the effect of nuclear quantum fluctuations on the structure and the water exchange dynamics of aqueous solutions of lithium and fluoride ions. While we obtain reasonably good agreement with experimental data for solutions of lithium by augmenting the Coulombic interactions between the ion and the water molecules with a standard Lennard-Jones ion-oxygen potential, the same is not true for solutions of fluoride, for which we find that a potential with a softer repulsive wall gives much better agreement. A small degree of destabilization of the first hydration shell is found in quantum simulations of both ions when compared with classical simulations, with the shell becoming less sharply defined and the mean residence time of the water molecules in the shell decreasing. In line with these modest differences, we find that the mechanisms of the exchange processes are unaffected by quantization, so a classical description of these reactions gives qualitatively correct and quantitatively reasonable results. We also find that the quantum effects in solutions of lithium are larger than in solutions of fluoride. This is partly due to the stronger interaction of lithium with water molecules, partly due to the lighter mass of lithium and partly due to competing quantum effects in the hydration of fluoride, which are absent in the hydration of lithium.

Extending the applicability of the Tkatchenko-Scheffler dispersion correction via iterative Hirshfeld partitioning

Recently we have demonstrated that the applicability of the Tkatchenko-Scheffler (TS) method for calculating dispersion corrections to density-functional theory can be extended to ionic systems if the Hirshfeld method for estimating effective volumes and charges of atoms in molecules or solids (AIM's) is replaced by its iterative variant [T. Bučko, S. Lebègue, J. Hafner, and J. Ángyán, J. Chem. Theory Comput. 9, 4293 (2013)]. The standard Hirshfeld method uses neutral atoms as a reference, whereas in the iterative Hirshfeld (HI) scheme the fractionally charged atomic reference states are determined self-consistently. We show that the HI method predicts more realistic AIM charges and that the TS/HI approach leads to polarizabilities and C6 dispersion coefficients in ionic or partially ionic systems which are, as expected, larger for anions than for cations (in contrast to the conventional TS method). For crystalline materials, the new algorithm predicts polarizabilities per unit cell in better agreement with the values derived from the Clausius-Mosotti equation. The applicability of the TS/HI method has been tested for a wide variety of molecular and solid-state systems. It is demonstrated that for systems dominated by covalent interactions and/or dispersion forces the TS/HI method leads to the same results as the conventional TS approach. The difference between the TS/HI and TS approaches increases with increasing ionicity. A detailed comparison is presented for isoelectronic series of octet compounds, layered crystals, complex intermetallic compounds, and hydrides, and for crystals built of molecules or containing molecular anions. It is demonstrated that only the TS/HI method leads to accurate results for systems where both electrostatic and dispersion interactions are important, as illustrated for Li-intercalated graphite and for molecular adsorption on the surfaces in ionic solids and in the cavities of zeolites.

Role of anion polarizability in fluorescence sensitization of DNA-templated silver nanoclusters

Fluorescent silver nanoclusters (Ag NCs) as novel fluorophores have received much attention because of their high brightness, good photostability and widely tunable emissions from the visible to the near-infrared range as a result of their size and existing environment. However, efforts are still needed to find the factors that tune the emission of Ag NCs. In this work, Ag NCs that were size-selectively grown on DNA were used to investigate the effect of the electronic properties of coordinating ligands. Halogen anions were used as the paradigm because of their periodicity in element properties. We found that addition of halogen anions did not alter the emission wavelength of Ag NCs, but the fluorescence intensity showed an initial increase at low concentrations of Cl(-), Br(-) and I(-) followed by a gradual decrease at high concentrations. No increase in fluorescence was observed for F(-) at either low or high concentration. Such specific halogen-anion sensitization of the fluorescence of Ag NCs suggests that the binding strength/manner and dipole polarizability of these anions synergistically tune the emission behavior of Ag NCs. Less fluorescence sensitization occurred for the anion having high enough polarizability to form a covalent bond with Ag NCs. The anion polarizability-sensitized fluorescence indicates the role of anion electronic properties in tuning the emission behavior of Ag NCs, which should be seriously considered in designing Ag NC-based sensors and devices.

A continuum solvent model of the multipolar dispersion solvation energy

The dispersion energy is an important contribution to the total solvation energies of ions and neutral molecules. Here, we present a new continuum model calculation of these energies, based on macroscopic quantum electrodynamics. The model uses the frequency dependent multipole polarizabilities of molecules in order to accurately calculate the dispersion interaction of a solute particle with surrounding water molecules. It includes the dipole, quadrupole, and octupole moment contributions. The water is modeled via a bulk dielectric susceptibility with a spherical cavity occupied by the solute. The model invokes damping functions to account for solute-solvent wave function overlap. The assumptions made are very similar to those used in the Born model. This provides consistency and additivity of electrostatic and dispersion (quantum mechanical) interactions. The energy increases in magnitude with cation size, but decreases slightly with size for the highly polarizable anions. The higher order multipole moments are essential, making up more than 50% of the dispersion solvation energy of the fluoride ion. This method provides an accurate and simple way of calculating the notoriously problematic dispersion contribution to the solvation energy. The result establishes the importance of using accurate calculations of the dispersion energy for the modeling of solvation.

Polarizable Poisson-Boltzmann equation: the study of polarizability effects on the structure of a double layer

We incorporate ion polarizabilities into the Poisson-Boltzmann equation by modifying the effective dielectric constant and the Boltzmann distribution of ions. The extent of the polarizability effects is controlled by two parameters, γ(1) and γ(2); γ(1) determines the polarization effects in a dilute system and γ(2) regulates the dependence of the polarizability effects on the concentration of ions. For a polarizable ion in an aqueous solution γ(1) ≈ 0.01 and the polarizability effects are negligible. The conditions where γ(1) and/or γ(2) are large and the polarizability is relevant involve the low dielectric constant media, high surface charge, and/or large ionic concentrations.

Ions in solutions: Determining their polarizabilities from first-principles

Dipole polarizabilities of a series of ions in aqueous solutions are computed from first-principles. The procedure is based on the study of the linear response of the maximally localized Wannier functions to an applied external field, within density functional theory. For most monoatomic cations (Li(+), Na(+), K(+), Rb(+), Mg(2+), Ca(2+) and Sr(2+)) the computed polarizabilities are the same as in the gas phase. For Cs(+) and a series of anions (F(-), Cl(-), Br(-) and I(-)), environmental effects are observed, which reduce the polarizabilities in aqueous solutions with respect to their gas phase values. The polarizabilities of H((aq)) (+), OH((aq)) (-) have also been determined along an ab initio molecular dynamics simulation. We observe that the polarizability of a molecule instantaneously switches upon proton transfer events. Finally, we also computed the polarizability tensor in the case of a strongly anisotropic molecular ion, UO(2) (2+). The results of these calculations will be useful in building interaction potentials that include polarization effects.

The polarizable point dipoles method with electrostatic damping: implementation on a model system

Recently, the use of polarizable force fields in Molecular Dynamics simulations has been gaining importance, since they allow a better description of heterogeneous systems compared to simple point charges force fields. Among the various techniques developed in the last years the one based on polarizable point dipoles represents one of the most used. In this paper, we review the basic technical issues of the method, illustrating the way to implement intramolecular and intermolecular damping of the electrostatic interactions, either with and without the Ewald summation method. We also show how to reduce the computational overhead for evaluating the dipoles, introducing to the state-of-the-art methods: the extended Lagrangian method and the always stable predictor corrector method. Finally we discuss the importance of screening the electrostatic interactions at short range, defending this technique against simpler approximations usually made. We compare results of density functional theory and classical force field-based Molecular Dynamics simulations of chloride in water.

Radii of atomic ions determined from diatomic ion-He bond lengths

We propose a new definition of the effective radius of an atomic ion: the bond distance (R(e)) of the ion/He diatomic complex minus the van der Waals radius of the helium atom. Our rationale is that He is the most chemically inert and least polarizable atom, so that its interaction with the outer portions of the electron cloud causes the smallest perturbation of it. We show that such radii, which we denote R(XHe), make good qualitative sense. We also compare our R(XHe) values to more traditional ionic radii from solid crystal X-ray measurements, as well as estimates of such radii from "ionic" gas-phase MF, MOM, MF(+), and MO molecules, where M is a metal atom. Such comparisons lead to interesting conclusions about bonding in ionic crystals and in simple gas-phase oxide and fluoride molecules. The definition is shown to be reasonable for -1, +1, and even for many of the larger +2 atomic ions. Another advantage of the R(XHe) definition is that it is also consistently valid for ground states and excited states of both neutral atoms and atomic ions, even for open-shell np and nd cases where the electron clouds of the ions are not spherically symmetric and R(XHe) thus depends on the "approach" direction of the He atom. Finally, we note that when there is a contribution from covalent bonding with the He atom, and/or in cases where the ion is small and has a very high charge, so that there is distortion even of the He 1s electrons, R(XHe) is not expected to be representative of the size of the ion. We then suggest that in these cases small, and sometimes unphysical, values of R(XHe) are diagnostic of the fact that simple "physical" interactions have been supplemented by a "chemical" component.

Simulating Monovalent and Divalent Ions in Aqueous Solution Using a Drude Polarizable Force Field

An accurate representation of ion solvation in aqueous solution is critical for meaningful computer simulations of a broad range of physical and biological processes. Polarizable models based on classical Drude oscillators are introduced and parametrized for a large set of monoatomic ions including cations of the alkali metals (Li(+), Na(+), K(+), Rb(+) and Cs(+)) and alkaline earth elements (Mg(2+), Ca(2+), Sr(2+) and Ba(2+)) along with Zn(2+) and halide anions (F(-), Cl(-), Br(-) and I(-)). The models are parameterized, in conjunction with the polarizable SWM4-NDP water model [Lamoureux et al., Chem. Phys. Lett. 418, 245 (2006)], to be consistent with a wide assortment of experimentally measured aqueous bulk thermodynamic properties and the energetics of small ion-water clusters. Structural and dynamic properties of the resulting ion models in aqueous solutions at infinite dilution are presented.

Absolute hydration free energy scale for alkali and halide ions established from simulations with a polarizable force field

A polarizable potential function for the hydration of alkali and halide ions is developed on the basis of the recent SWM4-DP water model [Lamoureux, G.; MacKerell, A. D., Jr.; Roux, B. J. Chem. Phys. 2003, 119, 5185]. Induced polarization is incorporated using classical Drude oscillators that are treated as auxiliary dynamical degrees of freedom. The ions are represented as polarizable Lennard-Jones centers, whose parameters are optimized to reproduce the binding energies of gas-phase monohydrates and the hydration free energies in the bulk liquid. Systematic exploration of the parameters shows that the monohydrate binding energies can be consistent with a unique hydration free energy scale if the computed hydration free energies incorporate the contribution from the air/water interfacial electrostatic potential (-540 mV for SWM4-DP). The final model, which can satisfyingly reproduce both gas and bulk-phase properties, corresponds to an absolute scale in which the intrinsic hydration free energy of the proton is -247 kcal/mol.

Excess compressibility in binary liquid mixtures

Brillouin scattering experiments have been carried out on some mixtures of molecular liquids. From the measurement of the hypersonic velocities we have evaluated the adiabatic compressibility as a function of the volume fraction. We show how the quadratic form of the excess compressibility dependence on the solute volume fraction can be derived by simple statistical effects and does not imply any interaction among the components of the system other than excluded volume effects. This idea is supported by the comparison of the experimental results with a well-established prototype model, consisting of a binary mixture of hard spheres with a nonadditive interaction potential. This naive model turns out to be able to produce a very wide spectrum of structural and thermodynamic features depending on values of its parameters. An attempt has made to understand what kind of structural information can be gained through the analysis of the volume fraction dependence of the compressibility.

Interaction potentials for Br−–Rg (Rg=He–Rn): Spectroscopy and transport coefficients

High-level ab initio CCSD(T) calculations are performed in order to obtain accurate interaction potentials for the Br(-) anion interacting with each rare gas (Rg) atom. For the Rg atoms from He to Ar, two approaches are taken. The first one implements a relativistic core potential and an aug-cc-pVQZ basis set for bromine, an aug-cc-pV5Z basis set for Rg, and a set of bond functions placed at the midpoint of the Rg-Br distance. The second one uses the all-electron approximation with aug-cc-pV5Z bases further augmented by an extra diffuse function in each shell. Comparison reveals close similarity between both sets of results, so for Rg atoms from Kr to Rn only the second approach is exploited. Calculated potentials are assessed against the previous empirical, semiempirical, and ab initio potentials, and against available beam scattering data, zero electron kinetic energy spectroscopic data, and various sets of the measured ion mobilities and diffusion coefficients. This multiproperty analysis leads to the conclusion that the present potentials are consistently good for the whole series of Br(-)-Rg pairs over the whole range of internuclear distances covered.

Solvation of N3- at the Water Surface: The Polarizable Continuum Model Approach

We present a new quantum mechanical model to introduce Pauli repulsion interaction between a molecular solute and the surrounding solvent in the framework of the Polarizable Continuum Model. The new expression is derived in a way to allow naturally for a position-dependent solvent density. This development makes it possible to employ the derived expression for the calculation of molecular properties at the interface between two different dielectrics. The new formulation has been tested on the azide anion (N3-) for which we have calculated the solvation energy, the dipole moment, and the static polarizability at the interface as a function of the ion position. The calculations have been carried out for different ion-surface orientations, and the results have also been compared with the parallel electrostatic-only solvation model.

Interaction potentials of the RG-I anions, neutrals, and cations (RG = He, Ne, Ar)

Interaction potentials of the iodine atom, atomic cation, and anion with light rare-gas atoms from He to Ar are calculated within the unified ab initio approach using the unrestricted coupled-cluster with singles and doubles and perturbative treatment of triples correlation treatment, relativistic small-core pseudopotential, and an extended basis set. Ab initio points are fit to a flexible analytical function. The calculated potentials are compared with available literature data, assessed in the I(-)-and I+-ion mobility calculations and the Ar-I(-)-anion zero electron kinetic-energy spectra simulations, and analyzed using the correlation rules. The results indicate a high precision of the reported potentials.

Ion Solvation Thermodynamics from Simulation with a Polarizable Force Field

Thermodynamic measurements of the solvation of salts and electrolytes are relatively straightforward, but it is not possible to separate total solvation free energies into distinct cation and anion contributions without reference to an additional extrathermodynamic assumption. The present work attempts to resolve this difficulty using molecular dynamics simulations with the AMOEBA polarizable force field and perturbation techniques to directly compute absolute solvation free energies for potassium, sodium, and chloride ions in liquid water and formamide. Corresponding calculations are also performed with two widely used nonpolarizable force fields. The simulations with the polarizable force field accurately reproduce in vacuo quantum mechanical results, experimental ion-cluster solvation enthalpies, and experimental solvation free energies for whole salts, while the other force fields do not. The results indicate that calculations with a polarizable force field can capture the thermodynamics of ion solvation and that the solvation free energies of the individual ions differ by several kilocalories from commonly cited values.

An environmental pseudopotential approach to molecular interactions: Implementation in MOLPRO

We present the implementation into the MOLPRO package of a model for the interaction of a central system with its surrounding environment. The properties of a target system enclosed by a noncovalently bound environment or solvent are modeled as those of a system embedded into the effective pseudopotential arising from the exact electrostatic Coulomb potential and the approximated exchange-repulsion potential. For the latter we use the charge-density overlap model, which relates the exchange-repulsion interaction energy between two species with the overlap of their ground-state electron charge densities. The solutions of the modified Hartree-Fock equations for the target system are obtained self-consistently. This way the exchange-induction effects arising from the converged electron-charge density of the embedded system are implicitly included. Inclusion of the correlation effects is provided by the use of post-Hartree-Fock and density-functional techniques available in the MOLPRO package. The computational and conceptual advantages provided by this approach are shown in the calculation of the dipole polarizabilities of halide and chalcogenide anions in different environments.