We use the kinematics of ∼200 000 giant stars that lie within ∼1.5 kpc of the plane to measure the vertical profile of mass density near the Sun. We find that the dark mass contained within the ...isodensity surface of the dark halo that passes through the Sun ((6 ± 0.9) × 1010 M⊙), and the surface density within 0.9 kpc of the plane ((69 ± 10) M⊙ pc−2) are almost independent of the (oblate) halo's axis ratio q. If the halo is spherical, 46 per cent of the radial force on the Sun is provided by baryons, and only 4.3 per cent of the Galaxy's mass is baryonic. If the halo is flattened, the baryons contribute even less strongly to the local radial force and to the Galaxy's mass. The dark matter density at the location of the Sun is 0.0126 q
−0.89 M⊙ pc−3 = 0.48 q
−0.89 GeV cm−3. When combined with other literature results we find hints for a mildly oblate dark halo with q ≃ 0.8. Our value for the dark mass within the solar radius is larger than that predicted by cosmological dark-matter-only simulations but in good agreement with simulations once the effects of baryonic infall are taken into account. Our mass models consist of three double-exponential discs, an oblate bulge and a Navarro–Frenk–White dark matter halo, and we model the dynamics of the RAVE (RAdial Velocity Experiment) stars in the corresponding gravitational fields by finding distribution functions f J
that depend on three action integrals. Statistical errors are completely swamped by systematic uncertainties, the most important of which are the distance to the stars in the photometric and spectroscopic samples and the solar distance to the Galactic Centre. Systematics other than the flattening of the dark halo yield overall uncertainties ∼15 per cent.
New distances to RAVE stars Binney, J; Burnett, B; Kordopatis, G ...
Monthly Notices of the Royal Astronomical Society,
2014, Letnik:
437, Številka:
1
Journal Article
Recenzirano
Odprti dostop
Probability density functions (pdfs) are determined from new stellar parameters for the distance moduli of stars for which the RAdial Velocity Experiment (RAVE) has obtained spectra with S/N ≥ 10. ...Single-Gaussian fits to the pdf in distance modulus suffice for roughly half the stars, with most of the other half having satisfactory two-Gaussian representations. As expected, early-type stars rarely require more than one Gaussian. The expectation value of distance is larger than the distance implied by the expectation of distance modulus; the latter is itself larger than the distance implied by the expectation value of the parallax. Our parallaxes of Hipparcos stars agree well with the values measured by Hipparcos, so the expectation of parallax is the most reliable distance indicator. The latter are improved by taking extinction into account. The effective temperature-absolute magnitude diagram of our stars is significantly improved when these pdfs are used to make the diagram. We use the method of kinematic corrections devised by Schönrich, Binney and Asplund to check for systematic errors for general stars and confirm that the most reliable distance indicator is the expectation of parallax. For cool dwarfs and low-gravity giants, 〈ϖ〉 tends to be larger than the true distance by up to 30 per cent. The most satisfactory distances are for dwarfs hotter than 5500 K. We compare our distances to stars in 13 open clusters with cluster distances from the literature and find excellent agreement for the dwarfs and indications that we are overestimating distances to giants, especially in young clusters.
We present the stellar atmospheric parameters (effective temperature, surface gravity, overall metallicity), radial velocities, individual abundances, and distances determined for 425,561 stars, ...which constitute the fourth public data release of the RAdial Velocity Experiment (RAVE). The stellar atmospheric parameters are computed using a new pipeline, based on the algorithms of MATISSE and DEGAS. The spectral degeneracies and the Two Micron All Sky Survey photometric information are now better taken into consideration, improving the parameter determination compared to the previous RAVE data releases. The individual abundances for six elements (magnesium, aluminum, silicon, titanium, iron, and nickel) are also given, based on a special-purpose pipeline that is also improved compared to that available for the RAVE DR3 and Chemical DR1 data releases. Together with photometric information and proper motions, these data can be retrieved from the RAVE collaboration Web site and the Vizier database.
This paper presents the current state of the global gyrokinetic code Orb5 as an update of the previous reference (Jolliet et al., 2007). The Orb5 code solves the electromagnetic Vlasov–Maxwell system ...of equations using a PIC scheme and also includes collisions and strong flows. The code assumes multiple gyrokinetic ion species at all wavelengths for the polarization density and drift-kinetic electrons. Variants of the physical model can be selected for electrons such as assuming an adiabatic response or a “hybrid” model in which passing electrons are assumed adiabatic and trapped electrons are drift-kinetic. A Fourier filter as well as various control variates and noise reduction techniques enable simulations with good signal-to-noise ratios at a limited numerical cost. They are completed with different momentum and zonal flow-conserving heat sources allowing for temperature-gradient and flux-driven simulations. The code, which runs on both CPUs and GPUs, is well benchmarked against other similar codes and analytical predictions, and shows good scalability up to thousands of nodes.
Mass models of the Milky Way McMillan, Paul J.
Monthly Notices of the Royal Astronomical Society,
July 2011, Letnik:
414, Številka:
3
Journal Article
Recenzirano
Odprti dostop
We present a simple method for fitting parametrized mass models of the Milky Way to observational constraints. We take a Bayesian approach which allows us to take into account input from photometric ...and kinematic data, and expectations from theoretical modelling. This provides us with a best-fitting model, which is a suitable starting point for dynamical modelling. We also determine a probability density function on the properties of the model, which demonstrates that the mass distribution of the Galaxy remains very uncertain. For our choices of parametrization and constraints, we find disc scalelengths of 3.00 ± 0.22 and 3.29 ± 0.56 kpc for the thin and thick discs, respectively, a solar radius of 8.29 ± 0.16 kpc and a circular speed at the Sun of 239 ± 5 km s−1, a total stellar mass of 6.43 ± 0.63 × 1010 M⊙, a virial mass of 1.26 ± 0.24 × 1012 M⊙ and a local dark matter density of 0.40 ± 0.04 GeV cm−3. We find some correlations between the best-fitting parameters of our models (for example, between the disc scalelengths and the solar radius), which we discuss. The chosen disc scaleheights are shown to have little effect on the key properties of the model.
There is mounting evidence that associates brain injury and offending behaviour, and there is a need to understand the epidemiology of head injury in prisoners in order to plan interventions to ...reduce associated disability and risk of reoffending. This is the first study to determine the lifetime prevalence of hospitalised head injury (HHI) in a national population of current prison inmates. In addition characteristics of prisoners with HHI and were compared to prisoners without HHI to discover whether those with HI differed demographically.
Whole life hospital records of everyone aged 35 years or younger and resident in a prison in Scotland on a census date in 2015 were electronically linked via their unique NHS identifier and checked for ICD-9 and 10 codes for head injury. Using a case-control design, these data were compared with a sample from the general population matched 3:1 for age, gender and area-based social deprivation. Comparison of demographic variables was made between prisoners with and without HHI.
HHI was found in 24.7% (1,080/4,374) of prisoners and was significantly more prevalent than found in the matched general population sample (18.2%; 2394/13122; OR 2.10; 95%CI 1.87, 2.16). The prevalence of HHI in prisoners and controls was similar with the exception of a higher risk of HHI in prisoners in lower deprivation quintiles. Having three or more HHI was more common in prisoners (OR 3.04; 95%CI 2.33, 3.97) as were HHI with ICD codes for intracranial injuries (OR 1.81; 95% CI 1.54, 2.11), suggesting that more severe HHI is more prevalent in prisoners than the general population. The distributions within demographic variables and the characteristics of HHI admissions in prisoners with and without a history of HHI were similar.
Prisoners in Scotland aged 35 years or younger have a higher lifetime prevalence of HHI than the general population and are more likely to have had repeated HI or intracranial injuries. Further work is required to elucidate the correspondence between self-report of HI and hospitalised records and to ascertain persisting effects of HI in prisoners and the need for services to reduce associated disability and risk of reoffending.
Context.
With data releases from the astrometric space mission
Gaia
, exploration of the structure of the Milky Way is now possible in unprecedented detail, and has unveiled many previously unknown ...structures in the Galactic disc and halo. One such feature is the
Gaia
phase spiral where the stars in the Galactic disc form a spiral density pattern in the
Z
−
V
Z
plane. Many questions regarding the phase spiral remain, particularly how its amplitude and rotation change with position in the Galaxy.
Aims.
We aim to characterize the shape, rotation, amplitude, and metallicity of the phase spiral in the outer disc of the Milky Way. This will allow us to better understand which physical processes caused the phase spiral and may provide further clues as to the Milky Way’s past and the events that contributed to its current state.
Methods.
We use
Gaia
data release 3 (DR3) to get full position and velocity data on approximately 31.5 million stars, and metallicity for a subset of them. We then compute the angular momenta of the stars and develop a model to characterise the phase spiral in terms of amplitude and rotation at different locations in the disc.
Results.
We find that the rotation angle of the phase spiral changes with Galactic azimuth and galactocentric radius, making the phase spiral appear to rotate about 3° per degree in Galactic azimuth. Furthermore, we find that the phase spiral in the 2200 − 2400 kpc km s
−1
range of angular momentum is particularly strong compared to the phase spiral that can be observed in the solar neighbourhood. The metallicity of the phase spiral appears to match that of the field stars of the Milky Way disc.
Conclusions.
We created a new model capable of fitting several key parameters of the
Gaia
phase spiral. We have been able to determine the rotation rate of the phase spiral to be about 3° per degree in Galactic azimuth. We find a maximum in the amplitude of the phase spiral at
L
Z
≈ 2300 km kpc s
−1
, which makes the phase spiral clearly visible. These results provide insights into the physical processes that led to the formation of the phase spiral and contribute to our understanding of the Milky Way’s past and present state.
Context.
The physical processes driving the formation of Galactic spiral arms are still under debate. Studies using open clusters favour the description of the Milky Way spiral arms as long-lived ...structures following the classical density wave theory. Current studies comparing the
Gaia
DR2 field stars kinematic information of the solar neighbourhood to simulations, find a better agreement with short-lived arms with a transient behaviour.
Aims.
Our aim is to provide an observational, data-driven view of the Milky Way spiral structure and its dynamics using open clusters as the main tracers, and to contrast it with simulation-based approaches. We used the most complete catalogue of Milky Way open clusters, with astrometric
Gaia
EDR3 updated parameters, estimated astrophysical information, and radial velocities, to revisit the nature of the spiral pattern of the Galaxy.
Methods.
We used a Gaussian mixture model to detect overdensities of open clusters younger than 30 Myr that correspond to the Perseus, Local, Sagittarius, and Scutum spiral Arms, respectively. We used the birthplaces of the open cluster population younger than 80 Myr to trace the evolution of the different spiral arms and compute their pattern speed. We analysed the age-distribution of the open clusters across the spiral arms to explore the differences in the rotational velocity of stars and spiral arms.
Results.
We are able to increase the range in Galactic azimuth where present-day spiral arms are described, better estimating its parameters by adding 264 young open clusters to the 84 high-mass star-forming regions used so far, thus increasing the number of tracers by 314%. We used the evolution of the open clusters from their birth positions to find that spiral arms nearly co-rotate with field stars at any given radius, discarding a common spiral pattern speed for the spiral arms explored.
Conclusions.
The derivation of different spiral pattern speeds for the different spiral arms disfavours classical density waves as the main drivers for the formation of the Milky Way spiral structure, and it is in better agreement with simulation-based approaches that tend to favour transient spirals. The increase in the number of known open clusters, as well as in their derived properties, allows us to use them as effective spiral structure tracers and homogenise the view from open clusters and field stars on the nature of the Galactic spiral arms.
Using the RAdial Velocity Experiment fourth data release (RAVE DR4), and a new metallicity calibration that will be also taken into account in the future RAVE DR5, we investigate the existence and ...the properties of supersolar metallicity stars (M/H ... +0.1 dex) in the sample, and in particular in the solar neighbourhood. We find that RAVE is rich in supersolar metallicity stars, and that the local metallicity distribution function declines remarkably slowly up to +0.4 dex. Our results show that the kinematics and height distributions of the supersolar metallicity stars are identical to those of the M/H ... 0 thin-disc giants that we presume were locally manufactured. The eccentricities of the supersolar metallicity stars indicate that half of them are on a roughly circular orbit (e = 0.15), so under the assumption that the metallicity of the interstellar medium at a given radius never decreases with time, they must have increased their angular momenta by scattering at corotation resonances of spiral arms from regions far inside the solar annulus. The likelihood that a star will migrate radially does not seem to decrease significantly with increasing amplitude of vertical oscillations within range of oscillation amplitudes encountered in the disc. (ProQuest: ... denotes formulae/symbols omitted.)
Kinematics with Gaia DR2: the force of a dwarf Carrillo, I; Minchev, I; Steinmetz, M ...
Monthly Notices of the Royal Astronomical Society,
11/2019, Letnik:
490, Številka:
1
Journal Article
Recenzirano
Odprti dostop
ABSTRACT
We use Gaia DR2 astrometric and line-of-sight velocity information combined with two sets of distances obtained with a Bayesian inference method to study the 3D velocity distribution in the ...Milky Way disc. We search for variations in all Galactocentric cylindrical velocity components (Vϕ, VR, and Vz) with Galactic radius, azimuth, and distance from the disc mid-plane. We confirm recent work showing that bulk vertical motions in the R–z plane are consistent with a combination of breathing and bending modes. In the x–y plane, we show that, although the amplitudes change, the structure produced by these modes is mostly invariant as a function of distance from the plane. Comparing to two different Galactic disc models, we demonstrate that the observed patterns can drastically change in short time intervals, showing the complexity of understanding the origin of vertical perturbations. A strong radial VR gradient was identified in the inner disc, transitioning smoothly from 16 km s−1 kpc−1 at an azimuth of 30° < ϕ < 45° ahead of the Sun-Galactic centre line to −16 km s−1 kpc−1 at an azimuth of −45° < ϕ < −30° lagging the solar azimuth. We use a simulation with no significant recent mergers to show that exactly the opposite trend is expected from a barred potential, but overestimated distances can flip this trend to match the data. Alternatively, using an N-body simulation of the Sagittarius dwarf–Milky Way interaction, we demonstrate that a major recent perturbation is necessary to reproduce the observations. Such an impact may have strongly perturbed the existing bar or even triggered its formation in the last 1–2 Gyr.