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

Jianyuan Xiao

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

ResearcherUnclaimed source record
physics.plasm-phphysics.comp-phquant-phhep-latmath-phmath.MPphysics.acc-phphysics.optics

Catalog footprint

What is connected

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

13 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Explicit high-order non-canonical symplectic particle-in-cell algorithms for Vlasov-Maxwell systems

Explicit high-order non-canonical symplectic particle-in-cell algorithms for classical particle-field systems governed by the Vlasov-Maxwell equations are developed. The algorithm conserves a discrete non-canonical symplectic structure derived from the Lagrangian of the particle-field system, which is naturally discrete in particles. The electromagnetic field is spatially-discretized using the method of discrete exterior calculus with high-order interpolating differential forms for a cubic grid. The resulting time-domain Lagrangian assumes a non-canonical symplectic structure. It is also gauge invariant and conserves charge. The system is then solved using a splitting method discovered by He et al., which produces five exactly-soluable sub-systems, and high-order structure- preserving algorithms follow by combinations. The explicit, high-order, and conservative nature of the algorithms is especially suitable for long-term simulations of particle-field systems with extremely large number of degrees of freedom on massively parallel supercomputers. The algorithms have been tested and verified by the two physics problems, i.e., the nonlinear Landau damping and the electron Bernstein wave.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Explicit high-order noncanonical symplectic algorithms for ideal two-fluid systems

An explicit high-order noncanonical symplectic algorithm for ideal two-fluid systems is developed. The fluid is discretized as particles in the Lagrangian description, while the electromagnetic fields and internal energy are treated as discrete differential form fields on a fixed mesh. With the assistance of Whitney interpolating forms, this scheme preserves the gauge symmetry of the electromagnetic field, and the pressure field is naturally derived from the discrete internal energy. The whole system is solved using the Hamiltonian splitting method discovered by He et al., which was been successfully adopted in constructing symplectic particle-in-cell schemes. Because of its structure preserving and explicit nature, this algorithm is especially suitable for large-scale simulations for physics problems that are multi-scale and require long-term fidelity and accuracy. The algorithm is verified via two tests: studies of the dispersion relation of waves in a two-fluid plasma system and the oscillating two-stream instability.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Explicit symplectic algorithms based on generating functions for charged particle dynamics

Dynamics of a charged particle in the canonical coordinates is a Hamiltonian system, and the well-known symplectic algorithm has been regarded as the de facto method for numerical integration of Hamiltonian systems due to its long-term accuracy and fidelity. For long-term simulations with high efficiency, explicit symplectic algorithms are desirable. However, it is widely accepted that explicit symplectic algorithms are only available for sum-separable Hamiltonians, and that this restriction severely limits the application of explicit symplectic algorithms to charged particle dynamics. To overcome this difficulty, we combine the familiar sum-split method and a generating function method to construct second and third order explicit symplectic algorithms for dynamics of charged particle. The generating function method is designed to generate explicit symplectic algorithms for product-separable Hamiltonian with form of $H(\mathbf{p},\mathbf{q})=\mathbf{p}_{i}f(\mathbf{q})$ or $H(\mathbf{p},\mathbf{q})=\mathbf{q}_{i}f(\mathbf{p})$. Applied to the simulations of charged particle dynamics, the explicit symplectic algorithms based on generating functions demonstrate superiorities in conservation and efficiency.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the structure of the two-stream instability -- complex G-Hamiltonian structure and Krein collisions between positive- and negative-action modes

The two-stream instability is probably the most important elementary example of collective instabilities in plasma physics and beam-plasma systems. For a warm plasma with two charged particle species based on a 1D warm-fluid model, the instability diagram of the two-stream instability exhibits an interesting band structure that has not been explained. We show that the band structure for this instability is the consequence of the Hamiltonian nature of the warm two-fluid system. Interestingly, the Hamiltonian nature manifests as a complex G-Hamiltonian structure in wave-number space, which directly determines the instability diagram. Specifically, it is shown that the boundaries between the stable and unstable regions are locations for Krein collisions between eigenmodes with different Krein signatures. In terms of physics, this rigorously implies that the system is destabilized when a positive-action mode resonates with a negative-action mode, and that this is the only mechanism by which the system can be destabilized. It is anticipated that this physical mechanism of destabilization is valid for other collective instabilities in conservative systems in plasma physics, accelerator physics, and fluid dynamics systems, which admit infinite-dimensional Hamiltonian structures.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Canonical symplectic particle-in-cell method for long-term large-scale simulations of the Vlasov-Maxwell system

Particle-in-Cell (PIC) simulation is the most important numerical tool in plasma physics. However, its long-term accuracy has not been established. To overcome this difficulty, we developed a canonical symplectic PIC method for the Vlasov-Maxwell system by discretizing its canonical Poisson bracket. A fast local algorithm to solve the symplectic implicit time advance is discovered without root searching or global matrix inversion, enabling applications of the proposed method to very large-scale plasma simulations with many, e.g., $10^{9}$, degrees of freedom. The long-term accuracy and fidelity of the algorithm enables us to numerically confirm Mouhot and Villani's theory and conjecture on nonlinear Landau damping over several orders of magnitude using the PIC method, and to calculate the nonlinear evolution of the reflectivity during the mode conversion process from extraordinary waves to Bernstein waves.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Comment on "Hamiltonian splitting for the Vlasov-Maxwell equations"

The paper [1] by Crouseilles, Einkemmer, and Faou used an incorrect Poisson bracket for the Vlasov-Maxwell equations. If the correct Poisson bracket is used, the solution of one of the subsystems cannot be computed exactly in general. As a result, one cannot construct a symplectic scheme for the Vlasov-Maxwell equations using the splitting Hamiltonian method proposed in Ref [1].

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hamiltonian integration methods for Vlasov-Maxwell equations

Hamiltonian integration methods for the Vlasov-Maxwell equations are developed by a Hamiltonian splitting technique. The Hamiltonian functional is split into five parts, i.e., the electrical energy, the magnetic energy, and the kinetic energy in three Cartesian components. Each of the subsystems is a Hamiltonian system with respect to the Morrison-Marsden-Weinstein Poisson bracket and can be solved exactly. Compositions of the exact solutions yield Poisson structure preserving, or Hamiltonian, integration methods for the Vlasov-Maxwell equations, which have superior long-term fidelity and accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Variational Symplectic Particle-in-cell Simulation of Nonlinear Mode Conversion from Extraordinary waves to Bernstein Waves

In this paper, the nonlinear mode conversion of extraordinary waves in nonuniform magnetized plasmas is studied using the variational symplectic particle-in-cell simulation. The accuracy of the nonlinear simulation is guaranteed by the long-term accuracy and conservativeness of the symplectic algorithm. The spectra of the electromagnetic wave, the evolution of the wave reflectivity, the energy deposition profile, and the parameter-dependent properties of radio-frequency waves during the nonlinear mode conversion are investigated. It is illustrated that nonlinear effects significantly modify the physics of the radio-frequency injection in magnetized plasmas. The evolutions of the radio-frequency wave reflectivity and the energy deposition are observed, as well as the self-interaction of the Bernstein waves and mode excitations. Even for waves with small magnitude, nonlinear effects can also become important after continuous wave injections, which are common in the realistic radio-frequency wave heating and current drive experiments.

0 citations0 reviewsNo code linkedOpen access