Energy landscapes for water clusters in a uniform electric field

Understanding the Underlying Field Evaporation Mechanism of Pure Water Tips in High Electrical Fields

We investigated the field evaporation process of frozen water in atom probe tomography (APT) by density functional simulations. In previous experiments, a strong tailing effect was observed for peaks caused by the molecular structure (H₂O)_nH⁺, in contrast to other peaks. In purely field-induced and thermally assisted evaporation simulations, we found that chains of protonated water molecules were pulled out of the dielectric surface by up to 6 Å, which are stable over a wide range of field strengths. Therefore, the resulting water clusters experience only part of the acceleration after evaporation compared to molecules evaporating directly from the surface and, thus, exhibit an energy deficit, which explains the tailing effect. Our simulations provide new insight into the complex evaporation behavior of water in high electrical fields and reveal possibilities for adapting the existing reconstruction algorithms.

Solvent Organization and Rate Regulation of a Menshutkin Reaction by Oriented External Electric Fields are Revealed by Combined MD and QM/MM Calculations

When and how do external electric fields (EEFs) lead to catalysis in the presence of a (polar or nonpolar) solvent? This is the question that is addressed here using a combination of molecular dynamics (MD) simulations, quantum mechanical/molecular mechanical calculations with EEF, and quantum mechanical/(local) electric field calculations. The paper focuses on a model reaction, the Menshutkin reaction between CH3I and pyridine in three solvents of varying polarity. Using MD simulations, we find that the EEF causes the solvent to undergo organization; the solvent molecules gradually align with the applied field as the field strength increases. The collective orientation of the solvent molecules modifies the electrostatic environment around the Menshutkin species and induces a global electric field pointing in the opposite direction of the applied EEF. The combination of these two entangled effects leads to partial or complete screening of the EEF, with the extent of screening being proportional to the polarity/polarizability of the solvent. Nevertheless, we find that catalysis of the Menshutkin reaction inevitably emerges once the EEF exceeds the opposing field of the organizing solvent, i.e., once polarization of the Menshutkin complex is observed to set in. Overall, our analysis provides a lucid and pictorial interpretation of the behavior of solutions in the presence of EEFs and indicates that EEF-mediated catalysis should, in principle, be feasible in bulk setups, especially for nonpolar and mildly polar solvents. By application of the charge-transfer paradigm, it is shown that the emergence of OEEF catalysis in solution can be generalized to other reactions as well.

Anhydrous octyl-glucoside phase transition from lamellar to isotropic induced by electric and magnetic fields

A static deuterium nuclear magnetic resonance (²HNMR) technique (magnetic field, B = 7.05 T) was employed to monitor the thermotropic lamellar phase of the anhydrous 1:1 mixture sample of octyl-b-D-glucoside (βOG) and that of partially deuterium labelled at the alpha position on the chain, i.e.,βOG-d₂ In the absence of an electric field, the ²H NMR spectrum of the mixture gives a typical quadrupolar doublet representing the aligned lamellar phase. Upon heating to beyond the clearing temperature at 112 °C, this splitting converts to a single line expected for an isotropic phase. Simultaneous application of magnetic and electric fields (E = 0.4 MV/m) at 85 °C in the lamellar phase, whose direction was set to be parallel or perpendicular to the magnetic field, resulted in the change of the doublet into a single line and this recovers to the initial doublet with time for both experimental geometries. This implies E- and B-field-induced phase transitions from the lamellar to an isotropic phase and a recovery to the lamellar phase again with time. Moreover, these phase transformations are accompanied by a transient current. A similar observation was made in a computational study when an electric field was applied to a water cluster system. Increasing the field strength distorts the water cluster and weakens its hydrogen bonds leading to a structural breakdown beyond a threshold field-strength. Therefore, we suggest the observed field-induced transition is likely due to a structure change of the βOG lamellar assembly caused by the field effect and not due to Joule heating.

Exploring experimental fitness landscapes for chemical synthesis and property optimization

The topology of experimental fitness landscapes for chemical optimization objectives is assessed through svr-based HDMR modeling.

Internal electric fields in small water clusters [(H2O)n; n = 2-6]

The electric field experienced by a water molecule within a water cluster depends on its position relative to the rest of the water molecules. The stabilization energies and the red-shifts in the donor O-H stretching vibrations in the water clusters increase with the cluster size concomitant with the increase in the electric field experienced by the donor O-H of a particular water molecule due to the hydrogen bonding network. The red-shifts in O-H stretching frequencies show a spread of about ±100 cm(-1) against the corresponding electric fields. Deviations from linearity were marked in the region of 100-160 MV cm(-1), which can be attributed to the strain in the hydrogen bonding network, especially for structures with DDAA and DDA motifs. The linear Stark effect holds up to 200 MV cm(-1) of internal electric field for the average red-shifts in the O-H stretching frequencies, with a Stark tuning rate of 2.4 cm(-1) (MV cm(-1))(-1), suggesting the validity of the classical model in small water clusters.

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.

A theoretical study of the dissociation of the sI methane hydrate induced by an external electric field

Molecular dynamics simulations in the equilibrium isobaric-isothermal (NPT) ensemble were used to examine the strength of an external electric field required to dissociate the methane hydrate sI structure. The water molecules were modeled using the four-site TIP4P/Ice analytical potential and methane was described as a simple Lennard-Jones interaction site. A series of simulations were performed at T = 260 K with P = 80 bars and at T = 285 K with P = 400 bars with an applied electric field ranging from 1.0 V nm(-1) to 5.0 V nm(-1). For both (T,P) conditions, applying a field greater than 1.5 V nm(-1) resulted in the orientation of the water molecules such that an ice Ih-type structure was formed, from which the methane was segregated. When the simulations were continued without the external field, the ice-like structures became disordered, resulting in two separate phases: gas methane and liquid water.

Exploring Energy Landscapes: Metrics, Pathways, and Normal-Mode Analysis for Rigid-Body Molecules

We present new methodology for exploring the energy landscapes of molecular systems, using angle-axis variables for the rigid-body rotational coordinates. The key ingredient is a distance measure or metric tensor, which is invariant to global translation and rotation. The metric is used to formulate a generalized nudged elastic band method for calculating pathways, and a full prescription for normal-mode analysis is described. The methodology is tested by mapping the potential energy and free energy landscape of the water octamer, described by the TIP4P potential.

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.

Quasi-Continuous Interpolation Scheme for Pathways between Distant Configurations

A quasi-continuous interpolation (QCI) scheme is introduced for characterizing physically realistic initial pathways from which to initiate transition state searches and construct kinetic transition networks. Applications are presented for peptides, proteins, and a morphological transformation in an atomic cluster. The first step in each case involves end point alignment, and we describe the use of a shortest augmenting path algorithm for optimizing permutational isomers. The QCI procedure then employs an interpolating potential, which preserves the covalent bonding framework for the biomolecules and includes repulsive terms between unconstrained atoms. This potential is used to identify an interpolating path by minimizing contributions from a connected set of images, including terms corresponding to minima in the interatomic distances between them. This procedure detects unphysical geometries in the line segments between images. The most difficult cases, where linear interpolation would involve chain crossings, are treated by growing the structure an atom at a time using the interpolating potential. To test the QCI procedure, we carry through a series of benchmark calculations where the initial interpolation is coupled to explicit transition state searches to produce complete pathways between specified local minima.

The Activation-Relaxation Technique: ART Nouveau and Kinetic ART

The evolution of many systems is dominated by rare activated events that occur on timescale ranging from nanoseconds to the hour or more. For such systems, simulations must leave aside the full thermal description to focus specifically on mechanisms that generate a configurational change. We present here the activation relaxation technique (ART), an open-ended saddle point search algorithm, and a series of recent improvements to ART nouveau and kinetic ART, an ART-based on-the-fly off-lattice self-learning kinetic Monte Carlo method.

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.

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.

Multiscale approach to explore the potential energy surface of water clusters (H2O)nn <or= 8

We propose a multiscale method to explore the energy landscape of water clusters. An asynchronous genetic algorithm is employed to explore the potential energy surface (PES) of OSS2 and TTM2.1-F models. Local minimum structures are collected on the fly, and the ultrafast shape recognition algorithm was used to remove duplicate structures. These structures are then refined at the B3LYP/6-31+G* level. The number of distinct local minima we found (21, 76, 369, 1443, and 3563 isomers for n = 4-8, respectively) reflects the complexity of the PES of water clusters.