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

Zhenning Cai

Zhenning Cai contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.flu-dynmath.NAphysics.comp-phquant-phmath-phmath.MPNumerical Analysiscond-mat.stat-mechhep-lat

Trust snapshot

Quick read

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

11 published item(s)

preprint2022arXiv

Asymptotic analysis of diabatic surface hopping algorithm in the adiabatic and non-adiabatic limits

Surface hopping algorithms, as an important class of quantum dynamics simulation algorithms for non-adiabatic dynamics, are typically performed in the adiabatic representation, which can break down in the presence of ill-defined adiabatic potential energy surfaces (PESs) and adiabatic coupling term. Another issue of surface hopping algorithms is the difficulty in capturing the correct scaling of the transition rate in the Marcus (weak-coupling/non-adiabatic) regime. Though the first issue can be circumvented by exploiting the diabatic representation, diabatic surface hopping algorithms usually lack justification on the theoretical level. We consider the diabatic surface hopping algorithm proposed in [Fang, Lu. Multiscale Model. Simul. 16:4, 1603-1622, 2018] and provide the asymptotic analysis of the transition rate in the Marcus regime that justifies the correct scaling for the spin-boson model. We propose two conditions that guarantee the correctness for general potentials. In the opposite (strong-coupling/adiabatic) regime, we derive the asymptotic behavior of the algorithm that interestingly matches a type of mean-field description. The techniques used here may shed light on the analysis for other diabatic-based algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fast Algorithms of Bath Calculations in Simulations of Quantum System-Bath Dynamics

We present fast algorithms for the summation of Dyson series and the inchworm Monte Carlo method for quantum systems that are coupled with harmonic baths. The algorithms are based on evolving the integro-differential equations where the most expensive part comes from the computation of bath influence functionals. To accelerate the computation, we design fast algorithms based on reusing the bath influence functionals computed in the previous time steps to reduce the number of calculations. It is proven that the proposed fast algorithms reduce the number of such calculations by a factor of $O(N)$, where $N$ is the total number of time steps. Numerical experiments are carried out to show the efficiency of the method and to verify the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Differential Equation Based Path Integral for System-Bath Dynamics

We propose the differential equation based path integral (DEBPI) method to simulate the real-time evolution of open quantum systems. In this method, a system of partial differential equations is derived based on the continuation of a classical numerical method called iterative quasi-adiabatic propagator path integral (i-QuAPI). While the resulting system has infinite equations, we introduce a reasonable closure to obtain a series of finite systems. New numerical schemes can be derived by discretizing these differential equations. It is numerically verified that in certain cases, by selecting appropriate systems and applying suitable numerical schemes, the memory cost required in the i-QuAPI method can be significantly reduced.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Regularization of Complex Langevin Method

The complex Langevin method, a numerical method used to compute the ensemble average with a complex partition function, often suffers from runaway instability. We study the regularization of the complex Langevin method via augmenting the action with a stabilization term. Since the regularization introduces biases to the numerical result, two approaches, named 2R and 3R methods, are introduced to recover the unbiased result. The 2R method supplements the regularization with regression to estimate the unregularized ensemble average, and the 3R method reduces the computational cost by coupling the regularization with a reweighting strategy before regression. Both methods can be generalized to the SU(n) theory and are assessed from several perspectives. Several numerical experiments in the lattice field theory are carried out to show the effectiveness of our approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Flows Between Parallel Plates: Analytical Solutions of Regularized 13-Moment Equations for Inverse-Power-Law Models

We study the structure of stationary channel flows predicted by the regularized 13-moment equations. Compared with the previous work [P. Taheri et al., Phys. Fluids, 21 (2009), 017102], we focus on gases whose molecules satisfy the general inverse power law. The analytical solutions are obtained for the semi-linear equations, and the structures of Couette, Fourier, and Poiseuille flows are solved by coupling the general solutions with newly derived boundary conditions. The results show excellent agreement with the reference solution in the slip-flow regime. Our results also show that the R13 equations derived from inverse power law models can have better accuracy than the R13 equations of Maxwell molecules with altered viscosity.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Hyperbolic Model Reduction for Kinetic Equations

We make a brief historical review to the moment model reduction to the kinetic equations, particularly the Grad's moment method for Boltzmann equation. The focus is on the hyperbolicity of the reduced model, which is essential to the existence of its classical solution as a Cauchy problem. The theory of the framework we developed in last years is then introduced, which may preserve the hyperbolic nature of the kinetic equations with high universality. Some lastest progress on the comparison between models with/without hyperbolicity is presented to validate the hyperbolic moment models for rarefied gases.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical Analysis of a High-Order Scheme for Nonlinear Fractional Differential Equations with Uniform Accuracy

We introduce a high-order numerical scheme for fractional ordinary differential equations with the Caputo derivative. The method is developed by dividing the domain into a number of subintervals, and applying the quadratic interpolation on each subinterval. The method is shown to be unconditionally stable, and for general nonlinear equations, the uniform sharp numerical order $3-ν$ can be rigorously proven for sufficiently smooth solutions at all time steps. The proof provides a general guide for proving the sharp order for higher-order schemes in the nonlinear case. Some numerical examples are given to validate our theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

How does Gauge Cooling Stabilize Complex Langevin?

We study the mechanism of the gauge cooling technique to stabilize the complex Langevin method in the one-dimensional periodic setting. In this case, we find the exact solutions for the gauge transform which minimizes the Frobenius norm of link variables. Thereby, we derive the underlying stochastic differential equations by continuing the numerical method with gauge cooling, and thus provide a number of insights on the effects of gauge cooling. A specific case study is carried out for the Polyakov loop model in $SU(2)$ theory, in which we show that the gauge cooling may help form a localized distribution to guarantee there is no excursion too far away from the real axis.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

On the Holway-Weiss Debate: Convergence of the Grad-Moment-Expansion in Kinetic Gas Theory

Moment expansions are used as model reduction technique in kinetic gas theory to approximate the Boltzmann equation. Rarefied gas models based on so-called moment equations became increasingly popular recently. However, in a seminal paper by Holway [Phys. Fluids 7/6, (1965)] a fundamental restriction on the existence of the expansion was used to explain sub-shock behavior of shock profile solutions obtained by moment equations. Later, Weiss [Phys. Fluids 8/6, (1996)] argued that this restriction does not exist. We will revisit and discuss their findings and explain that both arguments have a correct and incorrect part. While a general convergence restriction for moment expansions does exist, it cannot be attributed to sub-shock solutions. We will also discuss the implications of the restriction and give some numerical evidence for our considerations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Regularized 13-Moment Equations for Inverse Power Law Models

We propose a systematic methodology to derive the regularized thirteen-moment equations in the rarefied gas dynamics for a general class of linearized collision models. Detailed expressions of the moment equations are written down for all inverse power law models as well as the hard-sphere model. By linear analysis, we show that the equations are stable near the equilibrium. The models are tested for shock structure problems to show its capability to capture the correct flow structure in strong nonequilibirum.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Numerical Regularized Moment Method of Arbitrary Order for Boltzmann-BGK Equation

We introduce a numerical method for solving Grad's moment equations or regularized moment equations for arbitrary order of moments. In our algorithm, we do not need explicitly the moment equations. As an instead, we directly start from the Boltzmann equation and perform Grad's moment method \cite{Grad} and the regularization technique \cite{Struchtrup2003} numerically. We define a conservative projection operator and propose a fast implementation which makes it convenient to add up two distributions and provides more efficient flux calculations compared with the classic method using explicit expressions of flux functions. For the collision term, the BGK model is adopted so that the production step can be done trivially based on the Hermite expansion. Extensive numerical examples for one- and two-dimensional problems are presented. Convergence in moments can be validated by the numerical results for different number of moments.

0 citations0 reviewsNo code linkedOpen access