Search

We derive a Sturm-Liouville system of equations for the exact calculation of the survival probability in first-passage time problems. This system is the one associated with the Wiener-Hopf integral equation obtained from the theory of random walks. The derived approach is an alternative to the existing literature and we tested it against direct calculations from both discrete- and continuous-time random walks in a manageable, but meaningful, example. Within this framework, the Sparre Andersen theorem results to be a boundary condition for the system., Comment: 10 pages

The nature of the first Pop III stars is still a mystery and the energy distribution of the first supernovae is completely unexplored. For the first time we account simultaneously for the unknown initial mass function (IMF), stellar mixing, and energy distribution function (EDF) of Pop III stars in the context of a cosmological model for the formation of a MW-analogue. Our data-calibrated semi-analytic model is based on a N-body simulation and follows the formation and evolution of both Pop III and Pop II/I stars in their proper timescales. We discover degeneracies between the adopted Pop III unknowns, in the predicted metallicity and carbonicity distribution functions and the fraction of C-enhanced stars. Nonetheless, we are able to provide the first available constraints on the EDF, $dN/dE_\star \propto E_{\star}^{-\alpha_e}$ with $1\leq \alpha_e \leq2.5$. In addition, the characteristic mass of the Pop III IMF should be $m_{\rm ch}<100\:{\rm M_\odot}$, assuming a mass range consistent with hydrodynamical simulations (0.1-1000$\:{\rm M_\odot}$). Independent of the assumed Pop III properties, we find that all [C/Fe]>+0.7 stars (with [Fe/H]<-2.8) have been enriched by Pop III supernovae at a $>20\%$ level, and all [C/Fe]>+2 stars at a $>95\%$ level. All very metal-poor stars with $\rm [C/Fe]<0$ are predicted to be predominantly enriched by Pop III hypernovae and/or pair instabillity supernovae. To better constrain the primordial EDF, it is absolutely crucial to have a complete and accurate determination of the metallicity distribution function, and the properties of C-enhanced metal-poor stars (frequency and [C/Fe]) in the Galactic halo., Comment: 22 pages, 13 figures

We study the mean first-passage time (MFPT) for asymmetric continuous-time random walks in continuous-space characterised by waiting-times with finite mean and by jump-sizes with both finite mean and finite variance. In the asymptotic limit, this well-controlled process is governed by an advection-diffusion equation and the MFPT results to be finite when the advecting velocity is in the direction of the boundary. We derive a nonhomogeneous Wiener-Hopf integral equation that allows for the exact calculation of the MFPT by avoiding asymptotic limits and it emerges to depend on the whole distribution of the jump-sizes and on the mean-value only of the waiting-times, thus it holds for general non-Markovian random walks. Through the case study of a quite general family of asymmetric distributions of the jump-sizes that is exponential towards the boundary and arbitrary in the opposite direction, we show that the MFPT is indeed independent of the jump-sizes distribution in the opposite direction to the boundary. Moreover, we show also that there exists a length-scale, which depends only on the features of the distribution of jumps in the direction of the boundary, such that for starting points near the boundary the MFPT depends on the specific whole distribution of jump-sizes, in opposition to the universality emerging for starting points far-away from the boundary., Comment: 20 pages, 3 figures

The various types of generalized Cattaneo, called also telegrapher's equation, are studied. We find conditions under which solutions of the equations considered so far can be recognized as probability distributions, \textit{i.e.} are normalizable and non-negative on their domains. Analysis of the relevant mean squared displacements enables us to classify diffusion processes described by such obtained solutions and to identify them with either ordinary or anomalous super- or subdiffusion. To complete our study we analyse derivations of just considered examples the generalized Cattaneo equations using the continuous time random walk and the persistent random walk approaches.

Fractional Brownian motion is a Gaussian stochastic process with stationary, long-time correlated increments and is frequently used to model anomalous diffusion processes. We study numerically fractional Brownian motion confined to a finite interval with reflecting boundary conditions. The probability density function of this reflected fractional Brownian motion at long times converges to a stationary distribution showing distinct deviations from the fully flat distribution of amplitude $1/L$ in an interval of length $L$ found for reflected normal Brownian motion. While for superdiffusion, corresponding to a mean squared displacement $\langle X^2(t)\rangle\simeq t^{\alpha}$ with $1<\alpha<2$, the probability density function is lowered in the centre of the interval and rises towards the boundaries, for subdiffusion ($0<\alpha<1$) this behaviour is reversed and the particle density is depleted close to the boundaries. The mean squared displacement in these cases at long times converges to a stationary value, which is, remarkably, monotonically increasing with the anomalous diffusion exponent $\alpha$. Our a priori surprising results may have interesting consequences for the application of fractional Brownian motion for processes such as molecule or tracer diffusion in the confined of living biological cells or organelles, or other viscoelastic environments such as dense liquids in microfluidic chambers., Comment: 19 pages, 8 figures, IOP LaTeX

A considerable number of systems have recently been reported in which Brownian yet non-Gaussian dynamics was observed. These are processes characterised by a linear growth in time of the mean squared displacement, yet the probability density function of the particle displacement is distinctly non-Gaussian, and often of exponential (Laplace) shape. This apparently ubiquitous behaviour observed in very different physical systems has been interpreted as resulting from diffusion in inhomogeneous environments and mathematically represented through a variable, stochastic diffusion coefficient. Indeed different models describing a fluctuating diffusivity have been studied. Here we present a new view of the stochastic basis describing time dependent random diffusivities within a broad spectrum of distributions. Concretely, our study is based on the very generic class of the generalised Gamma distribution. Two models for the particle spreading in such random diffusivity settings are studied. The first belongs to the class of generalised grey Brownian motion while the second follows from the idea of diffusing diffusivities. The two processes exhibit significant characteristics which reproduce experimental results from different biological and physical systems. We promote these two physical models for the description of stochastic particle motion in complex environments., Comment: 33 pages, 13 figures, IOP LaTeX

The emerging diffusive dynamics in many complex systems shows a characteristic crossover behaviour from anomalous to normal diffusion which is otherwise fitted by two independent power-laws. A prominent example for a subdiffusive-diffusive crossover are viscoelastic systems such as lipid bilayer membranes, while superdiffusive-diffusive crossovers occur in systems of actively moving biological cells. We here consider the general dynamics of a stochastic particle driven by so-called tempered fractional Gaussian noise, that is noise with Gaussian amplitude and power-law correlations, which are cut off at some mesoscopic time scale. Concretely we consider such noise with built-in exponential or power-law tempering, driving an overdamped Langevin equation (fractional Brownian motion) and fractional Langevin equation motion. We derive explicit expressions for the mean squared displacement and correlation functions, including different shapes of the crossover behaviour depending on the concrete tempering, and discuss the physical meaning of the tempering. In the case of power-law tempering we also find a crossover behaviour from faster to slower superdiffusion and slower to faster subdiffusion. As a direct application of our model we demonstrate that the obtained dynamics quantitatively described the subdiffusion-diffusion and subdiffusion-subdiffusion crossover in lipid bilayer systems. We also show that a model of tempered fractional Brownian motion recently proposed by Sabzikar and Meerschaert leads to physically very different behaviour with a seemingly paradoxical ballistic long time scaling., Comment: 40 pages, 14 figures, IOP LaTeX

A Voigt profile function emerges in several physical investigations (e.g. atmospheric radiative transfer, astrophysical spectroscopy, plasma waves and acoustics) and it turns out to be the convolution of the Gaussian and the Lorentzian densities. Its relation with a number of special functions has been widely derived in literature starting from its Fourier type integral representation. The main aim of the present paper is to introduce the Mellin-Barnes integral representation as a useful tool to obtain new analytical results. Here, starting from the Mellin-Barnes integral representation, the Voigt function is expressed in terms of the Fox H-function which includes representations in terms of the Meijer G-function and previously well-known representations with other special functions., Comment: Submitted to: J. Phys. A: Math. Gen

This paper has been withdrawn by the author due to a crucial error in the formulation., Comment: This paper has been withdrawn

The aim of the article is to investigate the relative dispersion properties of the Well Mixed class of Lagrangian Stochastic Models. Dimensional analysis shows that given a model in the class, its properties depend solely on a non-dimensional parameter, which measures the relative weight of Lagrangian-to-Eulerian scales. This parameter is formulated in terms of Kolmogorov constants, and model properties are then studied by modifying its value in a range that contains the experimental variability. Large variations are found for the quantity $g^*=2gC_0^{-1}$, where $g$ is the Richardson constant, and for the duration of the $t^3$ regime. Asymptotic analysis of model behaviour clarifies some inconsistencies in the literature and excludes the Ornstein-Uhlenbeck process from being considered a reliable model for relative dispersion., Comment: 15 pages, 10 figures, LaTeX source