Source author record

Raymond Kapral

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

25works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

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

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

preprint2022arXiv

Control of active polymeric filaments by chemically-powered nanomotors

Active materials with distinctive nonequilibrium properties have diverse materials science applications. Active systems are common in living matter, such as the filament network in the cell that is activated by molecular motors, and in materials science as exemplified by hydrogels activated by chemical reactions. Here we describe another class of active polymeric filament systems where the filaments are activated by embedded chemically-powered nanomotors that have catalytic and noncatalytic parts. Chemical reactions on the catalytic surfaces produce forces that act on the polymeric filaments. By changing the nonequilibrium conditions these forces can be made to change sign and thereby compress or expand the filaments. The embedded motors provide both the source of activity and the means to control the filament conformational structure. As an example of control, we show that oscillatory variations of the chemical constraints yield gel-like networks that oscillate between expanded or compressed forms, much like those of hydrogels.

preprint2020arXiv

Active matter, microreversibility, and thermodynamics

Active matter, comprising many active agents interacting and moving in fluids or more complex environments, is a commonly occurring state of matter in biological and physical systems. By its very nature active matter systems exist in nonequilibrium states. In this paper the active agents are small Janus colloidal particles that use chemical energy provided by chemical reactions occurring on their surfaces for propulsion through a diffusiophoretic mechanism. As a result of interactions among these colloids, either directly or through fluid velocity and concentration fields, they may act collectively to form structures such as dynamic clusters. A general nonequilibrium thermodynamics framework for the description of such systems is presented that accounts for both self-diffusiophoresis and diffusiophoresis due to external concentration gradients, and is consistent with microreversibility. It predicts the existence of a reciprocal effect of diffusiophoresis back onto the reaction rate for the entire collection of colloids in the system, as well as the existence of a clustering instability that leads to nonequilibrium inhomogeneous system states.

preprint2018arXiv

Thermodynamics and statistical mechanics of chemically-powered synthetic nanomotors

Colloidal motors without moving parts can be propelled by self-diffusiophoresis, coupling molecular concentration gradients generated by surface chemical reactions to the velocity slip between solid Janus particles and the surrounding fluid solution. The interfacial properties involved in this propulsion mechanism can be described by nonequilibrium thermodynamics and statistical mechanics, disclosing the fundamental role of microreversibility in the coupling between motion and reaction. Among other phenomena, the approach predicts that propulsion by fuel consumption has the reciprocal effect of fuel synthesis by mechanical action.

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.

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.

preprint2016arXiv

Quantum Dynamics in Open Quantum-Classical Systems

Often quantum systems are not isolated and interactions with their environments must be taken into account. In such open quantum systems these environmental interactions can lead to decoherence and dissipation, which have a marked influence on the properties of the quantum system. In many instances the environment is well-approximated by classical mechanics, so that one is led to consider the dynamics of open quantum-classical systems. Since a full quantum dynamical description of large many-body systems is not currently feasible, mixed quantum-classical methods can provide accurate and computationally tractable ways to follow the dynamics of both the system and its environment. This review focuses on quantum-classical Liouville dynamics, one of several quantum-classical descriptions, and discusses the problems that arise when one attempts to combine quantum and classical mechanics, coherence and decoherence in quantum-classical systems, nonadiabatic dynamics, surface-hopping and mean-field theories and their relation to quantum-classical Liouville dynamics, as well as methods for simulating the dynamics.

preprint2016arXiv

Surface hopping from the perspective of quantum-classical Liouville dynamics

Fewest-switches surface hopping is studied in the context of quantum-classical Liouville dynamics. Both approaches are mixed quantum-classical theories that provide a way to describe and simulate the nonadiabatic quantum dynamics of many-body systems. Starting from a surface-hopping solution of the quantum-classical Liouville equation, it is shown how fewest-switches dynamics can be obtained by dropping terms that are responsible for decoherence and restricting the nuclear momentum changes that accompany electronic transitions to those events that occur between population states. The analysis provides information on some of the elements that are essential for the construction of accurate and computationally tractable algorithms for nonadiabatic processes.

preprint2015arXiv

Hydrodynamic Collective Effects of Active Protein Machines in Solution and Lipid Bilayers

The cytoplasm and biomembranes in biological cells contain large numbers of proteins that cyclically change their shapes. They are molecular machines that can function as molecular motors or carry out many other tasks in the cell. We analyze the effects that hydrodynamic flows induced by active proteins have on other passive molecules in solution or membranes. We show that the diffusion constants of passive particles are enhanced substantially. Furthermore, when gradients of active proteins are present, a chemotaxis-like drift of passive particles takes place. In lipid bilayers, the effects are strongly nonlocal, so that active inclusions in the membrane contribute to diffusion enhancement and the drift. The results indicate that the transport properties of passive particles in systems containing active proteins machines operating under nonequilibrium conditions differ from their counterparts in systems at thermal equilibrium.

preprint2015arXiv

Nanoconfined catalytic Ångström-size motors

Self-propelled chemically powered synthetic micron and nano-scale motors are being intensively studied because of the wide range of potential applications that exploit their directed motion. This paper considers even smaller Ångström-size synthetic motors. Such very small motors in bulk solution display effects arising from their self-propulsion. Recent experiments have shown that small-molecule catalysts and single enzyme molecules exhibit properties that have been attributed to their chemical activity. Molecular dynamics is used to investigate the properties of very small Ångström-size synthetic chemically powered sphere-dimer motors in a simple atomic-like solvent confined between walls separated by distances of tens of nanometers. Evidence for strong structural ordering of the motors between the walls, which reflects the finite size of solvent molecules and depends on solvent depletion forces, is provided. Dynamical properties, such as average motor velocity, orientational relaxation, and mean square displacement, are anisotropic and depend on the distance from the walls. This research provides information needed for potential applications that use molecular-scale motors in the complex confined geometries encountered in biology and the laboratory.

preprint2014arXiv

Ångström-scale chemically powered motors

Like their larger micron-scale counterparts, Ångström-scale chemically self-propelled motors use asymmetric catalytic activity to produce self-generated concentration gradients that lead to directed motion. Unlike their micron-scale counterparts, the sizes of Ångström-scale motors are comparable to the solvent molecules in which they move, they are dominated by fluctuations, and they operate on very different time scales. These new features are studied using molecular dynamics simulations of small sphere dimer motors. We show that the ballistic regime is dominated by the thermal speed but the diffusion coefficients of these motors are orders of magnitude larger than inactive dimers. Such small motors may find applications in nano-confined systems or perhaps eventually in the cell.

preprint2014arXiv

Chemistry in Motion: Tiny Synthetic Motors

In this Account, we describe how synthetic motors that operate by self-diffusiophoresis make use of a self-generated concentration gradient to drive motor motion. A description of propulsion by self-diffusiophoresis is presented for Janus particle motors comprising catalytic and noncatalytic faces. The properties of the dynamics of chemically powered motors are illustrated by presenting the results of particle-based simulations of sphere-dimer motors constructed from linked catalytic and noncatalytic spheres. The geometries of both Janus and sphere-dimer motors with asymmetric catalytic activity support the formation of concentration gradients around the motors. Because directed motion can occur only when the system is not in equilibrium, the nature of the environment and the role it plays in motor dynamics are described. Rotational Brownian motion also acts to limit directed motion, and it has especially strong effects for very small motors. We address the following question: how small can motors be and still exhibit effects due to propulsion, even if only to enhance diffusion? Synthetic motors have the potential to transform the manner in which chemical dynamical processes are carried out for a wide range of applications.

preprint2014arXiv

Self-propulsion through symmetry breaking

In addition to self-propulsion by phoretic mechanisms that arises from an asymmetric distribution of reactive species around a catalytic motor, spherical particles with a uniform distribution of catalytic activity may also propel themselves under suitable conditions. Reactive fluctuation-induced asymmetry can give rise to transient concentration gradients which may persist under certain conditions, giving rise to a bifurcation to self-propulsion. The nature of this phenomenon is analyzed in detail, and particle-level simulations are carried out to demonstrate its existence.

preprint2013arXiv

Analysis of geometric phase effects in the quantum-classical Liouville formalism

We analyze two approaches to the quantum-classical Liouville (QCL) formalism that differ in the order of two operations: Wigner transformation and projection onto adiabatic electronic states. The analysis is carried out on a two-dimensional linear vibronic model where geometric phase (GP) effects arising from a conical intersection profoundly affect nuclear dynamics. We find that the Wigner-then-Adiabatic (WA) QCL approach captures GP effects, whereas the Adiabatic-then-Wigner (AW) QCL approach does not. Moreover, the Wigner transform in AW-QCL leads to an ill-defined Fourier transform of double-valued functions. The double-valued character of these functions stems from the nontrivial GP of adiabatic electronic states in the presence of a conical intersection. In contrast, WA-QCL avoids this issue by starting with the Wigner transform of single-valued quantities of the full problem. Since the WA-QCL approach uses solely the adiabatic potentials and non-adiabatic derivative couplings as an input, our results indicate that WA-QCL can capture GP effects in general two-state crossing problems using first-principles electronic structure calculations without prior diabatization or introduction of explicit phase factors.

preprint2013arXiv

Analysis of the Forward-Backward Trajectory Solution for the Mixed Quantum-Classical Liouville Equation

Mixed quantum-classical methods provide powerful algorithms for the simulation of quantum processes in large and complex systems. The forward-backward trajectory solution of the mixed quantum-classical Liouville equation in the mapping basis [J. Chem. Phys. 137, 22A507 (2012)] is one such scheme. It simulates the dynamics via the propagation of forward and backward trajectories of quantum coherent state variables, and the propagation of bath trajectories on a mean-field potential determined jointly by the forward and backward trajectories. An analysis of the properties of this solution, numerical tests of its validity and an investigation of its utility for the study of nonadiabtic quantum processes are given. In addition, we present an extension of this approximate solution that allows one to systematically improve the results. This extension, termed the jump forward-backward trajectory solution, is analyzed and tested in detail and its various implementations are discussed.

preprint2012arXiv

Coarse-Grain Model for Lipid Bilayer Self-Assembly and Dynamics: Multiparticle Collision Description of the Solvent

A mesoscopic coarse-grain model for computationally-efficient simulations of biomembranes is presented. It combines molecular dynamics simulations for the lipids, modeled as elastic chains of beads, with multiparticle collision dynamics for the solvent. Self-assembly of a membrane from a uniform mixture of lipids is observed. Simulations at different temperatures demonstrate that it reproduces the gel and liquid phases of lipid bilayers. Investigations of lipid diffusion in different phases reveals a crossover from subdiffusion to normal diffusion at long times. Macroscopic membrane properties, such as stretching and bending elastic moduli, are determined directly from the mesoscopic simulations. Velocity correlation functions for membrane flows are determined and analyzed.

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.

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.

preprint2012arXiv

Nonadiabatic Dynamics in Open Quantum-Classical Systems: Forward-Backward Trajectory Solution

A new approximate solution to the quantum-classical Liouville equation is derived starting from the formal solution of this equation in forward-backward form. The time evolution of a mixed quantum-classical system described by this equation is obtained in a coherent state basis using the mapping representation, which expresses $N$ quantum degrees of freedom in a 2N-dimensional phase space. The solution yields a simple non-Hamiltonian dynamics in which a set of $N$ coherent state coordinates evolve in forward and backward trajectories while the bath coordinates evolve under the influence of the mean potential that depends on these forward and backward trajectories. It is shown that the solution satisfies the differential form of the quantum-classical Liouville equation exactly. Relations to other mixed quantum-classical and semi-classical schemes are discussed.

preprint2012arXiv

Phoretic self-propulsion: a mesoscopic description of reaction dynamics that powers motion

The fabrication of synthetic self-propelled particles and the experimental investigations of their dynamics have stimulated interest in self-generated phoretic effects that propel nano- and micron-scale objects. Theoretical modeling of these phenomena is often based on a continuum description of the solvent for different phoretic propulsion mechanisms, including, self-electrophoresis, self-diffusiophoresis and self-thermophoresis. The work in this paper considers various types of catalytic chemical reaction at the motor surface and in the bulk fluid that come into play in mesoscopic descriptions of the dynamics. The formulation is illustrated by developing the mesoscopic reaction dynamics for exothermic and dissociation reactions that are used to power motor motion. The results of simulations of the self-propelled dynamics of composite Janus particles by these mechanisms are presented.

preprint2010arXiv

Analysis of the quantum-classical Liouville equation in the mapping basis

The quantum-classical Liouville equation provides a description of the dynamics of a quantum subsystem coupled to a classical environment. Representing this equation in the mapping basis leads to a continuous description of discrete quantum states of the subsystem and may provide an alternate route to the construction of simulation schemes. In the mapping basis the quantum-classical Liouville equation consists of a Poisson bracket contribution and a more complex term. By transforming the evolution equation, term-by-term, back to the subsystem basis, the complex term (excess coupling term) is identified as being due to a fraction of the back reaction of the quantum subsystem on its environment. A simple approximation to quantum-classical Liouville dynamics in the mapping basis is obtained by retaining only the Poisson bracket contribution. This approximate mapping form of the quantum-classical Liouville equation can be simulated easily by Newtonian trajectories. We provide an analysis of the effects of neglecting the presence of the excess coupling term on the expectation values of various types of observables. Calculations are carried out on nonadiabatic population and quantum coherence dynamics for curve crossing models. For these observables, the effects of the excess coupling term enter indirectly in the computation and good estimates are obtained with the simplified propagation.

preprint2009arXiv

Mapping Approach for Quantum-Classical Time Correlation Functions

The calculation of quantum canonical time correlation functions is considered in this paper. Transport properties, such as diffusion and reaction rate coefficients, can be determined from time integrals of these correlation functions. Approximate, quantum-classical expressions for correlation functions, which are amenable to simulation, are derived. These expressions incorporate the full quantum equilibrium structure of the system but approximate the dynamics by quantum-classical evolution where a quantum subsystem is coupled to a classical environment. The main feature of the formulation is the use of a mapping basis where the subsystem quantum states are represented by fictitious harmonic oscillator states. This leads to a full phase space representation of the dynamics that can be simulated without appeal to surface-hopping methods. The results in this paper form the basis for new simulation algorithms for the computation of quantum transport properties of large many-body systems.

preprint2009arXiv

Quantum-Classical Liouville Dynamics in the Mapping Basis

The quantum-classical Liouville equation describes the dynamics of a quantum subsystem coupled to a classical environment. It has been simulated using various methods, notably, surface-hopping schemes. A representation of this equation in the mapping Hamiltonian basis for the quantum subsystem is derived. The resulting equation of motion, in conjunction with expressions for quantum expectation values in the mapping basis, provide another route to the computation of the nonadiabatic dynamics of observables that does not involve surface-hopping dynamics. The quantum-classical Liouville equation is exact for the spin-boson system. This well-known model is simulated using an approximation to the evolution equation in the mapping basis and close agreement with exact quantum results is found.

preprint2007arXiv

Decoherence and Quantum-Classical Master Equation Dynamics

The conditions under which quantum-classical Liouville dynamics may be reduced to a master equation are investigated. Systems that can be partitioned into a quantum-classical subsystem interacting with a classical bath are considered. Starting with an exact non-Markovian equation for the diagonal elements of the density matrix, an evolution equation for the subsystem density matrix is derived. One contribution to this equation contains the bath average of a memory kernel that accounts for all coherences in the system. It is shown to be a rapidly decaying function, motivating a Markovian approximation on this term in the evolution equation. The resulting subsystem density matrix equation is still non-Markovian due to the fact that bath degrees of freedom have been projected out of the dynamics. Provided the computation of non-equilibrium average values or correlation functions is considered, the non-Markovian character of this equation can be removed by lifting the equation into the full phase space of the system. This leads to a trajectory description of the dynamics where each fictitious trajectory accounts for decoherence due to the bath degrees of freedom. The results are illustrated by computations of the rate constant of a model nonadiabatic chemical reaction.

preprint1999arXiv

Spiral Waves in Media with Complex Excitable Dynamics

The structure of spiral waves is investigated in super-excitable reaction-diffusion systems where the local dynamics exhibits multi-looped phase space trajectories. It is shown that such systems support stable spiral waves with broken symmetry and complex temporal dynamics. The main structural features of such waves, synchronization defect lines, are demonstrated to be similar to those of spiral waves in systems with complex-oscillatory dynamics.