In the mathematical modeling of many processes and phenomena in the Sciences and Technology one employs systems with many random particles and interactions; see the Application Theme Particlebased modelling in the Sciences. Here we mean many different types of such systems, which, depending on the application, contain moving or static particles, have interactions with each other or with a (possibly random) surrounding medium, have attracting or repelling forces, etc. Also the mathematical questions that we have about the system are quite diverse, e.g., for an emerging cluster structure, for regularities within the clusters, for properties like percolation, crystallization or condensation, for temporal development of the contacts of the particles with each other, for dependencies of the global behavior on parameters, to the point of phase transitions, and more.
In most of the models considered there is one or more order parameters, in terms of which one well expresses the macroscopic behavior, like the average size of a cluster, the occupation probabilities of the particles at a given site at a given time, the empirical mean of all the particles, etc. In some cases, these order parameters satisfy deterministic equations (e.g. differential equations) in the limit considered (often large number of particles or late times), which can afterwards be studied with the help of analytic or numeric methods. In other cases, formulae are derived for the limiting free energy or for laws of large numbers or for ergodic theorems.
For answering these questions and for finding and derivation of formulae, at WIAS mathematical tools and concepts are usually employed and (if not yet available) developed. Stochastic and analytic methods are combined, and in the best case a mathematical solution is derived. Accompanying simulations give visualizations and produce explicit data. Examples of mathematical theories employed are, depending on the model, Gibbs measures, percolation, stochastic (partial) differential equations, Markov processes, and large deviations on the stochastic side, and weak convergence, calculus of variations, convex analysis and partial differential equations on the analytic side.
Contribution of the Institute
Here is a selection of some of the mathematical achievements of the WIAS, see also the Application Theme Particlebased modelling in the Sciences.
A (static) interacting manybody system is given if a large number of points are randomly distributed in a large box, such that they do not accumulate and such that a certain energy is given to the configuration as an exponential probability weight. In this way, a probability measure on configurations is given. The energy term carries a prefactor, the inverse temperature. In the FG5, the thermodynamic limit at low temperatures in particularly large boxes was studied, such that the entropy (the part of the probability that comes from the spatial distribution of the particle cloud) has a particular relation with the energy and the free energy depends on just one parameter. Interesting phase transitions could be obtained. In the future, however, the usual thermodynamic limit will be studied, where the box volume has a fixed relation to the particle number. For a system in which each particle also carries a kinetic energy and the cloud is subject to some symmetry condition, an ansatz with large deviations for a meanfield version of the entire particle system was developed and was used for expressing the free energy in terms of a variational formula. However, the desired effect of a condensation could not be proved in this way, and future investigations will be made.
In telecommunication, the locations of many users are usually modeled as the points of a spatial Poisson point process, see the Application Theme Mobile Communication Networks. Each two of these points interact if their distance is small enough. The question of the probability with which many messages, which are sent through the system via a system of relays, actually reach their intended receiver was studied. That is, a global connectivity property for the system was studied.
The connectivity property was studied in the limit of a high spatial density of the users per volume, in which case the probabilities decay exponentially. With the help of the theory of large deviations for highly dense point processes, the decay rate was expressed in term of an entropy. Afterwards, those configurations were analyzed that minimize the entropy, as these carry the interpretation of the (random) situations of best connectivity under the assumptions made. Similar investigations were made for interference and capacity properties. Further research currently concerns optimal trajectories of the messages, subject to interference and under avoidance of congestion of the relays, as well as the implementation of realistic movement schemes for the users.
A key task in the case of dynamic models is the establishment of a hydrodynamic limit. Such limits are typically first proved for complete applications in order to find an evolution equation for the macroscopic model properties. For examples see the application oriented themes "Coagulation" and Particlebased modelling in the Sciences. A fundamental part of these proofs is the demonstration of compactness in distribution of the Markov Processes that make up the model. In a number of works, properties of application specific problems have been abstracted and the results generalised so that they can be used in further applications.
In Biology the definition of useful stochastic models is an active topic of research that is far from complete. Established models for populations and their movements include spatial branching processes with random motions, which the WIAS studies in random environments; see the mathematical theme Spectral theory of random operators.
Publications
Articles in Refereed Journals

CH. Hirsch, B. Jahnel, Large deviations for the capacity in dynamic spatial relay networks, Markov Processes and Related Fields, 25 (2019), pp. 3373.
Abstract
We derive a large deviation principle for the spacetime evolution of users in a relay network that are unable to connect due to capacity constraints. The users are distributed according to a Poisson point process with increasing intensity in a bounded domain, whereas the relays are positioned deterministically with given limiting density. The preceding work on capacity for relay networks by the authors describes the highly simplified setting where users can only enter but not leave the system. In the present manuscript we study the more realistic situation where users leave the system after a random transmission time. For this we extend the point process techniques developed in the preceding work thereby showing that they are not limited to settings with strong monotonicity properties. 
L. Andreis, A. Asselah, P. Dai Pra , Ergodicity of a system of interacting random walks with asymmetric interaction, Annales de l'Institut Henri Poincare. Probabilites et Statistiques, 55 (2019), pp. 590606.
Abstract
We study N interacting random walks on the positive integers. Each particle has drift delta towards infinity, a reflection at the origin, and a drift towards particles with lower positions. This inhomogeneous mean field system is shown to be ergodic only when the interaction is strong enough. We focus on this latter regime, and point out the effect of piles of particles, a phenomenon absent in models of interacting diffusion in continuous space. 
L. Andreis, P. Dai Pra, M. Fischer, McKeanVlasov limit for interacting systems with simultaneous jumps, Stochastic Analysis and Applications, 36 (2018), pp. 960995, DOI 10.1080/07362994.2018.1486202 .
Abstract
Motivated by several applications, including neuronal models, we consider the McKeanVlasov limit for meanfield systems of interacting diffusions with simultaneous jumps. We prove propagation of chaos via a coupling technique that involves an intermediate process and that gives a rate of convergence for the W1 Wasserstein distance between the empirical measures of the two systems on the space of trajectories D([0,T],R^d). 
O. Gün, A. Yilmaz, The stochastic encountermating model, Acta Applicandae Mathematicae. An International Survey Journal on Applying Mathematics and Mathematical Applications, 148 (2017), pp. 71102.

A. Mielke, R.I.A. Patterson, M.A. Peletier, D.R.M. Renger, Nonequilibrium thermodynamical principles for chemical reactions with massaction kinetics, SIAM Journal on Applied Mathematics, 77 (2017), pp. 15621585, DOI 10.1137/16M1102240 .
Abstract
We study stochastic interacting particle systems that model chemical reaction networks on the micro scale, converging to the macroscopic Reaction Rate Equation. One abstraction level higher, we study the ensemble of such particle systems, converging to the corresponding Liouville transport equation. For both systems, we calculate the corresponding large deviations and show that under the condition of detailed balance, the large deviations induce a nonlinear relation between thermodynamic fluxes and free energy driving force. 
R.I.A. Patterson, S. Simonella, W. Wagner, A kinetic equation for the distribution of interaction clusters in rarefied gases, Journal of Statistical Physics, 169 (2017), pp. 126167.

J. Blath, A. González Casanova Soberón, B. Eldon, N. Kurt, M. WilkeBerenguer, Genetic variability under the seedbank coalescent, Genetics, 200 (2015), pp. 921934.
Abstract
We analyze patterns of genetic variability of populations in the presence of a large seedbank with the help of a new coalescent structure called the seedbank coalescent. This ancestral process appears naturally as a scaling limit of the genealogy of large populations that sustain seedbanks, if the seedbank size and individual dormancy times are of the same order as those of the active population. Mutations appear as Poisson processes on the active lineages and potentially at reduced rate also on the dormant lineages. The presence of "dormant" lineages leads to qualitatively altered times to the most recent common ancestor and nonclassical patterns of genetic diversity. To illustrate this we provide a WrightFisher model with a seedbank component and mutation, motivated from recent models of microbial dormancy, whose genealogy can be described by the seedbank coalescent. Based on our coalescent model, we derive recursions for the expectation and variance of the time to most recent common ancestor, number of segregating sites, pairwise differences, and singletons. Estimates (obtained by simulations) of the distributions of commonly employed distance statistics, in the presence and absence of a seedbank, are compared. The effect of a seedbank on the expected sitefrequency spectrum is also investigated using simulations. Our results indicate that the presence of a large seedbank considerably alters the distribution of some distance statistics, as well as the sitefrequency spectrum. Thus, one should be able to detect from genetic data the presence of a large seedbank in natural populations.
Contributions to Collected Editions

B. Jahnel, W. König, Probabilistic methods for spatial multihop communication systems, in: Topics in Applied Analysis and Optimisation, M. Hintermüller, J.F. Rodrigues, eds., CIM Series in Mathematical Sciences, Springer, Springer International Publishing, 2019, pp. 239268.
Preprints, Reports, Technical Reports

B. Jahnel, A. Tóbiás, E. Cali, Phase transitions for the Boolean model of continuum percolation for Cox point processes, Preprint no. 2704, WIAS, Berlin, 2020, DOI 10.20347/WIAS.PREPRINT.2704 .
Abstract, PDF (389 kByte)
We consider the Boolean model with random radii based on Cox point processes. Under a condition of stabilization for the random environment, we establish existence and nonexistence of subcritical regimes for the size of the cluster at the origin in terms of volume, diameter and number of points. Further, we prove uniqueness of the infinite cluster for sufficiently connected environments. 
A. Hinsen, B. Jahnel, E. Cali, J.P. Wary, Malware propagation in urban D2D networks, Preprint no. 2674, WIAS, Berlin, 2020, DOI 10.20347/WIAS.PREPRINT.2674 .
Abstract, PDF (3133 kByte)
We introduce and analyze models for the propagation of malware in pure D2D networks given via stationary CoxGilbert graphs. Here, the devices form a Poisson point process with random intensity measure λ, Λ where Λ is stationary and given, for example, by the edgelength measure of a realization of a PoissonVoronoi tessellation that represents an urban street system. We assume that, at initial time, a typical device at the center of the network carries a malware and starts to infect neighboring devices after random waiting times. Here we focus on Markovian models, where the waiting times are exponential random variables, and nonMarkovian models, where the waiting times feature strictly positive minimal and finite maximal waiting times. We present numerical results for the speed of propagation depending on the system parameters. In a second step, we introduce and analyze a counter measure for the malware propagation given by special devices called white knights, which have the ability, once attacked, to eliminate the malware from infected devices and turn them into white knights. Based on simulations, we isolate parameter regimes in which the malware survives or is eliminated, both in the Markovian and nonMarkovian setting. 
B. Jahnel, A. Tóbiás, SINR percolation for Cox point processes with random powers, Preprint no. 2659, WIAS, Berlin, 2019, DOI 10.20347/WIAS.PREPRINT.2659 .
Abstract, PDF (356 kByte)
Signaltointerference plus noise ratio (SINR) percolation is an infiniterange dependent variant of continuum percolation modeling connections in a telecommunication network. Unlike in earlier works, in the present paper the transmitted signal powers of the devices of the network are assumed random, i.i.d. and possibly unbounded. Additionally, we assume that the devices form a stationary Cox point process, i.e., a Poisson point process with stationary random intensity measure, in two or higher dimensions. We present the following main results. First, under suitable moment conditions on the signal powers and the intensity measure, there is percolation in the SINR graph given that the device density is high and interferences are sufficiently reduced, but not vanishing. Second, if the interference cancellation factor γ and the SINR threshold τ satisfy γ ≥ 1/(2τ), then there is no percolation for any intensity parameter. Third, in the case of a Poisson point process with constant powers, for any intensity parameter that is supercritical for the underlying Gilbert graph, the SINR graph also percolates with some small but positive interference cancellation factor. 
A. Mielke, A. Stephan, Coarsegraining via EDPconvergence for linear fastslow reaction systems, Preprint no. 2643, WIAS, Berlin, 2019, DOI 10.20347/WIAS.PREPRINT.2643 .
Abstract, PDF (426 kByte)
We consider linear reaction systems with slow and fast reactions, which can be interpreted as master equations or Kolmogorov forward equations for Markov processes on a finite state space. We investigate their limit behavior if the fast reaction rates tend to infinity, which leads to a coarsegrained model where the fast reactions create microscopically equilibrated clusters, while the exchange mass between the clusters occurs on the slow time scale. Assuming detailed balance the reaction system can be written as a gradient flow with respect to the relative entropy. Focusing on the physically relevant coshtype gradient structure we show how an effective limit gradient structure can be rigorously derived and that the coarsegrained equation again has a coshtype gradient structure. We obtain the strongest version of convergence in the sense of the EnergyDissipation Principle (EDP), namely EDPconvergence with tilting. 
A. Hinsen, B. Jahnel, E. Cali, J.P. Wary, Phase transitions for chaseescape models on Gilbert graphs, Preprint no. 2642, WIAS, Berlin, 2019, DOI 10.20347/WIAS.PREPRINT.2642 .
Abstract, PDF (219 kByte)
We present results on phase transitions of local and global survival in a twospecies model on Gilbert graphs. At initial time there is an infection at the origin that propagates on the Gilbert graph according to a continuoustime nearestneighbor interacting particle system. The Gilbert graph consists of susceptible nodes and nodes of a second type, which we call white knights. The infection can spread on susceptible nodes without restriction. If the infection reaches a white knight, this white knight starts to spread on the set of infected nodes according to the same mechanism, with a potentially different rate, giving rise to a competition of chase and escape. We show welldefinedness of the model, isolate regimes of global survival and extinction of the infection and present estimates on local survival. The proofs rest on comparisons to the process on trees, percolation arguments and finitedegree approximations of the underlying random graphs. 
D. Heydecker, R.I.A. Patterson, Bilinear coagulation equations, Preprint no. 2637, WIAS, Berlin, 2019, DOI 10.20347/WIAS.PREPRINT.2637 .
Abstract, PDF (453 kByte)
We consider coagulation equations of Smoluchowski or Flory type where the total merge rate has a bilinear form π(y) · Aπ (x) for a vector of conserved quantities π, generalising the multiplicative kernel. For these kernels, a gelation transition occurs at a finite time t_{g} ∈ (0,∞), which can be given exactly in terms of an eigenvalue problem in finite dimensions. We prove a hydrodynamic limit for a stochastic coagulant, including a corresponding phase transition for the largest particle, and exploit a coupling to random graphs to extend analysis of the limiting process beyond the gelation time. 
CH. Hirsch, B. Jahnel, A. Tóbiás, Lower large deviations for geometric functionals, Preprint no. 2632, WIAS, Berlin, 2019, DOI 10.20347/WIAS.PREPRINT.2632 .
Abstract, PDF (1268 kByte)
This work develops a methodology for analyzing largedeviation lower tails associated with geometric functionals computed on a homogeneous Poisson point process. The technique applies to characteristics expressed in terms of stabilizing score functions exhibiting suitable monotonicity properties. We apply our results to clique counts in the random geometric graph, intrinsic volumes of PoissonVoronoi cells, as well as powerweighted edge lengths in the random geometric, κnearest neighbor and relative neighborhood graph. 
A. Stephan, Combinatorial considerations on the invariant measure of a stochastic matrix, Preprint no. 2627, WIAS, Berlin, 2019, DOI 10.20347/WIAS.PREPRINT.2627 .
Abstract, PDF (225 kByte)
The invariant measure is a fundamental object in the theory of Markov processes. In finite dimensions a Markov process is defined by transition rates of the corresponding stochastic matrix. The Markov tree theorem provides an explicit representation of the invariant measure of a stochastic matrix. In this note, we given a simple and purely combinatorial proof of the Markov tree theorem. In the symmetric case of detailed balance, the statement and the proof simplifies even more. 
S. Jansen, W. König, B. Schmidt, F. Theil, Surface energy and boundary layers for a chain of atoms at low temperature, Preprint no. 2589, WIAS, Berlin, 2019, DOI 10.20347/WIAS.PREPRINT.2589 .
Abstract, PDF (529 kByte)
We analyze the surface energy and boundary layers for a chain of atoms at low temperature for an interaction potential of LennardJones type. The pressure (stress) is assumed small but positive and bounded away from zero, while the temperature goes to zero. Our main results are: (1) As the temperature goes to zero and at fixed positive pressure, the Gibbs measures for infinite chains and semiinfinite chains satisfy path large deviations principles. The rate functions are bulk and surface energy functionals. The minimizer of the surface functional corresponds to zero temperature boundary layers. (2) The surface correction to the Gibbs free energy converges to the zero temperature surface energy, characterized with the help of the minimum of the surface energy functional. (3) The bulk Gibbs measure and Gibbs free energy can be approximated by their Gaussian counterparts. (4) Bounds on the decay of correlations are provided, some of them uniform in the inverse temperature. 
A. Tóbiás, B. Jahnel, Exponential moments for planar tessellations, Preprint no. 2572, WIAS, Berlin, 2019, DOI 10.20347/WIAS.PREPRINT.2572 .
Abstract, PDF (276 kByte)
In this paper we show existence of all exponential moments for the total edge length in a unit disc for a family of planar tessellations based on Poisson point processes. Apart from classical such tessellations like the PoissonVoronoi, PoissonDelaunay and Poisson line tessellation, we also treat the JohnsonMehl tessellation, Manhattan grids, nested versions and Palm versions. As part of our proofs, for some planar tessellations, we also derive existence of exponential moments for the number of cells and the number of edges intersecting the unit disk. 
L. Andreis, W. König, R.I.A. Patterson, A largedeviations approach to gelation, Preprint no. 2568, WIAS, Berlin, 2019, DOI 10.20347/WIAS.PREPRINT.2568 .
Abstract, PDF (338 kByte)
A largedeviations principle (LDP) is derived for the state, at fixed time, of the multiplicative coalescent in the large particle number limit. The rate function is explicit and describes each of the three parts of the state: microscopic, mesoscopic and macroscopic. In particular, it clearly captures the well known gelation phase transition given by the formation of a particle containing a positive fraction of the system mass at time t=1. Via a standard map of the multiplicative coalescent onto a timedependent version of the ErdősRényi random graph, our results can also be rephrased as an LDP for the component sizes in that graph. Our proofs rely on estimates and asymptotics for the probability that smaller ErdősRényi graphs are connected. 
M. Mittnenzweig, Hydrodynamic limit and large deviations of reactiondiffusion master equations, Preprint no. 2521, WIAS, Berlin, 2018, DOI 10.20347/WIAS.PREPRINT.2521 .
Abstract, PDF (389 kByte)
We derive the hydrodynamic limit of a reactiondiffusion master equation, that combines an exclusion process with a reversible chemical master equation expression for the reaction rates. The crucial assumption is that the associated macroscopic reaction network has a detailed balance equilibrium. The hydrodynamic limit is given by a system of reactiondiffusion equations with a modified mass action law for the reaction rates. We provide the upper bound for large deviations of the empirical measure from the hydrodynamic limit.
Talks, Poster

A. Stephan, Coarsegraining via EDPconvergence for linear fastslow reaction systems, Seminar ``Applied Analysis'', Eindhoven University of Technology, Centre for Analysis, Scientific Computing, and Applications  Mathematics and Computer Science, Netherlands, January 20, 2020.

A. Stephan, On mathematical coarsegraining for linear reaction systems, 8th BMS Student Conference, February 19  21, 2020, Technische Universität Berlin, February 21, 2020.

A. Hinsen, Data mobility in adhoc networks: Vulnerability and security, KEIN öffentlicher Vortrag (Orange), Telecom Orange Paris, France, December 12, 2019.

A. Stephan, Rigorous derivation of the effective equation of a linear reaction system with different time scales, 90th Annual Meeting of the International Association of Applied Mathematics and Mechanics (GAMM 2019), Section S14 ``Applied Analysis'', February 18  22, 2019, Universität Wien, Technische Universität Wien, Austria, February 21, 2019.

B. Jahnel, Continuum percolation in random environment I  IV, Workshop on Probability, Analysis and Applications (PAA), September 23  October 4, 2019, African Institute for Mathematical Sciences  Ghana (AIMS Ghana), Accra, Ghana.

R.I.A. Patterson, Interaction clusters for the Kac process, BerlinLeipzig workshop in analysis and stochastics, January 16  18, 2019, Max Planck Institute ür Mathematik in den Naturwissenschaften, Leipzig, January 18, 2019.

R.I.A. Patterson, Fluctuations and confidence intervals for stochastic particle simulations, First Berlin  Leipzig Workshop on Fluctuating Hydrodynamics, August 26  30, 2019, Max Planck Institute for Mathematics in the Sciences, Leipzig, August 29, 2019.

R.I.A. Patterson, Flux large deviations, Workshop on Chemical Reaction Networks, July 1  3, 2019, Politecnico di Torino, Dipartimento di Scienze Matematiche ''G. L. Lagrange``, Italy, July 2, 2019.

R.I.A. Patterson, Flux large deviations, Seminar Statistical Laboratory, University of Cambridge, Faculty of Mathematics, UK, May 7, 2019.

R.I.A. Patterson, Interaction clusters for the Kac process, Workshopon Effective equations: frontiers in classical and quantum systems, June 24  28, 2019, Hausdorff Research Institute for Mathematics, Bonn, June 28, 2019.

R.I.A. Patterson, The role of fluctuating hydrodynamics in the CRC 1114, CRC 1114 School 2019: Fluctuating Hydrodynamics, Zuse Institute Berlin (ZIB), October 28, 2019.

R.I.A. Patterson, Kinetic interaction clusters, Oberseminar, MartinLutherUniversität HalleWittenberg, Naturwissenschaftliche Fakultät II  Chemie, Physik und Mathematik, April 17, 2019.

L. Taggi, Critical density in activated random walks, Horowitz Seminar on Probability, Ergodic Theory and Dynamical Systems, Tel Aviv University, School of Mathematical Sciences, Israel, May 20, 2019.

M. Maurelli, McKeanVlasov SDEs with irregular drift: Large deviations for particle approximation, University of Oxford, Mathematical Institute, UK, March 5, 2018.

M. Maurelli , A McKeanVlasov SDE with reflecting boundaries, CASA Colloquium, Eindhoven University of Technology, Department of Mathematics and Computer Science, Netherlands, January 10, 2018.

L. Andreis, A largedeviations approach to the multiplicative coagulation process, Probability Seminar, Università degli Studi di Padova, Dipartimento di Matematica ``Tullio LeviCivita'', Italy, October 12, 2018.

L. Andreis, A largedeviations approach to the multiplicative coagulation process, Seminar ''Theory of Complex Systems and Neurophysics  Theory of Statistical Physics and Nonlinear Dynamics``, HumboldtUniversität zu Berlin, Institut für Physik, October 30, 2018.

L. Andreis, Ergodicity of a system of interacting random walks with asymmetric interaction, 13th German Probability and Statistics Days 2018  Freiburger StochastikTage, February 27  March 2, 2018, AlbertLudwigsUniversität Freiburg, Abteilung für Mathematische Stochastik, Freiburg, February 1, 2018.

L. Andreis, Networks of interacting components with macroscopic selfsustained periodic behavior, Neural Coding 2018, September 9  14, 2018, University of Torino, Department of Mathematics, Italy, September 10, 2018.

L. Andreis, Selfsustained periodic behavior in interacting systems, Random Structures in Neuroscience and Biology, March 26  29, 2018, LudwigMaximilians Universität München, Fakultät für Mathematik, Informatik und Statistik, Herrsching, March 26, 2018.

L. Andreis, System of interacting random walks with asymmetric interaction, 48th Probability Summer School, July 8  20, 2018, Clermont Auvergne University, Saint Flour, France, July 17, 2018.

W. Dreyer, Thermodynamics and kinetic theory of nonNewtonian fluids, Technische Universität Darmstadt, Mathematische Modellierung und Analysis, June 13, 2018.

W. Dreyer, J. Fuhrmann, P. Gajewski, C. Guhlke, M. Landstorfer, M. Maurelli, R. Müller, Stochastic model for LiFePO4electrodes, ModVal14  14th Symposium on Fuel Cell and Battery Modeling and Experimental Validation, Karlsruhe, March 2  3, 2017.

D.R.M. Renger, From large deviations to Wasserstein gradient flows in multiple dimensions, Workshop on Gradient Flows, Large Deviations and Applications, November 22  29, 2015, EURANDOM, Mathematics and Computer Science Department, Eindhoven, Netherlands, November 23, 2015.

D.R.M. Renger, The inverse problem: From gradient flows to large deviations, Workshop ``Analytic Approaches to Scaling Limits for Random System'', January 26  30, 2015, Universität Bonn, Hausdorff Research Institute for Mathematics, January 26, 2015.
External Preprints

D. Heydecker , R.I.A. Patterson, Kac interaction clusters: A bilinear coagulation equation and phase transition, Preprint no. arXiv:1902.07686, Cornell University Library, 2019.
Abstract
We consider the interaction clusters for Kac's model of a gas with quadratic interaction rates, and show that they behave as coagulating particles with a bilinear coagulation kernel. In the large particle number limit the distribution of the interaction cluster sizes is shown to follow an equation of Smoluchowski type. Using a coupling to random graphs, we analyse the limiting equation, showing wellposedness, and a closed form for the time of the gelation phase transition tg when a macroscopic cluster suddenly emerges. We further prove that the second moment of the cluster size distribution diverges exactly at tg. Our methods apply immediately to coagulating particle systems with other bilinear coagulation kernels. 
A. González Casanova Soberón, J.C. Pardo, J.L. Perez, Branching processes with interactions: The subcritical cooperative regime, Preprint no. arXiv:1704.04203, Cornell University Library, arXiv.org, 2017.
Abstract
In this paper, we introduce a particular family of processes with values on the nonnegative integers that model the dynamics of populations where individuals are allow to have different types of inter actions. The types of interactions that we consider include pairwise: competition, annihilation and cooperation; and interaction among several individuals that can be consider as catastrophes. We call such families of processes branching processes with interactions. In particular, we prove that a process in this class has a moment dual which turns out to be a jumpdiffusion that can be thought as the evolution of the frequency of a trait or phenotype. The aim of this paper is to study the long term behaviour of branching processes with interac tions under the assumption that the cooperation parameter satisfies a given condition that we called subcritical cooperative regime. The moment duality property is useful for our purposes. 
J. Blath, E. Buzzoni, A. Casanova Soberón, M.W. Berenguer, The seed bank diffusion, and its relation to the twoisland model, Preprint no. arXiv:1710.08164, Cornell University Library, arXiv.org, 2017.
Abstract
In this paper, we introduce a particular family of processes with values on the nonnegative integers that model the dynamics of populations where individuals are allow to have different types of inter actions. The types of interactions that we consider include pairwise: competition, annihilation and cooperation; and interaction among several individuals that can be consider as catastrophes. We call such families of processes branching processes with interactions. In particular, we prove that a process in this class has a moment dual which turns out to be a jumpdiffusion that can be thought as the evolution of the frequency of a trait or phenotype. The aim of this paper is to study the long term behaviour of branching processes with interac tions under the assumption that the cooperation parameter satisfies a given condition that we called subcritical cooperative regime. The moment duality property is useful for our purposes. 
K.F. Lee, M. Dosta, A.D. Mcguire, S. Mosbach, W. Wagner, S. Heinrich, M. Kraft, Development of a multicompartment population balance model for highshear wet granulation with Discrete Element Method, Technical report no. 170, c4ePreprint Series, 2016.
Abstract
This paper presents a multicompartment population balance model for wet granulation coupled with DEM (Discrete Element Method) simulations. Methodologies are developed to extract relevant data from the DEM simulations to inform the population balance model. First, compartmental residence times are calculated for the population balance model from DEM. Then, a suitable collision kernel is chosen for the population balance model based on particleparticle collision frequencies extracted from DEM. It is found t hat the population balance model is able to predict the trends exhibited by the experimental size and porosity distributions by utilising the information provided by the DEM simulations.