Stellar Number Density Research Articles

First spiral arm detection using dynamical mass measurements of the Milky Way disk

We applied the vertical Jeans equation to the Milky Way disk in order to study non-axisymmetric variations in the thin disk surface density. We divided the disk plane into area cells with a 100 pc grid spacing and used four separate subsets of the Gaia DR3 stars, defined by cuts in absolute magnitude, that reach distances up to 3 kpc. The vertical Jeans equation is informed by the stellar number density field and the vertical velocity field; for the former, we used maps produced via Gaussian process regression; for the latter, we used Bayesian neural network radial velocity predictions, which allowed us to utilise the full power of the Gaia DR3 proper motion sample. For the first time, we find evidence of a spiral arm in the form of an over-density in the dynamically measured disk surface density, detected in all four data samples, which agrees very well with the spiral arm as traced by stellar age and chemistry. We fitted a simple spiral arm model to this feature and infer a relative over-density of roughly 20% and a width of roughly 400 pc. We also infer a thin disk surface density scale length of 3.3–4.2 kpc when restricting the analysis to stars within a distance of 2 kpc.

Astrometry in crowded fields towards the Galactic bulge

Context. The astrometry towards the Galactic bulge is hampered by high stellar crowding and patchy extinction. This effect is particularly severe for optical surveys such as the European Space Agency satellite Gala. Aims. In this study, we assess the consistency of proper motion measurements between optical (Gaia DR3) and near-infrared (VIRAC2) catalogues in comparison with proper motions measured with the Hubble Space Telescope (HST) observations in several crowded fields towards the Galactic bulge and in Galactic globular clusters. Methods. Assuming that the proper motion measurements are well characterised, the uncertainty-normalised proper motion differences between pairs of catalogues are expected to follow a normal distribution. A deviation from a normal distribution defines the inflation factor r. By multiplying the proper motion uncertainties with the appropriate inflation factor values, the Gaia (VIRAC2) proper motion measurements are brought into a 1σ agreement with the HST proper motions. Results. The inflation factor (r) depends on stellar surface density. For the brightest stars in our sample (G < 18), the dependence on G-band magnitude is strong, corresponding to the most precise Gaia DR3 proper motions. We used the number of observed Gaia DR3 sources as a proxy for the stellar surface density. Assuming that the HST proper motion measurements are well determined and free from systematic errors, we find that Gaia DR3 proper motion uncertainties are better characterised, having r < 1.5 in fields with a stellar number density with fewer than 200 Gaia DR3 sources per arcmin2, and are underestimated by up to a factor of 4 in fields with stellar densities higher than 300 sources per arcmin2. For the most crowded fields in VIRAC2, the proper motion uncertainties are underestimated by a factor of 1.1 up to 1.5, with a dependence on J-band magnitude. In all fields, the brighter sources have the higher r value. At the faint end (G > 19), the inflation factor is close to 1, meaning that the proper motions already fully agree with the HST measurements within 1σ. Conclusions. In the crowded fields common to both catalogues, VIRAC2 proper motions agree with HST proper motions and do not need an inflation factor for their uncertainties. Because of the depth and completeness of VIRAC2 in these fields, it is an ideal complement to Gaia DR3 for proper motion studies towards the Galactic bulge.

Radial distribution of distant trans-Neptunian objects points to Sun’s formation in a stellar cluster

The Scattered Disk Objects (SDOs) are a population of trans-Neptunian bodies with semimajor axes 50<a≲1000 au and perihelion distances q≳30 au. The detached SDOs with orbits beyond the reach of Neptune (roughly q>35 au) are of special interest here as an important constraint on the early evolution of the outer Solar System. The semimajor axis profile of detached SDOs at 50–500 au, as characterized from the Dark Energy Survey (DES), is radially extended, but previous dynamical models of Neptune’s early migration produce a relatively compact profile. This problem is most likely related to Sun’s birth environment in a stellar cluster. We perform new dynamical simulations that account for cluster effects and show that the orbital distribution of SDOs can be explained if a particularly close stellar encounter occurred early on (e.g., M dwarf with the mass ≃0.2M⊙ approaching the Sun at ≃200 au). For such an encounter to happen with a reasonably high probability the Sun must have formed in a stellar cluster with ηT≳104 Myr pc−3, where η is the stellar number density and T is the Sun’s residence time in the cluster.

S5: Probing the Milky Way and Magellanic Clouds potentials with the 6D map of the Orphan–Chenab stream

ABSTRACTWe present a 6D map of the Orphan–Chenab (OC) stream by combining the data from Southern Stellar Stream Spectroscopic Survey (S5) and Gaia. We reconstruct the proper motion, radial velocity, distance, on-sky track, and stellar density along the stream with spline models. The stream has a total luminosity of MV = −8.2 and metallicity of [Fe/H] = −1.9, similar to classical Milky Way (MW) satellites like Draco. The stream shows drastic changes in its physical width varying from 200 pc to 1 kpc, but a constant line-of-sight velocity dispersion of 5 $\mathrm{km\, s^{-1}}$. Despite the large apparent variation in the stellar number density along the stream, the flow rate of stars along the stream is remarkably constant. We model the 6D stream track by a Lagrange-point stripping method with a flexible MW potential in the presence of a moving extended Large Magellanic Cloud (LMC). This allows us to constrain the mass profile of the MW within the distance range 15.6 &amp;lt; r &amp;lt; 55.5 kpc, with the best measured enclosed mass of $(2.85\pm 0.1)\times 10^{11}\, \mathrm{\, M_\odot }$ within 32.4 kpc. Our stream measurements are highly sensitive to the LMC mass profile with the most precise measurement of its enclosed mass made at 32.8 kpc, $(7.02\pm 0.9)\times 10^{10}\, {\rm M}_\odot$. We also detect that the LMC dark matter halo extends to at least 53 kpc. The fitting of the OC stream allows us to constrain the past LMC trajectory and the degree of dynamical friction it experienced. We demonstrate that the stars in the OC stream show large energy and angular momentum spreads caused by LMC perturbation.

The Fifth Catalogue of Nearby Stars (CNS5)

Context. We present the compilation of the Fifth Catalogue of Nearby Stars (CNS5), based on astrometric and photometric data from Gaia EDR3 and HIPPARCOS and supplemented with parallaxes from ground-based astrometric surveys carried out in the infrared. Aims. The aim of the CNS5 is to provide the most complete sample of objects in the solar neighbourhood. For all known stars and brown dwarfs in the 25 pc sphere around the Sun, basic astrometric and photometric parameters are given. Furthermore, we provide the colour-magnitude diagram (CMD) and various luminosity functions of the stellar content in the solar neighbourhood, and characterise the completeness of the CNS5 catalogue. Methods. We compiled a sample of stars and brown dwarfs that are most likely located within 25 pc of the Sun, taking space-based parallaxes from Gaia EDR3 and HIPPARCOS as well as ground-based parallaxes from Best et al. (2021, AJ, 161, 42), Kirkpatrick et al. (2021, ApJS, 253, 7), and from the CNS4 into account. We developed a set of selection criteria to clean the sample from spurious sources. Furthermore, we show that effects of blending in the Gaia photometry, which mainly affect the faint and red sources in Gaia, can be mitigated to reliably place those objects in a CMD. We also assessed the completeness of the CNS5 using a Kolmogorov-Smirnov test and derive observational optical and mid-infrared (MIR) luminosity functions for the main-sequence stars and white dwarfs (WDs) in the solar neighbourhood. Results. The CNS5 contains 5931 objects, including 5230 stars (4946 main-sequence stars, 20 red giants and 264 white dwarfs) and 701 brown dwarfs. We find that the CNS5 catalogue is statistically complete down to 19.7 mag in the G-band and 11.8 mag in W1-band absolute magnitudes, corresponding to a spectral type of L8. The stellar number density in the solar neighbourhood is (7.99 ± 0.11) × 10−2 stars pc−3, and about 72% of stars in the solar neighbourhood are M dwarfs. Furthermore, we show that the WD sample in CNS5 is statistically complete within 25 pc. The derived number density of WDs is (4.03 ± 0.25) × 10−3 stars pc−3. The ratio between stars and brown dwarfs within 15 pc is 4.6 ± 0.4, whereas within 25 pc it is 7.5 ± 0.3. Thus, we estimate that about one third of brown dwarfs are still missing within 25 pc, particularly those with spectral types later than L8 and distances close to the 25 pc limit.

Mapping Milky Way disk perturbations in stellar number density and vertical velocity using Gaia DR3

We have mapped the number density and mean vertical velocity of the Milky Way’s stellar disk out to roughly two kiloparsecs from the Sun using Gaia Data Release 3 (DR3) and complementary photo-astrometric distance information from StarHorse. For the number counts, we carefully masked spatial regions that are compromised by open clusters, great distances, or dust extinction and used Gaussian processes to arrive at a smooth, non-parametric estimate for the underlying number density field. We find that the number density and velocity fields depart significantly from an axisymmetric and mirror-symmetric model. These departures, which include projections of the Gaia phase-space spiral, signal the presence of local disturbances in the disk. We identify two features that are present in both stellar number density and mean vertical velocity. One of these features appears to be associated with the Local Spiral Arm. It is most prominent at small heights and is largely symmetric across the mid-plane of the disk. The density and velocity field perturbations are phase-shifted by roughly a quarter wavelength, suggesting a breathing mode that is propagating in the direction of Galactic longitude l ∼ 270 deg. The second feature is a gradient in the stellar number density and mean vertical velocity with respect to galactocentric radius. This feature, which extends across the entire region of our analysis, may be associated with the extension of the Galactic warp into the solar neighbourhood in combination with more localised bending waves.

Assessing the physical reality of Milky Way open cluster candidates

ABSTRACT We report results on the analysis of 11 new Milky Way open cluster candidates, recently discovered from the detection of stellar overdensities in the Vector Point diagram, by employing Xtreme deconvolution Gaussian mixture models. We treated these objects as real open clusters and derived their fundamental properties with their associated intrinsic dispersions by exploring the parameter space through the minimization of likelihood functions on the generated synthetic colour–magnitude diagrams. The intrinsic dispersions of the resulting ages turned out to be much larger than those usually obtained for open clusters. Indeed, they resemble the ages and metallicities of composite star field populations. We also traced their stellar number density profiles and mass functions and derived their total masses and Jacobi and tidal radii, which helped us as criteria while assessing their physical nature as real open clusters. Because the 11 candidates show a clear gathering of stars in the proper-motion plane and some hint for similar distances, we concluded that they are possibly sparse groups of stars.

KiDS and Euclid: Cosmological implications of a pseudo angular power spectrum analysis of KiDS-1000 cosmic shear tomography

We present a tomographic weak lensing analysis of the Kilo Degree Survey Data Release 4 (KiDS-1000), using a new pseudo angular power spectrum estimator (pseudo-Cℓ) under development for the ESA Euclid mission. Over 21 million galaxies with shape information are divided into five tomographic redshift bins, ranging from 0.1 to 1.2 in photometric redshift. We measured pseudo-Cℓ using eight bands in the multipole range 76 &lt; ℓ &lt; 1500 for auto- and cross-power spectra between the tomographic bins. A series of tests were carried out to check for systematic contamination from a variety of observational sources including stellar number density, variations in survey depth, and point spread function properties. While some marginal correlations with these systematic tracers were observed, there is no evidence of bias in the cosmological inference. B-mode power spectra are consistent with zero signal, with no significant residual contamination from E/B-mode leakage. We performed a Bayesian analysis of the pseudo-Cℓ estimates by forward modelling the effects of the mask. Assuming a spatially flat ΛCDM cosmology, we constrained the structure growth parameter S8 = σ8(Ωm/0.3)1/2 = 0.754−0.029+0.027. When combining cosmic shear from KiDS-1000 with baryon acoustic oscillation and redshift space distortion data from recent Sloan Digital Sky Survey (SDSS) measurements of luminous red galaxies, as well as the Lyman-α forest and its cross-correlation with quasars, we tightened these constraints to S8 = 0.771−0.032+0.006. These results are in very good agreement with previous KiDS-1000 and SDSS analyses and confirm a ∼3σ tension with early-Universe constraints from cosmic microwave background experiments.

NIR jets from a clustered region of massive star formation

Context. Massive stars play crucial roles in determining the physical and chemical evolution of galaxies. However, they form deeply embedded in their parental clouds, making it challenging to directly observe these stars and their immediate environments. It is known that accretion and ejection processes are intrinsically related, thus observing the massive protostellar outflows can provide crucial information about the processes governing massive star formation very close to the central engine. Aims. We aim to probe the IRAS 18264-1152 (also known as G19.88-0.53) high-mass star-forming complex in the near infrared (NIR) through its molecular hydrogen (H2) jets to analyse the morphology and composition of the line emitting regions and to compare with other outflow tracers. Methods. We observed the H2 NIR jets via K-band (1.9–2.5 μm) observations obtained with the integral field units VLT/SINFONI and VLT/KMOS. VLT/SINFONI provides the highest NIR angular resolution achieved so far for the central region of IRAS 18264-1152 (~0.2′′). We compared the geometry of the NIR outflows with that of the associated molecular outflow, probed by CO (2-1) emission mapped with the Submillimeter Array. Results. We identify nine point sources in the SINFONI and KMOS fields of view. Four of these display a rising continuum in the K-band and are Brγ emitters, revealing that they are young, potentially jet-driving sources. The spectro-imaging analysis focusses on the H2 jets, for which we derived visual extinction, temperature, column density, area, and mass. The intensity, velocity, and excitation maps based on H2 emission strongly support the existence of a protostellar cluster in this region, with at least two (and up to four) different large-scale outflows, found through the NIR and radio observations. We compare our results with those found in the literature and find good agreement in the outflow morphology. This multi-wavelength comparison also allows us to derive a stellar density of ~4000 stars pc−3. Conclusions. Our study reveals the presence of several outflows driven by young sources from a forming cluster of young, massive stars, demonstrating the utility of such NIR observations for characterising massive star-forming regions. Moreover, the derived stellar number density together with the geometry of the outflows suggest that stars can form in a relatively ordered manner in this cluster.

Weighing the Galactic disk using phase-space spirals

We present a new method for inferring the gravitational potential of the Galactic disk, using the time-varying structure of a phase-space spiral in the (z, w)-plane (where z and w represent vertical position and vertical velocity). Our method of inference extracts information from the shape of the spiral and disregards the bulk density distribution that is usually used to perform dynamical mass measurements. In this manner, it is complementary to traditional methods that are based on the assumption of a steady state. Our method consists of fitting an analytical model for the phase-space spiral to data, where the spiral is seen as a perturbation of the stellar number density in the (z, w)-plane. We tested our method on one-dimensional simulations, which were initiated in a steady state and then perturbed by an external force similar to that of a passing satellite. We were able to retrieve the true gravitational potentials of the simulations with high accuracy. The gravitational potential at 400–500 parsec distances from the disk mid-plane was inferred with an error of only a few percent. This is the first paper of a series in which we plan to test and refine our method on more complex simulations, as well as apply our method to Gaia data.

Mapping the Galactic Disk with the LAMOST and Gaia Red Clump Sample. VII. The Stellar Disk Structure Revealed by the Mono-abundance Populations

Abstract Using a sample of 96,201 primary red clump stars selected from the LAMOST and Gaia surveys, we investigate the stellar structure of the Galactic disk. The sample stars show two separated sequences of high-[α/Fe] and low-[α/Fe] in the [α/Fe]–[Fe/H] plane. We divide the sample stars into five mono-abundance populations (MAPs) with different ranges of [α/Fe] and [Fe/H], named as the high-[α/Fe], high-[α/Fe] and high-[Fe/H], low-[Fe/H], solar, high-[Fe/H] MAPs, respectively. We present the stellar number density distributions in the R–Z plane, and the scale heights and scale lengths of the individual MAPs by fitting their vertical and radial density profiles. The vertical profiles, the variation trend of scale height with the Galactocentric radius, indicate that there is a clear disk flare in the outer disk both for the low-[α/Fe] and the high-[α/Fe] MAPs. While the radial surface-density profiles show a peak radius of 7 kpc and 8 kpc for the high-[α/Fe] and low-[α/Fe] MAPs, respectively. We also investigate the correlation between the mean rotation velocity and metallicity of the individual MAPs, and find that the mean rotation velocities are well separated and show different trends between the high-[α/Fe] and the low-[α/Fe] MAPs. Finally, we discuss the character of the high-[α/Fe] and high-[Fe/H] MAP and find that it is more similar to the high-[α/Fe] MAP either in the radial and vertical density profiles or in the rotation velocity.

Dynamics of Planetary Systems within Star Clusters: Aspects of the Solar System’s Early Evolution

Abstract Most planetary systems—including our own—are born within stellar clusters, where interactions with neighboring stars can help shape the system architecture. This paper develops an orbit-averaged formalism to characterize the cluster’s mean-field effects, as well as the physics of long-period stellar encounters. Our secular approach allows for an analytic description of the dynamical consequences of the cluster environment on its constituent planetary systems. We analyze special cases of the resulting Hamiltonian, corresponding to eccentricity evolution driven by planar encounters, as well as hyperbolic perturbations upon dissipative disks. We subsequently apply our results to the early evolution of our solar system, where the cluster’s collective potential perturbs the solar system’s plane, and stellar encounters act to increase the velocity dispersion of the Kuiper Belt. Our results are twofold. First, we find that cluster effects can alter the mean plane of the solar system by ≲1° and are thus insufficient to explain the ψ ≈ 6° obliquity of the Sun. Second, we delineate the extent to which stellar flybys excite the orbital dispersion of the cold classical Kuiper Belt and show that while stellar flybys may grow the cold belt’s inclination by the observed amount, the resulting distribution is incompatible with the data. Correspondingly, our calculations place an upper limit on the product of the stellar number density and residence time of the Sun in its birth cluster, η τ ≲ 2 × 104 Myr pc−3.

New constraints on the structure of the nuclear stellar cluster of the Milky Way from star counts and MIR imaging

Context. The Milky Way nuclear star cluster (MWNSC) is a crucial laboratory for studying the galactic nuclei of other galaxies, but its properties have not been determined unambiguously until now. Aims. We aim to study the size and spatial structure of the MWNSC. Methods. This study uses data and methods that address potential shortcomings of previous studies on the topic. We use 0.2″ angular resolution Ks data to create a stellar density map in the central 86.4 pc × 21 pc at the Galactic center. We include data from selected adaptive-optics-assisted images obtained for the inner parsecs. In addition, we use Spitzer/IRAC mid-infrared (MIR) images. We model the Galactic bulge and the nuclear stellar disk in order to subtract them from the MWNSC. Finally, we fit a Sérsic model to the MWNSC and investigate its symmetry. Results. Our results are consistent with previous work. The MWNSC is flattened with an axis ratio of q = 0.71 ± 0.10, an effective radius of Re = (5.1 ± 1.0) pc, and a Sérsic index of n = 2.2 ± 0.7. Its major axis may be tilted out of the Galactic plane by up to −10°. The distribution of the giants brighter than the Red Clump (RC) is found to be significantly flatter than the distribution of the faint stars. We investigate the 3D structure of the central stellar cusp using our results on the MWNSC structure on large scales to constrain the deprojection of the measured stellar surface number density, obtaining a value of the 3D inner power law of γ = 1.38 ± 0.06sys ± 0.01stat. Conclusions. The MWNSC shares its main properties with other extragalactic NSCs found in spiral galaxies. The differences in the structure between bright giants and RC stars might be related to the existence of not completely mixed populations of different ages. This may hint at recent growth of the MWNSC through star formation or cluster accretion.

Exchange of material between solar systems by random stellar encounters

AbstractIt is shown that a mechanism involving only random motion of the sun with respect to the surrounding star field can account for the ~1 per 25 Myr characteristic frequency of large cometary impacts on Earth. In the proposed mechanism, the sun travels through the Oort Cloud of an encounter star, most typically a Type M dwarf, while the dwarf flies through the Oort cloud of our Sun. As a result, Oort Cloud objects from our Solar System are precipitated in large numbers to impact planets in the dwarf star system, while the dwarf's Oort Cloud objects are destabilized to impact planets in our Solar System. It is shown that it is this exchange of Oort cloud object between stellar systems, rather than the precipitation of Oort Cloud objects within a stellar system, that can account for the apparent periodicity of mass extinctions. Because the sun is more massive than ~90% of stars, its Oort cloud extends further, resulting in it delivering about a factor of three more bombardments on other solar systems than our Solar System receives. About 60% of the bombardments on our Solar System are found to be delivered by Type M dwarfs, about 20% by type K dwarfs, with the remaining 20% being delivered by stars of type G or larger. Foreign star Oort cloud objects can be captured by our Sun at typical ranges of 10 AU, resulting in a cometary approach to perihelion within a few years. It is found that assuming an effective Oort Cloud radius of 40 000 AU for a star of solar mass, increasing in size with the square root of the mass, accounts for the observed characteristic frequency of mass extinction events on Earth, given the local stellar number density of 0.003 stars per cubic light year. The frequency of mass extinction events in other solar systems would increase or decrease in linear proportion to the local stellar number density. It is shown that this exchange of materials between solar systems during close stellar encounters could be an important mechanism for spreading life throughout the galaxy. Implications for the evolution of life on Earth and in other solar systems are discussed.

Asymmetric Mean Metallicity Distribution of the Milky Way’s Disk

Abstract I present the mean metallicity distribution of stars in the Milky Way based on photometry from the Sloan Digital Sky Survey. I utilize an empirically calibrated set of stellar isochrones developed in previous work to estimate the metallicities of individual stars to a precision of 0.2 dex for reasonably bright stars across the survey area. I also obtain more precise metallicity estimates using priors from the Gaia parallaxes for relatively nearby stars. Close to the Galactic mid-plane ( kpc), a mean metallicity map reveals deviations from the mirror symmetry between the northern and southern hemispheres, displaying wave-like oscillations. The observed metallicity asymmetry structure is almost parallel to the Galactic mid-plane, and coincides with the previously known asymmetry in the stellar number density distribution. This result reinforces the previous notion of the plane-parallel vertical waves propagating through the disk, in which a local metallicity perturbation from the mean vertical metallicity gradient is induced by the phase-space wrapping of stars in the Z–V Z plane. The maximum amplitude of the metallicity asymmetry (Δ[Fe/H] ∼ 0.05) implies that these stars have been pulled away from the Galactic mid-plane by an order of pc as a massive halo substructure such as the Sagittarius dwarf galaxy plunged through the Milky Way. This work provides evidence that the Gaia phase-space spiral may continue out to .

Galactic Stellar Populations from Photometric Metallicity Distribution Functions

Abstract Based on Sloan Digital Sky Survey photometric data, Gu developed a new Monte-Carlo-based method for estimating the stellar metallicity distribution functions (MDFs). This method enables a more reliable determination of MDFs compared with the conventional polynomial-based methods. In this work, MDF determined from the method are well fit by a three-Gaussian model, with peaks at [Fe/H] = −0.68, −1.38, and −1.90, associated with the thick disk, the inner halo, and the outer halo, respectively. The vertical metallicity gradient within 1 &lt; Z &lt; 5 kpc is around R = 8.25 kpc. But the mean radial gradient is almost negligible. The density profile of the thick disk is fitted with a modified double exponential law decaying to a constant at far distance. The scale height and scale length thus estimated are H ≈ 1.13 kpc and L ≈ 3.63 kpc, which are consistent with the results determined from star-count methods in previous studies. The halos are described with a two-axial power-law ellipsoid, and the axis ratios of both the inner halo and the outer halo, inferred from stellar number density in the R–Z plane, are q ih ≈ 0.49 and q oh ≈ 0.61, respectively. It also manifests that the outer halo is more spherical than the inner halo. Moreover, the halo power-law indices estimated are n ih ≈ 3.4 and n oh ≈ 3.1, indicating that the stellar number density of the inner halo changes more steeply than that of outer halo.

Measuring the local matter density usingGaiaDR2

Aims.We determine the total dynamical matter density in the solar neighbourhood using the secondGaiadata release (DR2).Methods.The dynamical matter density distribution is inferred in a framework of a Bayesian hierarchical model, which accounts for position and velocity of all individual stars, as well as the full error covariance matrix of astrometric observables, in a joint fit of the vertical velocity distribution and stellar number density distribution. This was done for eight separate data samples, with different cuts in observed absolute magnitude, each containing about 25 000 stars. The model for the total matter density does not rely on any underlying baryonic model, although we assumed that it is symmetrical, smooth, and monotonically decreasing with distance from the mid-plane.Results.We infer a density distribution which is strongly peaked in the region close to the Galactic plane (≲60 pc), for all eight stellar samples. Assuming a baryonic model and a dark matter halo of constant density, this corresponds to a surplus surface density of approximately 5–9M⊙pc−2. For the Sun’s position and vertical velocity with respect to the Galactic plane, we inferZ⊙ = 4.76 ± 2.27 pc andW⊙ = 7.24 ± 0.19 km s−1.Conclusions.These results suggest a surplus of matter close to the Galactic plane, possibly explained by an underestimated density of cold gas. We discuss possible systematic effects that could bias our result, for example unmodelled non-equilibrium effects, and how to account for such effects in future extensions of this work.

The dynamical matter density in the solar neighbourhood inferred from Gaia DR1

We determine the total dynamical density in the solar neighbourhood using the Tycho-Gaia Astrometric Solution (TGAS) catalogue. Astrometric measurements of proper motion and parallax of stars inform us of both the stellar number density distribution and the velocity distribution of stars close to the plane. Assuming equilibrium, these distributions are interrelated through the local dynamical density. For the first time, we do a full joint fit of the velocity and stellar number density distributions while accounting for the astrometric error of individual stars, in the framework of Bayesian Hierarchical Modelling. We use a sample of stars whose distance extends to approximately 160 pc from the Sun. We find a local matter density of $\rho_0 = 0.119^{+0.015}_{-0.012}~M_\odot$pc$^{-3}$, where the result is presented as the median to the posterior distribution, plus/minus its 16th and 84th percentiles. We find the Sun's position above the Galactic plane to be $z_\odot = 15.3^{+2.24}_{-2.16}$ pc, and the Sun's velocity perpendicular to the Galactic plane to be $W_\odot = 7.19^{+0.18}_{-0.18}$ km/s.

Galactic cartography with SkyMapper – I. Population substructure and the stellar number density of the inner halo

The stars within our Galactic halo presents a snapshot of its ongoing growth and evolution, probing galaxy formation directly. Here, we present our first analysis of the stellar halo from detailed maps of Blue Horizontal Branch (BHB) stars drawn from the SkyMapper Southern Sky Survey. To isolate candidate BHB stars from the overall population, we develop a machine-learning approach through the application of an Artificial Neural Network (ANN), resulting in a relatively pure sample of target stars. From this, we derive the absolute $u$ magnitude for the BHB sample to be $\sim2\ mag$, varying slightly with $(v-g)_0$ and $(u-v)_0$ colours. We examine the BHB number density distribution from 5272 candidate stars, deriving a double power-law with a break radius of $r_s = 11.8\pm0.3\ kpc$, and inner and outer slopes of $\alpha_{in} = -2.5\pm0.1$ and $\alpha_{out} = -4.5\pm0.3$ respectively. Through isochrone fitting of simulated BHB stars, we find a colour-age/metallicity correlation, with older/more metal-poor stars being bluer, and establish a parameter to indicate this age (or metallicity) variation. Using this, we construct the three-dimensional population distribution of BHB stars in the halo and identify significant substructure. Finally, in agreement with previous studies, we also identify a systemic age/metallicity shift spanning $\sim3\ kpc$ to $ \sim20\ kpc$ in galactocentric distance.

Protoplanetary disc truncation mechanisms in stellar clusters: comparing external photoevaporation and tidal encounters

Most stars form and spend their early life in regions of enhanced stellar density. Therefore the evolution of protoplanetary discs (PPDs) hosted by such stars are subject to the influence of other members of the cluster. Physically, PPDs might be truncated either by photoevaporation due to ultraviolet flux from massive stars, or tidal truncation due to close stellar encounters. Here we aim to compare the two effects in real cluster environments. In this vein we first review the properties of well studied stellar clusters with a focus on stellar number density, which largely dictates the degree of tidal truncation, and far ultraviolet (FUV) flux, which is indicative of the rate of external photoevaporation. We then review the theoretical PPD truncation radius due to an arbitrary encounter, additionally taking into account the role of eccentric encounters that play a role in hot clusters with a 1D velocity dispersion $\sigma_v > 2$ km/s. Our treatment is then applied statistically to varying local environments to establish a canonical threshold for the local stellar density ($n_{c} > 10^4$ pc$^{-3}$) for which encounters can play a significant role in shaping the distribution of PPD radii over a timescale $\sim 3$ Myr. By combining theoretical mass loss rates due to FUV flux with viscous spreading in a PPD we establish a similar threshold for which a massive disc is completely destroyed by external photoevaporation. Comparing these thresholds in local clusters we find that if either mechanism has a significant impact on the PPD population then photoevaporation is always the dominating influence.