Source author record

Hong Qin

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

51works
23topics
4close collaborators

Actions

Connect this record

Log in to claim

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

51 published item(s)

preprint2026arXiv

Hidden free energy released by explicit parity-time-symmetry breaking

It is shown that the familiar two-stream instability is the result of spontaneous parity-time (PT)-symmetry breaking in a conservative system, and more importantly, explicit PT-symmetry breaking by viscosity can destabilize the system in certain parameter regimes that are stable when viscosity vanishes. This reveals that complex systems may possess hidden free energies protected by PT-symmetry and viscosity, albeit dissipative, can expose the systems to these freed energies by breaking PT-symmetry explicitly. Such a process is accompanied by instability and total variation growth.

preprint2022arXiv

A gauge-compatible Hamiltonian splitting algorithm for particle-in-cell simulations using finite element exterior calculus

A particle-in-cell algorithm is derived with a canonical Poisson structure in the formalism of finite element exterior calculus. The resulting method belongs to the class of gauge-compatible splitting algorithms, which exactly preserve gauge symmetries and their associated conservation laws via the momentum map. We numerically demonstrate this time invariance of the momentum map and its usefulness in establishing precise initial conditions with a desired initial electric field and fixed background charge. The restriction of this canonical, finite element Poisson structure to the 1X2P phase space is also considered and simulated numerically.

preprint2022arXiv

Cointegration of SARS-CoV-2 Transmission with Weather Conditions and Mobility during the First Year of the COVID-19 Pandemic in the United States

Correlation between weather and the transmission of SARS-CoV-2 may suggest its seasonality. Cointegration analysis can avoid spurious correlation among time series data. We examined the cointegration of virus transmission with daily temperature, dewpoint, and confounding factors of mobility measurements during the first year of the pandemic in the United States. We examined the cointegration of the effective reproductive rate, Rt, of the virus with the dewpoint at two meters, the temperature at two meters, Apple driving mobility, and Google workplace mobility measurements. We found that dewpoint and Apple driving mobility are the best factors to cointegrate with Rt, although temperature and Google workplace mobility also cointegrate with Rt at substantial levels. We found that the optimal lag is two days for cointegration between Rt and weather variables, and three days for Rt and mobility. We observed clusters of states that share similar cointegration results of Rt, weather, and mobility, suggesting regional patterns. Our results support the correlation of weather with the spread of SARS-CoV-2 and its potential seasonality.

preprint2022arXiv

Collective Effects and Intense Beam-Plasma Interactions in Ion-Beam-Driven High Energy Density Matter and Inertial Fusion Energy

For the successful generation of ion-beam-driven high energy density matter and heavy ion fusion energy, intense ion beams must be transported and focused onto a target with small spot size. One of the successful approaches to achieve this goal is to accelerate and transport intense ion charge bunches in an accelerator and then focus the charge bunches ballistically in a section of the accelerator that contains a neutralizing background plasma. This requires the ability to control space-charge effects during un-neutralized (non-neutral) beam transport in the accelerator and transport sections, and the ability to effectively neutralize the space charge and current by propagating the beam through background plasma. As the beam intensity and energy are increased in future heavy ion fusion (HIF) drivers and Fast Ignition (FI) approaches, it is expected that nonlinear processes and collective effects will become much more pronounced than in previous experiments. Making use of 3D electromagnetic particle-in-cell simulation (PIC) codes (BEST, WARP-X, and LTP-PIC, etc.), the theory and modelling studies will be validated by comparing with experimental data on the 100kV Princeton Advanced Test Stand, and future experiments at the FAIR facility. The theoretical predictions that are developed will be scaled to the beam and plasma parameters relevant to heavy ion fusion drivers and Fast Ignition scenarios. Therefore, the theoretical results will also contribute significantly toward the long-term goal of fusion energy production by ion-beam-driven inertial confinement fusion.

preprint2022arXiv

Quantum Instability

The physics of many closed, conservative systems can be described by both classical and quantum theories. The dynamics according to classical theory is symplectic and admits linear instabilities which would initially seem at odds with a unitary quantum description. Using the example of three-wave interactions, we describe how a time-independent, finite-dimensional quantum system, which is Hermitian with all real eigenvalues, can give rise to a linear instability corresponding to that in the classical system. We show that the instability is realized in the quantum theory as a cascade of the wave function in the space of occupation number states, and an unstable quantum system has a richer spectrum and a much longer recurrence time than a stable quantum system. The conditions for quantum instability are described.

preprint2022arXiv

Shifting Trends of COVID-19 Tweet Sentiment with Respect to Voting Preferences in the 2020 Election Year of the United States

COVID-19 related policies were extensively politicized during the 2020 election year of the United States, resulting in polarizing viewpoints. Twitter users were particularly engaged during the 2020 election year. Here we investigated whether COVID-19 related tweets were associated with the overall election results at the state level during the period leading up to the election day. We observed weak correlations between the average sentiment of COVID-19 related tweets and popular votes in two-week intervals, and the trends gradually become opposite. We then compared the average sentiments of COVID-19 related tweets between states called in favor of Republican (red states) or Democratic parties (blue states). We found that at the beginning of lockdowns sentiments in the blue states were much more positive than those in the red states. However, sentiments in the red states gradually become more positive during the summer of 2020 and persisted until the election day.

preprint2022arXiv

The dispersion and propagation of topological Langmuir-cyclotron waves in cold magnetized plasmas

Topological Langmuir-Cyclotron Wave (TLCW) is a recently identified topological surface excitation in magnetized plasmas. We show that TLCW originates from the topological phase transition at the Langmuir wave-cyclotron wave resonance. By isofrequency surface analysis and 2D and 3D time-dependent simulations, we demonstrate that the TLCW can propagate robustly along complex phase transition interfaces in a unidirectional manner and without scattering. Because of these desirable features, the TLCW could be explored as an effective mechanism to drive current and flow in magnetized plasmas. The analysis also establishes a close connection between the newly instituted topological phase classification of plasmas and the classical CMA diagram of plasma waves.

preprint2022arXiv

Weakly Supervised Visual-Auditory Fixation Prediction with Multigranularity Perception

Thanks to the rapid advances in deep learning techniques and the wide availability of large-scale training sets, the performance of video saliency detection models has been improving steadily and significantly. However, deep learning-based visualaudio fixation prediction is still in its infancy. At present, only a few visual-audio sequences have been furnished, with real fixations being recorded in real visual-audio environments. Hence, it would be neither efficient nor necessary to recollect real fixations under the same visual-audio circumstances. To address this problem, this paper promotes a novel approach in a weakly supervised manner to alleviate the demand of large-scale training sets for visual-audio model training. By using only the video category tags, we propose the selective class activation mapping (SCAM) and its upgrade (SCAM+). In the spatial-temporal-audio circumstance, the former follows a coarse-to-fine strategy to select the most discriminative regions, and these regions are usually capable of exhibiting high consistency with the real human-eye fixations. The latter equips the SCAM with an additional multi-granularity perception mechanism, making the whole process more consistent with that of the real human visual system. Moreover, we distill knowledge from these regions to obtain complete new spatial-temporal-audio (STA) fixation prediction (FP) networks, enabling broad applications in cases where video tags are not available. Without resorting to any real human-eye fixation, the performances of these STA FP networks are comparable to those of fully supervised networks. The code and results are publicly available at https://github.com/guotaowang/STANet.

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.

preprint2020arXiv

$\mathbb{1}$-Loop Theory

A new formalism for lattice gauge theory is developed that preserves Poincaré symmetry in a discrete universe. We define the $\mathbb{1}$-loop, a generalization of the Wilson loop that reformulates classical differential equations of motion as identity-valued multiplicative loops of Lie group elements of the form ${[g_1\cdots g_n]=\mathbb{1}}$. A lattice Poincaré gauge theory of gravity is thus derived that employs a novel matter field construction and recovers Einstein's vacuum equations in the appropriate limit.

preprint2020arXiv

A Deeper Look at Salient Object Detection: Bi-stream Network with a Small Training Dataset

Compared with the conventional hand-crafted approaches, the deep learning based methods have achieved tremendous performance improvements by training exquisitely crafted fancy networks over large-scale training sets. However, do we really need large-scale training set for salient object detection (SOD)? In this paper, we provide a deeper insight into the interrelationship between the SOD performances and the training sets. To alleviate the conventional demands for large-scale training data, we provide a feasible way to construct a novel small-scale training set, which only contains 4K images. Moreover, we propose a novel bi-stream network to take full advantage of our proposed small training set, which is consisted of two feature backbones with different structures, achieving complementary semantical saliency fusion via the proposed gate control unit. To our best knowledge, this is the first attempt to use a small-scale training set to outperform state-of-the-art models which are trained on large-scale training sets; nevertheless, our method can still achieve the leading state-of-the-art performance on five benchmark datasets.

preprint2020arXiv

A Plug-and-play Scheme to Adapt Image Saliency Deep Model for Video Data

With the rapid development of deep learning techniques, image saliency deep models trained solely by spatial information have occasionally achieved detection performance for video data comparable to that of the models trained by both spatial and temporal information. However, due to the lesser consideration of temporal information, the image saliency deep models may become fragile in the video sequences dominated by temporal information. Thus, the most recent video saliency detection approaches have adopted the network architecture starting with a spatial deep model that is followed by an elaborately designed temporal deep model. However, such methods easily encounter the performance bottleneck arising from the single stream learning methodology, so the overall detection performance is largely determined by the spatial deep model. In sharp contrast to the current mainstream methods, this paper proposes a novel plug-and-play scheme to weakly retrain a pretrained image saliency deep model for video data by using the newly sensed and coded temporal information. Thus, the retrained image saliency deep model will be able to maintain temporal saliency awareness, achieving much improved detection performance. Moreover, our method is simple yet effective for adapting any off-the-shelf pre-trained image saliency deep model to obtain high-quality video saliency detection. Additionally, both the data and source code of our method are publicly available.

preprint2020arXiv

Deep Patch-based Human Segmentation

3D human segmentation has seen noticeable progress in re-cent years. It, however, still remains a challenge to date. In this paper, weintroduce a deep patch-based method for 3D human segmentation. Wefirst extract a local surface patch for each vertex and then parameterizeit into a 2D grid (or image). We then embed identified shape descriptorsinto the 2D grids which are further fed into the powerful 2D Convolu-tional Neural Network for regressing corresponding semantic labels (e.g.,head, torso). Experiments demonstrate that our method is effective inhuman segmentation, and achieves state-of-the-art accuracy.

preprint2020arXiv

Depth Quality Aware Salient Object Detection

The existing fusion based RGB-D salient object detection methods usually adopt the bi-stream structure to strike the fusion trade-off between RGB and depth (D). The D quality usually varies from scene to scene, while the SOTA bi-stream approaches are depth quality unaware, which easily result in substantial difficulties in achieving complementary fusion status between RGB and D, leading to poor fusion results in facing of low-quality D. Thus, this paper attempts to integrate a novel depth quality aware subnet into the classic bi-stream structure, aiming to assess the depth quality before conducting the selective RGB-D fusion. Compared with the SOTA bi-stream methods, the major highlight of our method is its ability to lessen the importance of those low-quality, no-contribution, or even negative-contribution D regions during the RGB-D fusion, achieving a much improved complementary status between RGB and D.

preprint2020arXiv

Energy Transformations in a Relativistic Engine

In a previous paper we have shown that Newton's third law cannot strictly hold in a distributed system of which the different parts are at a finite distance from each other. This is due to the finite speed of signal propagation which cannot exceed the speed of light at vacuum, which in turn means that when summing the total force in the system the force does not add up to zero. This was demonstrated in a specific example of two current loops with time dependent currents, the above analysis led to suggestion of a relativistic engine. Since the system is effected by a total force for a finite period of time this means that the system acquires mechanical momentum and energy, the question then arises how can we accommodate the law of momentum and energy conservation. The subject of momentum conversation was discussed in a previous paper, while preliminary results regarding energy conservation were discussed in previous papers. Here we give a complete analysis of the exchange of energy between the mechanical part of the relativistic engine and the field part, the energy radiated from the relativistic engine is also discussed. We show that the relativistic engine effect on the energy is 4th-order in 1/c and no lower order relativistic engine effect on the energy exist.

preprint2020arXiv

GmFace: A Mathematical Model for Face Image Representation Using Multi-Gaussian

Establishing mathematical models is a ubiquitous and effective method to understand the objective world. Due to complex physiological structures and dynamic behaviors, mathematical representation of the human face is an especially challenging task. A mathematical model for face image representation called GmFace is proposed in the form of a multi-Gaussian function in this paper. The model utilizes the advantages of two-dimensional Gaussian function which provides a symmetric bell surface with a shape that can be controlled by parameters. The GmNet is then designed using Gaussian functions as neurons, with parameters that correspond to each of the parameters of GmFace in order to transform the problem of GmFace parameter solving into a network optimization problem of GmNet. The face modeling process can be described by the following steps: (1) GmNet initialization; (2) feeding GmNet with face image(s); (3) training GmNet until convergence; (4) drawing out the parameters of GmNet (as the same as GmFace); (5) recording the face model GmFace. Furthermore, using GmFace, several face image transformation operations can be realized mathematically through simple parameter computation.

preprint2020arXiv

Knowing Depth Quality In Advance: A Depth Quality Assessment Method For RGB-D Salient Object Detection

Previous RGB-D salient object detection (SOD) methods have widely adopted deep learning tools to automatically strike a trade-off between RGB and D (depth), whose key rationale is to take full advantage of their complementary nature, aiming for a much-improved SOD performance than that of using either of them solely. However, such fully automatic fusions may not always be helpful for the SOD task because the D quality itself usually varies from scene to scene. It may easily lead to a suboptimal fusion result if the D quality is not considered beforehand. Moreover, as an objective factor, the D quality has long been overlooked by previous work. As a result, it is becoming a clear performance bottleneck. Thus, we propose a simple yet effective scheme to measure D quality in advance, the key idea of which is to devise a series of features in accordance with the common attributes of high-quality D regions. To be more concrete, we conduct D quality assessments for each image region, following a multi-scale methodology that includes low-level edge consistency, mid-level regional uncertainty and high-level model variance. All these components will be computed independently and then be assembled with RGB and D features, applied as implicit indicators, to guide the selective fusion. Compared with the state-of-the-art fusion schemes, our method can achieve a more reasonable fusion status between RGB and D. Specifically, the proposed D quality measurement method achieves steady performance improvements for almost 2.0\% in general.

preprint2020arXiv

Recursive Multi-model Complementary Deep Fusion forRobust Salient Object Detection via Parallel Sub Networks

Fully convolutional networks have shown outstanding performance in the salient object detection (SOD) field. The state-of-the-art (SOTA) methods have a tendency to become deeper and more complex, which easily homogenize their learned deep features, resulting in a clear performance bottleneck. In sharp contrast to the conventional ``deeper'' schemes, this paper proposes a ``wider'' network architecture which consists of parallel sub networks with totally different network architectures. In this way, those deep features obtained via these two sub networks will exhibit large diversity, which will have large potential to be able to complement with each other. However, a large diversity may easily lead to the feature conflictions, thus we use the dense short-connections to enable a recursively interaction between the parallel sub networks, pursuing an optimal complementary status between multi-model deep features. Finally, all these complementary multi-model deep features will be selectively fused to make high-performance salient object detections. Extensive experiments on several famous benchmarks clearly demonstrate the superior performance, good generalization, and powerful learning ability of the proposed wider framework.

preprint2020arXiv

Rethinking of the Image Salient Object Detection: Object-level Semantic Saliency Re-ranking First, Pixel-wise Saliency Refinement Latter

The real human attention is an interactive activity between our visual system and our brain, using both low-level visual stimulus and high-level semantic information. Previous image salient object detection (SOD) works conduct their saliency predictions in a multi-task manner, i.e., performing pixel-wise saliency regression and segmentation-like saliency refinement at the same time, which degenerates their feature backbones in revealing semantic information. However, given an image, we tend to pay more attention to those regions which are semantically salient even in the case that these regions are perceptually not the most salient ones at first glance. In this paper, we divide the SOD problem into two sequential tasks: 1) we propose a lightweight, weakly supervised deep network to coarsely locate those semantically salient regions first; 2) then, as a post-processing procedure, we selectively fuse multiple off-the-shelf deep models on these semantically salient regions as the pixel-wise saliency refinement. In sharp contrast to the state-of-the-art (SOTA) methods that focus on learning pixel-wise saliency in "single image" using perceptual clues mainly, our method has investigated the "object-level semantic ranks between multiple images", of which the methodology is more consistent with the real human attention mechanism. Our method is simple yet effective, which is the first attempt to consider the salient object detection mainly as an object-level semantic re-ranking problem.

preprint2020arXiv

The geometric theory of charge conservation in particle-in-cell simulations

In recent years, several gauge-symmetric particle-in-cell (PIC) methods have been developed whose simulations of particles and electromagnetic fields exactly conserve charge. While it is rightly observed that these methods' gauge symmetry gives rise to their charge conservation, this causal relationship has generally been asserted via ad hoc derivations of the associated conservation laws. In this work, we develop a comprehensive theoretical grounding for charge conservation in gauge-symmetric Lagrangian and Hamiltonian PIC algorithms. For Lagrangian variational PIC methods, we apply Noether's second theorem to demonstrate that gauge symmetry gives rise to a local charge conservation law as an off-shell identity. For Hamiltonian splitting methods, we show that the momentum map establishes their charge conservation laws. We define a new class of algorithms -- gauge-compatible splitting methods -- that exactly preserve the momentum map associated with a Hamiltonian system's gauge symmetry -- even after time discretization. This class of algorithms affords splitting schemes a decided advantage over alternative Hamiltonian integrators. We apply this general technique to design a novel, explicit, symplectic, gauge-compatible splitting PIC method, whose momentum map yields an exact local charge conservation law. Our study clarifies the appropriate initial conditions for such schemes and examines their symplectic reduction.

preprint2020arXiv

The Physics of Spontaneous Parity-Time Symmetry Breaking in the Kelvin-Helmholtz Instability

We show that the dynamics, in particular the Kelvin-Helmholtz (KH) instability, of an inviscid fluid with velocity shear admits Parity-Time (PT) symmetry, which provides a physical explanation to the well-known observation that the spectrum of the perturbation eigenmodes of the system is symmetric with respect to the real axis. It is found that the KH instability is triggered when and only when the PT symmetry is spontaneously broken. The analysis of PT symmetry also reveals that the relative phase between parallel velocity and pressure perturbations needs to be locked at $π/2$ when the instability is suppressed.

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

Linear Beam Stability in Periodic Focusing Systems:Krein Signature and Band Structure

The general question of how a beam becomes unstable has been one of the fundamental research topics among beam and accelerator physicists for several decades. In this study, we revisited the general problem of linear beam stability in periodic focusing systems by applying the concepts of Krein signature and band structure. We numerically calculated the eigenvalues and other associated characteristics of one-period maps, and discussed the stability properties of single-particle motions with skew quadrupoles and envelope perturbations in high-intensity beams on an equal footing. In particular, an application of the Krein theory to envelope instability analysis was newly attempted in this study. The appearance of instabilities is interpreted as the result of the collision between eigenmodes of opposite Krein signatures and the formation of a band gap.

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.

preprint2016arXiv

Approach to Chandrasekhar-Kendall-Woltjer State in a Chiral Plasma

We study the time evolution of the magnetic field in a plasma with a chiral magnetic current. The Vector Spherical Harmonic functions (VSH) are used to expand all fields. We define a measure for the Chandrasekhar-Kendall-Woltjer (CKW) state, which has a simple form in VSH expansion. We propose the conditions for a general class of initial momentum spectra that will evolve into the CKW state. For this class of initial conditions, to approach the CKW state, (i) a non-vanishing chiral magnetic conductivity is necessary, and (ii) the time integration of the product of the electric resistivity and chiral magnetic conductivity must grow faster than the time integration of the resistivity. We give a few examples to test these conditions numerically which work very well.

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.

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.

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.

preprint2016arXiv

Formation of current singularity in a topologically constrained plasma

Recently a variational integrator for ideal magnetohydrodynamics in Lagrangian labeling has been developed. Its built-in frozen-in equation makes it optimal for studying current sheet formation. We use this scheme to study the Hahm-Kulsrud-Taylor problem, which considers the response of a 2D plasma magnetized by a sheared field under sinusoidal boundary forcing. We obtain an equilibrium solution that preserves the magnetic topology of the initial field exactly, with a fluid mapping that is non-differentiable. Unlike previous studies that examine the current density output, we identify a singular current sheet from the fluid mapping. These results are benchmarked with a constrained Grad-Shafranov solver. The same signature of current singularity can be found in other cases with more complex magnetic topologies.

preprint2016arXiv

Generalized Kapchinskij-Vladimirskij Distribution and Beam Matrix for Phase-Space Manipulations of High-Intensity Beams

In an uncoupled linear lattice system, the Kapchinskij-Vladimirskij (KV) distribution, formulated on the basis of the single-particle Courant-Snyder (CS) invariants, has served as a fundamental theoretical basis for the analyses of the equilibrium, stability, and transport properties of high-intensity beams for the past several decades. Recent applications of high-intensity beams, however, require beam phase-space manipulations by intentionally introducing strong coupling. In this Letter, we report the full generalization of the KV model by including all of the linear (both external and space-charge) coupling forces, beam energy variations, and arbitrary emittance partition, which all form essential elements for phase-space manipulations. The new generalized KV model yields spatially uniform density profiles and corresponding linear self-field forces as desired. The corresponding matrix envelope equations and beam matrix for the generalized KV model provide important new theoretical tools for the detailed design and analysis of high-intensity beam manipulations, for which previous theoretical models are not easily applicable.

preprint2016arXiv

Hamiltonian particle-in-cell methods for Vlasov-Maxwell equations

In this paper, we develop Hamiltonian particle-in-cell methods for Vlasov-Maxwell equations by applying conforming finite element methods in space and splitting methods in time. For the spatial discretisation, the criteria for choosing finite element spaces are presented such that the semi-discrete system possesses a discrete non-canonical Poisson structure. We apply a Hamiltonian splitting method to the semi-discrete system in time, then the resulting algorithm is Poisson preserving and explicit. The conservative properties of the algorithm guarantee the efficient and accurate numerical simulation of the Vlasov-Maxwell equations over long-time.

preprint2016arXiv

High order volume-preserving algorithms for relativistic charged particles in general electromagnetic fields

We construct high order symmetric volume-preserving methods for the relativistic dynamics of a charged particle by the splitting technique with processing. Via expanding the phase space to include time $t$, we give a more general construction of volume-preserving methods that can be applied to systems with time-dependent electromagnetic fields. The newly derived methods provide numerical solutions with good accuracy and conservative properties over long time of simulation. Furthermore, because of the use of processing technique the high order methods are explicit, and cost less than the methods derived from standard compositions, thus are more efficient. The results are verified by the numerical experiments. Linear stability analysis of the methods show that the high order processed method allows larger time step size during integration.

preprint2016arXiv

Largest Particle Simulations Downgrade the Runaway Electron Risk for ITER

Fusion energy will be the ultimate clean energy source for mankind. One of the most visible concerns of the future fusion device is the threat of deleterious runaway electrons (REs) produced during unexpected disruptions of the fusion plasma. Both efficient long-term algorithms and super-large scale computing power are necessary to reveal the complex dynamics of REs in a realistic fusion reactor. In the present study, we deploy the world's fastest supercomputer, Sunway TaihuLight, and the newly developed relativistic volume-preserving algorithm to carry out long-term particle simulations of 10^7 sampled REs in 6D phase space, which involves simulation scale of 10^18 particle-steps, the largest ever achieved in fusion research. Our simulations show that in a realistic fusion reactor, the concern of REs is not as serious as previously thought. Specifically, REs are confined much better than previously predicted and the maximum average energy is in the range of 150MeV, less than half of previous estimate.

preprint2016arXiv

Multi-scale Full-orbit Analysis on Phase-space Behavior of Runaway Electrons in Tokamak Fields with Synchrotron Radiation

In this paper, the secular full-orbit simulations of runaway electrons with synchrotron radiation in tokamak fields are carried out using a relativistic volume-preserving algorithm. Detailed phase-space behaviors of runaway electrons are investigated in different dynamical timescales spanning 11 orders. When looking into the small timescale, i.e., the characteristic timescale imposed by Lorentz force, the severely deformed helical trajectory of energetic runaway electron is witnessed. A qualitative analysis of the neoclassical scattering, a kind of collisionless pitch-angle scattering phenomena, is provided when considering the coupling between the rotation of momentum vector and the background magnetic field. In large timescale up to one second, it is found that the initial condition of runaway electrons in phase space globally influences the pitch-angle scattering, the momentum evolution, and the loss-gain ratio of runaway energy evidently. However, the initial value has little impact on the synchrotron energy limit. It is also discovered that the parameters of tokamak device, such as the toroidal magnetic field, the loop voltage, the safety factor profile, and the major radius, can modify the synchrotron energy limit as well as the strength of neoclassical scattering. The maximum runaway energy is also proved to be lower than the synchrotron limit when the magnetic field ripple is considered.

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.

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.

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

preprint2015arXiv

Explicit non-canonical symplectic algorithms for charged particle dynamics

We study the non-canonical symplectic structure, or K-symplectic structure inherited by the charged particle dynamics. Based on the splitting technique, we construct non-canonical symplectic methods which is explicit and stable for the long-term simulation. The key point of splitting is to decompose the Hamiltonian as four parts, so that the resulting four subsystems have the same structure and can be solved exactly. This guarantees the K-symplectic preservation of the numerical methods constructed by composing the exact solutions of the subsystems. The error convergency of numerical solutions is analyzed by means of the Darboux transformation. The numerical experiment display the long-term stability and efficiency for these methods.

preprint2015arXiv

Field theory and weak Euler-Lagrange equation for classical particle-field systems

It is commonly believed as a fundamental principle that energy-momentum conservation of a physical system is the result of space-time symmetry. However, for classical particle-field systems, e.g., charged particles interacting through self-consistent electromagnetic or electrostatic fields, such a connection has only been cautiously suggested. It has not been formally established. The difficulty is due to the fact that the dynamics of particles and the electromagnetic fields reside on different manifolds. We show how to overcome this difficulty and establish the connection by generalizing the Euler-Lagrange equation, the central component of a field theory, to a so-called weak form. The weak Euler-Lagrange equation induces a new type of flux, called the weak Euler-Lagrange current, which enters conservation laws. Using field theory together with the weak Euler-Lagrange equation developed here, energy-momentum conservation laws that are difficult to find otherwise can be systematically derived from the underlying space-time symmetry.

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.

preprint2015arXiv

Neoclassical Pitch-Angle Scattering of Runaway Electrons

It is discovered that the tokamak field geometry generates a pitch-angle scattering effect for runaway electrons. This neoclassical pitch-angle scattering is much stronger than the collisional scattering and invalidates the gyro-center model for runaway electrons. As a result, the energy limit of runaway electrons is found to be larger than the prediction of the gyro-center model and to depend heavily on the background magnetic field.

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.

preprint2015arXiv

What is the fate of runaway positrons in tokamaks?

Massive runaway positrons are generated by runaway electrons in tokamaks. The fate of these positrons encodes valuable information about the runaway dynamics. The phase space dynamics of a runaway position is investigated using a Lagrangian that incorporates the tokamak geometry, loop voltage, radiation and collisional effects. It is found numerically that runaway positrons will drift out of the plasma to annihilate on the first wall, with an in-plasma annihilation possibility less than 0.1\%. The dynamics of runaway positrons provides signatures that can be observed as diagnostic tools.

preprint2014arXiv

Asymptotic normality in the maximum entropy models on graphs with an increasing number of parameters

Maximum entropy models, motivated by applications in neuron science, are natural generalizations of the $β$-model to weighted graphs. Similar to the $β$-model, each vertex in maximum entropy models is assigned a potential parameter, and the degree sequence is the natural sufficient statistic. Hillar and Wibisono (2013) has proved the consistency of the maximum likelihood estimators. In this paper, we further establish the asymptotic normality for any finite number of the maximum likelihood estimators in the maximum entropy models with three types of edge weights, when the total number of parameters goes to infinity. Simulation studies are provided to illustrate the asymptotic results.

preprint2014arXiv

Gauge properties of the guiding center variational symplectic integrator

Variational symplectic algorithms have recently been developed for carrying out long-time simulation of charged particles in magnetic fields. As a direct consequence of their derivation from a discrete variational principle, these algorithms have very good long-time energy conservation, as well as exactly preserving discrete momenta. We present stability results for these algorithms, focusing on understanding how explicit variational integrators can be designed for this type of system. It is found that for explicit algorithms an instability arises because the discrete symplectic structure does not become the continuous structure in the $t \rightarrow 0$ limit. We examine how a generalized gauge transformation can be used to put the Lagrangian in the "antisymmetric discretization gauge," in which the discrete symplectic structure has the correct form, thus eliminating the numerical instability. Finally, it is noted that the variational guiding center algorithms are not electromagnetically gauge invariant. By designing a model discrete Lagrangian, we show that the algorithms are approximately gauge invariant as long as $\boldsymbol{A}$ and $ϕ$ are relatively smooth. A gauge invariant discrete Lagrangian is very important in a variational particle-in-cell algorithm where it ensures current continuity and preservation of Gauss's law

preprint2014arXiv

Geometric integration of the Vlasov-Maxwell system with a variational particle-in-cell scheme

A fully variational, unstructured, electromagnetic particle-in-cell integrator is developed for integration of the Vlasov-Maxwell equations. Using the formalism of Discrete Exterior Calculus, the field solver, interpolation scheme and particle advance algorithm are derived through minimization of a single discrete field theory action. As a consequence of ensuring that the action is invariant under discrete electromagnetic gauge transformations, the integrator exactly conserves Gauss's law.

preprint2014arXiv

Variational integration for ideal magnetohydrodynamics with built-in advection equations

Newcomb's Lagrangian for ideal magnetohydrodynamics (MHD) in Lagrangian labeling is discretized using discrete exterior calculus. Variational integrators for ideal MHD are derived thereafter. Besides being symplectic and momentum-preserving, the schemes inherit built-in advection equations from Newcomb's formulation, and therefore avoid solving them and the accompanying error and dissipation. We implement the method in 2D and show that numerical reconnection does not take place when singular current sheets are present. We then apply it to studying the dynamics of the ideal coalescence instability with multiple islands. The relaxed equilibrium state with embedded current sheets is obtained numerically.

preprint2014arXiv

VEST: abstract vector calculus simplification in Mathematica

We present a new package, VEST (Vector Einstein Summation Tools), that performs abstract vector calculus computations in Mathematica. Through the use of index notation, VEST is able to reduce three-dimensional scalar and vector expressions of a very general type to a well defined standard form. In addition, utilizing properties of the Levi-Civita symbol, the program can derive types of multi-term vector identities that are not recognized by reduction, subsequently applying these to simplify large expressions. In a companion paper (Burby et al., 2013), we employ VEST in the automation of the calculation of high-order Lagrangians for the single particle guiding center system in plasma physics, a computation which illustrates its ability to handle very large expressions. VEST has been designed to be simple and intuitive to use, both for basic checking of work and more involved computations.

preprint2013arXiv

4-Dimensional Geometry Lens: A Novel Volumetric Magnification Approach

We present a novel methodology that utilizes 4-Dimensional (4D) space deformation to simulate a magnification lens on versatile volume datasets and textured solid models. Compared with other magnification methods (e.g., geometric optics, mesh editing), 4D differential geometry theory and its practices are much more flexible and powerful for preserving shape features (i.e., minimizing angle distortion), and easier to adapt to versatile solid models. The primary advantage of 4D space lies at the following fact: we can now easily magnify the volume of regions of interest (ROIs) from the additional dimension, while keeping the rest region unchanged. To achieve this primary goal, we first embed a 3D volumetric input into 4D space and magnify ROIs in the 4th dimension. Then we flatten the 4D shape back into 3D space to accommodate other typical applications in the real 3D world. In order to enforce distortion minimization, in both steps we devise the high dimensional geometry techniques based on rigorous 4D geometry theory for 3D/4D mapping back and forth to amend the distortion. Our system can preserve not only focus region, but also context region and global shape. We demonstrate the effectiveness, robustness, and efficacy of our framework with a variety of models ranging from tetrahedral meshes to volume datasets.

preprint2013arXiv

A Lagrangian perspective on the stability of ideal MHD equilibria with flow

We take a careful look at two approaches to deriving stability criteria for ideal MHD equilibria. One is based on a tedious analysis of the linearized equations of motion, while the other examines the second variation of the MHD Hamiltonian computed with proper variational constraints. For equilibria without flow, the two approaches are known to be fully consistent. However, for equilibria with flow, the stability criterion obtained from the constrained variation approach was claimed to be stronger than that derived using the linearized equations of motion. We show this claim is incorrect by deriving and comparing both criteria within the same framework. It turns out that the criterion obtained from the constrained variation approach has stricter requirements on the initial perturbations than the other. Such requirements naturally emerge in our new treatment of the constrained variation approach using the Euler-Poincaré structure of ideal MHD, which is more direct and simple than the previous derivation from the Poisson perspective.