Christian J. Burnham

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.

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.

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.

Dynamics of hydrogen guests in ice XVII nanopores

The present high-resolution inelastic neutron scattering experiment on ice XVII, containing molecular hydrogen with different ortho/para ratio, allows to assign the H

Free-Energy Calculations of the Intercage Hopping Barriers of Hydrogen Molecules in Clathrate Hydrates

_2

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

motion spectral bands to rotational and center-of-mass translational transitions of either {\it para}- or {\it ortho}-H

Vibrational Modes of Hydrogen Hydrates: A First-Principles Molecular Dynamics and Raman Spectra Study

_2

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

. Due to its structure, ice XVII confines H

Ice-Amorphization of Supercooled Water Nanodroplets in No Man’s Land

_2

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

molecules to move in spiral channels of molecular size. Reported data demonstrate that H