The stochastic thermodynamics of a rotating Brownian particle in a gradient flow.
Ontology highlight
ABSTRACT:
SUBMITTER: Lan Y
PROVIDER: S-EPMC4508657 | biostudies-literature | 2015 Jul
REPOSITORIES: biostudies-literature
Publications
The stochastic thermodynamics of a rotating Brownian particle in a gradient flow.
Scientific reports 20150721
We compute the entropy production engendered in the environment from a single Brownian particle which moves in a gradient flow, and show that it corresponds in expectation to classical near-equilibrium entropy production in the surrounding fluid with specific mesoscopic transport coefficients. With temperature gradient, extra terms are found which result from the nonlinear interaction between the particle and the non-equilibrated environment. The calculations are based on the fluctuation relatio ...[more]
Similar Datasets
Project description:A fundamental and intrinsic property of any device or natural system is its relaxation time relax, which is the time it takes to return to equilibrium after the sudden change of a control parameter [1]. Reducing ?relax, is frequently necessary, and is often obtained by a complex feedback process. To overcome the limitations of such an approach, alternative methods based on driving have been recently demonstrated [2, 3], for isolated quantum and classical systems [4-9]. Their extension to open systems in contact with a thermostat is a stumbling block for applications. Here, we design a protocol, named Engineered Swift Equilibration (ESE), that shortcuts time-consuming relaxations, and we apply it to a Brownian particle trapped in an optical potential whose properties can be controlled in time. We implement the process experimentally, showing that it allows the system to reach equilibrium times faster than the natural equilibration rate. We also estimate the increase of the dissipated energy needed to get such a time reduction. The method paves the way for applications in micro and nano devices, where the reduction of operation time represents as substantial a challenge as miniaturization [10].
| S-EPMC5011424 | biostudies-literature
Project description:Some developmental events rely on an electrophoretic force to produce morphogenetic gradients. To quantitatively explore the dynamics of this process, we constructed a stochastic model of an early phase of left-right patterning: serotonin movement through the gap junction-coupled blastomeres of the Xenopus embryo. Particle-tracking simulations showed that a left-right gradient is formed rapidly, quickly reaching a final stable level. The voltage difference was critical for producing a morphogen gradient of the right steepness; gap junctional connectivity and morphogen mass determined the timing of the gradient. Endogenous electrophoresis drives approximately 50% of the particles across more than one cell width, and approximately 20% can travel across half the embryo. The stochastic behavior of the resulting gradients exhibited unexpected complexity among blastomeres' morphogen content, and showed how spatiotemporal variability within individual cells resulted in robust and consistent gradients across the embryonic left-right axis. Analysis of the distribution profile of gradient gain values made quantitative predictions about the conditions that result in the observed background level of laterality defects in unperturbed frog embryos. This work provides a general model that can be used to quantitatively analyze the unexpectedly complex dynamics of morphogens in a wide variety of systems.
| S-EPMC2915568 | biostudies-other
Project description:Thermal fluctuations significantly affect the behavior of microscale systems rotating in shear flow, such as microvortexes, microbubbles, rotating micromotors, microactuators and other elements of lab-on-a-chip devices. The influence of Brownian torque on the motion of individual magnetic microparticles in a rotating magnetic field is experimentally determined using optical tweezers. Rotational Brownian motion induces the flattening of the breakdown transition between the synchronous and asynchronous modes of microparticle rotation. The experimental findings regarding microparticle rotation in the presence of Brownian torque are compared with the results of numerical Brownian dynamics simulations.
| S-EPMC4753417 | biostudies-other
Project description:Small particles attach to liquid-fluid interfaces due to capillary forces. The influence of rotation on the capillary force is largely unexplored, despite being relevant whenever particles roll at a liquid-fluid interface or on a moist solid. Here, we demonstrate that due to contact angle hysteresis, a particle needs to overcome a resistive capillary torque to rotate at an interface. We derive a general model for the capillary torque on a spherical particle. The capillary torque is given by M = γRLk(cos ΘR - cos ΘA), where γ is the interfacial tension, R is the radius of the particle, L is the diameter of the contact line, k = 24/π3 is a geometrical constant, and ΘR and ΘA are the receding and advancing contact angles, respectively. The expression for the capillary torque (normalized by the radius of the particle) is equivalent to the expression for the friction force that a drop experiences when moving on a flat surface. Our theory predicts that capillary torque reduces the mobility of wet granular matter and prevents small (nano/micro) particles from rotating when they are in Brownian motion at an interface.
| S-EPMC8280771 | biostudies-literature
Project description:The Ornstein-Uhlenbeck stochastic process is an exact mathematical model providing accurate representations of many real dynamic processes in systems in a stationary state. When applied to the description of random motion of particles such as that of Brownian particles, it provides exact predictions coinciding with those of the Langevin equation but not restricted to systems in thermal equilibrium but only conditioned to be stationary. Here, we investigate experimentally single particle motion in a two-dimensional granular system in a stationary state, consisting of 1 mm stainless balls on a plane circular surface. The motion of the particles is produced by an alternating magnetic field applied perpendicular to the surface of the container. The mean square displacement of the particles is measured for a range of low concentrations and it is found that following an appropriate scaling of length and time, the short-time experimental curves conform a master curve covering the range of particle motion from ballistic to diffusive in accordance with the description of the Ornstein-Uhlenbeck model.
| S-EPMC5626708 | biostudies-literature
Project description:The structure and free energy of multistranded linear polymer ends evolves as individual subunits are added and lost. Thus, the energetic state of the polymer end is not constant, as assembly theory has assumed. Here we utilize a Brownian dynamics approach to simulate the addition and loss of individual subunits at the polymer tip. Using the microtubule as a primary example, we examined how the structure of the polymer tip dictates the rate at which units are added to and lost from individual protofilaments. We find that freely diffusing subunits arrive less frequently to lagging protofilaments but bind more efficiently, such that there is no kinetic difference between leading and lagging protofilaments within a tapered tip. However, local structure at the nanoscale has up to an order-of-magnitude effect on the rate of addition. Thus, the kinetic on-rate constant, integrated across the microtubule tip (kon,MT), is an ensemble average of the varying individual protofilament on-rate constants (kon,PF). Our findings have implications for both catastrophe and rescue of the dynamic microtubule end, and provide a subnanoscale framework for understanding the mechanism of action of microtubule-associated proteins and microtubule-directed drugs. Although we utilize the specific example of the microtubule here, the findings are applicable to multistranded polymers generally.
| S-EPMC3853320 | biostudies-literature
Project description:BackgroundStochastic effects can be important for the behavior of processes involving small population numbers, so the study of stochastic models has become an important topic in the burgeoning field of computational systems biology. However analysis techniques for stochastic models have tended to lag behind their deterministic cousins due to the heavier computational demands of the statistical approaches for fitting the models to experimental data. There is a continuing need for more effective and efficient algorithms. In this article we focus on the parameter inference problem for stochastic kinetic models of biochemical reactions given discrete time-course observations of either some or all of the molecular species.ResultsWe propose an algorithm for inference of kinetic rate parameters based upon maximum likelihood using stochastic gradient descent (SGD). We derive a general formula for the gradient of the likelihood function given discrete time-course observations. The formula applies to any explicit functional form of the kinetic rate laws such as mass-action, Michaelis-Menten, etc. Our algorithm estimates the gradient of the likelihood function by reversible jump Markov chain Monte Carlo sampling (RJMCMC), and then gradient descent method is employed to obtain the maximum likelihood estimation of parameter values. Furthermore, we utilize flux balance analysis and show how to automatically construct reversible jump samplers for arbitrary biochemical reaction models. We provide RJMCMC sampling algorithms for both fully observed and partially observed time-course observation data. Our methods are illustrated with two examples: a birth-death model and an auto-regulatory gene network. We find good agreement of the inferred parameters with the actual parameters in both models.ConclusionsThe SGD method proposed in the paper presents a general framework of inferring parameters for stochastic kinetic models. The method is computationally efficient and is effective for both partially and fully observed systems. Automatic construction of reversible jump samplers and general formulation of the likelihood gradient function makes our method applicable to a wide range of stochastic models. Furthermore our derivations can be useful for other purposes such as using the gradient information for parametric sensitivity analysis or using the reversible jump samplers for full Bayesian inference. The software implementing the algorithms is publicly available at http://cbcl.ics.uci.edu/sgd.
| S-EPMC2914651 | biostudies-literature
Project description:Stochastic gradient-based Monte Carlo methods such as stochastic gradient Langevin dynamics are useful tools for posterior inference on large scale datasets in many machine learning applications. These methods scale to large datasets by using noisy gradients calculated using a mini-batch or subset of the dataset. However, the high variance inherent in these noisy gradients degrades performance and leads to slower mixing. In this paper, we present techniques for reducing variance in stochastic gradient Langevin dynamics, yielding novel stochastic Monte Carlo methods that improve performance by reducing the variance in the stochastic gradient. We show that our proposed method has better theoretical guarantees on convergence rate than stochastic Langevin dynamics. This is complemented by impressive empirical results obtained on a variety of real world datasets, and on four different machine learning tasks (regression, classification, independent component analysis and mixture modeling). These theoretical and empirical contributions combine to make a compelling case for using variance reduction in stochastic Monte Carlo methods.
| S-EPMC5508544 | biostudies-literature
Project description:BackgroundBrownian Dynamics (BD) is a coarse-grained implicit-solvent simulation method that is routinely used to investigate binary protein association dynamics, but due to its efficiency in handling large simulation volumes and particle numbers it is well suited to also describe many-protein scenarios as they often occur in biological cells.ResultsHere we introduce our "brownmove" simulation package which was designed to handle many-particle problems with varying particle numbers and allows for a very flexible definition of rigid and flexible protein and polymer models. Both a Brownian and a Langevin dynamics (LD) propagation scheme can be used and hydrodynamic interactions are treated efficiently with our recently introduced TEA-HI ansatz [Geyer, Winter, JCP 130 (2009) 114905]. With simulations of constrained polymers and flexible models of spherical proteins we demonstrate that it is crucial to include hydrodynamics when multi-bead models are used in BD or LD simulations. Only then both the translational and the rotational diffusion coefficients and the timescales of the internal dynamics can be reproduced correctly. In the third example project we show how constant density boundary conditions [Geyer et al, JCP 120 (2004) 4573] can be used to set up a non-equilibrium simulation of diffusional transport across an array of fixed obstacles. Finally, we demonstrate how the agglomeration dynamics of multiple particles with attractive patches can be analysed conveniently with the help of a dynamic interaction network.ConclusionsCombining BD and LD propagation, fast hydrodynamics, a flexible protein model, and interfaces for "open" simulation settings, our freely available "brownmove" simulation package constitutes a new platform for coarse-grained many-particle simulations of biologically relevant diffusion and transport processes.
| S-EPMC3098214 | biostudies-literature
Project description:Interacting Particle Systems (IPSs) are used to model spatio-temporal stochastic systems in many disparate areas of science. We design an algorithmic framework that reduces IPS simulation to simulation of well-mixed Chemical Reaction Networks (CRNs). This framework minimizes the number of associated reaction channels and decouples the computational cost of the simulations from the size of the lattice. Decoupling allows our software to make use of a wide class of techniques typically reserved for well-mixed CRNs. We implement the direct stochastic simulation algorithm in the open source programming language Julia. We also apply our algorithms to several complex spatial stochastic phenomena. including a rock-paper-scissors game, cancer growth in response to immunotherapy, and lipid oxidation dynamics. Our approach aids in standardizing mathematical models and in generating hypotheses based on concrete mechanistic behavior across a wide range of observed spatial phenomena.
| S-EPMC7924777 | biostudies-literature