On the surface of water ice, a quasi-liquid layer (QLL) has been extensively reported at temperatures below its bulk melting point at 273 K. Approaching the bulk melting temperature from below, the thickness of the QLL is known to increase. To elucidate the precise temperature variation of the QLL, and its nature, we investigate the surface melting of hexagonal ice by combining noncontact, surface-specific vibrational sum frequency generation (SFG) spectroscopy and spectra calculated from molecular dynamics simulations. Using SFG, we probe the outermost water layers of distinct single crystalline ice faces at different temperatures. For the basal face, a stepwise, sudden weakening of the hydrogen-bonded structure of the outermost water layers occurs at 257 K. The spectral calculations from the molecular dynamics simulations reproduce the experimental findings; this allows us to interpret our experimental findings in terms of a stepwise change from one to two molten bilayers at the transition temperature.

We introduce a simple "patchy particle" model to study the thermodynamics and dynamics of self-assembly of homomeric protein complexes. Our calculations allow us to rationalize recent results for dihedral complexes. Namely, why evolution of such complexes naturally takes the system into a region of interaction space where (i) the evolutionarily newer interactions are weaker, (ii) subcomplexes involving the stronger interactions are observed to be thermodynamically stable on destabilization of the protein-protein interactions and (iii) the self-assembly dynamics are hierarchical with these same subcomplexes acting as kinetic intermediates.

We discuss recent advances of the VOTCA package for systematic coarse-graining. Two methods have been implemented, namely the downhill simplex optimization and the relative entropy minimization. We illustrate the new methods by coarse-graining SPC/E bulk water and more complex water-methanol mixture systems. The CG potentials obtained from both methods are then evaluated by comparing the pair distributions from the coarse-grained to the reference atomistic simulations. In addition to the newly implemented methods, we have also added a parallel analysis framework to improve the computational efficiency of the coarse-graining process.

Structure-based coarse-graining relies on matching the pair correlation functions of a reference (atomistic) and a coarse-grained system. As such, it is designed for systems with uniform density distributions. Here, we demonstrate how it can be generalized for inhomogeneous systems by coarse-graining slabs of liquid water and methanol in vacuum, as well as a single benzene molecule at the water-vacuum interface. Our conclusion is that coarse-graining performed in inhomogeneous systems improves thermodynamic properties and the structure of interfaces without significant alterations to the local structure of the bulk liquid.

Biomineralization is the intricate, biomedically highly relevant process by which living organisms deposit minerals on biological matrices to stiffen tissues and build skeletal structures and shells. Rapaport and coworkers ( J. Am. Chem. Soc. 2000 , 122 , 12523 ; Adv. Funct. Mater. 2008 , 18 , 2889 ; Acta Biomater. 2012 , 8 , 2466 ) have designed a class of self-assembling amphiphilic peptides that are capable of forming hydrogels and attracting ions from the environment, generating structures akin to the extracellular matrix and promoting bone regeneration. The air-water interface serves both in experiment and in simulations as a model hydrophobic surface to mimic the cell's organic-aqueous interface and to investigate the organization of the peptide matrix into ordered β-pleated monolayers and the subsequent onset of biomineral formation. To obtain insight into the underlying molecular mechanism, we have used molecular dynamics simulations to study the effect of peptide sequence on aggregate stability and ion-peptide interactions. We find-in excellent agreement with experimental observations-that the nature of the peptide termini (proline vs phenylalanine) affect the aggregate order, while the nature of the acidic side chains (aspartic vs glutamic acid) affect the aggregate's stability in the presence of ions. These simulations provide valuable microscopic insight into the way ions and peptide templates mutually affect each other during the early stages of biomineralization preceding nucleation.