Researcher profile

Andrea Zoia

Andrea Zoia contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2022arXiv

Space and time correlations for diffusion models with prompt and delayed birth-and-death events

Understanding the statistical properties of a collection of individuals subject to random displacements and birth-and-death events is key to several applications in physics and life sciences, encompassing the diagnostic of nuclear reactors and the analysis of epidemic patterns. Previous investigations of the critical regime, where births and deaths balance on average, have shown that highly non-Poissonian fluctuations might occur in the population, leading to spontaneous spatial clustering, and eventually to a critical catastrophe, where fluctuations can result in the extinction of the population. A milder behaviour is observed when the population size is kept constant: thefluctuations asymptotically level off and the critical catastrophe is averted. In this paper, we shall extend these results by considering the broader class of models with prompt and delayed birth-and-death events, which mimic the presence of precursors in nuclear reactor physics or incubation inepidemics. We shall consider models with and without population control mechanisms. Analytical or semi-analytical results for the density, the two-point correlation function and the mean-squared pair distance will be derived and compared to Monte Carlo simulations, which will be used as a reference.

preprint2022arXiv

Unbiasedness and Optimization of Regional Weight Cancellation

The Monte Carlo method is often used to simulate systems which can be modeled by random walks. In order to calculate observables, in many implementations the "walkers" carry a statistical weight which is generally assumed to be positive. Some random walk simulations, however, may require walkers to have positive or negative weights: it has been shown that the presence of a mixture of positive and negative weights can impede the statistical convergence, and special weight-cancellation techniques must be adopted in order to overcome these issues. In a recent work we demonstrated the usefulness of one such method, exact regional weight cancellation, to solve eigenvalue problems in nuclear reactor physics in three spatial dimensions. The method previously exhibited had several limitations (including multi-group transport and isotropic scattering) and needed homogeneous cuboid cancellation regions. In this paper we lift the previous limitations, in view of applying exact regional cancellation to more realistic continuous-energy neutron transport problems. This extended regional cancellation framework is used to optimize the efficiency of the weight cancellation. Our findings are illustrated on a benchmark configuration for reactor physics.

preprint2019arXiv

Strongly constrained stochastic processes: the multi-ends Brownian bridge

In a recent article, Krapivsky and Redner (J. Stat. Mech. 093208 (2018)) established that the distribution of the first hitting times for a diffusing particle subject to hitting an absorber is independent of the direction of the external flow field. In the present paper, we build upon this observation and investigate when the conditioning on the diffusion leads to a process that is totally independent of the flow field. For this purpose, we adopt the Langevin approach, or more formally the theory of conditioned stochastic differential equations. This technique allows us to derive a large variety of stochastic processes: in particular, we introduce a new kind of Brownian bridge ending at two different final points and calculate its fundamental probabilities. This method is also very well suited for generating statistically independent paths. Numerical simulations illustrate our findings.

preprint2013arXiv

Spatial extent of an outbreak in animal epidemics

Characterizing the spatial extent of epidemics at the outbreak stage is key to controlling the evolution of the disease. At the outbreak, the number of infected individuals is typically small, so that fluctuations around their average are important: then, it is commonly assumed that the susceptible-infected-recovered (SIR) mechanism can be described by a stochastic birth-death process of Galton-Watson type. The displacements of the infected individuals can be modelled by resorting to Brownian motion, which is applicable when long-range movements and complex network interactions can be safely neglected, as in case of animal epidemics. In this context, the spatial extent of an epidemic can be assessed by computing the convex hull enclosing the infected individuals at a given time. We derive the exact evolution equations for the mean perimeter and the mean area of the convex hull, and compare them with Monte Carlo simulations.

preprint2012arXiv

Branching exponential flights: travelled lengths and collision statistics

The evolution of several physical and biological systems, ranging from neutron transport in multiplying media to epidemics or population dynamics, can be described in terms of branching exponential flights, a stochastic process which couples a Galton-Watson birth-death mechanism with random spatial displacements. Within this context, one is often called to assess the length $\ell_V$ that the process travels in a given region $V$ of the phase space, or the number of visits $n_V$ to this same region. In this paper, we address this issue by resorting to the Feynman-Kac formalism, which allows characterizing the full distribution of $\ell_V$ and $n_V$ and in particular deriving explicit moment formulas. Some other significant physical observables associated to $\ell_V $ and $n_V$, such as the survival probability, are discussed as well, and results are illustrated by revisiting the classical example of the rod model in nuclear reactor physics.

preprint2012arXiv

Discrete Feynman-Kac formulas for branching random walks

Branching random walks are key to the description of several physical and biological systems, such as neutron multiplication, genetics and population dynamics. For a broad class of such processes, in this Letter we derive the discrete Feynman-Kac equations for the probability and the moments of the number of visits $n_V$ of the walker to a given region $V$ in the phase space. Feynman-Kac formulas for the residence times of Markovian processes are recovered in the diffusion limit.

preprint2012arXiv

Properties of branching exponential flights in bounded domains

Branching random flights are key to describing the evolution of many physical and biological systems, ranging from neutron multiplication to gene mutations. When their paths evolve in bounded regions, we establish a relation between the properties of trajectories starting on the boundary and those starting inside the domain. Within this context, we show that the total length travelled by the walker and the number of performed collisions in bounded volumes can be assessed by resorting to the Feynman-Kac formalism. Other physical observables related to the branching trajectories, such as the survival and escape probability, are derived as well.

preprint2011arXiv

Collision densities and mean residence times for $d$-dimensional exponential flights

In this paper we analyze some aspects of {\em exponential flights}, a stochastic process that governs the evolution of many random transport phenomena, such as neutron propagation, chemical/biological species migration, or electron motion. We introduce a general framework for $d$-dimensional setups, and emphasize that exponential flights represent a deceivingly simple system, where in most cases closed-form formulas can hardly be obtained. We derive a number of novel exact (where possible) or asymptotic results, among which the stationary probability density for 2d systems, a long-standing issue in Physics, and the mean residence time in a given volume. Bounded or unbounded, as well as scattering or absorbing domains are examined, and Monte Carlo simulations are performed so as to support our findings.

preprint2011arXiv

Collision number statistics for transport processes

Many physical observables can be represented as a particle spending some random time within a given domain. For a broad class of transport-dominated processes, we detail how it is possible to express the moments of the number of particle collisions in an arbitrary volume in terms of repeated convolutions of the ensemble equilibrium distribution. This approach is shown to generalize the celebrated Kac formula for the moments of residence times, which is recovered in the diffusion limit. Some practical applications are illustrated for bounded, unbounded and absorbing domains.

preprint2011arXiv

Collision statistics for random flights with anisotropic scattering and absorption

For a broad class of random walks with anisotropic scattering kernel and absorption, we derive explicit formulas that allow expressing the moments of the collision number $n_V$ performed in a volume $V$ as a function of the particle equilibrium distribution. Our results apply to arbitrary domains $V$ and boundary conditions, and allow assessing the hitting statistics for systems where the typical displacements are comparable to the domain size, so that the diffusion limit is possibly not attained. An example is discussed for one-dimensional (1d) random flights with exponential displacements, where analytical calculations can be carried out.

preprint2011arXiv

Counting statistics: a Feynman-Kac perspective

By building upon a Feynman-Kac formalism, we assess the distribution of the number of hits in a given region for a broad class of discrete-time random walks with scattering and absorption. We derive the evolution equation for the generating function of the number of hits, and complete our analysis by examining the moments of the distribution, and their relation to the walker equilibrium density. Some significant applications are discussed in detail: in particular, we revisit the gambler's ruin problem and generalize to random walks with absorption the arcsine law for the number of hits on the half-line.

preprint2011arXiv

Residence time and collision statistics for exponential flights: the rod problem revisited

Many random transport phenomena, such as radiation propagation, chemical/biological species migration, or electron motion, can be described in terms of particles performing {\em exponential flights}. For such processes, we sketch a general approach (based on the Feynman-Kac formalism) that is amenable to explicit expressions for the moments of the number of collisions and the residence time that the walker spends in a given volume as a function of the particle equilibrium distribution. We then illustrate the proposed method in the case of the so-called {\em rod problem} (a 1d system), and discuss the relevance of the obtained results in the context of Monte Carlo estimators.

preprint2009arXiv

A continuous time random walk model of transport in variably saturated heterogeneous porous media

We propose a unified physical framework for transport in variably saturated porous media. This approach allows fluid flow and solute migration to be treated as ensemble averages of fluid and solute particles, respectively. We consider the cases of homogeneous and heterogeneous porous materials. Within a fractal mobile-immobile (MIM) continuous time random walk framework, the heterogeneity will be characterized by algebraically decaying particle retention-times. We derive the corresponding (nonlinear) continuum limit partial differential equations and we compare their solutions to Monte Carlo simulation results. The proposed methodology is fairly general and can be used to track fluid and solutes particles trajectories, for a variety of initial and boundary conditions.

preprint2009arXiv

Hitting probability for anomalous diffusion processes

We present the universal features of the hitting probability $Q(x,L)$, the probability that a generic stochastic process starting at $x$ and evolving in a box $[0,L]$ hits the upper boundary $L$ before hitting the lower boundary at 0. For a generic self-affine process (describing, for instance, the polymer translocation through a nanopore) we show that $Q(x,L)=Q(x/L)$ and the scaling function $Q(z)\sim z^ϕ$ as $z\to 0$ with $ϕ=θ/H$ where $H$ and $θ$ are respectively the Hurst exponent and the persistence exponent of the process. This result is verified in several exact calculations including when the process represents the position of a particle diffusing in a disordered potential. We also provide numerical supports for our analytical results.

preprint2009arXiv

Mass transport subject to time-dependent flow with non-uniform sorption in porous media

We address the description of solutes flow with trapping processes in porous media. Starting from a small-scale model for tracer particles trajectories, we derive the corresponding governing equations for the concentration of the mobile and immobile phases. We show that this formulation is fairly general and can easily take into account non-constant coefficients and in particular space-dependent sorption rates. The transport equations are solved numerically and a comparison with Monte Carlo particle-tracking simulations of spatial contaminant profiles and breakthrough curves is proposed, so to illustrate the obtained results.