A number of periodic lattices have historically been used to represent ice-1h in computer simulations. These vary in size, shape, and method of generation, and while they have served their intended purposes, their properties have rarely been documented in detail and their intercompatibility is unknown. We develop a method for generating sets of internally consistent lattices and apply it to determine eight unit cells containing from 96 to 768 water molecules in both near-cubic and slab arrangements. It can easily be applied to generate additional (larger) cells or representations of specific crystal faces. Each unit cell in this set has zero net dipole moment and minimal net quadrupole moment and is optimized using four different criteria to measure the randomness of the hydrogen bonding; if required, these criteria can easily be modified to suit the intended application and alternate sets thus generated. We find that Cota and Hoover’s much used constraint for selecting unit cells with zero dipole moment is too restrictive, not permitting a fully random hydrogen-bonding network; also, unit-cell generation methods based on potential-energy minimization are found to prefer unrepresentative, highly ordered structures.
The version in the Kent Academic Repository may differ from the final published version. Users are advised to check http://kar.kent.ac.uk for the status of the paper. Users should always cite the published version of record.
The temperature dependence of the internal dynamics of an isolated protein, bovine pancreatic trypsin inhibitor, is examined using normal mode analysis and molecular dynamics (MD) simulation. It is found that the protein exhibits marked anharmonic dynamics at temperatures of approximately 100-120 K, as evidenced by departure of the MD-derived average mean square displacement from that of the harmonic model. This activation of anharmonic dynamics is at lower temperatures than previously detected in proteins and is found in the absence of solvent molecules. The simulation data are also used to investigate neutron scattering properties. The effects are determined of instrumental energy resolution and of approximations commonly used to extract mean square displacement data from elastic scattering experiments. Both the presence of a distribution of mean square displacements in the protein and the use of the Gaussian approximation to the dynamic structure factor lead to quantified underestimation of the mean square displacement obtained.
The structure and dynamics of the {0001} (basal), {101̄0} (prism), {202̄1}, and {21̄1̄0} ice Ih/water interfaces have been investigated using molecular dynamics and the flexible CF1 central force model of water. The translational order profile, the average density profile, and the diffusion profile have been calculated for all four interfaces as a function of distance normal to the interface. Dynamical molecular trajectories have been used to explore the loss of translational order from within the crystal region, through the interface, and into the liquid region. The thickness of the interfaces has been determined from each order parameter and compared with results from rigid models of water and experiment. The high index faces have thinner interfacial regions than the basal and prism interfaces. All interfacial regions contain molecules that are neither ice-like nor water-like.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.