Astrophysics of Galaxies (astro-ph.GA)

  • PDF
    Over 80% of the cosmic baryon density is likely to be distributed in the diffuse, $\gtrsim10^{4}$ K circum- and inter-galactic medium (CGM and IGM respectively). We demonstrate that the dispersion measures (DMs) of samples of localized Fast Radio Bursts (FRBs) can be used to measure the distribution of baryons between the CGM and IGM. We propose to separate the CGM and IGM contributions to FRB DMs by including redshift and mass measurements of intervening galaxies. Using simulated samples of FRB sightlines through intervening galaxy halos and an illustrative model for the CGM, and including realistic observational uncertainties, we show that small samples ($O(10^{1})-O(10^{2})$) of localized FRBs are sensitive to the presence of CGM gas. The fractions of baryons in the CGM and IGM can be accurately estimated with 100 FRBs at $z<1$, and the characteristic radial density profiles of CGM halos may also be possible to constrain. The required samples of localized FRBs are expected to be assembled in the coming few years by instruments such as the Australian Square Kilometre Array Pathfinder, the Very Large Array, and the Deep Synoptic Array.
  • PDF
    The behavior of light scalar dark matter has been a subject of much interest recently as it can lead to interesting small scale behavior. In particular, this can lead to the formation of gravitationally bound clumps for these light scalars, including axions. In Ref. [1] we analyzed the behavior of these clumps, assuming spherical symmetry, allowing for both attractive and repulsive self-interactions. There is a maximum allowed mass for the clumps in the case of attractive interactions, and a minimum radius for the clumps in the case of repulsive interactions, which is saturated at large mass. Here we extend this work to include non-spherically symmetric clumps. Since the system tries to re-organize into a BEC ground state, we consider configurations with a conserved non-zero angular momentum, and construct minimum energy configurations at fixed particle number and fixed angular momentum. We find generalizations of the previous spherically symmetric results. In particular, while there is still a maximum mass for the attractive case, its value increases with angular momentum. Also, the minimum radius in the repulsive case is raised to higher radii. In a forthcoming paper we shall investigate the possibility of resonance of axion clumps into both visible and hidden sector photons, and analyze how the altered mass and radius from non-zero angular momentum affects the resonance.
  • PDF
    We present the results of our automatic search for proximate damped Ly$\alpha$ absorption (PDLA) systems in the quasar spectra from the Sloan Digital Sky Survey Data Release 12. We constrain our search to those PDLAs lying within 1500 km s$^{-1}$ from the quasar to make sure that the broad DLA absorption trough masks most of the strong Ly$\alpha$ emission from the broad line region (BLR) of the quasar. When the Ly$\alpha$ emission from the BLR is blocked by these so-called eclipsing DLAs, narrow Ly$\alpha$ emission from the host galaxy could be revealed as a narrow emission line (NEL) in the DLA trough. We define a statistical sample of 399 eclipsing DLAs with log$N$(HI)$\ge$21.10. We divide our statistical sample into three subsamples based on the strength of the NEL detected in the DLA trough. By studying the stacked spectra of these subsamples, we found that absorption from high ionization species are stronger in DLAs with stronger NEL in their absorption core. Moreover, absorption from the excited states of species like SiII are also stronger in DLAs with stronger NEL. We also found no correlation between the luminosity of the Ly$\alpha$ NEL and the quasar luminosity. These observations are consistent with a scenario in which the DLAs with stronger NEL are denser and physically closer to the quasar. We propose that these eclipsing DLAs could be the product of the interaction between infalling and outflowing gas. High resolution spectroscopic observation would be needed to shed some light on the nature of these eclipsing DLAs.
  • PDF
    If the Galaxy is axisymmetric and in dynamical equilibrium, we expect negligible fluctuations in the residual line-of-sight velocity field. However, non-axisymmetric structures like a bar, spiral arms and merger events can generate velocity fluctuations. Recent results using the APOGEE survey find significant fluctuations in velocity for stars in the midplane (|z|< 0.25 kpc) and out to 5 kpc, which suggests that the dynamical influence of the Milky Way's bar extends out to the Solar neighborhood. Their measured power spectrum has a characteristic amplitude of 11 km/s on a scale of ~ 2.5 kpc. The existence of large streaming motions on these scales has important implications for determining the Sun's motion about the Galactic Centre. Using red clump stars from the GALAH and APOGEE surveys, we map the line-of-sight velocity field around the Sun out to distances of 5 kpc and up to 1.25 kpc from the Galactic Plane. By subtracting a smooth axisymmetric model for the velocity field, we study the residual velocity fluctuations and compare our findings with mock survey generated by Galaxia based on an axisymmetric, steady state model. We find negligible large-scale fluctuations away from the plane. In the mid-plane, we reproduce the earlier APOGEE power spectrum results but with 20\% smaller amplitude (9.5 km/s) after taking a few systematic effects into account (e.g. volume completeness). The amplitude power is further reduced to 6.7 km/s if a flexible axisymmetric model is used. Additionally, our mock simulations show that, in the plane, the distances are underestimated for high mass red clump stars and this can lead to spurious power with amplitude of about 5.5 km/s. Taking this into account, we estimate the amplitude of real fluctuations to be less than 4.2 km/s, about a factor of three less than the previous result from APOGEE.
  • PDF
    We cross-match objects from several different astronomical catalogs to determine the absolute proper proper motions of stars within the 30-arcmin radius fields of 115 Milky-Way globular clusters with the accuracy of 1--2~mas/yr. The proper motions are based on positional data recovered from the USNO-B1, 2MASS, URAT1, ALLWISE, UCAC5, and GAIA DR1 surveys with up to 10 positions spanning an epoch difference of up to $\sim$~65~years, and reduced to GAIA DR1 TGAS frame using UCAC5 as the reference catalog. Cluster members are photometrically identified by selecting horizontal- and red-giant branch stars on color-magnitude diagrams, and the mean absolute proper motions of the clusters with a typical formal error of $\sim$~0.4~mas/yr are computed by averaging the proper motions of selected members. The inferred absolute proper motions of clusters are combined with available radial-velocity data and heliocentric distance estimates to compute the cluster orbits in terms of the Galactic potential models based on Miyamoto and Nagai disk, Hernquist spheroid, and modified isothermal dark-matter halo (axisymmetric model without a bar) and the same model + rotating Ferre's bar (non-axisymmetric). Five distant clusters have higher-than-escape velocities, most likely due to large errors of computed transversal velocities, whereas the computed orbits of all other clusters remain bound to the Galaxy. Unlike previously published results, we find the bar to affect substantially the orbits of most of the clusters, even those at large Galactocentric distances, bringing appreciable chaotization, especially in the portions of the orbits close to the Galactic center, and stretching out the orbits of some of the thick-disk clusters.
  • PDF
    We present the relation between stellar specific angular momentum $j_*$, stellar mass $M_*$, and bulge-to-total light ratio $\beta$ for THINGS, CALIFA and Romanowsky \& Fall datasets, exploring the existence of a fundamental plane between these parameters as first suggested by Obreschkow \& Glazebrook. Our best-fit $M_*-j_*$ relation yields a slope of $\alpha = 1.03 \pm 0.11$ with a trivariate fit including $\beta$. When ignoring the effect of $\beta$, the exponent $\alpha = 0.56 \pm 0.06$ is consistent with $\alpha = 2/3$ predicted for dark matter halos. There is a linear $\beta - j_*/M_*$ relation for $\beta \lesssim 0.4$, exhibiting a general trend of increasing $\beta$ with decreasing $j_*/M_*$. Galaxies with $\beta \gtrsim 0.4$ have higher $j_*$ than predicted by the relation. Pseudobulge galaxies have preferentially lower $\beta$ for a given $j_*/M_*$ than galaxies that contain classical bulges. Pseudobulge galaxies follow a well-defined track in $\beta - j_*/M_*$ space, consistent with Obreschkow \& Glazebrook, while galaxies with classical bulges do not. These results are consistent with the hypothesis that while growth in either bulge type is linked to a decrease in $j_*/M_*$, the mechanisms that build pseudobulges seem to be less efficient at increasing bulge mass per decrease in specific angular momentum than those that build classical bulges.
  • PDF
    Continuum emissions from dust grains are used as a general probe to constrain the initial physical conditions of molecular dense cores where new stars may born. To get as much information as possible from dust emissions, we have developed a tool, named as $COREGA$, which is capable of identifying positions of dense cores, optimizing a three-dimensional model for the dense cores with well characterized uncertainties. $COREGA$ can also estimate the physical properties of dense cores, such as density, temperature, and dust emissivity, through analyzing multi-wavelength dust continuum data sets. In the numerical tests on $COREGA$, the results of fitting simulated data are consistent with initial built-in parameters. We also demonstrate $COREGA$ by adding random gaussian noises with Monte Carlo methods and show that the results are stable against varying observational noise intensities within certain levels. A beam size $<$ 3 arcsec and rms $<$ 0.2mJy/pixel (1 pixel = 0.1") is needed for ALMA to distinguish different collapse models, such as power law and Bonner-Ebert sphere, during continuum observations of massive dense cores in Orion molecular cloud. Based on its advanced algorithm, $COREGA$ is capable of giving a quick and deep analysis on dust cores.
  • PDF
    We present optical photometry of Hubble Space Telescope (HST) ACS/WFC data of the resolved stellar populations in the outer disc of the dwarf irregular galaxy DDO 154. The photometry reveals that young main sequence stars are almost absent from the outermost HI disc. Instead, most are clustered near the main stellar component of the galaxy. We constrain the stellar initial mass function (IMF) by comparing the luminosity function of the main sequence stars to simulated stellar populations assuming a constant star formation rate over the dynamical timescale. The best-fitting IMF is deficient in high mass stars compared to a canonical Kroupa IMF, with a best-fit slope $\alpha = -2.45$ and upper mass limit $M_U = 16\ M_{\odot}$. This top-light IMF is consistent with predictions of the Integrated Galaxy-wide IMF theory. Combining the HST images with HI data from The HI Nearby Galaxy Survey Treasury (THINGS) we determine the star formation law (SFL) in the outer disc. The fit has a power law exponent $N = 2.92 \pm0.22$ and zero point $A=4.47 \pm 0.65 \times 10^{-7} \ M_{\odot} \ \text{yr}^{-1} \ \text{kpc}^{-2}$. This is depressed compared to the Kennicutt-Schmidt Star Formation Law, but consistent with weak star formation observed in diffuse HI environments. Extrapolating the SFL over the outer disc implies that there could be significant star formation occurring that is not detectable in H$\alpha$. Last, we determine the Toomre stability parameter $Q$ of the outer disc of DDO 154 using the THINGS HI rotation curve and velocity dispersion map. 72% of the HI in our field has $Q\leq 4$ and this incorporates 96% of the observed MS stars. Hence 28% of the HI in the field is largely dormant.
  • PDF
    We use the SDSS-Gaia catalogue to search for material stripped off the atypical globular cluster $\omega$ Centauri. As actions are conserved under slow changes of the potential, this allows identification of groups of stars with a common accretion history. We devise a method to assess the significance of halo substructures based on their clustering in action space, using metallicity as a secondary check. This is validated against smooth models and numerically constructed stellar halos. We identify a number of new substructures in the SDSS-Gaia catalogue, including 7 high significance, high energy and retrograde ones. Using a simple model of the accretion of the progenitor of the $\omega$ Centauri, we argue for the possible association of up to 4 of our new substructures (labelled Rg1, Rg3, Rg4 and Rg6) with this event. This sets a minimum mass of $5 \times 10^8 M_\odot$ for the progenitor, so as to bring $\omega$ Centauri to its current location in action -- energy space. Our proposal can be tested by high resolution spectroscopy of the candidates to look for the unusual abundance patterns possessed by $\omega$ Centauri stars.
  • PDF
    We used the Tycho-Gaia Astrometric Solution catalogue, part of the Gaia Data Release 1, to search for candidate astrometric microlensing events expected to occur within the remaining lifetime of the Gaia satellite. Our search yielded one promising candidate. We predict that the nearby DQ type white dwarf LAWD 37 (WD 1142-645) will lens a background star and will reach closest approach on November 11th 2019 ($\pm$ 4 days) with impact parameter $380\pm10$ mas. This will produce an apparent maximum deviation of the source position of $2.8\pm0.1$ mas. In the most propitious circumstance, Gaia will be able to determine the mass of LAWD 37 to $\sim3\%$. This mass determination will provide an independent check on atmospheric models of white dwarfs with helium rich atmospheres, as well as tests of white dwarf mass radius relationships and evolutionary theory.
  • PDF
    We present new atomic hydrogen (HI) observations of the collisional ring galaxy NGC 922 obtained using the Australia Telescope Compact Array. Our observations reveal for the first time the vast extent of the HI disc of this galaxy. The HI morphology and kinematics of NGC 922 show that this galaxy is not the product of a simple drop-through interaction, but has a more complex interaction history. The integrated HI flux density of NGC 922 from our observations is 24.7 Jy km s$^{-1}$, which is within the error of the flux value obtained using the $64$-m Parkes radio telescope. This flux density translates to a total HI mass of $1.1*10^{10}$ M$_{\circ}$ and corresponds to an HI to total mass fraction (M$_{HI}$/M$_{tot}$) of approximately $0.11$. The gaseous structures of NGC 922 are more extended to the north and include an HI tail that has a projected physical length of $8$ kpc. Gas warps are also evident in the velocity field of NGC 922 and are more prominent on the approaching and the western side of the disc. In comparison with a large sample of star-forming galaxies in the local Universe, NGC 922 possesses a high gas fraction relative to galaxies with a similar stellar mass of ~$10^{10.4}$ M$_{\circ}$, and exhibits a high specific star formation rate.
  • PDF
    Previous radio recombination line (RRL) observations of dust clumps identified in the APEX Telescope Large Area Survey of the Galaxy (ATLASGAL) have led to the detection of a large number of RRLs in the 3mm range. Here, we aim to study their excitation with shorter wavelength (sub)millimeter radio recombination line (submm-RRL) observations. We made observations of submm-RRLs with low principal quantum numbers ($n$ $\leq$ 30) using the APEX 12 m telescope, toward 104 HII regions associated with massive dust clumps from ATLASGAL. The observations covered the H25$\alpha$, H28$\alpha$, and H35$\beta$ transitions. Toward a small subsample the H26$\alpha$, H27$\alpha$, H29$\alpha$, and H30$\alpha$ lines were observed to avoid contamination by molecular lines at adjacent frequencies. We have detected submm-RRLs (signal-to-noise $\geq$ 3 $\sigma$) from compact HII regions embedded within 93 clumps. The submm-RRLs are approximately a factor of two brighter than the mm-RRLs and consistent with optically thin emission in local thermodynamic equilibrium (LTE). The average ratio (0.31) of the measured H35$\beta$/H28$\alpha$ fluxes is close to the LTE value of 0.28. No indication of RRL maser emission has been found. The Lyman photon flux, bolometric, and submm-RRL luminosities toward the submm-RRL detected sources present significant correlations. The trends of dust temperature and the ratio of bolometric luminosity to clump mass, $L_{\rm bol}/M_{\rm clump}$, indicate that the HII regions are related to the most massive and luminous clumps. By estimating the production rate of ionizing photons, $Q$, from the submm-RRL flux, we find that the $Q$(H28$\alpha$) measurements provide estimates of the Lyman continuum photon flux consistent with those determined from 5 GHz radio continuum emission.
  • PDF
    We analyse the geometry of the Pipe Nebula, drawn by the distribution ($Q$-spatial parameter) and hierarchy ($\Lambda$ spatial segregation) of column density peaks previously detected and catalogued. By analysing the mass and volume density of the cores, we determine that both variables shown to be spatially segregated and with a high degree of substructure. Given the early evolutionary state of the Pipe Nebula, our results suggest that both, mass and volume-density segregations, may be primordial, in the sense of appearing in the early phases of the chain of physical mechanisms which conform the cluster-formation process. We also propose that volume density, and not mass, is the parameter that most clearly determines the initial spatial distribution of pre-stellar cores.
  • PDF
    We simultaneously study the growth of errors and the question of the faithfulness of simulations of $N$-body systems. The errors are quantified through the numerical reversibility of trajectories of small-$N$ spherical systems integrated to high accuracy. Initially, the errors add randomly, before exponential divergence sets in. Though the exponentiation rate is virtually independent of $N$, the instability saturates at scales $1/\sqrt{N}$. This is interpreted by adopting a model due to Goodman, Heggie \& Hut (1993). In the third phase, the (diminished) growth is initially driven by multiplicative enhancement of errors as in the exponential stage. It is then qualitatively different for the errors in the phase space variables and the mean field conserved quantities (energy and momentum); the former grow systematically through phase mixing while the latter grow diffusively. For energy, the $N$-variation of the `relaxation time' of error growth follows expectations of two-body relaxation theory. This is not the case for angular momentum, at least up to the particle numbers and timescales considered, and even less so for the velocities. Due to increasingly smaller saturation scales, the information loss associated with the exponential instability decreases with $N$, especially when viewed in terms of the mean-field conserved quantities. Indeed, the dynamical entropy vanishes at any finite resolution as $N \rightarrow \infty$. In this sense there is convergence to the collisionless limit and confidence that numerical simulations may faithfully represent it, despite the exponential instability and loss of information on phase space trajectories. Nevertheless, the rapid initial growth of errors. and the relatively slow $N$-variation in its saturation, point to the slowness of the convergence.
  • PDF
    CONTEXT.The first Gaia Data Release (DR1) significantly improved the previously available proper motions for the majority of the Tycho-2 stars. AIMS. We want to detect runaway stars using Gaia DR1 proper motions and compare our results with previous searches. METHODS. Runaway O stars and BA supergiants are detected using a 2-D proper-motion method. The sample is selected using Simbad, spectra from our GOSSS project, literature spectral types, and photometry processed using CHORIZOS. RESULTS. We detect 76 runaway stars, 17 (possibly 19) of them with no prior identification as such, with an estimated detection rate of approximately one half of the real runaway fraction. An age effect appears to be present, with objects of spectral subtype B1 and later having travelled for longer distances than runaways of earlier subtypes. We also tentatively propose that the fraction of runaways is lower among BA supergiants that among O stars but further studies using future Gaia data releases are needed to confirm this. The frequency of fast rotators is high among runaway O stars, which indicates that a significant fraction of them (and possibly a majority) is produced in supernova explosions.
  • PDF
    Using visible, radio, microwave, and sub-mm data, we study several lines of sight toward stars generally closer than 1 kpc on a component-by-component basis. We derive the component structure seen in absorption at visible wavelengths from Ca II, Ca I, K I, CH, CH$^{+}\!,$ and CN and compare it to emission from H I, CO and its isotopologues, and C$^{+}$ from the GOT C+ survey. The correspondence between components in emission and absorption help create a more unified picture of diffuse atomic and molecular gas in the interstellar medium. We also discuss how these tracers are related to the CO-dark H$_{2}$ gas probed by C$^{+}$ emission and discuss the kinematic connections among the species observed.
  • PDF
    The majority of low-ionization nuclear emission-line regions (LINERs) harbor supermassive black holes with very low accretion rates. However, the accretion flows do not produce enough ionizing photons to power the emission lines emitted on scales of ~100 pc, therefore additional sources of power are required. We present and analyze Hubble Space Telescope spectra of three nearby luminous LINERs that are spatially resolved on scales of ~9 pc. The targets have multiple indicators of an accreting black hole, as well as a deficient ionizing photon budget. We measure diagnostic emission line ratios as a function of distance from the nucleus and compare them to models for different excitation mechanisms: shocks, photoionization by the accreting black hole, and photoionization by young or old hot stars. The physical model that best describes our targets comprises a low-luminosity, accretion-powered active nucleus that photoionizes the gas within ~20 pc of the galaxy center and shocks that excite the gas at larger distances. We conclude that, in LINERs with low-luminosity active nuclei, shocks by jets or other outflows are crucial in exciting the gas in and around the nucleus, as suggested by other authors.
  • PDF
    The light curves of Type IIn supernovae are dominated by the radiative energy released through the interaction of the supernova shockwaves with their dense circumstellar medium (CSM). The ultraluminous Type IIn supernova SN 2010jl exhibits an infrared emission component that is in excess of the extrapolated UV-optical spectrum as early as a few weeks post-explosion. This emission has been attributed by some as evidence for rapid formation of dust in the cooling postshock CSM. We investigate the physical processes that may inhibit or facilitate the formation of dust in the CSM. When only radiative cooling is considered, the temperature of the dense shocked gas rapidly drops below the dust condensation temperature. However, by accounting for the heating of the postshock gas by the downstream radiation from the shock, we show that dust formation is inhibited until the radiation from the shock weakens, as the shock propagates into the less dense outer regions of the CSM. In SN 2010jl dust formation can therefore only commence after day 380. Only the IR emission since that epoch can be attributed to the newly formed CSM dust. Observations on day 460 and later show that the IR luminosity exceeds the UV-optical luminosity. The post-shock dust cannot extinct the radiation emitted by the expanding SN shock. Therefore, its IR emission must be powered by an interior source, which we identify with the reverse shock propagating through the SN ejecta. IR emission before day 380 must therefore be an IR echo from preexisting CSM dust.
  • PDF
    We present a semi-empirical, largely model-independent approach for estimating Galactic birth radii, r_birth, for Milky Way disk stars. The technique relies on the justifiable assumption that a negative radial metallicity gradient in the interstellar medium (ISM) existed for most of the disk lifetime. Stars are projected back to their birth positions according to the observationally derived age and [Fe/H] with no kinematical information required. Applying our approach to the AMBRE:HARPS and HARPS-GTO local samples, we show that we can constrain the ISM metallicity evolution with Galactic radius and cosmic time, [Fe/H]_ISM(r, t), by requiring a physically meaningful r_birth distribution. We find that the data are consistent with an ISM radial metallicity gradient that flattens with time from ~-0.15 dex/kpc at the beginning of disk formation, to its measured present-day value (-0.07 dex/kpc). We present several chemo-kinematical relations in terms of mono-r_birth populations. One remarkable result is that the kinematically hottest stars would have been born locally or in the outer disk, consistent with thick disk formation from the nested flares of mono-age populations and predictions from cosmological simulations. This phenomenon can be also seen in the observed age-velocity dispersion relation, in that its upper boundary is dominated by stars born at larger radii. We also find that the flatness of the local age-metallicity relation (AMR) is the result of the superposition of the AMRs of mono-r_birth populations, each with a well-defined negative slope. The solar birth radius is estimated to be 7.3+-0.6 kpc, for a current Galactocentric radius of 8 kpc.
  • PDF
    Cold stellar streams---produced by tidal disruptions of clusters---are long-lived, coherent dynamical features in the halo of the Milky Way. Due to their different ages and different positions in phase space, different streams tell us different things about the Galaxy. Here we employ a Cramer--Rao (CRLB) or Fisher-matrix approach to understand the quantitative information content in eleven known streams (ATLAS, GD-1, Hermus, Kwando, Orinoco, PS1A, PS1C, PS1D, PS1E, Sangarius and Triangulum). This approach depends on a generative model, which we have developed previously, and which permits calculation of derivatives of predicted stream properties with respect to Galaxy and stream parameters. We find that in simple analytic models of the Milky Way, streams on eccentric orbits contain the most information about the halo shape. For each stream, there are near-degeneracies between dark-matter-halo properties and parameters of the bulge, the disk, and the stream progenitor, but simultaneous fitting of multiple streams will constrain all parameters at the percent level. At this precision, simulated dark matter halos deviate from simple analytic parametrizations, so we add an expansion of basis functions to give the gravitational potential more freedom. As freedom increases, the information about the halo reduces overall, and it becomes more localized to the current position of the stream. In the limit of high model freedom, a stellar stream appears to measure the local acceleration at its current position; this motivates thinking about future non-parametric approaches. The CRLB formalism also permits us to assess the value of future measurements of stellar velocities, distances, and proper motions. We show that kinematic measurements of stream stars are essential for producing competitive constraints on the distribution of dark matter, which bodes well for stream studies in the age of Gaia.
  • PDF
    We implement a state-of-the-art treatment of the processes affecting the production and Interstellar Medium (ISM) evolution of carbonaceous and silicate dust grains within SPH simulations. We trace the dust grain size distribution by means of a two-size approximation. We test our method on zoom-in simulations of four massive ($M_{200} \geq 3 \times 10^{14} M_{\odot}$) galaxy clusters. We predict that during the early stages of assembly of the cluster at $z \gtrsim 3$, where the star formation activity is at its maximum in our simulations, the proto-cluster regions are rich of dusty gas. Compared to the case in which only dust production in stellar ejecta is active, if we include processes occurring in the cold ISM,the dust content is enhanced by a factor $2-3$. However, the dust properties in this stage turn out to be significantly different than those observationally derived for the \it average Milky Way dust, and commonly adopted in calculations of dust reprocessing. We show that these differences may have a strong impact on the predicted spectral energy distributions. At low redshift in star forming regions our model reproduces reasonably well the trend of dust abundances over metallicity as observed in local galaxies. However we under-produce by a factor of 2 to 3 the total dust content of clusters estimated observationally at low redshift, $z \lesssim 0.5$ using IRAS, Planck and Herschel satellites data. This discrepancy can be solved by decreasing the efficiency of sputtering which erodes dust grains in the hot Intracluster Medium (ICM).
  • PDF
    Samples of reliably identified halo wide binaries are scarce. If reasonably free from selection effects and with a small degree of contamination by chance alignments, these wide binaries become a powerful dynamical tool, having provided one of the very few experiments capable of constraining the nature of dark matter in the Milky Way halo. Currently, however, the best available sample of halo wide binaries is confined to the solar neighborhood, and is plagued by small number statistics at the widest separations. We present the results of a program aimed to probe the wide binary population of the Galactic halo at significantly larger distances, and which informs future searches that could improve the statistics by orders of magnitude. Halo stars were taken from the Sloan Digital Sky Survey after analyzing the Galactic orbits of stars in the reduced proper motion diagram. We then select candidate binaries by searching for pairs with small differences in proper motion and small projected separation on the sky. Using medium-resolution spectroscopy, a subsample of candidates is validated via radial velocities, finding that about 68% of candidate pairs up to 20$^{\prime\prime}$ separation can be considered genuine halo wide binaries, to the limits of the available data. Photometric distance estimates show that the components of our validated candidates lie at the same distances, independently confirming the robustness of the selection method. These results should prove valuable to guide the optimal assembly of larger catalogs of halo wide binaries from upcoming large databases such as Gaia and LSST.
  • PDF
    We use ALMA observations of four sub-millimetre galaxies (SMGs) at $z\sim2-3$ to investigate the spatially resolved properties of the inter-stellar medium (ISM) at scales of 1--5 kpc (0.1--0.6$''$). The velocity fields of our sources, traced by the $^{12}$CO($J$=3-2) emission, are consistent with disk rotation to first order, implying average dynamical masses of $\sim$3$\times10^{11}$M$_{\odot}$ within two half-light radii. Through a Bayesian approach we investigate the uncertainties inherent to dynamically constraining total gas masses. We explore the covariance between the stellar mass-to-light ratio and CO-to-H$_{2}$ conversion factor, $\alpha_{\rm CO}$, finding values of $\alpha_{\rm CO}=1.1^{+0.8}_{-0.7}$ for dark matter fractions of 15 \%. We show that the resolved spatial distribution of the gas and dust continuum can be uncorrelated to the stellar emission, challenging energy balance assumptions in global SED fitting. Through a stacking analysis of the resolved radial profiles of the CO(3-2), stellar and dust continuum emission in SMG samples, we find that the cool molecular gas emission in these sources (radii $\sim$5--14 kpc) is clearly more extended than the rest-frame $\sim$250 $\mu$m dust continuum by a factor $>2$. We propose that assuming a constant dust-to-gas ratio, this apparent difference in sizes can be explained by temperature and optical-depth gradients alone. Our results suggest that caution must be exercised when extrapolating morphological properties of dust continuum observations to conclusions about the molecular gas phase of the ISM.
  • PDF
    Recent work has suggested extreme [OIII] emitting star-forming galaxies are important to reionization. Relatedly, [OIII]/[OII] has been put forward as an indirect estimator of the Lyman Continuum (LyC) escape fraction ($f_{esc}$) at $z\gtrsim4.5$ when the opaque IGM renders LyC photons unobservable. Using deep archival U-band (VLT/VIMOS) imaging of a recently confirmed overdensity at $z\sim3.5$ we calculate tight constraints on $f_{esc}$ for a sample (N=73) dominated by extreme [OIII] emitters. We find no Lyman Continuum signal ($f_{esc}^{rel} < 6.3^{+0.7}_{-0.7} \%$ at $1\sigma$) in a deep U-band stack of our sample (31.98 mag at 1$\sigma$). This constraint is in agreement with recent studies of star-forming galaxies spanning $z\sim1-4$ that have found very low average $f_{esc}$. Despite the galaxies in our study having an estimated average rest-frame EW([OIII]$\lambda5007$)$\sim400\AA$ and [OIII]/[OII]$\sim 4$ from composite SED-fitting, we find no LyC detection, which brings into question the potential of [OIII]/[OII] as an effective probe of the LyC--a majority of LyC emitters have [OIII]/[OII]$>3$, but we establish here that [OIII]/[OII]$>3$ does not guarantee significant LyC leakage for a population. Since even extreme star-forming galaxies are unable to produce the $f_{esc}\sim10-15\%$ required by most theoretical calculations for star-forming galaxies to drive reionization, there must either be a rapid evolution of $f_{esc}$ between $z\sim3.5$ and the Epoch of Reionization, or hitherto observationally unstudied sources (e.g. ultra-faint low-mass galaxies with $\log(M/M_\odot)\sim7-8.5$) must make an outsized contribution to reionization.
  • PDF
    The degree of fractal substructure in molecular clouds can be quantified by comparing them with Fractional Brownian Motion (FBM) surfaces or volumes. These fields are self-similar over all length scales and characterised by a drift exponent $H$, which describes the structural roughness. Given that the structure of molecular clouds and the initial structure of star clusters are almost certainly linked, it would be advantageous to also apply this analysis to clusters. Currently, the structure of star clusters is often quantified by applying $\mathcal{Q}$ analysis. $\mathcal{Q}$ values from observed targets are interpreted by comparing them with those from artificial clusters. These are typically generated using a Box-Fractal (BF) or Radial Density Profile (RDP) model. We present a single cluster model, based on FBM, as an alternative to these models. Here, the structure is parameterised by $H$, and the standard deviation of the log-surface/volume density $\sigma$. The FBM model is able to reproduce both centrally concentrated and substructured clusters, and is able to provide a much better match to observations than the BF model. We show that $\mathcal{Q}$ analysis is unable to estimate FBM parameters. Therefore, we develop and train a machine learning algorithm which can estimate values of $H$ and $\sigma$, with uncertainties. This provides us with a powerful method for quantifying the structure of star clusters in terms which relate to the structure of molecular clouds. We use the algorithm to estimate the $H$ and $\sigma$ for several young star clusters, some of which have no measurable BF or RDP analogue.