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

Pingwen Zhang

Pingwen Zhang appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.APcond-mat.softmath.NANumerical Analysisphysics.comp-phmath-phmath.MPApplicationsastro-ph.GAcond-mat.mtrl-scicond-mat.stat-mechphysics.chem-phPopulations and Evolution

Catalog footprint

What is connected

30works
13topics
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

30 published item(s)

preprint2024arXiv

Accurately recover global quasiperiodic systems by finite points

Quasiperiodic systems, related to irrational numbers, are space-filling structures without decay nor translation invariance. How to accurately recover these systems, especially for non-smooth cases, presents a big challenge in numerical computation. In this paper, we propose a new algorithm, finite points recovery (FPR) method, which is available for both smooth and non-smooth cases, to address this challenge. The FPR method first establishes a homomorphism between the lower-dimensional definition domain of the quasiperiodic function and the higher-dimensional torus, then recovers the global quasiperiodic system by employing interpolation technique with finite points in the definition domain without dimensional lifting. Furthermore, we develop accurate and efficient strategies of selecting finite points according to the arithmetic properties of irrational numbers. The corresponding mathematical theory, convergence analysis, and computational complexity analysis on choosing finite points are presented. Numerical experiments demonstrate the effectiveness and superiority of FPR approach in recovering both smooth quasiperiodic functions and piecewise constant Fibonacci quasicrystals. While existing spectral methods encounter difficulties in accurately recovering non-smooth quasiperiodic functions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A model-free shrinking-dimer saddle dynamics for finding saddle point and solution landscape

We propose a model-free shrinking-dimer saddle dynamics for finding any-index saddle points and constructing the solution landscapes, in which the force in the standard saddle dynamics is replaced by a surrogate model trained by the Gassian process learning. By this means, the exact form of the model is no longer necessary such that the saddle dynamics could be implemented based only on some observations of the force. This data-driven approach not only avoids the modeling procedure that could be difficult or inaccurate, but also significantly reduces the number of queries of the force that may be expensive or time-consuming. We accordingly develop a sequential learning saddle dynamics algorithm to perform a sequence of local saddle dynamics, in which the queries of the training samples and the update or retraining of the surrogate force are performed online and around the latent trajectory in order to improve the accuracy of the surrogate model and the value of each sampling. Numerical experiments are performed to demonstrate the effectiveness and efficiency of the proposed algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Discretization and index-robust error analysis for constrained high-index saddle dynamics on high-dimensional sphere

We develop and analyze numerical discretization to the constrained high-index saddle dynamics, the dynamics searching for the high-index saddle points confined on the high-dimensional unit sphere. Compared with the saddle dynamics without constraints, the constrained high-index saddle dynamics has more complex dynamical forms, and additional operations such as the retraction and vector transport are required due to the constraint, which significantly complicate the numerical scheme and the corresponding numerical analysis. Furthermore, as the existing numerical analysis results usually depend on the index of the saddle points implicitly, the proved numerical accuracy may be reduced if the index is high in many applications, which indicates the lack of robustness with respect to the index. To address these issues, we derive the error estimates for numerical discretization of the constrained high-index saddle dynamics on high-dimensional sphere, and then improve it by providing an index-robust error analysis in an averaged norm by adjusting the relaxation parameters. The developed results provide mathematical supports for the accuracy of numerical computations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Error estimates for Euler discretization of high-index saddle dynamics

High-index saddle dynamics provides an effective means to compute the any-index saddle points and construct the solution landscape. In this paper we prove error estimates for Euler discretization of high-index saddle dynamics with respect to the time step size, which remains untreated in the literature. We overcome the main difficulties that lie in the strong nonlinearity of the saddle dynamics and the orthonormalization procedure in the numerical scheme that is uncommon in standard discretization of differential equations. The derived methods are further extended to study the generalized high-index saddle dynamics for non-gradient systems and provide theoretical support for the accuracy of numerical implementations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Mathematical and numerical analysis to shrinking-dimer saddle dynamics with local Lipschitz conditions

We present a mathematical and numerical investigation to the shrinkingdimer saddle dynamics for finding any-index saddle points in the solution landscape. Due to the dimer approximation of Hessian in saddle dynamics, the local Lipschitz assumptions and the strong nonlinearity for the saddle dynamics, it remains challenges for delicate analysis, such as the the boundedness of the solutions and the dimer error. We address these issues to bound the solutions under proper relaxation parameters, based on which we prove the error estimates for numerical discretization to the shrinking-dimer saddle dynamics by matching the dimer length and the time step size. Furthermore, the Richardson extrapolation is employed to obtain a high-order approximation. The inherent reason of requiring the matching of the dimer length and the time step size lies in that the former serves a different mesh size from the later, and thus the proposed numerical method is close to a fully-discrete numerical scheme of some spacetime PDE model with the Hessian in the saddle dynamics and its dimer approximation serving as a "spatial operator" and its discretization, respectively, which in turn indicates the PDE nature of the saddle dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Computing solution landscape of nonlinear space-fractional problems via fast approximation algorithm

The nonlinear space-fractional problems often allow multiple stationary solutions, which can be much more complicated than the corresponding integer-order problems. In this paper, we systematically compute the solution landscapes of nonlinear constant/variable-order space-fractional problems. A fast approximation algorithm is developed to deal with the variable-order spectral fractional Laplacian by approximating the variable-indexing Fourier modes, and then combined with saddle dynamics to construct the solution landscape of variable-order space-fractional phase field model. Numerical experiments are performed to substantiate the accuracy and efficiency of fast approximation algorithm and elucidate essential features of the stationary solutions of space-fractional phase field model. Furthermore, we demonstrate that the solution landscapes of spectral fractional Laplacian problems can be reconfigured by varying the diffusion coefficients in the corresponding integer-order problems.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Transition pathways connecting crystals and quasicrystals

Due to structural incommensurability, the emergence of a quasicrystal from a crystalline phase represents a challenge to computational physics. Here the nucleation of quasicrystals is investigated by using an efficient computational method applied to a Landau free-energy functional. Specifically, transition pathways connecting different local minima of the Lifshitz-Petrich model are obtained by using the high-index saddle dynamics. Saddle points on these paths are identified as the critical nuclei of the 6-fold crystals and 12-fold quasicrystals. The results reveal that phase transitions between the crystalline and quasicrystalline phases could follow two possible pathways, corresponding to a one-stage phase transition and a two-stage phase transition involving a metastable lamellar quasicrystalline state, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Detecting Suspected Epidemic Cases Using Trajectory Big Data

Emerging infectious diseases are existential threats to human health and global stability. The recent outbreaks of the novel coronavirus COVID-19 have rapidly formed a global pandemic, causing hundreds of thousands of infections and huge economic loss. The WHO declares that more precise measures to track, detect and isolate infected people are among the most effective means to quickly contain the outbreak. Based on trajectory provided by the big data and the mean field theory, we establish an aggregated risk mean field that contains information of all risk-spreading particles by proposing a spatio-temporal model named HiRES risk map. It has dynamic fine spatial resolution and high computation efficiency enabling fast update. We then propose an objective individual epidemic risk scoring model named HiRES-p based on HiRES risk maps, and use it to develop statistical inference and machine learning methods for detecting suspected epidemic-infected individuals. We conduct numerical experiments by applying the proposed methods to study the early outbreak of COVID-19 in China. Results show that the HiRES risk map has strong ability in capturing global trend and local variability of the epidemic risk, thus can be applied to monitor epidemic risk at country, province, city and community levels, as well as at specific high-risk locations such as hospital and station. HiRES-p score seems to be an effective measurement of personal epidemic risk. The accuracy of both detecting methods are above 90\% when the population infection rate is under 20\%, which indicates great application potential in epidemic risk prevention and control practice.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Nonlocal Diffusion Operators for Normal and Anomalous Dynamics

The Laplacian $Δ$ is the infinitesimal generator of isotropic Brownian motion, being the limit process of normal diffusion, while the fractional Laplacian $Δ^{β/2}$ serves as the infinitesimal generator of the limit process of isotropic Lévy process. Taking limit, in some sense, means that the operators can approximate the physical process well after sufficient long time. We introduce the nonlocal operators (being effective from the starting time), which describe the general processes undergoing normal diffusion. For anomalous diffusion, we extend to the anisotropic fractional Laplacian $Δ_m^{β/2}$ and the tempered one $Δ_m^{β/2,λ}$ in $\mathbb{R}^n$. Their definitions are proved to be equivalent to an alternative one in Fourier space. Based on these new nonlocal diffusion operators, we further derive the deterministic governing equations of some interesting statistical observables of the very general jump processes with multiple internal states. Finally, we consider the associated initial and boundary value problems and prove their well-posedness of the Galerkin weak formulation in $\mathbb{R}^n$. To obtain the coercivity, we claim that the probability density function $m(Y)$ should be nondegenerate.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Application of Self-Consistent Field Theory to Self-Assembled Bilayer Membranes

Bilayer membranes self-assembled from amphiphilic molecules such as lipids, surfactants and block copolymers are ubiquitous in biological and physiochemical systems. The shape and structure of bilayer membranes depend crucially on their mechanical properties such surface tension, bending moduli and line tension. Understanding how the molecular property of the amphiphiles determine the structure and mechanics of the self-assembled bilayers requires a molecularly detailed theoretical framework. The self-consistent field theory provides such a theoretical framework, which is capable of accurately predicting mechanical parameters of self-assembled bilayer membranes. In this mini review we summarize the formulation of the self-consistent field theory, as exemplified by a model system composed of flexible amphiphilic chains dissolved in hydrophilic polymeric solvents, and its application to the study of self-assembled bilayer membranes.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On a molecular based Q-tensor model for liquid crystals with density variations

In this article, we study the new Q-tensor model previously derived from Onsager's molecular theory by Han \textit{et al.} [Arch. Rational Mech. Anal., 215.3 (2014), pp. 741-809] for static liquid crystal modeling. Taking density and Q-tensor as order parameters, the new Q-tensor model not only characterizes important phases while capturing density variation effects, but also remains computationally tractable and efficient. We report the results of two numerical applications of the model, namely the isotropic--nematic--smectic-A--smectic-C phase transitions and the isotropic--nematic interface problem, in which density variations are indispensable. Meanwhile, we show the connections of the new Q-tensor model with classical models including generalized Landau-de Gennes models, generalized McMillan models, and the Chen-Lubensky model. The new Q-tensor model is the pivot and an appropriate trade-off between the classical models in three scales.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stability of Soft Quasicrystals in a Coupled-Mode Swift-Hohenberg Model for Three-Component Systems

In this article, we discuss the stability of soft quasicrystalline phases in a coupled-mode Swift-Hohenberg model for three-component systems, where the characteristic length scales are governed by the positive-definite gradient terms. Classic two-mode approximation method and direct numerical minimization are applied to the model. In the latter approach, we apply the projection method to deal with the potentially quasiperiodic ground states. A variable cell method of optimizing the shape and size of higher-dimensional periodic cell is developed to minimize the free energy with respect to the order parameters. Based on the developed numerical methods, we rediscover decagonal and dodecagonal quasicrystalline phases, and find diverse periodic phases and complex modulated phases. Furthermore, phase diagrams are obtained in various phase spaces by comparing the free energies of different candidate structures. It does show not only the important roles of system parameters, but also the effect of optimizing computational domain. In particular, the optimization of computational cell allows us to capture the ground states and phase behavior with higher fidelity. We also make some discussions on our results and show the potential of applying our numerical methods to a larger class of mean-field free energy functionals.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stability of Two-Dimensional Soft Quasicrystals

The relative stability of two-dimensional soft quasicrystals is examined using a recently developed projection method which provides a unified numerical framework to compute the free energy of periodic crystal and quasicrystals. Accurate free energies of numerous ordered phases, including dodecagonal, decagonal and octagonal quasicrystals, are obtained for a simple model, i.e. the Lifshitz-Petrich free energy functional, of soft quasicrystals with two length-scales. The availability of the free energy allows us to construct phase diagrams of the system, demonstrating that, for the Lifshitz-Petrich model, the dodecagonal and decagonal quasicrystals can become stable phases, whereas the octagonal quasicrystal stays as a metastable phase.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Transition of defect patterns from 2D to 3D in liquid crystals

Defects arise when nematic liquid crystals are under topological constraints at the boundary. Recently the study of defects has drawn a lot of attention. In this paper, we investigate the relationship between two-dimensional defects and three-dimensional defects within nematic liquid crystals confined in a shell under the Landau-de Gennes model. We use a highly accurate spectral method to numerically solve the Landau- de Gennes model to get the detailed static structures of defects. Interestingly, the solution is radial-invariant when the thickness of the shell is sufficiently small. As the shell thickness increase, the solution undergo symmetry break to reconfigure the disclination lines. We study this three-dimensional reconfiguration of disclination lines in detail under different boundary conditions. We also discuss the topological charge of defects in two- and three-dimensional spaces within the tensor model.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Local well-posedness and small Deborah limit of a molecule-based $Q$-tensor system

In this paper, we consider a hydrodynamic $Q$-tensor system for nematic liquid crystal flow, which is derived from Doi-Onsager molecular theory by the Bingham closure. We first prove the existence and uniqueness of local strong solution. Furthermore, by taking Deborah number goes to zero and using the Hilbert expansion method, we present a rigorous derivation from the molecule-based $Q$-tensor theory to the Ericksen-Leslie theory.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Analytic Structure of the SCFT Energy Functional of Multicomponent Block Copolymers

This paper concerns the analytic structure of the self-consistent field theory (SCFT) energy functional of multicomponent block copolymer systems which contain more than two chemically distinct blocks. The SCFT has enjoyed considered success and wide usage in investigation of the complex phase behavior of block copolymers. It is well-known that the physical solutions of the SCFT equations are saddle points, however, the analytic structure of the SCFT energy functional has received little attention over the years. A recent work by Fredrickson and collaborators [see the monograph by Fredrickson, The Equilibrium Theory of Inhomogeneous Polymers, (2006), pp. 203-209] has analysed the mathematical structure of the field energy functional for polymeric systems, and clarified the index-1 saddle point nature of the problem produced by the incompressibility constraint. In this paper, our goals are to draw further attention to multicomponent block copolymers utilizing the Hubbard-Stratonovich transformation used by Fredrickson and co-workers. We first show that the saddle point character of the SCFT energy functional of multicomponent block copolymer systems may be high index, not only produced by the incompressibility constraint, but also by the Flory-Huggins interaction parameters. Our analysis will be beneficial to many theoretical studies, such as the nucleation theory of ordered phases, the mesoscopic dynamics. As an application, we then utilize the discovery to develop the gradient-based iterative schemes to solve the SCFT equations, and illustrate its performance through several numerical experiments taking ABC star triblock copolymers as an example.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

From microscopic theory to macroscopic theory: a systematic study on static modeling for liquid crystals

In this paper, we propose a systematic way of liquid crystal modeling to build connection between microscopic theory and macroscopic theory. A new Q-tensor theory based on Onsager's molecular theory which leads to liquid crystals with certain shape has been proposed. Making uniaxial assumption, we can recover the Oseen-Frank theory from the derived $Q$-tensor theory, and the Oseen-Frank model coefficients can be examined. In addition, the smectic-A phase can also be characterized by the derived macroscopic model.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Is there a third order phase transition for supercritical fluids?

We prove that according to Molecular Dynamics (MD) simulations of liquid mixtures of Lennard-Jones (L-J) particles, there is no third order phase transition in the supercritical regime beyond Andrew's critical point. This result is in open contrast with recent theoretical studies and experiments which instead suggest not only its existence but also its universality regarding the chemical nature of the fluid. We argue that our results are solid enough to go beyond the limitations of MD and the generic character of L-J models, thus suggesting a rather smooth liquid-vapor thermodynamic behavior of fluids in supercritical regime.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Numerical Methods for Quasicrystals

Quasicrystals are one kind of space-filling structures. The traditional crystalline approximant method utilizes periodic structures to approximate quasicrystals. The errors of this approach come from two parts: the numerical discretization, and the approximate error of Simultaneous Diophantine Approximation which also determines the size of the domain necessary for accurate solution. As the approximate error decreases, the computational complexity grows rapidly, and moreover, the approximate error always exits unless the computational region is the full space. In this work we focus on the development of numerical method to compute quasicrystals with high accuracy. With the help of higher-dimensional reciprocal space, a new projection method is developed to compute quasicrystals. The approach enables us to calculate quasicrystals rather than crystalline approximants. Compared with the crystalline approximant method, the projection method overcomes the restrictions of the Simultaneous Diophantine Approximation, and can also use periodic boundary conditions conveniently. Meanwhile, the proposed method efficiently reduces the computational complexity through implementing in a unit cell and using pseudospectral method. For illustrative purpose we work with the Lifshitz-Petrich model, though our present algorithm will apply to more general systems including quasicrystals. We find that the projection method can maintain the rotational symmetry accurately. More significantly, the algorithm can calculate the free energy density to high precision.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Unstable Galaxy Models

The dynamics of collisionless galaxy can be described by the Vlasov-Poisson system. By the Jean's theorem, all the spherically symmetric steady galaxy models are given by a distribution of Φ(E,L), where E is the particle energy and L the angular momentum. In a celebrated Doremus-Feix-Baumann Theorem, the galaxy model Φ(E,L) is stable if the distribution Φ is monotonically decreasing with respect to the particle energy E. On the other hand, the stability of Φ(E,L) remains largely open otherwise. Based on a recent abstract instability criterion of Guo-Lin, we constuct examples of unstable galaxy models of f(E,L) and f(E) in which f fails to be monotone in E.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Optimizing working parameters of the twin-range cutoff method in terms of accuracy and efficiency

We construct a priori error estimation for the force error of the twin-range cutoff method, which is widely used to treat the short-range non-bonded interactions in molecular simulations. Based on the error and cost estimation, we develop a work flow that can automatically determine the nearly most efficient twin-range cutoff parameters (i.e. the cutoff radii and the neighbor list updating frequency) prior to a simulation for a predetermined accuracy. Both the error estimate and the parameter tuning method are demonstrated to be effective by testing simulations of the standard Lennard-Jones 6-12 fluid in gas, liquid as well as supercritical state. We recommend the tuned twin-range cutoff method that can save precious user time and computational resources.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The small Deborah number limit of the Doi-Onsager equation to the Ericksen-Leslie equation

We present a rigorous derivation of the Ericksen-Leslie equation starting from the Doi-Onsager equation. As in the fluid dynamic limit of the Boltzmann equation, we first make the Hilbert expansion for the solution of the Doi-Onsager equation. The existence of the Hilbert expansion is connected to an open question whether the energy of the Ericksen-Leslie equation is dissipated. We show that the energy is dissipated for the Ericksen-Leslie equation derived from the Doi-Onsager equation. The most difficult step is to prove a uniform bound for the remainder in the Hilbert expansion. This question is connected to the spectral stability of the linearized Doi-Onsager operator around a critical point. By introducing two important auxiliary operators, the detailed spectral information is obtained for the linearized operator around all critical points. However, these are not enough to justify the small Deborah number limit for the inhomogeneous Doi-Onsager equation, since the elastic stress in the velocity equation is also strongly singular. For this, we need to establish a precise lower bound for a bilinear form associated with the linearized operator. In the bilinear form, the interactions between the part inside the kernel and the part outside the kernel of the linearized operator are very complicated. We find a coordinate transform and introduce a five dimensional space called the Maier-Saupe space such that the interactions between two parts can been seen explicitly by a delicate argument of completing the square. However, the lower bound is very weak for the part inside the Maier-Saupe space. In order to apply them to the error estimates, we have to analyze the structure of the singular terms and introduce a suitable energy functional.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Well-posedness of Hydrodynamics on the Moving Elastic Surface

The dynamics of a membrane is a coupled system comprising a moving elastic surface and an incompressible membrane fluid. We will consider a reduced elastic surface model, which involves the evolution equations of the moving surface, the dynamic equations of the two-dimensional fluid, and the incompressible equation, all of which operate within a curved geometry. In this paper, we prove the local existence and uniqueness of the solution to the reduced elastic surface model by reformulating the model into a new system in the isothermal coordinates. One major difficulty is that of constructing an appropriate iterative scheme such that the limit system is consistent with the original system.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Well-posedness of the Ericksen-Leslie system

In this paper, we prove the local well-posedness of the Ericksen-Leslie system, and the global well-posednss for small initial data under the physical constrain condition on the Leslie coefficients, which ensures that the energy of the system is dissipated. Instead of the Ginzburg-Landau approximation, we construct an approximate system with the dissipated energy based on a new formulation of the system.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dynamic Transitions and Pattern Formations for Cahn-Hilliard Model with Long-Range Repulsive Interactions

The main objective of this article is to study the order-disorder phase transition and pattern formation for systems with long-range repulsive interactions. The main focus is on the Cahn-Hilliard model with a nonlocal term in the corresponding energy functional, representing the long-range repulsive interaction. First, we show that as soon as the linear problem loses stability, the system always undergoes a dynamic transition to one of the three types, forming different patterns/structures. The types of transition are then dictated by a nondimensional parameter, measuring the interactions between the long-range repulsive term and the quadratic and cubic nonlinearities in the model. The derived explicit form of this parameter offers precise information for the phase diagrams. Second, we obtain a novel and explicit pattern selection mechanism associated with the competition between the long-range repulsive interaction and the short-range attractive interactions. In particular, the hexagonal pattern is unique to the long-range interaction, and is associated with a novel two-dimensional reduced transition equations on the center manifold generated by the unstable modes, consisting of (degenerate) quadratic terms and non-degenerate cubic terms. Finally, explicit information on the metastability and basin of attraction of different disordered/ordered states and patterns are derived as well.

0 citations0 reviewsNo code linkedOpen access