Crystal Structure Prediction via Basin-Hopping Global Optimization Employing Tiny Periodic Simulation Cells, with Application to Water-Ice.

A crystal structure prediction algorithm for use in periodic boundary conditions with empirical rigid models is presented, which employs (i) unrestricted cutoff radii for the real-space interactions, thus allowing the treatment of even very small unit cells, and (ii) a global-optimization algorithm based on the basin-hopping method of Wales et al. (D. J. Wales and J. P. K. Doye, J. Phys. Chem. A 1997, 101, 5111). The algorithm is then applied to the TIP4P model of water (W. L. Jorgensen et al., J. Chem. Phys. 1983, 79, 926.) in order to find the lowest enthalpy water-ice crystalline structures in the pressure region 0-8000 bar, in unit cells holding in the range of 1-16 molecules, and a database of the 10 lowest enthalpy structures found at pressures 0, 4000, and 8000 bar is presented. The algorithm finds many of the ice polymorphs and, in particular, finds that the lowest energy structure at zero pressure is almost exactly tied between an ice Ic (cubic ice) and ice Ih (hexagonal ice) structure, having near-identical energies.

Controlling ionic conductivity through transprotein electropores in human aquaporin 4: a non-equilibrium molecular-dynamics study.

Electroporation is a matter of intensive ongoing research interest, and a much-neglected topic in trans-membrane proteins, particularly in view of such promising potential applications in medicine and biotechnology. In particular, selected such novel and exciting applications are predicated on controlling ionic conductivity through electro-pores. Here, we scrutinise the mechanisms of ions' electric conductivity, by means of structural rearrangements, through quasi-stable electro-pores through human-AQP4 as a well-representative prototype of trans-membrane ionic conduction, achieving exquisite control over ionic permeability manipulated by the application of intense static electric fields.

A New Relatively Simple Approach to Multipole Interactions in Either Spherical Harmonics or Cartesians, Suitable for Implementation into Ewald Sums.

We present a novel derivation of the multipole interaction (energies, forces and fields) in spherical harmonics, which results in an expression that is able to exactly reproduce the results of earlier Cartesian formulations. Our method follows the derivations of Smith (W. Smith, CCP5 Newsletter 1998, 46, 18.) and Lin (D. Lin, J. Chem. Phys. 2015, 143, 114115), who evaluate the Ewald sum for multipoles in Cartesian form, and then shows how the resulting expressions can be converted into spherical harmonics, where the conversion is performed by establishing a relation between an inner product on the space of symmetric traceless Cartesian tensors, and an inner product on the space of harmonic polynomials on the unit sphere. We also introduce a diagrammatic method for keeping track of the terms in the multipole interaction expression, such that the total electrostatic energy can be viewed as a 'sum over diagrams', and where the conversion to spherical harmonics is represented by 'braiding' subsets of Cartesian components together. For multipoles of maximum rank n, our algorithm is found to have scaling of n 3.7 vs. n 4.5 for our most optimised Cartesian implementation.

Transprotein-Electropore Characterization: A Molecular Dynamics Investigation on Human AQP4.

Electroporation characterization is a topic of intensive interest probed by extensive ongoing research efforts. Usually, these studies are carried out on lipid-bilayer electroporation. Surprisingly, the possibility of water-channel electropore formation across transmembrane proteins themselves, particularly in view of such a promising application, has not yet been elucidated. The present work examines the geometrical and kinetic aspects of electropores and their stability in such a protein milieux (as opposed through the phospholipid membranes) in depth, by means of scrutiny of such a process in human-AQP4 as a well-representative prototype. The residues forming the electropore's walls, organized in loops, reveal the formation mechanism by their dipole alignment and translational response in response to applied axial electric fields in nonequilibrium molecular dynamics simulation. The magnitude of sustaining electric fields (keeping a stable electropore open) were determined. This suggests that transmembrane proteins could play a central role in electroporation applications, e.g., in medicine and biotechnology.

Does Local Structure Bias How a Crystal Nucleus Evolves?

The broad scientific and technological importance of crystallization has led to significant research probing and rationalizing crystal nucleation processes. Previous work has generally neglected the possibility of the molecular-level dynamics of individual crystal nuclei coupling to local structures. However, recent experimental work has conjectured that this can occur. Therefore, to address a deficiency in scientific understanding of crystallization, we have probed the nucleation of prototypical single and multicomponent crystals (specifically, ice and mixed gas hydrates). We establish that local structures can bias the evolution of nascent crystal phases on a nanosecond time scale by, for example, promoting the appearance or disappearance of specific crystal motifs and thus reveal a new facet of crystallization behavior. Moreover, we demonstrate structural biases are likely present during crystallization processes beyond ice and gas hydrate formation. Structurally biased dynamics are a lens for understanding existing computational and experimental results while pointing to future opportunities.

Electro-nucleation of water nano-droplets in No Man's Land to fault-free ice Ic.

Elucidating water-to-ice freezing, especially in "No Man's Land" (150 K < T < 235 K), is fundamentally important (e.g., predicting upper-troposphere cirrus-cloud formation) - and elusive. An oft-neglected aspect of tropospheric ice-crystallite formation lies in inevitably-present electric fields' role. Exploring nucleation in No Man's Land is technically demanding, owing to rapid nucleation rates, to mention nothing of difficulties of applying relevant electric fields thereto. Here, we tackle these intriguing open questions, via non-equilibrium molecular-dynamics simulations of sub-microsecond formation of rhombus-shaped ice Ic nano-crystallites from aggressively-quenched supercooled water nano-droplets in the gas phase, in external static electric fields. We explore droplets' nano-confined geometries and the entropic-ordering agent of external electric fields as a means of realising cubic-ice formation, especially with very few stacking faults and defects.

Study of hydrogen-molecule guests in type II clathrate hydrates using a force-matched potential model parameterised from ab initio molecular dynamics.

The force-matching method has been applied to parameterise an empirical potential model for water-water and water-hydrogen intermolecular interactions for use in clathrate-hydrate simulations containing hydrogen guest molecules. The underlying reference simulations constituted ab initio molecular dynamics (AIMD) of clathrate hydrates with various occupations of hydrogen-molecule guests. It is shown that the resultant model is able to reproduce AIMD-derived free-energy curves for the movement of a tagged hydrogen molecule between the water cages that make up the clathrate, thus giving us confidence in the model. Furthermore, with the aid of an umbrella-sampling algorithm, we calculate barrier heights for the force-matched model, yielding the free-energy barrier for a tagged molecule to move between cages. The barrier heights are reasonably large, being on the order of 30 kJ/mol, and are consistent with our previous studies with empirical models [C. J. Burnham and N. J. English, J. Phys. Chem. C 120, 16561 (2016) and C. J. Burnham et al., Phys. Chem. Chem. Phys. 19, 717 (2017)]. Our results are in opposition to the literature, which claims that this system may have very low barrier heights. We also compare results to that using the more ad hoc empirical model of Alavi et al. [J. Chem. Phys. 123, 024507 (2005)] and find that this model does very well when judged against the force-matched and ab initio simulation data.

Electro-suppression of water nano-droplets' solidification in no man's land: Electromagnetic fields' entropic trapping of supercooled water.

Understanding water solidification, especially in "No Man's Land" (NML) (150 K < T < 235 K) is crucially important (e.g., upper-troposphere cloud processes) and challenging. A rather neglected aspect of tropospheric ice-crystallite formation is inevitably present electromagnetic fields' role. Here, we employ non-equilibrium molecular dynamics of aggressively quenched supercooled water nano-droplets in the gas phase under NML conditions, in externally applied electromagnetic (e/m) fields, elucidating significant differences between effects of static and oscillating fields: although static fields induce "electro-freezing," e/m fields exhibit the contrary - solidification inhibition. This anti-freeze action extends not only to crystal-ice formation but also restricts amorphisation, i.e., suppression of low-density amorphous ice which forms otherwise in zero-field NML environments. E/m-field applications maintain water in the deeply supercooled state in an "entropic trap," which is ripe for industrial impacts in cryo-freezing, etc.

Human aquaporin 4 gating dynamics under axially oriented electric-field impulses: A non-equilibrium molecular-dynamics study.

Human aquaporin 4 has been studied using non-equilibrium molecular dynamics simulations in the absence and presence of pulses of external electric fields. The pulses were 100 ns in duration and 0.005-0.015 V/Å in intensity acting along the pores' axes. Water diffusivity and the dipolar response of various residues of interest within the pores have been studied. Results show relatively little change in levels of water permeability per se within aquaporin channels during axially oriented field impulses, although care must be taken with regard to statistical certainty. However, the spatial variation of water permeability vis-à-vis electric-field intensity within the milieu of the channels, as revealed by heterogeneity in diffusivity-map gradients, indicates the possibility of somewhat enhanced diffusivity, owing to several residues being affected substantially by external fields, particularly for HIS 201 and 95 and ILE 93. This has the effect of increasing slightly intra-pore water diffusivity in the "pore-mouths" locale, albeit rendering it more spatially uniform overall vis-à-vis zero-field conditions (via manipulation of the selectivity filter).

Electropumping of Water Through Human Aquaporin 4 by Circularly Polarized Electric Fields: Dramatic Enhancement and Control Revealed by Non-Equilibrium Molecular Dynamics.

An extensive suite of nonequilibrium molecular-dynamics (NEMD) simulations have been performed for ∼60 ns of human aquaporin 4 in externally applied circularly polarized (CP) electric fields, applied axially along channels. These external fields were 0.05 V/Å in intensity and 100 GHz in frequency. This has the effect of "electro-pumping" the water through the pores as prototypical biochannels, from conversion of molecules' spin angular momentum to linear momentum in the asymmetric heterogeneous-frictional environment of the pores, thus inducing overall net flow. Water's osmotic permeability was enhanced very substantially (doubled) vis-à-vis the zero-field case. This raises the tantalizing possibility of CP-field-mediated control of water permeability in aquaporins, or other biological (or biomimetic) channels as a potential viable and competitive water-treatment technology.

Quantum and classical inter-cage hopping of hydrogen molecules in clathrate hydrate: temperature and cage-occupation effects.

Building on our previous work (J. Phys. Chem. C, 2016, 120, 16561), using an empirical model we run both classical and path-integral molecular dynamics simulations for a type II clathrate hydrate containing different amounts of guest H2 molecules from 1 to 5 molecules per large cage, with results presented at temperatures of 50, 100 K and 200 K. We present results for the density isosurfaces of the guest molecules at all different occupations and temperatures, showing how the density approaches the perfect tetrahedral structure which has been found for the n = 4 case in which each molecule sits on the vertex of a tetrahedron about the centre of each large cage. We calculate free-energy profiles of the molecules over the volume interior to the cage, and using umbrella sampling, we also calculate the free energy barrier for the molecule to hop between cages. We show that this barrier reduces almost linearly for n = 1-3 molecules per large cage, but becomes larger than would be expected (from extrapolation) for the n = 4 case, with this departure from linearity becoming larger at lower temperatures. We show that, perhaps counter-intuitively, these barriers tend to increase with raising the temperature, and also counter-intuitively that quantisation of the nuclei acts to increase the barriers. Finally, for the n = 4 case, a comparison is made between the empirical model results and those from an ab initio molecular dynamics calculation, which shows that qualitative agreement exists between the two models.

Human Aquaporin 4 Gating Dynamics under Perpendicularly-Oriented Electric-Field Impulses: A Molecular Dynamics Study.

Human aquaporin 4 has been studied using molecular dynamics (MD) simulations in the absence and presence of pulses of external static electric fields. The pulses were 10 ns in duration and 0.012-0.065 V/Å in intensity acting along both directions perpendicular to the pores. Water permeability and the dipolar response of all residues of interest (including the selectivity filter) within the pores have been studied. Results showed decreased levels of water osmotic permeability within aquaporin channels during orthogonally-oriented field impulses, although care must be taken with regard to statistical certainty. This can be explained observing enhanced "dipolar flipping" of certain key residues, especially serine 211, histidine 201, arginine 216, histidine 95 and cysteine 178. These residues are placed at the extracellular end of the pore (serine 211, histidine 201, and arginine 216) and at the cytoplasm end (histidine 95 and cysteine 178), with the key role in gating mechanism, hence influencing water permeability.

Study of clathrate hydrates via equilibrium molecular-dynamics simulation employing polarisable and non-polarisable, rigid and flexible water models.

Equilibrium molecular-dynamics (MD) simulations have been performed on metastable sI and sII polymorphs of empty hydrate lattices, in addition to liquid water and ice Ih. The non-polarisable TIP4P-2005, simple point charge model (SPC), and polarisable Thole-type models (TTM): TTM2, TTM3, and TTM4 water models were used in order to survey the differences between models and to see what differences can be expected when polarisability is incorporated. Rigid and flexible variants were used of each model to gauge the effects of flexibility. Power spectra are calculated and compared to density-of-states spectra inferred from inelastic neutron scattering (INS) measurements. Thermodynamic properties were also calculated, as well as molecular-dipole distributions. It was concluded that TTM models offer optimal fidelity vis-à-vis INS spectra, together with thermodynamic properties, with the flexible TTM2 model offering optimal placement of vibrational modes.

Communication: Librational dynamics in water, sI and sII clathrate hydrates, and ice Ih: Molecular-dynamics insights.

Equilibrium molecular-dynamics simulations have been performed for liquid water, and on metastable sI and sII polymorphs of empty hydrate lattices, in addition to ice Ih, in order to study the dynamical properties of librational motion (rotation oscillation) depicted by protons in water molecules. In particular, hydrate lattices were found to display prominent "bifurcated" features, or peaks, at circa 70 and 80-95 meV (or ∼560 and 640-760 cm(-1), respectively), also displayed by ice, in essentially quantitative agreement with experimental neutron-scattering data. However, observed differences in dispersion between these librational modes between these two structures (both hydrate polymorphs vis-à-vis ice), owing primarily to density effects, have been decomposed into contributions arising from angular-velocity dynamics about axes in the local molecular frame of water molecules, with in-plane "wagging" and "twisting" rationalising one mode at ∼70 meV, and out-of-plane motion for the higher-frequency band. This was confirmed explicitly by a type of de facto normal-mode analysis, in which only immediate layers of water molecules about the one under consideration were allowed to move. In contrast, liquid water displayed no marked preference for such local in- or out-of-plane modes characterising librational motion, owing to the marked absence of rigid, pentamers or hexamers therein.

The properties of ion-water clusters. II. Solvation structures of Na+, Cl-, and H+ clusters as a function of temperature.

Ion-water-cluster properties are investigated both through the multistate empirical valence bond potential and a polarizable model. Equilibrium properties of the ion-water clusters H+(H2O)100, Na+(H2O)100, Na+(H2O)20, and Cl-(H2O)17 in the temperature region 100-450 K are explored using a hybrid parallel basin-hopping and tempering algorithm. The effect of the solid-liquid phase transition in both caloric curves and structural distribution functions is investigated. It is found that sodium and chloride ions largely reside on the surface of water clusters below the cluster melting temperature but are solvated into the interior of the cluster above the melting temperature, while the solvated proton was found to have significant propensity to reside on or near the surface in both the liquid- and solid-state clusters.

The properties of ion-water clusters. I. The protonated 21-water cluster.

The ab initio atom-centered density-matrix propagation approach and the multistate empirical valence bond method have been employed to study the structure, dynamics, and rovibrational spectrum of a hydrated proton in the "magic" 21 water cluster. In addition to the conclusion that the hydrated proton tends to reside on the surface of the cluster, with the lone pair on the protonated oxygen pointing "outwards," it is also found that dynamical effects play an important role in determining the vibrational properties of such clusters. This result is used to analyze and complement recent experimental and theoretical studies.

Anomalously soft dynamics of water in a nanotube: a revelation of nanoscale confinement.

Quasi-one-dimensional water encapsulated inside single-walled carbon nanotubes, here referred to as nanotube water, was studied by neutron scattering. The results reveal an anomalously soft dynamics characterized by pliable hydrogen bonds, anharmonic intermolecular potentials, and large-amplitude motions in nanotube water. Molecular dynamics simulations consistently describe the observed phenomena and propose the structure of nanotube water, which comprises a square-ice sheet wrapped into a cylinder inside the carbon nanotube and interior molecules in a chainlike configuration.

Effective force fields for condensed phase systems from ab initio molecular dynamics simulation: a new method for force-matching.

A novel least-squares fitting approach is presented to obtain classical force fields from trajectory and force databases produced by ab initio (e.g., Car-Parrinello) molecular dynamics (MD) simulations. The method was applied to derive effective nonpolarizable three-site force fields for liquid water at ambient conditions from Car-Parrinello MD simulations in the Becke-Lee-Yang-Parr approximation to the electronic density functional theory. The force-matching procedure includes a fit of short-ranged nonbonded forces, bonded forces, and atomic partial charges. The various parameterizations of the water force field differ by an enforced smooth cut-off applied to the short-ranged interaction term. These were obtained by fitting to the trajectory and force data produced by Car-Parrinello MD simulations of systems of 32 and 64 H(2)O molecules. The new water force fields were developed assuming both flexible or rigid molecular geometry. The simulated structural and self-diffusion properties of liquid water using the fitted force fields are in close agreement with those observed in the underlying Car-Parrinello MD simulations. The resulting empirical models compare to experiment much better than many conventional simple point charge (SPC) models. The fitted potential is also shown to combine well with more sophisticated intramolecular potentials. Importantly, the computational cost of the new models is comparable to that for SPC-like potentials.