Sodium Chloride, NaCl/ϵ: New Force Field

Electrokinetic properties of NaCl solution via molecular dynamics simulations with scaled-charge electrolytes

Scaling ionic charges has become an alternative to polarizable force fields for representing indirect charge transfer effects in molecular simulations. In our work, we apply molecular dynamics simulations to investigate the properties of NaCl aqueous solutions in homogeneous and confined media. We compare classical integer- and scaled-charge force fields for the ions. In the bulk, we validate the force fields by computing equilibrium and transport properties and comparing them with experimental data. Integer-charge ions overestimate dielectric saturation and ionic association. Both force fields present an excess in ion-ion correlation, which leads to a deviation in the ionic conductivity at higher ionic strengths. Negatively charged quartz is used to simulate the confinement effect. Electrostatic interactions dominate counter-ion adsorption. Full-charge ions have stronger and more defined adsorption planes. We obtain the electroosmotic mobility of the solution by combining the shear plane location from non-equilibrium simulations with the ionic distribution from equilibrium simulations. From the Helmholtz-Smoluchowski equation, the zeta potential and the streaming potential coupling coefficient are computed. From an atomic-scale perspective, our molecular dynamics simulations corroborate the hypothesis of maximum packing of the Stern layer, which results in a stable and non-zero zeta potential at high salinity. The scaled-charge model representation of both properties is in excellent qualitative and quantitative agreement with experimental data. With our work, we demonstrate how useful and precise simple scaled-charge models for electrolytes can be to represent complex systems, such as the electrical double layer.

Accuracy limit of non-polarizable four-point water models: TIP4P/2005 vs OPC. Should water models reproduce the experimental dielectric constant?

The last generation of four center non-polarizable models of water can be divided into two groups: those reproducing the dielectric constant of water, as OPC, and those significantly underestimating its value, as TIP4P/2005. To evaluate the global performance of OPC and TIP4P/2005, we shall follow the test proposed by Vega and Abascal in 2011 evaluating about 40 properties to fairly address this comparison. The liquid-vapor and liquid-solid equilibria are computed, as well as the heat capacities, isothermal compressibilities, surface tensions, densities of different ice polymorphs, the density maximum, equations of state at high pressures, and transport properties. General aspects of the phase diagram are considered by comparing the ratios of different temperatures (namely, the temperature of maximum density, the melting temperature of hexagonal ice, and the critical temperature). The final scores are 7.2 for TIP4P/2005 and 6.3 for OPC. The results of this work strongly suggest that we have reached the limit of what can be achieved with non-polarizable models of water and that the attempt to reproduce the experimental dielectric constant deteriorates the global performance of the water force field. The reason is that the dielectric constant depends on two surfaces (potential energy and dipole moment surfaces), whereas in the absence of an electric field, all properties can be determined simply from just one surface (the potential energy surface). The consequences of the choice of the water model in the modeling of electrolytes in water are also discussed.

Molecular dynamics simulations of the dielectric constants of salt-free and salt-doped polar solvents

We develop a Stockmayer fluid model that accounts for the dielectric responses of polar solvents (water, MeOH, EtOH, acetone, 1-propanol, DMSO, and DMF) and NaCl solutions. These solvent molecules are represented by Lennard-Jones (LJ) spheres with permanent dipole moments and the ions by charged LJ spheres. The simulated dielectric constants of these liquids are comparable to experimental values, including the substantial decrease in the dielectric constant of water upon the addition of NaCl. Moreover, the simulations predict an increase in the dielectric constant when considering the influence of ion translations in addition to the orientation of permanent dipoles.

A Nanoconfined Water-Ion Coordination Network for Flexible Energy-Dissipation Devices

Water-ion interaction in a nanoconfined environment that deeply constrains spatial freedoms of local atomistic motion with unconventional coupling mechanisms beyond that in a free, bulk state is essential to spark designs of a broad spectrum of nanofluidic devices with unique properties and functionalities. Here, it is reported that the interaction between ions and water molecules in a hydrophobic nanopore forms a coordination network with an interaction density that is nearly fourfold that of the bulk counterpart. Such strong interaction facilitates the connectivity of the water-ion network and is uncovered by corroborating the formation of ion clusters and the reduction of particle dynamics. A liquid-nanopore energy-dissipation system is designed and demonstrated in both molecular simulations and experiments that the formed coordination network controls the outflow of confined electrolytes along with a pressure reduction, capable of providing flexible protection for personnel and devices and instrumentations against external mechanical impact and attack.

Computation of Electrical Conductivities of Aqueous Electrolyte Solutions: Two Surfaces, One Property

In this work, we computed electrical conductivities under ambient conditions of aqueous NaCl and KCl solutions by using the Einstein-Helfand equation. Common force fields (charge q = ±1 e) do not reproduce the experimental values of electrical conductivities, viscosities, and diffusion coefficients. Recently, we proposed the idea of using different charges to describe the potential energy surface (PES) and the dipole moment surface (DMS). In this work, we implement this concept. The equilibrium trajectories required to evaluate electrical conductivities (within linear response theory) were obtained by using scaled charges (with the value q = ±0.75 e) to describe the PES. The potential parameters were those of the Madrid-Transport force field, which accurately describe viscosities and diffusion coefficients of these ionic solutions. However, integer charges were used to compute the conductivities (thus describing the DMS). The basic idea is that although the scaled charge describes the ion-water interaction better, the integer charge reflects the value of the charge that is transported due to the electric field. The agreement obtained with experiments is excellent, as for the first time electrical conductivities (and the other transport properties) of NaCl and KCl electrolyte solutions are described with high accuracy for the whole concentration range up to their solubility limit. Finally, we propose an easy way to obtain a rough estimate of the actual electrical conductivity of the potential model under consideration using the approximate Nernst-Einstein equation, which neglects correlations between different ions.

Why Dissolving Salt in Water Decreases Its Dielectric Permittivity

The dielectric permittivity of salt water decreases on dissolving more salt. For nearly a century, this phenomenon has been explained by invoking saturation in the dielectric response of the solvent water molecules. Herein, we employ an advanced deep neural network (DNN), built using data from density functional theory, to study the dielectric permittivity of sodium chloride solutions. Notably, the decrease in the dielectric permittivity as a function of concentration, computed using the DNN approach, agrees well with experiments. Detailed analysis of the computations reveals that the dominant effect, caused by the intrusion of ionic hydration shells into the solvent hydrogen-bond network, is the disruption of dipolar correlations among water molecules. Accordingly, the observed decrease in the dielectric permittivity is mostly due to increasing suppression of the collective response of solvent waters.

Scaled charges for ions: An improvement but not the final word for modeling electrolytes in water

In this work, we discuss the use of scaled charges when developing force fields for NaCl in water. We shall develop force fields for Na⁺ and Cl^- using the following values for the scaled charge (in electron units): ±0.75, ±0.80, ±0.85, and ±0.92 along with the TIP4P/2005 model of water (for which previous force fields were proposed for q = ±0.85 and q = ±1). The properties considered in this work are densities, structural properties, transport properties, surface tension, freezing point depression, and maximum in density. All the developed models were able to describe quite well the experimental values of the densities. Structural properties were well described by models with charges equal to or larger than ±0.85, surface tension by the charge ±0.92, maximum in density by the charge ±0.85, and transport properties by the charge ±0.75. The use of a scaled charge of ±0.75 is able to reproduce with high accuracy the viscosities and diffusion coefficients of NaCl solutions for the first time. We have also considered the case of KCl in water, and the results obtained were fully consistent with those of NaCl. There is no value of the scaled charge able to reproduce all the properties considered in this work. Although certainly scaled charges are not the final word in the development of force fields for electrolytes in water, its use may have some practical advantages. Certain values of the scaled charge could be the best option when the interest is to describe certain experimental properties.

Where Lennard-Jones Potentials Fail: Iterative Optimization of Ion-Water Pair Potentials Based on Ab Initio Molecular Dynamics Data

The use of the Lennard-Jones (LJ) potential in computer simulations of aqueous electrolyte solutions is widespread. The standard approach is to parametrize LJ potential parameters against thermodynamic solution properties, but problems in representing the local structural and dynamic properties of ion hydration shells remain. The r^-12-term in the LJ potential is responsible for this as it leads to overly repulsive ion-water interactions at short range. As a result, the LJ potential predicts blue-shifted vibrational peaks of the cations' rattling mode and too large negative ion hydration entropies. We demonstrate that cation-water effective pair potentials derived from ab initio MD data have softer short-range repulsions and represent hydration shell properties significantly better. Our findings indicate that replacing the LJ potential with these effective pair potentials offers a promising route to represent thermodynamic solution properties and local interactions of specific ions with nonpolarizable force field models.

Maximum in density of electrolyte solutions: Learning about ion-water interactions and testing the Madrid-2019 force field

In this work, we studied the effect of Li⁺, Na⁺, K⁺, Mg²⁺, and Ca²⁺ chlorides and sulfates on the temperature of maximum density (TMD) of aqueous solutions at room pressure. Experiments at 1 molal salt concentration were carried out to determine the TMD of these solutions. We also performed molecular dynamics simulations to estimate the TMD at 1 and 2 m with the Madrid-2019 force field, which uses the TIP4P/2005 water model and scaled charges for the ions, finding an excellent agreement between experiment and simulation. All the salts studied in this work shift the TMD of the solution to lower temperatures and flatten the density vs temperature curves (when compared to pure water) with increasing salt concentration. The shift in the TMD depends strongly on the nature of the electrolyte. In order to explore this dependence, we have evaluated the contribution of each ion to the shift in the TMD concluding that Na⁺, Ca²⁺, and SO₄ ^2- seem to induce the largest changes among the studied ions. The volume of the system has been analyzed for salts with the same anion and different cations. These curves provide insight into the effect of different ions upon the structure of water. We claim that the TMD of electrolyte solutions entails interesting physics regarding ion-water and water-water interactions and should, therefore, be considered as a test property when developing force fields for electrolytes. This matter has been rather unnoticed for almost a century now and we believe it is time to revisit it.

Freezing point depression of salt aqueous solutions using the Madrid-2019 model

Salt aqueous solutions are relevant in many fields, ranging from biological systems to seawater. Thus, the availability of a force-field that is able to reproduce the thermodynamic and dynamic behavior of salt aqueous solutions would be of great interest. Unfortunately, this has been proven challenging, and most of the existing force-fields fail to reproduce much of their behavior. In particular, the diffusion of water or the salt solubility are often not well reproduced by most of the existing force-fields. Recently, the Madrid-2019 model was proposed, and it was shown that this force-field, which uses the TIP4P/2005 model for water and non-integer charges for the ions, provides a good description of a large number of properties, including the solution densities, viscosities, and the diffusion of water. In this work, we assess the performance of this force-field on the evaluation of the freezing point depression. Although the freezing point depression is a colligative property that at low salt concentrations depends solely on properties of pure water, a good model for the electrolytes is needed to accurately predict the freezing point depression at moderate and high salt concentrations. The coexistence line between ice and several salt aqueous solutions (NaCl, KCl, LiCl, MgCl₂, and Li₂SO₄) up to the eutectic point is estimated from direct coexistence molecular dynamics simulations. Our results show that this force-field reproduces fairly well the experimentally measured freezing point depression with respect to pure water freezing for all the salts and at all the compositions considered.

Anomalous Phase Behaviors of Monolayer NaCl Aqueous Solutions Induced by Effective Coulombic Interactions within Angstrom-Scale Slits

Interests in subnanofluidic devices have called for molecular dynamics (MD) simulation studies of the thermodynamic behavior of monolayer salt solution within angstrom-scale slits. However, it still remains a grand challenge to accurately describe the Coulombic interactions by incorporating the effects of charge transfer and electronic dielectric screening. Herein, by using the electronic continuum model, where the effective ion charges are fine-tuned with a scaling factor of λ, we present simulation evidence that the effective Coulombic interactions among Na⁺/Cl^- ions can strongly affect the behavior of monolayer ionic aqueous solution. Our microsecond-scale MD simulations show that only the counterions with moderate effective charges (0.3 ≤ λ ≤ 0.8) can dissolve in monolayer water, whereas the high effective charges (λ ≥ 0.85) induce ions to assemble into monolayer nanocrystals, and ions with the low effective charges (λ ≤ 0.2) exhibit gas-like nanobubble. These findings could provide deeper insights into the physical chemistry behind subnanofluidic iontronic devices.

Activity coefficients of aqueous electrolytes from implicit-water molecular dynamics simulations

We obtain activity coefficients in NaCl and KCl solutions from implicit-water molecular dynamics simulations, at 298.15 K and 1 bar, using two distinct approaches. In the first approach, we consider ions in a continuum with constant relative permittivity (ɛ_r) equal to that of pure water; in the other approach, we take into account the concentration-dependence of ɛ_r, as obtained from explicit-water simulations. Individual ion activity coefficients (IIACs) are calculated using gradual insertion of single ions with uniform neutralizing backgrounds to ensure electroneutrality. Mean ionic activity coefficients (MIACs) obtained from the corresponding IIACs in simulations with constant ɛ_r show reasonable agreement with experimental data for both salts. Surprisingly, large systematic negative deviations are observed for both IIACs and MIACs in simulations with concentration-dependent ɛ_r. Our results suggest that the absence of hydration structure in implicit-water simulations cannot be compensated by correcting for the concentration-dependence of the relative permittivity ɛ_r. Moreover, even in simulations with constant ɛ_r for which the calculated MIACs are reasonable, the relative positioning of IIACs of anions and cations is incorrect for NaCl. We conclude that there are severe inherent limitations associated with implicit-water simulations in providing accurate activities of aqueous electrolytes, a finding with direct relevance to the development of electrolyte theories and to the use and interpretation of implicit-solvent simulations.

Structural and dynamic properties of some aqueous salt solutions

Aqueous salt solutions are utilized and encountered in wide-ranging technological applications and natural settings. Towards improved understanding of the effect of salts on the dynamic properties of such systems, dilute aqueous salt solutions (up to 1 molar concentration) are investigated here, via experiments and molecular simulations. Four salts are considered: sodium chloride, for which published results are readily available for comparison, ammonium acetate, barium acetate and barium nitrate, for which published data are scarce. In the present work, molecular dynamics (MD) simulations are conducted to quantify viscosity and water self-diffusion coefficients, together with rheometry and Pulsed Field Gradient Spin Echo (PFGSE)-NMR experiments for validation. Simulation predictions are consistent with experimental observations in terms of trend and magnitude of salt-specific effects. Combining insights from the approaches considered, an interpretation of the results is proposed whereby the capacity of salts to influence bulk dynamics arises from their molecular interfacial area and strength of interaction with first hydration-shell water molecules. For the concentration range investigated, the interpretation could be useful in formulating aqueous systems for applications including the manufacturing of advanced catalysts.

Molecular fluid flow in MoS2 nanoporous membranes and hydrodynamics interactions

We study the impact of the induced pressure fields on the water flow and salt rejection in nanopores produced in MoS₂ membranes. We observe that the water permeability and the salt rejection are not impacted by the distance between the pores. This result contradicts the continuous fluid mechanics calculations in microfilters, which indicates the existence of hydrodynamic interactions between adjacent pores that increase the water mobility. Our results suggest that at this nanoscale, the hydrodynamic interactions do not affect the water mobility through nanopores.

A practical guide to biologically relevant molecular simulations with charge scaling for electronic polarization

Molecular simulations can elucidate atomistic-level mechanisms of key biological processes, which are often hardly accessible to experiment. However, the results of the simulations can only be as trustworthy as the underlying simulation model. In many of these processes, interactions between charged moieties play a critical role. Current empirical force fields tend to overestimate such interactions, often in a dramatic way, when polyvalent ions are involved. The source of this shortcoming is the missing electronic polarization in these models. Given the importance of such biomolecular systems, there is great interest in fixing this deficiency in a computationally inexpensive way without employing explicitly polarizable force fields. Here, we review the electronic continuum correction approach, which accounts for electronic polarization in a mean-field way, focusing on its charge scaling variant. We show that by pragmatically scaling only the charged molecular groups, we qualitatively improve the charge-charge interactions without extra computational costs and benefit from decades of force field development on biomolecular force fields.

Salt parameterization can drastically affect the results from classical atomistic simulations of water desalination by MoS2 nanopores

Water scarcity is a reality in our world, and scenarios predicted by leading scientists in this area indicate that it will worsen in the next decades. However, new technologies based on low-cost seawater desalination can prevent the worst scenarios, providing fresh water for humanity. With this goal, membranes based on nanoporous materials have been suggested in recent years. One of the materials suggested is MoS₂, and classical Molecular Dynamics (MD) simulation is one of the most powerful tools to explore these nanomaterials. However, distinct force fields employed in MD simulations are parameterized based on distinct experimental quantities. In this paper, we compare two models of salt that were built based on distinct properties of water-salt mixtures. One model fits the hydration free energy and lattice properties, and the second fits the crystal density and the density and the dielectric constant of water and salt mixtures. To compare the models, MD simulations for salty water flow through nanopores of two sizes were used - one pore big enough to accommodate hydrated ions, and one smaller in which the ion has to dehydrate to enter - and two rigid water models from the TIP4P family - TIP4P/2005 and TIP4P/ε. Our results indicate that the water permeability and salt rejection by the membrane are more influenced by the salt model than by the water model, especially for the narrow pore. In fact, completely distinct mechanisms were observed, and they are related to the characteristics employed in the ion model parameterization. The results show that not only can the water model influence the outcomes, but the ion model plays a crucial role when the pore is small enough.

Dielectric constant and density of aqueous alkali halide solutions by molecular dynamics: A force field assessment

The concentration dependence of the dielectric constant and the density of 11 aqueous alkali halide solutions (LiCl, NaCl, KCl, RbCl, CsCl, LiI, NaI, KI, CsI, KF, and CsF) is investigated by molecular simulation. Predictions using eight non-polarizable ion force fields combined with the TIP4P/ε water model are compared to experimental data. The influence of the water model and the temperature on the results for the NaCl brine are also addressed. The TIP4P/ε water model improves the accuracy of dielectric constant predictions compared to the SPC/E water model. The solution density is predicted well by most ion models. Almost all ion force fields qualitatively capture the decline of the dielectric constant with the increase of concentration for all solutions and with the increase of temperature for NaCl brine. However, the sampled dielectric constant is mostly in poor quantitative agreement with experimental data. These results are related to the microscopic solution structure, ion pairing, and ultimately the force field parameters. Ion force fields with excessive contact ion pairing and precipitation below the experimental solubility limit generally yield higher dielectric constant values. An adequate reproduction of the experimental solubility limit should therefore be a prerequisite for further investigations of the dielectric constant of aqueous electrolyte solutions by molecular simulation.

On the transferability of ion parameters to the TIP4P/2005 water model using molecular dynamics simulations

Countless molecular dynamics studies have relied on available ion and water force field parameters to model aqueous electrolyte solutions. The TIP4P/2005 model has proven itself to be among the best rigid water force fields, whereas many of the most successful ion parameters were optimized in combination with SPC/E, TIP3P, or TIP4P/Ew water. Many researchers have combined these ions with TIP4P/2005, hoping to leverage the strengths of both parameter sets. To assess if this widely used approach is justified and to provide a guide in selecting ion parameters, we investigated the transferability of various commonly used monovalent and multivalent ion parameters to the TIP4P/2005 water model. The transferability is evaluated in terms of ion hydration free energy, hydration radius, coordination number, and self-diffusion coefficient at infinite dilution. For selected ion parameters, we also investigated density, ion pairing, chemical potential, and mean ionic activity coefficients at finite concentrations. We found that not all ions are equally transferable to TIP4P/2005 without compromising their performance. In particular, ions optimized for TIP3P water were found to be poorly transferable to TIP4P/2005, whereas ions optimized for TIP4P/Ew water provided nearly perfect transferability. The latter ions also showed good overall agreement with experimental values. The one exception is that no combination of ion parameters and water model considered here was found to accurately reproduce experimental self-diffusion coefficients. Additionally, we found that cations optimized for SPC/E and TIP3P water displayed consistent underpredictions in the hydration free energy, whereas anions consistently overpredicted the hydration free energy.

Correction: Salt parameterization can drastically affect the results from classical atomistic simulations of water desalination by MoS2 nanopores

Correction for ‘Salt parameterization can drastically affect the results from classical atomistic simulations of water desalination by MoS₂ nanopores’ by João P. K. Abal et al., Phys. Chem. Chem. Phys., 2020, 22, 11053–11061, DOI: 10.1039/d0cp00484g.

Dielectric Decrement for Aqueous NaCl Solutions: Effect of Ionic Charge Scaling in Nonpolarizable Water Force Fields

We investigate the dielectric constant and the dielectric decrement of aqueous NaCl solutions by means of molecular dynamic simulations. We thereby compare the performance of four different force fields and focus on disentangling the origin of the dielectric decrement and the influence of scaled ionic charges, as often used in nonpolarizable force fields to account for the missing dynamic polarizability in the shielding of electrostatic ion interactions. Three of the force fields showed excessive contact ion pair formation, which correlates with a reduced dielectric decrement. In spite of the fact that the scaling of charges only weakly influenced the average polarization of water molecules around an ion, the rescaling of ionic charges did influence the dielectric decrement, and a close-to-linear relation of the slope of the dielectric constant as a function of concentration with the ionic charge was found.

A force field of Li+, Na+, K+, Mg2+, Ca2+, Cl-, and SO4 2- in aqueous solution based on the TIP4P/2005 water model and scaled charges for the ions

In this work, a force field for several ions in water is proposed. In particular, we consider the cations Li⁺, Na⁺, K⁺, Mg²⁺, and Ca²⁺ and the anions Cl^- and SO₄ ^2-. These ions were selected as they appear in the composition of seawater, and they are also found in biological systems. The force field proposed (denoted as Madrid-2019) is nonpolarizable, and both water molecules and sulfate anions are rigid. For water, we use the TIP4P/2005 model. The main idea behind this work is to further explore the possibility of using scaled charges for describing ionic solutions. Monovalent and divalent ions are modeled using charges of 0.85 and 1.7, respectively (in electron units). The model allows a very accurate description of the densities of the solutions up to high concentrations. It also gives good predictions of viscosities up to 3 m concentrations. Calculated structural properties are also in reasonable agreement with the experiment. We have checked that no crystallization occurred in the simulations at concentrations similar to the solubility limit. A test for ternary mixtures shows that the force field provides excellent performance at an affordable computer cost. In summary, the use of scaled charges, which could be regarded as an effective and simple way of accounting for polarization (at least to a certain extend), improves the overall description of ionic systems in water. However, for purely ionic systems, scaled charges will not adequately describe neither the solid nor the melt.

Systematically improved melting point prediction: a detailed physical simulation model is required

Our detailed physical simulation model allows for an unprecedented and systematically improved prediction of melting points of alkali halides.

Incremental update of electrostatic interactions in adaptively restrained particle simulations

The computation of long-range potentials is one of the demanding tasks in Molecular Dynamics. During the last decades, an inventive panoply of methods was developed to reduce the CPU time of this task. In this work, we propose a fast method dedicated to the computation of the electrostatic potential in adaptively restrained systems. We exploit the fact that, in such systems, only some particles are allowed to move at each timestep. We developed an incremental algorithm derived from a multigrid-based alternative to traditional Fourier-based methods. Our algorithm was implemented inside LAMMPS, a popular molecular dynamics simulation package. We evaluated the method on different systems. We showed that the new algorithm's computational complexity scales with the number of active particles in the simulated system, and is able to outperform the well-established Particle Particle Particle Mesh (P3M) for adaptively restrained simulations. © 2018 Wiley Periodicals, Inc.

Collision-Induced Dissociation of Electrosprayed NaCl Clusters: Using Molecular Dynamics Simulations to Visualize Reaction Cascades in the Gas Phase

Infusion of NaCl solutions into an electrospray ionization (ESI) source produces [Na_(n+1)Cl _n ]⁺ and other gaseous clusters. The n = 4, 13, 22 magic number species have cuboid ground state structures and exhibit elevated abundance in ESI mass spectra. Relatively few details are known regarding the mechanisms whereby these clusters undergo collision-induced dissociation (CID). The current study examines to what extent molecular dynamics (MD) simulations can be used to garner insights into the sequence of events taking place during CID. Experiments on singly charged clusters reveal that the loss of small neutrals is the dominant fragmentation pathway. MD simulations indicate that the clusters undergo extensive structural fluctuations prior to decomposition. Consistent with the experimentally observed behavior, most of the simulated dissociation events culminate in ejection of small neutrals ([NaCl] _i , with i = 1, 2, 3). The MD data reveal that the prevalence of these dissociation channels is linked to the presence of short-lived intermediates where a relatively compact core structure carries a small [NaCl] _i protrusion. The latter can separate from the parent cluster via cleavage of a single Na-Cl contact. Fragmentation events of this type are kinetically favored over other dissociation channels that would require the quasi-simultaneous rupture of multiple electrostatic contacts. The CID behavior of NaCl cluster ions bears interesting analogies to that of collisionally activated protein complexes. Overall, it appears that MD simulations represent a valuable tool for deciphering the dissociation of noncovalently bound systems in the gas phase. Graphical Abstract ᅟ.