Related papers: SCALAR: an AMR code to simulate axion-like dark ma…
We extend our previous work on a semi-Lagrangian adaptive rank (SLAR) integrator, in the finite difference framework for nonlinear Vlasov-Poisson systems, to the general high-order tensor setting. The proposed scheme retains the high-order…
Obtainable computational efficiency is evaluated when using an Adaptive Mesh Refinement (AMR) strategy in time accurate simulations governed by sets of conservation laws. For a variety of 1D, 2D, and 3D hydro- and magnetohydrodynamic…
We present a new massively parallel code for N-body and cosmological hydrodynamical simulations of modified gravity models. The code employs a multigrid-accelerated Newton-Gauss-Seidel relaxation solver on an adaptive mesh to efficiently…
A covariant non-local extention if the stationary Schr\"odinger equation is presented and it's solution in terms of Heisenbergs's matrix quantum mechanics is proposed. For the special case of the Riesz fractional derivative, the calculation…
We develop a quantum algorithm for solving high-dimensional fractional Poisson equations. By applying the Caffarelli-Silvestre extension, the $d$-dimensional fractional equation is reformulated as a local partial differential equation in…
Axions are candidates for dark matter in the universe.We develop an accurate Boltzmann code to calculate the linear growth of the plasma. As an interesting example, we investigate a mixed dark matter model consisting of cold dark matter…
Scalar, tensor waves induce oscillatory perturbations in Keplerian systems that can be probed with measurements of pulsar timing residuals. In this paper, we consider the imprint of coherent oscillations produced by ultralight axion dark…
The era of precision cosmology allows us to test the composition of the dark matter. Mixed ultralight or fuzzy dark matter (FDM) is a cosmological model with dark matter composed of a combination of particles of mass $m\leq…
We demonstrate, for the first time, that neural scaling laws can deliver near-exact solutions to the many-electron Schr\"odinger equation across a broad range of realistic molecules. This progress is enabled by the Lookahead Variational…
High-order solvers for compressible flows are vital in scientific applications. Adaptive mesh refinement (AMR) is a key technique for reducing computational cost by concentrating resolution in regions of interest. In this work, we develop…
We study the Schr\"odinger-Poisson (SP) method in the context of cosmological large-scale structure formation in an expanding background. In the limit $\hbar \to 0$, the SP technique can be viewed as an effective method to sample the phase…
For simulations that deal only with dark matter or stellar systems, the conventional N-body technique is fast, memory efficient, and relatively simple to implement. However when including the effects of gas physics, mesh codes are at a…
This paper presents a linear, decoupled, mass- and energy-conserving numerical scheme for the multi-dimensional coupled nonlinear Schr\"odinger (CNLS) system. The scheme combines the fourth-order compact difference approximation in space…
In this paper, we revisit an auxiliary space preconditioning method proposed by Xu [Computing 56, 1996], in which low-order finite element spaces are employed as auxiliary spaces for solving linear algebraic systems arising from high-order…
Numerical simulations of two-phase flow and fluid structure interaction problems are of great interest in many environmental problems and engineering applications. To capture the complex physical processes involved in these problems, a high…
We derive non-relativistic equations of motion for the formation of cosmological structure in a Scalar Field Dark Matter (SFDM) model corresponding to a complex scalar field endowed with a quadratic scalar potential. Starting with the full…
We study a nonlinear Schr\"odinger equation which arises as an effective single particle model in X-ray Free Electron Lasers (XFEL). This equation appears as a first-principles model for the beam-matter interactions that would take place in…
We report on a series of tests of agreement between three types of N-body simulations: PM, P$^3$M, and Tree codes. We find good agreement in both the individual and the statistical properties only on scales larger than the mean…
Efforts to achieve better accuracy in numerical relativity have so far focused either on implementing second order accurate adaptive mesh refinement or on defining higher order accurate differences and update schemes. Here, we argue for the…
To date, fully cosmological hydrodynamic disk simulations to redshift zero have only been undertaken with particle-based codes, such as GADGET, Gasoline, or GCD+. In light of the (supposed) limitations of traditional implementations of…