A Positivity-Preserving Numerical Method for the Free Cox-Ingersoll-Ross Equation
Abstract
We introduce a positivity-preserving drift-implicit square-root Euler method for the free Cox--Ingersoll--Ross (CIR) equation driven by free Brownian motion. The method is obtained by discretizing the square-root equation implicitly in the drift and preserves strict positivity of the numerical solution without any step-size restriction. Under the Feller condition, we establish the uniqueness of the global strictly positive solution and derive an upper bound for its operator norm. A key ingredient for the error analysis is a free stochastic differential representation of the inverse square-root process. Another key ingredient is an a priori grid-point error estimate for sufficiently small step sizes that yields uniform upper and strictly positive lower spectral bounds for the numerical solution on any fixed finite time interval. These results support the strong convergence of order $1$ and improve the classical results. Moreover, both the exact and numerical solutions are exponentially stable with respect to perturbations of the initial value. The numerical stability holds for arbitrary step sizes, and the associated discrete decay rate converges to its continuous counterpart as the step size tends to zero. Numerical experiments confirm the theoretical results.
Related
- Uniqueness and structure-preserving numerical simulations for an Ultrahyperbolic Schrödinger equation
- A second-order structure- and positivity-preserving convex limiting method for the Vlasov equations
- A compact high-order and positivity-preserving active flux method for compressible Navier--Stokes equations on Cartesian meshes
- Structure-Preserving Numerical Schemes for a Time-Fractional Porous Medium Equation with Nonlocal Pressure
- The nodally bound-preserving finite element method for the transport equation
BibTeX
@misc{niu2026positivitypreserving,
title = {{A Positivity-Preserving Numerical Method for the Free Cox-Ingersoll-Ross Equation}},
author = {Y. Niu and Z. Peng and J. Wei and Z. Yin and D. Zeng},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.39631},
eprint = {2610.39631},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.39631},
}
RIS
TY - UNPB
TI - A Positivity-Preserving Numerical Method for the Free Cox-Ingersoll-Ross Equation
AU - Y. Niu
AU - Z. Peng
AU - J. Wei
AU - Z. Yin
AU - D. Zeng
PY - 2026
DA - 2026/10/29
PB - arXiv
DO - 10.48550/arxiv.2610.39631
UR - https://arxiv.org/abs/2610.39631
AB - We introduce a positivity-preserving drift-implicit square-root Euler method for the free Cox--Ingersoll--Ross (CIR) equation driven by free Brownian motion. The method is obtained by discretizing the square-root equation implicitly in the drift and preserves strict positivity of the numerical solution without any step-size restriction. Under the Feller condition, we establish the uniqueness of the global strictly positive solution and derive an upper bound for its operator norm. A key ingredient for the error analysis is a free stochastic differential representation of the inverse square-root process. Another key ingredient is an a priori grid-point error estimate for sufficiently small step sizes that yields uniform upper and strictly positive lower spectral bounds for the numerical solution on any fixed finite time interval. These results support the strong convergence of order $1$ and improve the classical results. Moreover, both the exact and numerical solutions are exponentially stable with respect to perturbations of the initial value. The numerical stability holds for arbitrary step sizes, and the associated discrete decay rate converges to its continuous counterpart as the step size tends to zero. Numerical experiments confirm the theoretical results.
ER -