Thermodynamics of water octamer in a uniform electric field

de Almeida M, Marques JMC, Prudente FV. A thermodynamic view on the microsolvation of ions by rare gas: application to Li+ with argon

Parallel tempering Monte Carlo calculations on the Li⁺Ar_n microsolvation clusters have shown that the two peaks appearing in the heat capacity curve as a function of temperature correspond to the melting of the second and first solvation shells.

Perspectives on external electric fields in molecular simulation: progress, prospects and challenges

In this review, the application of a wide variety of external electric fields in molecular simulation shall be discussed, including time-varying and electromagnetic, as well as the utility and potential impact and prospects for exploitation of such simulations for real-world and industrial end use. In particular, non-equilibrium molecular dynamics will be discussed, as well as challenges in addressing adequate thermostatting and scaling field amplitudes to more experimentally relevant levels. Attention shall be devoted to recent progress and advances in external fields in ab initio molecular simulation and dynamics, as well as elusive challenges thereof (and, to some extent, for molecular dynamics from empirical potentials), such as timescales required to observe low-frequency and intensity field effects. The challenge of deterministic molecular dynamics in external fields in sampling phase space shall be discussed, along with prospects for application of fields in enhanced-sampling simulations. Finally, the application of external electric fields to a wide variety of aqueous, nanoscale and biological systems will be discussed, often motivated by the possibility of exploitation in real-world applications, which serve to underpin our molecular-level understanding of field effects in terms of microscopic mechanisms, and possibly with a view to control thereof.

Exploring energy landscapes: from molecular to mesoscopic systems

We review a comprehensive computational framework to survey the potential energy landscape for systems composed of rigid or partially rigid molecules. Illustrative case studies relevant to a wide range of molecular clusters and soft and condensed matter systems are discussed.

Structure of the floating water bridge and water in an electric field

The floating water bridge phenomenon is a freestanding rope-shaped connection of pure liquid water, formed under the influence of a high potential difference (approximately 15 kV). Several recent spectroscopic, optical, and neutron scattering studies have suggested that the origin of the bridge is associated with the formation of anisotropic chains of water molecules in the liquid. In this work, high energy X-ray diffraction experiments have been performed on a series of floating water bridges as a function of applied voltage, bridge length, and position within the bridge. The two-dimensional X-ray scattering data showed no direction-dependence, indicating that the bulk water molecules do not exhibit any significant preferred orientation along the electric field. The only structural changes observed were those due to heating, and these effects were found to be the same as for bulk water. These X-ray scattering measurements are supported by molecular dynamics (MD) simulations which were performed under electric fields of 10(6) V/m and 10(9) V/m. Directional structure factor calculations were made from these simulations parallel and perpendicular to the E-field. The 10(6) V/m model showed no significant directional-dependence (anisotropy) in the structure factors. The 10(9) V/m model however, contained molecules aligned by the E-field, and had significant structural anisotropy.

Physical properties of small water clusters in low and moderate electric fields

Likely candidates for the lowest minima of water clusters (H(2)O)(N) for N ≤ 20 interacting with a uniform electric field strength in the range E ≤ 0.6 V/Å have been identified using basin-hopping global optimization. Two water-water model potentials were considered, namely TIP4P and the polarizable Dang-Chang potential. The two models produce some consistent results but also exhibit significant differences. The cluster internal energy and dipole moment indicate two varieties of topological transition in the structure of the global minimum as the field strength is increased. The first takes place at low field strengths (0.1 V/Å<E < 0.2 V/Å) and reorganizes the hydrogen-bonds to orient the water permanent dipoles along the field. The second type of transition occurs at larger field strengths (0.3 V/Å<E < 0.5 V/Å) and corresponds to an extensive structural reorganization, where several hydrogen-bonds break as the cluster stretches along the field direction, the larger clusters (N > 10) usually forming helical structures.

Ionization of doped helium nanodroplets: complexes of C60 with water clusters

Water clusters are known to undergo an autoprotonation reaction upon ionization by photons or electron impact, resulting in the formation of (H(2)O)(n)H(3)O(+). Ejection of OH cannot be quenched by near-threshold ionization; it is only partly quenched when clusters are complexed with inert gas atoms. Mass spectra recorded by electron ionization of water-doped helium droplets show that the helium matrix also fails to quench OH loss. The situation changes drastically when helium droplets are codoped with C(60). Charged C(60)-water complexes are predominantly unprotonated; C(60)(H(2)O)(4)(+) and (C(60))(2)(H(2)O)(4)(+) appear with enhanced abundance. Another intense ion series is due to C(60)(H(2)O)(n)OH(+); dehydrogenation is proposed to be initiated by charge transfer between the primary He(+) ion and C(60). The resulting electronically excited C(60)(+*) leads to the formation of a doubly charged C(60)-water complex either via emission of an Auger electron from C(60)(+*), or internal Penning ionization of the attached water complex, followed by charge separation within {C(60)(H(2)O)(n)}(2+). This mechanism would also explain previous observations of dehydrogenation reactions in doped helium droplets. Mass-analyzed ion kinetic energy scans reveal spontaneous (unimolecular) dissociation of C(60)(H(2)O)(n)(+). In addition to the loss of single water molecules, a prominent reaction channel yields bare C(60)(+) for sizes n=3, 4, or 6. Ab initio Hartree-Fock calculations for C(60)-water complexes reveal negligible charge transfer within neutral complexes. Cationic complexes are well described as water clusters weakly bound to C(60)(+). For n=3, 4, or 6, fissionlike desorption of the entire water complex from C(60)(H(2)O)(n)(+) energetically competes with the evaporation of a single water molecule.

Folding of gas-phase polyalanines in a static electric field: alignment, deformations, and polarization effects

Monte Carlo simulations of the temperature-induced unfolding of small gas-phase polyalanines in a static, homogeneous electric field are reported, based on the AMBER ff96 force field. The peptides exhibit a structural transition from the native alpha-helix state to entropically favored beta-sheet conformations, before eventually turning to extended coil at higher temperatures. Upon switching the electric field, the molecules undergo preferential alignment of their dipole moment vector toward the field axis and a shift of the alpha-beta transition to higher temperatures. At higher field strengths (>10(8) V/m) the molecules stretch and the alpha-beta and beta-coil transitions merge. A simple three-state model is shown to account for the observed behavior. Under even higher fields, density functional theory calculations and a polarizable force field both show that electronic rearrangements tend to further increase the dipole moment, polarization effects being approximately half in magnitude with respect to stretching effect. Finally a tentative (temperature, field-strength) phase diagram is sketched.

Water absorption in polyethylene under external electric fields

Monte Carlo simulations of the solubility and structure of water in polyethylene in thermodynamic equilibrium with liquid water were performed in external fields ranging from 2 x 10(5) to 4 x 10(9) V/m. For a given equilibrium temperature and pressure, the water solubility decreases at higher fields. This occurs since it is energetically favorable for water molecules to be in the pure water phase than in the polyethylene matrix at high field strengths, and results in an increased density in the water phase. However, fields relevant to high voltage conduction (in the absence of defects that can lead to large local field strengths) do not change the solubility. In addition, at large fields the number of water clusters decreases for all cluster sizes. The rate of decrease is highest for large clusters, and a larger fraction of water molecules exist as monomers in the polyethylene matrix at high fields. Large fields also cause alignment of the water molecules, which leads to more clusters with linear topologies and hence an increase in the cluster radius of gyration.

Energy landscapes for water clusters in a uniform electric field

The behavior of small water clusters, (H2O)n, n=2-5 and n=8, in a uniform electric field is investigated for three related rigid-body models. Changes in the properties of the low-lying potential energy minima and the rearrangement pathways between them are examined. Results for certain structural transitions are compared with recent ab initio calculations. The models are found to give qualitatively similar trends, and there is some evidence that as the applied field strength is increased the quantitative differences between the models are also reduced.