Related papers: Basis function expansions for galactic dynamics: S…
We devise a Hybrid High-Order (HHO) method for highly oscillatory elliptic problems that is capable of handling general meshes. The method hinges on discrete unknowns that are polynomials attached to the faces and cells of a coarse mesh;…
Tremendous progress has been made recently in modelling the morphology and kinematics of centers of galaxies. Increasingly realistic models are built for central bar, bulge, nucleus and black hole of galaxies, including our own. The newly…
Using one-center expansion relations for the Slater type orbitals (STOs) of noninteger principal quantum numbers in terms of integer n STOs derived in this study with the help of - exponential type orbitals (-ETOs, the general formulas are…
(Abridged:) Schwarzschild's method was used to construct equilibrium solutions to the collisionless Boltzmann equation corresponding to a Plummer sphere, which were compared with analytical results to test the robustness of the numerical…
Gravitational features are a fundamental source of information to learn more about the interior structure and composition of planets, moons, asteroids and comets. Gravitational field modeling typically approximates the target body with a…
The dynamics and evolution of any galactic structure are strongly influenced by the properties of the orbits that constitute it. In this paper, we compare two orbit classification schemes, one by Laskar (NAFF) , and the other by Carpintero…
This paper explores orbits in extended mass distributions and develops an analytic approximation scheme based on epicycloids (spirograph patterns). We focus on the Hernquist potential which provides a good model for many astrophysical…
I examine differences in non-linear structure formation between cosmological models that share a $z=0$ linear power spectrum in both shape and amplitude, but that differ via their growth history. $N$-body simulations of these models display…
We present GPU implementations of two fast force calculation methods, based on series expansions of the Poisson equation. One is the Self-Consistent Field (SCF) method, which is a Fourier-like expansion of the density field in some basis…
This paper presents two implementations of the recently developed motion planning algorithms HyRRT arXiv:2210.1508(2) and HySST arXiv:2305.1864(9). Specifically, cHyRRT, an implementation of the HyRRT algorithm, generates solutions to…
We introduce a closed-form expansion for the transition density of elliptic and hypo-elliptic multivariate Stochastic Differential Equations (SDEs), over a period $\Delta\in (0,1)$, in terms of powers of $\Delta^{j/2}$, $j\ge 0$. Our…
The analytical phase-space distribution function (DF) of spherical self--consistent galaxy (or cluster) models, embedded in a dark matter halo, where both density distributions follow the Hernquist profile, with different total masses and…
Numerical integration on spheres, including the computation of the areas of spherical triangles, is a core computation in geomathematics. The commonly used techniques sometimes suffer from instabilities and significant loss of accuracy. We…
Multi-revolution elliptic Halo (ME-Halo) orbits are a special class of symmetric and periodic solutions within the framework of the elliptic restricted three-body problem (ERTBP). During a single period, an M:N ME-Halo orbit completes $M$…
Imaginary-time path integral (PI) is a rigorous tool to compute static properties at finite temperatures. However, the stiff PI internal modes poses a sampling challenge. This is commonly tackled using staging coordinates, in which the free…
We show that generalized spherical harmonics are well suited for representing the space and orientation molecular density in the resolution of the molecular density functional theory. We consider the common system made of a rigid solute of…
We describe a new algorithm for including the dynamical effects of ionizing radiation in SPH simulations, and we present several examples of how the algorithm can be applied to problems in star formation. We use the HEALPix software to…
We devise and analyze two hybrid high-order (HHO) methods for the numerical approximation of the biharmonic problem. The methods support polyhedral meshes, rely on the primal formulation of the problem, and deliver $O(h^{k+1})$ $H^2$-error…
Near-rectilinear halo orbits (NRHOs) around Earth-Moon L2 in the Circular Restricted 3-Body Problem (CR3BP) exhibit a complex dynamical landscape, featuring a band of normally elliptic orbits embedded within regions of strong instability.…
We take into account the dynamics of three types of models of rotating galaxies in polar coordinates in a rotating frame. Due to non-axisymmetric potential perturbations, the angular momentum varies with time, and the kinetic energy depends…