Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userCharting Galactic Accelerations: When and How to Extract a Unique Potential from the Distribution Function
613
0
0.0
(
0
)
Authors
J. An
Ask ChatGPT about the research
No Arabic abstract
The advent of datasets of stars in the Milky Way with six-dimensional phase-space information makes it possible to construct empirically the distribution function (DF). Here, we show that the accelerations can be uniquely determined from the DF using the collisionless Boltzmann equation, providing the Hessian determinant of the DF with respect to the velocities is non-vanishing. We illustrate this procedure and requirement with some analytic examples. Methods to extract the potential from datasets of discrete positions and velocities of stars are then discussed. Following Green & Ting (arXiv:2011.04673), we advocate the use of normalizing flows on a sample of observed phase-space positions to obtain a differentiable approximation of the DF. To then derive gravitational accelerations, we outline a semi-analytic method involving direct solutions of the over-constrained linear equations provided by the collisionless Boltzmann equation. Testing our algorithm on mock datasets derived from isotropic and anisotropic Hernquist models, we obtain excellent accuracies even with added noise. Our method represents a new, flexible and robust means of extracting the underlying gravitational accelerations from snapshots of six-dimensional stellar kinematics of an equilibrium system.
rate research
Read More
We develop a novel method to simultaneously determine the vertical potential, force and stellar $z-v_z$ phase space distribution function (DF) in our local patch of the Galaxy. We assume that the Solar Neighborhood can be treated as a one-dimensional system in dynamical equilibrium and directly fit the number density in the $z-v_z$ plane to what we call the Rational Linear DF (RLDF) model. This model can be regarded as a continuous sum of isothermal DFs though it has only one more parameter than the isothermal model. We apply our method to a sample of giant stars from Gaia Data Release 2 and show that the RLDF provides an excellent fit to the data. The well-known phase space spiral emerges in the residual map of the $z-v_z$ plane. We use the best-fit potential to plot the residuals in terms of the frequency and angle of vertical oscillations and show that the spiral maps into a straight line. From its slope, we estimate that the phase spirals were generated by a perturbation $sim540$ Myr years ago. We also determine the differential surface density as a function of vertical velocity dispersion, a.k.a. the vertical temperature distribution. The result is qualitatively similar to what was previously found for SDSS/SEGUE G dwarfs. Finally, we address parameter degeneracies and the validity of the 1D approximation. Particularly, the mid-plane density derived from a cold subsample, where the 1D approximation is more secure, is closer to literature values than that derived from the sample as a whole.
Astronomy has evolved almost exclusively by the use of spectroscopic and imaging techniques, operated separately. With the development of modern technologies it is possible to obtain datacubes in which one combines both techniques simultaneously, producing images with spectral resolution. To extract information from them can be quite complex, and hence the development of new methods of data analysis is desirable. We present a method of analysis of datacube (data from single field observations, containing two spatial and one spectral dimension) that uses PCA (Principal Component Analysis) to express the data in the form of reduced dimensionality, facilitating efficient information extraction from very large data sets. PCA transforms the system of correlated coordinates into a system of uncorrelated coordinates ordered by principal components of decreasing variance. The new coordinates are referred to as eigenvectors, and the projections of the data onto these coordinates produce images we will call tomograms. The association of the tomograms (images) to eigenvectors (spectra) is important for the interpretation of both. The eigenvectors are mutually orthogonal and this information is fundamental for their handling and interpretation. When the datacube shows objects that present uncorrelated physical phenomena, the eigenvectors orthogonality may be instrumental in separating and identifying them. By handling eigenvectors and tomograms one can enhance features, extract noise, compress data, extract spectra, etc. We applied the method, for illustration purpose only, to the central region of the LINER galaxy NGC 4736, and demonstrate that it has a type 1 active nucleus, not known before. Furthermore we show that it is displaced from the centre of its stellar bulge.
We use compiled high-precision pulsar timing measurements to directly measure the Galactic acceleration of binary pulsars relative to the Solar System barycenter. Given the vertical accelerations, we use the Poisson equation to derive the Oort limit, i.e., the total volume mass density in the Galactic mid-plane. Our best-fitting model gives an Oort limit of $0.08^{0.05}_{-0.02} M_{odot}/rm pc^{3}$, which is close to estimates from recent Jeans analyses. Given the accounting of the baryon budget from McKee et al. (2015), we obtain a local dark matter density of $-0.004^{0.05}_{-0.02}~M_{odot}/rm pc^{3}$, which is slightly below other modern estimates but consistent within the current uncertainties of our method. While this first measurement of the Oort limit (and other Galactic parameters) has error bars that are currently several times larger than kinematical estimates, they should improve in the future. We also constrain the oblateness of the potential, finding it consistent with that expected from the disk and inconsistent with a potential dominated by a spherical halo, as is appropriate for our sample which is within a $sim$ kpc of the Sun. We find that the slope of the rotation curve is not constrained by current measurements of binary pulsar accelerations. We give a fitting function for the vertical acceleration $a_{z}$: $a_{z} = -alpha_{1}z$; $log_{10} (alpha_{1}/{rm Gyr}^{-2})=3.69^{0.19}_{-0.12}$. By analyzing interacting simulations of the Milky Way, we find that large asymmetric variations in $da_{z}/dz$ as a function of vertical height may be a signature of sub-structure. We end by discussing the power of combining constraints from pulsar timing and high-precision radial velocity (RV) measurements towards lines-of-sight near pulsars, to test theories of gravity and constrain dark matter sub-structure.
We develop an approach for fitting the results of modeling of wriggling radial large scale iron pattern along the Galactic disk, derived over young (high massive) Cepheids, with the metallicity distribution, obtained using low mass long living dwarf stars in the close solar vicinity. For this, at the step of computing of the theoretical abundance distribution over low mass stars in the solar vicinity we propose to redefine the initial mass function so as the resulting theoretical stellar distribution over masses would be close to the distribution in the observed sample. By means of the above algorithm and subsequent corrections of the theoretical metallicity distribution function, described in literature, we have achieved fairly well agreement of the theoretical and observed metallicity distribution functions for low mass stars in the local solar vicinity.
We present maps of the stellar streams detected in the Gaia Data Release 2 (DR2) and Early Data Release 3 (EDR3) catalogs using the STREAMFINDER algorithm. We also report the spectroscopic follow-up of the brighter DR2 stream members obtained with the high-resolution CFHT/ESPaDOnS and VLT/UVES spectrographs as well as with the medium-resolution NTT/EFOSC2 spectrograph. Two new stellar streams that do not have a clear progenitor are detected in DR2 (named Hrid and Gunnthra), and seven are detected in EDR3 (named Gaia-6 to Gaia-12). Several candidate streams are also identified. The software also finds very long tidal tails associated with the 15 globular clusters NGC 288, NGC 1261, NGC 1851, NGC 2298, NGC 2808, NGC 3201, M 68, $omega$Cen, NGC 5466, Palomar 5, M 5, NGC 6101, M 92, NGC 6397 and NGC 7089. These stellar streams will be used in subsequent contributions in this series to chart the properties of the Galactic acceleration field on $sim$ 100 pc to $sim$ 100 kpc scales.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
