BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

M. V. Tretyakov

M. V. Tretyakov contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NAmath.PRq-fin.PRphysics.comp-phq-fin.CPComputationcond-mat.mes-hallcond-mat.mtrl-sciq-fin.MF

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
9topics
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

12 published item(s)

preprint2021arXiv

Pricing FX Options under Intermediate Currency

We suggest an intermediate currency approach that allows us to price options on all FX markets simultaneously under the same risk-neutral measure which ensures consistency of FX option prices across all markets. In particular, it is sufficient to calibrate a model to the volatility smile on the domestic market as, due to the consistency of pricing formulas, the model automatically reproduces the correct smile for the inverse pair (the foreign market). We first consider the case of two currencies and then the multi-currency setting. We illustrate the intermediate currency approach by applying it to the Heston and SABR stochastic volatility models, to the model in which exchange rates are described by an extended skewed normal distribution, and also to the model-free approach of option pricing

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Block Circulant Embedding Method for Simulation of Stationary Gaussian Random Fields on Block-regular Grids

We propose a new method for sampling from stationary Gaussian random field on a grid which is not regular but has a regular block structure which is often the case in applications. The introduced block circulant embedding method (BCEM) can outperform the classical circulant embedding method (CEM) which requires a regularization of the irregular grid before its application. Comparison of BCEM vs CEM is performed on some typical model problems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Recursive Sparse Grid Collocation Method for Differential Equations with White Noise

We consider a sparse grid collocation method in conjunction with a time discretization of the differential equations for computing expectations of functionals of solutions to differential equations perturbed by time-dependent white noise. We first analyze the error of Smolyak's sparse grid collocation used to evaluate expectations of functionals of solutions to stochastic differential equations discretized by the Euler scheme. We show theoretically and numerically that this algorithm can have satisfactory accuracy for small magnitude of noise or small integration time, however it does not converge neither with decrease of the Euler scheme's time step size nor with increase of Smolyak's sparse grid level. Subsequently, we use this method as a building block for proposing a new algorithm by combining sparse grid collocation with a recursive procedure. This approach allows us to numerically integrate linear stochastic partial differential equations over longer times, which is illustrated in numerical tests on a stochastic advection-diffusion equation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

New Langevin and Gradient Thermostats for Rigid Body Dynamics

We introduce two new thermostats, one of Langevin type and one of gradient (Brownian) type, for rigid body dynamics. We formulate rotation using the quaternion representation of angular coordinates; both thermostats preserve the unit length of quaternions. The Langevin thermostat also ensures that the conjugate angular momenta stay within the tangent space of the quaternion coordinates, as required by the Hamiltonian dynamics of rigid bodies. We have constructed three geometric numerical integrators for the Langevin thermostat and one for the gradient thermostat. The numerical integrators reflect key properties of the thermostats themselves. Namely, they all preserve the unit length of quaternions, automatically, without the need of a projection onto the unit sphere. The Langevin integrators also ensure that the angular momenta remain within the tangent space of the quaternion coordinates. The Langevin integrators are quasi-symplectic and of weak order two. The numerical method for the gradient thermostat is of weak order one. Its construction exploits ideas of Lie-group type integrators for differential equations on manifolds. We numerically compare the discretization errors of the Langevin integrators, as well as the efficiency of the gradient integrator compared to the Langevin ones when used in the simulation of rigid TIP4P water model with smoothly truncated electrostatic interactions. We observe that the gradient integrator is computationally less efficient than the Langevin integrators. We also compare the relative accuracy of the Langevin integrators in evaluating various static quantities and give recommendations as to the choice of an appropriate integrator.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the long-time integration of stochastic gradient systems

This article addresses the weak convergence of numerical methods for Brownian dynamics. Typical analyses of numerical methods for stochastic differential equations focus on properties such as the weak order which estimates the asymptotic (stepsize $h \to 0$) convergence behavior of the error of finite time averages. Recently it has been demonstrated, by study of Fokker-Planck operators, that a non-Markovian numerical method [Leimkuhler and Matthews, 2013; Leimkuhler et al., 2013] generates approximations in the long time limit with higher accuracy order (2nd order) than would be expected from its weak convergence analysis (finite-time averages are 1st order accurate). In this article we describe the transition from the transient to the steady-state regime of this numerical method by estimating the time-dependency of the coefficients in an asymptotic expansion for the weak error, demonstrating that the convergence to 2nd order is exponentially rapid in time. Moreover, we provide numerical tests of the theory, including comparisons of the efficiencies of the Euler-Maruyama method, the popular 2nd order Heun method, and the non-Markovian method.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Layer methods for Navier-Stokes equations with additive noise

We propose and study a number of layer methods for stochastic Navier-Stokes equations (SNSE) with spatial periodic boundary conditions and additive noise. The methods are constructed using conditional probabilistic representations of solutions to SNSE and exploiting ideas of the weak sense numerical integration of stochastic differential equations. We prove some convergence results for the proposed methods. Results of numerical experiments on two model problems are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A fundamental mean-square convergence theorem for SDEs with locally Lipschitz coefficients and its applications

A version of the fundamental mean-square convergence theorem is proved for stochastic differential equations (SDE) which coefficients are allowed to grow polynomially at infinity and which satisfy a one-sided Lipschitz condition. The theorem is illustrated on a number of particular numerical methods, including a special balanced scheme and fully implicit methods. Some numerical tests are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Application of simplest random walk algorithms for pricing barrier options

We demonstrate effectiveness of the first-order algorithm from [Milstein, Tretyakov. Theory Prob. Appl. 47 (2002), 53-68] in application to barrier option pricing. The algorithm uses the weak Euler approximation far from barriers and a special construction motivated by linear interpolation of the price near barriers. It is easy to implement and is universal: it can be applied to various structures of the contracts including derivatives on multi-asset correlated underlyings and can deal with various type of barriers. In contrast to the Brownian bridge techniques currently commonly used for pricing barrier options, the algorithm tested here does not require knowledge of trigger probabilities nor their estimates. We illustrate this algorithm via pricing a barrier caplet, barrier trigger swap and barrier swaption.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Numerical integration of Heath-Jarrow-Morton model of interest rates

We propose and analyze numerical methods for the Heath-Jarrow-Morton (HJM) model. To construct the methods, we first discretize the infinite dimensional HJM equation in maturity time variable using quadrature rules for approximating the arbitrage-free drift. This results in a finite dimensional system of stochastic differential equations (SDEs) which we approximate in the weak and mean-square sense using the general theory of numerical integration of SDEs. The proposed numerical algorithms are computationally highly efficient due to the use of high-order quadrature rules which allow us to take relatively large discretization steps in the maturity time without affecting overall accuracy of the algorithms. Convergence theorems for the methods are proved. Results of some numerical experiments with European-type interest rate derivatives are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Convergence of Numerical Time-Averaging and Stationary Measures via Poisson Equations

Numerical approximation of the long time behavior of a stochastic differential equation (SDE) is considered. Error estimates for time-averaging estimators are obtained and then used to show that the stationary behavior of the numerical method converges to that of the SDE. The error analysis is based on using an associated Poisson equation for the underlying SDE. The main advantage of this approach is its simplicity and universality. It works equally well for a range of explicit and implicit schemes including those with simple simulation of random variables, and for hypoelliptic SDEs. To simplify the exposition, we consider only the case where the state space of the SDE is a torus and we study only smooth test functions. However we anticipate that the approach can be applied more widely. An analogy between our approach and Stein's method is indicated. Some practical implications of the results are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stable and fast semi-implicit integration of the stochastic Landau-Lifshitz equation

We propose new semi-implicit numerical methods for the integration of the stochastic Landau-Lifshitz equation with built-in angular momentum conservation. The performance of the proposed integrators is tested on the 1D Heisenberg chain. For this system, our schemes show better stability properties and allow us to use considerably larger time steps than standard explicit methods. At the same time, these semi-implicit schemes are also of comparable accuracy to and computationally much cheaper than the standard midpoint implicit method. The results are of key importance for atomistic spin dynamics simulations and the study of spin dynamics beyond the macro spin approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Langevin Thermostat for Rigid Body Dynamics

We present a new method for isothermal rigid body simulations using the quaternion representation and Langevin dynamics. It can be combined with the traditional Langevin or gradient (Brownian) dynamics for the translational degrees of freedom to correctly sample the NVT distribution in a simulation of rigid molecules. We propose simple, quasi-symplectic second-order numerical integrators and test their performance on the TIP4P model of water. We also investigate the optimal choice of thermostat parameters.

0 citations0 reviewsNo code linkedOpen access