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

Guido Germano

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

ResearcherUnclaimed source record
cond-mat.stat-mechcond-mat.softcond-mat.mtrl-sciq-fin.CPphysics.soc-phq-fin.GNmath-phmath.MPmath.PRMathematical SoftwareNumerical Analysisphysics.comp-phq-fin.MF

Catalog footprint

What is connected

18works
13topics
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

18 published item(s)

preprint2020arXiv

Hilbert transform, spectral filters and option pricing

We show how spectral filters can improve the convergence of numerical schemes which use discrete Hilbert transforms based on a sinc function expansion, and thus ultimately on the fast Fourier transform. This is relevant, for example, for the computation of fluctuation identities, which give the distribution of the maximum or the minimum of a random path, or the joint distribution at maturity with the extrema staying below or above barriers. We use as examples the methods by Feng and Linetsky (2008) and Fusai, Germano and Marazzina (2016) to price discretely monitored barrier options where the underlying asset price is modelled by an exponential Lévy process. Both methods show exponential convergence with respect to the number of grid points in most cases, but are limited to polynomial convergence under certain conditions. We relate these rates of convergence to the Gibbs phenomenon for Fourier transforms and achieve improved results with spectral filtering.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Full and fast calibration of the Heston stochastic volatility model

This paper presents an algorithm for a complete and efficient calibration of the Heston stochastic volatility model. We express the calibration as a nonlinear least squares problem. We exploit a suitable representation of the Heston characteristic function and modify it to avoid discontinuities caused by branch switchings of complex functions. Using this representation, we obtain the analytical gradient of the price of a vanilla option with respect to the model parameters, which is the key element of all variants of the objective function. The interdependency between the components of the gradient enables an efficient implementation which is around ten times faster than a numerical gradient. We choose the Levenberg-Marquardt method to calibrate the model and do not observe multiple local minima reported in previous research. Two-dimensional sections show that the objective function is shaped as a narrow valley with a flat bottom. Our method is the fastest calibration of the Heston model developed so far and meets the speed requirement of practical trading.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stability of calibration procedures: fractals in the Black-Scholes model

Usually, in the Black-Scholes pricing theory the volatility is a positive real parameter. Here we explore what happens if it is allowed to be a complex number. The function for pricing a European option with a complex volatility has essential singularities at zero and infinity. The singularity at zero reflects the put-call parity. Solving for the implied volatility that reproduces a given market price yields not only a real root, but also infinitely many complex roots in a neighbourhood of the origin. The Newton-Raphson calculation of the complex implied volatility has a chaotic nature described by fractals.

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
preprint2010arXiv

Efficiency of linked cell algorithms

The linked cell list algorithm is an essential part of molecular simulation software, both molecular dynamics and Monte Carlo. Though it scales linearly with the number of particles, there has been a constant interest in increasing its efficiency, because a large part of CPU time is spent to identify the interacting particles. Several recent publications proposed improvements to the algorithm and investigated their efficiency by applying them to particular setups. In this publication we develop a general method to evaluate the efficiency of these algorithms, which is mostly independent of the parameters of the simulation, and test it for a number of linked cell list algorithms. We also propose a combination of linked cell reordering and interaction sorting that shows a good efficiency for a broad range of simulation setups.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

First-passage and first-exit times of a Bessel-like stochastic process

We study a stochastic process $X_t$ related to the Bessel and the Rayleigh processes, with various applications in physics, chemistry, biology, economics, finance and other fields. The stochastic differential equation is $dX_t = (nD/X_t) dt + \sqrt{2D} dW_t$, where $W_t$ is the Wiener process. Due to the singularity of the drift term for $X_t = 0$, different natures of boundary at the origin arise depending on the real parameter $n$: entrance, exit, and regular. For each of them we calculate analytically and numerically the probability density functions of first-passage times or first-exit times. Nontrivial behaviour is observed in the case of a regular boundary.

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

Fluctuating Interfaces in Liquid Crystals

We review and compare recent work on the properties of fluctuating interfaces between nematic and isotropic liquid-crystalline phases. Molecular dynamics and Monte Carlo simulations have been carried out for systems of ellipsoids and hard rods with aspect ratio 15:1, and the fluctuation spectrum of interface positions (the capillary wave spectrum) has been analyzed. In addition, the capillary wave spectrum has been calculated analytically within the Landau-de Gennes theory. The theory predicts that the interfacial fluctuations can be described in terms of a wave vector dependent interfacial tension, which is anisotropic at small wavelengths (stiff director regime) and becomes isotropic at large wavelengths (flexible director regime). After determining the elastic constants in the nematic phase, theory and simulation can be compared quantitatively. We obtain good agreement for the stiff director regime. The crossover to the flexible director regime is expected at wavelengths of the order of several thousand particle diameters, which was not accessible to our simulations.

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
preprint2006arXiv

Expressions for forces and torques in molecular simulations using rigid bodies

Expressions for intermolecular forces and torques, derived from pair potentials between rigid non-spherical units, are presented. The aim is to give compact and clear expressions, which are easily generalised, and which minimise the risk of error in writing molecular dynamics simulation programs. It is anticipated that these expressions will be useful in the simulation of liquid crystalline systems, and in coarse-grained modelling of macromolecules.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Influence of saving propensity on the power law tail of wealth distribution

Some general features of kinetic multi-agent models are reviewed, with particular attention to the relation between the agent saving propensities and the form of the equilibrium wealth distribution. The effect of a finite cutoff of the saving propensity distribution on the corresponding wealth distribution is studied. Various results about kinetic multi-agent models are collected and used to construct a realistic wealth distribution with zero limit for small values of wealth, an exponential form at intermediate and a power law tail at larger values of wealth.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Kinetic theory models for the distribution of wealth: power law from overlap of exponentials

Various multi-agent models of wealth distributions defined by microscopic laws regulating the trades, with or without a saving criterion, are reviewed. We discuss and clarify the equilibrium properties of the model with constant global saving propensity, resulting in Gamma distributions, and their equivalence to the Maxwell-Boltzmann kinetic energy distribution for a system of molecules in an effective number of dimensions $D_λ$, related to the saving propensity $λ$ [M. Patriarca, A. Chakraborti, and K. Kaski, Phys. Rev. E 70 (2004) 016104]. We use these results to analyze the model in which the individual saving propensities of the agents are quenched random variables, and the tail of the equilibrium wealth distribution exhibits a Pareto law $f(x) \propto x^{-α-1}$ with an exponent $α=1$ [A. Chatterjee, B. K. Chakrabarti, and S. S. Manna, Physica Scripta T106 (2003) 367]. Here, we show that the observed Pareto power law can be explained as arising from the overlap of the Maxwell-Boltzmann distributions associated to the various agents, which reach an equilibrium state characterized by their individual Gamma distributions. We also consider the influence of different types of saving propensity distributions on the equilibrium state.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Nematic-Isotropic Interfaces Under Shear: A Molecular Dynamics Simulation

We present a large-scale molecular dynamics study of nematic-paranematic interfaces under shear. We use a model of soft repulsive ellipsoidal particles with well-known equilibrium properties, and consider interfaces which are oriented normal to the direction of the shear gradient (common stress case). The director at the interface is oriented parallel to the interface (planar). A fixed average shear rate is imposed with Lees-Edwards boundary conditions, and the heat is dissipated with a profile-unbiased thermostat. First we study the properties of the interface at one particular shear rate in detail. The local interfacial profiles and the capillary wave fluctuations of the interfaces are calculated and compared with those of the corresponding equilibrium interface. Under shear, the interfacial width broadens and the capillary wave amplitudes at large wavelengths increase. The strain is distributed inhomogeneously in the system (shear banding), the local shear rate in the nematic region being distinctly higher than in the paranematic region. Surprisingly, we also observe (symmetry breaking) flow in the {\em vorticity} direction, with opposite direction in the nematic and the paranematic state. Finally, we investigate the stability of the interface for other shear rates and construct a nonequilibrium phase diagram.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Elastic constants from direct correlation functions in nematic liquid crystals: a computer simulation study

Density functional theories such as the Poniewierski-Stecki theory relate the elastic properties of nematic liquid crystals with their local liquid structure, i.e., with the direct correlation function (DCF) of the particles. We propose a way to determine the DCF in the nematic state from simulations without any approximations, taking into account the dependence of pair correlations on the orientation of the director explicitly. Using this scheme, we evaluate the Frank elastic constants K1, K2 and K3 in a system of soft ellipsoids. The values are in good agreement with those obtained directly from an analysis of order fluctuations. Our method thus establishes a reliable way to calculate elastic constants from pair distributions in computer simulations.

0 citations0 reviewsNo code linkedOpen access