Related papers: Implementing Few-Body Algorithmic Regularization w…
Numerical evolution of the spherically symmetric, massive Klein-Gordon field is presented using a new adaptive mesh refinement (AMR) code with fourth order discretization in space and time, along with compactification in space. The system…
The two-dimensional n-body problem of classical mechanics is a non-integrable Hamiltonian system for n > 2. Traditional numerical integration algorithms, which are polynomials in the time step, typically lead to systematic drifts in the…
The theory of the post-Newtonian (PN) planar circular restricted three-body problem is used for numerically investigating the orbital dynamics of a test particle (e.g., a comet, asteroid, meteor or spacecraft) in the planar Sun-Jupiter…
We describe novel subgradient methods for a broad class of matrix optimization problems involving nuclear norm regularization. Unlike existing approaches, our method executes very cheap iterations by combining low-rank stochastic…
In the present paper, using the first-order approximation of the $n$-body Lagrangian (derived on the basis of the post-Newtonian gravitational theory of Einstein, Infeld, and Hoffman), we explicitly write down the equations of motion for…
Aims: We study the impact of Post-Newtonian correction terms on the energetic interaction between a gravitational wave (GW)-emitting supermassive black hole (SMBH) binary system and incoming stars via three-body scattering experiments.…
We consider an analytic way to make the interacting N-body problem tractable by using harmonic oscillators in place of the relevant two-body interactions. The two body terms of the N-body Hamiltonian are approximated by considering the…
The post-Newtonian (PN) perturbative framework has been successful in understanding the slow-motion, weak field limit of Einstein's theory of gravity on solar system scales, and for isolated astrophysical systems. The parameterized…
The numerical integration of particle trajectories in curved spacetimes is fundamental for obtaining realistic models of the particle dynamics around massive compact objects such as black holes and neutron stars. Generalized algorithms…
In this paper, we utilize stochastic optimization to reduce the space complexity of convex composite optimization with a nuclear norm regularizer, where the variable is a matrix of size $m \times n$. By constructing a low-rank estimate of…
Black hole binaries on non-eccentric orbits form an important subclass of gravitational wave sources, but it is a non-trivial issue to construct numerical initial data with minimal initial eccentricity for numerical simulations. We compute…
We have developed a simulation code with the techniques which enhance both spatial and time resolution of the PM method for which the spatial resolution is restricted by the spacing of structured mesh. The adaptive mesh refinement (AMR)…
The $N$-body problem is of historical significance because it was the first implementation of the Newtonian dynamical laws for the description of our Solar System. Motivated by this, the project's goal is to revisit this problem for small…
Overdetermined systems of first kind integral equations appear in many applications. When the right-hand side is discretized, the resulting finite-data problem is ill-posed and admits infinitely many solutions. We propose a numerical method…
We initiate a systematic study of utilizing predictions to improve over approximation guarantees of classic algorithms, without increasing the running time. We propose a systematic method for a wide class of optimization problems that ask…
We introduce a numerically exact real-time evolution scheme for quantum impurities in a macroscopically large bath. The algorithm is few-body revealing, namely it identifies the electronic orbitals that can be made inactive (in a trivial…
In this work, we consider the matrix completion problem, where the objective is to reconstruct a low-rank matrix from a few observed entries. A commonly employed approach involves nuclear norm minimization. For this method to succeed, the…
We discuss a new N-body simulation method for studying black hole binary dynamics. This method avoids previous numerical problems due to large mass ratios and trapped orbits with short periods. A treatment of relativistic effects is…
The Similarity Renormalization Group (SRG) is investigated as a powerful yet practical method to modify nuclear potentials so as to reduce computational requirements for calculations of observables. The key feature of SRG transformations…
One of the central challenges in the study of quantum many-body systems is the complexity of simulating them on a classical computer. A recent advance of Landau et al. gave a polynomial time algorithm to actually compute a succinct…