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

Jian-Guo Liu

Jian-Guo Liu contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.APmath.NANumerical Analysismath.PRmath-phmath.DSmath.MPMachine Learningmath.CAPopulations and Evolution

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
10topics
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

16 published item(s)

preprint2024arXiv

Some Grönwall inequalities for a class of discretizations of time fractional equations on nonuniform meshes

We consider the completely positive discretizations of fractional ordinary differential equations (FODEs) on nonuniform meshes. Making use of the resolvents for nonuniform meshes, we first establish comparison principles for the discretizations. Then we prove some discrete Grönwall inequalities using the comparison principles and careful analysis of the solutions to the time continuous FODEs. Our results do not have any restrictions on the step size ratio. The Grönwall inequalities for dissipative equations can be used to obtain the uniform-in-time error control and decay estimates of the numerical solutions. The Grönwall inequalities are then applied to subdiffusion problems and the time fractional Allen-Cahn equations for illustration.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Diffusion approximations of Oja's online principal component analysis

Oja's algorithm of principal component analysis (PCA) has been one of the methods utilized in practice to reduce dimension. In this paper, we focus on the convergence property of the discrete algorithm. To realize that, we view the algorithm as a stochastic process on the parameter space and semi-group. We approximate it by SDEs, and prove large time convergence of the SDEs to ensure its performance. This process is completed in three steps. First, the discrete algorithm can be viewed as a semigroup: $S^kφ=\mathbb{E}[φ(\mathbf W(k))]$. Second, we construct stochastic differential equations (SDEs) on the Stiefel manifold, i.e. the diffusion approximation, to approximate the semigroup. By proving the weak convergence, we verify that the algorithm is 'close to' the SDEs. Finally, we use the reversibility of the SDEs to prove long-time convergence.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Convergence of Oja's online principal component flow

Online principal component analysis (PCA) has been an efficient tool in practice to reduce dimension. However, convergence properties of the corresponding ODE are still unknown, including global convergence, stable manifolds, and convergence rate. In this paper, we focus on the stochastic gradient ascent (SGA) method proposed by Oja. By regarding the corresponding ODE as a Landau-Lifshitz-Gilbert (LLG) equation on the Stiefel manifold, we proved global convergence of the ODE. Moreover, we developed a new technique to determine stable manifolds. This technique analyzes the rank of the initial datum. Using this technique, we derived the explicit expression of the stable manifolds. As a consequence, exponential convergence to stable equilibrium points was also proved. The success of this new technique should be attributed to the semi-decoupling property of the SGA method: iteration of previous components does not depend on that of later ones. As far as we know, our result is the first complete one on the convergence of an online PCA flow, providing global convergence, explicit characterization of stable manifolds, and closed formula of exponential convergence depending on the spectrum gap.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Homogenization of the Landau-Lifshitz-Gilbert equation with natural boundary condition

The full Landau-Lifshitz-Gilbert equation with periodic material coefficients and natural boundary condition is employed to model the magnetization dynamics in composite ferromagnets. In this work, we establish the convergence between the homogenized solution and the original solution via a Lax equivalence theorem kind of argument. There are a few technical difficulties, including: 1) it is proven the classic choice of corrector to homogenization cannot provide the convergence result in the $H^1$ norm; 2) a boundary layer is induced due to the natural boundary condition; 3) the presence of stray field give rise to a multiscale potential problem. To keep the convergence rates near the boundary, we introduce the Neumann corrector with a high-order modification. Estimates on singular integral for disturbed functions and boundary layer are deduced, to conduct consistency analysis of stray field. Furthermore, inspired by length conservation of magnetization, we choose proper correctors in specific geometric space. These, together with a uniform $W^{1,6}$ estimate on original solution, provide the convergence rates in the $H^1$ sense.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On the equilibrium of the Poisson-Nernst-Planck-Bikermann model equipping with the steric and correlation effects

The Poisson-Nernst-Planck-Bikermann (PNPB) model, in which the ions and water molecules are treated as different species with non-uniform sizes and valences with interstitial voids, can describe the steric and correlation effects in ionic solution neglected by the Poisson-Nernst-Planck and Poisson-Boltzmann theories with point charge assumption. In the PNPB model, the electric potential is governed by the fourth-order Poisson-Bikermann (4PBik) equation instead of the Poisson equation so that it can describe the correlation effect. What's more, the steric potential is included in the ionic and water fluxes as well as the equilibrium Fermi-like distributions which characterizes the steric effect quantitatively. In this work, after doing a nondimensionalization step, we analyze the self-adjointness and the kernel of the fourth-order operator of the 4PBik equation. Also, we show the positivity of the void volume function and the convexity of the free energy. Following these properties, the well-posedness of the PNPB model in equilibrium is given. Furthermore, because the PNPB model has an energy dissipated structure, we adopt a finite volume scheme which preserves the energy dissipated property at the semi-discrete level. After that, various numerical investigations are given to show the parameter dependence of the steric effect to the steady state.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Some random batch particle methods for the Poisson-Nernst-Planck and Poisson-Boltzmann equations

We consider in this paper random batch interacting particle methods for solving the Poisson-Nernst-Planck (PNP) equations, and thus the Poisson-Boltzmann (PB) equation as the equilibrium, in the external unbounded domain. To justify the simulation in a truncated domain, an error estimate of the truncation is proved in the symmetric cases for the PB equation. Then, the random batch interacting particle methods are introduced which are $O(N)$ per time step. The particle methods can not only be considered as a numerical method for solving the PNP and PB equations, but also can be used as a direct simulation approach for the dynamics of the charged particles in solution. The particle methods are preferable due to their simplicity and adaptivity to complicated geometry, and may be interesting in describing the dynamics of the physical process. Moreover, it is feasible to incorporate more physical effects and interactions in the particle methods and to describe phenomena beyond the scope of the mean-field equations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Investigating the integrate and fire model as the limit of a random discharge model: a stochastic analysis perspective

In the mean field integrate-and-fire model, the dynamics of a typical neuron within a large network is modeled as a diffusion-jump stochastic process whose jump takes place once the voltage reaches a threshold. In this work, the main goal is to establish the convergence relationship between the regularized process and the original one where in the regularized process, the jump mechanism is replaced by a Poisson dynamic, and jump intensity within the classically forbidden domain goes to infinity as the regularization parameter vanishes. On the macroscopic level, the Fokker-Planck equation for the process with random discharges (i.e. Poisson jumps) are defined on the whole space, while the equation for the limit process is on the half space. However, with the iteration scheme, the difficulty due to the domain differences has been greatly mitigated and the convergence for the stochastic process and the firing rates can be established. Moreover, we find a polynomial-order convergence for the distribution by a re-normalization argument in probability theory. Finally, by numerical experiments, we quantitatively explore the rate and the asymptotic behavior of the convergence for both linear and nonlinear models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Class of Functional Inequalities and their Applications to Fourth-Order Nonlinear Parabolic Equations

We study a class of fourth order nonlinear parabolic equations which include the thin-film equation and the quantum drift-diffusion model as special cases. We investigate these equations by first developing functional inequalities of the type $ \int_Ωu^{2γ-α-β}Δu^αΔu^βdx \geq c\int_Ω|Δu^γ|^2dx $, which seem to be of interest on their own right.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Proximal-Gradient Algorithm for Crystal Surface Evolution

As a counterpoint to recent numerical methods for crystal surface evolution, which agree well with microscopic dynamics but suffer from significant stiffness that prevents simulation on fine spatial grids, we develop a new numerical method based on the macroscopic partial differential equation, leveraging its formal structure as the gradient flow of the total variation energy, with respect to a weighted $H^{-1}$ norm. This gradient flow structure relates to several metric space gradient flows of recent interest, including 2-Wasserstein flows and their generalizations to nonlinear mobilities. We develop a novel semi-implicit time discretization of the gradient flow, inspired by the classical minimizing movements scheme (known as the JKO scheme in the 2-Wasserstein case). We then use a primal dual hybrid gradient (PDHG) method to compute each element of the semi-implicit scheme. In one dimension, we prove convergence of the PDHG method to the semi-implicit scheme, under general integrability assumptions on the mobility and its reciprocal. Finally, by taking finite difference approximations of our PDHG method, we arrive at a fully discrete numerical algorithm, with iterations that converge at a rate independent of the spatial discretization: in particular, the convergence properties do not deteriorate as we refine our spatial grid. We close with several numerical examples illustrating the properties of our method, including facet formation at local maxima, pinning at local minima, and convergence as the spatial and temporal discretizations are refined.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Convergence of Random Batch Method for interacting particles with disparate species and weights

We consider in this work the convergence of Random Batch Method proposed in our previous work [Jin et al., J. Comput. Phys., 400(1), 2020] for interacting particles to the case of disparate species and weights. We show that the strong error is of $O(\sqrtτ)$ while the weak error is of $O(τ)$ where $τ$ is the time step between two random divisions of batches. Both types of convergence are uniform in $N$, the number of particles. The proof of strong convergence follows closely the proof in [Jin et al., J. Comput. Phys., 400(1), 2020] for indistinguishable particles, but there are still some differences: since there is no exchangeability now, we have to use a certain weighted average of the errors; some refined auxiliary lemmas have to be proved compared with our previous work. To show that the weak convergence of empirical measure is uniform in $N$, certain sharp estimates for the derivatives of the backward equations have been used. The weak convergence analysis is also illustrating for the convergence of Random Batch Method for $N$-body Liouville equations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the mean-field limit for the Vlasov-Poisson-Fokker-Planck system

We rigorously justify the mean-field limit of a $N$-particle system subject to the Brownian motion and interacting through a Newtonian potential in $\mathbb{R}^3$. Our result leads to a derivation of the Vlasov-Poisson-Fokkker-Planck (VPFP) equation from the microscopic $N$-particle system. More precisely, we show that the maximal distance between the exact microscopic trajectories and trajectories following the the mean-field is bounded by $N^{-\frac{1}{3}+\varepsilon}$ ($\frac{1}{63}\leq\varepsilon<\frac{1}{36}$) for a system with blob size $N^{-δ}$ ($\frac{1}{3}\leqδ<\frac{19}{54}-\frac{2\varepsilon}{3}$) up to a probability $1-N^{-α}$ for any $α>0$. Moreover, we prove the convergence rate between the empirical measure associated to the particle system and the solution of the VPFP equations. The technical novelty of this paper is that our estimates crucially rely on the randomness coming from the initial data and from the Brownian motion.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A Markov jump process modelling animal group size statistics

We translate a coagulation-framentation model, describing the dynamics of animal group size distributions, into a model for the population distribution and associate the \blue{nonlinear} evolution equation with a Markov jump process of a type introduced in classic work of H.~McKean. In particular this formalizes a model suggested by H.-S. Niwa [J.~Theo.~Biol.~224 (2003)] with simple coagulation and fragmentation rates. Based on the jump process, we develop a numerical scheme that allows us to approximate the equilibrium for the Niwa model, validated by comparison to analytical results by Degond et al. [J.~Nonlinear Sci.~27 (2017)], and study the population and size distributions for more complicated rates. Furthermore, the simulations are used to describe statistical properties of the underlying jump process. We additionally discuss the relation of the jump process to models expressed in stochastic differential equations and demonstrate that such a connection is justified in the case of nearest-neighbour interactions, as opposed to global interactions as in the Niwa model.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A stochastic version of Stein Variational Gradient Descent for efficient sampling

We propose in this work RBM-SVGD, a stochastic version of Stein Variational Gradient Descent (SVGD) method for efficiently sampling from a given probability measure and thus useful for Bayesian inference. The method is to apply the Random Batch Method (RBM) for interacting particle systems proposed by Jin et al to the interacting particle systems in SVGD. While keeping the behaviors of SVGD, it reduces the computational cost, especially when the interacting kernel has long range. Numerical examples verify the efficiency of this new version of SVGD.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Large time behavior, bi-Hamiltonian structure and kinetic formulation for complex Burgers equation

We prove the existence and uniqueness of positive analytical solutions with positive initial data to the mean field equation (the Dyson equation) of the Dyson Brownian motion through the complex Burgers equation with a force term on the upper half complex plane. These solutions converge to a steady state given by Wigner&#39;s semicircle law. A unique global weak solution with nonnegative initial data to the Dyson equation is obtained and some explicit solutions are given by Wigner&#39;s semicircle laws. We also construct a bi-Hamiltonian structure for the system of the real and imaginary components of the complex Burgers equation (coupled Burgers system). We establish a kinetic formulation for the coupled Burgers system and prove the existence and uniqueness of entropy solutions. The coupled Burgers system in Lagrangian variable naturally leads to two interacting particle systems: Fermi-Pasta-Ulam-Tsingou model with nearest-neighbor interactions, and Calogero-Moser model. These two particle systems yield the same Lagrangian dynamics in the continuum limit.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

On mean field limit for Brownian particles with Coulomb interaction in 3D

In this paper, we consider the mean field limit of Brownian particles with Coulomb interaction in 3D space. In particular, using a symmetrization technique, we show that the limit measure almost surely is a weak solution to the limiting nonlinear Fokker-Planck equation. By proving that the energy almost surely is bounded by the initial energy, we improve the regularity of the weak solutions. Moreover, by a natural assumption, we establish the weak strong uniqueness principle, which is closely related to the propagation of chaos.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Well-posedness and derivative blow-up for a dispersionless regularized shallow water system

We study local-time well-posedness and breakdown for solutions of regularized Saint-Venant equations (regularized classical shallow water equations) recently introduced by Clamond and Dutykh. The system is linearly non-dispersive, and smooth solutions conserve an $H^1$-equivalent energy. No shock discontinuities can occur, but the system is known to admit weakly singular shock-profile solutions that dissipate energy. We identify a class of small-energy smooth solutions that develop singularities in the first derivatives in finite time.

0 citations0 reviewsNo code linkedOpen access