Search

How Supernova Remnant (SNR) shocks impact nearby molecular clouds is still poorly observationally constrained. It is unclear if SNRs can positively or negatively affect clouds star formation potential. We have studied the dense gas morphology and kinematics toward the Infrared Dark Cloud (IRDC) G034.77-00.55, shock-interacting with the SNR W44, to identify evidence of early stage star formation induced by the shock. We have used high-angular resolution N2H+(1-0) images across G034.77-00.55, obtained with ALMA. N2H+ is a well known tracer of dense and cold material, optimal to identify gas with the highest potential to harbour star formation. The N2H+ emission is distributed into two elongated structures, one toward the dense ridge at the edge of the source and one toward the inner cloud. Both elongations are spatially associated with well-defined mass-surface density features. The velocities of the gas in the two structures i.e., 38-41 km s-1 and 41-43 km s-1 are consistent with the lowest velocities of the J- and C-type parts of the SNR-driven shock, respectively. A third velocity component is present at 43-45.5 km s-1. The dense gas shows a fragmented morphology with core-like fragments of scales consistent with the Jeans lengths, masses $\sim$1-20 M$_{\odot}$, densities (n(H$_2$)$\geq$10$^5$ cm$^{-3}$) sufficient to host star formation in free-fall time scales (few 10$^4$ yr) and with virial parameters that hint toward possible collapse. The W44 driven shock may have swept up the encountered material which is now seen as a dense ridge, almost detached from the main cloud, and an elongation within the inner cloud, well constrained in both N2H+ emission and mass surface density. This shock compressed material may have then fragmented into cores that are either in a starless or pre-stellar stage. Additional observations are needed to confirm this scenario and the nature of the cores., Comment: 12 pages, 8 figures, accepted for publication in A&A

The increasing number of observations towards different environments in the Milky Way, as well as theoretical and experimental works, are improving our knowledge of the astrochemical processes in the interstellar medium (ISM). In this chapter we report some of the main projects to study the chemical complexity and isotopic ratios across the Galaxy. High-sensitivity spectral surveys covering broad bandwidths towards Galactic Center molecular clouds (e.g. G+0.693-0.027) and star-forming regions (e.g. the hot core G31.41+0.31) are revealing very rich astrochemical reservoirs, which include molecules of prebiotic interest. At the same time, isotopic ratios (e.g. $^{12}$C/$^{13}$C and $^{14}$N/$^{15}$N) can give important information on the Galactic chemical evolution, as well as on chemical local processes due to the physical conditions of the molecular clouds. We also highlight the role of cosmic rays as a key agent affecting the interstellar chemistry described above., Comment: 9 pages, 1 figures - Chapter II of the National Congress of (proto-)Planetary Astrochemistry (CNAP II) Proceedings - Accepted for publication on MemSAIt, in print

We use data from the ALMA Evolutionary Study of High Mass Protocluster Formation in the Galaxy (ALMAGAL) survey to study 100 ALMAGAL regions at $\sim$ 1 arsecond resolution located between $\sim$ 2 and 6 kpc distance. Using ALMAGAL $\sim$ 1.3mm line and continuum data we estimate flow rates onto individual cores. We focus specifically on flow rates along filamentary structures associated with these cores. Our primary analysis is centered around position velocity cuts in H$_2$CO (3$_{0,3}$ - 2$_{0,2}$) which allow us to measure the velocity fields, surrounding these cores. Combining this work with column density estimates we derive the flow rates along the extended filamentary structures associated with cores in these regions. We select a sample of 100 ALMAGAL regions covering four evolutionary stages from quiescent to protostellar, Young Stellar Objects (YSOs), and HII regions (25 each). Using dendrogram and line analysis, we identify a final sample of 182 cores in 87 regions. In this paper, we present 728 flow rates for our sample (4 per core), analysed in the context of evolutionary stage, distance from the core, and core mass. On average, for the whole sample, we derive flow rates on the order of $\sim$10$^{-4}$ M$_{sun}$yr$^{-1}$ with estimated uncertainties of $\pm$50%. We see increasing differences in the values among evolutionary stages, most notably between the less evolved (quiescent/protostellar) and more evolved (YSO/HII region) sources. We also see an increasing trend as we move further away from the centre of these cores. We also find a clear relationship between the flow rates and core masses $\sim$M$^{2/3}$ which is in line with the result expected from the tidal-lobe accretion mechanism. Overall, we see increasing trends in the relationships between the flow rate and the three investigated parameters; evolutionary stage, distance from the core, and core mass., Comment: 11 pages, 11 figures, accepted for publication in A&A

Ice-coated dust grains provide the main reservoir of volatiles that play an important role in planet formation processes and may become incorporated into planetary atmospheres. However, due to observational challenges, the ice abundance distribution in protoplanetary disks is not well constrained. We present JWST/MIRI observations of the edge-on disk HH 48 NE carried out as part of the IRS program Ice Age. We detect CO$_2$, NH$_3$, H$_2$O and tentatively CH$_4$ and NH$_4^+$. Radiative transfer models suggest that ice absorption features are produced predominantly in the 50-100 au region of the disk. The CO$_2$ feature at 15 micron probes a region closer to the midplane (z/r = 0.1-0.15) than the corresponding feature at 4.3 micron (z/r = 0.2-0.6), but all observations trace regions significantly above the midplane reservoirs where we expect the bulk of the ice mass to be located. Ices must reach a high scale height (z/r ~ 0.6; corresponding to modeled dust extinction Av ~ 0.1), in order to be consistent with the observed vertical distribution of the peak ice optical depths. The weakness of the CO$_2$ feature at 15 micron relative to the 4.3 micron feature and the red emission wing of the 4.3 micron CO$_2$ feature are both consistent with ices being located at high elevation in the disk. The retrieved NH$_3$ abundance and the upper limit on the CH$_3$OH abundance relative to H$_2$O are significantly lower than those in the interstellar medium (ISM), but consistent with cometary observations. Full wavelength coverage is required to properly study the abundance distribution of ices in disks. To explain the presence of ices at high disk altitudes, we propose two possible scenarios: a disk wind that entrains sufficient amounts of dust, thus blocking part of the stellar UV radiation, or vertical mixing that cycles enough ices into the upper disk layers to balance ice photodesorption., Comment: 16 pages, 11 figures, accepted in A&A

With two directly detected protoplanets, the PDS 70 system is a unique source in which to study the complex interplay between forming planets and their natal environment. The large dust cavity carved by the two giant planets can affect the disk chemistry, and therefore the molecular emission morphology. On the other hand, chemical properties of the gas component of the disk are expected to leave an imprint on the planetary atmospheres. In this work, we reconstruct the emission morphology of a rich inventory of molecular tracers in the PDS 70 disk, and we look for possible chemical signatures of the two actively accreting protoplanets, PDS b and c. We leverage Atacama Large Millimeter/submillimeter Array (ALMA) band 6 high-angular-resolution and deep-sensitivity line emission observations, together with image and $uv$-plane techniques, to boost the detection of faint lines. We robustly detect ring-shaped emission from $^{12}$CO, $^{13}$CO, C$^{18}$O, H$^{13}$CN, HC$^{15}$N, DCN, H$_2$CO, CS, C$_2$H, and H$^{13}$CO$^{+}$ lines in unprecedented detail. Most of the molecular tracers show a peak of the emission inside the millimeter dust peak. We interpret this as the direct impact of the effective irradiation of the cavity wall, as a result of the planet formation process. Moreover, we have found evidence of an O-poor gas reservoir in the outer disk, which is supported by the observations of bright C-rich molecules, the non-detection of SO, and a lower limit on the $\mathrm{CS/SO}$ ratio of $\sim1$. Eventually, we provide the first detection of the c-C$_3$H$_2$ transitions at 218.73 GHz, and the marginal detection of an azimuthal asymmetry in the higher-energy H$_2$CO (3$_{2,1}$-2$_{2,0}$) line, which could be due to accretion heating near PDS 70b., Comment: 25 pages, 10 figures, 3 tables, 8 figures and one table in appendix. Accepted for publication in A&A

Context. Dust polarization observations of the massive protocluster G31.41+0.31 carried out at ~1'' (~3750 au) resolution with the SMA at 870 microm have revealed one of the clearest examples to date of an hourglass-shaped magnetic field morphology in the high-massregime. Additionally, ~0.24'' (~900 au) resolution observations with ALMA at 1.3 mm have confirmed these results. The next step is to investigate whether the magnetic field maintains its hourglass-shaped morphology down to circumstellar scales. Aims. To study the magnetic field morphology toward the four (proto)stars A, B, C, and D contained in G31.41+0.31 and examine whether the self-similarity observed at core scales (1'' and 0.24'' resolution) still holds at circumstellar scales, we carried out ALMA observations of the polarized dust continuum emission at 1.3 mm and 3.1 mm at an angular resolution of ~0.068'' (~250 au), sufficient to resolve the envelope emission of the embedded protostars. Methods. We used ALMA to perform full polarization observations at 233 GHz (Band 6) and 97.5 GHz (Band 3) with a synthesized beam of 0.072'' x 0.064''. We carried out polarization observations at two different wavelengths to confirm that the polarization traces magnetically aligned dust grains and is not due to dust self-scattering. Results. The polarized emission and the direction of the magnetic field obtained at the two wavelengths are basically the same, except for an area between the embedded sources C and B. In such an area, the emission at 1.3 mm could be optically thick and affected by dichroic extinction. ..., Comment: 15 pages, 15 figures, Accepted for publication in A&A

Over the past few years, chemical studies have revealed multiple structures in the vicinity of young stellar objects (YSOs). It has become evident that specific physical conditions are associated with the emission of particular molecular lines, allowing us to use molecular probes of the YSO physics. Consequently, chemical surveys are now necessary to fully constrain the origin of the observed structures. Several surveys have been conducted to explore the chemistry of YSOs, focusing on Class 0 and Class II objects. However, our knowledge of intermediate objects, that are Class I objects, remains limited. To bridge the gap and establish the relationship between observed structures and molecular line emission at the Class I evolutionary stage, we investigate the spatial distribution of key molecular gas species in the low-mass Class I protostar L1489 IRS (IRAS 04016+2610), a source part of the ChemYSO survey. We performed a 3mm line survey at high spatial and high spectral resolution using the NOEMA interferometer and the IRAM-30m telescope. We present here the ten brightest lines of our survey, in which we identified a new ~ 3 000 au long streamer in HC3N, C2H, and c-C3H2 emission, likely associated with more localized accretion shocks probed in SO. In addition, two ~ 10 000 au bubbles are seen with the dense molecular tracers HCO+, CS, and HCN around the YSO. Additionally, potential indicators of a second outflow appear in CS and HCN emission, but its nature remains to be confirmed. The late infall identified at large scales may originate from the nearby prestellar core L1489 and is likely responsible for the formation of an external warped disk in this system. The detection of a potential second outflow could be the direct evidence of a binary system. Finally, we hypothesize that the bubbles may result from the magnetic pressure as observed in numerical simulations., Comment: 16 pages, 12 figures, accepted for publication in Astronomy & Astrophysics

Magnetic fields may play a crucial role in setting the initial conditions of massive star and star cluster formation. To investigate this, we report SOFIA-HAWC+ $214\:\mu$m observations of polarized thermal dust emission and high-resolution GBT-Argus C$^{18}$O(1-0) observations toward the massive Infrared Dark Cloud (IRDC) G28.37+0.07. Considering the local dispersion of $B$-field orientations, we produce a map of $B$-field strength of the IRDC, which exhibits values between $\sim0.03 - 1\:$mG based on a refined Davis-Chandrasekhar-Fermi (r-DCF) method proposed by Skalidis \& Tassis. Comparing to a map of inferred density, the IRDC exhibits a $B-n$ relation with a power law index of $0.51\pm0.02$, which is consistent with a scenario of magnetically-regulated anisotropic collapse. Consideration of the mass-to-flux ratio map indicates that magnetic fields are dynamically important in most regions of the IRDC. A virial analysis of a sample of massive, dense cores in the IRDC, including evaluation of magnetic and kinetic internal and surface terms, indicates consistency with virial equilibrium, sub-Alfv\'enic conditions and a dominant role for $B-$fields in regulating collapse. A clear alignment of magnetic field morphology with direction of steepest column density gradient is also detected. However, there is no preferred orientation of protostellar outflow directions with the $B-$field. Overall, these results indicate that magnetic fields play a crucial role in regulating massive star and star cluster formation and so need to be accounted for in theoretical models of these processes., Comment: Submitted to ApJ, comments welcome

We present the first compelling evidence of shock-heated molecular clouds associated with the supernova remnant (SNR) N49 in the Large Magellanic Cloud (LMC). Using $^{12}$CO($J$ = 2-1, 3-2) and $^{13}$CO($J$ = 2-1) line emission data taken with the Atacama Large Millimeter/Submillimeter Array, we derived the H$_2$ number density and kinetic temperature of eight $^{13}$CO-detected clouds using the large velocity gradient approximation at a resolution of 3.5$''$ (~0.8 pc at the LMC distance). The physical properties of the clouds are divided into two categories: three of them near the shock front show the highest temperatures of ~50 K with densities of ~500-700 cm$^{-3}$, while other clouds slightly distant from the SNR have moderate temperatures of ~20 K with densities of ~800-1300 cm$^{-3}$. The former clouds were heated by supernova shocks, but the latter were dominantly affected by the cosmic-ray heating. These findings are consistent with the efficient production of X-ray recombining plasma in N49 due to thermal conduction between the cold clouds and hot plasma. We also find that the gas pressure is roughly constant except for the three shock-engulfed clouds inside or on the SNR shell, suggesting that almost no clouds have evaporated within the short SNR age of ~4800 yr. This result is compatible with the shock-interaction model with dense and clumpy clouds inside a low-density wind bubble., Comment: 12 pages, 5 figures, 1 table, accepted for publication in The Astrophysical Journal (ApJ)

Ices are the main carriers of volatiles in protoplanetary disks and are crucial to our understanding of the chemistry that ultimately sets the organic composition of planets. The ERS program Ice Age on the JWST follows the ice evolution through all stages of star and planet formation. JWST/NIRSpec observations of the edge-on Class II protoplanetary disk HH~48~NE reveal spatially resolved absorption features of the major ice components H$_2$O, CO$_2$, CO, and multiple weaker signatures from less abundant ices NH$_3$, OCN$^-$, and OCS. Isotopologue $^{13}$CO$_2$ ice has been detected for the first time in a protoplanetary disk. Since multiple complex light paths contribute to the observed flux, the ice absorption features are filled in by ice-free scattered light. The $^{12}$CO$_2$/$^{13}$CO$_2$ ratio of 14 implies that the $^{12}$CO$_2$ feature is saturated, without the flux approaching 0, indicative of a very high CO$_2$ column density on the line of sight, and a corresponding abundance with respect to hydrogen that is higher than ISM values by a factor of at least a few. Observations of rare isotopologues are crucial, as we show that the $^{13}$CO$_2$ observation allows us to determine the column density of CO$_2$ to be at an order of magnitude higher than the lower limit directly inferred from the observed optical depth. Radial variations in ice abundance, e.g., snowlines, are significantly modified since all observed photons have passed through the full radial extent of the disk. CO ice is observed at perplexing heights in the disk, extending to the top of the CO-emitting gas layer. We argue that the most likely interpretation is that we observe some CO ice at high temperatures, trapped in less volatile ices like H$_2$O and CO$_2$. Future radiative transfer models will be required to constrain the implications on our current understanding of disk physics and chemistry., Comment: 16 pages, 8 figures, accepted for publication in A&A

Fast Radio Bursts (FRBs) are a powerful and mysterious new class of transient that are luminous enough to be detected at cosmological distances. By associating FRBs to host galaxies, we can measure intrinsic and environmental properties that test FRB origin models, in addition to using them as precise probes of distant cosmic gas. The Deep Synoptic Array (DSA-110) is a radio interferometer built to maximize the rate at which it can simultaneously detect and localize FRBs. Here, we present the first sample of FRBs and host galaxies discovered by the DSA-110. This sample of 11 FRBs is the largest, most uniform sample of localized FRBs to date, as it is selected based on association to host galaxies identified in optical imaging by Pan-STARRS1. These FRBs have not been observed to repeat and their radio properties (dispersion, temporal scattering, energy) are similar to that of the known non-repeating FRB population. Most host galaxies have ongoing star formation, as has been identified before for FRB hosts. Two hosts of the new sample are massive, quiescent galaxies. The distribution of star-formation history across this host-galaxy sample shows that the delay-time distribution is wide, with a powerlaw model that spans from $\sim100$\,Myr to $\gtrsim2$\,Gyr. This requires the existence of one or more progenitor formation channels associated with old stellar populations, such as the binary evolution of compact objects., Comment: 21 pages. Submitted to AAS Journals. Includes changes based on referee comments and improved host galaxy analysis

Embedded planets are potentially the cause of substructures like gaps and cavities observed in several protoplanetary disks. Thus, the substructures observed in the continuum and in line emission encode information about the presence of planets in the system and how they interact with the natal disk. The pre-transitional disk around the star PDS 70 is the first case of two young planets imaged within a dust depleted gap that was likely carved by themselves. We aim to determine the spatial distribution of the gas and dust components in the PDS 70 disk. The axisymmetric substructures observed in the resulting profiles are interpreted in the context of planet-disk interactions. We develop a thermo-chemical forward model for an axisymmetric disk to explain a subset of the Atacama Large Millimeter/Submillimeter Array (ALMA) band 6 observations of three CO isotopologues plus the continuum towards PDS 70. Combining the inferred gas and dust distributions, the model results in a variable gas-to-dust ratio profile throughout the disk that spans two orders of magnitude within the first $130$ au and shows a step gradient towards the outer disk, which is consistent with the presence of a pressure maxima driven by planet-disk interactions. We find a gas density drop factor of ${\sim} 19$ at the location of the planet PDS 70 c with respect to the peak gas density at $75$ au. Combining this value with literature results on the hydrodynamics of planet-disk interactions, we find this gas gap depth to be consistent with independent planet mass estimates from infrared observations. Our findings point towards gas stirring processes taking place in the common gap due to the gravitational perturbation of both planets., Comment: Accepted for publication in A&A

Supernova remnants (SNRs) may regulate star formation in galaxies. For example, SNR-driven shocks may form new molecular gas or compress pre-existing clouds and trigger the formation of new stars. To test this scenario, we measure the deuteration of $N_2H^+$, $D_{frac}^{N_2H^+}$, a well-studied tracer of pre-stellar cores, across the Infrared Dark Cloud (IRDC) G034.77-00.55, known to be experiencing a shock interaction with the SNR W44. We use N$_2$H$^+$ and N$_2$D$^+$ J=1-0 single pointing observations obtained with the 30m antenna at the Instituto de Radioastronomia Millimetrica to infer $D_{frac}^{N_2H^+}$ toward five positions across the cloud, namely a massive core, different regions across the shock front, a dense clump and ambient gas. We find $D_{frac}^{N_2H^+}$ in the range 0.03-0.1, several orders of magnitude larger than the cosmic D/H ratio ($\sim$10$^{-5}$). Across the shock front, $D_{frac}^{N_2H^+}$ is enhanced by more than a factor of 2 ($D_{frac}^{N_2H^+}\sim$0.05-0.07) with respect to the ambient gas ($\leq$0.03) and similar to that measured generally in pre-stellar cores. Indeed, in the massive core and dense clump regions of this IRDC we measure $D_{frac}^{N_2H^+}$}$\sim$0.1. We find enhanced deuteration of $N_2H^+$ across the region of the shock, at a level that is enhanced with respect to regions of unperturbed gas. It is possible that this has been induced by shock compression, which would then be indirect evidence that the shock is triggering conditions for future star formation. However, since unperturbed dense regions also show elevated levels of deuteration, further, higher-resolution studies are needed to better understand the structure and kinematics of the deuterated material in the shock region, e.g., if it still in relatively diffuse form or already organised in a population of low-mass pre-stellar cores., Comment: Accepted for publication on A&A; 8 pages, 5 figures

The abundance and distribution of ice in protoplanetary disks (PPD) is critical to understand the linkage between the composition of circumstellar matter and the composition of exoplanets. Edge-on PPDs are a useful tool to constrain such ice composition and its location in the disk, as ice spectral signatures can be observed in absorption against the continuum emission arising from the warmer central disk regions. The aim of this work is to model ice absorption features in PPDs and determine how well the abundance of the main ice species across the disk can be determined within the uncertainty of the physical parameter space. The edge-on PPD around HH 48 NE, a target of the JWST ERS program IceAge, is used as a reference system. We use RADMC-3D to raytrace the mid-infrared continuum. Using a constant parameterized ice abundance, ice opacities are added to the dust opacity in regions wherever the disk is cold enough for the main carbon, oxygen and nitrogen carriers to freeze out. The global abundance of the main ice carriers in HH 48 NE can be determined within a factor of 3, when taking the uncertainty of the physical parameters into account. Ice features in PPDs can be saturated at an optical depth <1, due to local saturation. Spatially observed ice optical depths cannot be directly related to column densities due to radiative transfer effects. Vertical snowlines will not be a clear transition due to the radially increasing height of the snowsurface, but their location may be constrained from observations using radiative transfer modeling. Radial snowlines are not really accesible. Not only the ice abundance, but also inclination, settling, grain size distribution and disk mass have strong impact on the observed ice absorption features in disks. Relative changes in ice abundance can be inferred from observations only if the source structure is well constrained, Comment: 18 pages, 12 figures, accepted for publication in Astronomy & Astrophysics

Context. Observations of edge-on disks are an important tool for constraining general protoplanetary disk properties that cannot be determined in any other way. However, most radiative transfer models cannot simultaneously reproduce the spectral energy distributions (SEDs) and resolved scattered light and submillimeter observations of these systems, due to the differences in geometry and dust properties at different wavelengths. Aims. We simultaneously constrain the geometry of the edge-on protoplanetary disk HH 48 NE and the characteristics of the host star. HH 48 NE is part of the JWST early release science program Ice Age. This work serves as a stepping stone towards a better understanding of the disk physical structure and icy chemistry in this particular source. This kind of modeling lays the groundwork for studying other edge-on sources to be observed with the JWST. Methods. We fit a parameterized dust model to HH 48 NE by coupling the radiative transfer code RADMC-3D and an MCMC framework. The dust structure was fitted independently to a compiled SED, a scattered light image at 0.8 ${\mu}$m and an ALMA dust continuum observation at 890 ${\mu}$m. Results. We find that 90% of the dust mass in HH 48 NE is settled to the disk midplane, less than in average disks, and that the atmospheric layers of the disk contain exclusively large grains (0.3-10 ${\mu}$m). The exclusion of small grains in the upper atmosphere likely has important consequences for the chemistry due to the deep penetration of high-energy photons. The addition of a relatively large cavity (ca. 50 au in radius) is necessary to explain the strong mid-infrared emission, and to fit the scattered light and continuum observations simultaneously., Comment: 16 pages, 8 figures, accepted for publication in Astronomy & Astrophysics

Context: Core accretion models of massive star formation require the existence of massive, starless cores within molecular clouds. Yet, only a small number of candidates for such truly massive, monolithic cores are currently known. Aims: Here we analyse a massive core in the well-studied infrared-dark cloud (IRDC) called the 'dragon cloud' (also known as G028.37+00.07 or 'Cloud C'). This core (C2c1) sits at the end of a chain of a roughly equally spaced actively star-forming cores near the centre of the IRDC. Methods: We present new high-angular resolution 1 mm ALMA dust continuum and molecular line observations of the massive core. Results: The high-angular resolution observations show that this region fragments into two cores C2c1a and C2c1b, which retain significant background-subtracted masses of 23 Msun and 2 Msun (31 Msun and 6 Msun without background subtraction), respectively. The cores do not appear to fragment further on the scales of our highest angular resolution images (0.200 arcsec, 0.005 pc ~ 1000 AU). We find that these cores are very dense (nH2 > 10^6 cm-3) and have only trans-sonic non-thermal motions (Ms ~ 1). Together the mass, density and internal motions imply a virial parameter of < 1, which suggests the cores are gravitationally unstable, unless supported by strong magnetic fields with strengths of ~ 1 - 10 mG. From CO line observations, we find that there is tentative evidence for a weak molecular outflow towards the lower-mass core, and yet the more massive core remains devoid of any star formation indicators. Conclusions: We present evidence for the existence of a massive, pre-stellar core, which has implications for theories of massive star formation. This source warrants follow-up higher-angular-resolution observations to further assess its monolithic and pre-stellar nature., Comment: 8+4 pages, 4+2 Figures, 2 Tables, accepted for publication in Astronomy & Astrophysics

Fast radio bursts (FRBs) are millisecond-timescale radio transients, the origins of which are predominantly extragalactic and likely involve highly magnetized compact objects. FRBs undergo multipath propagation, or scattering, from electron density fluctuations on sub-parsec scales in ionized gas along the line-of-sight. Scattering observations have located plasma structures within FRB host galaxies, probed Galactic and extragalactic turbulence, and constrained FRB redshifts. Scattering also inhibits FRB detection and biases the observed FRB population. We report the detection of scattering times from the repeating FRB 20190520B that vary by up to a factor of two or more on minutes to days-long timescales. In one notable case, the scattering time varied from $7.9\pm0.4$ ms to less than 3.1 ms ($95\%$ confidence) over 2.9 minutes at 1.45 GHz. The scattering times appear to be uncorrelated between bursts or with dispersion and rotation measure variations. Scattering variations are attributable to dynamic, inhomogeneous plasma in the circumsource medium, and analogous variations have been observed from the Crab pulsar. Under such circumstances, the frequency dependence of scattering can deviate from the typical power-law used to measure scattering. Similar variations may therefore be detectable from other FRBs, even those with inconspicuous scattering, providing a unique probe of small-scale processes within FRB environments., Comment: 10 pages, 9 figures, accepted to MNRAS

The environment of actively repeating fast radio bursts (FRBs) has been shown to be complex and varying. The recently localized FRB 20190520B is extremely active, has the largest confirmed host dispersion measure, and is only the second FRB source associated with a compact, persistent radio source (PRS). The main tracer of the magneto-ionic environments is the rotation measure (RM), a path-integral of the line-of-sight component of magnetic field strength (B) and electron density, which does not allow a direct probe of the B-field configuration. Here we report direct evidence for a B-field reversal based on the observed sign change and extreme variation of FRB 20190520B's RM, which changed from $\sim10000$ rad m$^{-2}$ to $\sim-16000$ rad m$^{-2}$ between June 2021 and January 2022. Such extreme RM reversal has never been observed before in any FRB nor in any astronomical object. The implied short-term change of the B-field configuration in or around the FRB could be due to the vicinity of massive black holes, or a magnetized companion star in binary systems, or a young supernova remnant along the line of sight., Comment: It was merged into arXiv:2202.11112. Science380,599-603(2023). DOI:10.1126/science.abo6526

The repeating FRB 20190520B is localized to a galaxy at $z=0.241$, much closer than expected given its dispersion measure $\rm DM=1205\pm4\ pc\ cm^{-3}$. Here we assess implications of the large DM and scattering observed from FRB 20190520B for the host galaxy's plasma properties. A sample of 75 bursts detected with the Five-hundred-meter Aperture Spherical radio Telescope shows scattering on two scales: a mean temporal delay $\tau(1.41\ {\rm GHz})=10.9\pm1.5$ ms, which is attributed to the host galaxy, and a mean scintillation bandwidth $\nu_{\rm d}(1.41\ {\rm GHz})=0.21\pm0.01$ MHz, which is attributed to the Milky Way. Balmer line measurements for the host imply an H$\alpha$ emission measure (galaxy frame) $\rm EM_s=620$ pc cm$^{-6} \times (T/10^4 {\rm K})^{0.9}$, implying $\rm DM_{\rm H\alpha}$ of order the value inferred from the FRB DM budget, $\rm DM_h=1121^{+89}_{-138}$ pc cm$^{-3}$ for plasma temperatures greater than the typical value $10^4$ K. Combining $\tau$ and $\rm DM_h$ yields a nominal constraint on the scattering amplification from the host galaxy $\tilde{F} G=1.5^{+0.8}_{-0.3}$ (pc$^2$ km)$^{-1/3}$, where $\tilde{F}$ describes turbulent density fluctuations and $G$ represents the geometric leverage to scattering that depends on the location of the scattering material. For a two-screen scattering geometry where $\tau$ arises from the host galaxy and $\Delta \nu_{\rm d}$ from the Milky Way, the implied distance between the FRB source and dominant scattering material is $\lesssim100$ pc. The host galaxy scattering and DM contributions support a novel technique for estimating FRB redshifts using the $\tau-\rm DM$ relation, and are consistent with previous findings that scattering of localized FRBs is largely dominated by plasma within host galaxies and the Milky Way., Comment: 14 pages, 9 figures, accepted to ApJ

Supernova remnants (SNRs) contribute to regulate the star formation efficiency and evolution of galaxies. As they expand into the interstellar medium (ISM), they transfer vast amounts of energy and momentum that displace, compress and heat the surrounding material. Despite the extensive work in galaxy evolution models, it remains to be observationally validated to what extent the molecular ISM is affected by the interaction with SNRs. We use the first results of the ESO-ARO Public Spectroscopic Survey SHREC, to investigate the shock interaction between the SNR IC443 and the nearby molecular clump G. We use high sensitivity SiO(2-1) and H$^{13}$CO$^+$(1-0) maps obtained by SHREC together with SiO(1-0) observations obtained with the 40m telescope at the Yebes Observatory. We find that the bulk of the SiO emission is arising from the ongoing shock interaction between IC443 and clump G. The shocked gas shows a well ordered kinematic structure, with velocities blue-shifted with respect to the central velocity of the SNR, similar to what observed toward other SNR-cloud interaction sites. The shock compression enhances the molecular gas density, n(H$_2$), up to $>$10$^5$ cm$^{-3}$, a factor of >10 higher than the ambient gas density and similar to values required to ignite star formation. Finally, we estimate that up to 50\% of the momentum injected by IC443 is transferred to the interacting molecular material. Therefore the molecular ISM may represent an important momentum carrier in sites of SNR-cloud interactions., Comment: Accepted for Publication in MNRAS; 11 pages, 7 figures

The first Fast Radio Burst (FRB) to be precisely localized was associated with a luminous persistent radio source (PRS). Recently, a second FRB/PRS association was discovered for another repeating source of FRBs. However, it is not clear what makes FRBs or PRS or how they are related. We compile FRB and PRS properties to consider the population of FRB/PRS sources. We suggest a practical definition for PRS as FRB associations with luminosity greater than $10^{29}$ erg s$^{-1}$ Hz$^{-1}$ that is not attributed to star-formation activity in the host galaxy. We model the probability distribution of the fraction of FRBs with PRS for repeaters and non-repeaters, showing there is not yet evidence for repeaters to be preferentially associated with PRS. We discuss how FRB/PRS sources may be distinguished by the combination of active repetition and an excess dispersion measure local to the FRB environment. We use CHIME/FRB event statistics to bound the mean per-source repetition rate of FRBs to be between 25 and 440 yr$^{-1}$. We use this to provide a bound on the density of FRB-emitting sources in the local universe of between $2.2\times10^2$ and $5.2\times10^4$ Gpc$^{-3}$ assuming a pulsar-like beam width for FRB emission. This density implies that PRS may comprise as much as 1\% of compact, luminous radio sources detected in the local universe. The cosmic density and phenomenology of PRS are similar to that of the newly-discovered, off-nuclear "wandering" AGN. We argue that it is likely that some PRS have already been detected and misidentified as AGN., Comment: Submitted. 12 pages, 3 figures

The dispersive sweep of fast radio bursts (FRBs) has been used to probe the ionized baryon content of the intergalactic medium, which is assumed to dominate the total extragalactic dispersion. While the host galaxy contributions to dispersion measure (DM) appear to be small for most FRBs, in at least one case there is evidence for an extreme magneto-ionic local environment and a compact persistent radio source. Here we report the detection and localization of the repeating FRB 20190520B, which is co-located with a compact, persistent radio source and associated with a dwarf host galaxy of high specific star formation rate at a redshift $z=0.241\pm0.001$. The estimated host galaxy DM $\approx 903^{+72}_{-111}$ pc cm$^{-3}$, nearly an order of magnitude higher than the average of FRB host galaxies, far exceeds the DM contribution of the intergalactic medium. Caution is thus warranted in inferring redshifts for FRBs without accurate host galaxy identifications., Comment: Accepted, Version 3

APOGEE is a high-resolution (R sim 22,000), near-infrared, multi-epoch, spectroscopic survey of the Milky Way. The second generation of the APOGEE project, APOGEE-2, includes an expansion of the survey to the Southern Hemisphere called APOGEE-2S. This expansion enabled APOGEE to perform a fully panoramic mapping of all the main regions of the Milky Way; in particular, by operating in the H-band, APOGEE is uniquely able to probe the dust-hidden inner regions of the Milky Way that are best accessed from the Southern Hemisphere. In this paper we present the targeting strategy of APOGEE-2S, with special attention to documenting modifications to the original, previously published plan. The motivation for these changes is explained as well as an assessment of their effectiveness in achieving their intended scientific objective. In anticipation of this being the last paper detailing APOGEE targeting, we present an accounting of all such information complete through the end of the APOGEE-2S project; this includes several main survey programs dedicated to exploration of major stellar populations and regions of the Milky Way, as well as a full list of programs contributing to the APOGEE database through allocations of observing time by the Chilean National Time Allocation Committee (CNTAC) and the Carnegie Institution for Science (CIS). This work was presented along with a companion article, R. Beaton et al. (submitted; AAS29028), presenting the final target selection strategy adopted for APOGEE-2 in the Northern Hemisphere., Comment: 43 pages; 5 figures; 6 Tables; 1 Appendix; Submitted to Journal and Under Review; Posting to accompany papers using the SDSS-IV/APOGEE-2 Data Release 17 scheduled for December 2021

We present new radio and optical data, including very long baseline interferometry, as well as archival data analysis, for the luminous decades-long radio transient FIRST J141918.9+394036. The radio data reveal a synchrotron self-absorption peak around 0.3 GHz and a radius of around 1.3 mas (0.5 pc) 26 years post-discovery, indicating a blastwave energy $\sim5 \times 10^{50}$ erg. The optical spectrum shows a broad [OIII]$\lambda$4959,5007 emission-line that may indicate collisional-excitation in the host galaxy, but its association with the transient cannot be ruled out. The properties of the host galaxy are suggestive of a massive stellar progenitor that formed at low metallicity. Based on the radio light curve, blastwave velocity, energetics, nature of the host galaxy and transient rates we find that the properties of FIRST J1419+39 are most consistent with long gamma-ray burst (LGRB) afterglows. Other classes of (optically-discovered) stellar explosions as well as neutron star mergers are disfavored, and invoking any exotic scenario may not be necessary. It is therefore likely that FIRST J1419+39 is an off-axis LGRB afterglow (as suggested by Law et al. and Marcote et al.), and under this premise the inverse beaming fraction is found to be $f_b^{-1}\simeq280^{+700}_{-200}$, corresponding to an average jet half-opening angle $<\theta_j>\simeq5^{+4}_{-2}$ degrees (68% confidence), consistent with previous estimates. From the volumetric rate we predict that surveys with the VLA, ASKAP and MeerKAT will find a handful of FIRST J1419+39-like events over the coming years., Comment: 23 pages, 7 figures, 2 tables. Accepted in ApJ

Fast radio bursts (FRBs) are exceptionally luminous flashes of unknown physical origin, reaching us from other galaxies (Petroff et al. 2019). Most FRBs have only ever been seen once, while others flash repeatedly, though sporadically (Spitler et al. 2016, CHIME/FRB Collaboration et al. 2021). Many models invoke magnetically powered neutron stars (magnetars) as the engines producing FRB emission (Margalit & Metzger 2018, CHIME/FRB Collaboration et al. 2020). Recently, CHIME/FRB announced the discovery (Bhardwaj et al. 2021) of the repeating FRB 20200120E, coming from the direction of the nearby grand design spiral galaxy M81. Four potential counterparts at other observing wavelengths were identified (Bhardwaj et al. 2021) but no definitive association with these sources, or M81, could be made. Here we report an extremely precise localisation of FRB 20200120E, which allows us to associate it with a globular cluster (GC) in the M81 galactic system and to place it ~2pc offset from the optical center of light of the GC. This confirms (Bhardwaj et al. 2021) that FRB 20200120E is 40 times closer than any other known extragalactic FRB. Because such GCs host old stellar populations, this association strongly challenges FRB models that invoke young magnetars formed in a core-collapse supernova as powering FRB emission. We propose, instead, that FRB 20200120E is a highly magnetised neutron star formed via either accretion-induced collapse of a white dwarf or via merger of compact stars in a binary system (Margalit et al. 2019). Alternative scenarios involving compact binary systems, efficiently formed inside globular clusters, could also be responsible for the observed bursts., Comment: Submitted. Comments welcome

Infrared dark clouds (IRDCs) are potential hosts of the elusive early phases of high-mass star formation (HMSF). Here we conduct an in-depth analysis of the fragmentation properties of a sample of 10 IRDCs, which have been highlighted as some of the best candidates to study HMSF within the Milky Way. To do so, we have obtained a set of large mosaics covering these IRDCs with ALMA at band 3 (or 3mm). These observations have a high angular resolution (~3arcsec or ~0.05pc), and high continuum and spectral line sensitivity (~0.15mJy/beam and ~0.2K per 0.1km/s channel at the N2H+(1-0) transition). From the dust continuum emission, we identify 96 cores ranging from low- to high-mass (M = 3.4 to 50.9Msun) that are gravitationally bound (alpha_vir = 0.3 to 1.3) and which would require magnetic field strengths of B = 0.3 to 1.0mG to be in virial equilibrium. We combine these results with a homogenised catalogue of literature cores to recover the hierarchical structure within these clouds over four orders of magnitude in spatial scale (0.01pc to 10pc). Using supplementary observations at an even higher angular resolution, we find that the smallest fragments (<0.02pc) within this hierarchy do not currently have the mass and/or the density required to form high-mass stars. Nonetheless, the new ALMA observations presented in this paper have facilitated the identification of 19 (6 quiescent and 13 star-forming) cores that retain >16Msun without further fragmentation. These high-mass cores contain trans-sonic non-thermal motions, are kinematically sub-virial, and require moderate magnetic field strengths for support against collapse. The identification of these potential sites of high-mass star formation represents a key step in allowing us to test the predictions from high-mass star and cluster formation theories., Comment: 21(+6 in appendix) pages, 12(+3) figures, 3(+3) tables. Machine-readable versions of Table A1, A2 and A3 are available online. Accepted for publication in MNRAS

The 14N/15N ratio in molecules exhibits a large variation in star-forming regions, especially when measured from N2H+ isotopologues. However, there are only a few studies performed at high-angular resolution. We present the first interferometric survey of the 14N/15N ratio in N2H+ obtained with the Atacama Large Millimeter Array towards four infrared-dark clouds harbouring 3~mm continuum cores associated with different physical properties. We detect N15NH+ (1-0) in about 20-40% of the cores, depending on the host cloud. The 14N/15N values measured towards the millimeter continuum cores range from a minimum of 80 up to a maximum of 400. The spread of values is narrower than that found in any previous single-dish survey of high-mass star-forming regions, and than that obtained using the total power data only. This suggests that the 14N/15N ratio is on average higher in the diffuse gaseous envelope of the cores, and stresses the need for high-angular resolution maps to measure correctly the 14N/15N ratio in dense cores embedded in IRDCs. The average 14N/15N ratio of 210 is also lower than the interstellar value at the Galactocentric distance of the clouds (300-330), although the sensitivity of our observations does not allow us to unveil 14N/15N ratios higher than 400. No clear trend is found between the 14N/15N ratio and the core physical properties. We find only a tentative positive trend between 14N/15N and the H2 column density. However, firmer conclusions can be drawn only with higher sensitivity measurements., Comment: 19 pages, 7 figures, 6 tables, 1 appendix, accepted for publication in MNRAS

[Abridged] Filaments and hubs have received special attention recently thanks to studies showing their role in star formation. While the column density and velocity structures of both filaments and hubs have been studied, their magnetic fields (B-field) are not yet characterized. We aim to understand the role of the B-field in the dynamical evolution of the NGC 6334 hub-filament network. We present new observations of the dust polarized emission at 850$\mu$m towards NGC 6334 obtained with the JCMT/POL-2. We study the distribution and dispersion of the polarized intensity ($PI$), the polarization fraction ($PF$), and the B-field angle ($\theta_{B}$). We derive the power spectrum of the intensity and $\theta_{B}$ along the ridge crest. Our analyses show a complex B-field structure when observed over the whole region ($\sim10$ pc), however, at smaller scales ($\sim1$ pc), $\theta_{B}$ varies coherently along the filaments. The observed power spectrum of $\theta_{B}$ can be well represented with a power law function with a slope $-1.33\pm0.23$, which is $\sim20\%$ shallower than that of $I$. This result is compatible with the properties of simulated filaments and may indicate the processes at play in the formation of filaments. $\theta_{B}$ rotates from being mostly perpendicular to the filament crests to mostly parallel as they merge with the hubs. This variation of $\theta_{B}$ may be tracing local velocity flows of matter in-falling onto the hubs. Our analysis suggests a variation of the energy balance along the crests of these filaments, from magnetically critical/supercritical at their far ends to magnetically subcritical near the hubs. We detect an increase of $PF$ towards the high-column density star cluster-forming hubs that may result from the increase of grain alignment efficiency due to stellar radiation from the newborn stars., Comment: Accepted for publication in Astronomy & Astrophysics

During the past decade the dynamical importance of magnetic fields in molecular clouds has been increasingly recognized, as observational evidence has accumulated. However, how a magnetic field affect star formation is still unclear. Typical star formation models still treat a magnetic fields as an isotropic pressure, ignoring the fundamental property of dynamically important magnetic fields: their direction. This study builds on our previous work which demonstrated how the mean magnetic field orientation relative to the global cloud elongation can affect cloud fragmentation. After the linear mass distribution reported earlier, we show here that the mass cumulative function (MCF) of a cloud is also regulated by the field orientation. A cloud elongated closer to the field direction tends to have a shallower MCF, in other words, a higher portion of the gas in high density. The evidence is consistent with our understanding of bimodal star formation efficiency discovered earlier, which is also correlated with the field orientations., Comment: Accepted by MNRAS on 11 August 2020

N132D is the brightest gamma-ray supernova remnant (SNR) in the Large Magellanic Cloud (LMC). We carried out $^{12}$CO($J$ = 1-0, 3-2) observations toward the SNR using the Atacama Large Millimeter/submillimeter Array (ALMA) and Atacama Submillimeter Telescope Experiment. We find diffuse CO emission not only at the southern edge of the SNR as previously known, but also inside the X-ray shell. We spatially resolved nine molecular clouds using ALMA with an angular resolution of $5''$, corresponding to a spatial resolution of $\sim$1 pc at the distance of the LMC. Typical cloud sizes and masses are $\sim$2.0 pc and $\sim$100 $M_\odot$, respectively. High-intensity ratios of CO $J$ = 3-2 / 1-0 $> 1.5$ are seen toward the molecular clouds, indicating that shock-heating has occurred. Spatially resolved X-ray spectroscopy reveals that thermal X-rays in the center of N132D are produced not only behind a molecular cloud, but also in front of it. Considering the absence of a thermal component associated with the forward shock towards one molecular cloud located along the line of sight to the center of the remnant, this suggests that this particular cloud is engulfed by shock waves and is positioned on the near side of remnant. If the hadronic process is the dominant contributor to the gamma-ray emission, the shock-engulfed clouds play a role as targets for cosmic-rays. We estimate the total energy of cosmic-ray protons accelerated in N132D to be $\sim$0.5-$3.8 \times 10^{49}$ erg as a conservative lower limit, which is similar to that observed in Galactic gamma-ray SNRs., Comment: 18 pages, 2 tables, 11 figures, accepted for publication in The Astrophysical Journal (ApJ)

We present the discovery and subarcsecond localization of a new Fast Radio Burst with the Karl G. Jansky Very Large Array and realfast search system. The FRB was discovered on 2019 June 14 with a dispersion measure of 959 pc/cm3. This is the highest DM of any localized FRB and its measured burst fluence of 0.6 Jy ms is less than nearly all other FRBs. The source is not detected to repeat in 15 hours of VLA observing and 153 hours of CHIME/FRB observing. We describe a suite of statistical and data quality tests we used to verify the significance of the event and its localization precision. Follow-up optical/infrared photometry with Keck and Gemini associate the FRB to a pair of galaxies with $\rm{r}\sim23$ mag. The false-alarm rate for radio transients of this significance that are associated with a host galaxy is roughly $3\times10^{-4}\ \rm{hr}^{-1}$. The two putative host galaxies have similar photometric redshifts of $z_{\rm{phot}}\sim0.6$, but different colors and stellar masses. Comparing the host distance to that implied by the dispersion measure suggests a modest (~ 50 pc/cm3) electron column density associated with the FRB environment or host galaxy/galaxies., Comment: Submitted to AAS Journals (ApJ) and revised for referee comments

We present low-radio-frequency follow-up observations of AT 2017gfo, the electromagnetic counterpart of GW170817, which was the first binary neutron star merger to be detected by Advanced LIGO-Virgo. These data, with a central frequency of 144 MHz, were obtained with LOFAR, the Low-Frequency Array. The maximum elevation of the target is just 13.7 degrees when observed with LOFAR, making our observations particularly challenging to calibrate and significantly limiting the achievable sensitivity. On time-scales of 130-138 and 371-374 days after the merger event, we obtain 3$\sigma$ upper limits for the afterglow component of 6.6 and 19.5 mJy beam$^{-1}$, respectively. Using our best upper limit and previously published, contemporaneous higher-frequency radio data, we place a limit on any potential steepening of the radio spectrum between 610 and 144 MHz: the two-point spectral index $\alpha^{610}_{144} \gtrsim -2.5$. We also show that LOFAR can detect the afterglows of future binary neutron star merger events occurring at more favourable elevations., Comment: 9 pages, 2 figures, accepted for publication in MNRAS

Fast radio bursts (FRBs) are brief, bright, extragalactic radio flashes. Their physical origin remains unknown, but dozens of possible models have been postulated. Some FRB sources exhibit repeat bursts. Though over a hundred FRB sources have been discovered to date, only four have been localised and associated with a host galaxy, with just one of the four known to repeat. The properties of the host galaxies, and the local environments of FRBs, provide important clues about their physical origins. However, the first known repeating FRB has been localised to a low-metallicity, irregular dwarf galaxy, and the apparently non-repeating sources to higher-metallicity, massive elliptical or star-forming galaxies, suggesting that perhaps the repeating and apparently non-repeating sources could have distinct physical origins. Here we report the precise localisation of a second repeating FRB source, FRB 180916.J0158+65, to a star-forming region in a nearby (redshift $z = 0.0337 \pm 0.0002$) massive spiral galaxy, whose properties and proximity distinguish it from all known hosts. The lack of both a comparably luminous persistent radio counterpart and a high Faraday rotation measure further distinguish the local environment of FRB 180916.J0158+65 from that of the one previously localised repeating FRB source, FRB 121102. This demonstrates that repeating FRBs have a wide range of luminosities, and originate from diverse host galaxies and local environments., Comment: 61 pages, 11 figures, 3 tables. Published in Nature

We present results of a search for late-time radio emission and Fast Radio Bursts (FRBs) from a sample of type-I superluminous supernovae (SLSNe-I). We used the Karl G. Jansky Very Large Array to observe ten SLSN-I more than 5 years old at a frequency of 3 GHz. We searched fast-sampled visibilities for FRBs and used the same data to perform a deep imaging search for late-time radio emission expected in models of magnetar-powered supernovae. No FRBs were found. One SLSN-I, PTF10hgi, is detected in deep imaging, corresponding to a luminosity of $1.2\times10^{28}$ erg s$^{-1}$. This luminosity, considered with the recent 6 GHz detection of PTF10hgi in Eftekhari et al (2019), supports the interpretation that it is powered by a young, fast-spinning ($\sim$ ms spin period) magnetar with $\sim$ 15 Msun of partially ionized ejecta. Broadly, our observations are most consistent with SLSNe-I being powered by neutron stars with fast spin periods, although most require more free-free absorption than is inferred for PTF10hgi. We predict that radio observations at higher frequencies or in the near future will detect these systems and begin constraining properties of the young pulsars and their birth environments., Comment: Accepted to ApJ

We present the DSA-2000: a world-leading radio survey telescope and multi-messenger discovery engine for the next decade. The array will be the first true radio camera, outputting science-ready image data over the 0.7 - 2 GHz frequency range with a spatial resolution of 3.5 arcsec. With 2000 x 5 m dishes, the DSA-2000 will have an equivalent point-source sensitivity to SKA1-mid, but with ten times the survey speed. The DSA-2000 is envisaged as an all-sky survey instrument complementary to the ngVLA, and as a counterpart to the LSST (optical), SPHEREx (near-infrared) and SRG/eROSITA (X-ray) all-sky surveys. Over a five-year prime phase, the DSA-2000 will image the entire sky above declination -30 degrees every four months, detecting > 1 unique billion radio sources in a combined full-Stokes sky map with 500 nJy/beam rms noise. This all-sky survey will be complemented by intermediate and deep surveys, as well as spectral and polarization image cubes. The array will be a cornerstone for multi-messenger science, serving as the principal instrument for the US pulsar timing array community, and by searching for radio afterglows of compact object mergers detected by LIGO and Virgo. The array will simultaneously detect and localize ~10,000 fast radio bursts each year, realizing their ultimate use as a cosmological tool. The DSA-2000 will be proposed to the NSF Mid-Scale Research Infrastructure-2 program with a view to first light in 2026, Comment: Submitted as a Project White Paper for the Astro2020 Decadal Survey (12 pages, 7 figures, 2 table)

The Very Large Array Sky Survey (VLASS) is a synoptic, all-sky radio sky survey with a unique combination of high angular resolution ($\approx$2.5"), sensitivity (a 1$\sigma$ goal of 70 $\mu$Jy/beam in the coadded data), full linear Stokes polarimetry, time domain coverage, and wide bandwidth (2-4 GHz). The first observations began in September 2017, and observing for the survey will finish in 2024. VLASS will use approximately 5500 hours of time on the Karl G. Jansky Very Large Array (VLA) to cover the whole sky visible to the VLA (Declination $>-40^{\circ}$), a total of 33,885 deg$^2$. The data will be taken in three epochs to allow the discovery of variable and transient radio sources. The survey is designed to engage radio astronomy experts, multi-wavelength astronomers, and citizen scientists alike. By utilizing an "on the fly" interferometry mode, the observing overheads are much reduced compared to a conventional pointed survey. In this paper, we present the science case and observational strategy for the survey, and also results from early survey observations., Comment: 34 pages, accepted by PASP (modified from prior version to address referee's and coauthor comments). (v2) Minor fixes to author list

Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New methodologies derived from advances in statistics, computer science, and machine learning are beginning to be employed in sophisticated investigations that are not only bringing forth new discoveries, but are placing them on a solid footing. Progress in wide-field sky surveys, interferometric imaging, precision cosmology, exoplanet detection and characterization, and many subfields of stellar, Galactic and extragalactic astronomy, has resulted in complex data analysis challenges that must be solved to perform scientific inference. Research in astrostatistics and astroinformatics will be necessary to develop the state-of-the-art methodology needed in astronomy. Overcoming these challenges requires dedicated, interdisciplinary research. We recommend: (1) increasing funding for interdisciplinary projects in astrostatistics and astroinformatics; (2) dedicating space and time at conferences for interdisciplinary research and promotion; (3) developing sustainable funding for long-term astrostatisics appointments; and (4) funding infrastructure development for data archives and archive support, state-of-the-art algorithms, and efficient computing., Comment: Submitted to the Astro2020 Decadal Survey call for science white papers

FRB 121102 is the only known repeating fast radio burst source. Here we analyze a wide-frequency-range (1-8 GHz) sample of high-signal-to-noise, coherently dedispersed bursts detected using the Arecibo and Green Bank telescopes. These bursts reveal complex time-frequency structures that include sub-bursts with finite bandwidths. The frequency-dependent burst structure complicates the determination of a dispersion measure (DM); we argue that it is appropriate to use a DM metric that maximizes frequency-averaged pulse structure, as opposed to peak signal-to-noise, and find DM = 560.57 +/- 0.07 pc/cc at MJD 57644. After correcting for dispersive delay, we find that the sub-bursts have characteristic frequencies that typically drift lower at later times in the total burst envelope. In the 1.1-1.7 GHz band, the ~ 0.5-1-ms sub-bursts have typical bandwidths ranging from 100-400 MHz, and a characteristic drift rate of ~ 200 MHz/ms towards lower frequencies. At higher radio frequencies, the sub-burst bandwidths and drift rate are larger, on average. While these features could be intrinsic to the burst emission mechanism, they could also be imparted by propagation effects in the medium local to the source. Comparison of the burst DMs with previous values in the literature suggests an increase of Delta(DM) ~ 1-3 pc/cc in 4 years, though this could be a stochastic variation as opposed to a secular trend. This implies changes in the local medium or an additional source of frequency-dependent delay. Overall, the results are consistent with previously proposed scenarios in which FRB 121102 is embedded in a dense nebula., Comment: Submitted to ApJ; comments welcome

We present a non-relativistic theory of quantum radiation generated by shaking a two-level atom in vacuum. Such radiation has the same origin of photon emission in dynamical Casimir effect. By performing a time-dependent "dressing" transformation to the Hamiltonian, we derive an interaction term that governs the radiation. In particular, we show that photon pairs can be generated, not only by shaking the position of the atom, but also by changing the internal states of the atom. As applications of our theory, we calculate the emission rate from an oscillating atom, and the multi-photon state generated in a single-photon scattering process., Comment: 7 pages

We present the discovery of a slowly-evolving, extragalactic radio transient, FIRST J141918.9+394036, identified by comparing a catalog of radio sources in nearby galaxies against new observations from the Very Large Array Sky Survey. Analysis of other archival data shows that FIRST J141918.9+394036 faded by a factor of ~50 over 23 years, from a flux of ~26 mJy at 1.4 GHz in 1993 to an upper limit of 0.4 mJy at 3 GHz in 2017. FIRST J141918.9+394036 is likely associated with the small star-forming galaxy SDSS J141918.81+394035.8 at a redshift z=0.01957 (d=87 Mpc), which implies a peak luminosity $\nu L_\nu \gtrsim 3\times10^{38}$ erg s$^{-1}$. If interpreted as an isotropic synchrotron blast wave, the source requires an explosion of kinetic energy ~10^{51} erg some time prior to our first detection in late 1993. This explosion could plausibly be associated with a long gamma-ray burst (GRB) or the merger of two neutron stars. Alternatively, FIRST J141918.9+394036 could be the nebula of a newly-born magnetar. The radio discovery of any of these phenomena would be unprecedented. Joint consideration of the event light curve, host galaxy, lack of a counterpart gamma-ray burst, and volumetric rate suggests that FIRST J141918.9+394036 is the afterglow of an off-axis (`orphan') long GRB. The long time baseline of this event offers the best available constraint in afterglow evolution as the bulk of shock-accelerated electrons become non-relativistic. The proximity, age, and precise localization of FIRST J141918.9+394036 make it a key object for understanding the aftermath of rare classes of stellar explosion., Comment: Submitted to AAS Journals (14 pages, 3 figures)

With their wide fields of view and often relatively long coverage of any position in the sky in imaging survey mode, modern radio telescopes provide a data stream that is naturally suited to searching for rare transients. However, Radio Frequency Interference (RFI) can show up in the data stream in similar ways to such transients, and thus the normal pre-treatment of filtering RFI (flagging) may also remove astrophysical transients from the data stream before imaging. In this paper we investigate how standard flagging affects the detectability of such transients by examining the case of transient detection in an observing mode used for Low Frequency Array (LOFAR; \citep{LOFAR}) surveys. We quantify the fluence range of transients that would be detected, and the reduction of their SNR due to partial flagging. We find that transients with a duration close to the integration sampling time, as well as bright transients with durations on the order of tens of seconds, are completely flagged. For longer transients on the order of several tens of seconds to minutes, the flagging effects are not as severe, although part of the signal is lost. For these transients, we present a modified flagging strategy which mitigates the effect of flagging on transient signals. We also present a script which uses the differences between the two strategies, and known differences between transient RFI and astrophysical transients, to notify the observer when a potential transient is in the data stream., Comment: 31 pages, 6 figures, accepted for publication in Astronomy and Computing

We report the first detections of the repeating fast radio burst source FRB 121102 above 5.2 GHz. Observations were performed using the 4$-$8 GHz receiver of the Robert C. Byrd Green Bank Telescope with the Breakthrough Listen digital backend. We present the spectral, temporal and polarization properties of 21 bursts detected within the first 60 minutes of a total 6-hour observations. These observations comprise the highest burst density yet reported in the literature, with 18 bursts being detected in the first 30 minutes. A few bursts clearly show temporal sub-structures with distinct spectral properties. These sub-structures superimpose to provide enhanced peak signal-to-noise ratio at higher trial dispersion measures. Broad features occur in $\sim 1$ GHz wide subbands that typically differ in peak frequency between bursts within the band. Finer-scale structures ($\sim 10-50$ MHz) within these bursts are consistent with that expected from Galactic diffractive interstellar scintillation. The bursts exhibit nearly 100% linear polarization, and a large average rotation measure of 9.359$\pm$0.012 $\times$ 10$^{\rm 4}$ rad m$^{\rm -2}$ (in the observer's frame). No circular polarization was found for any burst. We measure an approximately constant polarization position angle in the 13 brightest bursts. The peak flux densities of the reported bursts have average values (0.2$\pm$0.1 Jy), similar to those seen at lower frequencies ($<3$ GHz), while the average burst widths (0.64$\pm$0.46 ms) are relatively narrower., Comment: Accepted for publication in ApJ

We present SMA and NOEMA observations of the host galaxy of FRB 121102 in the CO 3-2 and 1-0 transitions, respectively. We do not detect emission from either transition. We set $3\sigma$ upper limits to the CO luminosity $L_{CO} < 2.5 \times 10^7\,{\rm K\,km\,s}^{-1} {\, \rm pc^{-2}}$ for CO 3-2 and $L_{CO} < 2.3 \times 10^9\, {\rm K\,km\,s}^{-1} {\, \rm pc^{-2}}$ for CO 1-0. For Milky-Way-like star formation properties, we set a $3\sigma$ upper limit on the $H_2$ mass of $2.5 \times 10^8 \rm\ M_{\odot}$, slightly less than the predictions for the $H_2$ mass based on the star formation rate. The true constraint on the $H_2$ mass may be significantly higher, however, because of the reduction in CO luminosity that is common forlow-metallicity dwarf galaxies like the FRB host galaxy. These results demonstrate the challenge of identifying the nature of FRB progenitors through study of the host galaxy molecular gas. We also place a limit of 42 $\mu$Jy ($3\sigma$) on the continuum flux density of the persistent radio source at 97 GHz, consistent with a power-law extrapolation of the low frequency spectrum, which may arise from an AGN or other nonthermal source., Comment: Accepted for publication in AJ

PySE is a Python software package for finding and measuring sources in radio telescope images. The software was designed to detect sources in the LOFAR telescope images, but can be used with images from other radio telescopes as well. We introduce the LOFAR Telescope, the context within which PySE was developed, the design of PySE, and describe how it is used. Detailed experiments on the validation and testing of PySE are then presented, along with results of performance testing. We discuss some of the current issues with the algorithms implemented in PySE and their inter- action with LOFAR images, concluding with the current status of PySE and its future development., Comment: 26 pages, 6 figures

We present LOFAR high-band data over the frequency range 115-189 MHz for the X-ray binary SS 433, obtained in an observing campaign from 2013 February - 2014 May. Our results include a deep, wide-field map, allowing a detailed view of the surrounding supernova remnant W 50 at low radio frequencies, as well as a light curve for SS 433 determined from shorter monitoring runs. The complex morphology of W 50 is in excellent agreement with previously published higher-frequency maps; we find additional evidence for a spectral turnover in the eastern wing, potentially due to foreground free-free absorption. Furthermore, SS 433 is tentatively variable at 150 MHz, with both a debiased modulation index of 11 per cent and a $\chi^2$ probability of a flat light curve of $8.2 \times 10^{-3}$. By comparing the LOFAR flux densities with contemporaneous observations carried out at 4800 MHz with the RATAN-600 telescope, we suggest that an observed $\sim$0.5-1 Jy rise in the 150-MHz flux density may correspond to sustained flaring activity over a period of approximately six months at 4800 MHz. However, the increase is too large to be explained with a standard synchrotron bubble model. We also detect a wealth of structure along the nearby Galactic plane, including the most complete detection to date of the radio shell of the candidate supernova remnant G 38.7-1.4. This further demonstrates the potential of supernova remnant studies with the current generation of low-frequency radio telescopes., Comment: 19 pages, 8 figures, 4 tables, accepted for publication in MNRAS

Radio interferometers have the ability to precisely localize and better characterize the properties of sources. This ability is having a powerful impact on the study of fast radio transients, where a few milliseconds of data is enough to pinpoint a source at cosmological distances. However, recording interferometric data at millisecond cadence produces a terabyte-per-hour data stream that strains networks, computing systems, and archives. This challenge mirrors that of other domains of science, where the science scope is limited by the computational architecture as much as the physical processes at play. Here, we present a solution to this problem in the context of radio transients: realfast, a commensal, fast transient search system at the Jansky Very Large Array. Realfast uses a novel architecture to distribute fast-sampled interferometric data to a 32-node, 64-GPU cluster for real-time imaging and transient detection. By detecting transients in situ, we can trigger the recording of data for those rare, brief instants when the event occurs and reduce the recorded data volume by a factor of 1000. This makes it possible to commensally search a data stream that would otherwise be impossible to record. This system will search for millisecond transients in more than 1000 hours of data per year, potentially localizing several Fast Radio Bursts, pulsars, and other sources of impulsive radio emission. We describe the science scope for realfast, the system design, expected outcomes, and ways real-time analysis can help in other fields of astrophysics., Comment: Accepted to ApJS Special Issue on Data; 11 pages, 4 figures

Fast radio bursts (FRBs) are millisecond-duration, extragalactic radio flashes of unknown physical origin. FRB 121102, the only known repeating FRB source, has been localized to a star-forming region in a dwarf galaxy at redshift z = 0.193, and is spatially coincident with a compact, persistent radio source. The origin of the bursts, the nature of the persistent source, and the properties of the local environment are still debated. Here we present bursts that show ~100% linearly polarized emission at a very high and variable Faraday rotation measure in the source frame: RM_src = +1.46 x 10^5 rad m^-2 and +1.33 x 10^5 rad m^-2 at epochs separated by 7 months, in addition to narrow (< 30 mus) temporal structure. The large and variable rotation measure demonstrates that FRB 121102 is in an extreme and dynamic magneto-ionic environment, while the short burst durations argue for a neutron star origin. Such large rotation measures have, until now, only been observed in the vicinities of massive black holes (M_BH > 10^4 MSun). Indeed, the properties of the persistent radio source are compatible with those of a low-luminosity, accreting massive black hole. The bursts may thus come from a neutron star in such an environment. However, the observed properties may also be explainable in other models, such as a highly magnetized wind nebula or supernova remnant surrounding a young neutron star., Comment: Published in Nature: DOI: 10.1038/nature25149

We present results of the coordinated observing campaign that made the first subarcsecond localization of a Fast Radio Burst, FRB 121102. During this campaign, we made the first simultaneous detection of an FRB burst by multiple telescopes: the VLA at 3 GHz and the Arecibo Observatory at 1.4 GHz. Of the nine bursts detected by the Very Large Array at 3 GHz, four had simultaneous observing coverage at other observatories. We use multi-observatory constraints and modeling of bursts seen only at 3 GHz to confirm earlier results showing that burst spectra are not well modeled by a power law. We find that burst spectra are characterized by a ~500 MHz envelope and apparent radio energy as high as $10^{40}$ erg. We measure significant changes in the apparent dispersion between bursts that can be attributed to frequency-dependent profiles or some other intrinsic burst structure that adds a systematic error to the estimate of DM by up to 1%. We use FRB 121102 as a prototype of the FRB class to estimate a volumetric birth rate of FRB sources $R_{FRB} \approx 5x10^{-5}/N_r$ Mpc$^{-3}$ yr$^{-1}$, where $N_r$ is the number of bursts per source over its lifetime. This rate is broadly consistent with models of FRBs from young pulsars or magnetars born in superluminous supernovae or long gamma-ray bursts, if the typical FRB repeats on the order of thousands of times during its lifetime., Comment: 17 pages, 7 figures. Submitted to AAS Journals

We present optical, near- and mid-infrared imaging of the host galaxy of FRB 121102 with the Gemini North telescope, the Hubble Space Telescope and the Spitzer Space Telescope. The FRB 121102 host galaxy is resolved, revealing a bright star forming region located in the outskirts of the irregular, low-metallicity dwarf galaxy. The star forming region has a half-light radius of 0.68 kpc (0.20 arcsec), encompassing the projected location of the compact (<0.7 pc), persistent radio source that is associated with FRB 121102. The half-light diameter of the dwarf galaxy is 5 to 7 kpc, and broadband spectral energy distribution fitting indicates that it has a total stellar mass of M*~10^8 Msun. The metallicity of the host galaxy is low, 12+log10 ([O/H])=8.0+-0.1. The properties of the host galaxy of FRB 121102 are comparable to those of extreme emission line galaxies, also known to host hydrogen-poor superluminous supernovae and long-duration gamma-ray bursts. The projected location of FRB 121102 within the star forming region supports the proposed connection of FRBs with newly born neutron stars or magnetars., Comment: 6 pages, 2 figures, submitted to ApJL; v2 has correct r' band magnitude

We undertook coordinated campaigns with the Green Bank, Effelsberg, and Arecibo radio telescopes during Chandra X-ray Observatory and XMM-Newton observations of the repeating fast radio burst FRB 121102 to search for simultaneous radio and X-ray bursts. We find 12 radio bursts from FRB 121102 during 70 ks total of X-ray observations. We detect no X-ray photons at the times of radio bursts from FRB 121102 and further detect no X-ray bursts above the measured background at any time. We place a 5$\sigma$ upper limit of $3\times10^{-11}$ erg cm$^{-2}$ on the 0.5--10 keV fluence for X-ray bursts at the time of radio bursts for durations $<700$ ms, which corresponds to a burst energy of $4\times10^{45}$ erg at the measured distance of FRB 121102. We also place limits on the 0.5--10 keV fluence of $5\times10^{-10}$ erg cm$^{-2}$ and $1\times10^{-9}$ erg cm$^{-2}$ for bursts emitted at any time during the XMM-Newton and Chandra observations, respectively, assuming a typical X-ray burst duration of 5 ms. We analyze data from the Fermi Gamma-ray Space Telescope Gamma-ray Burst Monitor and place a 5$\sigma$ upper limit on the 10--100 keV fluence of $4\times10^{-9}$ erg cm$^{-2}$ ($5\times10^{47}$ erg at the distance of FRB 121102) for gamma-ray bursts at the time of radio bursts. We also present a deep search for a persistent X-ray source using all of the X-ray observations taken to date and place a 5$\sigma$ upper limit on the 0.5--10 keV flux of $4\times10^{-15}$ erg s$^{-1}$ cm$^{-2}$ ($3\times10^{41}$ erg~s$^{-1}$ at the distance of FRB 121102). We discuss these non-detections in the context of the host environment of FRB 121102 and of possible sources of fast radio bursts in general., Comment: 13 pages, 5 figures, published in ApJ