-
Noise-cancellation algorithm for simulations of Brownian particles
Authors:
Regina Rusch,
Thomas Franosch,
Gerhard Jung
Abstract:
We investigate the usage of a recently introduced noise-cancellation algorithm for Brownian simulations to enhance the precision of measuring transport properties such as the mean-square displacement or the velocity-autocorrelation function. The algorithm is based on explicitly storing the pseudo-random numbers used to create the randomized displacements in computer simulations and subtracting the…
▽ More
We investigate the usage of a recently introduced noise-cancellation algorithm for Brownian simulations to enhance the precision of measuring transport properties such as the mean-square displacement or the velocity-autocorrelation function. The algorithm is based on explicitly storing the pseudo-random numbers used to create the randomized displacements in computer simulations and subtracting them from the simulated trajectories. The resulting correlation function of the reduced motion is connected to the target correlation function up to a cross-correlation term. Using analytical theory and computer simulations, we demonstrate that the cross-correlation term can be neglected in all three systems studied in this aper. We further expand the algorithm to Monte Carlo simulations and analyze the performance of the algorithm and rationalize that it works particularly well for unbounded, weakly interacting systems in which the precision of the mean-square displacement can be improved by orders of magnitude.
△ Less
Submitted 23 January, 2024;
originally announced January 2024.
-
Learning how to find targets in the micro-world: The case of intermittent active Brownian particles
Authors:
Michele Caraglio,
Harpreet Kaur,
Lukas J. Fiderer,
Andrea López-Incera,
Hans J. Briegel,
Thomas Franosch,
Gorka Muñoz-Gil
Abstract:
Finding the best strategy to minimize the time needed to find a given target is a crucial task both in nature and in reaching decisive technological advances. By considering learning agents able to switch their dynamics between standard and active Brownian motion, here we focus on developing effective target-search behavioral policies for microswimmers navigating a homogeneous environment and sear…
▽ More
Finding the best strategy to minimize the time needed to find a given target is a crucial task both in nature and in reaching decisive technological advances. By considering learning agents able to switch their dynamics between standard and active Brownian motion, here we focus on developing effective target-search behavioral policies for microswimmers navigating a homogeneous environment and searching for targets of unknown position. We exploit Projective Simulation, a reinforcement learning algorithm, to acquire an efficient stochastic policy represented by the probability of switching the phase, i.e. the navigation mode, in response to the type and the duration of the current phase. Our findings reveal that the target-search efficiency increases with the particle's self-propulsion during the active phase and that, while the optimal duration of the passive case decreases monotonically with the activity, the optimal duration of the active phase displays a non-monotonic behavior.
△ Less
Submitted 28 November, 2023;
originally announced November 2023.
-
Adaptive active Brownian particles searching for targets of unknown positions
Authors:
Harpreet Kaur,
Thomas Franosch,
Michele Caraglio
Abstract:
Developing behavioral policies designed to efficiently solve target-search problems is a crucial issue both in nature and in the nanotechnology of the 21st century. Here, we characterize the target-search strategies of simple microswimmers in a homogeneous environment containing sparse targets of unknown positions. The microswimmers are capable of controlling their dynamics by switching between Br…
▽ More
Developing behavioral policies designed to efficiently solve target-search problems is a crucial issue both in nature and in the nanotechnology of the 21st century. Here, we characterize the target-search strategies of simple microswimmers in a homogeneous environment containing sparse targets of unknown positions. The microswimmers are capable of controlling their dynamics by switching between Brownian motion and an active Brownian particle and by selecting the time duration of each of the two phases. The specific conduct of a single microswimmer depends on an internal decision-making process determined by a simple neural network associated with the agent itself. Starting from a population of individuals with random behavior, we exploit the genetic algorithm NeuroEvolution of Augmenting Topologies to show how an evolutionary pressure based on the target-search performances of single individuals helps to find the optimal duration of the two different phases. Our findings reveal that the optimal policy strongly depends on the magnitude of the particle's self-propulsion during the active phase and that a broad spectrum of network topology solutions exists, differing in the number of connections and hidden nodes.
△ Less
Submitted 24 July, 2023;
originally announced July 2023.
-
Active Brownian Particles in a Circular Disk with an Absorbing Boundary
Authors:
Francesco Di Trapani,
Thomas Franosch,
Michele Caraglio
Abstract:
We solve the time-dependent Fokker-Planck equation for a two-dimensional active Brownian particle exploring a circular region with an absorbing boundary. Using the passive Brownian particle as basis states and dealing with the activity as a perturbation, we provide a matrix representation of the Fokker-Planck operator and we express the propagator in terms of the perturbed eigenvalues and eigenfun…
▽ More
We solve the time-dependent Fokker-Planck equation for a two-dimensional active Brownian particle exploring a circular region with an absorbing boundary. Using the passive Brownian particle as basis states and dealing with the activity as a perturbation, we provide a matrix representation of the Fokker-Planck operator and we express the propagator in terms of the perturbed eigenvalues and eigenfunctions. Alternatively, we show that the propagator can be expressed as a combination of the equilibrium eigenstates with weights depending only on time and on the initial conditions, and obeying exact iterative relations. Our solution allows also obtaining the survival probability and the first-passage time distribution. These latter quantities exhibit peculiarities induced by the non-equilibrium character of the dynamics, in particular, they display a strong dependence on the activity of the particle and, to a less extent, also on its rotational diffusivity.
△ Less
Submitted 22 June, 2023;
originally announced June 2023.
-
The electric and magnetic disordered Maxwell equations as eigenvalue problem
Authors:
Walter Schirmacher,
Thomas Franosch,
Marco Leonetti,
Giancarlo Ruocco
Abstract:
We consider Maxwell's equations in a 3-dimensional material, in which both, the electric permittivity, as well as the magnetic permeability, fluctuate in space. Differently from all previous treatments of the disordered electromagnetic problem, we transform Maxwell's equations and the electric and magnetic fields in such a way that the linear operator in the resulting secular equations is manifest…
▽ More
We consider Maxwell's equations in a 3-dimensional material, in which both, the electric permittivity, as well as the magnetic permeability, fluctuate in space. Differently from all previous treatments of the disordered electromagnetic problem, we transform Maxwell's equations and the electric and magnetic fields in such a way that the linear operator in the resulting secular equations is manifestly Hermitian, in order to deal with a proper eigenvalue problem. As an application of our general formalism, we use an appropriate version of the Coherent-Potential approximation (CPA) to calculate the photon density of states and scattering-mean-free path. Applying standard localization theory, we find that in the presence of both electric and magnetic disorder the spectral range of Anderson localization appears to be much larger than in the case of electric (or magnetic) disorder only. Our result could explain the absence of experimental evidence of 3D Anderson localization of light (all the existing experiments has been performed with electric disorder only) and pave the way towards a successful search of this, up to now, elusive phenomenon.
△ Less
Submitted 10 February, 2023;
originally announced February 2023.
-
Characterization and Control of the Run-and-Tumble Dynamics of {\it Escherichia Coli}
Authors:
Christina Kurzthaler,
Yongfeng Zhao,
Nan Zhou,
Jana Schwarz-Linek,
Clemence Devailly,
Jochen Arlt,
Jian-Dong Huang,
Wilson C. K. Poon,
Thomas Franosch,
Julien Tailleur,
Vincent A. Martinez
Abstract:
We characterize the full spatiotemporal gait of populations of swimming {\it Escherichia coli} using renewal processes to analyze the measurements of intermediate scattering functions. This allows us to demonstrate quantitatively how the persistence length of an engineered strain can be controlled by a chemical inducer and to report a controlled transition from perpetual tumbling to smooth swimmin…
▽ More
We characterize the full spatiotemporal gait of populations of swimming {\it Escherichia coli} using renewal processes to analyze the measurements of intermediate scattering functions. This allows us to demonstrate quantitatively how the persistence length of an engineered strain can be controlled by a chemical inducer and to report a controlled transition from perpetual tumbling to smooth swimming. For wild-type {\it E.~coli}, we measure simultaneously the microscopic motility parameters and the large-scale effective diffusivity, hence quantitatively bridging for the first time small-scale directed swimming and macroscopic diffusion.
△ Less
Submitted 21 December, 2022;
originally announced December 2022.
-
Quantitative characterization of run-and-tumble statistics in bulk bacterial suspensions
Authors:
Yongfeng Zhao,
Christina Kurzthaler,
Nan Zhou,
Jana Schwarz-Linek,
Clemence Devailly,
Jochen Arlt,
Jian-Dong Huang,
Wilson C. K. Poon,
Thomas Franosch,
Vincent A. Martinez,
Julien Tailleur
Abstract:
We introduce a numerical method to extract the parameters of run-and-tumble dynamics from experimental measurements of the intermediate scattering function. We show that proceeding in Laplace space is unpractical and employ instead renewal processes to work directly in real time. We first validate our approach against data produced using agent-based simulations. This allows us to identify the leng…
▽ More
We introduce a numerical method to extract the parameters of run-and-tumble dynamics from experimental measurements of the intermediate scattering function. We show that proceeding in Laplace space is unpractical and employ instead renewal processes to work directly in real time. We first validate our approach against data produced using agent-based simulations. This allows us to identify the length and time scales required for an accurate measurement of the motility parameters, including tumbling frequency and swim speed. We compare different models for the run-and-tumble dynamics by accounting for speed variability at the single-cell and population level, respectively. Finally, we apply our approach to experimental data on wild-type Escherichia coli obtained using differential dynamic microscopy.
△ Less
Submitted 21 December, 2022;
originally announced December 2022.
-
Resonant diffusion of a gravitactic circle swimmer
Authors:
Oleksandr Chepizhko,
Thomas Franosch
Abstract:
We investigate the dynamics of a single chiral active particle subject to an external torque due to the presence of a gravitational field. Our computer simulations reveal an arbitrarily strong increase of the long-time diffusivity of the gravitactic agent when the external torque approaches the intrinsic angular drift. We provide analytic expressions for the mean-square displacement in terms of ei…
▽ More
We investigate the dynamics of a single chiral active particle subject to an external torque due to the presence of a gravitational field. Our computer simulations reveal an arbitrarily strong increase of the long-time diffusivity of the gravitactic agent when the external torque approaches the intrinsic angular drift. We provide analytic expressions for the mean-square displacement in terms of eigenfunctions and eigenvalues of the noisy-driven-pendulum problem. The pronounced maximum in the diffusivity is then rationalized by the vanishing of the lowest eigenvalues of the Fokker-Planck equation for the angular motion as the rotational diffusion decreases and the underlying classical bifurcation is approached. A simple harmonic-oscillator picture for the barrier-dominated motion provides a quantitative description for the onset of the resonance while its range of validity is determined by the crossover to a critical-fluctuation-dominated regime.
△ Less
Submitted 6 March, 2023; v1 submitted 29 November, 2022;
originally announced November 2022.
-
Crowding-Enhanced Diffusion: An Exact Theory for Highly Entangled Self-Propelled Stiff Filaments
Authors:
Suvendu Mandal,
Christina Kurzthaler,
Thomas Franosch,
Hartmut Löwen
Abstract:
We study a strongly interacting crowded system of self-propelled stiff filaments by event-driven Brownian dynamics simulations and an analytical theory to elucidate the intricate interplay of crowding and self-propulsion. We find a remarkable increase of the effective diffusivity upon increasing the filament number density by more than one order of magnitude. This counter-intuitive 'crowded is fas…
▽ More
We study a strongly interacting crowded system of self-propelled stiff filaments by event-driven Brownian dynamics simulations and an analytical theory to elucidate the intricate interplay of crowding and self-propulsion. We find a remarkable increase of the effective diffusivity upon increasing the filament number density by more than one order of magnitude. This counter-intuitive 'crowded is faster' behavior can be rationalized by extending the concept of a confining tube pioneered by Doi and Edwards for highly entangled crowded, passive to active systems. We predict a scaling theory for the effective diffusivity as a function of the Péclet number and the filament number density. Subsequently, we show that an exact expression derived for a single self-propelled filament with motility parameters as input can predict the non-trivial spatiotemporal dynamics over the entire range of length and time scales. In particular, our theory captures short-time diffusion, directed swimming motion at intermediate times, and the transition to complete orientational relaxation at long times.
△ Less
Submitted 21 September, 2022;
originally announced September 2022.
-
Target search of active agents crossing high energy barriers
Authors:
Luigi Zanovello,
Michele Caraglio,
Thomas Franosch,
Pietro Faccioli
Abstract:
Target search by active agents in rugged energy landscapes has remained a challenge because standard enhanced sampling methods do not apply to irreversible dynamics. We overcome this non-equilibrium rare-event problem by developing an algorithm generalizing transition-path sampling to active Brownian dynamics. This method is exemplified and benchmarked for a paradigmatic two-dimensional potential…
▽ More
Target search by active agents in rugged energy landscapes has remained a challenge because standard enhanced sampling methods do not apply to irreversible dynamics. We overcome this non-equilibrium rare-event problem by developing an algorithm generalizing transition-path sampling to active Brownian dynamics. This method is exemplified and benchmarked for a paradigmatic two-dimensional potential with a high barrier. We find that even in such a simple landscape the structure and kinetics of the ensemble of transition paths change drastically in the presence of activity. Indeed, active Brownian particles reach the target more frequently than passive Brownian particles, following longer and counterintuitive search patterns.
△ Less
Submitted 13 January, 2021; v1 submitted 24 July, 2020;
originally announced July 2020.
-
Random motion of a circle microswimmer in a random environment
Authors:
Oleksandr Chepizhko,
Thomas Franosch
Abstract:
We simulate the dynamics of a single circle microswimmer exploring a disordered array of fixed obstacles. The interplay of two different types of randomness, quenched disorder and stochastic noise, is investigated to unravel their impact on the transport properties. We compute lines of isodiffusivity as a function of the rotational diffusion coefficient and the obstacle density. We find that incre…
▽ More
We simulate the dynamics of a single circle microswimmer exploring a disordered array of fixed obstacles. The interplay of two different types of randomness, quenched disorder and stochastic noise, is investigated to unravel their impact on the transport properties. We compute lines of isodiffusivity as a function of the rotational diffusion coefficient and the obstacle density. We find that increasing noise or disorder tends to amplify diffusion, yet for large randomness the competition leads to a strong suppression of transport. We rationalize both the suppression and amplification of transport by comparing the relevant time scales of the free motion to the mean period between collisions with obstacles.
△ Less
Submitted 15 July, 2020;
originally announced July 2020.
-
An improved integration scheme for Mode-coupling-theory equations
Authors:
Michele Caraglio,
Lukas Schrack,
Gerhard Jung,
Thomas Franosch
Abstract:
Within the mode-coupling theory (MCT) of the glass transition, we reconsider the numerical schemes to evaluate the MCT functional. Here we propose nonuniform discretizations of the wave number, in contrast to the standard equidistant grid, in order to decrease the number of grid points without losing accuracy. We discuss in detail how the integration scheme on the new grids has to be modified from…
▽ More
Within the mode-coupling theory (MCT) of the glass transition, we reconsider the numerical schemes to evaluate the MCT functional. Here we propose nonuniform discretizations of the wave number, in contrast to the standard equidistant grid, in order to decrease the number of grid points without losing accuracy. We discuss in detail how the integration scheme on the new grids has to be modified from standard Riemann integration. We benchmark our approach by solving the MCT equations numerically for mono-disperse hard disks and hard spheres and by computing the critical packing fraction and the nonergodicity parameters. Our results show that significant improvements in performance can be obtained by employing a nonuniform grid.
△ Less
Submitted 15 July, 2020;
originally announced July 2020.
-
Ideal circle microswimmers in crowded media
Authors:
Oleksandr Chepizhko,
Thomas Franosch
Abstract:
Microswimmers are exposed in nature to crowded environments and their transport properties depend in a subtle way on the interaction with obstacles. Here, we investigate a model for a single ideal circle swimmer exploring a two-dimensional disordered array of impenetrable obstacles. The microswimmer moves on circular orbits in the freely accessible space and follows the surface of an obstacle for…
▽ More
Microswimmers are exposed in nature to crowded environments and their transport properties depend in a subtle way on the interaction with obstacles. Here, we investigate a model for a single ideal circle swimmer exploring a two-dimensional disordered array of impenetrable obstacles. The microswimmer moves on circular orbits in the freely accessible space and follows the surface of an obstacle for a certain time upon collision. Depending on the obstacle density and the radius of the circular orbits, the microswimmer displays either long-range transport or is localized in a finite region. We show that there are transitions from two localized states to a diffusive state each driven by an underlying static percolation transition. We determine the non-equilibrium state diagram and calculate the mean-square displacements and diffusivities by computer simulations. Close to the transition lines transport becomes subdiffusive which is rationalized as a dynamic critical phenomenon.
△ Less
Submitted 16 March, 2020;
originally announced March 2020.
-
Heterogeneities Shape Passive Intracellular Transport
Authors:
P. Witzel,
M. Götz,
Y. Lanoiselée,
T. Franosch,
D. S. Grebenkov,
D. Heinrich
Abstract:
A living cell's interior is one of the most complex and intrinsically dynamic systems, providing an elaborate interplay between cytosolic crowding and ATP-driven motion, which controls cellular functionality. Here, we investigated two distinct fundamental features of the merely passive, not-bio-motor shuttled material transport within the cytoplasm of Dictyostelium discoideum cells: the anomalous…
▽ More
A living cell's interior is one of the most complex and intrinsically dynamic systems, providing an elaborate interplay between cytosolic crowding and ATP-driven motion, which controls cellular functionality. Here, we investigated two distinct fundamental features of the merely passive, not-bio-motor shuttled material transport within the cytoplasm of Dictyostelium discoideum cells: the anomalous non-linear scaling of the mean-squared displacement of a 150nm-diameter particle and non-Gaussian distribution of increments. Relying on single-particle tracking data of 320,000 data points, we performed a systematic analysis of four possible origins for non-Gaussian transport: (1) sample-based variability, (2) rare occurring strong motion events, (3) ergodicity breaking/ageing, and (4) spatio-temporal heterogeneities of the intracellular medium. After excluding the first three reasons, we investigated the remaining hypothesis of a heterogeneous cytoplasm as cause for non-Gaussian transport. A novel fit model with randomly distributed diffusivities implementing medium heterogeneities suits the experimental data. Strikingly, the non-Gaussian feature is independent of the cytoskeleton condition and lag time. This reveals that efficiency and consistency of passive intracellular transport and the related anomalous scaling of the mean-squared displacement are regulated by cytoskeleton components, while cytoplasmic heterogeneities are responsible for the generic, non-Gaussian distribution of increments.
△ Less
Submitted 12 November, 2019;
originally announced November 2019.
-
Anomalous transport in the crowded world of biological cells
Authors:
Felix Höfling,
Thomas Franosch
Abstract:
A ubiquitous observation in cell biology is that diffusion of macromolecules and organelles is anomalous, and a description simply based on the conventional diffusion equation with diffusion constants measured in dilute solution fails. This is commonly attributed to macromolecular crowding in the interior of cells and in cellular membranes, summarising their densely packed and heterogeneous struct…
▽ More
A ubiquitous observation in cell biology is that diffusion of macromolecules and organelles is anomalous, and a description simply based on the conventional diffusion equation with diffusion constants measured in dilute solution fails. This is commonly attributed to macromolecular crowding in the interior of cells and in cellular membranes, summarising their densely packed and heterogeneous structures. The most familiar phenomenon is a power-law increase of the MSD, but there are other manifestations like strongly reduced and time-dependent diffusion coefficients, persistent correlations, non-gaussian distributions of the displacements, heterogeneous diffusion, and immobile particles. After a general introduction to the statistical description of slow, anomalous transport, we summarise some widely used theoretical models: gaussian models like FBM and Langevin equations for visco-elastic media, the CTRW model, and the Lorentz model describing obstructed transport in a heterogeneous environment. Emphasis is put on the spatio-temporal properties of the transport in terms of 2-point correlation functions, dynamic scaling behaviour, and how the models are distinguished by their propagators even for identical MSDs. Then, we review the theory underlying common experimental techniques in the presence of anomalous transport: single-particle tracking, FCS, and FRAP. We report on the large body of recent experimental evidence for anomalous transport in crowded biological media: in cyto- and nucleoplasm as well as in cellular membranes, complemented by in vitro experiments where model systems mimic physiological crowding conditions. Finally, computer simulations play an important role in testing the theoretical models and corroborating the experimental findings. The review is completed by a synthesis of the theoretical and experimental progress identifying open questions for future investigation.
△ Less
Submitted 29 January, 2013;
originally announced January 2013.
-
High-resolution detection of Brownian motion for quantitative Optical Tweezers experiments
Authors:
Matthias Grimm,
Thomas Franosch,
Sylvia Jeney
Abstract:
We have developed a new in situ method to calibrate optical tweezers experiments and simultaneously measure the size of the trapped particle or the viscosity of the surrounding fluid. The positional fluctuations of the trapped particle are recorded with a high-bandwidth photodetector. Next, we compute the mean-square displacement, as well as the velocity autocorrelation function of the sphere and…
▽ More
We have developed a new in situ method to calibrate optical tweezers experiments and simultaneously measure the size of the trapped particle or the viscosity of the surrounding fluid. The positional fluctuations of the trapped particle are recorded with a high-bandwidth photodetector. Next, we compute the mean-square displacement, as well as the velocity autocorrelation function of the sphere and compare it to the theory of Brownian motion including hydrodynamic memory effects. A careful measurement and analysis of the time scales characterizing the dynamics of the harmonically bound sphere fluctuating in a viscous medium then directly yields all relevant parameters. Finally, we test the method for different optical trap strengths, with different bead sizes and in different fluids, and we find excellent agreement with the values provided by the manufacturers. The proposed approach overcomes the most commonly encountered limitations in precision when analyzing the power spectrum of position fluctuations in the region around the corner frequency. These low frequencies are usually prone to errors due to drift, limitations in the detection and trap linearity as well as short acquisition times resulting in poor statistics. Furthermore, the strategy can be generalized to Brownian motion in more complex environments, provided the adequate theories are available.
△ Less
Submitted 6 August, 2012;
originally announced August 2012.
-
Driven Transport on Parallel Lanes with Particle Exclusion and Obstruction
Authors:
Anna Melbinger,
Tobias Reichenbach,
Thomas Franosch,
Erwin Frey
Abstract:
We investigate a driven two-channel system where particles on different lanes mutually obstruct each others motion extending an earlier model by Popkov and Peschel [1]. This obstruction may occur in biological contexts due to steric hinderance where motor proteins carry cargos by `walking' on microtubules. Similarly, the model serves as a description for classical spin transport where charged part…
▽ More
We investigate a driven two-channel system where particles on different lanes mutually obstruct each others motion extending an earlier model by Popkov and Peschel [1]. This obstruction may occur in biological contexts due to steric hinderance where motor proteins carry cargos by `walking' on microtubules. Similarly, the model serves as a description for classical spin transport where charged particles with internal states move unidirectionally on a lattice. Three regimes of qualitatively different behavior are identified depending on the strength of coupling between the lanes. For small and large coupling strengths the model can be mapped to a one-channel problem, whereas a new and rich phase behavior emerges for intermediate ones. We derive an approximate but quantitatively accurate theoretical description in terms of a one-site cluster approximation, and obtain insight into the phase behavior through the current-density relations combined with an extremal-current principle. Our results are confirmed by stochastic simulations.
[1] V. Popkov and I. Peschel, Phys. Rev. E 64, 026126 (2001).
△ Less
Submitted 8 February, 2011; v1 submitted 19 February, 2010;
originally announced February 2010.
-
Generic principles of active transport
Authors:
Mauro Mobilia,
Tobias Reichenbach,
Hauke Hinsch,
Thomas Franosch,
Erwin Frey
Abstract:
Nonequilibrium collective motion is ubiquitous in nature and often results in a rich collection of intringuing phenomena, such as the formation of shocks or patterns, subdiffusive kinetics, traffic jams, and nonequilibrium phase transitions. These stochastic many-body features characterize transport processes in biology, soft condensed matter and, possibly, also in nanoscience. Inspired by these…
▽ More
Nonequilibrium collective motion is ubiquitous in nature and often results in a rich collection of intringuing phenomena, such as the formation of shocks or patterns, subdiffusive kinetics, traffic jams, and nonequilibrium phase transitions. These stochastic many-body features characterize transport processes in biology, soft condensed matter and, possibly, also in nanoscience. Inspired by these applications, a wide class of lattice-gas models has recently been considered. Building on the celebrated {\it totally asymmetric simple exclusion process} (TASEP) and a generalization accounting for the exchanges with a reservoir, we discuss the qualitative and quantitative nonequilibrium properties of these model systems. We specifically analyze the case of a dimeric lattice gas, the transport in the presence of pointwise disorder and along coupled tracks.
△ Less
Submitted 20 December, 2006;
originally announced December 2006.
