Astrophysics

Since the discovery of hypervelocity stars in 2005, it has been widely believed that only the disruption of a binary system by a supermassive black hole at the Galactic center (GC), that is, the so-called Hills mechanism, is capable of accelerating stars to beyond the Galactic escape velocity. In the meantime, however, driven by the Gaia space mission, there is mounting evidence that many of the most extreme high-velocity early-type stars at high Galactic latitudes do originate in the Galactic disk and not in the GC. Moreover, the ejection velocities of these extreme disk-runaway stars exceed the predicted limits of the classical scenarios for the production of runaway stars. Based on proper motions from the Gaia early data release 3 and on recent and new spectrophotometric distances, we studied the kinematics of 30 such extreme disk-runaway stars, allowing us to deduce their spatial origins in and their ejection velocities from the Galactic disk with unprecedented precision. Only three stars in the sample have past trajectories that are consistent with an origin in the GC, most notably S5-HVS1, which is the most extreme object in the sample by far. All other program stars are shown to be disk runaways with ejection velocities that sharply contrast at least with classical ejection scenarios. They include HVS5 and HVS6, which are both gravitationally unbound to the Milky Way. While most stars originate from within a galactocentric radius of 15kpc, which corresponds to the observed extent of the spiral arms, a group of five stars stems from radii of about 21-29kpc. This indicates a possible link to outer Galactic rings and a potential origin from infalling satellite galaxies.

Hydride molecules lie at the base of interstellar chemistry, but the synthesis of sulfuretted hydrides is poorly understood. Motivated by new observations of the Orion Bar PDR - 1'' resolution ALMA images of SH+; IRAM 30m detections of H2S, H2S34, and H2S33; H3S+ (upper limits); and SOFIA observations of SH - we perform a systematic study of the chemistry of S-bearing hydrides. We determine their column densities using coupled excitation, radiative transfer as well as chemical formation and destruction models. We revise some of the key gas-phase reactions that lead to their chemical synthesis. This includes ab initio quantum calculations of the vibrational-state-dependent reactions SH+ + H2 <-> H2S+ + H and S + H2 <-> SH + H. We find that reactions of UV-pumped H2 (v>1) with S+ explain the presence of SH+ in a high thermal-pressure gas component, P_th~10^8 cm^-3 K, close to the H2 dissociation front. However, subsequent hydrogen abstraction reactions of SH+, H2S+, and S with vibrationally excited H2, fail to ultimately explain the observed H2S column density (~2.5x10^14 cm^-2, with an ortho-to-para ratio of 2.9+/-0.3). To overcome these bottlenecks, we build PDR models that include a simple network of grain surface reactions leading to the formation of solid H2S (s-H2S). The higher adsorption binding energies of S and SH suggested by recent studies imply that S atoms adsorb on grains (and form s-H2S) at warmer dust temperatures and closer to the UV-illuminated edges of molecular clouds. Photodesorption and, to a lesser extent, chemical desorption, produce roughly the same H2S column density (a few 10^14 cm-^2) and abundance peak (a few 10^-8) nearly independently of n_H and G_0. This agrees with the observed H2S column density in the Orion Bar as well as at the edges of dark clouds without invoking substantial depletion of elemental sulfur abundances.

An analytical chemical evolution model is constructed to investigate the radial distribution of gas-phase and stellar metallicity for star-forming galaxies. By means of the model, the gas-phase and stellar metallicity can be obtained from the stellar-to-gas mass ratio. Both the gas inflow and outflow processes play an important role in building the final gas-phase metallicity, and there exists degeneracy effect between the gas inflow and outflow rates for star-forming galaxies. On the other hand, stellar metallicity is more sensitive to the gas outflow rate than to the gas inflow rate, and this helps to break the parameter degeneracy for star-forming galaxies. We apply this analysis method to the nearby disc galaxy M\,101 and adopting the classical ? 2 methodology to explore the influence of model parameters on the resulted metallicity. It can be found that the combination of gas-phase and stellar metallicity is indeed more effective for constraining the gas inflow and outflow rates. Our results also show that the model with relatively strong gas outflows but weak gas inflow describes the evolution of M\,101 reasonably well.

The morphological classification of radio sources is important to gain a full understanding of galaxy evolution processes and their relation with local environmental properties. Furthermore, the complex nature of the problem, its appeal for citizen scientists and the large data rates generated by existing and upcoming radio telescopes combine to make the morphological classification of radio sources an ideal test case for the application of machine learning techniques. One approach that has shown great promise recently is Convolutional Neural Networks (CNNs). Literature, however, lacks two major things when it comes to CNNs and radio galaxy morphological classification. Firstly, a proper analysis of whether overfitting occurs when training CNNs to perform radio galaxy morphological classification using a small curated training set is needed. Secondly, a good comparative study regarding the practical applicability of the CNN architectures in literature is required. Both of these shortcomings are addressed in this paper. Multiple performance metrics are used for the latter comparative study, such as inference time, model complexity, computational complexity and mean per class accuracy. As part of this study we also investigate the effect that receptive field, stride length and coverage has on recognition performance. For the sake of completeness, we also investigate the recognition performance gains that we can obtain by employing classification ensembles. A ranking system based upon recognition and computational performance is proposed. MCRGNet, Radio Galaxy Zoo and ConvXpress (novel classifier) are the architectures that best balance computational requirements with recognition performance.

The chemistry of the diffuse interstellar medium is driven by the combined influences of cosmic rays, ultraviolet (UV) radiation, and turbulence. Previously detected at the outer edges of photodissociation regions (PDRs) and formed from the reaction of C+ and OH, CO+ is the main chemical precursor of HCO+ and CO in a thermal, cosmic-ray, and UV-driven chemistry. Our aim was to test whether the thermal cosmic-ray and UV-driven chemistry is producing CO in diffuse interstellar molecular gas through the intermediate formation of CO+ We searched for CO+ absorption with the Atacama Large Millimeter Array (ALMA) toward two quasars with known Galactic foreground absorption from diffuse interstellar gas, J1717-3342 and J1744-3116, targeting the two strongest hyperfine components of the J=2-1 transition near 236 GHz. We could not detect CO+ but obtained sensitive upper limits toward both targets. The derived upper limits on the CO+ column densities represent about 4% of the HCO+ column densities. The corresponding upper limit on the CO+ abundance relative to H2 is <1.2 x 10^{-10}. The non-detection of CO+ confirms that HCO+ is mainly produced in the reaction between oxygen and carbon hydrides, CH2+ or CH3+ , induced by suprathermal processes, while CO+ and HOC+ result from reactions of C+ with OH and H2O. The densities required to form CO molecules at low extinction are consistent with this scheme.

Some luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs) host extremely compact and dusty nuclei. The intense infrared radiation arising from warm dust in these sources is prone to excite vibrational levels of molecules such as HCN. This results in emission from the rotational transitions of vibrationally excited HCN (HCN-vib), with the brightest emission found in compact obscured nuclei (CONs). We aim to establish how common CONs are in the local Universe, and whether their prevalence depends on the luminosity or other properties of the host galaxy. We have conducted an Atacama Large Millimeter/submillimeter Array (ALMA) survey of the rotational J=3-2 transition of HCN-vib in a sample of 46 far-infrared luminous galaxies. Compact obscured nuclei are identified in 38 percent of ULIRGs, 21 percent of LIRGs, and 0 percent of lower luminosity galaxies. We find no dependence on the inclination of the host galaxy, but strong evidence of lower IRAS 25 to 60 {\mu}m flux density ratios (f25/f60) in CONs compared to the rest of the sample. Furthermore, we find that CONs have stronger silicate features (s9.7{\mu}m) but similar PAH equivalent widths (EQW6.2{\mu}m) compared to other galaxies. In the local Universe, CONs are primarily found in (U)LIRGs. High resolution continuum observations of the individual nuclei are required to determine if the CON phenomenon is related to the inclinations of the nuclear disks. The lower f25/f60 ratios in CONs as well as the results for the mid-infrared diagnostics investigated are consistent with large dust columns shifting the nuclear radiation to longer wavelengths, making the mid- and far-infrared "photospheres" significantly cooler than the interior regions. To assess the importance of CONs in the context of galaxy evolution, it is necessary to extend this study to higher redshifts where (U)LIRGs are more common.

We propose to use relative strengths of far-infrared fine structure lines from galaxies to characterise early phases of the inside-out quenching by massive black holes (BHs). The BH feedback is thought to quench star formation by evacuating the ambient gas. In order to quantify the feedback effect on the gas density in the galactic centres, we utilise the outputs of IllustrisTNG and Illustris simulations, which implement different BH feedback models. We devise a physical model of H II regions and compute the intensities of [O III ] 52 and 88 μm lines. The line intensity ratio is sensitive to the local electron density, and thus can be used to measure the strength and physical extent of the BH quenching. If the BH feedback abruptly operates and expel the gas when it grows to a certain mass, as modelled in IllustrisTNG, the low-density gas yields relatively weak [O III ] 52 line with respect to 88 μm . In contrast, if the feedback strength and hence the local gas density are not strongly correlated with the BH mass, as in Illustris, the line ratio is not expected to vary significantly among galaxies with different evolutionary stages. We find these features are reproduced in the simulations. We also show that the line ratios are not sensitive to the aperture size for measurement, and thus observations do not need to resolve the galactic centres. We argue that the integrated line ratios can be used to capture the onset of the inside-out quenching by BHs.

We present an infrared (IR) characterization of the born-again planetary nebulae (PNe) A30 and A78 using IR images and spectra. We demonstrate that the carbon-rich dust in A30 and A78 is spatially coincident with the H-poor ejecta and coexists with hot X-ray-emitting gas up to distances of 50 ?��?from the central stars (CSPNs). Dust forms immediately after the born-again event and survives for 1000 yr in the harsh environment around the CSPN as it is destroyed and pushed away by radiation pressure and dragged by hydrodynamical effects. Spitzer IRS spectral maps showed that the broad spectral features at 6.4 and 8.0 μ m, attributed to amorphous carbon formed in H-deficient environments, are associated with the disrupted disk around their CSPN, providing an optimal environment for charge exchange reactions with the stellar wind that produces the soft X-ray emission of these sources. Nebular and dust properties are modeled for A30 with Cloudy taking into account different carbonaceous dust species. Our models predict dust temperatures in the 40-230 K range, five times lower than predicted by previous works. Gas and dust masses for the born-again ejecta in A30 are estimated to be M gas =( 4.41 +0.55 ??.14 )? 10 ?? M ??and M dust =( 3.20 +3.21 ??.06 )? 10 ?? M ??, which can be used to estimate a total ejected mass and mass-loss rate for the born-again event of ( 7.61 +3.76 ??.20 )? 10 ?? M ??and M ? =[5??0]? 10 ?? M ??yr ?? , respectively. Taking into account the carbon trapped into dust grains, we estimate that the C/O mass ratio of the H-poor ejecta of A30 is larger than 1, which favors the very late thermal pulse model over the alternate hypothesis of a nova-like event.

We have analyzed ALMA Cycle 5 data in Band 4 toward three low-mass young stellar objects (YSOs), IRAS 03235+3004 (hereafter IRAS 03235), IRAS 03245+3002 (IRAS 03245), and IRAS 03271+3013 (IRAS 03271), in the Perseus region. The HC 3 N ( J=16??5 ; E up /k=59.4 K) line has been detected in all of the target sources, while four CH 3 OH lines ( E up /k=15.4??6.3 K) have been detected only in IRAS 03245. Sizes of the HC 3 N distributions ( ??930??230 au) in IRAS 03235 and IRAS 03245 are similar to those of the carbon-chain species in the warm carbon chain chemistry (WCCC) source L1527. The size of the CH 3 OH emission in IRAS 03245 is ??760 au, which is slightly smaller than that of HC 3 N in this source. We compare the CH 3 OH/HC 3 N abundance ratios observed in these sources with predictions of chemical models. We confirm that the observed ratio in IRAS 03245 agrees with the modeled values at temperatures around 30--35 K, which supports the HC 3 N formation by the WCCC mechanism. In this temperature range, CH 3 OH does not thermally desorb from dust grains. Non-thermal desorption mechanisms or gas-phase formation of CH 3 OH seem to work efficiently around IRAS 03245. The fact that IRAS 03245 has the highest bolometric luminosity among the target sources seems to support these mechanisms, in particular the non-thermal desorption mechanisms.

Current models of galaxy evolution are constrained by the analysis of catalogs containing the flux and size of galaxies extracted from multiband deep fields carrying inevitable observational and extraction-related biases which can be highly correlated. In practice, taking all of these effects simultaneously into account is difficult, and derived models are inevitably biased. To address this issue, we use robust likelihood-free methods for the inference of luminosity function parameters, made possible via massive compression of multiband images using artificial neural networks. This technique makes the use of catalogs unnecessary when comparing observed and simulated multiband deep fields and constraining model parameters. A forward modeling approach generates galaxies of multiple types depending on luminosity function parameters and paints them on photometric multiband deep fields including both the instrumental and observational characteristics. The simulated and the observed images present the same selection effects and can therefore be properly compared. We train a fully-convolutional neural network to extract the most model-parameter-sensitive summary statistics out of these realistic simulations, shrinking down the dimensionality of the summary space. Finally, using the trained network to compress both observed and simulated deep fields, the model parameter values are constrained through Population Monte Carlo likelihood-free inference. Using synthetic photometric multiband deep fields similar to the CFHTLS and D1/D2 deep fields and massively compressing them through the convolutional neural network, we demonstrate the robustness, accuracy and consistency of this new catalog-free inference method. We are able to constrain the parameters of luminosity functions of different types of galaxies and our results are fully compatible with the classic catalog extraction approaches.