相关论文: Two-scale exponential integrators with uniform acc…
We propose a family of reliable symplectic integrators adapted to the Discrete Non-Linear Schr\"odinger equation; based on an idea of Yoshida (H. Yoshida, Construction of higher order symplectic integrators, Physics Letters A, 150, 5,6,7,…
This article is a contribution to the study of superintegrable Hamiltonian systems with magnetic fields on the three-dimensional Euclidean space $\mathbb{E}_3$ in quantum mechanics. In contrast to the growing interest in complex…
We show how to increase the order of one-dimensional discrete gradient numerical integrator without losing its advantages, such as exceptional stability, exact conservation of the energy integral and exact preservation of the trajectories…
We derive the zero order approximation of a charged particle under the influence of a strong magnetic field in a mathematically rigorous manner and clarify in which sense this approximation is valid. We use this to further rigorously derive…
Recently, triple decomposition has attracted increasing attention for decomposing third-order tensors into three factor tensors. However, this approach is limited to third-order tensors and enforces uniformity in the lower dimensions across…
The nonlinear capacitors are important devices in modern technologies and applied physics. The aim of this paper is to calculate exactly the capacitance and the electrostatic potential energy of a nonlinear parallel-plate capacitor by using…
We optimize Hockney and Eastwood's Particle-Particle Particle-Mesh (P3M) algorithm to achieve maximal accuracy in the electrostatic energies (instead of forces) in 3D periodic charged systems. To this end we construct an optimal influence…
Due to the high computational cost of 3D particle-in-cell (PIC) simulations, lower-dimensional (2D or 1D) simulations are frequently used in their place. Our work shows that when modeling high-intensity laser ion acceleration, simulation…
We propose a second order exponential scheme suitable for two-component coupled systems of stiff evolutionary advection--diffusion--reaction equations in two and three space dimensions. It is based on a directional splitting of the involved…
A comprehensive linear stability analysis of force-gradient integrators and their Hessian-free variants is carried out by investigating the harmonic oscillator as a test equation. The analysis reveals that the linear stability of…
In this paper, we propose a $\mu$-mode integrator for computing the solution of stiff evolution equations. The integrator is based on a $d$-dimensional splitting approach and uses exact (usually precomputed) one-dimensional matrix…
By means of hybrid multi-particle collsion--particle-in-cell (MPC-PIC) simulations we study the dynamical scaling of energy and density correlations at equilibrium in moderately coupled 2D and quasi 1D plasmas. We find that the predictions…
We present the integrand decomposition of multiloop scattering amplitudes in parallel and orthogonal space-time dimensions, $d=d_\parallel+d_\perp$, being $d_\parallel$ the dimension of the parallel space spanned by the legs of the…
Integrators are fundamental instruments to recover differential signals from magnetic probes in Experimental Advanced Superconducting Tokamak (EAST) experiments. A kind of difference integrator is introduced which has the same structure as…
In the framework of a Particle-In-Cell scheme for some 1D Vlasov-Poisson system depending on a small parameter, we propose a time-stepping method which is numerically uniformly accurate when the parameter goes to zero. Based on an…
A recent proof-of-principle study proposes a nonlinear electrostatic implicit particle-in-cell (PIC) algorithm in one dimension (Chen, Chacon, Barnes, J. Comput. Phys. 230 (2011) 7018). The algorithm employs a kinetically enslaved…
Space plasma simulations have seen an increase in the use of magnetohydrodynamic (MHD) with embedded Particle-in-Cell (PIC) models. This combined MHD-EPIC algorithm simulates some regions of interest using the kinetic PIC method while…
High-fidelity modeling of plasma-based acceleration (PBA) requires the use of 3D fully nonlinear and kinetic descriptions based on the particle-in-cell (PIC) method. Three-dimensional PIC algorithms based on the quasi-static approximation…
A new family of methods involving complex coefficients for the numerical integration of differential equations is presented and analyzed. They are constructed as linear combinations of symmetric-conjugate compositions obtained from a basic…
In molecular dynamics, penalized overdamped Langevin dynamics are used to model the motion of a set of particles that follow constraints up to a parameter $\varepsilon$. The most used schemes for simulating these dynamics are the Euler…