BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

David Chandler

David Chandler appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechcond-mat.softphysics.chem-phcond-mat.mtrl-sciBiological Physicscond-mat.dis-nnmath.COastro-ph.IMBiomoleculesComputer Visioncond-mat.mes-halleess.IVMachine Learningphysics.comp-phphysics.hist-ph

Catalog footprint

What is connected

44works
15topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

44 published item(s)

preprint2022arXiv

Emu: A Case Study for TDI-like Imaging for Infrared Observation from Space

A wide-field zenith-looking telescope operating in a mode similar to Time-Delay-Integration (TDI) or drift scan imaging can perform an infrared sky survey without active pointing control but it requires a high-speed, low-noise infrared detector. Operating from a hosted payload platform on the International Space Station (ISS), the Emu space telescope employs the paradigm-changing properties of the Leonardo SAPHIRA electron avalanche photodiode array to provide powerful new observations of cool stars at the critical water absorption wavelength (1.4 $μ$m) largely inaccessible to ground-based telescopes due to the Earth's own atmosphere. Cool stars, especially those of spectral-type M, are important probes across contemporary astrophysics, from the formation history of the Galaxy to the formation of rocky exoplanets. Main sequence M-dwarf stars are the most abundant stars in the Galaxy and evolved M-giant stars are some of the most distant stars that can be individually observed. The Emu sky survey will deliver critical stellar properties of these cool stars by inferring oxygen abundances via measurement of the water absorption band strength at 1.4 $μ$m. Here we present the TDI-like imaging capability of Emu mission, its science objectives, instrument details and simulation results.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Very High Resolution Land Cover Mapping of Urban Areas at Global Scale with Convolutional Neural Networks

This paper describes a methodology to produce a 7-classes land cover map of urban areas from very high resolution images and limited noisy labeled data. The objective is to make a segmentation map of a large area (a french department) with the following classes: asphalt, bare soil, building, grassland, mineral material (permeable artificialized areas), forest and water from 20cm aerial images and Digital Height Model. We created a training dataset on a few areas of interest aggregating databases, semi-automatic classification, and manual annotation to get a complete ground truth in each class. A comparative study of different encoder-decoder architectures (U-Net, U-Net with Resnet encoders, Deeplab v3+) is presented with different loss functions. The final product is a highly valuable land cover map computed from model predictions stitched together, binarized, and refined before vectorization.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Pre-transition effects mediate forces of assembly between transmembrane proteins

We present a mechanism for a generic and powerful force of assembly and mobility for transmembrane proteins in lipid bilayers. This force is a pre-transition (or pre-melting) effect for the first-order transition between ordered and disordered phases in the host membrane. Using large scale molecular simulation, we show that a protein with hydrophobic thickness equal to that of the disordered phase embedded in an ordered bilayer stabilizes a microscopic order-disorder interface, and the stiffness of that interface is finite. When two such proteins approach each other, they assemble because assembly reduces the net interfacial free energy. In analogy with the hydrophobic effect, we refer to this phenomenon as the "orderphobic effect". The effect is mediated by proximity to the order-disorder phase transition and the size and hydrophobic mismatch of the protein. The strength and range of forces arising from the orderphobic effect are significantly larger than those that could arise from membrane elasticity for the membranes we examine.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The order-disorder transition in model lipid bilayers is a first-order hexatic to liquid phase transition

We characterize the order-disorder transition in a model lipid bilayer using molecular dynamics simulations. We find that the ordered phase is hexatic. In particular, in-plane structures possess a finite concentration of 5-7 disclination pairs that diffuse throughout the plane of the bilayer, and further, in-plane structures exhibit long-range orientational order and short-range translational order. In contrast, the disordered phase is liquid. The transition between the two phases is first order. Specifically, it exhibits hysteresis, and coexistence exhibits an interface with capillary scaling. The location of the interface and its spatial fluctuations are analyzed with a spatial field constructed from a rotational-invariant for local 6-fold orientational order. As a result of finite interfacial tension, there necessarily exist associated forces of assembly between membrane-bound solutes that pre-melt the ordered phase.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Time scales of supercooled water and implications for reversible polyamorphism

Deeply supercooled water exhibits complex dynamics with large density fluctuations, ice coarsening and characteristic time scales extending from picoseconds to milliseconds. Here, we discuss implications of these time scales as they pertain to two-phase coexistence and to molecular simulations of supercooled water. Specifically, we argue that it is possible to discount liquid-liquid criticality because the time scales imply that correlation lengths for such behavior would be bounded by no more than a few nanometers. Similarly, it is possible to discount two-liquid coexistence because the time scales imply a bounded interfacial free energy that cannot grow in proportion to a macroscopic surface area. From time scales alone, therefore, we see that coexisting domains of differing density in supercooled water can be no more than nano-scale transient fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Water exchange at a hydrated platinum electrode is rare and collective

We use molecular dynamics simulations to study the exchange kinetics of water molecules at a model metal electrode surface -- exchange between water molecules in the bulk liquid and water molecules bound to the metal. This process is a rare event, with a mean residence time of a bound water of about 40 ns for the model we consider. With analysis borrowed from the techniques of rare-event sampling, we show how this exchange or desorption is controlled by (1) reorganization of the hydrogen bond network within the adlayer of bound water molecules, and by (2) interfacial density fluctuations of the bulk liquid adjacent to the adlayer. We define collective coordinates that describe the desorption mechanism. Spatial and temporal correlations associated with a single event extend over nanometers and tens of picoseconds.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Comment on "Spontaneous liquid-liquid phase separation of water" by T. Yagasaki, M. Matsumoto and H. Tanaka, Phys. Rev. E 89, 020301 (2014)

Yagasaki et al. present results from a molecular dynamics trajectory illustrating coarsening of ice, which they interpret as evidence of transient coexistence between two distinct supercooled phases of liquid water. We point out that neither two distinct liquids nor criticality are demonstrated in this simulation study. Instead, the illustrated trajectory is consistent with coarsening behaviors analyzed and predicted in earlier works by others.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Excitations are localized and relaxation is hierarchical in glass-forming liquids

For several atomistic models of glass formers, at conditions below their glassy dynamics onset temperatures, ${T_\mathrm{o}}$, we use importance sampling of trajectory space to study the structure, statistics and dynamics of excitations responsible for structural relaxation. Excitations are detected in terms of persistent particle displacements of length $a$. At supercooled conditions, for $a$ of the order of or smaller than a particle diameter, we find that excitations are associated with correlated particle motions that are sparse and localized, occupying a volume with an average radius that is temperature independent and no larger than a few particle diameters. We show that the statistics and dynamics of these excitations are facilitated and hierarchical. Excitation energy scales grow logarithmically with $a$. Excitations at one point in space facilitate the birth and death of excitations at neighboring locations, and space-time excitation structures are microcosms of heterogeneous dynamics at larger scales. This nature of dynamics becomes increasingly dominant as temperature $T$ is lowered. We show that slowing of dynamics upon decreasing temperature below $T_\mathrm{o}$ is the result of a decreasing concentration of excitations and concomitant growing hierarchical length scales, and further that the structural relaxation time $τ$ follows the parabolic law, $\log(τ/ τ_\mathrm{o}) = J^2(1/T - 1/T_\mathrm{o})^2$, for $T<T_\mathrm{o}$, where $J$, $τ_\mathrm{o}$ and $T_\mathrm{o}$ can be predicted quantitatively from dynamics at short time scales. Particle motion is facilitated and directional, and we show this becomes more apparent with decreasing $T$. We show that stringlike motion is a natural consequence of facilitated, hierarchical dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Illusions of phase coexistence: Comments on "Metastable liquid-liquid transition ..." by J. C. Palmer et al., Nature 510, 385 (2014)

The recent paper cited above claims that a molecular simulation of one specific model of supercooled water establishes a stable interface separating two metastable liquid phases, which would imply the existence of metastable two-liquid criticality for that model. Here, we note that this claim conflicts with fundamental principles and with earlier work published in the \textit{Journal of Chemical Physics}, and we show that the claim is unjustified by the data put forward to support the conclusion. Other technical problems are also noted.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Premelting, fluctuations and coarse-graining of water-ice interfaces

Using statistical field theory supplemented with molecular dynamics simulations, we consider premelting on the surface of ice as a generic consequence of broken hydrogen bonds at the boundary between the condensed and gaseous phases. A procedure for coarse-graining molecular configurations onto a continuous scalar order parameter field is discussed, which provides a convenient representation of the interface between locally crystal-like and locally liquid-like regions. A number of interfacial properties are straightforwardly evaluated using this procedure such as the average premelting thickness and surface tension. The temperature and system size dependence of the premelting layer thickness calculated in this way confirms the characteristic logarithmic growth expected for the scalar field theory that the system is mapped onto through coarse-graining, though remains finite due to long-ranged interactions. Finally, from explicit simulations the existence of a premelting layer is shown to be insensitive to bulk lattice geometry, exposed crystal face and curvature.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The electric double layer has a life of its own

Using molecular dynamics simulations with recently developed importance sampling methods, we show that the differential capacitance of a model ionic liquid based double-layer capacitor exhibits an anomalous dependence on the applied electrical potential. Such behavior is qualitatively incompatible with standard mean-field theories of the electrical double layer, but is consistent with observations made in experiment. The anomalous response results from structural changes induced in the interfacial region of the ionic liquid as it develops a charge density to screen the charge induced on the electrode surface. These structural changes are strongly influenced by the out-of-plane layering of the electrolyte and are multifaceted, including an abrupt local ordering of the ions adsorbed in the plane of the electrode surface, reorientation of molecular ions, and the spontaneous exchange of ions between different layers of the electrolyte close to the electrode surface. The local ordering exhibits signatures of a first-order phase transition, which would indicate a singular charge-density transition in a macroscopic limit.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The molecular structure of the interface between water and a hydrophobic substrate is liquid-vapor like

With molecular simulation for water and a tunable hydrophobic substrate, we apply the instantaneous interface construction [A. P. Willard and D. Chandler, J. Phys. Chem. B, 114, 1954 (2010)] to examine the similarity between a water-vapor interface and a water-hydrophobic surface interface. The intrinsic interface refers to molecular structure in terms of distances from the instantaneous interface. We show that attractive interactions between a hydrophobic surface and water affect capillary wave fluctuations of the instantaneous liquid interface, but these attractive interactions have essentially no effect on the intrinsic interface. Further, the intrinsic interface of liquid water and a hydrophobic substrate differs little from that of water and its vapor.The same is not true, we show, for an interface between water and a hydrophilic substrate. In that case, strong directional substrate-water interactions disrupt the liquid-vapor-like interfacial hydrogen bonding network.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Theory of amorphous ices

We derive a phase diagram for amorphous solids and liquid supercooled water and explain why the amorphous solids of water exist in several different forms. Application of large-deviation theory allows us to prepare such phases in computer simulations. Along with nonequilibrium transitions between the ergodic liquid and two distinct amorphous solids, we establish coexistence between these two amorphous solids. The phase diagram we predict includes a nonequilibrium triple point where two amorphous phases and the liquid coexist. While the amorphous solids are long-lived and slowly-aging glasses, their melting can lead quickly to the formation of crystalline ice. Further, melting of the higher density amorphous solid at low pressures takes place in steps, transitioning to the lower density glass before accessing a nonequilibrium liquid from which ice coarsens.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Using the s-Ensemble to Probe Glasses Formed by Cooling and Aging

From length scale distributions characterizing frozen amorphous domains, we relate the s-ensemble method with standard cooling and aging protocols for forming glass. We show that in a general class of models, where space-time scaling is in harmony with that of experiment, the domain size distributions obtained with the s-ensemble are identical to those obtained through cooling or aging, but the computational effort for applying the s-ensemble is generally many orders of magnitude smaller than that of straightforward numerical simulation of cooling or aging.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Calorimetric glass transition explained by hierarchical dynamic facilitation

The glass transition refers to the non-equilibrium process by which an equilibrium liquid is transformed to a non-equilibrium disordered solid, or vice versa. Associated response functions, such as heat capacities, are markedly different on cooling than on heating, and the response to melting a glass depends markedly on the cooling protocol by which the glass was formed. This paper shows how this irreversible behavior can be interpreted quantitatively in terms of an East-model picture of localized excitations (or soft spots) in which molecules can move with a specific direction, and from which excitations with the same directionality of motion can appear or disappear in adjacent regions. As a result of this facilitated dynamics, excitations become correlated in a hierarchical fashion. These correlations are manifested in the dynamic heterogeneity of the supercooled liquid phase. While equilibrium thermodynamics is virtually featureless, a non-equilibrium glass phase emerges when the model is driven out of equilibrium with a finite cooling rate. The correlation length of this emergent phase is large and increases with decreasing cooling rate. A spatially and temporally resolved fictive temperature encodes memory of its preparation. Parameters characterizing the model can be determined from reversible transport data, and with these parameters, predictions of the model agree well with irreversible differential scanning calorimetry.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Characterizing heterogeneous dynamics at hydrated electrode surfaces

In models of Pt 111 and Pt 100 surfaces in water, motions of molecules in the first hydration layer are spatially and temporally correlated. To interpret these collective motions, we apply quantitative measures of dynamic heterogeneity that are standard tools for considering glassy systems. Specifically, we carry out an analysis in terms of mobility fields and distributions of persistence times and exchange times. In so doing, we show that dynamics in these systems is facilitated by transient disorder in frustrated two-dimensional hydrogen bonding networks. The frustration is the result of unfavorable geometry imposed by strong metal-water bonding. The geometry depends upon the structure of the underlying metal surface. Dynamic heterogeneity of water on the Pt 111 surface is therefore qualitatively different than that for water on the Pt 100 surface. In both cases, statistics of this adlayer dynamic heterogeneity responds asymmetrically to applied voltage.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Charge fluctuations in nano-scale capacitors

The fluctuations of the charge on an electrode contain information on the microscopic correlations within the adjacent fluid and their effect on the electronic properties of the interface. We investigate these fluctuations using molecular dynamics simulations in a constant-potential ensemble with histogram reweighting techniques. This approach offers in particular an efficient, accurate and physically insightful route to the differential capacitance that is broadly applicable. We demonstrate these methods with three different capacitors: pure water between platinum electrodes, and a pure as well as a solvent-based organic electrolyte each between graphite electrodes. The total charge distributions with the pure solvent and solvent-based electrolytes are remarkably Gaussian, while in the pure ionic liquid the total charge distribution displays distinct non-Gaussian features, suggesting significant potential-driven changes in the organization of the interfacial fluid.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Comment on "Is there a breakdown of the Stokes-Einstein relation in Kinetically Constrained Models at low temperature?" by O. Blondel and C. Toninelli, arXiv:1307.1651

We consider the diffusion constant, D, of a probe particle coupled to the East model, extending previous numerical results for this model to encompass a total of twelve orders of magnitude in relaxation time, τ. Our considerations thus cover the full range of relaxation times relevant to experimental studies of reversible glass formers. Throughout this range, we show that the results of the model are consistent with a fractional Stokes-Einstein law, D \propto τ^{-ξ} with ξ \approx 0.77. The same data can also be fit by relating Dτ to the East model excitation concentration q, specifically, Dτ \sim 1/q^α with α \approx 1.6. In either case, breakdown of the Stokes-Einstein relation in this kinetically constrained model is unquestionable and large.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Putative Liquid-Liquid Transition is a Liquid-Solid Transition in Atomistic Models of Water, Part II

This paper extends our earlier studies of free energy functions of density and crystalline order parameters for models of supercooled water, which allows us to examine the possibility of two distinct metastable liquid phases [J. Chem. Phys. 135, 134503 (2011) and arXiv:1107.0337v2]. Low-temperature reversible free energy surfaces of several different atomistic models are computed: mW water, TIP4P/2005 water, SW silicon and ST2 water, the last of these comparing three different treatments of long-ranged forces. In each case, we show that there is one stable or metastable liquid phase, and there is an ice-like crystal phase. The time scales for crystallization in these systems far exceed those of structural relaxation in the supercooled metastable liquid. We show how this wide separation in time scales produces an illusion of a low-temperature liquid-liquid transition. The phenomenon suggesting metastability of two distinct liquid phases is actually coarsening of the ordered ice-like phase, which we elucidate using both analytical theory and computer simulation. For the latter, we describe robust methods for computing reversible free energy surfaces, and we consider effects of electrostatic boundary conditions. We show that sensible alterations of models and boundary conditions produce no qualitative changes in low-temperature phase behaviors of these systems, only marginal changes in equations of state. On the other hand, we show that altering sampling time scales can produce large and qualitative nonequilibrium effects. Recent reports of evidence of a liquid-liquid critical point in computer simulations of supercooled water are considered in this light.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Lectures on Molecular- and Nano-scale Fluctuations in Water

This manuscript is the written form of three lectures delivered by David Chandler at the International School of Physics "Enrico Fermi", Course CLXXVI - "Complex materials in physics and biology", in Varenna, Italy in July 2010. It describes the physical properties of water from a molecular perspective and how these properties are reflected in the behaviors of water as a solvent. Theory of hydrophobicity and solvation of ions are topics included in the discussion.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Phase diagram of supercooled water confined to hydrophilic nanopores

We present a phase diagram for water confined to cylindrical silica nanopores in terms of pressure, temperature and pore radius. The confining cylindrical wall is hydrophilic and disordered, which has a destabilizing effect on ordered water structure. The phase diagram for this class of systems is derived from general arguments, with parameters taken from experimental observations and computer simulations and with assumptions tested by computer simulation. Phase space divides into three regions: a single liquid, a crystal-like solid, and glass. For large pores, radii exceeding 1 nm, water exhibits liquid and crystal-like behaviors, with abrupt crossovers between these regimes. For small pore radii, crystal-like behavior is unstable and water remains amorphous for all non-zero temperatures. At low enough temperatures, these states are glasses. Several experimental results for supercooled water can be understood in terms of the phase diagram we present.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Solvation at Aqueous Metal Electrodes

We present a study of the solvation properties of model aqueous electrode interfaces. The exposed electrodes we study strongly bind water and have closed packed crystalline surfaces, which template an ordered water adlayer adjacent to the interface. We find that these ordered water structures facilitate collective responses in the presence of solutes that are correlated over large lengthscales and across long timescales. Specifically, we show that the liquid water adjacent to the ordered adlayers forms a soft, liquid-vapor-like interface with concomitant manifestations of hydrophobicity. Temporal defects in the adlayer configurations create a dynamic heterogeneity in the degree to which different regions of the interface attract hydrophobic species. The structure and heterogeneous dynamics of the adlayer defects depend upon the geometry of the underlying ordered metal surface. For both 100 and 111 surfaces, the dynamical heterogeneity relaxes on times longer than nanoseconds. Along with analyzing time scales associated with these effects, we highlight implications for electrolysis and the particular catalytic efficiency of platinum.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Water Evaporation: A Transition Path Sampling Study

We use transition path sampling to study evaporation in the SPC/E model of liquid water. Based on thousands of evaporation trajectories, we characterize the members of the transition state ensemble (TSE), which exhibit a liquid-vapor interface with predominantly negative mean curvature at the site of evaporation. We also find that after evaporation is complete, the distributions of translational and angular momenta of the evaporated water are Maxwellian with a temperature equal to that of the liquid. To characterize the evaporation trajectories in their entirety, we find that it suffices to project them onto just two coordinates: the distance of the evaporating molecule to the instantaneous liquid-vapor interface, and the velocity of the water along the average interface normal. In this projected space, we find that the TSE is well-captured by a simple model of ballistic escape from a deep potential well, with no additional barrier to evaporation beyond the cohesive strength of the liquid. Equivalently, they are consistent with a near-unity probability for a water molecule impinging upon a liquid droplet to condense. These results agree with previous simulations and with some, but not all, recent experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Extended surfaces modulate and can catalyze hydrophobic effects

Interfaces are a most common motif in complex systems. To understand how the presence of interfaces affect hydrophobic phenomena, we use molecular simulations and theory to study hydration of solutes at interfaces. The solutes range in size from sub-nanometer to a few nanometers. The interfaces are self-assembled monolayers with a range of chemistries, from hydrophilic to hydrophobic. We show that the driving force for assembly in the vicinity of a hydrophobic surface is weaker than that in bulk water, and decreases with increasing temperature, in contrast to that in the bulk. We explain these distinct features in terms of an interplay between interfacial fluctuations and excluded volume effects---the physics encoded in Lum-Chandler-Weeks theory [J. Phys. Chem. B 103, 4570--4577 (1999)]. Our results suggest a catalytic role for hydrophobic interfaces in the unfolding of proteins, for example, in the interior of chaperonins and in amyloid formation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Molecular explanation for why talc surfaces can be both hydrophilic and hydrophobic

While individual water molecules adsorb strongly on a talc surface (hydrophilic behavior), a droplet of water beads up on the same surface (hydrophobic behavior). To rationalize this dichotomy, we investigate the influence of the microscopic structure of the surface and the strength of adhesive (surface-water) interactions on surface hydrophobicity. We show that at low relative humidity, the competition between adhesion and the favorable entropy of being in the vapor phase determines the surface coverage. However, at saturation, it is the competition between adhesion and cohesion (water-water interactions) that determines surface hydrophobicity. The adhesive interactions in talc are strong enough to overcome the unfavorable entropy, and water adsorbs strongly on talc surfaces. However, they are too weak to overcome the cohesive interactions, and water thus beads up on talc surfaces. Surprisingly, even (talc-like) surfaces that are highly adhesive, do not fully wet at saturation. Instead, a water droplet forms on top of a strongly adsorbed monolayer of water. Our results imply that the interior of hydrophobic zeolites suspended in water may contain adsorbed water molecules at pressures much smaller than the intrusion pressure.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Preparation and relaxation of very stable glassy states of a simulated liquid

We prepare metastable glassy states in a model glass-former made of Lennard-Jones particles by sampling biased ensembles of trajectories with low dynamical activity. These trajectories form an inactive dynamical phase whose `fast' vibrational degrees of freedom are maintained at thermal equilibrium by contact with a heat bath, while the `slow' structural degrees of freedom are located in deep valleys of the energy landscape. We examine the relaxation to equilibrium and the vibrational properties of these metastable states. The glassy states we prepare by our trajectory sampling method are very stable to thermal fluctuations and also more mechanically rigid than low-temperature equilibrated configurations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Quantifying density fluctuations in volumes of all shapes and sizes using indirect umbrella sampling

Water density fluctuations are an important statistical mechanical observable that is related to many-body correlations, as well as hydrophobic hydration and interactions. Local water density fluctuations at a solid-water surface have also been proposed as a measure of its hydrophobicity. These fluctuations can be quantified by calculating the probability, $P_v(N)$, of observing $N$ waters in a probe volume of interest $v$. When $v$ is large, calculating $P_v(N)$ using molecular dynamics simulations is challenging, as the probability of observing very few waters is exponentially small, and the standard procedure for overcoming this problem (umbrella sampling in $N$) leads to undesirable impulsive forces. Patel et al. [J. Phys. Chem. B, 114, 1632 (2010)] have recently developed an indirect umbrella sampling (INDUS) method, that samples a coarse-grained particle number to obtain $P_v(N)$ in cuboidal volumes. Here, we present and demonstrate an extension of that approach to other basic shapes, like spheres and cylinders, as well as to collections of such volumes. We further describe the implementation of INDUS in the NPT ensemble and calculate $P_v(N)$ distributions over a broad range of pressures. Our method may be of particular interest in characterizing the hydrophobicity of interfaces of proteins, nanotubes and related systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Sitting at the edge: How biomolecules use hydrophobicity to tune their interactions and function

Water near hydrophobic surfaces is like that at a liquid-vapor interface, where fluctuations in water density are substantially enhanced compared to that in bulk water. Here we use molecular simulations with specialized sampling techniques to show that water density fluctuations are similarly enhanced, even near hydrophobic surfaces of complex biomolecules, situating them at the edge of a dewetting transition. Consequently, water near these surfaces is sensitive to subtle changes in surface conformation, topology, and chemistry, any of which can tip the balance towards or away from the wet state, and thus significantly alter biomolecular interactions and function. Our work also resolves the long-standing puzzle of why some biological surfaces dewet and other seemingly similar surfaces do not.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Putative Liquid-Liquid Transition is a Liquid-Solid Transition in Atomistic Models of Water

We use numerical simulation to examine the possibility of a reversible liquid-liquid transition in supercooled water and related systems. In particular, for two atomistic models of water, we have computed free energies as functions of multiple order parameters, where one is density and another distinguishes crystal from liquid. For a range of temperatures and pressures, separate free energy basins for liquid and crystal are found, conditions of phase coexistence between these phases are demonstrated, and time scales for equilibration are determined. We find that at no range of temperatures and pressures is there more than a single liquid basin, even at conditions where amorphous behavior is unstable with respect to the crystal. We find a similar result for a related model of silicon. This result excludes the possibility of the proposed liquid-liquid critical point for the models we have studied. Further, we argue that behaviors others have attributed to a liquid-liquid transition in water and related systems are in fact reflections of transitions between liquid and crystal.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

An improved coarse-grained model of solvation and the hydrophobic effect

We present a coarse-grained lattice model of solvation thermodynamics and the hydrophobic effect that implements the ideas of Lum-Chandler-Weeks (LCW) theory [J. Phys. Chem. B 103, 4570 (1999)] and improves upon previous lattice models based on it. Through comparison with molecular simulation, we show that our model captures the length-scale and curvature dependence of solvation free energies with near-quantitative accuracy and two to three orders of magnitude less computational effort, and further, correctly describes the large but rare solvent fluctuations that are involved in dewetting, vapor tube formation and hydrophobic assembly. Our model is intermediate in detail and complexity between implicit-solvent models and explicit-water simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Corresponding States of Structural Glass Formers. II

The earlier paper of this same title demonstrated a collapse of relaxation data of fragile supercooled glass forming liquids [\textit{J. Phys. Chem. B} \textbf{113}, 5563-5567 (2009)]. For temperature $T$ below that of the onset to supercooled behavior, $T_{\mathrm{o}}$, the logarithm of structural relaxation time, $\log τ$, is given by the parabolic form $\log (τ/τ_{\mathrm{o}}) \, =\, J^2(1/T - 1/T_{\mathrm{o}})^{2}$, where $J$ and $τ_{\mathrm{o}}$ are temperature independent. This paper presents further applications of this formula. In particular, it is shown that the effects of attractive forces in numerical simulation of glass forming liquids can be logically organized in terms of $J$ and $T_{\mathrm{o}}$. Further, analysis of experimental data for several systems suggests that $J$ and $T_{\mathrm{o}}$ are material properties. In contrast, values of similar parameters for other fitting formulas are shown to depend not only upon the material but also upon the range of data used in fitting these formulas. Expressions demonstrated to fail in this way include the Vogel-Fulcher-Tammann formula, a double-exponential formula, and a fractional exponential formula.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Finite-temperature critical point of a glass transition

We generalize the simplest kinetically constrained model of a glass-forming liquid by softening kinetic constraints, allowing them to be violated with a small finite rate. We demonstrate that this model supports a first-order dynamical (space-time) phase transition, similar to those observed with hard constraints. In addition, we find that the first-order phase boundary in this softened model ends in a finite-temperature dynamical critical point, which we expect to be present in natural systems. We discuss links between this critical point and quantum phase transitions, showing that dynamical phase transitions in $d$ dimensions map to quantum transitions in the same dimension, and hence to classical thermodynamic phase transitions in $d+1$ dimensions. We make these links explicit through exact mappings between master operators, transfer matrices, and Hamiltonians for quantum spin chains.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Dynamics on the Way to Forming Glass: Bubbles in Space-time

We review a theoretical perspective of the dynamics of glass forming liquids and the glass transition. It is a perspective we have developed with our collaborators during this decade. It is based upon the structure of trajectory space. This structure emerges from spatial correlations of dynamics that appear in disordered systems as they approach non-ergodic or jammed states. It is characterized in terms of dynamical heterogeneity, facilitation and excitation lines. These features are associated with a newly discovered class of non-equilibrium phase transitions. Equilibrium properties have little if anything to do with it. The broken symmetries of these transitions are obscure or absent in spatial structures, but they are vivid in space-time (i.e., trajectory space). In our view, the glass transition is an example of this class of transitions. The basic ideas and principles we review were originally developed through the analysis of idealized and abstract models. Nevertheless, the central ideas are easily illustrated with reference to molecular dynamics of more realistic atomistic models, and we use that illustrative approach here.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Fluctuations of water near extended hydrophobic and hydrophilic surfaces

We use molecular dynamics simulations of the SPC-E model of liquid water to derive probability distributions for water density fluctuations in probe volumes of different shapes and sizes, both in the bulk as well as near hydrophobic and hydrophilic surfaces. To obtain our results, we introduce a biased sampling of coarse-grained densities, which in turn biases the actual solvent density. The technique is easily combined with molecular dynamics integration algorithms. Our principal result is that the probability for density fluctuations of water near a hydrophobic surface, with or without surface-water attractions, is akin to density fluctuations at the water-vapor interface. Specifically, the probability of density depletion near the surface is significantly larger than that in bulk. In contrast, we find that the statistics of water density fluctuations near a model hydrophilic surface are similar to that in the bulk.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Instantaneous Liquid Interfaces

We describe and illustrate a simple procedure for identifying a liquid interface from atomic coordinates. In particular, a coarse grained density field is constructed, and the interface is defined as a constant density surface for this coarse grained field. In applications to a molecular dynamics simulation of liquid water, it is shown that this procedure provides instructive and useful pictures of liquid-vapor interfaces and of liquid-protein interfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Coarse Grained Modeling of The Interface Between Water and Heterogeneous Surfaces

Using coarse grained models we investigate the behavior of water adjacent to an extended hydrophobic surface peppered with various fractions of hydrophilic patches of different sizes. We study the spatial dependence of the mean interface height, the solvent density fluctuations related to drying the patchy substrate, and the spatial dependence of interfacial fluctuations. We find that adding small uniform attractive interactions between the substrate and solvent cause the mean position of the interface to be very close to the substrate. Nevertheless, the interfacial fluctuations are large and spatially heterogeneous in response to the underlying patchy substrate. We discuss the implications of these findings to the assembly of heterogeneous surfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Corresponding States of Structural Glass Formers

The variation with respect to temperature T of transport properties of 58 fragile structural glass forming liquids (68 data sets in total) are analyzed and shown to exhibit a remarkable degree of universality. In particular, super-Arrhenius behaviors of all super-cooled liquids appear to collapse to one parabola for which there is no singular behavior at any finite temperature. This behavior is bounded by an onset temperature To above which liquid transport has a much weaker temperature dependence. A similar collapse is also demonstrated, over the smaller available range, for existing numerical simulation data.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Water at an electrochemical interface - a simulation study

The results of molecular dynamics simulations of the properties of water in an aqueous ionic solution close to an interface with a model metallic electrode are described. In the simulations the electrode behaves as an ideally polarizable hydrophilic metal, supporting image charge interactions with charged species, and it is maintained at a constant electrical potential with respect to the solution so that the model is a textbook representation of an electrochemical interface through which no current is passing. We show how water is strongly attracted to and ordered at the electrode surface. This ordering is different to the structure that might be imagined from continuum models of electrode interfaces. Further, this ordering significantly affects the probability of ions reaching the surface. We describe the concomitant motion and configurations of the water and ions as functions of the electrode potential, and we analyze the length scales over which ionic atmospheres fluctuate. The statistics of these fluctuations depend upon surface structure and ionic strength. The fluctuations are large, sufficiently so that the mean ionic atmosphere is a poor descriptor of the aqueous environment near a metal surface. The importance of this finding for a description of electrochemical reactions is examined by calculating, directly from the simulation, Marcus free energy profiles for transfer of charge between the electrode and a redox species in the solution and comparing the results with the predictions of continuum theories. Significant departures from the electrochemical textbook descriptions of the phenomenon are found and their physical origins are characterized from the atomistic perspective of the simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

De-coupling of Exchange and Persistence Times in Atomistic Models of Glass Formers

With molecular dynamics simulations of a fluid mixture of classical particles interacting with pair-wise additive Weeks-Chandler-Andersen potentials, we consider the time series of particle displacements and thereby determine distributions for local persistence times and local exchange times. These basic characterizations of glassy dynamics are studied over a range of super-cooled conditions and shown to have behaviors, most notably de-coupling, similar to those found in kinetically constrained lattice models of structural glasses. Implications are noted.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Solvent coarse-graining and the string method applied to the hydrophobic collapse of a hydrated chain

Using computer simulations of over 100,000 atoms, the mechanism for the hydrophobic collapse of an idealized hydrated chain is obtained. This is done by coarse-graining the atomistic water molecule positions over 129,000 collective variables that represent the water density field and then using the string method in these variables to compute the minimum free energy pathway (MFEP) for the collapsing chain. The dynamical relevance of the MFEP (i.e. its coincidence with the mechanism of collapse) is validated a posteriori using conventional molecular dynamics trajectories. Analysis of the MFEP provides atomistic confirmation for the mechanism of hydrophobic collapse proposed by ten Wolde and Chandler. In particular, it is shown that lengthscale-dependent hydrophobic dewetting is the rate-limiting step in the hydrophobic collapse of the considered chain.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Model of a fluid at small and large length scales and the hydrophobic effect

We present a statistical field theory to describe large length scale effects induced by solutes in a cold and otherwise placid liquid. The theory divides space into a cubic grid of cells. The side length of each cell is of the order of the bulk correlation length of the bulk liquid. Large length scale states of the cells are specified with an Ising variable. Finer length scale effects are described with a Gaussian field, with mean and variance affected by both the large length scale field and by the constraints imposed by solutes. In the absence of solutes and corresponding constraints, integration over the Gaussian field yields an effective lattice gas Hamiltonian for the large length scale field. In the presence of solutes, the integration adds additional terms to this Hamiltonian. We identify these terms analytically. They can provoke large length scale effects, such as the formation of interfaces and depletion layers. We apply our theory to compute the reversible work to form a bubble in liquid water, as a function of the bubble radius. Comparison with molecular simulation results for the same function indicates that the theory is reasonably accurate. Importantly, simulating the large length scale field involves binary arithmetic only. It thus provides a computationally convenient scheme to incorporate explicit solvent dynamics and structure in simulation studies of large molecular assemblies.

0 citations0 reviewsNo code linkedOpen access