Refine
Year of publication
Language
- English (84)
Is part of the Bibliography
- yes (84)
Keywords
- anomalous diffusion (9)
- stochastic processes (8)
- ageing (5)
- diffusion (5)
- geometric Brownian motion (4)
- Debye screening (3)
- critical phenomena (3)
- electrostatic interactions (3)
- electrostatics (3)
- financial time series (3)
Institute
- Institut für Physik und Astronomie (84) (remove)
We study the adsorption–desorption transition of polyelectrolyte chains onto planar, cylindrical and spherical surfaces with arbitrarily high surface charge densities by massive Monte Carlo computer simulations. We examine in detail how the well known scaling relations for the threshold transition—demarcating the adsorbed and desorbed domains of a polyelectrolyte near weakly charged surfaces—are altered for highly charged interfaces. In virtue of high surface potentials and large surface charge densities, the Debye–Hückel approximation is often not feasible and the nonlinear Poisson–Boltzmann approach should be implemented. At low salt conditions, for instance, the electrostatic potential from the nonlinear Poisson–Boltzmann equation is smaller than the Debye–Hückel result, such that the required critical surface charge density for polyelectrolyte adsorption σc increases. The nonlinear relation between the surface charge density and electrostatic potential leads to a sharply increasing critical surface charge density with growing ionic strength, imposing an additional limit to the critical salt concentration above which no polyelectrolyte adsorption occurs at all. We contrast our simulations findings with the known scaling results for weak critical polyelectrolyte adsorption onto oppositely charged surfaces for the three standard geometries. Finally, we discuss some applications of our results for some physical–chemical and biophysical systems.
What are the physical laws of the diffusive search of proteins for their specific binding sites on DNA in the presence of the macromolecular crowding in cells? We performed extensive computer simulations to elucidate the protein target search on DNA. The novel feature is the viscoelastic non-Brownian protein bulk diffusion recently observed experimentally. We examine the influence of the protein-DNA binding affinity and the anomalous diffusion exponent on the target search time. In all cases an optimal search time is found. The relative contribution of intermittent three-dimensional bulk diffusion and one-dimensional sliding of proteins along the DNA is quantified. Our results are discussed in the light of recent single molecule tracking experiments, aiming at a better understanding of the influence of anomalous kinetics of proteins on the facilitated diffusion mechanism.
Bias plays an important role in the enhancement of diffusion in periodic potentials. Using the continuous-time random walk in the presence of a bias, we report on an interesting phenomenon for the enhancement of diffusion by the start of the measurement in a random energy landscape. When the variance of the waiting time diverges, in contrast to the bias-free case, the dynamics with bias becomes superdiffusive. In the superdiffusive regime, we find a distinct initial ensemble dependence of the diffusivity. Moreover, the diffusivity can be increased by the aging time when the initial ensemble is not in equilibrium. We show that the time-averaged variance converges to the corresponding ensemble-averaged variance; i.e., ergodicity is preserved. However, trajectory-to-trajectory fluctuations of the time-averaged variance decay unexpectedly slowly. Our findings provide a rejuvenation phenomenon in the superdiffusive regime, that is, the diffusivity for a nonequilibrium initial ensemble gradually increases to that for an equilibrium ensemble when the start of the measurement is delayed.
Wealth and income distributions are known to feature country-specific Pareto exponents for their long power-law tails. To propose a rationale for this, we introduce an agent-based dynamic model and use Monte Carlo simulations to unveil the wealth distributions in closed and open economical systems. The standard money-exchange scenario is supplemented with the position-exchange agent dynamics that vitally affects the Pareto law. Specifically, in closed systems with position-exchange dynamics the power law changes to an exponential shape, while for open systems with traps the Pareto law remains valid.
Native mucus is polymer-based soft-matter material of paramount biological importance. How non-Gaussian and non-ergodic is the diffusive spreading of pathogens in mucus? We study the passive, thermally driven motion of micron-sized tracers in hydrogels of mucins, the main polymeric component of mucus. We report the results of the Bayesian analysis for ranking several diffusion models for a set of tracer trajectories [C. E. Wagner et al., Biomacromolecules, 2017, 18, 3654]. The models with "diffusing diffusivity', fractional and standard Brownian motion are used. The likelihood functions and evidences of each model are computed, ranking the significance of each model for individual traces. We find that viscoelastic anomalous diffusion is often most probable, followed by Brownian motion, while the model with a diffusing diffusion coefficient is only realised rarely. Our analysis also clarifies the distribution of time-averaged displacements, correlations of scaling exponents and diffusion coefficients, and the degree of non-Gaussianity of displacements at varying pH levels. Weak ergodicity breaking is also quantified. We conclude that-consistent with the original study-diffusion of tracers in the mucin gels is most non-Gaussian and non-ergodic at low pH that corresponds to the most heterogeneous networks. Using the Bayesian approach with the nested-sampling algorithm, together with the quantitative analysis of multiple statistical measures, we report new insights into possible physical mechanisms of diffusion in mucin gels.
We perform a detailed statistical analysis of diffusive trajectories of membrane-enclosed vesicles (vacuoles) in the supercrowded cytoplasm of living Acanthamoeba castellanii cells. From the vacuole traces recorded in the center-of-area frame of moving amoebae, we examine the statistics of the time-averaged mean-squared displacements of vacuoles, their generalized diffusion coefficients and anomalous scaling exponents, the ergodicity breaking parameter, the non-Gaussian features of displacement distributions of vacuoles, the displacement autocorrelation function, as well as the distributions of speeds and positions of vacuoles inside the amoeba cells. Our findings deliver novel insights into the internal dynamics of cellular structures in these infectious pathogens. Published under license by AIP Publishing.
In this study we investigate, using all-atom molecular-dynamics computer simulations, the in-plane diffusion of a doxorubicin drug molecule in a thin film of water confined between two silica surfaces. We find that the molecule diffuses along the channel in the manner of a Gaussian diffusion process, but with parameters that vary according to its varying transversal position. Our analysis identifies that four Gaussians, each describing particle motion in a given transversal region, are needed to adequately describe the data. Each of these processes by itself evolves with time at a rate slower than that associated with classical Brownian motion due to a predominance of anticorrelated displacements. Long adsorption events lead to ageing, a property observed when the diffusion is intermittently hindered for periods of time with an average duration which is theoretically infinite. This study presents a simple system in which many interesting features of anomalous diffusion can be explored. It exposes the complexity of diffusion in nanoconfinement and highlights the need to develop new understanding.
Numerous examples for a priori unexpected non-Gaussian behaviour for normal and anomalous diffusion have recently been reported in single-particle tracking experiments. Here, we address the case of non-Gaussian anomalous diffusion in terms of a random-diffusivity mechanism in the presence of power-law correlated fractional Gaussian noise. We study the ergodic properties of this model via examining the ensemble- and time-averaged mean-squared displacements as well as the ergodicity breaking parameter EB quantifying the trajectory-to-trajectory fluctuations of the latter. For long measurement times, interesting crossover behaviour is found as function of the correlation time tau characterising the diffusivity dynamics. We unveil that at short lag times the EB parameter reaches a universal plateau. The corresponding residual value of EB is shown to depend only on tau and the trajectory length. The EB parameter at long lag times, however, follows the same power-law scaling as for fractional Brownian motion. We also determine a corresponding plateau at short lag times for the discrete representation of fractional Brownian motion, absent in the continuous-time formulation. These analytical predictions are in excellent agreement with results of computer simulations of the underlying stochastic processes. Our findings can help distinguishing and categorising certain nonergodic and non-Gaussian features of particle displacements, as observed in recent single-particle tracking experiments.
Quorum-sensing bacteria in a growing colony of cells send out signalling molecules (so-called “autoinducers”) and themselves sense the autoinducer concentration in their vicinity. Once—due to increased local cell density inside a “cluster” of the growing colony—the concentration of autoinducers exceeds a threshold value, cells in this clusters get “induced” into a communal, multi-cell biofilm-forming mode in a cluster-wide burst event. We analyse quantitatively the influence of spatial disorder, the local heterogeneity of the spatial distribution of cells in the colony, and additional physical parameters such as the autoinducer signal range on the induction dynamics of the cell colony. Spatial inhomogeneity with higher local cell concentrations in clusters leads to earlier but more localised induction events, while homogeneous distributions lead to comparatively delayed but more concerted induction of the cell colony, and, thus, a behaviour close to the mean-field dynamics. We quantify the induction dynamics with quantifiers such as the time series of induction events and burst sizes, the grouping into induction families, and the mean autoinducer concentration levels. Consequences for different scenarios of biofilm growth are discussed, providing possible cues for biofilm control in both health care and biotechnology.
Heterogeneous diffusion processes (HDPs) feature a space-dependent diffusivity of the form D(x) = D-0|x|(alpha). Such processes yield anomalous diffusion and weak ergodicity breaking, the asymptotic disparity between ensemble and time averaged observables, such as the mean-squared displacement. Fractional Brownian motion (FBM) with its long-range correlated yet Gaussian increments gives rise to anomalous and ergodic diffusion. Here, we study a combined model of HDPs and FBM to describe the particle dynamics in complex systems with position-dependent diffusivity driven by fractional Gaussian noise. This type of motion is, inter alia, relevant for tracer-particle diffusion in biological cells or heterogeneous complex fluids. We show that the long-time scaling behavior predicted theoretically and by simulations for the ensemble-and time-averaged mean-squared displacements couple the scaling exponents alpha of HDPs and the Hurst exponent H of FBM in a characteristic way. Our analysis of the simulated data in terms of the rescaled variable y similar to |x|(1/(2/(2-alpha)))/t(H) coupling particle position x and time t yields a simple, Gaussian probability density function (PDF), PHDP-FBM(y) = e(-y2)/root pi. Its universal shape agrees well with theoretical predictions for both uni- and bimodal PDF distributions.
We employ Bayesian statistics using the nested-sampling algorithm to compare and rank multiple models of ergodic diffusion (including anomalous diffusion) as well as to assess their optimal parameters for in silico-generated and real time-series. We focus on the recently-introduced model of Brownian motion with "diffusing diffusivity'-giving rise to widely-observed non-Gaussian displacement statistics-and its comparison to Brownian and fractional Brownian motion, also for the time-series with some measurement noise. We conduct this model-assessment analysis using Bayesian statistics and the nested-sampling algorithm on the level of individual particle trajectories. We evaluate relative model probabilities and compute best-parameter sets for each diffusion model, comparing the estimated parameters to the true ones. We test the performance of the nested-sampling algorithm and its predictive power both for computer-generated (idealised) trajectories as well as for real single-particle-tracking trajectories. Our approach delivers new important insight into the objective selection of the most suitable stochastic model for a given time-series. We also present first model-ranking results in application to experimental data of tracer diffusion in polymer-based hydrogels.
Employing extensive Monte Carlo computer simulations, we investigate in detail the properties of multichain adsorption of charged flexible polyelectrolytes (PEs) onto oppositely charged spherical nanoparticles (SNPs). We quantify the conditions of critical adsorption-the phase-separation curve between the adsorbed and desorbed states of the PEs-as a function of the SNP surface-charge density and the concentration of added salt. We study the degree of fluctuations of the PE-SNP electrostatic binding energy, which we use to quantify the emergence of the phase subtransitions, including a series of partially adsorbed PE configurations. We demonstrate how the phase-separation adsorption-desorption boundary shifts and splits into multiple subtransitions at low-salt conditions, thereby generalizing and extending the results for critical adsorption of a single PE onto the SNP. The current findings are relevant for finite concentrations of PEs around the attracting SNP, such as the conditions for PE adsorption onto globular proteins carrying opposite electric charges.
Proteins are capable of locating specific targets on DNA by employing a facilitated diffusion process with intermittent 1D and 3D search steps. Gene colocalisation and coregulation-i.e. the spatial proximity of two communicating genes-is one factor capable of accelerating the target search process along the DNA. We perform Monte Carlo computer simulations and demonstrate the benefits of gene colocalisation for minimising the search time in a model DNA-protein system. We use a simple diffusion model to mimic the search for targets by proteins, produced initially in bursts of multiple proteins and performing the first-passage search on the DNA chain. The behaviour of the mean first-passage times to the target is studied as a function of distance between the initial position of proteins and the DNA target position, as well as versus the concentration of proteins. We also examine the properties of bursty target search kinetics for varying physical-chemical protein-DNA binding affinity. Our findings underline the relevance of colocalisation of production and binding sites for protein search inside biological cells.
We consider diffusion processes with a spatially varying diffusivity giving rise to anomalous diffusion. Such heterogeneous diffusion processes are analysed for the cases of exponential, power-law, and logarithmic dependencies of the diffusion coefficient on the particle position. Combining analytical approaches with stochastic simulations, we show that the functional form of the space-dependent diffusion coefficient and the initial conditions of the diffusing particles are vital for their statistical and ergodic properties. In all three cases a weak ergodicity breaking between the time and ensemble averaged mean squared displacements is observed. We also demonstrate a population splitting of the time averaged traces into fast and slow diffusers for the case of exponential variation of the diffusivity as well as a particle trapping in the case of the logarithmic diffusivity. Our analysis is complemented by the quantitative study of the space coverage, the diffusive spreading of the probability density, as well as the survival probability.
Chromatin domains formed in vivo are characterized by different types of 3D organization of interconnected nucleosomes and architectural proteins. Here, we quantitatively test a hypothesis that the similarities in the structure of chromatin fibers (which we call "structural homology") can affect their mutual electrostatic and protein-mediated bridging interactions. For example, highly repetitive DNA sequences in heterochromatic regions can position nucleosomes so that preferred inter-nucleosomal distances are preserved on the surfaces of neighboring fibers. On the contrary, the segments of chromatin fiber formed on unrelated DNA sequences have different geometrical parameters and lack structural complementarity pivotal for stable association and cohesion. Furthermore, specific functional elements such as insulator regions, transcription start and termination sites, and replication origins are characterized by strong nucleosome ordering that might induce structure-driven iterations of chromatin fibers. We propose that shape-specific protein-bridging interactions facilitate long-range pairing of chromatin fragments, while for closely-juxtaposed fibers electrostatic forces can in addition yield fine-tuned structure-specific recognition and pairing. These pairing effects can account for some features observed for mitotic and inter-phase chromatins.
The semiconductor field-effect platform represents a powerful tool for detecting the adsorption and binding of charged macromolecules with direct electrical readout. In this work, a capacitive electrolyte-insulator-semiconductor (EIS) field-effect sensor consisting of an Al-p-Si-SiO2 structure has been applied for real-time in situ electrical monitoring of the layer-by-layer formation of polyelectrolyte (PE) multilayers (PEM). The PEMs were deposited directly onto the SiO2 surface without any precursor layer or drying procedures. Anionic poly(sodium 4-styrene sulfonate) and cationic weak polyelectrolyte poly(allylamine hydrochloride) have been chosen as a model system. The effect of the ionic strength of the solution, polyelectrolyte concentration, number and polarity of the PE layers on the characteristics of the PEM-modified EIS sensors have been studied by means of capacitance-voltage and constant-capacitance methods. In addition, the thickness, surface morphology, roughness and wettabilityof the PE mono- and multilayers have been characterised by ellipsometry, atomic force microscopy and water contact-angle methods, respectively. To explain potential oscillations on the gate surface and signal behaviour of the capacitive field-effect EIS sensor modified with a PEM, a simplified electrostatic model that takes into account the reduced electrostatic screening of PE charges by mobile ions within the PEM has been proposed and discussed.
We demonstrate the non-ergodicity of a simple Markovian stochastic process with space-dependent diffusion coefficient D(x). For power-law forms D(x) similar or equal to vertical bar x vertical bar(alpha), this process yields anomalous diffusion of the form < x(2)(t)> similar or equal to t(2/(2-alpha)). Interestingly, in both the sub- and superdiffusive regimes we observe weak ergodicity breaking: the scaling of the time-averaged mean-squared displacement <(delta(2)(Delta))over bar> remains linear in the lag time Delta and thus differs from the corresponding ensemble average < x(2)(t)>. We analyse the non-ergodic behaviour of this process in terms of the time-averaged mean- squared displacement (delta(2)) over bar and its random features, i.e. the statistical distribution of (delta(2)) over bar and the ergodicity breaking parameters. The heterogeneous diffusion model represents an alternative approach to non- ergodic, anomalous diffusion that might be particularly relevant for diffusion in heterogeneous media.
We model theoretically the electrostatic effects taking place upon DNA hybridization in dense DNA arrays immobilized on a layer of Au nano-particles deposited on the surface of a field-effect-based DNA capacitive biosensor. We consider the influence of separation of a charged analyte from the sensor surface and the salinity of electrolyte solution, in the framework of both linear and nonlinear Poisson-Boltzmann theories. The latter predicts a substantially weaker sensor signals due to electrostatic saturation effects that is the main conclusion of this paper. We analyze how different physical parameters of dense DNA brushes affect the magnitude of hybridization signals. The list includes the fraction of DNA charge neutralization, the length and spatial conformations of adsorbed DNA molecules, as well as the discreteness of DNA charges. We also examine the effect of Donnan ionic equilibrium in DNA lattices on the sensor response. The validity of theoretical models is contrasted against recent experimental observations on detection of DNA hybridization via its intrinsic electric charge. The sensitivity of such biochemical sensing devices, their detection limit, and DNA hybridization efficiency are briefly discussed in the end.
Polyelectrolytes are macromolecules composed of charged monomers and exhibit unique properties due to the interplay of their flexibility and electrostatic interactions. In solution, they are attracted to oppositely charged surfaces and interfaces and exhibit a transition to an adsorbed state when certain conditions are met concerning the charge densities of the polymer and surface and the properties of the solution. In this review, we discuss two limiting cases for adsorption of flexible polyelectrolytes on curved surfaces: weak and strong adsorption. In the first case, adsorption is strongly influenced by the entropic degrees of freedom of a flexible polyelectrolyte. By contrast, in the strong adsorption limit, electrostatic interactions dominate, which leads to particular adsorption patterns, specifically on spherical surfaces. We discuss the corresponding theoretical approaches, applying a mean-field description for the polymer and the polymer-surface interaction. For weak adsorption, we discuss the critical adsorption behavior by exactly solvable models for planar and spherical geometries and a generic approximation scheme, which is additionally applied to cylindrical surfaces. For strong adsorption, we investigate various polyelectrolyte patterns on cylinders and spheres and evaluate their stability. The results are discussed in the light of experimental results, mostly of DNA adsorption experiments.
Based on extensive Monte Carlo simulations and analytical considerations we study the electrostatically driven adsorption of flexible polyelectrolyte chains onto charged Janus nanospheres. These net-neutral colloids are composed of two equally but oppositely charged hemispheres. The critical binding conditions for polyelectrolyte chains are analysed as function of the radius of the Janus particle and its surface charge density, as well as the salt concentration in the ambient solution. Specifically for the adsorption of finite-length polyelectrolyte chains onto Janus nanoparticles, we demonstrate that the critical adsorption conditions drastically differ when the size of the Janus particle or the screening length of the electrolyte are varied. We compare the scaling laws obtained for the adsorption-desorption threshold to the known results for uniformly charged spherical particles, observing significant disparities. We also contrast the changes to the polyelectrolyte chain conformations close to the surface of the Janus nanoparticles as compared to those for simple spherical particles. Finally, we discuss experimentally relevant physicochemical systems for which our simulations results may become important. In particular, we observe similar trends with polyelectrolyte complexation with oppositely but heterogeneously charged proteins.