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

Felix Höfling

Felix Höfling 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-phphysics.comp-phBiological Physicsmath-phmath.MPmath.PRcond-mat.mes-hallDistributed, Parallel, and Cluster Computingnlin.AOnlin.PSPopulations and Evolution

Catalog footprint

What is connected

22works
13topics
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

22 published item(s)

preprint2022arXiv

A probabilistic framework for particle-based reaction-diffusion dynamics using classical Fock space representations

The modeling and simulation of stochastic reaction-diffusion processes is a topic of steady interest that is approached with a wide range of methods. \rev{At the level of particle-resolved descriptions, where chemical reactions are coupled to the spatial diffusion of individual particles, there exist comprehensive numerical simulation schemes, while the corresponding mathematical formalization is relatively underdeveloped. The aim of this paper is to provide a framework to systematically formulate the probabilistic evolution equation, termed chemical diffusion master equation (CDME), that governs particle-based stochastic reaction-diffusion processes. To account for the non-conserved and unbounded particle number of this type of open systems, we employ a classical analogue of the quantum mechanical Fock space that contains the symmetrized probability densities of the many-particle configurations in space. Following field-theoretical ideas of second quantization, we introduce creation and annihilation operators that act on single-particle densities and provide natural representations of symmetrized probability densities as well as of reaction and diffusion operators. These operators allow us to consistently and systematically formulate the CDME for arbitrary reaction schemes. The resulting form of the CDME further serves as the foundation to derive more coarse-grained descriptions of reaction-diffusion dynamics. In this regard, we show that a discretization of the evolution equation by projection onto a Fock subspace generated by a finite set of single-particle densities leads to a generalized form of the well-known reaction-diffusion master equation, which supports non-local reactions between grid cells and which converges properly in the continuum limit.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Chemical diffusion master equation: formulations of reaction--diffusion processes on the molecular level

The chemical diffusion master equation (CDME) describes the probabilistic dynamics of reaction--diffusion systems at the molecular level [del Razo et al., Lett. Math. Phys. 112:49, 2022]; it can be considered the master equation for reaction--diffusion processes. The CDME consists of an infinite ordered family of Fokker--Planck equations, where each level of the ordered family corresponds to a certain number of particles and each particle represents a molecule. The equations at each level describe the spatial diffusion of the corresponding set of particles, and they are coupled to each other via reaction operators --linear operators representing chemical reactions. These operators change the number of particles in the system, and thus transport probability between different levels in the family. In this work, we present three approaches to formulate the CDME and show the relations between them. We further deduce the non-trivial combinatorial factors contained in the reaction operators, and we elucidate the relation to the original formulation of the CDME, which is based on creation and annihilation operators acting on many-particle probability density functions. Finally we discuss applications to multiscale simulations of biochemical systems among other future prospects.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Spontaneous trail formation in populations of auto-chemotactic walkers

We study the formation of trails in populations of self-propelled agents that make oriented deposits of pheromones and also sense such deposits to which they then respond with gradual changes of their direction of motion. Based on extensive off-lattice computer simulations aiming at the scale of insects, e.g., ants, we identify a number of emerging stationary patterns and obtain qualitatively the non-equilibrium state diagram of the model, spanned by the strength of the agent--pheromone interaction and the number density of the population. In particular, we demonstrate the spontaneous formation of persistent, macroscopic trails, and highlight some behaviour that is consistent with a dynamic phase transition. This includes a characterisation of the mass of system-spanning trails as a potential order parameter. We also propose a dynamic model for a few macroscopic observables, including the sub-population size of trail-following agents, which captures the early phase of trail formation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Emergence of molecular friction in liquids: bridging between the atomistic and hydrodynamic pictures

Friction in liquids arises from conservative forces between molecules and atoms. Although the hydrodynamics at the nanoscale is subject of intense research and despite the enormous interest in the non-Markovian dynamics of single molecules and solutes, the onset of friction from the atomistic scale so far could not be demonstrated. Here, we fill this gap based on frequency-resolved friction data from high-precision simulations of three prototypical liquids, including water. Combining with rigorous theoretical arguments, we show that friction in liquids emerges abruptly at a characteristic frequency, beyond which viscous liquids appear as non-dissipative, elastic solids; as a consequence, its origin is non-local in time. Concomitantly, the molecules experience Brownian forces that display persistent correlations and long-lasting memory. A critical test of the generalised Stokes-Einstein relation, mapping the friction of single molecules to the viscoelastic response of the macroscopic sample, disproves the relation for Newtonian fluids, but substantiates it exemplarily for water and a moderately supercooled liquid. The employed approach is suitable to yield novel insights into vitrification mechanisms and the intriguing mechanical properties of soft materials.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Finite-size corrections for the static structure factor of a liquid slab with open boundaries

The presence of a confining boundary can modify the local structure of a liquid markedly. In addition, small samples of finite size are known to exhibit systematic deviations of thermodynamic quantities relative to their bulk values. Here, we consider the static structure factor of a liquid sample in slab geometry with open boundaries at the surfaces, which can be thought of as virtually cutting out the sample from a macroscopically large, homogeneous fluid. This situation is a relevant limit for the interpretation of grazing-incidence diffraction experiments at liquid interfaces and films. We derive an exact, closed expression for the slab structure factor, with the bulk structure factor as the only input. This shows that such free boundary conditions cause significant differences between the two structure factors, in particular at small wavenumbers. An asymptotic analysis of this result yields the scaling exponent and an accurate, useful approximation of these finite-size corrections. Furthermore, the open boundaries permit the interpretation of the slab as an open system, supporting particle exchange with a reservoir. We relate the slab structure factor to the particle number fluctuations and discuss conditions under which the subvolume of the slab represents a grand canonical ensemble with chemical potential $μ$ and temperature $T$. Thus, the open slab serves as a test-bed for the small-system thermodynamics in a $μT$ reservoir. We provide a microscopically justified and exact result for the size dependence of the isothermal compressibility. Our findings are corroborated by simulation data for Lennard-Jones liquids at two representative temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Open Systems out of Equilibrium: Theory and Simulation

We consider the theoretical model of Bergmann and Lebowitz for open systems out of equilibrium and translate its principles in the adaptive resolution molecular dynamics technique (AdResS). We simulate Lennard-Jones fluids with open boundaries in a thermal gradient and find excellent agreement of the stationary responses with results obtained from the simulation of a larger, locally forced closed system. The encouraging results pave the way for a computational treatment of open systems far from equilibrium framed in a well-established theoretical model that avoids possible numerical artifacts and physical misinterpretations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Characterizing anomalous diffusion in crowded polymer solutions and gels over five decades in time with variable-lengthscale fluorescence correlation spectroscopy

The diffusion of macromolecules in cells and in complex fluids is often found to deviate from simple Fickian diffusion. One explanation offered for this behavior is that molecular crowding renders diffusion anomalous, where the mean-squared displacement of the particles scales as $\langle r^2 \rangle \propto t^α$ with $α< 1$. Unfortunately, methods such as fluorescence correlation spectroscopy (FCS) or fluorescence recovery after photobleaching (FRAP) probe diffusion only over a narrow range of lengthscales and cannot directly test the dependence of the mean-squared displacement (MSD) on time. Here we show that variable-lengthscale FCS (VLS-FCS), where the volume of observation is varied over several orders of magnitude, combined with a numerical inversion procedure of the correlation data, allows retrieving the MSD for up to five decades in time, bridging the gap between diffusion experiments performed at different lengthscales. In addition, we show that VLS-FCS provides a way to assess whether the propagator associated with the diffusion is Gaussian or non-Gaussian. We used VLS-FCS to investigate two systems where anomalous diffusion had been previously reported. In the case of dense cross-linked agarose gels, the measured MSD confirmed that the diffusion of small beads was anomalous at short lengthscales, with a cross-over to simple diffusion around $\approx 1~μ$m, consistent with a caged diffusion process. On the other hand, for solutions crowded with marginally entangled dextran molecules, we uncovered an apparent discrepancy between the MSD, found to be linear, and the propagators at short lengthscales, found to be non-Gaussian. These contradicting features call to mind the "anomalous, yet Brownian" diffusion observed in several biological systems, and the recently proposed "diffusing diffusivity" model.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Splitting of the universality class of anomalous transport in crowded media

We investigate the emergence of subdiffusive transport by obstruction in continuum models for molecular crowding. While the underlying percolation transition for the accessible space displays universal behavior, the dynamic properties depend in a subtle non-universal way on the transport through narrow channels. At the same time, the different universality classes are robust with respect to introducing correlations in the obstacle matrix as we demonstrate for quenched hard-sphere liquids as underlying structures. Our results confirm that the microscopic dynamics can dominate the relaxational behavior even at long times, in striking contrast to glassy dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Tube Concept for Entangled Stiff Fibers Predicts Their Dynamics in Space and Time

We study dynamically crowded solutions of stiff fibers deep in the semidilute regime, where the motion of a single constituent becomes increasingly confined to a narrow tube. The spatiotemporal dynamics for wave numbers resolving the motion in the confining tube becomes accessible in Brownian dynamics simulations upon employing a geometry-adapted neighbor list. We demonstrate that in such crowded environments the intermediate scattering function, characterizing the motion in space and time, can be predicted quantitatively by simulating a single freely diffusing phantom needle only, yet with very unusual diffusion coefficients.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Anomalous magneto-transport in disordered structures: classical edge-state percolation

By event-driven molecular dynamics simulations we investigate magneto-transport in a two-dimensional model with randomly distributed scatterers close to the field-induced localization transition. This transition is generated by percolating skipping orbits along the edges of obstacle clusters. The dynamic exponents differ significantly from those of the conventional transport problem on percolating systems, thus establishing a new dynamic universality class. This difference is tentatively attributed to a weak-link scenario, which emerges naturally due to barely overlapping edge trajectories. We make predictions for the frequency-dependent conductivity and discuss implications for active colloidal circle swimmers in a heterogeneous environment.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

H5MD: a structured, efficient, and portable file format for molecular data

We propose a new file format named "H5MD" for storing molecular simulation data, such as trajectories of particle positions and velocities, along with thermodynamic observables that are monitored during the course of the simulation. H5MD files are HDF5 (Hierarchical Data Format) files with a specific hierarchy and naming scheme. Thus, H5MD inherits many benefits of HDF5, e.g., structured layout of multi-dimensional datasets, data compression, fast and parallel I/O, and portability across many programming languages and hardware platforms. H5MD files are self-contained and foster the reproducibility of scientific data and the interchange of data between researchers using different simulation programs and analysis software. In addition, the H5MD specification can serve for other kinds of data (e.g. experimental data) and is extensible to supplemental data, or may be part of an enclosing file structure.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Rounding of the localization transition in model porous media

The generic mechanisms of anomalous transport in porous media are investigated by computer simulations of two-dimensional model systems. In order to bridge the gap between the strongly idealized Lorentz model and realistic models of porous media, two models of increasing complexity are considered: a cherry-pit model with hard-core correlations as well as a soft-potential model. An ideal gas of tracer particles inserted into these structures is found to exhibit anomalous transport which extends up to several decades in time. Also, the self-diffusion of the tracers becomes suppressed upon increasing the density of the systems. These phenomena are attributed to an underlying percolation transition. In the soft potential model the transition is rounded, since each tracer encounters its own critical density according to its energy. Therefore, the rounding of the transition is a generic occurrence in realistic, soft systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Anomalous transport in the crowded world of biological cells

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dynamic arrest in model porous media -- intermediate scattering functions

Heterogeneous media constitute random disordered environments where transport is drastically hindered. Employing extensive molecular dynamics simulations, we investigate the spatio-temporal dynamics of tracer particles in the Lorentz model in the vicinity of the localization transition. There transport becomes anomalous and non-gaussian due to the presence of self-similar spatial structures, and dynamic scaling behavior is anticipated. The interplay of different time and length scales is revealed by the intermediate scattering functions, which are sensitive both to the underlying spatial fractal as well as the anomalous transport. We compare our numerical results in the transition regime to a mode-coupling approach, and find that certain aspects are surprisingly well predicted.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

GPU-accelerated simulation of colloidal suspensions with direct hydrodynamic interactions

Solvent-mediated hydrodynamic interactions between colloidal particles can significantly alter their dynamics. We discuss the implementation of Stokesian dynamics in leading approximation for streaming processors as provided by the compute unified device architecture (CUDA) of recent graphics processors (GPUs). Thereby, the simulation of explicit solvent particles is avoided and hydrodynamic interactions can easily be accounted for in already available, highly accelerated molecular dynamics simulations. Special emphasis is put on efficient memory access and numerical stability. The algorithm is applied to the periodic sedimentation of a cluster of four suspended particles. Finally, we investigate the runtime performance of generic memory access patterns of complexity $O(N^2)$ for various GPU algorithms relying on either hardware cache or shared memory.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Highly accelerated simulations of glassy dynamics using GPUs: caveats on limited floating-point precision

Modern graphics processing units (GPUs) provide impressive computing resources, which can be accessed conveniently through the CUDA programming interface. We describe how GPUs can be used to considerably speed up molecular dynamics (MD) simulations for system sizes ranging up to about 1 million particles. Particular emphasis is put on the numerical long-time stability in terms of energy and momentum conservation, and caveats on limited floating-point precision are issued. Strict energy conservation over 10^8 MD steps is obtained by double-single emulation of the floating-point arithmetic in accuracy-critical parts of the algorithm. For the slow dynamics of a supercooled binary Lennard-Jones mixture, we demonstrate that the use of single-floating point precision may result in quantitatively and even physically wrong results. For simulations of a Lennard-Jones fluid, the described implementation shows speedup factors of up to 80 compared to a serial implementation for the CPU, and a single GPU was found to compare with a parallelised MD simulation using 64 distributed cores.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Anomalous transport of a tracer on percolating clusters

We investigate the dynamics of a single tracer exploring a course of fixed obstacles in the vicinity of the percolation transition for particles confined to the infinite cluster. The mean-square displacement displays anomalous transport, which extends to infinite times precisely at the critical obstacle density. The slowing down of the diffusion coefficient exhibits power-law behavior for densities close to the critical point and we show that the mean-square displacement fulfills a scaling hypothesis. Furthermore, we calculate the dynamic conductivity as response to an alternating electric field. Last, we discuss the non-gaussian parameter as an indicator for heterogeneous dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Development of anomalous diffusion among crowding proteins

In cell membranes, proteins and lipids diffuse in a highly crowded and heterogeneous landscape, where aggregates and dense domains of proteins or lipids obstruct the path of diffusing molecules. In general, hindered motion gives rise to anomalous transport, though the nature of the onset of this behavior is still under debate and difficult to investigate experimentally. Here, we present a systematic study where proteins bound to supported lipid membranes diffuse freely in two dimensions, but are increasingly hindered by the presence of other like proteins. In our model system, the surface coverage of the protein avidin on the lipid bilayer is well controlled by varying the concentration of biotinylated lipid anchors. Using fluorescence correlation spectroscopy (FCS), we measure the time correlation function over long times and convert it to the mean-square displacement of the diffusing proteins. Our approach allows for high precision data and a clear distinction between anomalous and normal diffusion. It enables us to investigate the onset of anomalous diffusion, which takes place when the area coverage of membrane proteins increases beyond approximately 5%. This transition region exhibits pronounced spatial heterogeneities. Increasing the packing fraction further, transport becomes more and more anomalous, manifested in a decrease of the exponent of subdiffusion.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Localization phenomena in models of ion-conducting glass formers

The mass transport in soft-sphere mixtures of small and big particles as well as in the disordered Lorentz gas (LG) model is studied using molecular dynamics (MD) computer simulations. The soft-sphere mixture shows anomalous small-particle diffusion signifying a localization transition separate from the big-particle glass transition. Switching off small-particle excluded volume constraints slows down the small-particle dynamics, as indicated by incoherent intermediate scattering functions. A comparison of logarithmic time derivatives of the mean-squared displacements reveals qualitative similarities between the localization transition in the soft-sphere mixture and its counterpart in the LG. Nevertheless, qualitative differences emphasize the need for further research elucidating the connection between both models.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Persistent memory for a Brownian walker in a random array of obstacles

We show that for particles performing Brownian motion in a frozen array of scatterers long-time correlations emerge in the mean-square displacement. Defining the velocity autocorrelation function (VACF) via the second time-derivative of the mean-square displacement, power-law tails govern the long-time dynamics similar to the case of ballistic motion. The physical origin of the persistent memory is due to repeated encounters with the same obstacle which occurs naturally in Brownian dynamics without involving other scattering centers. This observation suggests that in this case the VACF exhibits these anomalies already at first order in the scattering density. Here we provide an analytic solution for the dynamics of a tracer for a dilute planar Lorentz gas and compare our results to computer simulations. Our result support the idea that quenched disorder provides a generic mechanism for persistent correlations irrespective of the microdynamics of the tracer particle.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Space-resolved dynamics of a tracer in a disordered solid

The dynamics of a tracer particle in a glassy matrix of obstacles displays slow complex transport as the free volume approaches a critical value and the void space falls apart. We investigate the emerging subdiffusive motion of the test particle by extensive molecular dynamics simulations and characterize the spatio-temporal transport in terms of two-time correlation functions, including the time-dependent diffusion coefficient as well as the wavenumber-dependent intermediate scattering function. We rationalize our findings within the framework of critical phenomena and compare our data to a dynamic scaling theory.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Localization Transition of the Two-Dimensional Lorentz Model

We investigate the dynamics of a single tracer particle performing Brownian motion in a two-dimensional course of randomly distributed hard obstacles. At a certain critical obstacle density, the motion of the tracer becomes anomalous over many decades in time, which is rationalized in terms of an underlying percolation transition of the void space. In the vicinity of this critical density the dynamics follows the anomalous one up to a crossover time scale where the motion becomes either diffusive or localized. We analyze the scaling behavior of the time-dependent diffusion coefficient D(t) including corrections to scaling. Away from the critical density, D(t) exhibits universal hydrodynamic long-time tails both in the diffusive as well as in the localized phase.

0 citations0 reviewsNo code linkedOpen access