We consider boundary-value problems for differential equations of second order containing a Brownian motion (random perturbation) and a small parameter and prove a special existence and uniqueness theorem for random solutions. We study the asymptotic behaviour of these solutions as the small parameter goes to zero and show the stochastic averaging theorem for such equations. We find the explicit limits for the solutions as the small parameter goes to zero.

We study a stochastic differential equation driven by a Poisson point process, which models the continuous change in a population's environment, as well as the stochastic fixation of beneficial mutations that might compensate for this change. The fixation probability of a given mutation increases as the phenotypic lag Xt between the population and the optimum grows larger, and successful mutations are assumed to fix instantaneously (leading to an adaptive jump). Our main result is that the process is transient (i.e. converges to -∞, so that continued adaptation is impossible) if the rate of environmental change v exceeds a parameter m, which can be interpreted as the rate of adaptation in case every beneficial mutation becomes fixed with probability 1. If v < m, the process is Harris recurrent and possesses a unique invariant probability measure, while in the limiting case m = v, Harris recurrence with an infinite invariant measure or transience depends upon additional technical conditions. We show how our results can be extended to a class of time varying rates of environmental change.

In this paper we consider a diffusive stochastic predator–prey model with a nonlinear functional response and the randomness is assumed to be of Gaussian nature. A large deviation principle is established for solution processes of the considered model by implementing the weak convergence technique.

We provide an elementary method for exploring pricing problems of one spread options within a fractional Wick–Itô–Skorohod integral framework. Its underlying assets come from two different interactive markets that are modelled by two mixed fractional Black–Scholes models with Hurst parameters, $H_{1}\neq H_{2}$, where $1/2\leq H_{i}<1$ for $i=1,2$. Pricing formulae of these options with respect to strike price $K=0$ or $K\neq 0$ are given, and their application to the real market is examined.

The deferred correction (DC) method is a classical method for solving ordinary differential equations; one of its key features is to iteratively use lower order numerical methods so that high-order numerical scheme can be obtained. The main advantage of the DC approach is its simplicity and robustness. In this paper, the DC idea will be adopted to solve forward backward stochastic differential equations (FBSDEs) which have practical importance in many applications. Noted that it is difficult to design high-order and relatively “clean” numerical schemes for FBSDEs due to the involvement of randomness and the coupling of the FSDEs and BSDEs. This paper will describe how to use the simplest Euler method in each DC step–leading to simple computational complexity–to achieve high order rate of convergence.

We investigate rare or small probability events in the context of large deviations of the stochastic Camassa–Holm equation. By the weak convergence approach and regularization, we get large deviations of the regularized equation. Then, by stochastic equations exponentially equivalent to the corresponding laws, we get large deviations of the stochastic Camassa–Holm equation.

This is one of our series papers on multistep schemes for solving forward backward stochastic differential equations (FBSDEs) and related problems. Here we extend (with non-trivial updates) our multistep schemes in [W. Zhao, Y. Fu and T. Zhou, SIAM J. Sci. Comput., 36 (2014), pp. A1731-A1751] to solve the second-order FBSDEs (2FBSDEs). The key feature of the multistep schemes is that the Euler method is used to discretize the forward SDE, which dramatically reduces the entire computational complexity. Moreover, it is shown that the usual quantities of interest (e.g., the solution tuple (Yt,Zt,At,Γt) of the 2FBSDEs) are still of high order accuracy. Several numerical examples are given to show the effectiveness of the proposed numerical schemes. Applications of our numerical schemes to stochastic optimal control problems are also presented.

In this paper, we propose a compact scheme to numerically study the coupled stochastic nonlinear Schrödinger equations. We prove that the compact scheme preserves the discrete stochastic multi-symplectic conservation law, discrete charge conservation law and discrete energy evolution law almost surely. Numerical experiments confirm well the theoretical analysis results. Furthermore, we present a detailed numerical investigation of the optical phenomena based on the compact scheme. By numerical experiments for various amplitudes of noise, we find that the noise accelerates the oscillation of the soliton and leads to the decay of the solution amplitudes with respect to time. In particular, if the noise is relatively strong, the soliton will be totally destroyed. Meanwhile, we observe that the phase shift is sensibly modified by the noise. Moreover, the numerical results present inelastic interaction which is different from the deterministic case.

We study the construction of symplectic Runge-Kutta methods for stochastic Hamiltonian systems (SHS). Three types of systems, SHS with multiplicative noise, special separable Hamiltonians and multiple additive noise, respectively, are considered in this paper. Stochastic Runge-Kutta (SRK) methods for these systems are investigated, and the corresponding conditions for SRK methods to preserve the symplectic property are given. Based on the weak/strong order and symplectic conditions, some effective schemes are derived. In particular, using the algebraic computation, we obtained two classes of high weak order symplectic Runge-Kutta methods for SHS with a single multiplicative noise, and two classes of high strong order symplectic Runge-Kutta methods for SHS with multiple multiplicative and additive noise, respectively. The numerical case studies confirm that the symplectic methods are efficient computational tools for long-term simulations.

We study the distributional properties of jumps in a continuous-state branching process with immigration. In particular, a representation is given for the distribution of the first jump time of the process with jump size in a given Borel set. From this result we derive a characterization for the distribution of the local maximal jump of the process. The equivalence of this distribution and the total Lévy measure is then studied. For the continuous-state branching process without immigration, we also study similar problems for its global maximal jump.

In this paper we uniformly approximate the trajectories of the Cox–Ingersoll–Ross (CIR) process. At a sequence of random times the approximate trajectories will be even exact. In between, the approximation will be uniformly close to the exact trajectory. From a conceptual point of view, the proposed method gives a better quality of approximation in a path-wise sense than standard, or even exact, simulation of the CIR dynamics at some deterministic time grid.

In this paper, we investigate the mean-square convergence of the split-step θ-scheme for nonlinear stochastic differential equations with jumps. Under some standard assumptions, we rigorously prove that the strong rate of convergence of the split-step θ-scheme in strong sense is one half. Some numerical experiments are carried out to assert our theoretical result.

By introducing a new Gaussian process and a new compensated Poisson random measure, we propose an explicit prediction-correction scheme for solving decoupled forward backward stochastic differential equations with jumps (FBSDEJs). For this scheme, we first theoretically obtain a general error estimate result, which implies that the scheme is stable. Then using this result, we rigorously prove that the accuracy of the explicit scheme can be of second order. Finally, we carry out some numerical experiments to verify our theoretical results.

This paper presents a computational method for solving stochastic Ito-Volterra integral equations. First, Haar wavelets and their properties are employed to derive a general procedure for forming the stochastic operational matrix of Haar wavelets. Then, application of this stochastic operational matrix for solving stochastic Ito-Volterra integral equations is explained. The convergence and error analysis of the proposed method are investigated. Finally, the efficiency of the presented method is confirmed by some examples.

Based on the martingale theory and large deviation techniques, we investigate the pth moment exponential stability criterion of the exact and numerical solutions to hybrid stochastic differential equations (SDEs) under the local Lipschitz condition. This new stability criterion shows that Markovian switching can serve as a stochastic stabilizing factor by its logarithmic moment-generating function. We also investigate the pth moment exponential stability of Euler–Maruyama (EM), backward EM (BEM) and split-step backward EM (SSBEM) approximations for hybrid SDEs and show that, under the additional linear growth condition, the EM method can share the mean-square exponential stability of the exact solution for sufficiently small step size. However, the BEM method can work without the linear growth condition. We further investigate the SSBEM method under a coupled condition.

We introduce exact methods for the simulation of sample paths of one-dimensional diffusions with a discontinuity in the drift function. Our procedures require the simulation of finite-dimensional candidate draws from probability laws related to those of Brownian motion and its local time, and are based on the principle of retrospective rejection sampling. A simple illustration is provided.

In this paper we consider linear functions constructed on two different weighted branching processes and provide explicit bounds for their Kantorovich–Rubinstein distance in terms of couplings of their corresponding generic branching vectors. Motivated by applications to the analysis of random graphs, we also consider a variation of the weighted branching process where the generic branching vector has a different dependence structure from the usual one. By applying the bounds to sequences of weighted branching processes, we derive sufficient conditions for the convergence in the Kantorovich–Rubinstein distance of linear functions. We focus on the case where the limits are endogenous fixed points of suitable smoothing transformations.

Various numerical methods have been developed in order to solve complex systems with uncertainties, and the stochastic collocation method using l1-minimisation on low discrepancy point sets is investigated here. Halton and Sobol' sequences are considered, and low discrepancy point sets and random points are compared. The tests discussed involve a given target function in polynomial form, high-dimensional functions and a random ODE model. Our numerical results show that the low discrepancy point sets perform as well or better than random sampling for stochastic collocation via l1-minimisation.

Upon a set of backward orthogonal polynomials, we propose a novel multi-step numerical scheme for solving the decoupled forward-backward stochastic differential equations (FBSDEs). Under Lipschtiz conditions on the coefficients of the FBSDEs, we first get a general error estimate result which implies zero-stability of the proposed scheme, and then we further prove that the convergence rate of the scheme can be of high order for Markovian FBSDEs. Some numerical experiments are presented to demonstrate the accuracy of the proposed multi-step scheme and to numerically verify the theoretical results.