Researcher profile

Jianyuan Xiao

Jianyuan Xiao contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
7topics
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

5 published item(s)

preprint2022arXiv

Explicit high-order symplectic integrators of coupled Schrodinger equations for pump-probe systems

Two-beam coupling within the field of nonlinear optics, which transfers energy from one light beam to the other under certain conditions, has received considerable attention in inertial confinement fusion (ICF) and plasma optics. To evaluate the coupling dynamics precisely, we modeled this process with full-wave coupled Schrodinger equations (CSEs) and a nonlinear refractive index. We found that the CSEs constituted a Hamiltonian system and proposed an arbitrary higher-order explicit symplectic algorithm to solve the CSEs numerically. The numerical results given by the developed BEAM code showed a good agreement with those from particle-in-cell simulations, which demonstrated the validity of the model and algorithm. The model and numerical algorithm presented in this work can be extended to more nonlinear optical interactions described by coupled-wave equations.

preprint2021arXiv

Gauge invariant canonical symplectic algorithms for real-time lattice strong-field quantum electrodynamics

A class of high-order canonical symplectic structure-preserving geometric algorithms are developed for high-quality simulations of the quantized Dirac-Maxwell theory based strong-field quantum electrodynamics (SFQED) and relativistic quantum plasmas (RQP) phenomena. The Lagrangian density of an interacting bispinor-gauge fields theory is constructed in a conjugate real fields form. The canonical symplectic form and canonical equations of this field theory are obtained by the general Hamilton's principle on cotangent bundle. Based on discrete exterior calculus, the gauge field components are discreted to form a cochain complex, and the bispinor components are naturally discreted on a staggered dual lattice as combinations of differential forms. With pull-back and push-forward gauge covariant derivatives, the discrete action is gauge invariant. A well-defined discrete canonical Poisson bracket generates a semi-discrete lattice canonical field theory (LCFT), which admits the canonical symplectic form, unitary property, gauge symmetry and discrete Poincaré subgroup. The Hamiltonian splitting method, Cayley transformation and symmetric composition technique are introduced to construct a class of high-order numerical schemes. These schemes involve two degenerate fermion flavors and are locally unconditional stable, which also preserve the geometric structures. Equipped with statistically quantization-equivalent ensemble models of the Dirac vacuum and non-trivial plasma backgrounds, the schemes are expected to have excellent performance in secular simulations of relativistic quantum effects. The algorithms are verified in detail by numerical energy spectra. Real-time LCFT simulations are successfully implemented for the nonlinear Schwinger mechanism induced $e$-$e^+$ pairs creation and vacuum Kerr effect, which open a new door toward high-quality simulations in SFQED and RQP.

preprint2021arXiv

Hybrid simulation of energetic particles interacting with magnetohydrodynamics using a slow manifold algorithm and GPU acceleration

The hybrid method combining particle-in-cell and magnetohydrodynamics can be used to study the interaction between energetic particles and global plasma modes. In this paper we introduce the M3D-C1-K code, which is developed based on the M3D-C1 finite element code solving the magnetohydrodynamics equations, with a newly developed kinetic module simulating energetic particles. The particle pushing is done using a new algorithm by applying the Boris pusher to the classical Pauli particles to simulate the slow-manifold of particle orbits, with long-term accuracy and fidelity. The particle pushing can be accelerated using GPUs with a significant speedup. The moments of the particles are calculated using the $δf$ method, and are coupled into the magnetohydrodynamics simulation through pressure or current coupling schemes. Several linear simulations of magnetohydrodynamics modes driven by energetic particles have been conducted using M3D-C1-K, including fishbone, toroidal Alfvén eigenmodes and reversed shear Alfvén eigenmodes. Good agreement with previous results from other eigenvalue, kinetic and hybrid codes has been achieved.

preprint2019arXiv

Field theory and structure-preserving geometric particle-in-cell algorithm for drift wave instability and turbulence

A field theory and the associated structure-preserving geometric Particle-In-Cell (PIC) algorithm are developed to study low frequency electrostatic perturbations with fully kinetic ions and adiabatic electrons in magnetized plasmas. The algorithm is constructed by geometrically discretizing the field theory using discrete exterior calculus, high-order Whitney interpolation forms, and non-canonical Hamiltonian splitting method. The discretization preserves the non-canonical symplectic structure of the particle-field system, as well as the electromagnetic gauge symmetry. As a result, the algorithm is charge-conserving and possesses long-term conservation properties. Because drift wave turbulence and anomalous transport intrinsically involve multi time-scales, simulation studies using fully kinetic particle demand algorithms with long-term accuracy and fidelity. The structure-preserving geometric PIC algorithm developed adequately servers this purpose. The algorithm has been implemented in the \textsl{SymPIC} code, tested and benchmarked using the examples of ion Bernstein waves and drift waves. We apply the algorithm to study the Ion Temperature Gradient (ITG) instability and turbulence in a 2D slab geometry. Simulation results show that at the early stage of the turbulence, the energy diffusion is between the Bohm scaling and gyro-Bohm scaling. At later time, the observed diffusion is closer to the gyro-Bohm scaling, and density blobs generated by the rupture of unstable modes are the prominent structures of the fully developed ITG turbulence.

preprint2019arXiv

PT-symmetry entails pseudo-Hermiticity regardless of diagonalizability

We prove that in finite dimensions, a Parity-Time (PT)-symmetric Hamiltonian is necessarily pseudo-Hermitian regardless of whether it is diagonalizable or not. This result is different from Mostafazadeh's, which requires the Hamiltonian to be diagonalizable. PT-symmetry breaking often occurs at exceptional points where the Hamiltonian is not diagonalizable. Our result implies that PT-symmetry breaking is equivalent to the onset of instabilities of pseudo-Hermitian systems, which was systematically studied by Krein et al. in 1950s. In particular, we show that the mechanism of PT-symmetry breaking is the resonance between eigenmodes with different Krein signatures.