BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Enrico Scalas

Enrico Scalas appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.PRcond-mat.stat-mechq-fin.STphysics.data-anphysics.soc-phq-fin.CPq-fin.TRcond-mat.dis-nnmath-phmath.MPMathematical SoftwareNumerical Analysisphysics.geo-phq-fin.PMq-fin.PR

Catalog footprint

What is connected

21works
15topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

21 published item(s)

preprint2022arXiv

A fractional approach to study the pure-temporal Epidemic Type Aftershock Sequence (ETAS) process for earthquakes modeling

In statistical seismology, the Epidemic Type Aftershocks Sequence (ETAS) model is a branching process used world-wide to forecast earthquake intensity rates and reproduce many statistical features observed in seismicity catalogs. In this paper, we describe a fractional differential equation that governs the earthquake intensity rate of the pure temporal ETAS model by using the Caputo fractional derivative and we solve it analytically. We highlight that the tools and special functions of fractional calculus simplify the classical methods employed to obtain the intensity rate and let us describe the change of solution decay for large times. We also apply and discuss the theoretical results to the Japanese catalog in the period 1965-2003.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A fractional generalization of the Dirichlet distribution and related distributions

This paper is devoted to a fractional generalization of the Dirichlet distribution. The form of the multivariate distribution is derived assuming that the $n$ partitions of the interval $[0,W_n]$ are independent and identically distributed random variables following the generalized Mittag-Leffler distribution. The expected value and variance of the one-dimensional marginal are derived as well as the form of its probability density function. A related generalized Dirichlet distribution is studied that provides a reasonable approximation for some values of the parameters. The relation between this distribution and other generalizations of the Dirichlet distribution is discussed. Monte Carlo simulations of the one-dimensional marginals for both distributions are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Generalization of the Space-Fractional Poisson Process and its Connection to some Lévy Processes

This paper introduces a generalization of the so-called space-fractional Poisson process by extending the difference operator acting on state space present in the associated difference-differential equations to a much more general form. It turns out that this generalization can be put in relation to a specific subordination of a homogeneous Poisson process by means of a subordinator for which it is possible to express the characterizing Lévy measure explicitly. Moreover, the law of this subordinator solves a one-sided first-order differential equation in which a particular convolution-type integral operator appears, called Prabhakar derivative. In the last section of the paper, a similar model is introduced in which the Prabhakar derivative also acts in time. In this case, too, the probability generating function of the corresponding process and the probability distribution are determined.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Solvable non-Markovian dynamic network

Non-Markovian processes are widespread in natural and human-made systems, yet explicit model- ling and analysis of such systems is underdeveloped. We consider a non-Markovian dynamic network with random link activation and deletion (RLAD) and heavy tailed Mittag-Leffler distribution for the inter-event times. We derive an analytically and computationally tractable system of Kolmogorov- like forward equations utilising the Caputo derivative for the probability of having a given number of active links in the network and solve them. Simulations for the RLAD are also studied for power-law inter-event times and we show excellent agreement with the Mittag-Leffler model. This agreement holds even when the RLAD network dynamics is coupled with the susceptible-infected-susceptible (SIS) spreading dynamics. Thus, the analytically solvable Mittag-Leffler model provides an excel- lent approximation to the case when the network dynamics is characterised by power-law distributed inter-event times. We further discuss possible generalizations of our result.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Velocity and energy distributions in microcanonical ensembles of hard spheres

In a microcanonical ensemble (constant $NVE$, hard reflecting walls) and in a molecular dynamics ensemble (constant $NVE\mathbf{PG}$, periodic boundary conditions) with a number $N$ of smooth elastic hard spheres in a $d$-dimensional volume $V$ having a total energy $E$, a total momentum $\mathbf{P}$, and an overall center of mass position $\mathbf{G}$, the individual velocity components, velocity moduli, and energies have transformed beta distributions with different arguments and shape parameters depending on $d$, $N$, $E$, the boundary conditions, and possible symmetries in the initial conditions. This can be shown marginalizing the joint distribution of individual energies, which is a symmetric Dirichlet distribution. In the thermodynamic limit the beta distributions converge to gamma distributions with different arguments and shape or scale parameters, corresponding respectively to the Gaussian, i.e., Maxwell-Boltzmann, Maxwell, and Boltzmann or Boltzmann-Gibbs distribution. These analytical results agree with molecular dynamics and Monte Carlo simulations with different numbers of hard disks or spheres and hard reflecting walls or periodic boundary conditions. The agreement is perfect with our Monte Carlo algorithm, which acts only on velocities independently of positions with the collision versor sampled uniformly on a unit half sphere in $d$ dimensions, while slight deviations appear with our molecular dynamics simulations for the smallest values of $N$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Ergodic transition in a simple model of the continuous double auction

We study a phenomenological model for the continuous double auction, equivalent to two independent $M/M/1$ queues. The continuous double auction defines a continuous-time random walk for trade prices. The conditions for ergodicity of the auction are derived and, as a consequence, three possible regimes in the behavior of prices and logarithmic returns are observed. In the ergodic regime, prices are unstable and one can observe an intermittent behavior in the logarithmic returns. On the contrary, non-ergodicity triggers stability of prices, even if two different regimes can be seen.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A parsimonious model for intraday European option pricing

A stochastic model for pure-jump diffusion (the compound renewal process) can be used as a zero-order approximation and as a phenomenological description of tick-by-tick price fluctuations. This leads to an exact and explicit general formula for the martingale price of a European call option. A complete derivation of this result is presented by means of elementary probabilistic tools.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the non-stationarity of financial time series: impact on optimal portfolio selection

We investigate the possible drawbacks of employing the standard Pearson estimator to measure correlation coefficients between financial stocks in the presence of non-stationary behavior, and we provide empirical evidence against the well-established common knowledge that using longer price time series provides better, more accurate, correlation estimates. Then, we investigate the possible consequences of instabilities in empirical correlation coefficient measurements on optimal portfolio selection. We rely on previously published works which provide a framework allowing to take into account possible risk underestimations due to the non-optimality of the portfolio weights being used in order to distinguish such non-optimality effects from risk underestimations genuinely due to non-stationarities. We interpret such results in terms of instabilities in some spectral properties of portfolio correlation matrices.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A class of CTRWs: Compound fractional Poisson processes

This chapter is an attempt to present a mathematical theory of compound fractional Poisson processes. The chapter begins with the characterization of a well-known Lévy process: The compound Poisson process. The semi-Markov extension of the compound Poisson process naturally leads to the compound fractional Poisson process, where the Poisson counting process is replaced by the Mittag-Leffler counting process also known as fractional Poisson process. This process is no longer Markovian and Lévy. However, several analytical results are available and some of them are discussed here. The functional limit of the compound Poisson process is an $α$-stable Lévy process, whereas in the case of the compound fractional Poisson process, one gets an $α$-stable Lévy process subordinated to the fractional Poisson process.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Full characterization of the fractional Poisson process

The fractional Poisson process (FPP) is a counting process with independent and identically distributed inter-event times following the Mittag-Leffler distribution. This process is very useful in several fields of applied and theoretical physics including models for anomalous diffusion. Contrary to the well-known Poisson process, the fractional Poisson process does not have stationary and independent increments. It is not a Lévy process and it is not a Markov process. In this letter, we present formulae for its finite-dimensional distribution functions, fully characterizing the process. These exact analytical results are compared to Monte Carlo simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Semi-Markov Graph Dynamics

In this paper, we outline a model of graph (or network) dynamics based on two ingredients. The first ingredient is a Markov chain on the space of possible graphs. The second ingredient is a semi-Markov counting process of renewal type. The model consists in subordinating the Markov chain to the semi-Markov counting process. In simple words, this means that the chain transitions occur at random time instants called epochs. The model is quite rich and its possible connections with algebraic geometry are briefly discussed. Moreover, for the sake of simplicity, we focus on the space of undirected graphs with a fixed number of nodes. However, in an example, we present an interbank market model where it is meaningful to use directed graphs or even weighted graphs.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On-line trading as a renewal process: Waiting time and inspection paradox

We briefly review our recent studies on stochastic processes modelling internet on-line trading. We present a way to evaluate the average waiting time between the observation of the price in financial markets and the next price change, especially in an on-line foreign exchange trading service for individual customers via the internet. The basic method of our approach depends on the so-called renewal-reward theorem. Assuming that the stochastic process modelling the price change is a renewal process, we use the theorem to calculate the average waiting time of the process. The so-called ``inspection paradox'' is discussed, which, in general, means that the average durations is shorter than the average waiting time.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Monte Carlo simulation of uncoupled continuous-time random walks yielding a stochastic solution of the space-time fractional diffusion equation

We present a numerical method for the Monte Carlo simulation of uncoupled continuous-time random walks with a Levy alpha-stable distribution of jumps in space and a Mittag-Leffler distribution of waiting times, and apply it to the stochastic solution of the Cauchy problem for a partial differential equation with fractional derivatives both in space and in time. The one-parameter Mittag-Leffler function is the natural survival probability leading to time-fractional diffusion equations. Transformation methods for Mittag-Leffler random variables were found later than the well-known transformation method by Chambers, Mallows, and Stuck for Levy alpha-stable random variables and so far have not received as much attention; nor have they been used together with the latter in spite of their mathematical relationship due to the geometric stability of the Mittag-Leffler distribution. Combining the two methods, we obtain an accurate approximation of space- and time-fractional diffusion processes almost as easy and fast to compute as for standard diffusion processes.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Random numbers from the tails of probability distributions using the transformation method

The speed of many one-line transformation methods for the production of, for example, Levy alpha-stable random numbers, which generalize Gaussian ones, and Mittag-Leffler random numbers, which generalize exponential ones, is very high and satisfactory for most purposes. However, for the class of decreasing probability densities fast rejection implementations like the Ziggurat by Marsaglia and Tsang promise a significant speed-up if it is possible to complement them with a method that samples the tails of the infinite support. This requires the fast generation of random numbers greater or smaller than a certain value. We present a method to achieve this, and also to generate random numbers within any arbitrary interval. We demonstrate the method showing the properties of the transform maps of the above mentioned distributions as examples of stable and geometric stable random numbers used for the stochastic solution of the space-time fractional diffusion equation.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Spectral densities of Wishart-Levy free stable random matrices: Analytical results and Monte Carlo validation

Random matrix theory is used to assess the significance of weak correlations and is well established for Gaussian statistics. However, many complex systems, with stock markets as a prominent example, exhibit statistics with power-law tails, that can be modelled with Levy stable distributions. We review comprehensively the derivation of an analytical expression for the spectra of covariance matrices approximated by free Levy stable random variables and validate it by Monte Carlo simulation.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Stochastic calculus for uncoupled continuous-time random walks

The continuous-time random walk (CTRW) is a pure-jump stochastic process with several applications in physics, but also in insurance, finance and economics. A definition is given for a class of stochastic integrals driven by a CTRW, that includes the Ito and Stratonovich cases. An uncoupled CTRW with zero-mean jumps is a martingale. It is proved that, as a consequence of the martingale transform theorem, if the CTRW is a martingale, the Ito integral is a martingale too. It is shown how the definition of the stochastic integrals can be used to easily compute them by Monte Carlo simulation. The relations between a CTRW, its quadratic variation, its Stratonovich integral and its Ito integral are highlighted by numerical calculations when the jumps in space of the CTRW have a symmetric Levy alpha-stable distribution and its waiting times have a one-parameter Mittag-Leffler distribution. Remarkably these distributions have fat tails and an unbounded quadratic variation. In the diffusive limit of vanishing scale parameters, the probability density of this kind of CTRW satisfies the space-time fractional diffusion equation (FDE) or more in general the fractional Fokker-Planck equation, that generalize the standard diffusion equation solved by the probability density of the Wiener process, and thus provides a phenomenologic model of anomalous diffusion. We also provide an analytic expression for the quadratic variation of the stochastic process described by the FDE, and check it by Monte Carlo.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The distribution of first-passage times and durations in FOREX and future markets

Possible distributions are discussed for intertrade durations and first-passage processes in financial markets. The view-point of renewal theory is assumed. In order to represent market data with relatively long durations, two types of distributions are used, namely, a distribution derived from the so-called Mittag-Leffler survival function and the Weibull distribution. For Mittag-Leffler type distribution, the average waiting time (residual life time) is strongly dependent on the choice of a cut-off parameter t_ max, whereas the results based on the Weibull distribution do not depend on such a cut-off. Therefore, a Weibull distribution is more convenient than a Mittag-Leffler type one if one wishes to evaluate relevant statistics such as average waiting time in financial markets with long durations. On the other side, we find that the Gini index is rather independent of the cut-off parameter. Based on the above considerations, we propose a good candidate for describing the distribution of first-passage time in a market: The Weibull distribution with a power-law tail. This distribution compensates the gap between theoretical and empirical results much more efficiently than a simple Weibull distribution. We also give a useful formula to determine an optimal crossover point minimizing the difference between the empirical average waiting time and the one predicted from renewal theory. Moreover, we discuss the limitation of our distributions by applying our distribution to the analysis of the BTP future and calculating the average waiting time. We find that our distribution is applicable as long as durations follow a Weibull-law for short times and do not have too heavy a tail.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

The Ehrenfest urn revisited: Playing the game on a realistic fluid model

The Ehrenfest urn process, also known as the dogs and fleas model, is realistically simulated by molecular dynamics of the Lennard-Jones fluid. The key variable is Delta z, i.e. the absolute value of the difference between the number of particles in one half of the simulation box and in the other half. This is a pure-jump stochastic process induced, under coarse graining, by the deterministic time evolution of the atomic coordinates. We discuss the Markov hypothesis by analyzing the statistical properties of the jumps and of the waiting times between jumps. In the limit of a vanishing integration time-step, the distribution of waiting times becomes closer to an exponential and, therefore, the continuous-time jump stochastic process is Markovian. The random variable Delta z behaves as a Markov chain and, in the gas phase, the observed transition probabilities follow the predictions of the Ehrenfest theory.

0 citations0 reviewsNo code linkedOpen access