BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Jinqiao Duan

Jinqiao Duan contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.DSmath.PRMachine Learningnlin.CDmath.APmath.STStatistics Theorymath.NANumerical Analysisphysics.comp-phmath-phmath.MPphysics.ao-phcond-mat.stat-mechMethodologyMolecular Networksphysics.data-anphysics.flu-dynPopulations and Evolution

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
31works
0followers
19topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

31 published item(s)

preprint2024arXiv

Fourier neural operator based fluid-structure interaction for predicting the vesicle dynamics

Solving complex fluid-structure interaction (FSI) problems, characterized by nonlinear partial differential equations, is crucial in various scientific and engineering applications. Traditional computational fluid dynamics (CFD) solvers are insufficient to meet the growing requirements for large-scale and long-period simulations. Fortunately, the rapid advancement in neural networks, especially neural operator learning mappings between function spaces, has introduced novel approaches to tackle these challenges via data-driven modeling. In this paper, we propose a Fourier neural operator-based fluid-structure interaction solver (FNO-based FSI solver) for efficient simulation of FSI problems, where the solid solver based on the finite difference method is seamlessly integrated with the Fourier neural operator to predict incompressible flow using the immersed boundary method. We analyze the performance of the FNO-based FSI solver in the following three situations: training data with or without the steady state, training method with one-step label or multi-step labels, and prediction in interpolation or extrapolation. We find that the best performance for interpolation is achieved by training the operator with multi-step labels using steady-state data. Finally, we train the FNO-based FSI solver using this optimal training method and apply it to vesicle dynamics. The results show that the FNO-based FSI solver is capable of capturing the variations in the fluid and the vesicle.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Learning effective dynamics from data-driven stochastic systems

Multiscale stochastic dynamical systems have been widely adopted to a variety of scientific and engineering problems due to their capability of depicting complex phenomena in many real world applications. This work is devoted to investigating the effective dynamics for slow-fast stochastic dynamical systems. Given observation data on a short-term period satisfying some unknown slow-fast stochastic systems, we propose a novel algorithm including a neural network called Auto-SDE to learn invariant slow manifold. Our approach captures the evolutionary nature of a series of time-dependent autoencoder neural networks with the loss constructed from a discretized stochastic differential equation. Our algorithm is also validated to be accurate, stable and effective through numerical experiments under various evaluation metrics.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An end-to-end deep learning approach for extracting stochastic dynamical systems with $α$-stable Lévy noise

Recently, extracting data-driven governing laws of dynamical systems through deep learning frameworks has gained a lot of attention in various fields. Moreover, a growing amount of research work tends to transfer deterministic dynamical systems to stochastic dynamical systems, especially those driven by non-Gaussian multiplicative noise. However, lots of log-likelihood based algorithms that work well for Gaussian cases cannot be directly extended to non-Gaussian scenarios which could have high error and low convergence issues. In this work, we overcome some of these challenges and identify stochastic dynamical systems driven by $α$-stable Lévy noise from only random pairwise data. Our innovations include: (1) designing a deep learning approach to learn both drift and diffusion coefficients for Lévy induced noise with $α$ across all values, (2) learning complex multiplicative noise without restrictions on small noise intensity, (3) proposing an end-to-end complete framework for stochastic systems identification under a general input data assumption, that is, $α$-stable random variable. Finally, numerical experiments and comparisons with the non-local Kramers-Moyal formulas with moment generating function confirm the effectiveness of our method.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An Onsager-Machlup approach to the most probable transition pathway for a genetic regulatory network

We investigate a quantitative network of gene expression dynamics describing the competence development in Bacillus subtilis. First, we introduce an Onsager-Machlup approach to quantify the most probable transition pathway for both excitable and bistable dynamics. Then, we apply a machine learning method to calculate the most probable transition pathway via the Euler-Lagrangian equation. Finally, we analyze how the noise intensity affects the transition phenomena.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Data-driven method to learn the most probable transition pathway and stochastic differential equations

Transition phenomena between metastable states play an important role in complex systems due to noisy fluctuations. In this paper, the physics informed neural networks (PINNs) are presented to compute the most probable transition pathway. It is shown that the expected loss is bounded by the empirical loss. And the convergence result for the empirical loss is obtained. Then, a sampling method of rare events is presented to simulate the transition path by the Markovian bridge process. And we investigate the inverse problem to extract the stochastic differential equation from the most probable transition pathway data and the Markovian bridge process data, respectively. Finally, several numerical experiments are presented to verify the effectiveness of our methods.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Learning the temporal evolution of multivariate densities via normalizing flows

In this work, we propose a method to learn multivariate probability distributions using sample path data from stochastic differential equations. Specifically, we consider temporally evolving probability distributions (e.g., those produced by integrating local or nonlocal Fokker-Planck equations). We analyze this evolution through machine learning assisted construction of a time-dependent mapping that takes a reference distribution (say, a Gaussian) to each and every instance of our evolving distribution. If the reference distribution is the initial condition of a Fokker-Planck equation, what we learn is the time-T map of the corresponding solution. Specifically, the learned map is a multivariate normalizing flow that deforms the support of the reference density to the support of each and every density snapshot in time. We demonstrate that this approach can approximate probability density function evolutions in time from observed sampled data for systems driven by both Brownian and Lévy noise. We present examples with two- and three-dimensional, uni- and multimodal distributions to validate the method.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Detection of bistable structures via the Conley index and applications to biological systems

Bistability is a ubiquitous phenomenon in life sciences. In this paper, two kinds of bistable structures in dynamical systems are studied: One is two one-point attractors, another is a one-point attractor accompanied by a cycle attractor. By the Conley index theory, we prove that there exist other isolated invariant sets besides the two attractors, and also obtain the possible components and their configuration. Moreover, we find that there is always a separatrix or cycle separatrix, which separates the two attractors. Finally, the biological meanings and implications of these structures are given and discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Extracting Governing Laws from Sample Path Data of Non-Gaussian Stochastic Dynamical Systems

Advances in data science are leading to new progresses in the analysis and understanding of complex dynamics for systems with experimental and observational data. With numerous physical phenomena exhibiting bursting, flights, hopping, and intermittent features, stochastic differential equations with non-Gaussian Lévy noise are suitable to model these systems. Thus it is desirable and essential to infer such equations from available data to reasonably predict dynamical behaviors. In this work, we consider a data-driven method to extract stochastic dynamical systems with non-Gaussian asymmetric (rather than the symmetric) Lévy process, as well as Gaussian Brownian motion. We establish a theoretical framework and design a numerical algorithm to compute the asymmetric Lévy jump measure, drift and diffusion (i.e., nonlocal Kramers-Moyal formulas), hence obtaining the stochastic governing law, from noisy data. Numerical experiments on several prototypical examples confirm the efficacy and accuracy of this method. This method will become an effective tool in discovering the governing laws from available data sets and in understanding the mechanisms underlying complex random phenomena.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Pointwise convergence problem of Ostrovsky equation with rough data and random data

In this paper, we consider the pointwise convergence problem of free Ostrovsky equation with rough data and random data. Firstly, we show the almost everywhere pointwise convergence of free Ostrovsky equation in $H^{s}(\mathbb{R})$ with $s\geq \frac{1}{4}$ with rough data. Secondly, we present counterexamples showing that the maximal function estimate related to the free Ostrovsky equation can fail if $s<\frac{1}{4}$. Finally, for every $x\in \mathbb{R}$, we show the almost surely pointwise convergence of free Ostrovsky equation in $L^{2}(\mathbb{R})$ with random data. The main tools are the density theorem, high-low frequency idea, Wiener decomposition and Lemmas 2.1-2.6 as well as the probabilistic estimates of some random series which are just Lemmas 3.2-3.4 in this paper. The main difficulty is that zero is the singular point of the phase functions of free Ostrovsky equation. We use high-low frequency idea to conquer the difficulties.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Quantifying model uncertainty for the observed non-Gaussian data by the Hellinger distance

Mathematical models for complex systems under random fluctuations often certain uncertain parameters. However, quantifying model uncertainty for a stochastic differential equation with an $α$-stable Lévy process is still lacking. Here, we propose an approach to infer all the uncertain non-Gaussian parameters and other system parameters by minimizing the Hellinger distance over the parameter space. The Hellinger distance measures the similarity between an empirical probability density of non-Gaussian observations and a solution (as a probability density) of the associated nonlocal Fokker-Planck equation. Numerical experiments verify that our method is feasible for estimating single and multiple parameters. Meanwhile, we find an optimal estimation interval of the estimated parameters. This method is beneficial for extracting governing dynamical system models under non-Gaussian fluctuations, as in the study of abrupt climate changes in the Dansgaard-Oeschger events.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Data-Driven Approach for Discovering Stochastic Dynamical Systems with Non-Gaussian Levy Noise

With the rapid increase of valuable observational, experimental and simulating data for complex systems, great efforts are being devoted to discovering governing laws underlying the evolution of these systems. However, the existing techniques are limited to extract governing laws from data as either deterministic differential equations or stochastic differential equations with Gaussian noise. In the present work, we develop a new data-driven approach to extract stochastic dynamical systems with non-Gaussian symmetric Lévy noise, as well as Gaussian noise. First, we establish a feasible theoretical framework, by expressing the drift coefficient, diffusion coefficient and jump measure (i.e., anomalous diffusion) for the underlying stochastic dynamical system in terms of sample paths data. We then design a numerical algorithm to compute the drift, diffusion coefficient and jump measure, and thus extract a governing stochastic differential equation with Gaussian and non-Gaussian noise. Finally, we demonstrate the efficacy and accuracy of our approach by applying to several prototypical one-, two- and three-dimensional systems. This new approach will become a tool in discovering governing dynamical laws from noisy data sets, from observing or simulating complex phenomena, such as rare events triggered by random fluctuations with heavy as well as light tail statistical features.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Logistic-Harvest Model with Allee Effect under Multiplicative Noise

This work is devoted to the study of a stochastic logistic growth model with and without the Allee effect. Such a model describes the evolution of a population under environmental stochastic fluctuations and is in the form of a stochastic differential equation driven by multiplicative Gaussian noise. With the help of the associated Fokker-Planck equation, we analyze the population extinction probability and the probability of reaching a large population size before reaching a small one. We further study the impact of the harvest rate, noise intensity, and the Allee effect on population evolution. The analysis and numerical experiments show that if the noise intensity and harvest rate are small, the population grows exponentially, and upon reaching the carrying capacity, the population size fluctuates around it. In the stochastic logistic-harvest model without the Allee effect, when noise intensity becomes small (or goes to zero), the stationary probability density becomes more acute and its maximum point approaches one. However, for large noise intensity and harvest rate, the population size fluctuates wildly and does not grow exponentially to the carrying capacity. So as far as biological meanings are concerned, we must catch at small values of noise intensity and harvest rate. Finally, we discuss the biological implications of our results.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Machine Learning Framework for Computing the Most Probable Paths of Stochastic Dynamical Systems

The emergence of transition phenomena between metastable states induced by noise plays a fundamental role in a broad range of nonlinear systems. The computation of the most probable paths is a key issue to understand the mechanism of transition behaviors. Shooting method is a common technique for this purpose to solve the Euler-Lagrange equation for the associated action functional, while losing its efficacy in high-dimensional systems. In the present work, we develop a machine learning framework to compute the most probable paths in the sense of Onsager-Machlup action functional theory. Specifically, we reformulate the boundary value problem of Hamiltonian system and design a neural network to remedy the shortcomings of shooting method. The successful applications of our algorithms to several prototypical examples demonstrate its efficacy and accuracy for stochastic systems with both (Gaussian) Brownian noise and (non-Gaussian) Lévy noise. This novel approach is effective in exploring the internal mechanisms of rare events triggered by random fluctuations in various scientific fields.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An averaging principle for fractional stochastic differential equations with Lévy noise

This paper is devoted to the study of an averaging principle for fractional stochastic differential equations in Rnwith Lévy motion, using an integral transform method. We obtain a time-averaged equation under suitable assumptions. Furthermore, we show that the solutions of averaged equation approach the solutions of the original equation. Our results in this paper provide better understanding for effective approximation of fractional dynamical systems with non-Gaussian Lévy noise.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Analysis of multiscale methods for stochastic dynamical systems driven by $α$-stable processes

In this paper, we first analyze the strong and weak convergence of projective integration methods for multiscale stochastic dynamical systems driven by $α$-stable processes, which are used to estimate the effect that the fast components have on slow ones. Then we obtain the $p$th moment error bounds between the solution of slow component produced by projective integration method and the solution of effective system with $p \in \left(1, α\right)$. Finally, we corroborate our analytical results through a specific numerical example.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Discovering transition phenomena from data of stochastic dynamical systems with Lévy noise

It is a challenging issue to analyze complex dynamics from observed and simulated data. An advantage of extracting dynamic behaviors from data is that this approach enables the investigation of nonlinear phenomena whose mathematical models are unavailable. The purpose of this present work is to extract information about transition phenomena (e.g., mean exit time and escape probability), from data of stochastic differential equations with non-Gaussian Lévy noise. As a tool in describing dynamical systems, the Koopman semigroup transforms a nonlinear system into a linear system, but at the cost of elevating a finite dimensional problem into an infinite dimensional one. In spite of this, using the relation between the stochastic Koopman semigroup and the infinitesimal generator of a stochastic differential equation, we learn the mean exit time and escape probability from data. Specifically, we first obtain a finite dimensional approximation of the infinitesimal generator by an extended dynamic mode decomposition algorithm. Then we identify the drift coefficient, diffusion coefficient and anomalous diffusion coefficient for the stochastic differential equation. Finally, we compute the mean exit time and escape probability by finite difference discretization of the associated nonlocal partial differential equations. This approach is applicable in extracting transition information from data of stochastic differential equations with either (Gaussian) Brownian motion or (non-Gaussian) Lévy motion. We present one - and two-dimensional examples to demonstrate the effectiveness of our approach.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Effective reduction for a nonlocal Zakai stochastic partial differential equation in data assimilation

We study the effective reduction for a nonlocal stochastic partial differential equation with oscillating coefficients. The nonlocal operator in this stochastic partial differential equation is the generator of non-Gaussian Lévy processes, with either \textbf{integrable} or \textbf{non-integrable} jump kernels. We examine the limiting behavior of this equation as a scaling parameter tends to zero, and derive a reduced (local or nonlocal) effective equation. In particular, this work leads to an effective reduction for a data assimilation system with Lévy noise, by examining the corresponding nonlocal Zakai stochastic partial differential equation. We show that the probability density for the reduced data assimilation system approximates that for the original system.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Effective Wave Factorization for a Stochastic Schrödinger Equation

We study the homogenization of a stochastic Schrödinger equation with a large periodic potential in solid state physics. Denoting by $\varepsilon$ the period, the potential is scaled as $\varepsilon^{-2}$. Under a generic assumption on the spectral properties of the associated cell problem, we prove that the solution can be approximately factorized as the product of a fast oscillating cell eigenfunction and of a slowly varying solution of an effective equation. Our method is based on two-scale convergence and Bloch waves theory.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Identifying stochastic governing equations from data of the most probable transition trajectories

Extracting governing stochastic differential equation models from elusive data is crucial to understand and forecast dynamics for complex systems. We devise a method to extract the drift term and estimate the diffusion coefficient of a governing stochastic dynamical system, from its time-series data of the most probable transition trajectory. By the Onsager-Machlup theory, the most probable transition trajectory satisfies the corresponding Euler-Lagrange equation, which is a second order deterministic ordinary differential equation involving the drift term and diffusion coefficient. We first estimate the coefficients of the Euler-Lagrange equation based on the data of the most probable trajectory, and then we calculate the drift and diffusion coefficients of the governing stochastic dynamical system. These two steps involve sparse regression and optimization. Finally, we illustrate our method with an example and some discussions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Linear Response Theory for Nonlinear Stochastic Differential Equations with $α$-stable Lévy Noises

We consider a nonlinear stochastic differential equation driven by an $α$-stable Lévy process ($1<α<2$). We first obtain some regularity results for the probability density of its invariant measure via establishing the a priori estimate of the corresponding stationary Fokker-Planck equation. Then by the a priori estimate of Kolmogorov backward equations and the perturbation property of Markov semigroup, we derive the response function and generalize the famous linear response theory in nonequilibrium statistical mechanics to non-Gaussian stochastic dynamic systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Maximal Likely Phase Lines for a Reduced Ice Growth Model

We study the impact of Brownian noise on transitions between metastable equilibrium states in a stochastic ice sheet model. Two methods to accomplish different objectives are employed. The maximal likely trajectory by maximizing the probability density function and numerically solving the Fokker-Planck equation shows how the system will evolve over time. We have especially studied the maximal likely trajectories starting near the ice-free metastable state, and examined whether they evolve to or near the ice-covered metastable state for certain parameters, in order to gain insights into how the ice sheet formed. Furthermore, for the transition from ice-covered metastable state to the ice-free metastable state, we study the most probable path for various noise parameters via the Onsager-Machlup least action principle. This enables us to predict and visualize the melting process of the ice sheet if such a rare event ever does take place.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Maximum Likelihood Estimation of Stochastic Differential Equations with Random Effects Driven by Fractional Brownian Motion

Stochastic differential equations and stochastic dynamics are good models to describe stochastic phenomena in real world. In this paper, we study N independent stochastic processes Xi(t) with real entries and the processes are determined by the stochastic differential equations with drift term relying on some random effects. We obtain the Girsanov-type formula of the stochastic differential equation driven by Fractional Brownian Motion through kernel transformation. Under some assumptions of the random effect, we estimate the parameter estimators by the maximum likelihood estimation and give some numerical simulations for the discrete observations. Results show that for the different H, the parameter estimator is closer to the true value as the amount of data increases.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Most Probable Dynamics of Stochastic Dynamical Systems with Exponentially Light Jump Fluctuations

The emergence of the exit events from a bounded domain containing a stable fixed point induced by non-Gaussian Lévy fluctuations plays a pivotal role in practical physical systems. In the limit of weak noise, we develop a Hamiltonian formalism under the Lévy fluctuations with exponentially light jumps for one- and two-dimensional stochastic dynamical systems. This formalism is based on a recently proved large deviation principle for dynamical systems under non-Gaussian Lévy perturbations. We demonstrate how to compute the most probable exit path and the quasi-potential by several examples. Meanwhile, we explore the impacts of the jump measure on the quasi-potential quantitatively and on the most probable exit path qualitatively. Results show that the quasi-potential can be well estimated by an approximate analytical expression. Moreover, we discover that although the most probable exit paths are analogous to the Gaussian case for the isotropic noise, the anisotropic noise leads to significant changes of the structure of the exit paths. These findings shed light on the underlying qualitative mechanism and quantitative feature of the exit phenomenon induced by non-Gaussian noise.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical analysis and applications of Fokker-Planck equations for stochastic dynamical systems with multiplicative $α$-stable noises

The Fokker-Planck equations (FPEs) for stochastic systems driven by additive symmetric $α$-stable noises may not adequately describe the time evolution for the probability densities of solution paths in some practical applications, such as hydrodynamical systems, porous media, and composite materials. As a continuation of previous works on additive case, the FPEs for stochastic dynamical systems with multiplicative symmetric $α$-stable noises are derived by the adjoint operator method, which satisfy the nonlocal partial differential equations. A finite difference method for solving the nonlocal Fokker-Planck equation (FPE) is constructed, which is shown to satisfy the discrete maximum principle and to be convergent. Moreover, an example is given to illustrate this method. For asymmetric case, general finite difference schemes are proposed, and some analyses of the corresponding numerical schemes are given. Furthermore, the corresponding result is successfully applied to the nonlinear filtering problem.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Recurrent Solutions of a Nonautonomous Modified Swift-Hohenberg Equation

We consider recurrent solutions of the nonautonomous modified Swift-Hohenberg equation $$u_t+Δ^2u+2Δu+au+b|\nabla u|^2+u^3=g(t,x).$$ We employ Conley index theory to show that, if the forcing $g:\mathbb{R}\rightarrow L^2(Ω)$ is a recurrent function, then there are at least two recurrent solutions in $H_0^2(Ω)$ under appropriate assumptions on the parameters $a$, $b$ and $g$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Solving Inverse Stochastic Problems from Discrete Particle Observations Using the Fokker-Planck Equation and Physics-informed Neural Networks

The Fokker-Planck (FP) equation governing the evolution of the probability density function (PDF) is applicable to many disciplines but it requires specification of the coefficients for each case, which can be functions of space-time and not just constants, hence requiring the development of a data-driven modeling approach. When the data available is directly on the PDF, then there exist methods for inverse problems that can be employed to infer the coefficients and thus determine the FP equation and subsequently obtain its solution. Herein, we address a more realistic scenario, where only sparse data are given on the particles&#39; positions at a few time instants, which are not sufficient to accurately construct directly the PDF even at those times from existing methods, e.g., kernel estimation algorithms. To this end, we develop a general framework based on physics-informed neural networks (PINNs) that introduces a new loss function using the Kullback-Leibler divergence to connect the stochastic samples with the FP equation, to simultaneously learn the equation and infer the multi-dimensional PDF at all times. In particular, we consider two types of inverse problems, type I where the FP equation is known but the initial PDF is unknown, and type II in which, in addition to unknown initial PDF, the drift and diffusion terms are also unknown. In both cases, we investigate problems with either Brownian or Levy noise or a combination of both. We demonstrate the new PINN framework in detail in the one-dimensional case (1D) but we also provide results for up to 5D demonstrating that we can infer both the FP equation and} dynamics simultaneously at all times with high accuracy using only very few discrete observations of the particles.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Symplectic Euler scheme for Hamiltonian stochastic differential equations driven by Levy noise

This paper proposes a general symplectic Euler scheme for a class of Hamiltonian stochastic differential equations driven by L$\acute{e}$vy noise in the sense of Marcus form. The convergence of the symplectic Euler scheme for this Hamiltonian stochastic differential equations is investigated. Realizable numerical implementation of this scheme is also provided in details. Numerical experiments are presented to illustrate the effectiveness and superiority of the proposed method by the simulations of its orbits, symplectic structure and Hamlitonian.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The role of slow manifolds in parameter estimation for a multiscale stochastic system with $α$-stable Lévy noise

This work is about parameter estimation for a fast-slow stochastic system with non-Gaussian $α$-stable Lévy noise. When the observations are only available for slow components, a system parameter is estimated and the accuracy for this estimation is quantified by $p$-moment with $p\in(1, α)$, with the help of a reduced system through random slow manifold approximation. This method provides an advantage in computational complexity and cost, due to the dimension reduction in stochastic systems. To numerically illustrate this method, and to corroborate that the parameter estimator based on the reduced slow system is a good approximation for the true parameter value of the original system, a prototypical example is present.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The tipping times in an Arctic sea ice system under influence of extreme events

In light of the rapid recent retreat of Arctic sea ice, the extreme weather events triggering the variability in Arctic ice cover has drawn increasing attention. A non-Gaussian $α$-stable Lévy process is thought to be an appropriate model to describe such extreme event. The maximal likely trajectory, based on the nonlocal Fokker-Planck equation, is applied to a nonautonomous Arctic sea ice system under $α$-stable Lévy noise. Two types of tipping times, the early-warning tipping time and the disaster-happening tipping time, are used to predict the critical time for the maximal likely transition from a perennially ice-covered state to a seasonally ice-free one, and from a seasonally ice-free state to a perennially ice-free one, respectively. We find that the increased intensity of extreme events results in shorter warning time for sea ice melting, and that an enhanced greenhouse effect will intensify this influence, making the arrival of warning time significantly earlier. Meanwhile, for the enhanced greenhouse effect, we discover that increased intensity and frequency of extreme events will advance the disaster-happening tipping time, in which an ice-free state is maintained throughout the year in the Arctic Ocean. Finally, we identify values of Lévy index $α$ and noise intensity $ε$ in $αε$-space that can trigger a transition between the Arctic sea ice state. These results provide an effective theoretical framework for studying Arctic sea ice variations under the influence of extreme events.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Onsager-Machlup Function as Lagrangian for the Most Probable Path of a Jump-diffusion Process

This work is devoted to deriving the Onsager-Machlup function for a class of stochastic dynamical systems under (non-Gaussian) Levy noise as well as (Gaussian) Brownian noise, and examining the corresponding most probable paths. This Onsager-Machlup function is the Lagrangian giving the most probable path connecting metastable states for jump-diffusion processes. This is done by applying the Girsanov transformation for measures induced by jump-diffusion processes. Moreover, we have found this Lagrangian function is consistent with the result in the special case of diffusion processes. Finally, we apply this new Onsager-Machlup function to investigate dynamical behaviors analytically and numerically in several examples. These include the transitions from one metastable state to another metastable state in a double-well system, with numerical experiments illustrating most probable transition paths for various noise parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Enstrophy Dynamics of Stochastically Forced Large-Scale Geophysical Flows

Enstrophy is an averaged measure of fluid vorticity. This quantity is particularly important in {\em rotating} geophysical flows. We investigate the dynamical evolution of enstrophy for large-scale quasi-geostrophic flows under random wind forcing. We obtain upper bounds on the enstrophy, as well as results establishing its Hölder continuity and describing the small-time asymptotics.

0 citations0 reviewsNo code linkedOpen access