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

Jeremy Schofield

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

ResearcherUnclaimed source record
cond-mat.softcond-mat.stat-mechphysics.chem-phBiological PhysicsBiomoleculescond-mat.dis-nnquant-ph

Catalog footprint

What is connected

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

12 published item(s)

preprint2024arXiv

Microscopic theory of a Janus motor in a non-equilibrium fluid: Surface hydrodynamics and boundary conditions

We present a derivation from first principles of the coupled equations of motion of an active self-diffusiophoretic Janus motor and the hydrodynamic densities of its fluid environment that are nonlinearly displaced from equilibrium. The derivation makes use of time-dependent projection operator techniques defined in terms of slowly varying coarse-grained microscopic densities of the fluid species number, total momentum, and energy. The exact equations of motion are simplified using time-scale arguments, resulting in Markovian equations for the Janus motor linear and angular velocities with average forces and torques that depend on the fluid densities. For a large colloid, the fluid equations are separated into bulk and interfacial contributions, and the conditions under which the dynamics of the fluid densities can be accurately represented by bulk hydrodynamic equations subject to boundary conditions on the colloid are determined. We show how the results for boundary conditions based on continuum theory can be obtained from the molecular description and provide Green-Kubo expressions for all transport coefficients, including the diffusiophoretic coupling and the slip coefficient.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A microscopic model for chemically-powered Janus motors

Very small synthetic motors that make use of chemical reactions to propel themselves in solution hold promise for new applications in the development of new materials, science and medicine. The prospect of such potential applications, along with the fact that systems with many motors or active elements display interesting cooperative phenomena of fundamental interest, has made the study of synthetic motors an active research area. Janus motors, comprising catalytic and noncatalytic hemispheres, figure prominently in experimental and theoretical studies of these systems. While continuum models of Janus particle systems are often used to describe motor dynamics, microscopic models that are able to account for intermolecular interactions, many-body concentration gradients, fluid flows and thermal fluctuations provide a way to explore the dynamical behavior of these complex out-of-equilibrium systems that does not rely on approximations that are often made in continuum theories. The analysis of microscopic models from first principles provides a foundation from which the range of validity and limitations of approximate theories of the dynamics may be assessed. In this paper, a microscopic model for the diffusiophoretic propulsion of Janus motors, where motor interactions with the environment occur only through hard collisions, is constructed, analyzed and compared to theoretical predictions. Microscopic simulations of both single-motor and many-motor systems are carried out to illustrate the results.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Microscopic and continuum descriptions of Janus motor fluid flow fields

Active media, whose constituents are able to move autonomously, display novel features that differ from those of equilibrium systems. In addition to naturally-occurring active systems such as populations of swimming bacteria, active systems of synthetic self-propelled nanomotors have been developed. These synthetic systems are interesting because of their potential applications in a variety of fields. Janus particles, synthetic motors of spherical geometry with one hemisphere that catalyzes the conversion of fuel to product and one noncatalytic hemisphere, can propel themselves in solution by self-diffusiophoresis. In this mechanism the concentration gradient generated by the asymmetric catalytic activity leads to a force on the motor that induces fluid flows in the surrounding medium. These fluid flows are studied in detail through microscopic simulations of Janus motor motion and continuum theory. It is shown that continuum theory is able to capture many but not all features of the dynamics of the Janus motor and the velocity fields of the fluid.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Derivation of a Markov state model of the dynamics of a protein-like chain immersed in an implicit solvent

A Markov state model of the dynamics of a protein-like chain immersed in an implicit hard sphere solvent is derived from first principles for a system of monomers that interact via discontinuous potentials designed to account for local structure and bonding in a coarse-grained sense. The model is based on the assumption that the implicit solvent interacts on a fast time scale with the monomers of the chain compared to the time scale for structural rearrangements of the chain and provides sufficient friction so that the motion of monomers is governed by the Smoluchowski equation. A microscopic theory for the dynamics of the system is developed that reduces to a Markovian model of the kinetics under well-defined conditions. Microscopic expressions for the rate constants that appear in the Markov state model are analyzed and expressed in terms of a temperature-dependent linear combination of escape rates that themselves are independent of temperature. The Markov model is studied by analyzing the eigenvalues and eigenvectors of the matrix of transition rates, and the equilibration process for a simple helix-forming system from an ensemble of initially extended configurations to mainly folded configurations is investigated as a function of temperature for a number of different chain lengths. For short chains, the relaxation is primarily single-exponential and becomes independent of temperature in the low-temperature regime. The profile is more complicated for longer chains, where multi-exponential relaxation behavior is seen at intermediate temperatures followed by a low temperature regime in which the folding becomes rapid and single exponential. It is demonstrated that the behavior of the equilibration profile as the temperature is lowered can be understood in terms of the number of relaxation modes or ``folding pathways'' that contribute to the evolution of the state populations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Free Energy Landscape of Protein-like Chains with Discontinuous Potentials

In this article the configurational space of two simple protein models consisting of polymers composed of a periodic sequence of four different kinds of monomers is studied as a function of temperature. In the protein models, hydrogen bond interactions, electrostatic repulsion, and covalent bond vibrations are modeled by discontinuous step, shoulder and square-well potentials, respectively. The protein-like chains exhibit a secondary alpha helix structure in their folded states at low temperatures, and allow a natural definition of a configuration by considering which beads are bonded. Free energies and entropies of configurations are computed using the parallel tempering method in combination with hybrid Monte Carlo sampling of the canonical ensemble of the discontinuous potential system. The probability of observing the most common configuration is used to analyze the nature of the free energy landscape, and it is found that the model with the least number of possible bonds exhibits a funnel-like free energy landscape at low enough temperature for chains with fewer than 30 beads. For longer proteins, the landscape consists of several minima, where the configuration with the lowest free energy changes significantly by lowering the temperature and the probability of observing the most common configuration never approaches one due to the degeneracy of the lowest accessible potential energy.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Mapping quantum-classical Liouville equation: projectors and trajectories

The evolution of a mixed quantum-classical system is expressed in the mapping formalism where discrete quantum states are mapped onto oscillator states, resulting in a phase space description of the quantum degrees of freedom. By defining projection operators onto the mapping states corresponding to the physical quantum states, it is shown that the mapping quantum-classical Liouville operator commutes with the projection operator so that the dynamics is confined to the physical space. It is also shown that a trajectory-based solution of this equation can be constructed that requires the simulation of an ensemble of entangled trajectories. An approximation to this evolution equation which retains only the Poisson bracket contribution to the evolution operator does admit a solution in an ensemble of independent trajectories but it is shown that this operator does not commute with the projection operators and the dynamics may take the system outside the physical space. The dynamical instabilities, utility and domain of validity of this approximate dynamics are discussed. The effects are illustrated by simulations on several quantum systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Modeling of solvent flow effects in enzyme catalysis under physiological conditions

A stochastic model for the dynamics of enzymatic catalysis in explicit, effective solvents under physiological conditions is presented. Analytically-computed first passage time densities of a diffusing particle in a spherical shell with absorbing boundaries are combined with densities obtained from explicit simulation to obtain the overall probability density for the total reaction cycle time of the enzymatic system. The method is used to investigate the catalytic transfer of a phosphoryl group in a phosphoglycerate kinase-ADP-bis phosphoglycerate system, one of the steps of glycolysis. The direct simulation of the enzyme-substrate binding and reaction is carried out using an elastic network model for the protein, and the solvent motions are described by multiparticle collision dynamics, which incorporates hydrodynamic flow effects. Systems where solvent-enzyme coupling occurs through explicit intermolecular interactions, as well as systems where this coupling is taken into account by including the protein and substrate in the multiparticle collision step, are investigated and compared with simulations where hydrodynamic coupling is absent. It is demonstrated that the flow of solvent particles around the enzyme facilitates the large-scale hinge motion of the enzyme with bound substrates, and has a significant impact on the shape of the probability densities and average time scales of substrate binding for substrates near the enzyme, the closure of the enzyme after binding, and the overall time of completion of the cycle.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Simulation of tethered oligomers in nanochannels using multi-particle collision dynamics

The effect of a high Reynold's number, pressure-driven flow of a compressible gas on the conformation of an oligomer tethered to the wall of a square-channel is studied under both ideal solvent and poor solvent conditions using a hybrid multiparticle collision dynamics and molecular dynamics algorithm. Unlike previous studies, the flow field contains an elongational component in addition to a shear component as well as fluid slip near the walls and results in a Schmidt number for the polymer beads that is less than unity. In both solvent regimes the oligomer is found to extend in the direction of flow. Under the ideal solvent conditions, torsional twisting of the chain and aperiodic cyclical dynamics are observed for the end of the oligomer. Under poor solvent conditions, a metastable helix forms in the end of the chain despite the lack of any attractive potential between beads in the oligomeric chain. The formation of the helix is postulated to be the result of a solvent induced chain collapse that has been confined to a single dimension by a strong flow field.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Free Energy Landscape of a Protein-Like Chain in a Fluid with Discontinuous Potentials

The free energy landscape of a protein-like chain in a fluid was studied by combining discontinuous molecular dynamics and parallel tempering. The model protein is a repeating sequence of four different beads, with interactions mimicking those in real proteins. Neighbor distances and angles are restricted to physical ranges and one out of the four kinds of beads can form hydrogen bonds with each other, except if they are too close in the chain. In contrast to earlier studies of this model, an explicit square-well solvent is included. Beads that can form intra-chain hydrogen bonds, can also form (weaker) hydrogen bonds with solvent molecules, while other beads are insoluble. By categorizing the protein configurations according to their intra-chain bonds, one can distinguish unfolded, helical, and collapsed helical structures. Simulations for chains of 15, 20 and 25 beads show that at low temperatures, the most likely structures are helical or collapsed helical, despite the low entropy of these structures. The temperature at which helical structures become dominant is higher than in the absence of a solvent. The cooperative effect of the solvent is attributed to the presence of hydrophobic beads. A phase transition of the solvent prevented the simulations of the 20-bead and 25-bead chains of reaching low enough temperatures to confirm whether the free energy landscape is funnel-shaped, although the results do not contradict that possibility.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Constructing smooth potentials of mean force, radial, distribution functions and probability densities from sampled data

In this paper a method of obtaining smooth analytical estimates of probability densities, radial distribution functions and potentials of mean force from sampled data in a statistically controlled fashion is presented. The approach is general and can be applied to any density of a single random variable. The method outlined here avoids the use of histograms, which require the specification of a physical parameter (bin size) and tend to give noisy results. The technique is an extension of the Berg-Harris method [B.A. Berg and R.C. Harris, Comp. Phys. Comm. 179, 443 (2008)], which is typically inaccurate for radial distribution functions and potentials of mean force due to a non-uniform Jacobian factor. In addition, the standard method often requires a large number of Fourier modes to represent radial distribution functions, which tends to lead to oscillatory fits. It is shown that the issues of poor sampling due to a Jacobian factor can be resolved using a biased resampling scheme, while the requirement of a large number of Fourier modes is mitigated through an automated piecewise construction approach. The method is demonstrated by analyzing the radial distribution functions in an energy-discretized water model. In addition, the fitting procedure is illustrated on three more applications for which the original Berg-Harris method is not suitable, namely, a random variable with a discontinuous probability density, a density with long tails, and the distribution of the first arrival times of a diffusing particle to a sphere, which has both long tails and short-time structure. In all cases, the resampled, piecewise analytical fit outperforms the histogram and the original Berg-Harris method.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Glassy dynamics and domains: exact results for the East model

A general matrix-based scheme for analyzing the long-time dynamics in kinetically constrained models such as the East model is presented. The treatment developed here is motivated by the expectation that slowly-relaxing spin domains of arbitrary size govern the highly cooperative events that lead to spin relaxation at long times. To account for the role of large spin domains in the dynamics, a complete basis expressed in terms of domains of all sizes is introduced. It is first demonstrated that accounting for single domains of all possible sizes leads to a simple analytical result for the two-time single-spin correlation function in the East model that is in excellent quantitative agreement with simulation data for equilibrium spin up density values c greater or equal to 0.6. It is then shown that including also two neighboring domains leads to a closed expression that describes the slow relaxation of the system down to c approximately 0.3. Ingredients of generalizing the method to lower values of c are also provided, as well as to other models. The main advantage of this approach is that it gives explicit analytical results and that it requires neither an arbitrary closure for the memory kernel nor the construction of an irreducible memory kernel. It also allows one to calculate quantities that measure heterogeneity in the same framework, as is illustrated on the neighbor-pair correlation function and the distribution of relaxation times.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Langevin equation for the extended Rayleigh model with an asymmetric bath

In this paper a one-dimensional model of two infinite gases separated by a movable heavy piston is considered. The non-linear Langevin equation for the motion of the piston is derived from first principles for the case when the thermodynamic parameters and/or the molecular masses of gas particles on left and right sides of the piston are different. Microscopic expressions involving time correlation functions of the force between bath particles and the piston are obtained for all parameters appearing in the non-linear Langevin equation. It is demonstrated that the equation has stationary solutions corresponding to directional fluctuation-induced drift in the absence of systematic forces. In the case of ideal gases interacting with the piston via a quadratic repulsive potential, the model is exactly solvable and explicit expressions for the kinetic coefficients in the non-linear Langevin equation are derived. The transient solution of the non-linear Langevin equation is analyzed perturbatively and it is demonstrated that previously obtained results for systems with the hard-wall interaction are recovered.

0 citations0 reviewsNo code linkedOpen access