No Arabic abstract
A field theoretical framework is developed for the Hawkes self-excited point process with arbitrary memory kernels by embedding the original non-Markovian one-dimensional dynamics onto a Markovian infinite-dimensional one. The corresponding Langevin dynamics of the field variables is given by stochastic partial differential equations that are Markovian. This is in contrast to the Hawkes process, which is non-Markovian (in general) by construction as a result of its (long) memory kernel. We derive the exact solutions of the Lagrange-Charpit equations for the hyperbolic master equations in the Laplace representation in the steady state, close to the critical point of the Hawkes process. The critical condition of the original Hawkes process is found to correspond to a transcritical bifurcation in the Lagrange-Charpit equations. We predict a power law scaling of the PDF of the intensities in an intermediate asymptotics regime, which crosses over to an asymptotic exponential function beyond a characteristic intensity that diverges as the critical condition is approached. We also discuss the formal relationship between quantum field theories and our formulation. Our field theoretical framework provides a way to tackle complex generalisation of the Hawkes process, such as nonlinear Hawkes processes previously proposed to describe the multifractal properties of earthquake seismicity and of financial volatility.
The Hawkes self-excited point process provides an efficient representation of the bursty intermittent dynamics of many physical, biological, geological and economic systems. By expressing the probability for the next event per unit time (called intensity), say of an earthquake, as a sum over all past events of (possibly) long-memory kernels, the Hawkes model is non-Markovian. By mapping the Hawkes model onto stochastic partial differential equations that are Markovian, we develop a field theoretical approach in terms of probability density functionals. Solving the steady-state equations, we predict a power law scaling of the probability density function (PDF) of the intensities close to the critical point $n=1$ of the Hawkes process, with a non-universal exponent, function of the background intensity $ u_0$ of the Hawkes intensity, the average time scale of the memory kernel and the branching ratio $n$. Our theoretical predictions are confirmed by numerical simulations.
The origin(s) of the ubiquity of Zipfs law (an inverse power law form for the probability density function (PDF) with exponent $1+1$) is still a matter of fascination and investigation in many scientific fields from linguistic, social, economic, computer sciences to essentially all natural sciences. In parallel, self-excited dynamics is a prevalent characteristic of many systems, from seismicity, financial volatility and financial defaults, to sociology, consumer behaviors, computer sciences, The Internet, neuronal discharges and spike trains in biological neuron networks, gene expression and even criminology. Motivated by financial and seismic modelling, we bring the two threads together by introducing a general class of nonlinear self-excited point processes with fast-accelerating intensities as a function of tension. Solving the corresponding master equations, we find that a wide class of such nonlinear Hawkes processes have the PDF of their intensities described by Zipfs law on the condition that (i) the intensity is a fast-accelerating function of tension and (ii) the distribution of the point fertilities is symmetric. This unearths a novel mechanism for Zipfs law, providing a new understanding of its ubiquity.
We study the coarse-graining approach to derive a generator for the evolution of an open quantum system over a finite time interval. The approach does not require a secular approximation but nevertheless generally leads to a Lindblad-Gorini-Kossakowski-Sudarshan generator. By combining the formalism with Full Counting Statistics, we can demonstrate a consistent thermodynamic framework, once the switching work required for the coupling and decoupling with the reservoir is included. Particularly, we can write the second law in standard form, with the only difference that heat currents must be defined with respect to the reservoir. We exemplify our findings with simple but pedagogical examples.
Master equations are common descriptions of mesoscopic systems. Analytical solutions to these equations can rarely be obtained. We here derive an analytical approximation of the time-dependent probability distribution of the master equation using orthogonal polynomials. The solution is given in two alternative formulations: a series with continuous and a series with discrete support both of which can be systematically truncated. While both approximations satisfy the system size expansion of the master equation, the continuous distribution approximations become increasingly negative and tend to oscillations with increasing truncation order. In contrast, the discrete approximations rapidly converge to the underlying non-Gaussian distributions. The theory is shown to lead to particularly simple analytical expressions for the probability distributions of molecule numbers in metabolic reactions and gene expression systems.
Traditionally, Hawkes processes are used to model time--continuous point processes with history dependence. Here we propose an extended model where the self--effects are of both excitatory and inhibitory type and follow a Gaussian Process. Whereas previous work either relies on a less flexible parameterization of the model, or requires a large amount of data, our formulation allows for both a flexible model and learning when data are scarce. We continue the line of work of Bayesian inference for Hawkes processes, and our approach dispenses with the necessity of estimating a branching structure for the posterior, as we perform inference on an aggregated sum of Gaussian Processes. Efficient approximate Bayesian inference is achieved via data augmentation, and we describe a mean--field variational inference approach to learn the model parameters. To demonstrate the flexibility of the model we apply our methodology on data from three different domains and compare it to previously reported results.