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

Phillip L. Geissler

Phillip L. Geissler 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.softBiological PhysicsBiomoleculescond-mat.mtrl-sciSubcellular Processesphysics.chem-ph

Catalog footprint

What is connected

21works
7topics
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

21 published item(s)

preprint2022arXiv

Dielectric response of thin water films: A thermodynamic perspective

The surface of a polar liquid presents a special environment for the solvation and organization of charged solutes, which differ from bulk behaviors in important ways. These differences have motivated many attempts to understand electrostatic response at aqueous interfaces in terms of a spatially varying dielectric permittivity, typically concluding that the dielectric constant of interfacial water is significantly lower than in the bulk liquid. Such analyses, however, are complicated by the potentially nonlocal nature of dielectric response over the short length scales of interfacial heterogeneity. Here we circumvent this problem for thin water films by adopting a thermodynamic approach. Using molecular simulations, we calculate the solvent's contribution to the reversible work of charging a parallel plate capacitor. We find good agreement with a simple dielectric continuum model that assumes bulk dielectric permittivity all the way up to the liquid's boundary, even for very thin ($\sim 1$ nm) films. This comparison requires careful attention to the placement of dielectric boundaries between liquid and vapor, which also resolves apparent discrepancies with dielectric imaging experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Local ice-like structure at the liquid water surface

Experiments and computer simulations have established that liquid water's surfaces can deviate in important ways from familiar bulk behavior. Even in the simplest case of an air-water interface, distinctive layering, orientational biases, and hydrogen bond arrangements have been reported, but an overarching picture of their origins and relationships has been incomplete. Here we show that a broad set of such observations can be understood through an analogy with the basal face of crystalline ice. Using simulations, we specifically demonstrate that water and ice surfaces share a set of structural features suggesting the presence of nanometer-scale ice-like domains at the air-water interface. Most prominent is a shared characteristic layering of molecular density and orientation perpendicular to the interface. Similarities in two-point correlations of hydrogen bond network geometry point to shared ice-like intermolecular structure in the parallel direction as well. Our results bolster and significantly extend previous conceptions of ice-like structure at the liquid's boundary, and suggest that the much-discussed quasi-liquid layer on ice evolves subtly above the melting point into a quasi-ice layer at the surface of liquid water.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Elastic Forces Drive Nonequilibrium Pattern Formation in a Model of Nanocrystal Ion Exchange

Chemical transformations, such as ion exchange, are commonly employed to modify nanocrystal compositions. Yet the mechanisms of these transformations, which often operate far from equilibrium and entail mixing diverse chemical species, remain poorly understood. Here, we explore an idealized model for ion exchange in which a chemical potential drives compositional defects to accumulate at a crystal's surface. These impurities subsequently diffuse inward. We find that the nature of interactions between sites in a compositionally impure crystal strongly impacts exchange trajectories. In particular, elastic deformations which accompany lattice-mismatched species promote spatially modulated patterns in the composition. These same patterns can be produced at equilibrium in core/shell nanocrystals, whose structure mimics transient motifs observed in nonequilibrium trajectories. Moreover, the core of such nanocrystals undergoes a phase transition - from modulated to unstructured - as the thickness or stiffness of the shell is decreased. Our results help explain the varied patterns observed in heterostructured nanocrystals produced by ion exchange and suggest principles for the rational design of compositionally-patterned nanomaterials.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Lattice models for protein organization throughout thylakoid membrane stacks

Proteins in photosynthetic membranes can organize into patterned arrays that span the membrane's lateral size. Attractions between proteins in different layers of a membrane stack can play a key role in this ordering, as was suggested by microscopy and fluorescence spectroscopy and demonstrated by computer simulations of a coarse-grained model. The architecture of thylakoid membranes, however, also provides opportunities for inter-layer interactions that instead disfavor the high protein densities of ordered arrangements. Here we explore the interplay between these opposing driving forces, and in particular the phase transitions that emerge in the periodic geometry of stacked thylakoid membrane discs. We propose a lattice model that roughly accounts for proteins' attraction within a layer and across the stromal gap, steric repulsion across the lumenal gap, and regulation of protein density by exchange with the stroma lamellae. Mean field analysis and computer simulation reveal rich phase behavior for this simple model, featuring a broken-symmetry striped phase that is disrupted at both high and low extremes of chemical potential. The resulting sensitivity of microscopic protein arrangement to the thylakoid's mesoscale vertical structure raises intriguing possibilities for regulation of photosynthetic function.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Origin of Mean-Field Behavior in an Elastic Ising Model

Simple elastic models of spin-crossover compounds are known empirically to exhibit classical critical behavior. We demonstrate how the long-ranged interactions responsible for this behavior arise naturally upon integrating out mechanical fluctuations of such a model. A mean field theory applied to the resulting effective Hamiltonian quantitatively accounts for both thermodynamics and kinetics observed in computer simulations, including a barrier to magnetization reversal that grows extensively with system size. For nanocrystals, which break translational symmetry, a straightforward extension of mean field theory yields similarly accurate results.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Near-optimal protocols in complex nonequilibrium transformations

The development of sophisticated experimental means to control nanoscale systems has motivated efforts to design driving protocols which minimize the energy dissipated to the environment. Computational models are a crucial tool in this practical challenge. We describe a general method for sampling an ensemble of finite-time, nonequilibrium protocols biased towards a low average dissipation. We show that this scheme can be carried out very efficiently in several limiting cases. As an application, we sample the ensemble of low-dissipation protocols that invert the magnetization of a 2D Ising model and explore how the diversity of the protocols varies in response to constraints on the average dissipation. In this example, we find that there is a large set of protocols with average dissipation close to the optimal value, which we argue is a general phenomenon.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the microscopic origin and macroscopic implications of lane formation in mixtures of oppositely-driven particles

Colloidal particles of two types, driven in opposite directions, can segregate into lanes [Vissers et al. Soft Matter 7, 2352 (2011)]. This phenomenon can be reproduced by two-dimensional Brownian dynamics simulations of model particles [Dzubiella et al. Phys. Rev. E 65, 021402 (2002)]. Here we use computer simulation to assess the generality of lane formation with respect to variation of particle type and dynamical protocol. We find that laning results from rectification of diffusion on the scale of a particle diameter: oppositely-driven particles must, in the time taken to encounter each other in the direction of the drive, diffuse in the perpendicular direction by about one particle diameter. This geometric constraint implies that the diffusion constant of a particle, in the presence of those of the opposite type, grows approximately linearly with Peclet number, a prediction confirmed by our numerics over a range of model parameters. Such environment-dependent diffusion is statistically similar to an effective interparticle attraction; consistent with this observation, we find that oppositely-driven non-attractive colloids display features characteristic of the simplest model system possessing both interparticle attractions and persistent motion, the driven Ising lattice gas [Katz, Leibowitz, Spohn, J. Stat. Phys. 34, 497 (1984)]. These features include long-ranged correlations in the disordered regime, and a critical regime characterized by a change in slope of the particle current with Peclet number and by fluctuations that grow with system size. By analogy, we suggest that lane formation in the driven colloid system is in the macroscopic limit a phase transition, but that macroscopic phase separation would not occur in finite time upon starting from disordered initial conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Preserving Correlations Between Trajectories for Efficient Path Sampling

Importance sampling of trajectories has proved a uniquely successful strategy for exploring rare dynamical behaviors of complex systems in an unbiased way. Carrying out this sampling, however, requires an ability to propose changes to dynamical pathways that are substantial, yet sufficiently modest to obtain reasonable acceptance rates. Satisfying this requirement becomes very challenging in the case of long trajectories, due to the characteristic divergences of chaotic dynamics. Here we examine schemes for addressing this problem, which engineer correlation between a trial trajectory and its reference path, for instance using artificial forces. Our analysis is facilitated by a modern perspective on Markov Chain Monte Carlo sampling, inspired by non-equilibrium statistical mechanics, which clarifies the types of sampling strategies that can scale to long trajectories. Viewed in this light, the most promising such strategy guides a trial trajectory by manipulating the sequence of random numbers that advance its stochastic time evolution, as done in a handful of existing methods. In cases where this "noise guidance" synchronizes trajectories effectively, as the Glauber dynamics of a two-dimensional Ising model, we show that efficient path sampling can be achieved even for very long trajectories.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dominance of extreme statistics in a prototype many-body Brownian ratchet

Many forms of cell motility rely on Brownian ratchet mechanisms that involve multiple stochastic processes. We present a computational and theoretical study of the nonequilibrium statistical dynamics of such a many-body ratchet, in the specific form of a growing polymer gel that pushes a diffusing obstacle. We find that oft-neglected correlations among constituent filaments impact steady-state kinetics and significantly deplete the gel's density within molecular distances of its leading edge. These behaviors are captured quantitatively by a self-consistent theory for extreme fluctuations in filaments' spatial distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dynamic phase transitions in simple driven kinetic networks

We analyze the probability distribution for entropy production rates of trajectories evolving on a class of out-of-equilibrium kinetic networks. These networks can serve as simple models for driven dynamical systems, which are of particular importance in biological processes, where energy fluxes typically result in non-equilibrium dynamics. By analyzing the fluctuations in the entropy production, we demonstrate the emergence, in a large system size limit, of a dynamic phase transition between two distinct dynamical regimes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Efficiency and Large Deviations in Time-Asymmetric Stochastic Heat Engines

In a stochastic heat engine driven by a cyclic non-equilibrium protocol, fluctuations in work and heat give rise to a fluctuating efficiency. Using computer simulations and tools from large deviation theory, we have examined these fluctuations in detail for a model two-state engine. We find in general that the form of efficiency probability distributions is similar to those described by Verley et al. [2014 Nat Comm, 5 4721], in particular featuring a local minimum in the long-time limit. In contrast to the time-symmetric engine protocols studied previously, however, this minimum need not occur at the value characteristic of a reversible Carnot engine. Furthermore, while the local minimum may reside at the global minimum of a large deviation rate function, it does not generally correspond to the least likely efficiency measured over finite time. We introduce a general approximation for the finite-time efficiency distribution, $P(η)$, based on large deviation statistics of work and heat, that remains very accurate even when $P(η)$ deviates significantly from its large deviation form.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Heterogeneity-induced large deviations in activity and (in some cases) entropy production

We solve a simple model that supports a dynamic phase transition and show conditions for the existence of the transition. Using methods of large deviation theory we analytically compute the probability distribution for activity and entropy production rates of the trajectories on a large ring with a single heterogeneous link. The corresponding joint rate function demonstrates two dynamical phases - one localized and the other delocalized, but the marginal rate functions do not always exhibit the underlying transition. Symmetries in dynamic order parameters influence the observation of a transition, such that distributions for certain dynamic order parameters need not reveal an underlying dynamical bistability. Solution of our model system furthermore yields the form of the effective Markov transition matrices that generate dynamics in which the two dynamical phases are at coexistence. We discuss the implications of the transition for the response of bacterial cells to antibiotic treatment, arguing that even simple models of a cell cycle lacking an explicit bistability in configuration space will exhibit a bistability of dynamical phases.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Coexistence between fluid and crystalline phases of proteins in photosynthetic membranes

Photosystem II (PSII) and its associated light-harvesting complex II (LHCII) are highly concentrated in the stacked grana regions of photosynthetic thylakoid membranes. Within the membrane, PSII-LHCII supercomplexes can be arranged in disordered packings, ordered arrays, or mixtures thereof. The physical driving forces underlying array formation are unknown, complicating attempts to determine a possible functional role for arrays in regulating light harvesting or energy conversion efficiency. Here we introduce a coarse-grained model of protein interactions in coupled photosynthetic membranes, focusing on just two particle types that feature simple shapes and potential energies motivated by structural studies. Reporting on computer simulations of the model's equilibrium fluctuations, we demonstrate its success in reproducing diverse structural features observed in experiments, including extended PSII-LHCII arrays. Free energy calculations reveal that the appearance of arrays marks a phase transition from the disordered fluid state to a system-spanning crystal, which can easily be arrested by thermodynamic constraints or slow dynamics. The region of fluid-crystal coexistence is broad, encompassing much of the physiologically relevant parameter regime. Our results suggest that grana membranes lie at or near phase coexistence, conferring significant structural and functional flexibility to this densely packed membrane protein system.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Patterns without patches: Hierarchical assembly of complex structures from simple building blocks

Nanoparticles with "sticky patches" have long been proposed as building blocks for the self-assembly of complex structures. The synthetic realizability of such patchy particles, however, greatly lags behind predictions of patterns they could form. Using computer simulations, we show that structures of the same genre can be obtained from a solution of simple isotropic spheres, provided control only over their sizes and a small number of binding affinities. In a first step, finite clusters of well-defined structure and composition emerge from natural dynamics with high yield. In effect a kind of patchy particle, these clusters can further assemble into a variety of complex superstructures, including filamentous networks, ordered sheets, and highly porous crystals.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Rare-event trajectory ensemble analysis reveals metastable dynamical phases in lattice proteins

We explore the dynamical large-deviations of a lattice heteropolymer model of a protein by means of path sampling of trajectories. We uncover the existence of non-equilibrium dynamical phase-transitions in ensembles of trajectories between active and inactive dynamical phases, whose nature depends on properties of the interaction potential. When the full heterogeneity of interactions due to the amino-acid sequence is preserved, as in a fully interacting model or in a heterogeneous version of the Gō model where only native interactions are considered, the transition is between the equilibrium native state and a highly native but kinetically trapped state. In contrast, for the homogeneous Gō model, where there is a single native energy and the sequence plays no role, the dynamical transition is a direct consequence of the static bi-stability between unfolded and native states. In the heterogeneous case the native-active and native-inactive states, despite their static similarity, have widely varying dynamical properties, and the transition between them occurs even in lattice proteins whose sequences are designed to make them optimal folders.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Transferable Pair Potentials for CdS and ZnS Crystals

A set of interatomic pair potentials is developed for CdS and ZnS crystals. We show that a simple energy function, which has been used to describe the properties of CdSe [J. Chem. Phys. 116, 258 (2002)], can be parametrized to accurately describe the lattice and elastic constants, and phonon dispersion relations of bulk CdS and ZnS in the wurtzite and rocksalt crystal structures. The predicted coexistence pressure of the wurtzite and rocksalt structures, as well as the equation of state are in good agreement with experimental observations. These new pair potentials enable the study of a wide range of processes in bulk and nanocrystalline II-VI semiconductor materials.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Consequences of local inter-strand dehybridization for large-amplitude bending fluctuations of double-stranded DNA

The wormlike chain (WLC) model of DNA bending accurately reproduces single-molecule force-extension profiles of long (kilobase) chains. These bending statistics over large scales do not, however, establish a unique microscopic model for elasticity at the 1-10 bp scale, which holds particular interest in biological contexts. Here we examine a class of microscopic models which allow for disruption of base pairing (i.e., a `melt' or `kink', generically an `excitation') and consequently enhanced local flexibility. We first analyze the effect on the excitation free energy of integrating out the spatial degrees of freedom in a wormlike chain. Based on this analysis, we present a formulation of these models that ensures consistency with the well-established thermodynamics of melting in long chains. Using a new method to calculate cyclization statistics of short chains from enhanced-sampling Monte Carlo simulations, we compute J-factors of a meltable wormlike chain (MWLC) over a broad range of chain lengths, including very short molecules (30 bp) that have not yet been explored experimentally. For chains longer than about 120 bp, including most molecules studied to date in the laboratory, we find that melting excitations have little impact on cyclization kinetics. Strong signatures of melting, which might be resolved within typical experimental scatter, emerge only for shorter chains.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Microscopic implications of S-DNA

Recent experiments [J. van Mameren et al. PNAS 106, 18231 (2009)] provide a detailed spatial picture of overstretched DNA, showing that under certain conditions the two strands of the double helix separate at about 65 pN. It was proposed that this observation rules out the existence of an elongated, hybridized form of DNA ('S-DNA'). Here we argue that the S-DNA picture is consistent with the observation of unpeeling during overstretching. We demonstrate that assuming the existence of S-DNA does not imply DNA overstretching to consist of the complete or near-complete conversion of the molecule from B- to S-form. Instead, this assumption implies in general a more complex dynamic coexistence of hybridized and unhybridized forms of DNA. We argue that such coexistence can rationalize several recent experimental observations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Faster strain fluctuation methods through partial volume updates

Elastic systems that are spatially heterogeneous in their mechanical response pose special challenges for molecular simulations. Standard methods for sampling thermal fluctuations of a system's size and shape proceed through a series of homogeneous deformations, whose magnitudes can be severely restricted by its stiffest parts. Here we present a Monte Carlo algorithm designed to circumvent this difficulty, which can be prohibitive in many systems of modern interest. By deforming randomly selected subvolumes alone, it naturally distributes the amplitude of spontaneous elastic fluctuations according to intrinsic heterogeneity. We describe in detail implementations of such "slice moves" that are consistent with detailed balance. Their practical application is illustrated for a random network of cross-linked polymers.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Large-scale simulations of fluctuating biological membranes

We present a simple, and physically motivated, coarse-grained model of a lipid bilayer, suited for micron scale computer simulations. Each ~25 nm^2 patch of bilayer is represented by a spherical particle. Mimicking forces of hydrophobic association, multi-particle interactions suppress the exposure of each sphere's equator to its implicit solvent surroundings. The requirement of high equatorial density stabilizes two-dimensional structures without necessitating crystalline order, allowing us to match both the elasticity and fluidity of natural lipid membranes. We illustrate the model's versatility and realism by characterizing membrane response to a prodding nanorod.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

There and (slowly) back again: Entropy-driven hysteresis in a model of DNA overstretching

When pulled along its axis, double-stranded DNA elongates abruptly at a force of about 65 pN. Two physical pictures have been developed to describe this overstretched state. The first proposes that strong forces induce a phase transition to a molten state consisting of unhybridized single strands. The second picture instead introduces an elongated hybridized phase, called S-DNA, structurally and thermodynamically distinct from standard B-DNA. Little thermodynamic evidence exists to discriminate directly between these competing pictures. Here we show that within a microscopic model of DNA we can distinguish between the dynamics associated with each. In experiment, considerable hysteresis in a cycle of stretching and shortening develops as temperature is increased. Since there are few possible causes of hysteresis in a system whose extent is appreciable in only one dimension, such behavior offers a discriminating test of the two pictures of overstretching. Most experiments are performed upon nicked DNA, permitting the detachment (`unpeeling') of strands. We show that the long-wavelength progression of the unpeeled front generates hysteresis, the character of which agrees with experiment only if we assume the existence of S-DNA. We also show that internal melting (distinct from unpeeling) can generate hysteresis, the degree of which is strongly dependent upon the nonextensive loop entropy of single-stranded DNA.

0 citations0 reviewsNo code linkedOpen access