Abstract We establish theoretical properties of the solution to a two-variance-driven interest rate model with super-linear coefficient terms. Since this model is not tractable analytically, we construct an implementable numerical method to approximate it and prove the finite-time strong convergence theory under the local Lipschitz condition. Finally, we provide simulation examples to demonstrate the theoretical results.
{"title":"Strong approximation of a two-factor stochastic volatility model under local Lipschitz condition","authors":"Emmanuel Coffie","doi":"10.1515/mcma-2023-2021","DOIUrl":"https://doi.org/10.1515/mcma-2023-2021","url":null,"abstract":"Abstract We establish theoretical properties of the solution to a two-variance-driven interest rate model with super-linear coefficient terms. Since this model is not tractable analytically, we construct an implementable numerical method to approximate it and prove the finite-time strong convergence theory under the local Lipschitz condition. Finally, we provide simulation examples to demonstrate the theoretical results.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-11-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135041948","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Evgeniya Kablukova, Karl K. Sabelfeld, Dmitry Protasov, Konstantin Zhuravlev
Abstract In this paper we develop a stochastic simulation algorithm for electron transport in a DA-pHEMT heterostructure. Mathematical formulation of the problem of electron gas transport in the heterostructure in the form of a coupled system of Poisson, Schrödinger and kinetic Boltzmann equations is given. A Monte Carlo model of electron transport in DA-pHEMT heterostructures which accounts for multivalley parabolic band structure, as well as relevant formulas for calculating electron scattering rates and scattering phase functions on polar optical, intervalley phonons and on impurities are developed. The results of a computational experiment involving the solution of the system of Poisson–Schrödinger–Boltzmann equations for the AlGaAs / GaAs / InGaAs / GaAs / AlGaAs heterostructure are presented. The distribution of electrons by energy subband in the main and satellite valleys and the field dependences of the electron drift velocity in each valley are calculated. It was discovered that there is no spatial transfer of electrons into wide-gap AlGaAs layers due to high barriers created by modulated-doped impurities. A comparative analysis of the electron drift velocities in the studied DA-pHEMT heterostructures and in the unstrained layer of the InGaAs is given.
{"title":"Stochastic simulation of electron transport in a strong electrical field in low-dimensional heterostructures","authors":"Evgeniya Kablukova, Karl K. Sabelfeld, Dmitry Protasov, Konstantin Zhuravlev","doi":"10.1515/mcma-2023-2019","DOIUrl":"https://doi.org/10.1515/mcma-2023-2019","url":null,"abstract":"Abstract In this paper we develop a stochastic simulation algorithm for electron transport in a DA-pHEMT heterostructure. Mathematical formulation of the problem of electron gas transport in the heterostructure in the form of a coupled system of Poisson, Schrödinger and kinetic Boltzmann equations is given. A Monte Carlo model of electron transport in DA-pHEMT heterostructures which accounts for multivalley parabolic band structure, as well as relevant formulas for calculating electron scattering rates and scattering phase functions on polar optical, intervalley phonons and on impurities are developed. The results of a computational experiment involving the solution of the system of Poisson–Schrödinger–Boltzmann equations for the AlGaAs / GaAs / InGaAs / GaAs / AlGaAs heterostructure are presented. The distribution of electrons by energy subband in the main and satellite valleys and the field dependences of the electron drift velocity in each valley are calculated. It was discovered that there is no spatial transfer of electrons into wide-gap AlGaAs layers due to high barriers created by modulated-doped impurities. A comparative analysis of the electron drift velocities in the studied DA-pHEMT heterostructures and in the unstrained layer of the InGaAs is given.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135216484","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract In this paper, under some suitable conditions, we prove existence of a strong solution and uniqueness for the perturbed stochastic differential equations with reflected boundary (PSDERB), that is, {x(t)=x(0)+∫0tσ(s,x(s))dB(s)+∫0tb(s,x(s))ds+α(t)H(max0≤u≤tx(u))+β(t)Lt0(x)
文摘本文在一些合适的条件下,我们证明强解的存在和唯一性的摄动随机微分方程反映边界(PSDERB) , { x(t ) = x(0)+∫0 tσ(x(年代 ) ) d B(s ) + ∫0 t b(x(年代 ) ) d s +α(t)Hu (max 0≤≤tx(u ) ) + β(t)L t 0(x)x(t ) 所有≥0 t≥0 , 左{{对齐}{}开始x (t)和= x (0) + int_ {0} ^ {t} σ(年代,x (s)) , dB (s) + int_ {0} ^ {t} b (s, x (s)) d + α(t) H bigl {(} max_ {0 leq u leq t} x (u) bigr{)} + β(t) L_ {t} ^ {0} (x) x (t)和 组0 四文本所有}{ t 组0 {对齐}正确的结束。其中𝐻为连续的r值函数,σ,b, α sigma,b, α alpha,和时延为可测函数,L t 0 L_{t}^{0}表示半鞅变量的时间在零点处的局部时间。
Abstract It is well known that shortened modeling of particle trajectories with the use multiplicative statistical weights, as a rule, increases the efficiency of the program (in terms of accuracy/time ratio). This trick is often used in non-branching schemes simulating transfer processes without multiplication (for example, the transfer of X-ray radiation), in which it is sufficient to confine ourselves to studying only the average values of the field characteristics. With an increase in energy, however, multiplication processes begin to play a significant role (the production of electron-photon pairs by gamma quanta with energies above 1.022 MeV, etc.), when the resulting trajectory is not just a broken curve in the phase space, but a branched tree. This technique is also applicable to this process, but only if the study of statistical fluctuations and correlations is not the purpose of the calculation. The present review contains the basic concepts of the Monte Carlo method as applied to the theory of particle transport, demonstration of the weighting method in non-branching processes, and ends with a discussion of unbiased estimates of the second moment and covariance of additive functionals.
{"title":"A weight Monte Carlo estimation of fluctuations in branching processes","authors":"Vladimir Uchaikin, Elena Kozhemiakina","doi":"10.1515/mcma-2023-2015","DOIUrl":"https://doi.org/10.1515/mcma-2023-2015","url":null,"abstract":"Abstract It is well known that shortened modeling of particle trajectories with the use multiplicative statistical weights, as a rule, increases the efficiency of the program (in terms of accuracy/time ratio). This trick is often used in non-branching schemes simulating transfer processes without multiplication (for example, the transfer of X-ray radiation), in which it is sufficient to confine ourselves to studying only the average values of the field characteristics. With an increase in energy, however, multiplication processes begin to play a significant role (the production of electron-photon pairs by gamma quanta with energies above 1.022 MeV, etc.), when the resulting trajectory is not just a broken curve in the phase space, but a branched tree. This technique is also applicable to this process, but only if the study of statistical fluctuations and correlations is not the purpose of the calculation. The present review contains the basic concepts of the Monte Carlo method as applied to the theory of particle transport, demonstration of the weighting method in non-branching processes, and ends with a discussion of unbiased estimates of the second moment and covariance of additive functionals.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-10-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135220233","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract There is no method of predicting the price of an option other than hedging strategies such as the binomial hedging strategy, the Black–Scholes hedging strategy and others. We will study these two basic hedging strategies in terms of their feasibility, and we will see that the Black–Scholes hedging strategy is not feasible because this strategy demands instantaneously rebuilding the replicating portfolio. Consequently, the real world prices of the options are not relevant at all with the Black–Scholes hedging strategy! We will suitably redefine the binomial hedging strategy so that it will be practically useful and present other feasible and generally more effective hedging strategies with some of them practically useful for options with no tradable underlying assets. Finally, we will mention some open questions related to the above.
{"title":"Option pricing: Examples and open problems","authors":"Nikolaos Halidias","doi":"10.1515/mcma-2023-2014","DOIUrl":"https://doi.org/10.1515/mcma-2023-2014","url":null,"abstract":"Abstract There is no method of predicting the price of an option other than hedging strategies such as the binomial hedging strategy, the Black–Scholes hedging strategy and others. We will study these two basic hedging strategies in terms of their feasibility, and we will see that the Black–Scholes hedging strategy is not feasible because this strategy demands instantaneously rebuilding the replicating portfolio. Consequently, the real world prices of the options are not relevant at all with the Black–Scholes hedging strategy! We will suitably redefine the binomial hedging strategy so that it will be practically useful and present other feasible and generally more effective hedging strategies with some of them practically useful for options with no tradable underlying assets. Finally, we will mention some open questions related to the above.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-10-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135220234","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Guilhem Balvet, Jean-Pierre Minier, Yelva Roustan, Martin Ferrand
Abstract Lagrangian stochastic methods are widely used to model turbulent flows. Scarce consideration has, however, been devoted to the treatment of the near-wall region and to the formulation of a proper wall-boundary condition. With respect to this issue, the main purpose of this paper is to present an in-depth analysis of such flows when relying on particle/mesh formulations of the probability density function (PDF) model. This is translated into three objectives. The first objective is to assess the existing an-elastic wall-boundary condition and present new validation results. The second objective is to analyse the impact of the interpolation of the mean fields at particle positions on their dynamics. The third objective is to investigate the spatial error affecting covariance estimators when they are extracted on coarse volumes. All these developments allow to ascertain that the key dynamical statistics of wall-bounded flows are properly captured even for coarse spatial resolutions.
{"title":"Analysis of wall-modelled particle/mesh PDF methods for turbulent parietal flows","authors":"Guilhem Balvet, Jean-Pierre Minier, Yelva Roustan, Martin Ferrand","doi":"10.1515/mcma-2023-2017","DOIUrl":"https://doi.org/10.1515/mcma-2023-2017","url":null,"abstract":"Abstract Lagrangian stochastic methods are widely used to model turbulent flows. Scarce consideration has, however, been devoted to the treatment of the near-wall region and to the formulation of a proper wall-boundary condition. With respect to this issue, the main purpose of this paper is to present an in-depth analysis of such flows when relying on particle/mesh formulations of the probability density function (PDF) model. This is translated into three objectives. The first objective is to assess the existing an-elastic wall-boundary condition and present new validation results. The second objective is to analyse the impact of the interpolation of the mean fields at particle positions on their dynamics. The third objective is to investigate the spatial error affecting covariance estimators when they are extracted on coarse volumes. All these developments allow to ascertain that the key dynamical statistics of wall-bounded flows are properly captured even for coarse spatial resolutions.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-10-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135666713","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract In this paper, we consider aggregated functional data composed by a linear combination of component curves and the problem of estimating these component curves. We propose the application of a bayesian wavelet shrinkage rule based on a mixture of a point mass function at zero and the logistic distribution as prior to wavelet coefficients to estimate mean curves of components. This procedure has the advantage of estimating component functions with important local characteristics such as discontinuities, spikes and oscillations for example, due the features of wavelet basis expansion of functions. Simulation studies were done to evaluate the performance of the proposed method, and its results are compared with a spline-based method. An application on the so-called Tecator dataset is also provided.
{"title":"A wavelet-based method in aggregated functional data analysis","authors":"Alex Rodrigo dos Santos Sousa","doi":"10.1515/mcma-2023-2016","DOIUrl":"https://doi.org/10.1515/mcma-2023-2016","url":null,"abstract":"Abstract In this paper, we consider aggregated functional data composed by a linear combination of component curves and the problem of estimating these component curves. We propose the application of a bayesian wavelet shrinkage rule based on a mixture of a point mass function at zero and the logistic distribution as prior to wavelet coefficients to estimate mean curves of components. This procedure has the advantage of estimating component functions with important local characteristics such as discontinuities, spikes and oscillations for example, due the features of wavelet basis expansion of functions. Simulation studies were done to evaluate the performance of the proposed method, and its results are compared with a spline-based method. An application on the so-called Tecator dataset is also provided.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-10-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135765975","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract In this paper we suggest randomized linear solvers with a focus on refinement issue to achieve a high precision while maintaining all the advantages of the Monte Carlo method for solving systems of large dimension with dense matrices. It is shown that each iterative refinement step reduces the error by one order of magnitude. The crucial point of the suggested method is, in contrast to the standard Monte Carlo method, that the randomized vector algorithm computes the entire solution column at once, rather than a single component. This makes it possible to efficiently construct the iterative refinement method. We apply the developed method for solving a system of elasticity equations.
{"title":"Randomized vector iterative linear solvers of high precision for large dense system","authors":"Karl K. Sabelfeld, Anastasiya Kireeva","doi":"10.1515/mcma-2023-2013","DOIUrl":"https://doi.org/10.1515/mcma-2023-2013","url":null,"abstract":"Abstract In this paper we suggest randomized linear solvers with a focus on refinement issue to achieve a high precision while maintaining all the advantages of the Monte Carlo method for solving systems of large dimension with dense matrices. It is shown that each iterative refinement step reduces the error by one order of magnitude. The crucial point of the suggested method is, in contrast to the standard Monte Carlo method, that the randomized vector algorithm computes the entire solution column at once, rather than a single component. This makes it possible to efficiently construct the iterative refinement method. We apply the developed method for solving a system of elasticity equations.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-10-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135548383","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2023-09-01DOI: 10.1515/mcma-2023-frontmatter3
{"title":"Frontmatter","authors":"","doi":"10.1515/mcma-2023-frontmatter3","DOIUrl":"https://doi.org/10.1515/mcma-2023-frontmatter3","url":null,"abstract":"","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"136310148","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract In this paper, we will consider periodic EGARCH ( p , p ) {operatorname{EGARCH}(p,p)} (exponential generalized autoregressive conditional heteroscedastic) processes denoted by PEGARCH ( p , p ) {operatorname{PEGARCH}(p,p)} . These processes are similar to the standard EGARCH processes, but include seasonally varying coefficients. We examine the probabilistic structure of an EGARCH-type stochastic difference equation with periodically-varying parameters. We propose necessary and sufficient conditions ensuring the existence of stationary solutions (in a periodic sense) based on a Markovian representation. The closed forms of higher moments are, under these conditions, established. Furthermore, the expressions for the Kurtosis coefficient and the autocorrelations of squared observations are derived. The general theory is illustrated by considering special cases such as the symmetric and the asymmetric cases of the second order PEGARCH model.
{"title":"On the stationarity and existence of moments of the periodic EGARCH process","authors":"Ines Lescheb, Walid Slimani","doi":"10.1515/mcma-2023-2011","DOIUrl":"https://doi.org/10.1515/mcma-2023-2011","url":null,"abstract":"Abstract In this paper, we will consider periodic EGARCH ( p , p ) {operatorname{EGARCH}(p,p)} (exponential generalized autoregressive conditional heteroscedastic) processes denoted by PEGARCH ( p , p ) {operatorname{PEGARCH}(p,p)} . These processes are similar to the standard EGARCH processes, but include seasonally varying coefficients. We examine the probabilistic structure of an EGARCH-type stochastic difference equation with periodically-varying parameters. We propose necessary and sufficient conditions ensuring the existence of stationary solutions (in a periodic sense) based on a Markovian representation. The closed forms of higher moments are, under these conditions, established. Furthermore, the expressions for the Kurtosis coefficient and the autocorrelations of squared observations are derived. The general theory is illustrated by considering special cases such as the symmetric and the asymmetric cases of the second order PEGARCH model.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":null,"pages":null},"PeriodicalIF":0.9,"publicationDate":"2023-08-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44711818","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}