Refine
Year of publication
- 2020 (17) (remove)
Document Type
- Article (7)
- Postprint (7)
- Doctoral Thesis (3)
Language
- English (17)
Is part of the Bibliography
- yes (17)
Keywords
- diffusion (17) (remove)
Institute
We study Brownian motion in a confining potential under a constant-rate resetting to a reset position x(0). The relaxation of this system to the steady-state exhibits a dynamic phase transition, and is achieved in a light cone region which grows linearly with time. When an absorbing boundary is introduced, effecting a symmetry breaking of the system, we find that resetting aids the barrier escape only when the particle starts on the same side as the barrier with respect to the origin. We find that the optimal resetting rate exhibits a continuous phase transition with critical exponent of unity. Exact expressions are derived for the mean escape time, the second moment, and the coefficient of variation (CV).
We study properties of magnetohydrodynamic (MHD) eigenmodes by decomposing the data of MHD simulations into linear MHD modes-namely, the Alfven, slow magnetosonic, and fast magnetosonic modes. We drive turbulence with a mixture of solenoidal and compressive driving while varying the Alfven Mach number (M-A), plasma beta, and the sonic Mach number from subsonic to transsonic. We find that the proportion of fast and slow modes in the mode mixture increases with increasing compressive forcing. This proportion of the magnetosonic modes can also become the dominant fraction in the mode mixture. The anisotropy of the modes is analyzed by means of their structure functions. The Alfven-mode anisotropy is consistent with the Goldreich-Sridhar theory. We find a transition from weak to strong Alfvenic turbulence as we go from low to high M-A. The slow-mode properties are similar to the Alfven mode. On the other hand, the isotropic nature of fast modes is verified in the cases where the fast mode is a significant fraction of the mode mixture. The fast-mode behavior does not show any transition in going from low to high M-A. We find indications that there is some interaction between the different modes, and the properties of the dominant mode can affect the properties of the weaker modes. This work identifies the conditions under which magnetosonic modes can be a major fraction of turbulent astrophysical plasmas, including the regime of weak turbulence. Important astrophysical implications for cosmic-ray transport and magnetic reconnection are discussed.
We study the extremal properties of a stochastic process xt defined by the Langevin equation ẋₜ =√2Dₜ ξₜ, in which ξt is a Gaussian white noise with zero mean and Dₜ is a stochastic‘diffusivity’, defined as a functional of independent Brownian motion Bₜ.We focus on threechoices for the random diffusivity Dₜ: cut-off Brownian motion, Dₜt ∼ Θ(Bₜ), where Θ(x) is the Heaviside step function; geometric Brownian motion, Dₜ ∼ exp(−Bₜ); and a superdiffusive process based on squared Brownian motion, Dₜ ∼ B²ₜ. For these cases we derive exact expressions for the probability density functions of the maximal positive displacement and of the range of the process xₜ on the time interval ₜ ∈ (0, T).We discuss the asymptotic behaviours of the associated probability density functions, compare these against the behaviour of the corresponding properties of standard Brownian motion with constant diffusivity (Dₜ = D0) and also analyse the typical behaviour of the probability density functions which is observed for a majority of realisations of the stochastic diffusivity process.
We study the extremal properties of a stochastic process xt defined by the Langevin equation ẋₜ =√2Dₜ ξₜ, in which ξt is a Gaussian white noise with zero mean and Dₜ is a stochastic‘diffusivity’, defined as a functional of independent Brownian motion Bₜ.We focus on threechoices for the random diffusivity Dₜ: cut-off Brownian motion, Dₜt ∼ Θ(Bₜ), where Θ(x) is the Heaviside step function; geometric Brownian motion, Dₜ ∼ exp(−Bₜ); and a superdiffusive process based on squared Brownian motion, Dₜ ∼ B²ₜ. For these cases we derive exact expressions for the probability density functions of the maximal positive displacement and of the range of the process xₜ on the time interval ₜ ∈ (0, T).We discuss the asymptotic behaviours of the associated probability density functions, compare these against the behaviour of the corresponding properties of standard Brownian motion with constant diffusivity (Dₜ = D0) and also analyse the typical behaviour of the probability density functions which is observed for a majority of realisations of the stochastic diffusivity process.
Research on novel and advanced biomaterials is an indispensable step towards their applications in desirable fields such as tissue engineering, regenerative medicine, cell culture, or biotechnology. The work presented here focuses on such a promising material: polyelectrolyte multilayer (PEM) composed of hyaluronic acid (HA) and poly(L-lysine) (PLL). This gel-like polymer surface coating is able to accumulate (bio-)molecules such as proteins or drugs and release them in a controlled manner. It serves as a mimic of the extracellular matrix (ECM) in composition and intrinsic properties. These qualities make the HA/PLL multilayers a promising candidate for multiple bio-applications such as those mentioned above. The work presented aims at the development of a straightforward approach for assessment of multi-fractional diffusion in multilayers (first part) and at control of local molecular transport into or from the multilayers by laser light trigger (second part).
The mechanism of the loading and release is governed by the interaction of bioactives with the multilayer constituents and by the diffusion phenomenon overall. The diffusion of a molecule in HA/PLL multilayers shows multiple fractions of different diffusion rate. Approaches, that are able to assess the mobility of molecules in such a complex system, are limited. This shortcoming motivated the design of a novel evaluation tool presented here.
The tool employs a simulation-based approach for evaluation of the data acquired by fluorescence recovery after photobleaching (FRAP) method. In this approach, possible fluorescence recovery scenarios are primarily simulated and afterwards compared with the data acquired while optimizing parameters of a model until a sufficient match is achieved. Fluorescent latex particles of different sizes and fluorescein in an aqueous medium are utilized as test samples validating the analysis results. The diffusion of protein cytochrome c in HA/PLL multilayers is evaluated as well.
This tool significantly broadens the possibilities of analysis of spatiotemporal FRAP data, which originate from multi-fractional diffusion, while striving to be widely applicable. This tool has the potential to elucidate the mechanisms of molecular transport and empower rational engineering of the drug release systems.
The second part of the work focuses on the fabrication of such a spatiotemporarily-controlled drug release system employing the HA/PLL multilayer. This release system comprises different layers of various functionalities that together form a sandwich structure. The bottom layer, which serves as a reservoir, is formed by HA/PLL PEM deposited on a planar glass substrate. On top of the PEM, a layer of so-called hybrids is deposited. The hybrids consist of thermoresponsive poly(N-isopropylacrylamide) (PNIPAM) -based hydrogel microparticles with surface-attached gold nanorods. The layer of hybrids is intended to serve as a gate that controls the local molecular transport through the PEM–solution-interface. The possibility of stimulating the molecular transport by near-infrared (NIR) laser irradiation is being explored.
From several tested approaches for the deposition of hybrids onto the PEM surface, the drying-based approach was identified as optimal. Experiments, that examine the functionality of the fabricated sandwich at elevated temperature, document the reversible volume phase transition of the PEM-attached hybrids while sustaining the sandwich stability. Further, the gold nanorods were shown to effectively absorb light radiation in the tissue- and cell-friendly NIR spectral region while transducing the energy of light into heat. The rapid and reversible shrinkage of the PEM-attached hybrids was thereby achieved. Finally, dextran was employed as a model transport molecule. It loads into the PEM reservoir in a few seconds with the partition constant of 2.4, while it spontaneously releases in a slower, sustained manner. The local laser irradiation of the sandwich, which contains the fluorescein isothiocyanate tagged dextran, leads to a gradual reduction of fluorescence intensity in the irradiated region.
The release system fabricated employs renowned photoresponsivity of the hybrids in an innovative setting. The results of the research are a step towards a spatially-controlled on-demand drug release system that paves the way to spatiotemporally controlled drug release.
The approaches developed in this work have the potential to elucidate the molecular dynamics in ECM and to foster engineering of multilayers with properties tuned to mimic the ECM. The work aims at spatiotemporal control over the diffusion of bioactives and their presentation to the cells.
The two hallmark features of Brownian motion are the linear growth < x2(t)> = 2Ddt of the mean squared displacement (MSD) with diffusion coefficient D in d spatial dimensions, and the Gaussian distribution of displacements. With the increasing complexity of the studied systems deviations from these two central properties have been unveiled over the years. Recently, a large variety of systems have been reported in which the MSD exhibits the linear growth in time of Brownian (Fickian) transport, however, the distribution of displacements is pronouncedly non-Gaussian (Brownian yet non-Gaussian, BNG). A similar behaviour is also observed for viscoelastic-type motion where an anomalous trend of the MSD, i.e., <x2(t)> ~ ta, is combined with a priori unexpected non-Gaussian distributions (anomalous yet non-Gaussian, ANG). This kind of behaviour observed in BNG and ANG diffusions has been related to the presence of heterogeneities in the systems and a common approach has been established to address it, that is, the random diffusivity approach.
This dissertation explores extensively the field of random diffusivity models. Starting from a chronological description of all the main approaches used as an attempt of describing BNG and ANG diffusion, different mathematical methodologies are defined for the resolution and study of these models. The processes that are reported in this work can be classified in three subcategories, i) randomly-scaled Gaussian processes, ii) superstatistical models and iii) diffusing diffusivity models, all belonging to the more general class of random diffusivity models. Eventually, the study focuses more on BNG diffusion, which is by now well-established and relatively well-understood. Nevertheless, many examples are discussed for the description of ANG diffusion, in order to highlight the possible scenarios which are known so far for the study of this class of processes.
The second part of the dissertation deals with the statistical analysis of random diffusivity processes. A general description based on the concept of moment-generating function is initially provided to obtain standard statistical properties of the models. Then, the discussion moves to the study of the power spectral analysis and the first passage statistics for some particular random diffusivity models. A comparison between the results coming from the random diffusivity approach and the ones for standard Brownian motion is discussed. In this way, a deeper physical understanding of the systems described by random diffusivity models is also outlined.
To conclude, a discussion based on the possible origins of the heterogeneity is sketched, with the main goal of inferring which kind of systems can actually be described by the random diffusivity approach.
We consider the first-passage problem for N identical independent particles that are initially released uniformly in a finite domain Ω and then diffuse toward a reactive area Γ, which can be part of the outer boundary of Ω or a reaction centre in the interior of Ω. For both cases of perfect and partial reactions, we obtain the explicit formulas for the first two moments of the fastest first-passage time (fFPT), i.e., the time when the first out of the N particles reacts with Γ. Moreover, we investigate the full probability density of the fFPT. We discuss a significant role of the initial condition in the scaling of the average fFPT with the particle number N, namely, a much stronger dependence (1/N and 1/N² for partially and perfectly reactive targets, respectively), in contrast to the well known inverse-logarithmic behaviour found when all particles are released from the same fixed point. We combine analytic solutions with scaling arguments and stochastic simulations to rationalise our results, which open new perspectives for studying the relevance of multiple searchers in various situations of molecular reactions, in particular, in living cells.
We consider the first-passage problem for N identical independent particles that are initially released uniformly in a finite domain Ω and then diffuse toward a reactive area Γ, which can be part of the outer boundary of Ω or a reaction centre in the interior of Ω. For both cases of perfect and partial reactions, we obtain the explicit formulas for the first two moments of the fastest first-passage time (fFPT), i.e., the time when the first out of the N particles reacts with Γ. Moreover, we investigate the full probability density of the fFPT. We discuss a significant role of the initial condition in the scaling of the average fFPT with the particle number N, namely, a much stronger dependence (1/N and 1/N² for partially and perfectly reactive targets, respectively), in contrast to the well known inverse-logarithmic behaviour found when all particles are released from the same fixed point. We combine analytic solutions with scaling arguments and stochastic simulations to rationalise our results, which open new perspectives for studying the relevance of multiple searchers in various situations of molecular reactions, in particular, in living cells.
Small moonlets or moons embedded in dense planetary rings create S-shaped density modulations called propellers if their masses are smaller than a certain threshold, alternatively they create a circumferential gap in the disk if the embedded body’s mass exceeds this threshold (Spahn and Sremčević, 2000). The gravitational perturber scatters the ring particles, depletes the disk’s density, and, thus, clears a gap, whereas counteracting viscous diffusion of the ring material has the tendency to close the created gap, thereby forming a propeller. Propeller objects were predicted by Spahn and Sremčević (2000) and Sremčević et al. (2002) and were later discovered by the Cassini space probe (Tiscareno et al., 2006, Sremčević et al., 2007, Tiscareno et al., 2008, and Tiscareno et al., 2010). The ring moons Pan and Daphnis are massive enough to maintain the circumferential Encke and Keeler gaps in Saturn’s A ring and were detected by Showalter (1991) and Porco (2005) in Voyager and Cassini images, respectively. In this thesis, a nonlinear axisymmetric diffusion model is developed to describe radial density profiles of circumferential gaps in planetary rings created by embedded moons (Grätz et al., 2018). The model accounts for the gravitational scattering of the ring particles by the embedded moon and for the counteracting viscous diffusion of the ring matter back into the gap. With test particle simulations it is shown that the scattering of the ring particles passing the moon is larger for small impact parameters than estimated by Goldreich and Tremaine (1980). This is especially significant for the modeling of the Keeler gap. The model is applied to the Encke and Keeler gaps with the aim to estimate the shear viscosity of the ring in their vicinities. In addition, the model is used to analyze whether tiny icy moons whose dimensions lie below Cassini’s resolution capabilities would be able to cause the poorly understood gap structure of the C ring and the Cassini Division. One of the most intriguing facets of Saturn’s rings are the extremely sharp edges of the Encke and Keeler gaps: UVIS-scans of their gap edges show that the optical depth drops from order unity to zero over a range of far less than 100 m, a spatial scale comparable to the ring’s vertical extent. This occurs despite the fact that the range over which a moon transfers angular momentum onto the ring material is much larger. Borderies et al. (1982, 1989) have shown that this striking feature is likely related to the local reversal of the usually outward-directed viscous transport of angular momentum in strongly perturbed regions. We have revised the Borderies et al. (1989) model using a granular flow model to define the shear and bulk viscosities, ν and ζ, in order to incorporate the angular momentum flux reversal effect into the axisymmetric diffusion model for circumferential gaps presented in this thesis (Grätz et al., 2019). The sharp Encke and Keeler gap edges are modeled and conclusions regarding the shear and bulk viscosities of the ring are discussed. Finally, we explore the question of whether the radial density profile of the central and outer A ring, recently measured by Tiscareno and Harris (2018) in the highest resolution to date, and in particular, the sharp outer A ring edge can be modeled consistently from the balance of gravitational scattering by several outer moons and the mass and momentum transport. To this aim, the developed model is extended to account for the inward drifts caused by multiple discrete and overlapping resonances with multiple outer satellites and is then used to hydrodynamically simulate the normalized surface mass density profile of the A ring. This section of the thesis is based on studies by Tajeddine et al. (2017a) who recently discussed the common misconception that the 7:6 resonance with Janus alone maintains the outer A ring edge, showing that the combined effort of several resonances with several outer moons is required to confine the A ring as observed by the Cassini spacecraft.
The passive and active motion of micron-sized tracer particles in crowded liquids and inside living biological cells is ubiquitously characterised by 'viscoelastic' anomalous diffusion, in which the increments of the motion feature long-ranged negative and positive correlations. While viscoelastic anomalous diffusion is typically modelled by a Gaussian process with correlated increments, so-called fractional Gaussian noise, an increasing number of systems are reported, in which viscoelastic anomalous diffusion is paired with non-Gaussian displacement distributions. Following recent advances in Brownian yet non-Gaussian diffusion we here introduce and discuss several possible versions of random-diffusivity models with long-ranged correlations. While all these models show a crossover from non-Gaussian to Gaussian distributions beyond some correlation time, their mean squared displacements exhibit strikingly different behaviours: depending on the model crossovers from anomalous to normal diffusion are observed, as well as a priori unexpected dependencies of the effective diffusion coefficient on the correlation exponent. Our observations of the non-universality of random-diffusivity viscoelastic anomalous diffusion are important for the analysis of experiments and a better understanding of the physical origins of 'viscoelastic yet non-Gaussian' diffusion.