Velocity Dispersion Of Stars Research Articles

The dual nature of the tidal tails of NGC 5904 (M5)

ABSTRACT The tangential velocity dispersion of stars belonging to the Milky Way globular cluster’s tidal tails has recently been found from N-body simulations to be a parameter that distinguishes between cored and cuspy profiles of low-mass dwarf galaxy dark matter subhaloes where that globular cluster formed, and the in situ formation scenario. In this context, we discovered that M5’s tidal tails are composed by stars at two different metallicity regimes ([Fe/H] ∼ −1.4 and −2.0 dex). The more metal-rich tidal tail stars are of the same metal content than M5’s members and have a tangential velocity dispersion that coincides with the predicted value for a cuspy formation scenario (subhalo mass ∼ 109 M⊙). The more metal-poor stars, that are found along the entire M5 tidal tails and have similar distributions to their more metal-rich counterparts in the M5 colour–magnitude diagram and orbit trajectory, have a tangential velocity dispersion that refers to a cored subhalo (mass ∼ 109 M⊙), or an in situ formation scenario. In order to reconcile the dual distribution of M5 tidal tail stars, in kinematics and chemistry, we propose that M5 collided with another more metal-poor and less-massive globular cluster anytime before or after it was accreted into the Milky Way.

A Statistical Relation between Mass, Age, and Velocity Dispersion in the Solar Neighborhood* *Based on the third data release of the Gaia space telescope.

The stellar kinematics of the Galactic disk are the main factors constraining disk formation and evolution processes in the Milky Way. In this paper, we investigate a statistical relation between stellar mass, age, and velocity dispersion for stars in the solar neighborhood. Age–velocity dispersion relations, with their applications, have been studied in detail before, but their correlation with mass was mostly neglected. To investigate this relation, we use proper motion data of more than 113,035 stars in the Galactic disk (with solar distances less than 150 pc) provided by the third data release of the Gaia mission, and for stellar mass and age, Gaia's Final Luminosity Age Mass Estimator is implemented. We analyze this data and the correlations between the parameters with random forest regression, which is an ensemble statistical-learning technique. Finally, we show that, by considering stellar mass alongside age, we can determine velocity dispersions with an average relative error and a mean absolute error of about 9% and 2.68 km s−1, respectively. We also find that the correlation of stellar age with velocity dispersion is 3–8 times more than mass, which varies due to the different stellar types and masses.

The Hubble Space Telescope UV Legacy Survey of Galactic Globular Clusters. XXIV. Differences in Internal Kinematics of Multiple Stellar Populations

Our understanding of the kinematic properties of multiple stellar populations (mPOPs) in Galactic globular clusters (GCs) is still limited compared to what we know about their chemical and photometric characteristics. Such limitation arises from the lack of a comprehensive observational investigation of this topic. Here we present the first homogeneous kinematic analysis of mPOPs in 56 GCs based on high-precision proper motions computed with Hubble Space Telescope data. We focused on red-giant-branch stars, for which the mPOP tagging is clearer, and measured the velocity dispersion of stars belonging to first (1G) and second generations (2G). We find that 1G stars are generally kinematically isotropic even at the half-light radius, whereas 2G stars are isotropic at the center and become radially anisotropic before the half-light radius. The radial anisotropy is induced by a lower tangential velocity dispersion of 2G stars with respect to the 1G population, while the radial component of the motion is comparable. We also show possible evidence that the kinematic properties of mPOPs are affected by the Galactic tidal field, corroborating previous observational and theoretical results suggesting a relation between the strength of the external tidal field and some properties of mPOPs. Although limited to the GCs’ central regions, our analysis leads to new insights into the mPOP phenomenon, and provides the motivation for future observational studies of the internal kinematics of mPOPs.

The origin of ultramassive white dwarfs: hints from Gaia EDR3

ABSTRACT Gaia Data Release 2 revealed a population of ultramassive white dwarfs on the Q branch that are moving anomalously fast for a local disc population with their young photometric ages. As the velocity dispersion of stars in the local disc increases with age, a proposed explanation of these white dwarfs is that they experience a cooling delay that causes current cooling models to infer photometric ages much younger than their true ages. To explore this explanation, we investigate the kinematics of ultramassive white dwarfs within 200 pc of the Sun using the improved Gaia Early Data Release 3 observations. We analyse the transverse motions of 0.95–1.25 M⊙ white dwarfs, subdivided by mass and age, and determine the distributions of the three-dimensional components of the transverse velocities. The results are compared to expectations based on observed kinematics of local main-sequence stars. We find a population of photometrically young (∼0.5–1.5 Gyr) ultramassive (∼1.15–1.25 M⊙) white dwarfs for which the transverse velocity component in the direction of Galactic rotation is more dispersed than for local disc stars of any age; thus, it is too dispersed to be explained by any cooling delay in white dwarfs originating from the local disc. Furthermore, the dispersion ratio of the velocity components in the Galactic plane for this population is also inconsistent with a local disc origin. We discuss some possible explanations of this kinematically anomalous population, such as a halo origin or production through dynamical effects of stellar triple systems.

Impact of Galaxy Dynamics on Modified Gravity Constraints from Strong Lensing Systems

AbstractGalaxy-galaxy strong lensing systems have been used in literature to test General Relativity by constraining the post-Newtonian parameter (γPPN). Nevertheless, these methods are prone to systematic errors, some of which arise from difficulties in modelling the dynamics of the lens galaxy. In this study, we address the systematic error related to the assumption of a constant anisotropy between the radial and tangential components of the velocity dispersion of stars within the lens galaxy, characterised by the parameter β. We considered two radial models for the anisotropy parameter, the Osipkov-Merritt and Mamon & Lokas, as well three Gaussian priors for constant β. Our analysis showed that the choice of β has a strong impact on the value of γPPN, with radial models leading to lower values of this parameter.

A Statistical Detection of Wide Binary Systems in the Ultrafaint Dwarf Galaxy Reticulum II

Binary stars can inflate the observed velocity dispersion of stars in dark matter-dominated systems such as ultrafaint dwarf galaxies (UFDs). However, the population of binaries in UFDs is poorly constrained by observations, with preferred binary fractions for individual galaxies ranging from a few percent to near unity. Searching for wide binaries through nearest neighbor (NN) statistics (or the two-point correlation function) has been suggested in the literature, and we apply this method for the first time to detect wide binaries in a UFD. By analyzing the positions of stars in Reticulum II (Ret II) from Hubble Space Telescope images, we search for angularly resolved wide binaries in Ret II. We find that the distribution of their NN distances shows an enhancement at projected separations of ≲1″ relative to a model containing no binaries. We show that such an enhancement can be explained by a wide binary fraction of at separations of more than 3000 au. Under the assumption that the binary separation distribution is similar to that in the Milky Way, the total binary fraction in Ret II may be on the order of 50%. We also use the observed magnitude distribution of stars in Ret II to constrain the initial mass function over the mass range 0.34–0.78 M ⊙, finding that a shallow power-law slope of 1.01 ≤ α ≤ 1.15 matches the data.

Hot Jupiters, cold kinematics

Context. The birth environments of planetary systems are thought to influence planet formation and orbital evolution through external photoevaporation and stellar flybys. Recent work has claimed observational support for this, in the form of a correlation between the properties of planetary systems and the local Galactic phase space density of the host star. In particular, hot Jupiters are overwhelmingly present around stars in regions of high phase space density, which may reflect a formation environment with high stellar density. Aims. We aim to investigate whether the high phase space density may have a Galactic kinematic origin: hot Jupiter hosts may be biased towards being young and therefore kinematically cold, because tidal inspiral leads to the destruction of the planets on gigayear timescales, and the velocity dispersion of stars in the Galaxy increases on similar timescales. Methods. We used 6D positions and kinematics from Gaia for the hot Jupiter hosts and their neighbours, and we constructed distributions of the phase space density. We investigated correlations between the stars’ local phase space density and peculiar velocity. Results. We find a strong anti-correlation between the phase space density and the host star’s peculiar velocity with respect to the Local Standard of Rest. Therefore, most stars in ‘high-density’ regions are kinematically cold, which may be caused by the aforementioned bias towards detecting hot Jupiters around young stars before the planets’ tidal destruction. Conclusions. We do not find evidence in the data for hot Jupiter hosts preferentially being in phase space overdensities compared to other stars of similar kinematics, nor therefore for their originating in birth environments of high stellar density.

Spiral structure of the galactic disk and its influence on the rotational velocity curve

AbstractMost spiral galaxies have a flat rotational velocity curve, according to the different observational techniques used in several wavelengths domain. In this work, we show that nonlinear terms are able to balance the dispersive effect of the wave, thus reviving the observed rotational curve profiles without the inclusion of any other but baryonic matter concentrated in the bulge and disk. To prove that the considered model is able to restore a flat rotational curve, Milky Way has been chosen as the best‐mapped galaxy, to apply on. Using the gravitational N‐body simulations with up to particles, we test this dynamical model in the case of the Milky Way with two different approaches. Within the direct approach, as an input condition in the simulation runs, we set the spiral surface density distribution which is previously obtained as an explicit solution to nonlinear Schrödinger equation (instead of a widely used exponential disk approximation). In the evolutionary approach, we initialize the runs with different initial mass and rotational velocity distributions, in order to capture the natural formation of spiral arms and to determine their role in the disk evolution. In both cases, we are able to reproduce the stable and nonexpanding disk structures at the simulation end times of years, with no halo inclusion. Although the given model does not take into account the velocity dispersion of stars and finite disk thickness, the results presented here still imply that nonlinear effects can significantly alter the amount of dark matter which is required to keep the galactic disk in a stable configuration.

Gamma Rays from Fast Black-hole Winds

Abstract Massive black holes at the centers of galaxies can launch powerful wide-angle winds that, if sustained over time, can unbind the gas from the stellar bulges of galaxies. These winds may be responsible for the observed scaling relation between the masses of the central black holes and the velocity dispersion of stars in galactic bulges. Propagating through the galaxy, the wind should interact with the interstellar medium creating a strong shock, similar to those observed in supernovae explosions, which is able to accelerate charged particles to high energies. In this work we use data from the Fermi Large Area Telescope to search for the γ-ray emission from galaxies with an ultrafast outflow (UFO): a fast (v ∼ 0.1 c), highly ionized outflow, detected in absorption at hard X-rays in several nearby active galactic nuclei (AGN). Adopting a sensitive stacking analysis we are able to detect the average γ-ray emission from these galaxies and exclude that it is due to processes other than UFOs. Moreover, our analysis shows that the γ-ray luminosity scales with the AGN bolometric luminosity and that these outflows transfer ∼0.04% of their mechanical power to γ-rays. Interpreting the observed γ-ray emission as produced by cosmic rays (CRs) accelerated at the shock front, we find that the γ-ray emission may attest to the onset of the wind–host interaction and that these outflows can energize charged particles up to the transition region between galactic and extragalactic CRs.

Revisiting pulsar velocities using Gaia Data Release 2

Precise measurements of neutron star (NS) velocities provide critical clues in regard to the supernova physics and evolution of binary systems. Based on Gaia Data Release 2 (DR2), we selected a sample of 24 young (< 3 Myr) pulsars with precise parallax measurements and measured the velocity of their local standard of rest (LSR) and the velocity dispersion among their respective local stellar groups. The median velocity difference between thus calculated LSRs and the Galactic rotation model is ∼ 7.6 km s−1, small compared to the typical velocity dispersion of ∼ 27.5 km s−1. For pulsars off the Galactic plane, such differences grow significantly to as large as ∼40 km s−1. More importantly, the velocity dispersion of stars in the local group of low-velocity pulsars can be comparable to their transverse velocities, suggesting that the intrinsic velocities of NS progenitors should be taken into account when we consider their natal kicks and binary evolution. We also examined the double NS system J0737–3039A/B, and measured its transverse velocity to be km s−1 assuming nearby Gaia sources are representative of their birth environment. This work demonstrated the feasibility and importance of using Gaia data to study the velocity of individual systems and velocity distribution of NSs.

Fundamental relations for the velocity dispersion of stars in the Milky Way

ABSTRACT We explore the fundamental relations governing the radial and vertical velocity dispersions of stars in the Milky Way, from combined studies of complementary surveys including GALAH, LAMOST, APOGEE, the NASA Kepler and K2 missions, and Gaia DR2. We find that different stellar samples, even though they target different tracer populations and employ a variety of age estimation techniques, follow the same set of fundamental relations. We provide the clearest evidence to date that, in addition to the well-known dependence on stellar age, the velocity dispersions of stars depend on orbital angular momentum Lz, metallicity, and height above the plane |z|, and are well described by a multiplicatively separable functional form. The dispersions have a power-law dependence on age with exponents of 0.441 ± 0.007 and 0.251 ± 0.006 for σz and σR, respectively, and the power law is valid even for the oldest stars. For the solar neighbourhood stars, the apparent break in the power law for older stars, as seen in previous studies, is due to the anticorrelation of Lz with age. The dispersions decrease with increasing Lz until we reach the Sun’s orbital angular momentum, after which σz increases (implying flaring in the outer disc) while σR flattens. For a given age, the dispersions increase with decreasing metallicity, suggesting that the dispersions increase with birth radius. The dispersions also increase linearly with |z|. The same set of relations that work in the solar neighbourhood also work for stars between 3 &amp;lt; R/kpc &amp;lt; 20. Finally, the high-[α/Fe] stars follow the same relations as the low-[α/Fe] stars.

Gyro-kinematic Ages for around 30,000 Kepler Stars

Abstract Estimating stellar ages is important for advancing our understanding of stellar and exoplanet evolution and investigating the history of the Milky Way. However, ages for low-mass stars are hard to infer as they evolve slowly on the main sequence. In addition, empirical dating methods are difficult to calibrate for low-mass stars as they are faint. In this work, we calculate ages for Kepler F, G, and crucially K and M dwarfs, using their rotation and kinematic properties. We apply the simple assumption that the velocity dispersion of stars increases over time and adopt an age–velocity-dispersion relation (AVR) to estimate average stellar ages for groupings of coeval stars. We calculate the vertical velocity dispersion of stars in bins of absolute magnitude, temperature, rotation period, and Rossby number and then convert velocity dispersion to kinematic age via an AVR. Using this method, we estimate gyro-kinematic ages for 29,949 Kepler stars with measured rotation periods. We are able to estimate ages for clusters and asteroseismic stars with an rms of 1.22 Gyr and 0.26 Gyr respectively. With our Astraea machine-learning algorithm, which predicts rotation periods, we suggest a new selection criterion (a weight of 0.15) to increase the size of the McQuillan et al. catalog of Kepler rotation periods by up to 25%. Using predicted rotation periods, we estimated gyro-kinematic ages for stars without measured rotation periods and found promising results by comparing 12 detailed age–element abundance trends with literature values.

How ‘cold’ are the stellar discs of superthin galaxies?

ABSTRACT Superthin galaxies are a class of bulgeless, low surface brightness galaxies with strikingly high values of planar-to-vertical axes ratio $\rm (b/a\gt 10\!-\!20)$, possibly indicating the presence of an ultracold stellar disc. Using the multicomponent galactic disc model of gravitationally coupled stars and gas in the force field of the dark matter halo as well as the stellar dynamical code AGAMA (Action-based Galaxy Modelling Architecture), we determine the vertical velocity dispersion of stars and gas as a function of galactocentric radius for five superthin galaxies (UGC 7321, IC 5249, FGC 1540, IC2233, and UGC00711) using observed stellar and atomic hydrogen (H i) scale heights as constraints, using a Markov chain Monte Carlo Method. We find that the central vertical velocity dispersion for the stellar disc in the optical band varies between σ0s ∼ 10.2−18.4 $\rm {kms}\,^{-1}$ and falls off with an exponential scale length of 2.6–3.2 Rd where Rd is the exponential stellar disc scale length. Interestingly, in the 3.6 μm, the same, averaged over the two components of the stellar disc, varies between 5.9 and 11.8 $\rm {kms}\,^{-1}$, both of which confirm the presence of ‘ultracold’ stellar discs in superthin galaxies. Interestingly, the global median of the multicomponent disc dynamical stability parameter QN of our sample superthins is found to be 5 ± 1.5, which higher than the global median value of 2.2 ± 0.6 for a sample of spiral galaxies.

Evaluating the merger rate of binary black holes from direct captures and third-body soft interactions using the Milky Way globular clusters

The multitude of binary black hole coalescence detections in gravitational waves has renewed our interest on environments that can be the cradle of these mergers. In this work we study merger rates of binary black holes in globular clusters that are among the most dense stellar environments and a natural place for the creation of black hole binaries. To model these systems with all their variations we rely on the observational properties of the known Milky Way globular clusters. We consider direct capture events between black holes, as well as soft interactions of black hole binaries with stars as third bodies that accelerate the evolution of these binaries. We find that binary black holes from direct captures merge at an averaged rate of $0.3-5 \times 10^{-11}$ yr$^{-1}$ per cluster. Third body soft interactions are a much more prominent channel giving an averaged rate of $2-4 \times 10^{-10}$ yr$^{-1}$ per cluster. Those rates in globular clusters can lead to a cumulative merger rate of about 100 mergers per year up to redshift of 1, i.e. a significant fraction of the detectable in the near future binary black hole coalescence events. Further observations of cluster properties both in terms of their masses, profile properties, velocity dispersion of stars and their cosmological distribution, will allow us to better constrain the contribution of these environments to the detectable coalescence events rate.

Resolving the Disc–Halo Degeneracy – II: NGC 6946

ABSTRACT The mass-to-light ratio (M/L) is a key parameter in decomposing galactic rotation curves into contributions from the baryonic components and the dark halo of a galaxy. One direct observational method to determine the disc M/L is by calculating the surface mass density of the disc from the stellar vertical velocity dispersion and the scale height of the disc. Usually, the scale height is obtained from near-IR studies of edge-on galaxies and pertains to the older, kinematically hotter stars in the disc, while the vertical velocity dispersion of stars is measured in the optical band and refers to stars of all ages (up to ∼10 Gyr) and velocity dispersions. This mismatch between the scale height and the velocity dispersion can lead to underestimates of the disc surface density and a misleading conclusion of the submaximality of galaxy discs. In this paper, we present the study of the stellar velocity dispersion of the disc galaxy NGC 6946 using integrated star light and individual planetary nebulae as dynamical tracers. We demonstrate the presence of two kinematically distinct populations of tracers that contribute to the total stellar velocity dispersion. Thus, we are able to use the dispersion and the scale height of the same dynamical population to derive the surface mass density of the disc over a radial extent. We find the disc of NGC 6946 to be closer to maximal with the baryonic component contributing most of the radial gravitational field in the inner parts of the galaxy (Vmax(bar) = 0.76(±0.14)Vmax).

Vertical stellar density distribution in a non-isothermal galactic disc

ABSTRACT The vertical density distribution of stars in a galactic disc is traditionally obtained by assuming an isothermal vertical velocity dispersion of stars. Recent observations from SDSS, LAMOST, RAVE, Gaia etc. show that this dispersion increases with height from the mid-plane. Here, we study the dynamical effect of such non-isothermal dispersion on the self-consistent vertical density distribution for the thin disc stars in the Galaxy, obtained by solving together the Poisson equation and the equation of hydrostatic equilibrium. We find that in the non-isothermal case the mid-plane density is lower and the scale height is higher than the corresponding values for the isothermal distribution, due to higher vertical pressure, hence the distribution is vertically more extended. The change is $\sim \! 35 {{\ \rm per\ cent}}$ at the solar radius for a stars-alone disc for the typical observed linear gradient of +6.7 km s−1 kpc−1 and becomes even higher with increasing radii and increasing gradients explored. The distribution shows a wing at high z, in agreement with observations, and is fitted well by a double $\operatorname{sech}^{2}$, which could be mis-interpreted as the existence of a second, thicker disc, specially in external galaxies. We also consider a more realistic disc consisting of gravitationally coupled stars and gas in the field of dark matter halo. The results show the same trend but the effect of non-isothermal dispersion is reduced due to the opposite, constraining effect of the gas and halo gravity. Further, the non-isothermal dispersion lowers the theoretical estimate of the total mid-plane density i.e. Oort limit value, by 16 per cent.

A Scenario for Ultradiffuse Satellite Galaxies with Low Velocity Dispersions: The Case of [KKS 2000]04

Abstract A scenario for achieving a low velocity dispersion for the galaxy [KKS 2000]04 (aka NGC 1052-DF2) and similar galaxies is presented. A progenitor halo corresponding to a z = 0 halo of mass and a low concentration parameter (but consistent with cosmological simulations) infalls onto a Milky Way–size host at early times. Substantial removal of cold gas from the inner regions by supernova feedback and ram pressure, assisted by tidal stripping of the dark matter in the outer regions, leads to a substantial reduction of the velocity dispersion of stars within one effective radius. In this framework, the observed stellar content of [KKS 2000]04 is associated with a progenitor mass close to that inferred from the global stellar-to-halo-mass ratio. As far as the implications of kinematics are concerned, even if at a ∼20 Mpc distance, it is argued that [KKS 2000]04 is no more peculiar than numerous early type galaxies with seemingly little total dark-matter content.

Clocking the formation of today’s largest galaxies: wide field integral spectroscopy of brightest cluster galaxies and their surroundings

ABSTRACT The formation and evolution of local brightest cluster galaxies (BCGs) is investigated by determining the stellar populations and dynamics from the galaxy core, through the outskirts and into the intracluster light (ICL). Integral spectroscopy of 23 BCGs observed out to $4\, r_{e}$ is collected and high signal-to-noise regions are identified. Stellar population synthesis codes are used to determine the age, metallicity, velocity, and velocity dispersion of stars within each region. The ICL spectra are best modelled with populations that are younger and less metal-rich than those of the BCG cores. The average BCG core age of the sample is $\rm 13.3\pm 2.8\, Gyr$ and the average metallicity is $\rm [Fe/H] = 0.30\pm 0.09$, whereas for the ICL the average age is $\rm 9.2\pm 3.5\, Gyr$ and the average metallicity is $\rm [Fe/H] = 0.18\pm 0.16$. The velocity dispersion profile is seen to be rising or flat in most of the sample (17/23), and those with rising values reach the value of the host cluster’s velocity dispersion in several cases. The most extended BCGs are closest to the peak of the cluster’s X-ray luminosity. The results are consistent with the idea that the BCG cores and inner regions formed quickly and long ago, with the outer regions and ICL forming more recently, and continuing to assemble through minor merging. Any recent star formation in the BCGs is a minor component, and is associated with the cluster cool core status.

Kinematic signatures of cluster formation from cool collapse in the Lagoon Nebula cluster NGC 6530

ABSTRACT We examine the mass dependence of the velocity dispersion of stars in the young cluster NGC 6530 to better understand how it formed. Using a large sample of members we find that the proper motion velocity dispersion increases with stellar mass. While this trend is the opposite of that predicted if the cluster were developing energy equipartition, it is in agreement with recent N-body simulations that find such a trend develops because of the Spitzer instability. In these simulations the massive stars sink to the centre of the cluster and form a self-gravitating system with a higher velocity dispersion. If the cluster has formed by the cool collapse of an initially substructured distribution, then this occurs within 1–2 Myr, in agreement with our observations of NGC 6530. We therefore conclude that NGC 6530 formed from much more extended initial conditions and has since collapsed to form the cluster we see now. This cluster formation model is inconsistent with the idea that all stars form in dense, compact clusters and provides the first dynamical evidence that star clusters can form by hierarchical mergers between subclusters.

Gravitational instability and star formation in NGC 628

The gas-stars instability criterion for infinitesimally thin disc was applied to the galaxy NGC 628. Instead of using the azimuthally averaged profiles of data the maps of the gas surface densities (THINGS, HERACLES), of the velocity dispersions of stars (VENGA) and gas (THINGS) and of the surface brightness of the galaxy (S$^4$G) were analyzed. All these maps were collected for the same region with a noticeable star formation rate and were superimposed on each other. Using the data on the rotation curve values of $Q_\mathrm{eff}$ were calculated for each pixel in the image. The areas within the contours $Q_\mathrm{eff}<3$ were compared with the ongoing star formation regions ($\Sigma_\mathrm{SFR}> 0.007 \, M_{\odot}$yr$^{-1}$kpc$^{-2}$) and showed a good coincidence between them. The Romeo-Falstad disc instability diagnostics taking into account the thickness of the stellar and gas layers does not change the result. If the one-fluid instability criterion is used, the coincidence is worse. The analysis was carried out for the area $r < $ 0.5$r_{25} $. Leroy et al. (2008) using azimuthally averaged data obtained $Q_\mathrm{eff} \approx 3-4$ for this area of the disc, which makes it stable against non-axisymmetric perturbations and gas dissipation, and does not predict the location of star forming regions. Since in the galaxies the distribution of hydrogen and the regions of star formation is often patchy, the relationship between gravitational instability and star formation should be sought using data maps rather than azimuthally averaged data.