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 T. Limmer

David T. Limmer 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-mechphysics.chem-phcond-mat.softcond-mat.mtrl-scicond-mat.mes-hallphysics.app-phcond-mat.dis-nnBiological PhysicsMachine Learningphysics.comp-phquant-ph

Catalog footprint

What is connected

31works
11topics
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

31 published item(s)

preprint2026arXiv

Atomic Alignment in PbS Nanocrystal Superlattices with Compact Inorganic Ligands via Reversible Oriented Attachment of Nanocrystals

Nanocrystals (NCs) serve as versatile building blocks for the creation of functional materials, with NC self-assembly offering opportunities to enable novel material properties. Here, we demonstrate that PbS NCs functionalized with strongly negatively charged metal chalcogenide complex (MCC) ligands, such as $Sn_2S_6^{4-}$ and $AsS_4^{3-}$, can self-assemble into all-inorganic superlattices with both long-range superlattice translational and atomic-lattice orientational order. Structural characterizations reveal that the NCs adopt unexpected edge-to-edge alignment, and numerical simulation clarifies that orientational order is thermodynamically stabilized by many-body ion correlations originating from the dense electrolyte. Furthermore, we show that the superlattices of $Sn_2S_6^{4-}$-functionalized PbS NCs can be fully disassembled back into the colloidal state, which is highly unusual for orientationally attached superlattices with atomic-lattice alignment. The reversible oriented attachment of NCs, enabling their dynamic assembly and disassembly into effectively single-crystalline superstructures, offers a pathway toward designing reconfigurable materials with adaptive and controllable electronic and optoelectronic properties.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Nanocrystal Geometry Governs Phase Transformation Pathways in Palladium Hydride

Pathways and structural dynamics of phase transformations impact performance of materials in energy and information storage technologies. Palladium hydride ($\mathrm{PdH}_x$) nanocrystals are an ideal model system for studying solute-induced phase transformations, where elastic energy from lattice mismatch between $α$-$\mathrm{PdH}_x$ and $β$-$\mathrm{PdH}_x$ phases is often considered a key to determining the transformation pathways. $α/β$-$\mathrm{PdH}_x$ interfacial elastic energy is affected by the confined geometry of a nanocrystal. However, how nanocrystal geometry influences phase transformation pathways is largely unknown. Using in situ liquid phase transmission electron microscopy, we directly visualize hydrogenation in Pd nanocrystals with two geometries -- a nanocube and a hexagonal nanoplate. Both follow similar sequences of an initially curved nucleus, interface flattening, and reverse-stage nucleation; however, their evolving $α/β$-$\mathrm{PdH}_x$ interfaces exhibit geometry-dependent crystallographic alignments. In nanocubes, $\{100\}$-aligned configurations conform to static elastic energy ordering, representing a pathway that maintains a local mechanical equilibrium, whereas nanoplates display both $\{110\}$- and $\{211\}$-aligned interfaces. Theoretical simulations show that geometry determines the accessibility of alternative phase transformation pathways as the system is driven far from equilibrium during hydrogenation. These findings identify geometry as a fundamental parameter for directing phase transformation pathways, offering design principles for accessing atypical configurations and improving properties of intercalation-based devices.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Nonlocal screening dictates the radiative lifetimes of excitations in lead halide perovskites

We use path integral molecular dynamics simulations and theory to elucidate the interactions between charge carriers, as mediated by a lead halide perovskite lattice. We find that the charge-lattice coupling of MAPbI$_3$ results in a repulsive interaction between electrons and holes at intermediate distances. The effective interaction is understood using a Gaussian field theory, whereby the underlying soft, polar lattice contributes a nonlocal screening between quasiparticles. Path integral calculations of this nonlocal screening model are used to rationalize the small exciton binding energy and low radiative recombination rate observed experimentally and are compared to traditional Wannier-Mott and Fröhlich models, which fail to do so. These results clarify the origin of the high power conversion efficiencies in lead halide perovskites. Emergent repulsive electron-hole interactions provide a design principle for optimizing soft, polar semiconductors.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Peptide Isomerization is Suppressed at the Air-Water Interface

We use molecular dynamics simulations to study the thermodynamics and kinetics of alanine dipeptide isomerization at the air-water interface. Thermodynamically, we find an affinity of the dipeptide to the interface. This affinity arises from stablizing intramolecular interactions that become unshielded as the dipeptide is desolvated. Kinetically, we consider the rate of transitions between the alpha-L and beta conformations of alanine dipeptide and evaluate it as a continuous function of the distance from the interface using a recent extension of transition path sampling, TPS+U. The rate of isomerization at the Gibbs dividing surface is suppressed relative to the bulk by a factor of 3. Examination of the ensemble of transition states elucidates the role of solvent degrees of freedom in mediating favorable intramolecular interactions along the reaction pathway of isomerization. Near the air-water interface, water is less effective at mediating these intramolecular interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Direct evaluation of rare events in active matter from variational path sampling

Active matter represents a broad class of systems that evolve far from equilibrium due to the local injection of energy. Like their passive analogues, transformations between distinct metastable states in active matter proceed through rare fluctuations, however their detailed balance violating dynamics renders these events difficult to study. Here, we present a simulation method for evaluating the rate and mechanism of rare events in generic nonequilibrium systems and apply it to study the conformational changes of a passive solute in an active fluid. The method employs a variational optimization of a control force that renders the rare event a typical one, supplying an exact estimate of its rate as a ratio of path partition functions. Using this method we find that increasing activity in the active bath can enhance the rate of conformational switching of the passive solute in a manner consistent with recent bounds from stochastic thermodynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Reinforcement learning of rare diffusive dynamics

We present a method to probe rare molecular dynamics trajectories directly using reinforcement learning. We consider trajectories that are conditioned to transition between regions of configuration space in finite time, like those relevant in the study of reactive events, as well as trajectories exhibiting rare fluctuations of time-integrated quantities in the long time limit, like those relevant in the calculation of large deviation functions. In both cases, reinforcement learning techniques are used to optimize an added force that minimizes the Kullback-Leibler divergence between the conditioned trajectory ensemble and a driven one. Under the optimized added force, the system evolves the rare fluctuation as a typical one, affording a variational estimate of its likelihood in the original trajectory ensemble. Low variance gradients employing value functions are proposed to increase the convergence of the optimal force. The method we develop employing these gradients leads to efficient and accurate estimates of both the optimal force and the likelihood of the rare event for a variety of model systems.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Variational control forces for enhanced sampling of nonequilibrium molecular dynamics simulations

We introduce a variational algorithm to estimate the likelihood of a rare event within a nonequilibrium molecular dynamics simulation through the evaluation of an optimal control force. Optimization of a control force within a chosen basis is made possible by explicit forms for the gradients of a cost function in terms of the susceptibility of driven trajectories to changes in variational parameters. We consider probabilities of time-integrated dynamical observables as characterized by their large deviation functions, and find that in many cases the variational estimate is quantitatively accurate. Additionally, we provide expressions to exactly correct the variational estimate that can be evaluated directly. We benchmark this algorithm against the numerically exact solution of a model of a driven particle in a periodic potential, where the control force can be represented with a complete basis. We then demonstrate the utility of the algorithm in a model of repulsive particles on a line, which undergo a dynamical phase transition, resulting in singular changes to the form of the optimal control force. In both systems, we find fast convergence and are able to evaluate large deviation functions with significant increases in statistical efficiency over alternative Monte Carlo approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Variational design principles for nonequilibrium colloidal assembly

Using large deviation theory and principles of stochastic optimal control, we show that rare molecular dynamics trajectories conditioned on assembling a specific target structure encode a set of interactions and external forces that lead to enhanced stability of that structure. Such a relationship can be formulated into a variational principle, for which we have developed an associated optimization algorithm and have used it to determine optimal forces for targeted self-assembly within nonequilibrium steady-states. We illustrate this perspective on inverse design in a model of colloidal cluster assembly within linear shear flow. We find that colloidal clusters can be assembled with high yield using specific short-range interactions of tunable complexity. Shear decreases the yields of rigid clusters, while small values of shear increase the yields of nonrigid clusters. The enhancement or suppression of the yield due to shear is rationalized with a generalized linear response theory. By studying 21 unique clusters made of 6, 7 or 8 particles, we uncover basic design principles for targeted assembly out of equilibrium.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Entropy production fluctuations encode collective behavior in active matter

We derive a general lower bound on distributions of entropy production in interacting active matter systems. The bound is tight in the limit that interparticle correlations are small and short-ranged, which we explore in four canonical active matter models. In all models studied, the bound is weak where collective fluctuations result in long-ranged correlations, which subsequently links the locations of phase transitions to enhanced entropy production fluctuations. We develop a theory for the onset of enhanced fluctuations and relate it to specific phase transitions in active Brownian particles. We also derive optimal control forces that realize the dynamics necessary to tune dissipation and manipulate the system between phases. In so doing, we uncover a general relationship between entropy production and pattern formation in active matter, as well as ways of controlling it.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Field-dependent ionic conductivities from generalized fluctuation-dissipation relations

We derive a relationship for the electric field dependent ionic conductivity in terms of fluctuations of time integrated microscopic variables. We demonstrate this formalism with molecular dynamics simulations of solutions of differing ionic strength with implicit solvent conditions and molten salts. These calculations are aided by a novel nonequilibrium statistical reweighting scheme that allows for the conductivity to be computed as a continuous function of the applied field. In strong electrolytes, we find the fluctuations of the ionic current are Gaussian and subsequently the conductivity is constant with applied field. In weaker electrolytes and molten salts, we find the fluctuations of the ionic current are strongly non-Gaussian and the conductivity increases with applied field. This nonlinear behavior, known phenomenologically for dilute electrolytes as the Onsager-Wien effect, is general and results from the suppression of ionic correlations at large applied fields, as we elucidate through both dynamic and static correlations within nonequilibrium steady-states.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fundamental limitations on photoisomerization from thermodynamic resource theories

Small, out-of-equilibrium, and quantum systems defy simple thermodynamic expressions. Such systems are exemplified by molecular switches, which exchange heat with a bath. These molecules can photoisomerize, or change conformation, or switch, upon absorbing light. The photoisomerization probability depends on kinetic details that couple the molecule's energetics to its dissipation. Therefore, a simple, general, thermodynamic-style bound on the photoisomerization probability seems out of reach. We derive such a bound using a resource theory. The resource-theory framework is a set of mathematical tools, developed in quantum information theory, used to generalize thermodynamics to small and quantum settings. From this toolkit has been derived a generalization of the second law, the thermomajorization preorder. We use thermomajorization to upper-bound the photoisomerization probability. Then, we compare the bound with an equilibrium prediction and with a Lindbladian model. We identify a realistic parameter regime in which the Lindbladian evolution saturates the thermomajorization bound. We also quantify the energy coherence in the electronic degree of freedom, and we argue that this coherence cannot promote photoisomerization. This work illustrates how quantum-information-theoretic thermodynamics can elucidate complex quantum processes in nature, experiments, and synthetics.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Photoinduced phase separation in the lead halides is a polaronic effect

We present a perspective on recent observations of the photoinduced phase separation of halides in multi-component lead-halide perovskites. The spontaneous phase separation of an initial homogeneous solid solution under steady-state illumination conditions is found experimentally to be reversible, stochastic, weakly dependent on morphology, yet strongly dependent on composition and thermodynamic state. Regions enriched in a specific halide species that form upon phase separation are self-limiting in size, pinned to specific compositions, and grow in number in proportion to the steady-state carrier concentration until saturation. These empirical observations of robustness rule out explanations based on specific defect structures and point to the local modulation of an existing miscibility phase transition in the presence of excess charge carriers. A model for rationalizing existing observations based on the coupling between composition, strain and charge density fluctuations through the formation of polarons is reviewed.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Heat current fluctuations and anomalous transport in low dimensional carbon lattices

Molecular dynamics simulations and nonequilibrium importance sampling are used to study the heat transport of low dimensional carbon lattices. For both carbon nanotubes and graphene sheets heat transport is found to be anomalous, violating Fourier's law of conduction with a system size dependent thermal conductivity and concomitant nonlinear temperature profiles. For carbon nanotubes, the thermal conductivity is found to increase as the square root of the length of the nanotube, while for graphene sheets the thermal conductivity is found to increase as the logarithm of the length of the sheet. The particular length dependence and nonlinear temperature profiles place carbon lattices into a universality class with nonlinear lattice models, and suggest that heat transport through carbon nano-structures is better described by a Levy walk rather than simple diffusion.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Liquid-Like Interfaces Mediate Structural Phase Transitions in Lead Halide Perovskites

Microscopic pathways of structural phase transitions are difficult to probe because they occur over multiple, disparate time and length scales. Using $in$ $situ$ nanoscale cathodoluminescence microscopy, we visualize the thermally-driven transition to the perovskite phase in hundreds of non-perovskite phase nanowires, resolving the initial nanoscale nucleation and subsequent mesoscale growth and quantifying the activation energy for phase propagation. In combination with molecular dynamics computer simulations, we reveal that the transformation does not follow a simple martensitic mechanism, and proceeds via ion diffusion through a liquid-like interface between the two structures. While cations are disordered in this liquid-like region, the halide ions retain substantial spatial correlations. We find that the anisotropic crystal structure translates to faster nucleation of the perovskite phase at nanowire ends and faster growth along the long nanowire axis. These results represent a significant step towards manipulating structural phases at the nanoscale for designer materials properties.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Microscopic dynamics of charge separation at the aqueous electrochemical interface

We have used molecular simulation and methods of importance sampling to study the thermodynamics and kinetics of ionic charge separation at a liquid water-metal interface. We have considered this process using canonical examples of two different classes of ions: a simple alkali-halide pair, Na$^+$I$^-$, or classical ions, and the products of water autoionization, H$_3$O$^+$OH$^-$, or water ions. We find that for both ion classes, the microscopic mechanism of charge separation, including water's collective role in the process, is conserved between the bulk liquid and the electrode interface. Despite this, the thermodynamic and kinetic details of the process differ between these two environments in a way that depends on ion type. In the case of the classical ion pairs, a higher free energy barrier to charge separation and a smaller flux over that barrier at the interface, results in a rate of dissociation that is 40x slower relative to the bulk. For water ions, a slightly higher free energy barrier is offset by a higher flux over the barrier from longer lived hydrogen bonding patters at the interface, resulting in a rate of association that is similar both at and away from the interface. We find that these differences in rates and stabilities of charge separation are due to the altered ability of water to solvate and reorganize in the vicinity of the metal interface.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Concentration Fluctuations and Capacitive Response in Dense Ionic Solutions

We use molecular dynamics simulations in a constant potential ensemble to study the effects of solution composition on the electrochemical response of a double layer capacitor. We find that the capacitance first increases with ion concentration following its expected ideal solution behavior, but decreases upon approaching a pure ionic liquid in agreement with recent experimental observations. The non-monotonic behavior of the capacitance as a function of ion concentration results from the competition between the independent motion of solvated ions in the dilute regime and solvation fluctuations in the concentrated regime. Mirroring the capacitance, we find that the characteristic decay length of charge density correlations away from the electrode is also non-monotonic. The correlation length first decreases with ion concentration as a result of better electrostatic screening but increases with ion concentration as a result of enhanced steric interactions. When charge fluctuations induced by correlated ion-solvent fluctuations are large relative to those induced by the pure ionic liquid, such capacitive behavior is expected to be generic.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Interfacial ordering and accompanying divergent capacitance at ionic liquid-metal interfaces

A theory is constructed for dense ionic solutions near charged planar walls that is valid for strong inter-ionic correlations. This theory predicts a fluctuation-induced, first-order transition and spontaneous charge density ordering at the interface, in the presence of an otherwise disordered bulk solution. The surface ordering is driven by applied voltage and results in an anomalous differential capacitance, in agreement with recent simulation results and consistent with experimental observations in a wide array of systems. Explicit forms for the charge density profile and capacitance are given. The theory is compared with numerical results for the charge frustrated Ising model, which is also found to exhibit a voltage driven first-order transition.

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

Nanoscale heterogeneity at the aqueous electrolyte-electrode interface

Using molecular dynamics simulations, we reveal emergent properties of hydrated electrode interfaces that while molecular in origin are integral to the behavior of the system across long times scales and large length scales. Specifically, we describe the impact of a disordered and slowly evolving adsorbed layer of water on the molecular structure and dynamics of the electrolyte solution adjacent to it. Generically, we find that densities and mobilities of both water and dissolved ions are spatially heterogeneous in the plane parallel to the electrode over nanosecond timescales. These and other recent results are analyzed in the context of available experimental literature from surface science and electrochemistry. We speculate on the implications of this emerging microscopic picture on the catalytic proficiency of hydrated electrodes, offering an new direction for study in heterogeneous catalysis at the nanoscale.

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 length and time scales of water's glass transitions

Using a general model for the equilibrium dynamics of supercooled liquids, I compute from molecular properties the emergent length and time scales that govern the nonequilibrium relaxation behavior of amorphous ice prepared by rapid cooling. Upon cooling, the liquid water falls out of equilibrium whereby the temperature dependence of its relaxation time is predicted to change from super-Arrhenius to Arrhenius. A consequence of this crossover is that the location of the apparent glass transition temperature depends logarithmically on cooling rate. Accompanying vitrification is the emergence of a dynamical length-scale, the size of which depends on the cooling rate and varies between angstroms and 10s of nanometers. While this protocol dependence clarifies a number of previous experimental observations for amorphous ice, the arguments are general and can be extended to other glass forming liquids.

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

Nature of the anomalies in the supercooled liquid state of the mW model of water

The thermodynamic properties of the supercooled liquid state of the mW model of water show anomalous behavior. Like in real water, the heat capacity and compressibility sharply increase upon supercooling. One of the possible explanations of these anomalies, the existence of a second (liquid-liquid) critical point, is not supported by simulations for this model. In this work, we reproduce the anomalies of the mW model with two thermodynamic scenarios: one based on a non-ideal "mixture" with two different types of local order of the water molecules, and one based on weak crystallization theory. We show that both descriptions accurately reproduce the model's basic thermodynamic properties. However, the coupling constant required for the power laws implied by weak crystallization theory is too large relative to the regular backgrounds, contradicting assumptions of weak crystallization theory. Fluctuation corrections outside the scope of this work would be necessary to fit the forms predicted by weak crystallization theory. For the two-state approach, the direct computation of the low-density fraction of molecules in the mW model is in agreement with the prediction of the phenomenological equation of state. The non-ideality of the "mixture" of the two states never becomes strong enough to cause liquid-liquid phase separation, also in agreement with simulation results.

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

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
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