Related papers: CheasePy
The current state-of-the-art methods for showing inapproximability in PPAD arise from the $\varepsilon$-Generalized-Circuit ($\varepsilon$-GCircuit) problem. Rubinstein (2018) showed that there exists a small unknown constant $\varepsilon$…
A python code (mvp.py) is presented for computing the mean-value point (MVP) in the Brillouin zone first introduced by Baldereschi [1]. The code allows calculations of the MVP for any input crystal structure. Having MVP allows approximating…
In this paper we present a computer-assisted procedure for proving the existence of transverse heteroclinic orbits connecting hyperbolic equilibria of polynomial vector fields. The idea is to compute high-order Taylor approximations of…
Energetic particle (EP) driven instabilities will be of strongly increased relevance in future burning plasmas as the EP pressure will be very large compared to the thermal plasma. Understanding the interaction of EPs and bulk plasma is…
An asymptotic preserving and energy stable scheme for the Euler-Poisson system under the quasineutral scaling is designed and analysed. Correction terms are introduced in the convective fluxes and the electrostatic potential, which lead to…
A necessary and sufficient set of conditions for a quasisymmetric magnetic field in the form of constraint equations is derived from first principles. Without any assumption regarding the magnetohydrodynamic (MHD) equilibrium of the plasma,…
Calculations of tearing mode stability in tokamaks split conveniently into an external region, where marginally stable ideal MHD is applicable, and a resonant layer around the rational surface where sophisticated kinetic physics is needed.…
QPot is an R package for analyzing two-dimensional systems of stochastic differential equations. It provides users with a wide range of tools to simulate, analyze, and visualize the dynamics of these systems. One of QPot's key features is…
QAssemble is a pure-Python package for the quantum many-body problem. It implements various functional approaches, such as tight-binding, Hartree-Fock, and GW approximations within a unified object-oriented architecture. Each physical…
In this paper, we introduce and discuss an exactly energy-conserving Particle-in-Cell method for arbitrary curvilinear coordinates. The flexibility provided by curvilinear coordinates enables the study of plasmas in complex-shaped domains…
In this article we expand and develop the authors' recent proposed methodology for efficient stochastic superparameterization (SP) algorithms for geophysical turbulence. Geophysical turbulence is characterized by significant intermittent…
A comprehensive numerical solution is developed for the transient thermo-elasto-hydrodynamic (TEHD) characteristics of mechanical face seals. Transient lubrication features of the fluid film, transient thermal deformation features of the…
Equations describing plasma equilibria are derived from the total energy of the system. The MHD equilibrium is shown to hold even for systems where the magnetic field may locally vanish. Through conservation of helicity, confinement is…
Motivated by the work of Cooper & Showman, we revisit the chemical relaxation method, which seeks to enhance the computational efficiency of chemical-kinetics calculations by replacing the chemical network with a handful of independent…
This article introduces a software package release for geometrically reasoning about the \textit{safety} desiderata of (complex) dynamical systems via level set methods. In emphasis, safety is analyzed with Hamilton-Jacobi equations. In…
We consider a topological stabilizer code on a honeycomb grid, the "XYZ$^2$" code. The code is inspired by the Kitaev honeycomb model and is a simple realization of a "matching code" discussed by Wootton [J. Phys. A: Math. Theor. 48, 215302…
We study the equilibrium instability of chiral quarkonia in a plasma in the presence of constant magnetic and electric field and at finite axial chemical potential using AdS/CFT duality. The model in use is a supersymmetric QCD at large…
The Vlasov equation is a nonlinear partial differential equation that provides a first-principles description of the dynamics of plasmas. Its linear limit is routinely used in plasma physics to investigate plasma oscillations and stability.…
In this work we propose a new algorithm for the computation of statistical equilibrium quantities on a cubic lattice when both an energy and a statistical temperature are involved. We demonstrate that the pivot algorithm used in situations…
3D equilibrium codes are vital for stellarator design and operation, and high-accuracy equilibria are also necessary for stability studies. This paper details comparisons of two three-dimensional equilibrium codes, VMEC, which uses a…