Related papers: Basis function expansions for galactic dynamics: S…
We construct numerical models of mildly triaxial elliptical galaxies with central density cusps. Using a technique we call ``adiabatic squeezing'', we begin with a spherical gamma=1 Hernquist model and apply a drag to the velocities of the…
The past years have witnessed impressive advances in electronic structure calculation, especially in the complexity and size of the systems studied, as well as in computation time. Linear scaling methods based on empirical tight-binding…
Given a point $A$ in the convex hull of a given adjoint orbit $\mathcal{O}(F)$ of a compact Lie group $G$, we give a polynomial time algorithm to compute the probability density supported on $\mathcal{O}(F)$ whose expectation is $A$ and…
We report an efficient algorithm using density fitting for the relativistic complete active space self-consistent field (CASSCF) method, which is significantly more stable than the algorithm previously reported by one of the authors [J. E.…
A number of methods for studying the large-scale cosmic matter distribution exist in the literature. One particularly common method employed to define the cosmic web is to examine the density, velocity or potential field. Such methods are…
We present a method to numerically estimate the densities of a discretely sampled data based on binary space partitioning tree. We start with a root node containing all the particles and then recursively divide each node into two nodes each…
Classification of galactic morphologies is a crucial task in galactic astronomy, and identifying fine structures of galaxies (e.g., spiral arms, bars, and clumps) is an essential ingredient in such a classification task. However, seeing…
We use techniques from nonparametric function estimation theory to extract the density profiles, and their derivatives, from a set of N-body dark matter halos. We consider halos generated from LCDM simulations of gravitational clustering,…
With the emergence of Artificial Intelligence, numerical algorithms are moving towards more approximate approaches. For methods such as PCA or diffusion maps, it is necessary to compute eigenvalues of a large matrix, which may also be dense…
Self-force methods can be applied in calculations of the scatter angle in two-body hyperbolic encounters, working order by order in the mass ratio (assumed small) but with no recourse to a weak-field approximation. This, in turn, can inform…
We construct time-evolving gravitational potential models for a Milky Way-mass galaxy from the FIRE-2 suite of cosmological-baryonic simulations using basis function expansions. These models capture the angular variation with spherical…
Cosmological tests based on cluster counts require accurate calibration of the space density of massive halos, but most calibrations to date have ignored complex gas physics associated with halo baryons. We explore the sensitivity of the…
This article describes the Schwarzschild orbit superposition method. It is the state-of-the-art dynamical modelling tool for early-type galaxies. Tests with analytic models show that masses and orbital anisotropies of not too face-on…
We present a new method that simultaneously solves for cosmology and galaxy bias on non-linear scales. The method uses the halo model to analytically describe the (non-linear) matter distribution, and the conditional luminosity function…
We revisit the J-matrix method for the one dimensional radial harmonic oscillator (RHO) and construct its tridiagonal matrix representation within an orthonormal basis phi(z)n of L2 (R+);parametrized by a fixed z in the complex unit disc D…
Current shallow granular flow models suited to arbitrary topography can be divided into two types, those formulated in bed-fitted curvilinear coordinates, and those formulated in global Cartesian coordinates. The shallow granular flow model…
Using methods from symplectic geometry, the second and fifth authors have provided theoretical groundwork and tools aimed at analyzing periodic orbits, their stability and their bifurcations in families, for the purpose of space mission…
Direct N-body simulations and symplectic integrators are effective tools to study the long-term evolution of planetary systems. The Wisdom-Holman (WH) integrator in particular has been used extensively in planetary dynamics as it allows for…
The Stochastic Series Expansion (SSE) quantum Monte Carlo method with directed loops is very efficient for spin and boson systems. The Heisenberg mode l and its generalizations, such as the $JQ_2$ model, are extensively simulated via this…
This report starts by describing the continuum model used by Lingle & Clark (1985) to approximate the deformation of the earth under changing ice sheet and ocean loads. That source considers a single ice stream, but we apply their…