pith. sign in

arxiv: 2606.31220 · v1 · pith:JSYZOAB3new · submitted 2026-06-30 · 🧮 math.PR

Efficient Computation Of Sensitivities For Derivatives In Energy Markets

Fred Espen Benth , Olfa Draouil , Farouk Hammami This is my paper

Pith reviewed 2026-07-01 04:41 UTC · model grok-4.3

classification 🧮 math.PR
keywords energy marketsdelta sensitivityMalliavin calculusdensity methodstochastic processesOrnstein-UhlenbeckCARMAjump diffusion
0
0 comments X

The pith

Delta sensitivities in energy markets can be computed using density or Malliavin methods on correlated price and volume processes

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a stochastic framework to calculate how the value of energy derivatives responds to changes in underlying prices when both prices and trading volumes follow linked random processes. It offers two approaches: a density method that requires knowing the probability distribution of the processes, and a Malliavin calculus method that works from the process dynamics alone without needing the density. This matters for risk management in energy markets, where volume uncertainty often correlates with price moves and hedging requires accurate sensitivities. The authors derive explicit formulas for models including Ornstein-Uhlenbeck and CARMA processes under the density approach, and for jump diffusions, time-changed Brownian motions, and NIG-driven processes under the Malliavin approach. Numerical checks on these models produce consistent delta values between the methods.

Core claim

In this study, we develop a stochastic framework for computing Delta sensitivities in energy markets, where both prices and traded volumes are modeled as correlated stochastic processes. Within this framework, we analyze two complementary approaches for sensitivity analysis: the density method, which is applicable when the density of the underlying process is known, and the Malliavin calculus method, which does not require any explicit knowledge of the density and relies only on the dynamics of the processes. We present illustrative examples for both methods. For the density-based approach, we consider Ornstein-Uhlenbeck and CARMA processes to model prices and energy volumes. For the Malliav

What carries the argument

Density method when the process density is known, and Malliavin calculus method using only process dynamics, applied to correlated price and volume stochastic processes

If this is right

  • Delta values follow directly from the derived formulas for Ornstein-Uhlenbeck and CARMA processes under the density method.
  • Malliavin formulas enable delta computation for jump-diffusion, time-changed Brownian motion, and NIG-driven processes without density knowledge.
  • Numerical implementations of both methods agree closely on overlapping models.
  • The framework supports efficient sensitivity calculation once the specific process dynamics are assumed.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The correlation structure between price and volume could be used to incorporate liquidity effects into hedging strategies.
  • Similar density or Malliavin derivations might apply to computing other sensitivities such as vega in the same market setting.
  • The methods could be applied to real traded volume and price time series from energy exchanges to check practical performance.

Load-bearing premise

The price and volume processes must follow the exact stochastic dynamics such as Ornstein-Uhlenbeck, CARMA, jump diffusion, or NIG for which the sensitivity formulas are derived.

What would settle it

A Monte Carlo simulation benchmark on one of the example models, such as an Ornstein-Uhlenbeck price-volume pair, that produces delta values differing substantially from the closed-form results would show the formulas do not hold.

read the original abstract

In this study, we develop a stochastic framework for computing Delta sensitivities in energy markets, where both prices and traded volumes are modeled as correlated stochastic processes. Within this framework, we analyze two complementary approaches for sensitivity analysis: the density method, which is applicable when the density of the underlying process is known, and the Malliavin calculus method, which does not require any explicit knowledge of the density and relies only on the dynamics of the processes. We present illustrative examples for both methods. For the density-based approach, we consider Ornstein-Uhlenbeck and CARMA processes to model prices and energy volumes. For the Malliavin calculus approach, we study Ornstein-Uhlenbeck processes, jump diffusion driven by a compound Poisson process, time-changed Brownian motion processes subordinated by an inverse Gaussian (IG) process, as well as Ornstein-Uhlenbeck processes driven by a normal inverse Gaussian (NIG) process. We provide some numerical examples illustrating the implementation of the proposed formulas and demonstrating a close agreement between the resulting delta estimates.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 3 minor

Summary. The manuscript develops a stochastic framework for computing Delta sensitivities of derivatives in energy markets, modeling both prices and traded volumes as correlated stochastic processes. It analyzes two complementary methods: the density method (applicable when the density is known), illustrated with Ornstein-Uhlenbeck and CARMA processes, and the Malliavin calculus method (not requiring explicit densities), applied to Ornstein-Uhlenbeck processes, compound-Poisson jump diffusions, inverse-Gaussian time-changed Brownian motion, and normal-inverse-Gaussian-driven Ornstein-Uhlenbeck processes. Numerical examples are included to illustrate implementation and demonstrate agreement between the two approaches.

Significance. If the derivations hold, the work supplies explicit sensitivity formulas for a range of models commonly used in energy markets, supporting practical risk-management computations under joint price-volume dynamics. The complementary density and Malliavin routes, together with the reported numerical agreement on simulated paths, constitute a concrete contribution to stochastic calculus applications in commodity derivatives.

minor comments (3)
  1. Abstract: the title emphasizes 'Efficient Computation' yet the abstract contains no discussion of computational cost, complexity, or runtime comparisons between the density and Malliavin implementations.
  2. The description of the joint correlation structure between price and volume processes is referenced but not stated explicitly in the provided abstract; a concise statement of the correlation assumption would improve readability.
  3. Numerical examples: the abstract reports 'close agreement' but supplies no information on simulation size, discretization scheme, or statistical error measures; adding these details would strengthen the validation section.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the careful summary of the manuscript and the positive assessment of its significance for stochastic calculus applications in commodity derivatives. The recommendation for minor revision is noted. However, the report lists no specific major comments, so we have no individual points to address at this time.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The manuscript derives Delta sensitivity formulas for energy market derivatives by applying the density method to OU and CARMA processes and Malliavin calculus to OU, compound-Poisson jump diffusion, IG-subordinated Brownian motion, and NIG-driven OU processes under a joint correlation structure. These derivations follow directly from standard stochastic calculus identities (e.g., integration by parts for Malliavin weights and explicit transition densities) applied to the given SDEs; no step reduces a claimed result to a fitted parameter, self-citation chain, or definitional tautology. Numerical illustrations compare the two independent methods on simulated paths and report agreement, which constitutes external validation rather than circular reinforcement. The central claims rest on process-specific but externally verifiable properties of the listed Lévy and Gaussian processes, rendering the derivation self-contained.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review provides no explicit free parameters, axioms, or invented entities; all modeling choices are implicit in the listed processes.

pith-pipeline@v0.9.1-grok · 5715 in / 1028 out tokens · 35750 ms · 2026-07-01T04:41:46.753188+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

27 extracted references · 27 canonical work pages

  1. [1]

    Approximations of small jumps of L´ evy processes with a view towards simulation.Journal of Applied Probability, 38(2):482–493, 2001

    Søren Asmussen and Jan Rosi´ nski. Approximations of small jumps of L´ evy processes with a view towards simulation.Journal of Applied Probability, 38(2):482–493, 2001

    work page 2001

  2. [2]

    Normal inverse Gaussian distributions and stochastic volatility modelling

    Ole E Barndorff-Nielsen. Normal inverse Gaussian distributions and stochastic volatility modelling. Scandinavian Journal of Statistics, 24(1):1–13, 1997

    work page 1997

  3. [3]

    Benth and P

    F.E. Benth and P. Kr¨ uhner.Stochastic Models for Price Dynamics in Energy and Commodity Markets: An Infinite-Dimensional Perspective. Springer Finance Series. Springer International Publishing AG, 2024

    work page 2024

  4. [4]

    Pricing and hedging quanto options in energy markets.Journal of Energy Markets, 8(1):—, 2015

    Fred Espen Benth, Nina Lange, and Tor Age Myklebust. Pricing and hedging quanto options in energy markets.Journal of Energy Markets, 8(1):—, 2015

    work page 2015

  5. [5]

    Robustness of option prices and their deltas in markets modelled by jump-diffusions.Communications on Stochastic Analysis, 5(2):285–307, 2011

    Fred Espen Benth, Giulia Di Nunno, and Asma Khedher. Robustness of option prices and their deltas in markets modelled by jump-diffusions.Communications on Stochastic Analysis, 5(2):285–307, 2011

    work page 2011

  6. [6]

    Computation of greeks in multifactor models with applications to power and commodity markets.Journal of Energy Markets, 5(4):3–31, 2012

    Fred Espen Benth, Giulia Di Nunno, and Asma Khedher. Computation of greeks in multifactor models with applications to power and commodity markets.Journal of Energy Markets, 5(4):3–31, 2012

    work page 2012

  7. [7]

    A non-Gaussian Ornstein–Uhlenbeck model for pricing wind power futures.Applied Mathematical Finance, 25(1):36–65, 2018

    Fred Espen Benth and Anca Pircalabu. A non-Gaussian Ornstein–Uhlenbeck model for pricing wind power futures.Applied Mathematical Finance, 25(1):36–65, 2018

    work page 2018

  8. [8]

    World Scientific, 2012

    Fred Espen Benth and Jurate Saltyte-Benth.Modelling and pricing in financial markets for weather derivatives, volume 17. World Scientific, 2012

    work page 2012

  9. [9]

    Brockwell

    Peter J. Brockwell. L´ evy-driven carma processes.Annals of the Institute of Statistical Mathematics, 53(1):113–124, 2001

    work page 2001

  10. [10]

    Ten things you should know about the dynamic conditional correlation representation.Econometrics, 1(1):115–126, 2013

    Massimiliano Caporin and Michael McAleer. Ten things you should know about the dynamic conditional correlation representation.Econometrics, 1(1):115–126, 2013

    work page 2013

  11. [11]

    Weather derivatives grow as risks intensify.OpenMarkets, 2024

    CME Group. Weather derivatives grow as risks intensify.OpenMarkets, 2024. Accessed: January, 2025

    work page 2024

  12. [12]

    On orthogonal polynomials and the Malliavin derivative for L´ evy stochastic measures

    Giulia Di Nunno. On orthogonal polynomials and the Malliavin derivative for L´ evy stochastic measures. InS´ eminaires et Congr` es, volume 16, pages 55–69. Soci´ et´ e Math´ ematique de France, 2007

    work page 2007

  13. [13]

    Stochastic integrals and adjoint derivatives

    Giulia Di Nunno and Yuri A Rozanov. Stochastic integrals and adjoint derivatives. InStochastic Analysis and Applications: The Abel Symposium 2005, pages 265–307. Springer, 2007

    work page 2005

  14. [14]

    Joseph L. Doob. The elementary Gaussian processes.Annals of Mathematical Statistics, 15(3):229–282, 1944

    work page 1944

  15. [15]

    Engle, Chowdhury Mustafa, and John Rice

    Robert F. Engle, Chowdhury Mustafa, and John Rice. Modelling peak electricity demand.Journal of Forecasting, 11(3):241–251, 1992

    work page 1992

  16. [16]

    Applications of Malliavin calculus to Monte-Carlo methods in finance

    Eric Fourni´ e, Jean-Michel Lasry, J´ erˆ ome Lebuchoux, and Pierre-Louis Lions. Applications of Malliavin calculus to Monte-Carlo methods in finance. II.Finance and Stochastics, 5:201–236, 2001

    work page 2001

  17. [17]

    Springer, 2003

    Paul Glasserman.Monte Carlo Methods in Financial Engineering, volume 53 ofApplications of Math- ematics. Springer, 2003

    work page 2003

  18. [18]

    El- sevier, 2014

    Nobuyuki Ikeda and Shinzo Watanabe.Stochastic Differential Equations and Diffusion Processes. El- sevier, 2014

    work page 2014

  19. [19]

    Spectral type of the shift transformation of differential processes with stationary increments

    Kiyosi Itˆ o. Spectral type of the shift transformation of differential processes with stationary increments. Transactions of the American Mathematical Society, 81(2):253–263, 1956

    work page 1956

  20. [20]

    Universitext

    Hui-Hsiung Kuo.Introduction to Stochastic Integration. Universitext. Springer, New York, 2006

    work page 2006

  21. [21]

    Lucia and Eduardo S

    Julio J. Lucia and Eduardo S. Schwartz. Electricity prices and power derivatives: Evidence from the nordic power exchange.The Journal of Derivatives, 9(3):5–50, 2002

    work page 2002

  22. [22]

    IMS Textbooks

    David Nualart and Eulalia Nualart.Introduction to Malliavin Calculus. IMS Textbooks. Cambridge University Press, 2002. 40 BENTH, DRAOUIL, AND HAMMAMI ENERGY DERIV ATIVES SENSITIVITIES

    work page 2002

  23. [23]

    Extreme weather ignites$25b industry.Governing Magazine, 2024

    Sam Potter. Extreme weather ignites$25b industry.Governing Magazine, 2024. Accessed: january, 2025

    work page 2024

  24. [24]

    The normal inverse Gaussian l´ evy process: simulation and approximation.Com- munications in statistics

    Tina Hviid Rydberg. The normal inverse Gaussian l´ evy process: simulation and approximation.Com- munications in statistics. Stochastic models, 13(4):887–910, 1997

    work page 1997

  25. [25]

    L´ evy Processes and Infinitely Divisible Distributions.Cambridge Studies in Advanced Mathematics, 68, 1999

    Ken-iti Sato. L´ evy Processes and Infinitely Divisible Distributions.Cambridge Studies in Advanced Mathematics, 68, 1999

    work page 1999

  26. [26]

    Canonical L´ evy process and Malliavin calculus

    Josep Llu´ ıs Sol´ e, Frederic Utzet, and Josep Vives. Canonical L´ evy process and Malliavin calculus. Stochastic Processes and Their Applications, 117(2):165–187, 2007

    work page 2007

  27. [27]

    Timmer and Peter J

    Reed P. Timmer and Peter J. Lamb. Relations between temperature and residential natural gas con- sumption in the central and eastern united states.Journal of Applied Meteorology and Climatology, 46(11):1993–2013, 2007. Department of Data Science and Analytics, BI Norwegian Business School, N-0484 Oslo, Nor- way & Department of Mathematics, University of O...

    work page 1993