Abstract
A univariate stochastic system driven by multiplicative Gaussian white noise is considered. The standard method for simulating its Langevin equation of motion involves incrementing the system's state variable by a biased Gaussian random number at each time step. It is shown that the efficiency of such simulations can be significantly enhanced by incorporating the skewness of the distribution of the updated state variable. A new algorithm based on this principle is introduced, and its superior performance is demonstrated using a model of free diffusion of a Brownian particle with a friction coefficient that decreases exponentially with the kinetic energy. The proposed simulation technique proves to be accurate over time steps that are an order of magnitude longer than those required by standard algorithms. The model used to test the new numerical technique is known to exhibit a transition from normal diffusion to superdiffusion as the environmental temperature rises above a certain critical value. A simple empirical formula for the time-dependent diffusion coefficient, which covers both diffusion regimes, is introduced, and its accuracy is confirmed through comparison with the simulation results.