Abstract:The present article aims to design and analyze efficient first-order strong schemes for a generalized Aït-Sahalia type model arising in mathematical finance and evolving in a positive domain $(0, \infty)$, which possesses a diffusion term with superlinear growth and a highly nonlinear drift that blows up at the origin. Such a complicated structure of the model unavoidably causes essential difficulties in the construction and convergence analysis of time discretizations. By incorporating implicitness in the term $\alpha_{-1} x^{-1}$ and a corrective mapping $\Phi_h$ in the recursion, we develop a novel class of explicit and unconditionally positivity-preserving (i.e., for any step-size $h>0$) Milstein-type schemes for the underlying model. In both non-critical and general critical cases, we introduce a novel approach to analyze mean-square error bounds of the novel schemes, without relying on a priori high-order moment bounds of the numerical approximations. The expected order-one mean-square convergence is attained for the proposed scheme. The above theoretical guarantee can be used to justify the optimal complexity of the Multilevel Monte Carlo method. Numerical experiments are finally provided to verify the theoretical findings.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to design and analyze an effective first - order strongly convergent numerical scheme applicable to Aït - Sahalia - type models. This type of model appears in mathematical finance and is characterized by a diffusion term with superlinear growth and a highly nonlinear drift term that diverges at the origin. This complex structure poses essential difficulties in the construction of time - discretization methods and the convergence analysis.
Specifically, the model is described as follows:
\[ dX_t=\left(\alpha_{ - 1}X_t^{- 1}-\alpha_0+\alpha_1X_t-\alpha_2X_t^r\right)dt+\sigma X_t^\rho dW_t,\quad t > 0,\quad X_0 = x_0>0. \]
where \(\alpha_{ - 1},\alpha_0,\alpha_1,\alpha_2>0\), \(r,\rho>1\) and satisfy \(r + 1\geq2\rho\). The drift term and the diffusion term in the model have polynomial growth and superlinear growth characteristics respectively, which makes it difficult for traditional numerical methods (such as the Euler - Maruyama method) to ensure the non - negativity and convergence of the solution.
To solve these problems, the author introduced a new explicit Milstein - type scheme, which has the following characteristics:
1. **Unconditional non - negativity**: For any step size \(h>0\), the numerical solution remains non - negative.
2. **First - order strong convergence**: In the non - critical case \(r + 1>2\rho\) and the general critical case \(r + 1 = 2\rho\), first - order strong convergence can be achieved.
3. **Easy to implement**: By introducing the correction mapping \(\Phi_h\) and partially implicit processing of \(\alpha_{ - 1}X_t^{- 1}\), the calculation process is simplified.
The main contribution of this method lies in proposing and analyzing for the first time an explicit, unconditionally non - negative numerical method with first - order strong convergence applicable to Aït - Sahalia - type models, thus filling the gap in existing research.