-
Comment on 'a slightly oblate dark matter halo revealed by a retrograde precessing Galactic disk warp' by Huang et al
Authors:
Walter Dehnen,
Ralph Schönrich,
Ronald Drimmel,
Žofia Chrobáková,
Eloisa Poggio,
Marcin Semczuk
Abstract:
Huang et al. (2024) measured the derivative of the phase $φ_{\mathrm{w}}$ of the Galactic warp traced by classical Cepheids with respect to their age $τ$ and interpreted it as the warp precession rate $ω\equiv\mathrm{d}φ_{\mathrm{w}}/\mathrm{d}t=-\mathrm{d}φ_{\mathrm{w}}/\mathrm{d}τ$. This interpretation is unfounded: young stars follow trajectories close to those of their parental gas and trace t…
▽ More
Huang et al. (2024) measured the derivative of the phase $φ_{\mathrm{w}}$ of the Galactic warp traced by classical Cepheids with respect to their age $τ$ and interpreted it as the warp precession rate $ω\equiv\mathrm{d}φ_{\mathrm{w}}/\mathrm{d}t=-\mathrm{d}φ_{\mathrm{w}}/\mathrm{d}τ$. This interpretation is unfounded: young stars follow trajectories close to those of their parental gas and trace the instantaneous gas warp, not its shape at their time of birth: $φ_{\mathrm{w}}$ should hardly depend on Cepheid age. We show that the measured $\mathrm{d}φ_{\mathrm{w}}/\mathrm{d}τ>0$ is consistent with an omitted-variable bias from neglecting the natural twist $\mathrm{d}φ_{\mathrm{w}}/\mathrm{d}R$ of the warp and the $R$-$τ$ correlation for Cepheids (originating from the Galactic metallicity gradient and the Cepheid metallicity-age correlation).
△ Less
Submitted 8 July, 2024;
originally announced July 2024.
-
Advancing Galaxy Analysis: AI-Powered Detection and Segmentation of Edge-On Galaxies
Authors:
Z. Chrobakova,
V. Kresnakova,
R. Nagy,
J. Gazdova,
P. Butka
Abstract:
Edge-on galaxies have many important applications in galactic astrophysics, but they can be difficult to identify in vast amounts of astronomical data. To facilitate the search for them, we have developed a deep learning algorithm designed to identify and extract edge-on galaxies from astronomical images. We utilised a sample of edge-on spiral galaxies from the Galaxy Zoo database, retrieving the…
▽ More
Edge-on galaxies have many important applications in galactic astrophysics, but they can be difficult to identify in vast amounts of astronomical data. To facilitate the search for them, we have developed a deep learning algorithm designed to identify and extract edge-on galaxies from astronomical images. We utilised a sample of edge-on spiral galaxies from the Galaxy Zoo database, retrieving the corresponding images from the Sloan Digital Sky Survey (SDSS). Our dataset comprises approximately $16,000$ galaxies, which we used to train the YOLOv5 algorithm for detection purposes. To isolate galaxies from their backgrounds, we trained the SCSS-Net neural network to generate segmentation masks. As a result, our algorithm detected $8,000$ edge-on galaxies, for which we compiled a catalogue including their parameters obtained from the SDSS database. We describe the basic properties of our sample, finding that most galaxies have redshifts $0.02<z<0.10$, have low values of $b/a$ and are mostly red, which is expected from edge-on galaxies and is consistent with our training sample, as well as other literature. The cutouts of the detected galaxies can be used for future studies and the algorithm can be applied to data from future surveys as well.
△ Less
Submitted 21 June, 2024;
originally announced June 2024.
-
Mass models of the Milky Way and estimation of its mass from the GAIA DR3 data-set
Authors:
Francesco Sylos Labini,
Zofia Chrobakova,
Roberto Capuzzo-Dolcetta,
Martin Lopez-Corredoira
Abstract:
We use data from the Gaia DR3 dataset to estimate the mass of the Milky Way (MW) by analyzing the rotation curve in the range of distances 5 kpc to 28 kpc. We consider three mass models: the first model adds a spherical dark matter (DM) halo, following the Navarro-Frenk-White (NFW) profile, to the known stellar components. The second model assumes that DM is confined to the Galactic disk, followin…
▽ More
We use data from the Gaia DR3 dataset to estimate the mass of the Milky Way (MW) by analyzing the rotation curve in the range of distances 5 kpc to 28 kpc. We consider three mass models: the first model adds a spherical dark matter (DM) halo, following the Navarro-Frenk-White (NFW) profile, to the known stellar components. The second model assumes that DM is confined to the Galactic disk, following the idea that the observed density of gas in the Galaxy is related to the presence of more massive DM disk (DMD), similar to the observed correlation between DM and gas in other galaxies. The third model only uses the known stellar mass components and is based on the Modified Newton Dynamics (MOND) theory. Our results indicate that the DMD model is comparable in accuracy to the NFW and MOND models and fits the data better at large radii where the rotation curve declines but has the largest errors. For the NFW model we obtain a virial mass $M_{vir}= (6.5 \pm 0.3) \times 10^{11} \; M_\odot$ with concentration parameter $c=14.5$, that is lower than what is typically reported. In the DMD case we find that the MW mass is $M_d = (1.6 \pm 0.5) \times 10^{11} \; M_\odot$ with a disk's characteristic radius of $R_d=17$ kpc.
△ Less
Submitted 2 February, 2023;
originally announced February 2023.
-
Evidence for Populations-dependent vertical motions and the Long-lived Non-Steady Lopsided Milky Way Warp
Authors:
Xiang Li,
Hai-Feng Wang,
Yang-Ping Luo,
Martín López-Corredoira,
Yuan-Sen Ting,
Žofia Chrobáková
Abstract:
We present the Galactic disk vertical velocity analysis using OB type stars (OB), Red Clump stars (RC), and Main-Sequence-Turn-Off stars (MSTO) with different average age populations crossed matched with LAMOST DR5 and Gaia DR3. We reveal the vertical velocities of the three populations varies clearly with the Galactocentric distance ($R$) and the younger stellar population has stronger increasing…
▽ More
We present the Galactic disk vertical velocity analysis using OB type stars (OB), Red Clump stars (RC), and Main-Sequence-Turn-Off stars (MSTO) with different average age populations crossed matched with LAMOST DR5 and Gaia DR3. We reveal the vertical velocities of the three populations varies clearly with the Galactocentric distance ($R$) and the younger stellar population has stronger increasing trend in general. The bending and breathing modes indicated by the vertical motions are dependent on the populations and they are varying with spatial locations. These vertical motions may be due to the Galactic warp, or minor mergers, or non-equilibrium of the disk. Assuming the warp is the dominant component, we find that the warp amplitude ($γ$, $Z_ω$) for OB (younger population) is larger than that for RC (medium population) and the later one is also larger than that for MSTO (older population), which is in agreement with other independent analyses of stellar density distribution, and supports the warp is long-lived, non-steady structure and has time evolution. This conclusion is robust whether or not the line-of-nodes $φ_w$ is fixed or as a free parameter (with $φ_w$ is around 3$-$8.5$^{\circ}$ as best fit). Furthermore, we find that warp is lopsided with asymmetries along azimuthal angle ($φ$).
△ Less
Submitted 12 December, 2022;
originally announced December 2022.
-
Mapping the Milky Way Disk with Gaia DR3: 3D extended kinematic maps and rotation curve to $\approx 30$ kpc
Authors:
Hai-Feng Wang,
Žofia Chrobáková,
Martín López-Corredoira,
Francesco Sylos Labini
Abstract:
We apply a statistical deconvolution of the parallax errors based on Lucy's inversion method (LIM) to the Gaia-DR3 sources to measure their three dimensional velocity components in the range of Galactocentric distances $R$ between 8 kpc and 30 kpc with their corresponding errors and root mean square values. We find results that are consistent with those obtained by applying LIM to the Gaia-DR2 sou…
▽ More
We apply a statistical deconvolution of the parallax errors based on Lucy's inversion method (LIM) to the Gaia-DR3 sources to measure their three dimensional velocity components in the range of Galactocentric distances $R$ between 8 kpc and 30 kpc with their corresponding errors and root mean square values. We find results that are consistent with those obtained by applying LIM to the Gaia-DR2 sources, and we conclude that the method gives convergent and more accurate results by improving the statistics of the data-set and lowering observational errors. The kinematic maps reconstructed with LIM up to $R \approx 30$ kpc show that the Milky Way is characterized by asymmetrical motions with significant gradients in all velocity components. Furthermore, we determine the Galaxy rotation curve $V_C(R)$ up to $\approx 27.5$ kpc with the cylindrical Jeans equation assuming an axisymmetric gravitational potential. We find that $V_C(R)$ is significantly declining up to the largest radius investigated. Finally, we also measure $V_C(R)$ at different vertical heights, showing that, for $R <15$ kpc, there is a marked dependence on $Z$, whereas at larger $R$ the dependence on $Z$ is negligible.
△ Less
Submitted 7 December, 2022; v1 submitted 10 November, 2022;
originally announced November 2022.
-
Virial theorem in clusters of galaxies with MOND
Authors:
M. Lopez-Corredoira,
J. E. Betancort-Rijo,
R. Scarpa,
Z. Chrobakova
Abstract:
A specific modification of Newtonian dynamics known as MOND has been shown to reproduce the dynamics of most astrophysical systems at different scales without invoking non-baryonic dark matter (DM). There is, however, a long-standing unsolved problem when MOND is applied to rich clusters of galaxies in the form of a deficit (by a factor around two) of predicted dynamical mass derived from the viri…
▽ More
A specific modification of Newtonian dynamics known as MOND has been shown to reproduce the dynamics of most astrophysical systems at different scales without invoking non-baryonic dark matter (DM). There is, however, a long-standing unsolved problem when MOND is applied to rich clusters of galaxies in the form of a deficit (by a factor around two) of predicted dynamical mass derived from the virial theorem with respect to observations. In this article we approach the virial theorem using the velocity dispersion of cluster members along the line of sight rather than using the cluster temperature from X-ray data and hydrostatic equilibrium. Analytical calculations of the virial theorem in clusters for Newtonian gravity+DM and MOND are developed, applying pressure (surface) corrections for non-closed systems. Recent calibrations of DM profiles, baryonic ratio and baryonic ($β$ model or others) profiles are used, while allowing free parameters to range within the observational constraints. It is shown that solutions exist for MOND in clusters that give similar results to Newton+DM -- particularly in the case of an isothermal $β$ model for $β=0.55-0.70$ and core radii $r_c$ between 0.1 and 0.3 times $r_{500}$ (in agreement with the known data). The disagreements found in previous studies seem to be due to the lack of pressure corrections (based on inappropriate hydrostatic equilibrium assumptions) and/or inappropriate parameters for the baryonic matter profiles.
△ Less
Submitted 25 October, 2022;
originally announced October 2022.
-
Analysing the structure of the bulge with Mira variables
Authors:
Ž. Chrobáková,
M. López-Corredoira,
F. Garzón
Abstract:
The Galactic bulge at latitude $4<|b|$(deg.)$<10$ was claimed to show an X-shape, which means that stellar density distributions along the line of sight have a double peak. However, this double peak is only observed with the red-clump population, and doubt has been cast on its use as a perfect standard candle.\ As such, a boxy bulge without an X-shape is not discarded. We aim to constrain the shap…
▽ More
The Galactic bulge at latitude $4<|b|$(deg.)$<10$ was claimed to show an X-shape, which means that stellar density distributions along the line of sight have a double peak. However, this double peak is only observed with the red-clump population, and doubt has been cast on its use as a perfect standard candle.\ As such, a boxy bulge without an X-shape is not discarded. We aim to constrain the shape of the bulge making use of a different population: Mira variables from the new Optical Gravitational Lensing Experiment data release, OGLE-IV, with an average age of 9 Gyr. We analysed an area of the bulge far from the plane, where we fitted the density of the Miras with boxy bulge and X-shaped bulge models and calculated the probability of each model. We find that the probability of a boxy bulge fitting the data is $p=0.19$, whereas the probability for the X-shaped bulge is only $p=2.85 \cdot 10^{-6}$ (equivalent to a tension of the model with the data of a 4.7$σ$ level). Therefore, the boxy bulge model seems to be more appropriate for describing the Galactic bulge, although we cannot exclude any model with complete certainty.
△ Less
Submitted 7 October, 2022;
originally announced October 2022.
-
Warp and flare of the Galactic disc revealed with supergiants by Gaia EDR3
Authors:
Ž. Chrobáková,
R. Nagy,
M. López-Corredoira
Abstract:
The outer Galactic disc contains some features such as the warp and flare, whose origin is still debated. The Gaia data provide an excellent opportunity to probe the Galactic disc at large distances and study these features. We derive the density distributions of the average (old) whole population and the supergiants (representative of a young population), and we use them to constrain their warp a…
▽ More
The outer Galactic disc contains some features such as the warp and flare, whose origin is still debated. The Gaia data provide an excellent opportunity to probe the Galactic disc at large distances and study these features. We derive the density distributions of the average (old) whole population and the supergiants (representative of a young population), and we use them to constrain their warp and flare. By comparing the results, we study how the properties of these phenomena depend on the studied population. We used Lucy's deconvolution method to recover corrected star counts as a function of distance, from which we derive the density distribution. We find that supergiants have an asymmetric warp, reaching a maximum amplitude of $z_w=0.658$ kpc and minimum amplitude of $z_w=-0.717$ kpc at a distance of $R=[19.5,20]$ kpc, which is almost twice as high as the amplitude of the whole population of the disc. We find a significant flare of the whole population, especially in the thick disc. The scale height increases from $h_{z,thick}\approx 0.8$ kpc and $h_{z,thin}\approx 0.3$ kpc in the solar neighbourhood, to $h_{z,thick} \approx 3$ kpc and $h_{z,thin}\approx 0.7$ kpc in the remote regions of the Milky Way ($R\approx 18$ kpc). The supergiants' population has only a small flare.
△ Less
Submitted 22 June, 2022; v1 submitted 16 June, 2022;
originally announced June 2022.
-
A Case against a Significant Detection of Precession in the Galactic Warp
Authors:
Žofia Chrobáková,
Martín López-Corredoira
Abstract:
Recent studies of warp kinematics using Gaia DR2 data have produced detections of warp precession for the first time, which greatly exceeds theoretical predictions of models. However, this detection assumes a warp model derived for a young population (few tens of megayears) to fit velocities of an average older stellar population of the thin disk (several gigayears) in Gaia-DR2 observations, which…
▽ More
Recent studies of warp kinematics using Gaia DR2 data have produced detections of warp precession for the first time, which greatly exceeds theoretical predictions of models. However, this detection assumes a warp model derived for a young population (few tens of megayears) to fit velocities of an average older stellar population of the thin disk (several gigayears) in Gaia-DR2 observations, which may lead to unaccounted systematic errors. Here, we recalculate the warp precession with the same approach and Gaia DR2 kinematic data, but using different warp parameters based on the fit of star counts of the Gaia DR2 sample, which has a much lower warp amplitude than the young population. When we take into account this variation of the warp amplitude with the age of the population, we find that there is no need for precession. We find the value of warp precession $β= 4^{+6}_{-4}$ km s$^{-1}$ kpc$^{-1}$, which does not exclude nonprecessing warp.
△ Less
Submitted 10 May, 2021;
originally announced May 2021.
-
Gaia-DR2 extended kinematical maps. Part III: Rotation curves analysis, dark matter, and Modified Newtonian dynamics tests
Authors:
Ž. Chrobáková,
M. López-Corredoira,
F. Sylos Labini,
H. -F. Wang,
R. Nagy
Abstract:
Recent statistical deconvolution methods have produced extended kinematical maps in a range of heliocentric distances that are a factor of two to three larger than those analysed in the Gaia Collaboration based on the same data. In this paper, we use such maps to derive the rotation curve both in the Galactic plane and in off-plane regions and to analyse the density distribution. By assuming stati…
▽ More
Recent statistical deconvolution methods have produced extended kinematical maps in a range of heliocentric distances that are a factor of two to three larger than those analysed in the Gaia Collaboration based on the same data. In this paper, we use such maps to derive the rotation curve both in the Galactic plane and in off-plane regions and to analyse the density distribution. By assuming stationary equilibrium and axisymmetry, we used the Jeans equation to derive the rotation curve. Then we fit it with density models that include both dark matter and predictions of the MOND (Modified Newtonian dynamics) theory. Since the Milky Way exhibits deviations from axisymmetry and equilibrium, we also considered corrections to the Jeans equation. To compute such corrections, we ran N-body experiments of mock disk galaxies where the departure from equilibrium becomes larger as a function of the distance from the centre. The rotation curve in the outer disk of the Milky Way that is constructed with the Jeans equation exhibits very low dependence on $R$ and $z$ and it is well-fitted both by dark matter halo and MOND models. The application of the Jeans equation for deriving the rotation curve, in the case of the systems that deviate from equilibrium and axisymmetry, introduces systematic errors that grow as a function of the amplitude of the average radial velocity. In the case of the Milky Way, we can observe that the amplitude of the radial velocity reaches $\sim 10\%$ that of the azimuthal one at $R\approx 20$ kpc. Based on this condition, using the rotation curve obtained from the Jeans equation to calculate the mass may overestimate its measurement.
△ Less
Submitted 29 July, 2020;
originally announced July 2020.
-
Mapping the Galactic disk with the LAMOST and Gaia Red clump sample. VI: An evidence for the long lived non-steady warp of non-gravitational scenarios
Authors:
H. -F. Wang,
M. López-Corredoira,
Y. Huang,
J. Chang,
H. -W. Zhang,
J. L. Carlin,
X. -D. Chen,
Ž. Chrobáková,
B. -Q. Chen
Abstract:
By combining LAMOST DR4 and Gaia DR2 common red clump stars with age and proper motion, we analyze the amplitude evolution of the stellar warp independently of any assumption with a simple model. The greatest height of the warp disk increases with Galactocentric distance in different populations and it is dependent on the age: the younger stellar populations exhibit stronger warp features than the…
▽ More
By combining LAMOST DR4 and Gaia DR2 common red clump stars with age and proper motion, we analyze the amplitude evolution of the stellar warp independently of any assumption with a simple model. The greatest height of the warp disk increases with Galactocentric distance in different populations and it is dependent on the age: the younger stellar populations exhibit stronger warp features than the old ones, accompanied with the warp amplitude $γ({\rm age})$ decreasing with age and its first derivative $\dotγ ({\rm age})$ is different from zero. The azimuth of line of nodes $φ_w$ is stable at $-$5 degree without clear time evolution, which perfectly confirms some of previous works. All these self-consistent evidences support that our Galactic warp should most likely be a long lived, but non-steady structure and not a transient one, which is supporting the warp is originated from gas infall onto the disk or other hypotheses that suppose that the warp mainly affects to the gas, and consequently younger populations tracing the gas are stronger than the older ones. In other words, the Galactic warp is induced by the non-gravitational interaction over the disk models.
△ Less
Submitted 14 May, 2020; v1 submitted 14 April, 2020;
originally announced April 2020.
-
Structure of the outer Galactic disc with Gaia-DR2
Authors:
Z. Chrobakova,
R. Nagy,
M. Lopez-Corredoira
Abstract:
AIMS. We calculate the stellar density using star counts obtained from Gaia DR2 up to a Galactocentric distance R=20 kpc with a deconvolution technique for the parallax errors. Then we analyse the density in order to study the structure of the outer Galactic disc, mainly the warp.
METHODS. In order to carry out the deconvolution, we used the Lucy inversion technique for recovering the corrected…
▽ More
AIMS. We calculate the stellar density using star counts obtained from Gaia DR2 up to a Galactocentric distance R=20 kpc with a deconvolution technique for the parallax errors. Then we analyse the density in order to study the structure of the outer Galactic disc, mainly the warp.
METHODS. In order to carry out the deconvolution, we used the Lucy inversion technique for recovering the corrected star counts. We also used the Gaia luminosity function of stars with $M_G<10$ to extract the stellar density from the star counts.
RESULTS. The stellar density maps can be fitted by an exponential disc in the radial direction $h_r=2.07\pm0.07$ kpc, with a weak dependence on the azimuth, extended up to 20 kpc without any cut-off. The flare and warp are clearly visible. The best fit of a symmetrical S-shaped warp gives $z_w= z_\odot+(37\pm 4.2(stat.)-0.91(syst.))$ pc $(R/R_\odot )^{2.42\pm 0.76(stat.) + 0.129 (syst.)} sin(φ+9.3\pm 7.37 (stat.) +4.48 (syst.))$ for the whole population. When we analyse the northern and southern warps separately, we obtain an asymmetry of an $\sim25\%$ larger amplitude in the north. This result may be influenced by extinction because the Gaia G band is quite prone to extinction biases. However, we tested the accuracy of the extinction map we used, which shows that the extinction is determined very well in the outer disc. Nevertheless, we recall that we do not know the full extinction error, and neither do we know the systematic error of the map, which may influence the final result.
The analysis was also carried out for very luminous stars alone ($M_G<-2$), which on average represents a younger population. We obtain similar scale-length values, while the maximum amplitude of the warp is $20-30\%$ larger than with the whole population. The north-south asymmetry is maintained.
△ Less
Submitted 7 April, 2020;
originally announced April 2020.
-
Gaia-DR2 extended kinematical maps. Part II: Dynamics in the Galactic disk explaining radial and vertical velocities
Authors:
M. Lopez-Corredoira,
F. Garzon,
H. -F. Wang,
F. Sylos Labini,
R. Nagy,
Z. Chrobakova,
J. Chang,
B. Villarroel
Abstract:
Context: In our Paper I, by using statistical deconvolution methods, extended kinematics maps of Gaia-DR2 data have been produced in a range of heliocentric distances that are a factor of two to three larger than those analyzed previously by the Gaia Collaboration with the same data. It added the range of Galactocentric distances between 13 kpc and 20 kpc to the previous maps.
Aims: Here, we inv…
▽ More
Context: In our Paper I, by using statistical deconvolution methods, extended kinematics maps of Gaia-DR2 data have been produced in a range of heliocentric distances that are a factor of two to three larger than those analyzed previously by the Gaia Collaboration with the same data. It added the range of Galactocentric distances between 13 kpc and 20 kpc to the previous maps.
Aims: Here, we investigate the dynamical effects produced by different mechanisms that can explain the radial and vertical components of these extended kinematic maps, including a decomposition of bending and breathing of the vertical components. This paper as a whole tries to be a compendium of different dynamical mechanisms whose predictions can be compared to the kinematic maps.
Methods: Using analytical methods or simulations, we are able to predict the main dynamical factors and compare them to the predictions of the extended kinematic maps of Gaia-DR2.
Results: The gravitational influence of Galactic components that are different from the disk, such as the long bar or bulge, the spiral arms, or a tidal interaction with Sagittarius dwarf galaxy, may explain some features of the velocity maps, especially in the inner parts of the disk. However, they are not sufficient in explaining the most conspicuous gradients in the outer disk. Vertical motions might be dominated by external perturbations or mergers, although a minor component may be due to a warp whose amplitude evolves with time. Here, we show with two different methods, which analyze the dispersion of velocities, that the mass distribution of the disk is flared. Despite these partial explanations, the main observed features can only be explained in terms of out-of-equilibrium models, which are either due to external perturbers or to the fact that the disk has not had time to reach equilibrium since its formation.
△ Less
Submitted 15 January, 2020;
originally announced January 2020.