BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jian-Guo Liu

Jian-Guo Liu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.APmath-phmath.MPphysics.soc-phInformation Retrievalmath.NAphysics.data-anSocial and Information Networksmath.PRNumerical Analysismath.DSphysics.flu-dynMachine LearningPopulations and Evolutionq-fin.GNcond-mat.mes-hallcond-mat.mtrl-scics.CYmath.CAmath.FAmath.OCnlin.AOnlin.CGphysics.plasm-ph

Catalog footprint

What is connected

55works
24topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

55 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'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'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
preprint2016arXiv

Continuum dynamics of the intention field under weakly cohesive social interactions

We investigate the long-time dynamics of an opinion formation model inspired by a work by Borghesi, Bouchaud and Jensen. Firstly, we derive a Fokker-Planck type equation under the assumption that interactions between individuals produce little consensus of opinion (grazing collision approximation). Secondly, we study conditions under which the Fokker-Planck equation has non-trivial equilibria and derive the macroscopic limit (corresponding to the long-time dynamics and spatially localized interactions) for the evolution of the mean opinion. Finally, we compare two different types of interaction rates: the original one given in the work of Borghesi, Bouchaud and Jensen (symmetric binary interactions) and one inspired from works by Motsch and Tadmor (non-symmetric binary interactions). We show that the first case leads to a conservative model for the density of the mean opinion whereas the second case leads to a non-conservative equation. We also show that the speed at which consensus is reached asymptotically for these two rates has fairly different density dependence.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Explicit and implicit TVD schemes for conservation laws with Caputo derivatives

In this paper, we investigate numerical approximations of the scalar conservation law with the Caputo derivative, which introduces the memory effect. We construct the first order and the second order explicit upwind schemes for such equations, which are shown to be conditionally $\ell^1$ contracting and TVD. However, the Caputo derivative leads to the modified CFL-type stability condition, $ (Δt)^α = O(Δx)$, where $α\in (0,1]$ is the fractional exponent in the derivative. When $α$ small, such strong constraint makes the numerical implementation extremely impractical. We have then proposed the implicit upwind scheme to overcome this issue, which is proved to be unconditionally $\ell^1$ contracting and TVD. Various numerical tests are presented to validate the properties of the methods and provide more numerical evidence in interpreting the memory effect in conservation laws.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Positivity-preserving and asymptotic preserving method for 2D Keller-Segal equations

We propose a semi-discrete scheme for 2D Keller-Segel equations based on a symmetrization reformation, which is equivalent to the convex splitting method and is free of any nonlinear solver. We show that, this new scheme is unconditionally stable as long as the initial condition does not exceed certain threshold, and it asymptotically preserves the quasi-static limit in the transient regime. Furthermore, we prove that the fully discrete scheme is conservative and positivity preserving, which makes it ideal for simulations. The analogical schemes for the radial symmetric cases and the subcritical degenerate cases are also presented and analyzed. With extensive numerical tests, we verify the claimed properties of the methods and demonstrate their superiority in various challenging applications.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An exact solution for Stokes flow in a channel with arbitrarily large wall permeability

We derive an exact solution for Stokes flow in an in a channel with permeable walls. We assume that at the channel walls, the normal component of the fluid velocity is described by Darcy's law and the tangential component of the fluid velocity is described by the no slip condition. The pressure exterior to the channel is assumed to be constant. Although this problem has been well studied, typical studies assume that the permeability of the wall is small relative to other non-dimensional parameters; this work relaxes this assumption and explores a regime in parameter space that has not yet been well studied. A consequence of this relaxation is that transverse velocity is no longer necessarily small when compared with the axial velocity. We use our result to explore how existing asymptotic theories break down in the limit of large permeability.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Coagulation-fragmentation model for animal group-size statistics

We study coagulation-fragmentation equations inspired by a simple model proposed in fisheries science to explain data for the size distribution of schools of pelagic fish. Although the equations lack detailed balance and admit no $H$-theorem, we are able to develop a rather complete description of equilibrium profiles and large-time behavior, based on recent developments in complex function theory for Bernstein and Pick functions. In the large-population continuum limit, a scaling-invariant regime is reached in which all equilibria are determined by a single scaling profile. This universal profile exhibits power-law behavior crossing over from exponent $-\frac23$ for small size to $-\frac32$ for large size, with an exponential cut-off.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Convergence of Stochastic Interacting Particle Systems in Probability under a Sobolev Norm

In this paper, we consider particle systems with interaction and Brownian motion. We prove that when the initial data is from the sampling of Chorin's method, i.e., the initial vertices are on lattice points $hi\in \mathbb{R}^d$ with mass $ρ_0(hi) h^d$, where $ρ_0$ is some initial density function, then the regularized empirical measure of the interacting particle system converges in probability to the corresponding mean-field partial differential equation with initial density $ρ_0$, under the Sobolev norm of $L^\infty(L^2)\cap L^2(H^1)$. Our result is true for all those systems when the interacting function is bounded, Lipschitz continuous and satisfies certain regular condition. And if we further regularize the interacting particle system, it also holds for some of the most important systems of which the interacting functions are not. For systems with repulsive Coulomb interaction, this convergence holds globally on any interval $[0,t]$. And for systems with attractive Newton force as interacting function, we have convergence within the largest existence time of the regular solution of the corresponding Keller-Segel equation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Emergence of step flow from atomistic scheme of epitaxial growth in 1+1 dimensions

The Burton-Cabrera-Frank (BCF) model for the flow of line defects (steps) on crystal surfaces has offered useful insights into nanostructure evolution. This model has rested on phenomenological grounds. Our goal is to show via scaling arguments the emergence of the BCF theory for non-interacting steps from a stochastic atomistic scheme of a simplified kinetic solid-on-solid model in one spatial dimension. Our main assumptions are: adsorbed atoms (adatoms) form a dilute system, and elastic effects of the crystal lattice are absent. The step edge is treated as a front that propagates via probabilistic rules for atom attachment and detachment at the step. We formally derive a quasistatic step flow description by averaging out the stochastic scheme when terrace diffusion, adatom desorption and deposition from above are present.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Empirical Analysis of the Online Rating Systems

This paper is to analyze the properties of evolving bipartite networks from four aspects, the growth of networks, the degree distribution, the popularity of objects and the diversity of user behaviours, leading a deep understanding on the empirical data. By empirical studies of data from the online bookstore Amazon and a question and answer site Stack Overflow, which are both rating bipartite networks, we could reveal the rules for the evolution of bipartite networks. These rules have significant meanings in practice for maintaining the operation of real systems and preparing for their future development. We find that the degree distribution of users follows a power law with an exponential cutoff. Also, according to the evolution of popularity for objects, we find that the large-degree objects tend to receive more new ratings than expected depending on their current degrees while the small-degree objects receive less ratings in terms of their degrees. Moreover, the user behaviours show such a trend that the larger degree the users have, the stronger purposes are with their behaviours except the initial periods when users choose a diversity of products to learn about what they want. Finally, we conclude with a discussion on how the bipartite network evolves, which provides guideline for meeting challenges brought by the growth of network.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Iterative resource allocation based on propagation feature of node for identifying the influential nodes

The Identification of the influential nodes in networks is one of the most promising domains. In this paper, we present an improved iterative resource allocation (IIRA) method by considering the centrality information of neighbors and the influence of spreading rate for a target node. Comparing with the results of the Susceptible Infected Recovered (SIR) model for four real networks, the IIRA method could identify influential nodes more accurately than the tradition IRA method. Specially, in the Erdos network, the Kendall's tau could be enhanced 23\% when the spreading rate is 0.12. In the Protein network, the Kendall's tau could be enhanced 24\% when the spreading rate is 0.08.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Locating influential nodes via dynamics-sensitive centrality

With great theoretical and practical significance, locating influential nodes of complex networks is a promising issues. In this paper, we propose a dynamics-sensitive (DS) centrality that integrates topological features and dynamical properties. The DS centrality can be directly applied in locating influential spreaders. According to the empirical results on four real networks for both susceptible-infected-recovered (SIR) and susceptible-infected (SI) spreading models, the DS centrality is much more accurate than degree, $k$-shell index and eigenvector centrality.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Optimal reservoir conditions for fluid extraction through permeable walls in the viscous limit

In biological transport mechanisms such as insect respiration and renal filtration, fluid travels along a leaky channel allowing exchange with systems exterior the the channel. The channels in these systems may undergo peristaltic pumping which is thought to enhance the material exchange. To date, little analytic work has been done to study the effect of pumping on material extraction across the channel walls. In this paper, we examine a fluid extraction model in which fluid flowing through a leaky channel is exchanged with fluid in a reservoir. The channel walls are allowed to contract and expand uniformly, simulating a pumping mechanism. In order to efficiently determine solutions of the model, we derive a formal power series solution for the Stokes equations in a finite channel with uniformly contracting/expanding permeable walls. This flow has been well studied in the case of weakly permeable channel walls in which the normal velocity at the channel walls is proportional to the wall velocity. In contrast we do not assume weakly driven flow, but flow driven by hydrostatic pressure, and we use Dacry's law to close our system for normal wall velocity. We use our flow solution to examine flux across the channel-reservoir barrier and demonstrate that pumping can either enhance or impede fluid extraction across channel walls. We find that associated with each set of physical flow and pumping parameters, there are optimal reservoir conditions that maximizes the amount of material flowing from the channel into the reservoir.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stability of similarity measurements for bipartite networks

Similarity is a fundamental measure in network analyses and machine learning algorithms, with wide applications ranging from personalized recommendation to socio-economic dynamics. We argue that an effective similarity measurement should guarantee the stability even under some information loss. With six bipartite networks, we investigate the stabilities of fifteen similarity measurements by comparing the similarity matrixes of two data samples which are randomly divided from original data sets. Results show that, the fifteen measurements can be well classified into three clusters according to their stabilities, and measurements in the same cluster have similar mathematical definitions. In addition, we develop a top-$n$-stability method for personalized recommendation, and find that the unstable similarities would recommend false information to users, and the performance of recommendation would be largely improved by using stable similarity measurements. This work provides a novel dimension to analyze and evaluate similarity measurements, which can further find applications in link prediction, personalized recommendation, clustering algorithms, community detection and so on.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Evolution of wealth in a nonconservative economy driven by local Nash equilibria

We develop a model for the evolution of wealth in a non-conservative economic environment, extending a theory developed earlier by the authors. The model considers a system of rational agents interacting in a game theoretical framework. This evolution drives the dynamic of the agents in both wealth and economic configuration variables. The cost function is chosen to represent a risk averse strategy of each agent. That is, the agent is more likely to interact with the market, the more predictable the market, and therefore the smaller its individual risk. This yields a kinetic equation for an effective single particle agent density with a Nash equilibrium serving as the local thermodynamic equilibrium. We consider a regime of scale separation where the large scale dynamics is given by a hydrodynamic closure with this local equilibrium. A class of generalized collision invariants (GCIs) is developed to overcome the difficulty of the non-conservative property in the hydrodynamic closure derivation of the large scale dynamics for the evolution of wealth distribution. The result is a system of gas dynamics-type equations for the density and average wealth of the agents on large scales. We recover the inverse Gamma distribution, which has been previously considered in the literature, as a local equilibrium for particular choices of the cost function.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Meanfield games and model predictive control

Mean-Field Games are games with a continuum of players that incorporate the time-dimension through a control-theoretic approach. Recently, simpler approaches relying on the Best Reply Strategy have been proposed. They assume that the agents navigate their strategies towards their goal by taking the direction of steepest descent of their cost function (i.e. the opposite of the utility function). In this paper, we explore the link between Mean-Field Games and the Best Reply Strategy approach. This is done by introducing a Model Predictive Control framework, which consists of setting the Mean-Field Game over a short time interval which recedes as time moves on. We show that the Model Predictive Control offers a compromise between a possibly unrealistic Mean-Field Game approach and the sub-optimal Best Reply Strategy.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Memory effect of the online user preference

The mechanism of the online user preference evolution is of great significance for understanding the online user behaviors and improving the quality of online services. Since users are allowed to rate on objects in many online systems, ratings can well reflect the users' preference. With two benchmark datasets from online systems, we uncover the memory effect in users' selecting behavior which is the sequence of qualities of selected objects and the rating behavior which is the sequence of ratings delivered by each user. Furthermore, the memory duration is presented to describe the length of a memory, which exhibits the power-law distribution, i.e., the probability of the occurring of long-duration memory is much higher than that of the random case which follows the exponential distribution. We present a preference model in which a Markovian process is utilized to describe the users' selecting behavior, and the rating behavior depends on the selecting behavior. With only one parameter for each of the user's selecting and rating behavior, the preference model could regenerate any duration distribution ranging from the power-law form (strong memory) to the exponential form (weak memory).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Phase transitions, hysteresis, and hyperbolicity for self-organized alignment dynamics

We provide a complete and rigorous description of phase transitions for kinetic models of self-propelled particles interacting through alignment. These models exhibit a competition between alignment and noise. Both the alignment frequency and noise intensity depend on a measure of the local alignment. We show that, in the spatially homogeneous case, the phase transition features (number and nature of equilibria, stability, convergence rate, phase diagram, hysteresis) are totally encoded in how the ratio between the alignment and noise intensities depend on the local alignment. In the spatially inhomogeneous case, we derive the macroscopic models associated to the stable equilibria and classify their hyperbolicity according to the same function.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Ranking the spreading influence in complex networks

Identifying the node spreading influence in networks is an important task to optimally use the network structure and ensure the more efficient spreading in information. In this paper, by taking into account the shortest distance between a target node and the node set with the highest $k$-core value, we present an improved method to generate the ranking list to evaluate the node spreading influence. Comparing with the epidemic process results for four real networks and the Barabási-Albert network, the parameterless method could identify the node spreading influence more accurately than the ones generated by the degree $k$, closeness centrality, $k$-shell and mixed degree decomposition methods. This work would be helpful for deeply understanding the node importance of a network.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Ultra accurate collaborative information filtering via directed user similarity

A key challenge of the collaborative filtering (CF) information filtering is how to obtain the reliable and accurate results with the help of peers' recommendation. Since the similarities from small-degree users to large-degree users would be larger than the ones opposite direction, the large-degree users' selections are recommended extensively by the traditional second-order CF algorithms. By considering the users' similarity direction and the second-order correlations to depress the influence of mainstream preferences, we present the directed second-order CF (HDCF) algorithm specifically to address the challenge of accuracy and diversity of the CF algorithm. The numerical results for two benchmark data sets, MovieLens and Netflix, show that the accuracy of the new algorithm outperforms the state-of-the-art CF algorithms. Comparing with the CF algorithm based on random-walks proposed in the Ref.7, the average ranking score could reach 0.0767 and 0.0402, which is enhanced by 27.3\% and 19.1\% for MovieLens and Netflix respectively. In addition, the diversity, precision and recall are also enhanced greatly. Without relying on any context-specific information, tuning the similarity direction of CF algorithms could obtain accurate and diverse recommendations. This work suggests that the user similarity direction is an important factor to improve the personalized recommendation performance.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A note on Aubin-Lions-Dubinskii lemmas

Strong compactness results for families of functions in seminormed nonnegative cones in the spirit of the Aubin-Lions-Dubinski\uı lemma are proven, refining some recent results in the literature. The first theorem sharpens slightly a result of Dubinski\uı (1965) for seminormed cones. The second theorem applies to piecewise constant functions in time and sharpens slightly the results of Dreher and Jüngel (2012) and Chen and Liu (2012). An application is given, which is useful in the study of porous-medium or fast-diffusion type equations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Evolution of the distribution of wealth in an economic environment driven by local Nash equilibria

We present and analyze a model for the evolution of the wealth distribution within a heterogeneous economic environment. The model considers a system of rational agents interacting in a game theoretical framework, through fairly general assumptions on the cost function. This evolution drives the dynamic of the agents in both wealth and economic configuration variables. We consider a regime of scale separation where the large scale dynamics is given by a hydrodynamic closure with a Nash equilibrium serving as the local thermodynamic equilibrium. The result is a system of gas dynamics-type equations for the density and average wealth of the agents on large scales. We recover the inverse gamma distribution as an equilibrium in the particular case of quadratic cost functions which has been previously considered in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Flow on sweeping networks

We introduce a cellular automaton model coupled with a transport equation for flows on graphs. The direction of the flow is described by a switching process where the switching probability dynamically changes according to the value of the transported quantity in the neighboring cells. A motivation is pedestrian dynamics in a small corridor where the propagation of people in a part of the corridor can be either left or rightgoing. Under the assumptions of propagation of chaos and mean-field limit, we derive a master equation and the corresponding meanfield kinetic and macroscopic models. Steady--states are computed and analyzed analytically and exhibit the possibility of multiple meta-stable states and hysteresis.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Macroscopic models of collective motion and self-organization

In this paper, we review recent developments on the derivation and properties of macroscopic models of collective motion and self-organization. The starting point is a model of self-propelled particles interacting with its neighbors through alignment. We successively derive a mean-field model and its hydrodynamic limit. The resulting macroscopic model is the Self-Organized Hydrodynamics (SOH). We review the available existence results and known properties of the SOH model and discuss it in view of its possible extensions to other kinds of collective motion.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Nash equilibrium macroscopic closure for kinetic models coupled with Mean-Field Games

We introduce a new mean field kinetic model for systems of rational agents interacting in a game theoretical framework. This model is inspired from non-cooperative anonymous games with a continuum of players and Mean-Field Games. The large time behavior of the system is given by a macroscopic closure with a Nash equilibrium serving as the local thermodynamic equilibrium. An application of the presented theory to a social model (herding behavior) is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A note on phase transitions for the Smoluchowski equation with dipolar potential

In this note, we study the phase transitions arising in a modified Smoluchowski equation on the sphere with dipolar potential. This equation models the competition between alignment and diffusion, and the modification consists in taking the strength of alignment and the intensity of the diffusion as functions of the order parameter. We characterize the stable and unstable equilibrium states. For stable equilibria, we provide the exponential rate of convergence. We detail special cases, giving rise to second order and first order phase transitions, respectively. We study the hysteresis diagram, and provide numerical illustrations of this phenomena.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Solving the accuracy-diversity dilemma via directed random walks

Random walks have been successfully used to measure user or object similarities in collaborative filtering (CF) recommender systems, which is of high accuracy but low diversity. A key challenge of CF system is that the reliably accurate results are obtained with the help of peers' recommendation, but the most useful individual recommendations are hard to be found among diverse niche objects. In this paper we investigate the direction effect of the random walk on user similarity measurements and find that the user similarity, calculated by directed random walks, is reverse to the initial node's degree. Since the ratio of small-degree users to large-degree users is very large in real data sets, the large-degree users' selections are recommended extensively by traditional CF algorithms. By tuning the user similarity direction from neighbors to the target user, we introduce a new algorithm specifically to address the challenge of diversity of CF and show how it can be used to solve the accuracy-diversity dilemma. Without relying on any context-specific information, we are able to obtain accurate and diverse recommendations, which outperforms the state-of-the-art CF methods. This work suggests that the random walk direction is an important factor to improve the personalized recommendation performance.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Blow-up, zero $α$ limit and the Liouville type theorem for the Euler-Poincaré equations

In this paper we study the Euler-Poincaré equations in $\Bbb R^N$. We prove local existence of weak solutions in $W^{2,p}(\Bbb R^N),$ $p>N$, and local existence of unique classical solutions in $H^k (\Bbb R^N)$, $k>N/2+3$, as well as a blow-up criterion. For the zero dispersion equation($α=0$) we prove a finite time blow-up of the classical solution. We also prove that as the dispersion parameter vanishes, the weak solution converges to a solution of the zero dispersion equation with sharp rate as $α\to0$, provided that the limiting solution belongs to $C([0, T);H^k(\Bbb R^N))$ with $k>N/2 +3$. For the {\em stationary weak solutions} of the Euler-Poincaré equations we prove a Liouville type theorem. Namely, for $α>0$ any weak solution $\mathbf{u}\in H^1(\Bbb R^N)$ is $\mathbf{u}=0$; for $α=0$ any weak solution $\mathbf{u}\in L^2(\Bbb R^N)$ is $\mathbf{u}=0$.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dynamics in a kinetic model of oriented particles with phase transition

Motivated by a phenomenon of phase transition in a model of alignment of self-propelled particles, we obtain a kinetic mean-field equation which is nothing else than the Doi equation (also called Smoluchowski equation) with dipolar potential. In a self-contained article, using only basic tools, we analyze the dynamics of this equation in any dimension. We first prove global well-posedness of this equation, starting with an initial condition in any Sobolev space. We then compute all possible steady-states. There is a threshold for the noise parameter: over this threshold, the only equilibrium is the uniform distribution, and under this threshold, there is also a family of non-isotropic equilibria. We give a rigorous prove of convergence of the solution to a steady-state as time goes to infinity. In particular we show that in the supercritical case, the only initial conditions leading to the uniform distribution in large time are those with vanishing momentum. For any positive value of the noise parameter, and any initial condition, we give rates of convergence towards equilibrium, exponentially for both supercritical and subcritical cases and algebraically for the critical case.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hydrodynamic models of self-organized dynamics: derivation and existence theory

This paper is concerned with the derivation and analysis of hydrodynamic models for systems of self-propelled particles subject to alignment interaction and attraction-repulsion. The starting point is the kinetic model considered in earlier work of Degond & Motsch with the addition of an attraction-repulsion interaction potential. Introducing different scalings than in Degond & Motsch, the non-local effects of the alignment and attraction-repulsion interactions can be kept in the hydrodynamic limit and result in extra pressure, viscosity terms and capillary force. The systems are shown to be symmetrizable hyperbolic systems with viscosity terms. A local-in-time existence result is proved in the 2D case for the viscous model and in the 3D case for the inviscid model. The proof relies on the energy method.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hydrodynamics of self-alignment interactions with precession and derivation of the Landau-Lifschitz-Gilbert equation

We consider a kinetic model of self-propelled particles with alignment interaction and with precession about the alignment direction. We derive a hydrodynamic system for the local density and velocity orientation of the particles. The system consists of the conservative equation for the local density and a non-conservative equation for the orientation. First, we assume that the alignment interaction is purely local and derive a first order system. However, we show that this system may lose its hyperbolicity. Under the assumption of weakly non-local interaction, we derive diffusive corrections to the first order system which lead to the combination of a heat flow of the harmonic map and Landau-Lifschitz-Gilbert dynamics. In the particular case of zero self-propelling speed, the resulting model reduces to the phenomenological Landau-Lifschitz-Gilbert equations. Therefore the present theory provides a kinetic formulation of classical micromagnetization models and spin dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Information filtering via biased heat conduction

Heat conduction process has recently found its application in personalized recommendation [T. Zhou \emph{et al.}, PNAS 107, 4511 (2010)], which is of high diversity but low accuracy. By decreasing the temperatures of small-degree objects, we present an improved algorithm, called biased heat conduction (BHC), which could simultaneously enhance the accuracy and diversity. Extensive experimental analyses demonstrate that the accuracy on MovieLens, Netflix and Delicious datasets could be improved by 43.5%, 55.4% and 19.2% compared with the standard heat conduction algorithm, and the diversity is also increased or approximately unchanged. Further statistical analyses suggest that the present algorithm could simultaneously identify users' mainstream and special tastes, resulting in better performance than the standard heat conduction algorithm. This work provides a creditable way for highly efficient information filtering.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Information Filtering via Implicit Trust-based Network

Based on the user-item bipartite network, collaborative filtering (CF) recommender systems predict users' interests according to their history collections, which is a promising way to solve the information exploration problem. However, CF algorithm encounters cold start and sparsity problems. The trust-based CF algorithm is implemented by collecting the users' trust statements, which is time-consuming and must use users' private friendship information. In this paper, we present a novel measurement to calculate users' implicit trust-based correlation by taking into account their average ratings, rating ranges, and the number of common rated items. By applying the similar idea to the items, a item-based CF algorithm is constructed. The simulation results on three benchmark data sets show that the performances of both user-based and item-based algorithms could be enhanced greatly. Finally, a hybrid algorithm is constructed by integrating the user-based and item-based algorithms, the simulation results indicate that hybrid algorithm outperforms the state-of-the-art methods. Specifically, it can not only provide more accurate recommendations, but also alleviate the cold start problem.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Kinetic formulation and global existence for the Hall-Magneto-hydrodynamics system

This paper deals with the derivation and analysis of the the Hall Magneto-Hydrodynamic equations. We first provide a derivation of this system from a two-fluids Euler-Maxwell system for electrons and ions, through a set of scaling limits. We also propose a kinetic formulation for the Hall-MHD equations which contains as fluid closure different variants of the Hall-MHD model. Then, we prove the existence of global weak solutions for the incompressible viscous resistive Hall-MHD model. We use the particular structure of the Hall term which has zero contribution to the energy identity. Finally, we discuss particular solutions in the form of axisymmetric purely swirling magnetic fields and propose some regularization of the Hall equation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Macroscopic limits and phase transition in a system of self-propelled particles

We investigate systems of self-propelled particles with alignment interaction. Compared to previous work, the force acting on the particles is not normalized and this modification gives rise to phase transitions from disordered states at low density to aligned states at high densities. This model is the space inhomogeneous extension of a previous work by Frouvelle and Liu in which the existence and stability of the equilibrium states were investigated. When the density is lower than a threshold value, the dynamics is described by a non-linear diffusion equation. By contrast, when the density is larger than this threshold value, the dynamics is described by a hydrodynamic model for self-alignment interactions previously derived in Degond and Motsch. However, the modified normalization of the force gives rise to different convection speeds and the resulting model may lose its hyperbolicity in some regions of the state space.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Solving the apparent diversity-accuracy dilemma of recommender systems

Recommender systems use data on past user preferences to predict possible future likes and interests. A key challenge is that while the most useful individual recommendations are to be found among diverse niche objects, the most reliably accurate results are obtained by methods that recommend objects based on user or object similarity. In this paper we introduce a new algorithm specifically to address the challenge of diversity and show how it can be used to resolve this apparent dilemma when combined in an elegant hybrid with an accuracy-focused algorithm. By tuning the hybrid appropriately we are able to obtain, without relying on any semantic or context-specific information, simultaneous gains in both accuracy and diversity of recommendations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Degree correlation effect of bipartite network on personalized recommendation

In this paper, by introducing a new user similarity index base on the diffusion process, we propose a modified collaborative filtering (MCF) algorithm, which has remarkably higher accuracy than the standard collaborative filtering. In the proposed algorithm, the degree correlation between users and objects is taken into account and embedded into the similarity index by a tunable parameter. The numerical simulation on a benchmark data set shows that the algorithmic accuracy of the MCF, measured by the average ranking score, is further improved by 18.19% in the optimal case. In addition, two significant criteria of algorithmic performance, diversity and popularity, are also taken into account. Numerical results show that the presented algorithm can provide more diverse and less popular recommendations, for example, when the recommendation list contains 10 objects, the diversity, measured by the hamming distance, is improved by 21.90%.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Effect of user tastes on personalized recommendation

In this paper, based on a weighted projection of the user-object bipartite network, we study the effects of user tastes on the mass-diffusion-based personalized recommendation algorithm, where a user's tastes or interests are defined by the average degree of the objects he has collected. We argue that the initial recommendation power located on the objects should be determined by both of their degree and the users' tastes. By introducing a tunable parameter, the user taste effects on the configuration of initial recommendation power distribution are investigated. The numerical results indicate that the presented algorithm could improve the accuracy, measured by the average ranking score, more importantly, we find that when the data is sparse, the algorithm should give more recommendation power to the objects whose degrees are close to the users' tastes, while when the data becomes dense, it should assign more power on the objects whose degrees are significantly different from user's tastes.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Highly accurate recommendation algorithm based on high-order similarities

In this Letter, we introduce a modified collaborative filtering (MCF) algorithm, which has remarkably higher accuracy than the standard collaborative filtering. In the MCF, instead of the standard Pearson coefficient, the user-user similarities are obtained by a diffusion process. Furthermore, by considering the second order similarities, we design an effective algorithm that depresses the influence of mainstream preferences. The corresponding algorithmic accuracy, measured by the ranking score, is further improved by 24.9% in the optimal case. In addition, two significant criteria of algorithmic performance, diversity and popularity, are also taken into account. Numerical results show that the algorithm based on second order similarity can outperform the MCF simultaneously in all three criteria.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Improved Collaborative Filtering Algorithm via Information Transformation

In this paper, we propose a spreading activation approach for collaborative filtering (SA-CF). By using the opinion spreading process, the similarity between any users can be obtained. The algorithm has remarkably higher accuracy than the standard collaborative filtering (CF) using Pearson correlation. Furthermore, we introduce a free parameter $β$ to regulate the contributions of objects to user-user correlations. The numerical results indicate that decreasing the influence of popular objects can further improve the algorithmic accuracy and personality. We argue that a better algorithm should simultaneously require less computation and generate higher accuracy. Accordingly, we further propose an algorithm involving only the top-$N$ similar neighbors for each target user, which has both less computational complexity and higher algorithmic accuracy.

0 citations0 reviewsNo code linkedOpen access