Search

We present Chandra X-ray observations of four optically-selected tidal disruption events (TDEs) obtained 4-9 years after discovery. Three sources were detected with luminosities between 9X10^40 and 3X10^42 erg/s. The spectrum of PTF09axc is consistent with a power law of index 2.5+-0.1, whereas the spectrum of PTF09ge is very soft. The power law spectrum of PTF09axc and prior literature findings, provide evidence that TDEs transition from an early-time soft state to a late-time hard state many years after disruption. We propose that the time to peak luminosity for optical and X-ray emission may differ substantially in TDEs, with X-rays being produced or becoming observable later. This delay helps explain the differences in observed properties such as L_opt/L_ X of optically and X-ray selected TDEs. We update TDE rate predictions for the eROSITA instrument: it ranges from 3 per yr to 990 per yr, depending sensitively on the distribution of black hole spins and the time delay between disruption and peak X-ray brightness. We further predict an asymmetry in the number of retrograde and prograde disks in samples of optically and X-ray selected TDEs. The details of the observational biases can contribute to observed differences between optically and X-ray selected TDEs (with optically selected TDEs being fainter in X-rays for retrograde TDE disks)., Comment: Accepted for publication to ApJ, 18 pages, 4 figures

High-energy gamma rays of interstellar origin are produced by the interaction of cosmic-ray (CR) particles with the diffuse gas and radiation fields in the Galaxy. The main features of this emission are well understood and are reproduced by existing CR propagation models employing 2D Galactocentric cylindrically symmetrical geometry. However, the high-quality data from instruments like the Fermi Large Area Telescope reveal significant deviations from the model predictions on few to tens of degree scales indicating the need to include the details of the Galactic spiral structure and thus require 3D spatial modelling. In this paper the high-energy interstellar emissions from the Galaxy are calculated using the new release of the GALPROP code employing 3D spatial models for the CR source and interstellar radiation field (ISRF) densities. Three models for the spatial distribution of CR sources are used that are differentiated by their relative proportion of input luminosity attributed to the smooth disc or spiral arms. Two ISRF models are developed based on stellar and dust spatial density distributions taken from the literature that reproduce local near- to far-infrared observations. The interstellar emission models that include arms and bulges for the CR source and ISRF densities provide plausible physical interpretations for features found in the residual maps from high-energy gamma-ray data analysis. The 3D models for CR and ISRF densities provide a more realistic basis that can be used for the interpretation of the non-thermal interstellar emissions from the Galaxy., 23 pages, 14 figures, 1 appendix. ApJ in press

We quantify the luminosity contribution of active galactic nuclei (AGN) to the 12 $\mu$m, mid-infrared (MIR; 5-38 $\mu$m), and the total IR (5-1000 $\mu$m) emission in the local AGN detected in the all-sky 70-month Swift/Burst Alert Telescope (BAT) ultra hard X-ray survey. We decompose the IR spectral energy distributions (SEDs) of 587 objects into AGN and starburst components using AGN torus and star-forming galaxy templates. This enables us to recover the AGN torus emission also for low-luminosity end, down to $\log (L_{14-150}/{\rm erg}~{\rm s}^{-1}) \simeq 41$, which typically have significant host galaxy contamination. We find that the luminosity contribution of the AGN to the 12 $\mu$m, the MIR, and the total IR band is an increasing function of the 14-150 keV luminosity. We also find that for the most extreme cases, the IR pure-AGN emission from the torus can extend up to 90 $\mu$m. The obtained total IR AGN luminosity through the IR SED decomposition enables us to estimate the fraction of the sky obscured by dust, i.e., the dust covering factor. We demonstrate that the median of the dust covering factor is always smaller than that of the X-ray obscuration fraction above the AGN bolometric luminosity of $\log (L_{\rm bol}/{\rm erg}~{\rm s}^{-1}) \simeq 42.5$. Considering that X-ray obscuration fraction is equivalent to the covering factor coming from both the dust and gas, it indicates that an additional neutral gas component, along with the dusty torus, is responsible for the absorption of X-ray emission., Comment: 21 pages, 15 figures, accepted for publication in ApJ. The full list of Table 1 is available at http://www.kusastro.kyoto-u.ac.jp/~ichikawa/Table1_MR_20181107.txt

Core-collapse supernovae (SNe) expand into a medium created by winds from the pre-SN progenitor. The SN explosion and resulting shock wave(s) heat up the surrounding plasma, giving rise to thermal X-ray emission, which depends on the density of the emitting material. Tracking the variation of the X-ray luminosity over long periods of time thus allows for investigation of the kinematics of the SN shock waves, the structure of the surrounding medium, and the nature of the progenitor star. In this paper X-ray observations of five of the oldest known X-ray supernovae - SN 1970G, SN 1968D, SN 1959D, SN 1957D and SN 1941C - are analyzed, with the aim of reconstructing their light curves over several decades. For those supernovae for which we can extract multi-epoch data, the X-ray luminosity appears to decline with time, although with large error bars. No increase in the X-ray emission from SN 1970G is found at later epochs, contrary to previous reports. All five SNe show X-ray luminosities that are of comparable magnitude. We compare the late-time X-ray luminosities of these SNe to those of supernova remnants (SNRs) in the Galaxy which are a few hundred years old, and find that when the tentative decline is taken into account, the luminosity of the old SNe studied herein could fall below the luminosity of some of the younger SNRs within a few hundred years. However, the X-ray luminosity should begin to increase as the SNe expand in the Sedov phase, thus reaching that of the observed SNRs., Accepted to the Astrophysical Journal. 24 pages, 12 tables, 16 figures

We present a detailed study of emission-line systems in the GAMA G23 region, making use of $\textit{WISE}$ photometry that includes carefully measured resolved sources. After applying several cuts to the initial catalogue of $\sim$41,000 galaxies, we extract a sample of 9,809 galaxies. We then compare the spectral diagnostic (BPT) classification of 1154 emission-line galaxies (38$\%$ resolved in W1) to their location in the $\textit{WISE}$ colour-colour diagram, leading to the creation of a new zone for mid-infrared "warm" galaxies located 2$\sigma$ above the star-forming sequence, below the standard $\textit{WISE}$ AGN region. We find that the BPT and $\textit{WISE}$ diagrams agree on the classification for 85$\%$ and 8$\%$ of the galaxies as non-AGN (star forming = SF) and AGN, respectively, and disagree on $\sim$7$\%$ of the entire classified sample. 39$\%$ of the AGN (all types) are broad-line systems for which the [\ion{N}{ii}] and [H$\alpha$] fluxes can barely be disentangled, giving in most cases spurious [\ion{N}{ii}]/[H$\alpha$] flux ratios. However, several optical AGN appear to be completely consistent with SF in $\textit{WISE}$. We argue that these could be low power AGN, or systems whose hosts dominate the IR emission. Alternatively, given the sometimes high [\ion{O}{iii}] luminosity in these galaxies, the emission lines may be generated by shocks coming from super-winds associated with SF rather than the AGN activity. Based on our findings, we have created a new diagnostic: [W1-W2] vs [\ion{N}{ii}]/[H$\alpha$], which has the virtue of separating SF from AGN and high-excitation sources. It classifies 3$\sim$5 times more galaxies than the classic BPT, Comment: 43 pages, 32 figures, 4 tables

The size and structure of the dusty circumnuclear torus in active galactic nuclei (AGN) can be investigated by analyzing the temporal response of the torus's infrared (IR) dust emission to variations in the AGN ultraviolet/optical luminosity. This method, reverberation mapping, is applicable over a wide redshift range, but the IR response is sensitive to several poorly constrained variables relating to the dust distribution and its illumination, complicating the interpretation of measured reverberation lags. We have used an enhanced version of our torus reverberation mapping code (TORMAC) to conduct a comprehensive exploration of the torus response functions at selected wavelengths, for the standard interstellar medium grain composition. The shapes of the response functions vary widely over the parameter range covered by our models, with the largest variations occurring at shorter wavelengths ($\leq 4.5\,\mu$m). The reverberation lag, quantified as the response-weighted delay (RWD), is most affected by the radial depth of the torus, the steepness of the radial cloud distribution, the degree of anisotropy of the AGN radiation field, and the volume filling factor. Nevertheless, we find that the RWD provides a reasonably robust estimate, to within a factor of $\sim 3$, of the luminosity-weighted torus radius, confirming the basic assumption underlying reverberation mapping. However, overall, the models predict radii at $2.2\,\mu$m that are a typically factor of $\sim 2$ larger than those derived from K-band reverberation mapping. This is likely an indication that the innermost region of the torus is populated by clouds dominated by large graphite grains., Comment: 38 pages, 36 figures (13 in appendices), published in the Astrophysical Journal

Extreme outbursts in young stars may be a common stage of pre-main-sequence stellar evolution. These outbursts, caused by enhanced accretion and accompanied by increased luminosity, can also strongly impact the evolution of the circumstellar environment. We present ALMA and VLA observations of EX Lupi, a prototypical outburst system, at 100 GHz, 45 GHz, and 15 GHz. We use these data, along with archival ALMA 232 GHz data, to fit radiative transfer models to EX Lupi's circumstellar disk in its quiescent state following the extreme outburst in 2008. The best fit models show a compact disk with a characteristic dust radius of 45 au and a total mass of 0.01 M$_{\odot}$. Our modeling suggests grain growth to sizes of at least 3 mm in the disk, possibly spurred by the recent outburst, and an ice line that has migrated inward to $0.2-0.3$ au post-outburst. At 15 GHz, we detected significant emission over the expected thermal disk emission which we attribute primarily to stellar (gyro)synchrotron and free-free disk emission. Altogether, these results highlight what may be a common impact of outbursts on the circumstellar dust., Accepted to ApJ, 15 pages, 8 figures

We present optical photometry and spectroscopy of the Type Ia supernova SN2018cqj/ATLAS18qtd. The supernova exploded in an isolated region at $\sim 65$~kpc from the S0 galaxy IC~550 at $z=0.0165$ ($D\approx 74$~Mpc) and has a redshift consistent with a physical association to this galaxy. Multicolor photometry show that SN2018cqj/ATLAS18qtd is a low-luminosity ($M_{B_{max}}\approx -17.9$ mag), fast-declining Type Ia with color stretch $s_{BV} \approx 0.6$ and $B$-band decline rate $\Delta m_{15}(B) \approx 1.77$ mag. Two nebular-phase spectra obtained as part of the 100IAS survey at +193 and +307 days after peak show the clear detection of a narrow H$\alpha$ line in emission that is resolved in the first spectrum with $\rm FWHM \approx 1200$ km s$^{-1}$ and $L_{H\alpha} \approx 3.8\times 10^{37}$ erg s$^{-1}$. The detection of a resolved H$\alpha$ line with a declining luminosity is broadly consistent with recent models where hydrogen is stripped from the non-degenerate companion in a single-degenerate progenitor system. However, the amount of hydrogen consistent with the luminosities of the H$\alpha$ line would be $\sim 10^{-3}$ M$_{\odot}$, significantly less than theoretical model predictions in the classical single-degenerate progenitor systems. SN2018cqj/ATLAS18qtd is the second low-luminosity, fast-declining Type Ia SN after SN2018fhw/ASASSN-18tb that shows narrow H$\alpha$ in emission in its nebular-phase spectra., Comment: Accepted for publication in ApJ

We present the results of a detailed investigation of the prompt and afterglow emission in the HESS detected GRB 190829A. Swift and Fermi observations of the prompt phase of this GRB reveal two isolated sub-bursts or episodes, separated by a quiescent phase. The energetic and the spectral properties of the first episode are in stark contrast to the second. The first episode, which has a higher spectral peak of $\sim 120\:\text{keV}$ and a low isotropic energy $\sim 10^{50}\:\text{erg}$ is an outlier to the Amati correlation and marginally satisfies the Yonetoku correlation. However, the energetically dominant second episode has lower peak energy and is consistent with the above correlations. We compared this GRB to other low luminosity GRBs (LLGRBs). Prompt emission of LLGRBs also indicates a relativistic shock breakout origin of the radiation. For GRB 190829A, some of the properties of a shock breakout origin are satisfied. However, the absence of an accompanying thermal component and energy above the shock breakout critical limit precludes a shock breakout origin. In the afterglow, an unusual long-lasting late time flare of duration $\sim 10^4\:\text{s}$ is observed. We also analyzed the late-time \fermi-LAT emission that encapsulates the H.E.S.S. detection. Some of the LAT photons are likely to be associated with the source. All the above observational facts suggest GRB 190829A is a peculiar low luminosity GRB that is not powered by a shock breakout, and with an unusual rebrightening due to a patchy emission or a refreshed shock during the afterglow. Furthermore, our results show that TeV energy photons seem common in both high luminosity GRBs and LLGRBs., Comment: Accepted for publication in ApJ. The shock breakout origin of the prompt emission and central engine activity for the flare are revised with more counter-evidences. Major changes are marked in bold-face

We present optical and near-infrared ($ugriYJH$) photometry of host galaxies of Type Ia supernovae (SN~Ia) observed by the \textit{Carnegie Supernova Project-I}. We determine host galaxy stellar masses and, for the first time, study their correlation with SN~Ia standardized luminosity across optical and near-infrared ($uBgVriYJH$) bands. In the individual bands, we find that SNe~Ia are more luminous in more massive hosts with luminosity offsets ranging between $-0.07 \pm0.03$ mag to $-0.15\pm0.04$ mag after light-curve standardization. The slope of the SN~Ia Hubble residual-host mass relation is negative across all $uBgVriYJH$ bands with values ranging between $-0.036\pm 0.025$ mag/dex to $-0.097\pm 0.027$ mag/dex -- implying that SNe~Ia in more massive galaxies are brighter than expected. The near-constant observed correlations across optical and near-infrared bands indicate that dust may not play a significant role in the observed luminosity offset--host mass correlation. We measure projected separations between SNe~Ia and their host centers, and find that SNe~Ia that explode beyond a projected 10 kpc have a $\rm 30\% \ to \ 50\%$ reduction of the dispersion in Hubble residuals across all bands -- making them a more uniform subset of SNe~Ia. Dust in host galaxies, peculiar velocities of nearby SN~Ia, or a combination of both may drive this result as the color excesses of SNe~Ia beyond 10 kpc are found to be generally lower than those interior, but there is also a diminishing trend of the dispersion as we exclude nearby events. We do not find that SN~Ia average luminosity varies significantly when they are grouped in various host morphological types. Host galaxy data from this work will be useful, in conjunction with future high-redshift samples, in constraining cosmological parameters., Accepted to The Astrophysical Journal

We present CO(J=1-0) and/or CO(J=2-1) spectroscopy for 31 galaxies selected from the ongoing MaNGA survey, obtained with multiple telescopes. This sample is combined with CO observations from the literature to study the correlation of the CO luminosities ($L_{\rm CO(1-0)}$) with the mid-infrared luminosities at 12 ($L_{12 \mu m}$) and 22 $\mu$m ($L_{\rm 22 \mu m}$), as well as the dependence of the residuals on a variety of galaxy properties. The correlation with $L_{\rm 12 \mu m}$ is tighter and more linear, but galaxies with relatively low stellar masses and blue colors fall significantly below the mean $L_{\rm CO(1-0)}-L_{\rm 12\mu m}$ relation. We propose a new estimator of the CO(1-0) luminosity (and thus the total molecular gas mass) that is a linear combination of three parameters: $L_{\rm 12 \mu m}$, $M_\ast$ and $g-r$. We show that, with a scatter of only 0.18 dex in log $(L_{\rm CO(1-0)})$, this estimator provides unbiased estimates for galaxies of different properties and types. An immediate application of this estimator to a compiled sample of galaxies with only CO(J=2-1) observations yields a distribution of the CO(J=2-1) to CO(J=1-0) luminosity ratios ($R21$) that agrees well with the distribution of real observations, in terms of both the median and the shape. Application of our estimator to the current MaNGA sample reveals a gas-poor population of galaxies that are predominantly early-type and show no correlation between molecular gas-to-stellar mass ratio and star formation rate, in contrast to gas-rich galaxies. We also provide alternative estimators with similar scatters, based on $r$ and/or $z$ band luminosities instead of $M_\ast$. These estimators serve as cheap and convenient $M_{\rm mol}$ proxies to be potentially applied to large samples of galaxies, thus allowing statistical studies of gas-related processes of galaxies., Comment: 24 pages, 12 figures, accepted for publication in ApJ

We present a catalog of physical properties for galaxies hosting active galactic nuclei (AGN) detected by the Wide-field Infrared Survey Explorer (WISE). By fitting broadband spectral energy distributions of sources in the WISE AGN Catalog (Assef et al. 2018) with empirical galaxy and AGN templates, we derive photometric redshifts, AGN bolometric luminosities, measures of AGN obscuration, host galaxy stellar masses, and host galaxy star formation rates (SFRs) for 695,273 WISE AGN. The wide-area nature of this catalog significantly augments the known number of obscured AGN out to redshifts of z~3 and will be useful for studies focused on AGN or their host galaxy physical properties. We first show that the most likely non-AGN contaminants are galaxies at redshifts of z=0.2-0.3, with relatively blue W1-W2 colors, and with high specific SFRs for which the dust continuum emission is elevated in the W2 filter. Toward increasingly lower redshifts, WISE AGN host galaxies have systematically lower specific SFRs, relative to those of normal star forming galaxies, likely due to decreased cold gas fractions and the time delay between global star formation and AGN triggering. Finally, WISE AGN obscuration is not strongly correlated with AGN bolometric luminosity but shows a significant negative correlation with Eddington ratio. This result is consistent with a version of the `receding torus' model in which the obscuring material is located within the supermassive black hole gravitational sphere of influence and the dust inner radius increases due to radiation pressure., 21 pages, 17 figures. Accepted for publication in the Astrophysical Journal. The full catalog is available from the publisher or from the corresponding author upon request

We report the discovery of a massive protostar M17~MIR embedded in a hot molecular core in M17. The multiwavelength data obtained during 1993--2019 show significant mid-IR (MIR) variations, which can be split into three stages: the decreasing phase during 1993.03--mid-2004, the quiescent phase from mid-2004 to mid-2010, and the rebrightening phase from mid-2010 until now. The variation of the 22\,GHz H$_2$O maser emission, together with the MIR variation, indicates an enhanced disk accretion rate onto M17~MIR during the decreasing and rebrightening phases. Radiative transfer modeling of the spectral energy distributions of M17~MIR in the 2005 epoch (quiescent) and 2017 epoch (accretion outburst) constrains the basic stellar parameters of M17~MIR, which is an intermediate-mass protostar (M~5.4 Msun) with accretion rate ~1.1x10^-5 Msun in the 2005 epoch and ~1.7x10^-3 Msun/yr in the 2017 epoch. The enhanced accretion rate during outburst induces the luminosity outburst $\Delta L\approx7600 $Lsun. In the accretion outburst, a larger stellar radius is required to produce accretion rate consistent with the value estimated from the kinematics of water masers. M17 MIR shows two accretion outbursts ($\Delta t\sim 9-20$ yr) with outburst magnitudes of 2 mag, separated by a 6 yr quiescent phase. The accretion outbusrt occupies 83\% of the time over 26 yr. The accretion rate in outburst is variable with amplitude much lower than the contrast between quiescent and outburst phases. The extreme youth of M17 MIR suggests that minor accretion bursts are frequent in the earliest stages of massive star formation., Comment: Typos corrected. The version most close to the publication

We consider whether the subset of carbon-rich asymptotic giant branch (AGB) stars that exhibit detached, expanding circumstellar shells may reveal the past histories of these stars as having undergone helium shell flashes (thermal pulses) on the AGB. We exploit newly available Gaia parallaxes and photometry, along with archival infrared photometry, to obtain refined estimates of the luminosities of all (12) known detached shell carbon stars. We examine the relationship between these luminosities and the estimated dynamical ages (ejection times) of the detached shells associated with the 12 stars, which range from $\sim$1000 to $\sim$30000 yr. When arranged according to detached shell dynamical age, the (implied) luminosity evolution of the known detached shell carbon stars closely follows the predicted "light curves" of individual thermal pulses obtained from models of AGB stars. The comparison between data and models suggests that detached shell carbon stars are descended from $\sim$2.5-4.0 $M_\odot$ progenitors. We conclude that detached shell carbon stars may serve as effective tracers of the luminosity evolution of AGB thermal pulses., Comment: 14 pages, 4 figures; accepted for publication in The Astrophysical Journal

Warm absorbers (WAs) located approximately in the region of $1-1000$ parsecs are common phenomena in many active galactic nuclei (AGNs). The driving mechanism of WAs is still under debate. Ultra-fast outflows (UFOs) which are launched very close to the central black hole are also frequently observed in AGNs. When UFOs move outwards, they will collide with the interstellar medius (ISM) gas. In this paper, we study the possibility that whether WAs can be generated by the interaction between ISM gas and the UFOs. We find that under some ISM gas conditions, WAs can be generated. However, the covering factor of WAs is much smaller than that given by observations. This indicates that other mechanisms should also be at work. We also find that the properties of the WAs mainly depend on the density of the ISM injected into the computational domain from the outer radial boundary (1000 parsec). The higher the density of the ISM is, the higher the mass flux and kinetic power of the WAs will be. The kinetic power of the UFO driven WAs is much less than $1\%$ of the bolometric luminosity of its host AGNs. Therefore, the UFO driven WAs might not contribute sufficient feedback to its host galaxy., Comment: 10 pages, 6 figures, published in ApJ

Millimeter continuum imaging of protoplanetary disks reveals the distribution of solid particles and the presence of substructures (gaps and rings) beyond 5-10 au, while infrared (IR) spectra provide access to abundances of gaseous species at smaller disk radii. Building on recent observational findings of an anti-correlation between the inner disk water luminosity and outer dust disk radius, we aim here at investigating the dynamics of icy solids that drift from the outer disk and sublimate their ice inside the snow line, enriching the water vapor that is observed in the IR. We use a volatile-inclusive disk evolution model to explore a range of conditions (gap location, particle size, disk mass, and alpha-viscosity) under which gaps in the outer disk efficiently block the inward drift of icy solids. We find that inner-disk vapor enrichment is highly sensitive to the location of a disk gap, yielding for each particle size a radial "sweet spot" that reduces the inner-disk vapor enrichment to a minimum. For pebbles of 1-10 mm in size, which carry the most mass, this sweet spot is at 7-15 au, suggesting that inner gaps may have a key role in reducing ice delivery to the inner disk and may not allow the formation of Earths and super-Earths. This highlights the importance of observationally determining the presence and properties of inner gaps in disks. Finally, we argue that the inner water vapor abundance can be used as a proxy for estimating the pebble drift efficiency and mass-flux entering the inner disk., Comment: 24 pages, 14 figures; Accepted for publication in ApJ

For over 60 years, the scientific community has studied actively growing central super-massive black holes (active galactic nuclei -- AGN) but fundamental questions on their genesis remain unanswered. Numerical simulations and theoretical arguments show that black hole growth occurs during short-lived periods ($\sim$ 10$^{7}$ -10$^{8}$ yr) of powerful accretion. Major mergers are commonly invoked as the most likely dissipative process to trigger the rapid fueling of AGN. If the AGN-merger paradigm is true, we expect galaxy mergers to coincide with black hole accretion during a heavily obscured AGN phase (N$_H$ $ > 10^{23}$ cm$^{-2}$). Starting from one of the largest samples of obscured AGN at 0.5 $, Comment: 19 pages, 11 figures, accepted in ApJ

We present a systematic study of the diffuse hot gas around early-type galaxies (ETGs) residing in the Virgo cluster, based on archival {\it Chandra} observations. Our representative sample consists of 79 galaxies with low-to-intermediate stellar masses ($M_* \approx 10^{9-11}\rm~M_\odot$), a mass range that has not been extensively explored with X-ray observations thus far. We detect diffuse X-ray emission in only eight galaxies and find that in five cases a substantial fraction of the detected emission can be unambiguously attributed to truly diffuse hot gas, based on their spatial distribution and spectral properties. For the individually non-detected galaxies, we constrain their average X-ray emission by performing a stacking analysis, finding a specific X-ray luminosity of $L_{\rm X}/M_* \sim 10^{28}{\rm~erg~s^{-1}~M_{\odot}^{-1}}$, which is consistent with unresolved stellar populations. The apparent paucity of truly diffuse hot gas in these low- and intermediate-mass ETGs may be the result of efficient ram pressure stripping by the hot intra-cluster medium. However, we also find no significant diffuse hot gas in a comparison sample of 57 field ETGs of similar stellar masses, for which archival {\it Chandra} observations with similar sensitivity are available. This points to the alternative possibility that galactic winds evacuate the hot gas from the inner region of low- and intermediate-mass ETGs, regardless of the galactic environment. Nevertheless, we do find strong morphological evidence for on-going ram pressure stripping in two galaxies (NGC 4417 and NGC 4459). A better understanding of the roles of ram pressure stripping and galactic winds in regulating the hot gas content of ETGs, invites sensitive X-ray observations for a large galaxy sample., 25 pages, 7 figures. Accepted for publication in ApJ. Comments Welcome

The detection of TeV gamma-ray bursts (GRBs) brought new opportunities for studying the physics of particle acceleration at relativistic shocks. The \hess telescopes recently observed very-high-energy (VHE) emission from a nearby low-luminosity GRB, GRB 190829A. Follow-up observations with, e.g., Swift-XRT, revealed unusual flare activities at $\sim 10^3\rm~s$, which can be caused by a long-lasting central engine. We show that the VHE emission during the H.E.S.S. observation time is naturally produced in the external inverse-Compton (EIC) scenario, where seed photons supplied by the flares or other late-time dissipation are upscattered to VHE energies by the non-thermal electrons accelerated at the external forward shock. Our calculations show that the EIC flare nearly coincides with the late-prompt flare, but extends $\sim3-4$ times longer than the duration of the late-prompt flare. The preferred kinetic energy and initial Lorentz factor used in our model are $\sim 10^{52}\rm~erg$ and $\sim 20$, respectively. Understanding the mechanisms of the VHE emission from low-luminosity GRBs will help us constrain the properties of the outflow and the central engine activities, as well as the particle acceleration mechanism., Comment: 12 pages, 6 figures, accepted for publication in ApJ. HESS data are added, figures are updated, appendix is extended, main conclusions unchanged

We present the mid-infrared (mid-IR) period$-$luminosity relations (PLRs) using over 1000 Mira variables in the Large Magellanic Cloud (LMC), for the four Wide-Field Infrared Survey Explorer (WISE) and the four Spitzer bands. These PLRs cover a mid-IR wavelength range from $3.4$ $\mu$m to $22$ $\mu$m and are presented separately for the oxygen-rich (O-rich) and carbon-rich (C-rich) Miras. These relations can be used to measure distances to individual O-rich and/or C-rich Mira stars with the accuracy of $5\%$ and $12\%$, respectively. They are the most accurate Mira PLRs in the mid-IR to date., Comment: Accepted for publication in The Astrophysical Journal. 18 pages, 8 figures, 4 tables. The full machine readable Table 3 will be available in the online journal just after publication

We present 0\farcs{14}-resolution Atacama Large Millimeter/submillimeter Array (ALMA) CO(2$-$1) observations of the circumnuclear gas disk in UGC 2698, a local compact galaxy. The disk exhibits regular rotation with projected velocities rising to 450 km s$^{-1}$ near the galaxy center. We fit gas-dynamical models to the ALMA data cube, assuming the CO emission originates from a dynamically cold, thin disk, and measured the mass of the supermassive black hole (BH) in UGC 2698 to be $M_{\mathrm{BH}} = (2.46 \pm{0.07}$ [$1\sigma$ stat] $^{+0.70}_{-0.78}$ [sys])$\times 10^9$ $M_\odot$. UGC 2698 is part of a sample of nearby early-type galaxies that are plausible $z\sim2$ red nugget relics. Previous stellar-dynamical modeling for three galaxies in the sample found BH masses consistent with the BH mass$-$stellar velocity dispersion ($M_{\mathrm{BH}}-\sigma_\star$) relation but over-massive relative to the BH mass$-$bulge luminosity ($M_{\mathrm{BH}}-L_{\mathrm{bul}}$) correlation, suggesting that BHs may gain the majority of their mass before their host galaxies. However, UGC 2698 is consistent with both $M_{\mathrm{BH}}-\sigma_\star$ and $M_{\mathrm{BH}}-L_{\mathrm{bul}}$. As UGC 2698 has the largest stellar mass and effective radius in the local compact galaxy sample, it may have undergone more recent mergers that brought it in line with the BH scaling relations. Alternatively, given that the three previously-measured compact galaxies are outliers from $M_{\mathrm{BH}}-L_{\mathrm{bul}}$, while UGC 2698 is not, there may be significant scatter at the poorly sampled high-mass end of the relation. Additional gas-dynamical $M_{\mathrm{BH}}$ measurements for the compact galaxy sample will improve our understanding of BH$-$galaxy co-evolution., Comment: 18 pages, 7 figures. Submitted

We present extensive optical photometric and spectroscopic observations of the high-velocity (HV) Type Ia supernova (SN Ia) 2017fgc, covering the phase from $\sim$ 12 d before to $\sim 389$ d after maximum brightness. SN 2017fgc is similar to normal SNe Ia, with an absolute peak magnitude of $M_{\rm max}^{B} \approx$ $-19.32 \pm 0.13$ mag and a post-peak decline of ${\Delta}m_{15}(B)$ = $1.05 \pm 0.07$ mag. Its peak bolometric luminosity is derived as $1.32 \pm 0.13) \times 10^{43} $erg s$^{-1}$, corresponding to a $^{56}$Ni mass of $0.51 \pm 0.03 M_{\odot}$. The light curves of SN 2017fgc are found to exhibit excess emission in the $UBV$ bands in the early nebular phase and pronounced secondary shoulder/maximum features in the $RrIi$ bands. Its spectral evolution is similar to that of HV SNe Ia, with a maximum-light Si II velocity of $15,000 \pm 150 $km s$^{-1}$ and a post-peak velocity gradient of $\sim$ $120 \pm 10 $km s$^{-1} $d$^{-1}$. The Fe II and Mg II lines blended near 4300 {\AA} and the Fe II, Si II, and Fe III lines blended near 4800 {\AA} are obviously stronger than those of normal SNe Ia. Inspecting a large sample reveals that the strength of the two blends in the spectra, and the secondary peak in the $i/r$-band light curves, are found to be positively correlated with the maximum-light Si II velocity. Such correlations indicate that HV SNe~Ia may experience more complete burning in the ejecta and/or that their progenitors have higher metallicity. Examining the birthplace environment of SN 2017fgc suggests that it likely arose from a stellar environment with young and high-metallicity populations., Comment: 30 pages, 21 figures, 5 tables

Young stellar objects are observed to have large X-ray fluxes and are thought to produce commensurate luminosities in energetic particles (cosmic rays). This particle radiation, in turn, can synthesize short-lived radioactive nuclei through spallation. With a focus on $^{26}$Al, this paper estimates the expected abundances of radioactive nulcei produced by spallation during the epoch of planet formation. In this model, cosmic rays are accelerated near the inner truncation radii of circumstellar disks, $r_{\scriptstyle X}\approx0.1$ AU, where intense magnetic activity takes place. For planets forming in this region, radioactive abundances can be enhanced over the values inferred for the early solar system (from meteoritic measurements) by factors of $\sim10-20$. These short-lived radioactive nuclei influence the process of planet formation and the properties of planets in several ways. The minimum size required for planetesimals to become fully molten decreases with increasing levels of radioactive enrichment, and such melting leads to loss of volatile components including water. Planets produced with an enhanced radioactive inventory have significant internal luminosity which can be comparable to that provided by the host star; this additional heating affects both atmospheric mass loss and chemical composition. Finally, the habitable zone of red dwarf stars is coincident with the magnetic reconnection region, so that planets forming at those locations will experience maximum exposure to particle radiation, and subsequent depletion of volatiles., Comment: 29 pages, 4 figures, accepted to ApJ

A newly born millisecond magnetar is thought to be the central engine of some gamma-ray bursts (GRBs), especially those that present long-lasting X-ray plateau emissions. By solving the field equations, we find that when the rotational speed of the magnetar is approaching the breakup limit, its radius $R$ and moment of inertia $I$ would undergo an obvious evolution as the magnetar spins down. Meanwhile, the values of $R$ and $I$ would sensitively depend on the adoption of neutron star (NS) equation of state (EoS) and the NS baryonic mass. With different EoSs and baryonic masses considered, the magnetic dipole radiation luminosity ($L_{\rm dip}$) could be variant within one to two orders of magnitude. We thus suggest that when using the X-ray plateau data of GRBs to diagnose the properties of the nascent NSs, EoS and NS mass information should be invoked as simultaneously constrained parameters. On the other hand, due to the evolution of $R$ and $I$, the temporal behavior of $L_{\rm dip}$ would become more complicated. For instance, if the spin-down process is dominated by gravitational wave emission due to the NS asymmetry caused by magnetic field distortion ($\epsilon\propto B_{p}^{2}$), the segment $L_{\rm dip}\propto t^{0}$ could be followed by $L_{\rm dip}\propto t^{-\gamma}$ with $\gamma$ larger than 3. This case could naturally interpret the so-called internal X-ray plateau feature shown in some GRB afterglows, which means the sharp decay following the plateau is unnecessarily corresponding to the NS collapsing. This may explain why some internal X-ray plateaus are followed by late time central engine activity, manifested through flares and second shallow plateaus.

PKS\,1351$-$018 at a redshift of $z=3.71$ is one of the most luminous, steady synchrotron sources with a luminosity $> 10^{47}$\,erg~s$^{-1}$. The synchrotron luminosity does not seem to vary by more than $\sim 25\%$ over 35 years. In order to appreciate this remarkable behavior, if it were at $z=0.5$, it would have a flux density at 15 GHz in a range of $110 - 137$\,Jy over 11 yrs. In spite of this steady behavior, two strong $\gamma$-ray flares $\lesssim 10^{49}$\,erg~s$^{-1}$ were detected in 2011 and 2016. There is blazar-like behavior coexisting with the steady behavior. This study is aimed at elucidating the dual nature of this source. We find that the radio source is extremely compact with a bright core and a steep spectrum secondary component, 12\,mas away, that appears to be constant in position and flux density in six epochs from 1995 to 2018. We estimate that a jet with a time averaged power of $(5.2 \pm 3.2) \times 10^{45}$\,erg~s$^{-1}$ terminates in this lobe that is advancing $\gtrsim 0.9 c$ at a deprojected distance of $1-3$\,kpc from the central engine. This is the rare case of a young ($\sim 6000$\,yr), very powerful radio source that is viewed a few degrees from the jet axis. We find evidence of a high velocity (4000\,km~s$^{-1}$), high ionization wind emanating form a luminous quasar. The young radio jet appears to experience modest bending as it navigates through the intense quasar environment., Comment: To appear in ApJ

We present new observations of the odd $z=0.96$ weak-line quasar PG1407+265, and report the discovery of CXOU J140927.9+261813, a $z=0.68$ X-ray cluster. Archival X-ray photometry spanning nearly four decades reveals that PG1407+265 is variable at the 1 dex level on a timescale of years. V-band variability is present with an amplitude less than 0.1 mag. The emission-line properties of PG1407+265 also reveal clear evidence for a powerful inflow or outflow due to near- or super-Eddington accretion, having a mechanical luminosity of order $10^{48}$ erg s$^{-1}$. Our follow-up {\sl Chandra} exposure centered on this object reveal a foreground $z=0.68$ cluster roughly 1' x 1'.5 in extent, offset to the east of PG1407+265, roughly coincident with the $z=0.68$ radio galaxy FIRST J140927.8+261818. This non-cool-core cluster contributes about 10\% of the X-ray flux of PG1407+265, has a mass of $(0.6- 5.5)\times10^{14} M_\odot$, and an X-ray gas temperature of ($2.2-4.3$) keV. Because the projected position of the quasar lies at about twice that of the cluster's inferred Einstein radius, lensing by the cluster is unlikely to explain the quasar's unusual properties. We also discuss the evidence for a second cluster centered on and at the redshift of the quasar., Comment: 15 pages, 11 figures

Very-high-energy (VHE) emission is usually interpreted in the synchrotron-self Compton (SSC) scenario, and expected from the low-redshift and high-luminosity gamma-ray bursts (GRBs), as GRB 180720B and GRB 190114C. Recently, VHE emission was detected by the H.E.S.S. telescopes from one of the closest burst GRB 190829A which was associated with the supernova (SN) 2019oyw. In this paper, we present a temporal and spectral analysis from optical bands to Fermi-LAT energy range over multiple observational periods beginning just after the BAT trigger time and extending for almost three months. We show that the X-ray and optical observations are consistent with synchrotron forward-shock emission evolving between the characteristic and cooling spectral breaks during the early and late afterglow in a uniform-density medium. Modeling the light curves together with its spectral energy distribution, it is shown that the outflow expands with an initial bulk Lorentz factor of $\Gamma\sim 30$, which is high for a low-luminosity GRBs and low for a high-luminosity GRBs. The values of the initial bulk Lorentz factor and the isotropic equivalent energy suggest that GRB 190829A is classified as an intermediate-luminosity burst and consequently, it becomes the first burst of this class in being detected in the VHE gamma-ray band by an imaging atmospheric Cherenkov telescope, and, in turn, the first event without being simultaneously observed by the Fermi-LAT instrument. Analyzing the intermediate-luminosity bursts with $z\lesssim 0.2$ such as GRB 130702A, we show that bursts with intermediate luminosity are potential candidates to be detected in very-high energies., Comment: 26 pages, 7 figures, Submitted to ApJ

SN 2010jl is a Type IIn core collapse supernova whose radiative output is powered by the interaction of the SN shock wave with its surrounding dense circumstellar medium (CSM). After day ~60, its light curve developed a NIR excess emission from dust. This excess could be a thermal IR echo from pre-existing CSM dust, or emission from newly-formed dust either in the cooling postshock region of the CSM, or in the cooling SN ejecta. Recent analysis has shown that dust formation in the CSM can commence only after day ~380, and has also ruled out newly-formed ejecta dust as the source of the NIR emission. The early (< 380 d) NIR emission can therefore only be attributed to an IR echo. The H-K color temperature of the echo is about 1250 K. The best fitting model requires the presence of about 1.6e-4 Msun of amorphous carbon dust at a distance of 2.2e16 cm from the explosion. The CSM-powered luminosity is preceded by an intense burst of hard radiation generated by the breakout of the SN shock through the stellar surface. The peak burst luminosity seen by the CSM dust is significantly reduced by Thomson scattering in the CSM, but still has the potential of evaporating the dust needed to produce the echo. We show that the survival of the echo-producing dust provides important constraints on the intensity, effective temperature, and duration of the burst., Accepted for publication in the ApJ, 20 pages, 1 table, and 17 figures

One way an Active Galactic Nucleus (AGN) influences the evolution of their host galaxy is by generating a large-scale (kpc-scale) outflow. The content, energetics, and impact of such outflows depend on the properties of both the AGN and host galaxy, and understanding the relationship between them requires measuring the properties of all three. In this paper, we do so by analyzing recent radio and optical integral field unit (IFU) spectroscopic observations of MaNGA 1-166919. Our results indicate that the bi-conical outflow in this galaxy is powered by a low-luminosity, low-Eddington ratio AGN ejecting material that drives ~100-200 km/s shocks into the surrounding interstellar medium (ISM) -- producing the hot, ionized gas and relativistic particles associated with the observed outflow. The energetics of the relativistic and ionized gas material produced at this shock are comparable, and both the mass outflow and kinetic power of the ionized gas in this outflow are higher than other AGN with similar bolometric luminosities. Lastly, while the host galaxy's total star formation rate is comparable to that of other star-forming galaxies with a similar stellar mass, there is evidence that the outflow both suppresses and enhances star formation in its immediate surroundings., 29 pages, 28 figures, 11 tables, accepted by The Astrophysical Journal for publication

Model-independent estimation of H 0 and Ω K can provide clues about the origin of the intractable Hubble constant tension and the current cosmic-curvature crisis. Strongly lensed fast radio bursts (FRBs) have been proposed as precision probes of the universe since the time delay ∼(10 days) between images could be precisely measured due to short durations ∼ of this kind of bright radio pulses. Here, on the basis of the distance sum rule, we investigate the capacity of model-independently estimating these two parameters from time-delay distances of strongly lensed FRBs and luminosity distances of the upcoming Wide Field InfraRed Survey Telescope (WFIRST) type Ia supernovae observations. Considering the expected FRB detection rate of upcoming facilities, we find that H 0 could be determined to a ∼1% precision and Ω K could be constrained to ∼0.1 simultaneously from 10 lensed FRBs. These estimations, which are independent of the energy contents of the universe and validity of Einstein’s equation on cosmological scales, will be complementary to popular probes and of great importance for clarifying the current crises of cosmology.

We have carried out an extensive X-ray spectral study of the bare Seyfert-1 galaxy MCG--02--58--22 to ascertain the nature of the X-ray reprocessing media, using observations from Suzaku (2009) and simultaneous observations from XMM-Newton and NuSTAR (2016) . The most significant results of our investigation are: 1. The primary X-ray emission from the corona is constant in these observations, both in terms of the power law slope ($\Gamma=1.80$) and luminosity ($L_{2-10 \rm keV}= 2.55\times 10^{44} $ erg/s). 2. The soft excess flux decreased by a factor of two in 2016, the Compton hump weakened/vanished in 2016, and the narrow FeK$\alpha$ emission line became marginally broad ($\sigma=0.35\pm0.08$ keV) and its flux doubled in 2016. 3. From physical model fits we find that the normalization of the narrow component of the FeK$\alpha$ line does not change in the two epochs, although the Compton hump vanishes in the same time span. Since the primary X-ray continuum does not change, we presume that any changes in the reprocessed emission must arise due to changes in the reprocessing media. Our primary conclusions are: A. The vanishing of the Compton hump in 2016 can probably be explained by a dynamic clumpy torus which is infalling/outflowing, or by a polar torus wind. B. The torus in this AGN possibly has two structures: an equatorial toroidal disk (producing the narrow FeK$\alpha$ emission) and a polar component (producing the variable Compton hump), C. The reduction of the soft-excess flux by half and increase in the FeK$\alpha$ flux by a factor of two in the same period cannot be adequately explained by ionized disk reflection model alone., Comment: Accepted for publication in ApJ

Low-frequency radio selection finds radio-bright galaxies regardless of the amount of obscuration by gas and dust. We report Chandra observations of a complete 178 MHz–selected, and so orientation-unbiased, sample of 44 0.5 < z < 1 3CRR sources. The sample is comprised of quasars and narrow-line radio galaxies (NLRGs) with similar radio luminosities, and the radio structure serves as both an age and an orientation indicator. Consistent with unification, intrinsic obscuration (measured by N H, X-ray hardness ratio, and X-ray luminosity) generally increases with inclination. However, the sample includes a population not seen in high-z 3CRR sources: NLRGs viewed at intermediate inclination angles with N H < 1022 cm−2. Multiwavelength analysis suggests that these objects have lower L/L Edd than typical NLRGs at similar orientation. Thus, both orientation and L/L Edd are important, and a “radiation-regulated unification” provides a better explanation of the sample’s observed properties. In comparison with the 3CRR sample at 1 < z < 2, our lower-redshift sample shows a higher fraction of Compton-thin NLRGs (45% versus 29%) but a similar Compton-thick fraction (20%), implying a larger covering factor of Compton-thin material at intermediate viewing angles and thus a more “puffed-up” torus atmosphere. We posit that this is due to a range of L/L Edd extending to lower values in this sample. In contrast, at high redshifts, the narrower range and high L/L Edd values allowed orientation (and so simple unification) to dominate the sample’s observed properties.

We present theoretical predictions for the free-free emission at cm wavelengths obtained from photoevaporation and MHD wind disk models adjusted to the case of the TW Hydrae young stellar object. For this system, disk photoevaporation with heating due to the high-energy photons from the star has been proposed as a possible mechanism to open the gap observed in the dust emission with ALMA. We show that the photoevaporation disk model predicts a radial profile for the free-free emission that is made of two main spatial components, one originated from the bound disk atmosphere at 0.5-1 au from the star, and another more extended component from the photoevaporative wind at larger disk radii. We also show that the stellar X-ray luminosity has a significant impact on both these components. The predicted radio emission from the MHD wind model has a smoother radial distribution which extends to closer distances to the star than the photoevaporation case. We also show that a future radio telescope such as the \textit{Next Generation Very Large Array} (ngVLA) would have enough sensitivity and angular resolution to spatially resolve the main structures predicted by these models., 11 pages, 5 figures, accepted for publication in The Astrophysical Journal

The bulk propagation speed of GeV-energy cosmic rays is limited by frequent scattering off hydromagnetic waves. Most galaxy evolution simulations that account for this confinement assume the gas is fully ionized and cosmic rays are well-coupled to Alfv\'en waves; however, multiphase density inhomogeneities, frequently under-resolved in galaxy evolution simulations, induce cosmic ray collisions and ionization-dependent transport driven by cosmic ray decoupling and elevated streaming speeds in partially neutral gas. How do cosmic rays navigate and influence such a medium, and can we constrain this transport with observations? In this paper, we simulate cosmic ray fronts impinging upon idealized, partially neutral clouds and lognormally-distributed clumps, with and without ionization-dependent transport. With these high-resolution simulations, we identify cloud interfaces as crucial regions where cosmic ray fronts can develop a stair-step pressure gradient sufficient to collisionlessly generate waves, overcome ion-neutral damping, and exert a force on the cloud. We find that the acceleration of cold clouds is hindered by only a factor of a few when ionization-dependent transport is included, with additional dependencies on magnetic field strength and cloud dimensionality. We also probe how cosmic rays sample the background gas and quantify collisional losses. Hadronic gamma-ray emission maps are qualitatively different when ionization-dependent transport is included, but the overall luminosity varies by only a small factor, as the short cosmic ray residence times in cold clouds are offset by the higher densities that cosmic rays sample., Comment: Accepted to ApJ. Revised version includes 3 tables and additional text that clarify the simulation setups and parameters. Figures 3 and 20 have been modified to show results from higher resolution simulations

Using the Monte Carlo code, SEDONA, multiband photometry and spectra are calculated for supernovae derived from stripped helium stars with presupernova masses from 2.2 to 10.0 $M_\odot$. The models are representative of evolution in close binaries and have previously been exploded using a parametrized one-dimensional model for neutrino-transport. A subset, those with presupernova masses in the range 2.2 - 5.6 $M_\odot$, have many properties in common with observed Type Ib and Ic supernovae, including a median ejected mass near 2 $M_\odot$, explosion energies near $1 \times 10^{51}$ erg, typical $^{56}$Ni masses 0.07 - 0.09 $M_\odot$, peak times of about 20 days, and a narrow range for the $V$-$R$ color index 10 days post $V$-maximum near 0.3 mag. The median peak bolometric luminosity, near 10$^{42.3}$ erg s$^{-1}$, is fainter, however, than for several observational tabulations and the brightest explosion has a bolometric luminosity of only 10$^{42.50}$ erg s$^{-1}$. The brightest absolute $B$, $V$, and $R$ magnitudes at peak are $-17.2$, $-17.8$, and $-18.0$. These limits are fainter than some allegedly typical Type Ib and Ic supernovae and could reflect problems in our models or the observational analysis. Helium stars with lower and higher masses also produce interesting transients that may have been observed including fast, faint, blue transients and long, red, faint Type Ic supernovae. New models are specifically presented for SN 2007Y, SN 2007gr, SN 2009jf, LSQ13abf, SN 2008D, and SN 2010X., Comment: submitted to ApJ

We present the analytical framework for converting projected light distributions with a S\'ersic profile into three-dimensional light distributions for stellar systems of arbitrary triaxial shape. The main practical result is the definition of a simple yet robust measure of intrinsic galaxy size: the median radius $r_\mathrm{med}$, defined as the radius of a sphere that contains 50% of the total luminosity or mass, that is, the median distance of a star to the galaxy center. We examine how $r_\mathrm{med}$ depends on projected size measurements as a function of S\'ersic index and intrinsic axis ratios, and demonstrate its relative independence of these parameters. As an application we show that the projected semi-major axis length of the ellipse enclosing 50% of the light is an unbiased proxy for $r_\mathrm{med}$, with small galaxy-to-galaxy scatter of $\sim$10% (1$\sigma$), under the condition that the variation in triaxiality within the population is small. For galaxy populations with unknown or a large range in triaxiality an unbiased proxy for $r_\mathrm{med}$ is $1.3\times R_{e}$, where $R_{e}$ is the circularized half-light radius, with galaxy-to-galaxy scatter of 20-30% (1$\sigma$). We also describe how inclinations can be estimated for individual galaxies based on the measured projected shape and prior knowledge of the intrinsic shape distribution of the corresponding galaxy population. We make the numerical implementation of our calculations available., Comment: 7 pages, 4 figures; accepted for publication in ApJ

We matched the 1.4 GHz local luminosity functions of star-forming galaxies (SFGs) and active galactic nuclei to the 1.4 GHz differential source counts from $0.25 \ \mu\mathrm{Jy}$ to 25 Jy using combinations of luminosity and density evolution. We present the most robust and complete local far-infrared (FIR)/radio luminosity correlation to date in a volume-limited sample of $\approx 4.3 \times 10^3$ nearby SFGs, finding that it is very tight but distinctly sub-linear: $L_\mathrm{FIR} \propto L_\mathrm{1.4\,GHz}^{0.85}$. If the local FIR/radio correlation does not evolve, the evolving 1.4 GHz luminosity function of SFGs yields the evolving star-formation rate density (SFRD) $\psi (M_\odot \ \mathrm{year}^{-1} \ \mathrm{Mpc}^{-3}$) as a function of time since the big bang. The SFRD measured at 1.4 GHz grows rapidly at early times, peaks at "cosmic noon" when $t \approx 3 \ \mathrm{Gyr}$ and $z \approx 2$, and subsequently decays with an $e$-folding time scale $\tau = 3.2 \ \mathrm{Gyr}$. This evolution is similar to, but somewhat stronger than, SFRD evolution estimated from UV and FIR data., Comment: 19 pages, 13 figures; accepted for publication in ApJ

We present observations of the Trumpler 14/Carina I region carried out using the Stratospheric Terahertz Observatory 2 (STO2). The Trumpler 14/Carina I region is in the west part of the Carina Nebula Complex, which is one of the most extreme star-forming regions in the Milky Way. We observed Trumpler 14/Carina I in the 158 $\mu$m transition of [C\,{\sc ii}] with a spatial resolution of 48$''$ and a velocity resolution of 0.17 km s$^{-1}$. The observations cover a 0.25$^\circ$ by 0.28$^\circ$ area with central position {\it l} = 297.34$^\circ$, {\it b} = -0.60$^\circ$. The kinematics show that bright [C\,{\sc ii}] structures are spatially and spectrally correlated with the surfaces of CO clouds, tracing the photodissociation region and ionization front of each molecular cloud. Along 7 lines of sight that traverse Tr 14 into the dark ridge to the southwest, we find that the [C\,{\sc ii}] luminosity from the HII region is 3.7 times that from the PDR. In same los we find in the PDRs an average ratio of 1:4.1:5.6 for the mass in atomic gas:dark-CO gas: molecular gas traced by CO. Comparing multiple gas tracers including HI 21cm, [C\,{\sc ii}], CO, and radio recombination lines, we find that the HII regions of the Carina Nebula Complex are well-described as HII regions with one-side freely expanding towards us, consistent with the champagne model of ionized gas evolution. The dispersal of the GMC in this region is dominated by EUV photoevaporation; the dispersal timescale is 20-30 Myr., Comment: ApJ accepted

We measure the average deflection of cosmic microwave background photons by quasars at $\langle z \rangle =1.7$. Our sample is selected from the Sloan Digital Sky Survey to cover the redshift range $0.9\leq z\leq2.2$ with absolute i-band magnitudes of $M_i\leq-24$ (K-corrected to z=2). A stack of nearly 200,000 targets reveals an 8$\sigma$ detection of Planck's estimate of the lensing convergence towards the quasars. We fit the signal with a model comprising a Navarro-Frenk-White density profile and a 2-halo term accounting for correlated large scale structure, which dominates the observed signal. The best-fitting model is described by an average halo mass $\log_{10}(M_{\rm h}/h^{-1}M_\odot)=12.6\pm0.2$ and linear bias $b=2.7\pm0.3$ at $z=1.7$, in excellent agreement with clustering studies. We also report of a hint, at a 90% confidence level, of a correlation between the convergence amplitude and luminosity, indicating that quasars brighter than $M_i\lesssim -26$ reside in halos of typical mass ${M_{\rm h}\approx 10^{13}\,h^{-1}M_\odot}$, scaling roughly as ${M_{\rm h}\propto L_{\rm opt}^{3/4}}$ at ${M_i\lesssim-24}$, in good agreement with physically-motivated quasar demography models. Although we acknowledge this luminosity dependence is a marginal result, the observed $M_{\rm h}$-$L_{\rm opt}$ relationship could be interpreted as a reflection of the cutoff in the distribution of black hole accretion rates towards high Eddington ratios: the weak trend of $M_{\rm h}$ with $L_{\rm opt}$ observed at low luminosity becomes stronger for the most powerful quasars, which tend to be accreting close to the Eddington limit., Comment: Accepted for publication in ApJ

We present Zwicky Transient Facility (ZTF) observations of the tidal disruption flare AT2018zr/PS18kh reported by Holoien et al. and detected during ZTF commissioning. The ZTF light curve of the tidal disruption event (TDE) samples the rise-to-peak exceptionally well, with 50 days of g- and r-band detections before the time of maximum light. We also present our multi-wavelength follow-up observations, including the detection of a thermal (kT~100 eV) X-ray source that is two orders of magnitude fainter than the contemporaneous optical/UV blackbody luminosity, and a stringent upper limit to the radio emission. We use observations of 128 known active galactic nuclei (AGN) to assess the quality of the ZTF astrometry, finding a median host-flare distance of 0.2" for genuine nuclear flares. Using ZTF observations of variability from known AGN and supernovae we show how these sources can be separated from TDEs. A combination of light-curve shape, color, and location in the host galaxy can be used to select a clean TDE sample from multi-band optical surveys such as ZTF or LSST., accepted to ApJ, updated with recent Swift data

In the hierarchical theory of galaxy formation, a galaxy overdensity is a hallmark of a massive cosmic structure. However, it is less well understood how different types of galaxies trace the underlying large-scale structure. Motivated by the discovery of a z=3.13 protocluster, we examine how the same structure is populated by Ly$\alpha$-emitting galaxies (LAEs). To this end, we have undertaken a deep narrow-band imaging survey sampling Ly$\alpha$ emission at this redshift. Of the 93 LAE candidates within a 36'x36'~(70x70~Mpc^2) field, 21 galaxies form a significant surface overdensity (delta_g=3.3+/-0.9), which is spatially segregated from the Lyman break galaxy (LBG) overdensity. One possible interpretation is that they trace two separate structures of comparable masses (~ 10^{15}M_sun) where the latter is hosted by a halo assembled at an earlier time. We speculate that the dearth of LAEs in the LBG overdensity region may signal the role of halo assembly bias in galaxy formation, which would suggest that different search techniques may be biased accordingly to the formation age or dynamical state of the host halo. The median Ly$\alpha$- and UV luminosity is 30--70\% higher for the protocluster LAEs relative to the field. This difference cannot be explained by the galaxy overdensity alone, and may require a top-heavy mass function, higher star formation efficiency for protocluster halos, or suppression of galaxy formation in low-mass halos. A luminous Ly$\alpha$ blob and an ultramassive galaxy found in this region paint a picture consistent with the expected early growth of galaxies in clusters., Comment: Accepted to ApJ

We use the results from the ALMA large program ASPECS, the spectroscopic survey in the Hubble Ultra Deep Field (HUDF), to constrain CO luminosity functions of galaxies and the resulting redshift evolution of $\rho$(H$_2$). The broad frequency range covered enables us to identify CO emission lines of different rotational transitions in the HUDF at $z>1$. We find strong evidence that the CO luminosity function evolves with redshift, with the knee of the CO luminosity function decreasing in luminosity by an order of magnitude from $\sim$2 to the local universe. Based on Schechter fits, we estimate that our observations recover the majority (up to $\sim$90%, depending on the assumptions on the faint end) of the total cosmic CO luminosity at $z$=1.0-3.1. After correcting for CO excitation, and adopting a Galactic CO-to-H$_2$ conversion factor, we constrain the evolution of the cosmic molecular gas density $\rho$(H$_2$): this cosmic gas density peaks at $z\sim1.5$ and drops by factor of $6.5_{-1.4}^{+1.8}$ to the value measured locally. The observed evolution in $\rho$(H$_2$) therefore closely matches the evolution of the cosmic star formation rate density $\rho_{\rm SFR}$. We verify the robustness of our result with respect to assumptions on source inclusion and/or CO excitation. As the cosmic star formation history can be expressed as the product of the star formation efficiency and the cosmic density of molecular gas, the similar evolution of $\rho$(H$_2$) and $\rho_{\rm SFR}$ leaves only little room for a significant evolution of the average star formation efficiency in galaxies since $z\sim 3$ (85% of cosmic history)., Comment: 22 pages, 11 figures. Paper re-submitted to ApJ after addressing the first round of comments by the referee

A pair of non-thermal radio bubbles recently discovered in the inner few hundred parsecs of the Galactic center bears a close spatial association with elongated, thermal X-ray features called the X-ray chimneys. While their morphology, position, and orientation vividly point to an outflow from the Galactic center, the physical processes responsible for the outflow remain to be understood. We use three-dimensional magnetohydrodynamic simulations to test the hypothesis that the radio bubbles/X-ray chimneys are the manifestation of an energetic outflow driven by multiple core-collapsed supernovae in the nuclear stellar disk, where numerous massive stars are known to be present. Our simulations are run with different combinations of two main parameters, the supernova birth rate and the strength of a global magnetic field being vertically oriented with respect to the disk. The simulation results show that a hot gas outflow can naturally form and acquire a vertically elongated shape due to collimation by the magnetic pressure. In particular, the simulation with an initial magnetic field strength of 80 $\mu$G and a supernova rate of 1 $kyr^{-1}$ can well reproduce the observed morphology, internal energy and X-ray luminosity of the bubbles after an evolutionary time of 330 kyr. On the other hand, a magnetic field strength of 200 $\mu$G gives rise to an overly elongated outflow that is inconsistent with the observed bubbles. The simulations also reveal that, inside the bubbles, mutual collisions between the shock waves of individual supernovae produce dense filaments of locally amplified magnetic field. Such filaments may account for a fraction of the synchrotron-emitting radio filaments known to exist in the Galactic center., Comment: 20 pages, 10 figures, accepted by ApJ

The Atacama Large Millimeter/submillimeter Array (ALMA) Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) Band 6 scan (212–272 GHz) covers potential [C ii] emission in galaxies at 6 ≤ z ≤ 8 throughout a 2.9 arcmin2 area. By selecting on known Lyα emitters (LAEs) and photometric dropout galaxies in the field, we perform targeted searches down to a 5σ [C ii] luminosity depth L [C II] ∼ 2.0 × 108 L ⊙, corresponding roughly to star formation rates (SFRs) of 10–20 M ⊙ yr−1 when applying a locally calibrated conversion for star-forming galaxies, yielding zero detections. While the majority of galaxies in this sample are characterized by lower SFRs, the resulting upper limits on [C ii] luminosity in these sources are consistent with the current literature sample of targeted ALMA observations of z = 6–7 LAEs and Lyman-break galaxies (LBGs), as well as the locally calibrated relations between L [C ii] and SFR—with the exception of a single [C ii]-deficient, UV-luminous LBG. We also perform a blind search for [C ii]-bright galaxies that may have been missed by optical selections, resulting in an upper limit on the cumulative number density of [C ii] sources with L [C II] > 2.0 × 108 L ⊙ (5σ) to be less than 1.8 × 10−4 Mpc−3 (90% confidence level). At this luminosity depth and volume coverage, we present an observed evolution of the [C ii] luminosity function from z = 6–8 to z ∼ 0 by comparing the ASPECS measurement to literature results at lower redshift.

We probe the dusty environment of the archetypical type-1 AGN in NGC 4151 with comprehensive IR reverberation mapping over several decades, in J ($\sim1.22\mu$m), H ($\sim1.63\mu$m), K ($\sim2.19\mu$m), L ($\sim3.45\mu$m) and N-band ($\sim10.6\mu m$), plus multiple measurements at 20--40 $\mu$m. At 1--4$\mu$m, the hot dust reverberation signals come from two distinct dust populations at separate radii ($\sim$0.033 pc and $\sim$0.076 pc), with temperatures of $\sim$1500--2500 K and $\sim$900--1000 K, consistent with the expected properties of sublimating graphite and silicate dust grains. The domination of the torus infrared output by carbon and silicate grains near their sublimation temperatures and radii may account for the general similarity of AGN near-IR spectral energy distributions. The torus inner edge defined by the hottest dust remains at roughly the same radius independent of the AGN optical luminosity over $\sim$25 years. The emission by hot dust warmed directly by the optical/UV AGN output has increased gradually by $\sim$4\%/year, indicating a possibly growing torus. A third dust component at $\sim$700 K does not seem to participate directly in the IR reverberation behavior and its emission may originate deep in the circumnuclear torus. We find a reverberation signal at $\sim10\mu$m with an inferred radius for the warm dust of $\sim$2.2--3.1 pc. The lack of variability at 20--40$\mu$m indicates the far-IR emission comes from even more extended regions. The torus properties revealed by dust reverberation analysis are inconsistent with predictions from pure clumpy torus models. Instead, the longer wavelength emission possibly originates in a flared torus or the polar wind., Comment: 34 pages, 23 figures, 2 tables; ApJ in press

Recent studies of gamma-ray, cosmic-ray, and radio data put stringent constraints on the fraction of primordial black holes (PBHs) in our universe. In this article, we propose a new indirect method in using the X-ray luminosity data of cool-core clusters to constrain the evaporating PBH fraction for the monochromatic, log-normal and power-law mass distributions. The present results show that the amount of evaporating PBHs only constitutes a minor component of dark matter for a large parameter space. The constraints are consistent with and close to that obtained from other cosmic-ray and multiwavelength observations.

Building on recent work by Chandar et al. (2020), we construct X-ray luminosity functions (XLFs) for different classes of X-ray binary (XRB) donors in the nearby star-forming galaxy M83 through a novel methodology: rather than classifying low- vs. high-mass XRBs based on the scaling of the number of X-ray sources with stellar mass and star formation rate, respectively, we utilize multi-band Hubble Space Telescope imaging data to classify each Chandra-detected compact X-ray source as a low-mass (i.e. donor mass ~ 8 solar masses) or intermediate-mass XRB based on either the location of its candidate counterpart on optical color-magnitude diagrams or the age of its host star cluster. In addition to the the standard (single and/or truncated) power-law functional shape, we approximate the resulting XLFs with a Schechter function. We identify a marginally significant (at the 1-to-2 sigma level) exponential downturn for the high-mass XRB XLF, at logLx ~ 38.48^{+0.52}_{-0.33} (in log CGS units). In contrast, the low- and intermediate-mass XRB XLFs, as well as the total XLF of M83, are formally consistent with sampling statistics from a single power-law. Our method suggests a non-negligible contribution from low- and possibly intermediate-mass XRBs to the total XRB XLF of M83, i.e. between 20 and 50%, in broad agreement with X-ray based XLFs. More generally, we caution against considerable contamination from X-ray emitting supernova remnants to the published, X-ray based XLFs of M83, and possibly all actively star-forming galaxies., Comment: 18 pages (text, 44 with figures), 8 figures + 1 long figure, 1 table + 1 long table. Accepted by ApJ

Outflows of ionized gas driven by active galactic nuclei (AGN) may significantly impact the evolution of their host galaxies. However, determining the energetics of these outflows is difficult with spatially unresolved observations that are subject to strong global selection effects. We present part of an ongoing study using Hubble Space Telescope (HST) and Apache Point Observatory (APO) spectroscopy and imaging to derive spatially-resolved mass outflow rates and energetics for narrow line region (NLR) outflows in nearby AGN that are based on multi-component photoionization models to account for spatial variations in the gas ionization, density, abundances, and dust content. This expanded analysis adds Mrk 3, Mrk 78, and NGC 1068, doubling the sample in Revalski (2019). We find that the outflows contain total ionized gas masses of $M \approx 10^{5.5} - 10^{7.5}$ $M_{\odot}$ and reach peak velocities of $v \approx 800 - 2000$ km s$^{-1}$. The outflows reach maximum mass outflow rates of $\dot M_{out} \approx 3 - 12$ $M_{\odot}$ yr$^{-1}$ and encompass total kinetic energies of $E \approx 10^{54} - 10^{56}$ erg. The outflows extend to radial distances of $r \approx 0.1 - 3$ kpc from the nucleus, with the gas masses, outflow energetics, and radial extents positively correlated with AGN luminosity. The outflow rates are consistent with in-situ ionization and acceleration where gas is radiatively driven at multiple radii. These radial variations indicate that spatially-resolved observations are essential for localizing AGN feedback and determining the most accurate outflow parameters., Accepted for Publication in ApJ on January 12, 2021. The paper has 40 pages and 15 figures, with results tabulated in the Appendix. Version two includes minor corrections to match the journal publication

The detailed observations of GW170817 proved for the first time directly that neutron star mergers are a major production site of heavy elements. The observations could be fit by a number of simulations that qualitatively agree, but can quantitatively differ (e.g. in total r-process mass) by an order of magnitude. We categorize kilonova ejecta into several typical morphologies motivated by numerical simulations, and apply a radiative transfer Monte Carlo code to study how the geometric distribution of the ejecta shapes the emitted radiation. We find major impacts on both spectra and light curves. The peak bolometric luminosity can vary by two orders of magnitude and the timing of its peak by a factor of five. These findings provide the crucial implication that the ejecta masses inferred from observations around the peak brightness are uncertain by at least an order of magnitude. Mixed two-component models with lanthanide-rich ejecta are particularly sensitive to geometric distribution. A subset of mixed models shows very strong viewing angle dependence due to lanthanide "curtaining," which persists even if the relative mass of lanthanide-rich component is small. The angular dependence is weak in the rest of our models, but different geometric combinations of the two components lead to a highly diverse set of light curves. We identify geometry-dependent {P Cygni} features in late spectra that directly map out strong lines in the simulated opacity of neodymium, which can help to constrain the ejecta geometry and to directly probe the r-process abundances., Comment: 23 pages, 18 figures; accepted to ApJ

We present the results of a single-dish survey toward 95 very low luminosity objects (VeLLOs) in optically thick (HCN 1−0) and thin (N2H+ 1−0) lines performed for the purpose of understanding the physical processes of inward motions in the envelopes of the VeLLOs and characterizing their true nature. The normalized velocity differences ( δ V HCN ) between the peak velocities of the two lines were derived for 41 VeLLOs detected in both lines. The δ V distribution of these VeLLOs is found to be significantly skewed to the blue, indicating the dominance of infalling motions in their envelopes. The infall speeds were derived for 15 infall candidates by using the HILL5 radiative transfer model. The speeds were in the range of 0.03−0.3 km s−1, with a median value of 0.16 km s−1, consistent with the gravitational freefall speeds from pressure-free envelopes. The mass infall rates calculated from the infall speeds are mostly of the order of 10−6 M ⊙ yr−1, with a median value of (3.4 ± 1.5) × 10−6 M ⊙ yr−1. These are found to be also consistent with the values predicted with the inside-out collapse model and show a fairly good correlation with the internal luminosities of the VeLLOs. This again indicates that the infall motions observed toward the VeLLOs are likely to be due to the gravitational infall motions in their envelopes. Our study suggests that most of the VeLLOs are potentially faint protostars, while two of the VeLLOs could possibly be proto−brown dwarf candidates.