Variational Schemes and Geometric Simulations for a Hydrodynamic-Electrodynamic Model of Surface Plasmon Polaritons
arXiv:1810.06205 · doi:10.1103/PhysRevE.99.023313
Abstract
A class of variational schemes for the hydrodynamic-electrodynamic model of lossless free-electron gas in a quasineutral background is developed for high-quality simulations of surface plasmon polaritons. The Lagrangian density of lossless free-electron gas with a self-consistent electromagnetic field is established, and the dynamical equations with the associated constraints are obtained via a variational principle. Based on discrete exterior calculus, the action functional of this system is discretized and minimized to obtain the discrete dynamics. Newton-Raphson iteration and the biconjugate gradient stabilized method are equipped as a hybrid nonlinear-linear algebraic solver. Instead of discretizing the partial differential equations, the variational schemes have better numerical properties in secular simulations, as they preserve the discrete Lagrangian symplectic structure, gauge symmetry, and general energy-momentum density. Two numerical experiments were performed. The numerical results reproduce characteristic dispersion relations of bulk plasmons and surface plasmon polaritons, and the numerical errors of conserved quantities in all experiments are bounded by a small value after long-term simulations.
arXiv comment9 figures
References in corpus (12)
- Finite element exterior calculus: from Hodge theory to numerical stability
- GEMPIC: Geometric ElectroMagnetic Particle-In-Cell Methods
- Structure and structure-preserving algorithms for plasma physics
- Explicit symplectic approximation of nonseparable Hamiltonians: algorithm and long time performance
- Geometric integration of the Vlasov-Maxwell system with a variational particle-in-cell scheme
- Canonical symplectic particle-in-cell method for long-term large-scale simulations of the Vlasov-Maxwell system
- Geometric variational crimes: Hilbert complexes, finite element exterior calculus, and problems on hypersurfaces
- Explicit high-order noncanonical symplectic algorithms for ideal two-fluid systems
- Maxwell-Hydrodynamic Model for Simulating Nonlinear Terahertz Generation from Plasmonic Metasurfaces
- Simulations of relativistic-quantum plasmas using real-time lattice scalar QED
- Canonical symplectic structure and structure-preserving geometric algorithms for Schrödinger-Maxwell systems
- Gauge properties of the guiding center variational symplectic integrator