skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Spectral stacking of radio-interferometric data
Context.Mapping molecular line emission beyond the bright low-JCO transitions is still challenging in extragalactic studies, even with the latest generation of (sub-)millimetre interferometers, such as ALMA and NOEMA. Aims.We summarise and test a spectral stacking method that has been used in the literature to recover low-intensity molecular line emission, such as HCN(1−0), HCO+(1−0), and even fainter lines in external galaxies. The goal is to study the capabilities and limitations of the stacking technique when applied to imaged interferometric observations. Methods.The core idea of spectral stacking is to align spectra of the low S/N spectral lines to a known velocity field calculated from a higher S/N line expected to share the kinematics of the fainter line (e.g. CO(1−0) or 21 cm emission). Then these aligned spectra can be coherently averaged to produce potentially high S/N spectral stacks. Here we used imaged simulated interferometric and total power observations at different S/N levels, based on real CO observations. Results.For the combined interferometric and total power data, we find that the spectral stacking technique is capable of recovering the integrated intensities even at low S/N levels across most of the region where the high S/N prior is detected. However, when stacking interferometer-only data for low S/N emission, the stacks can miss up to 50% of the emission from the fainter line. Conclusions.A key result of this analysis is that the spectral stacking method is able to recover the true mean line intensities in low S/N cubes and to accurately measure the statistical significance of the recovered lines. To facilitate the application of this technique we provide a public Python package, called PYSTACKER.  more » « less
Award ID(s):
1653300
PAR ID:
10469143
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
A&A
Date Published:
Journal Name:
Astronomy & Astrophysics
Volume:
675
ISSN:
0004-6361
Page Range / eLocation ID:
A104
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Existing star-forming vs. active galactic nucleus (AGN) classification schemes using optical emission-line diagnostics mostly fail for low-metallicity and/or highly star-forming galaxies, missing AGN in typicalz∼ 0 dwarfs. To recover AGN in dwarfs with strong emission lines (SELs), we present a classification scheme optimizing the use of existing optical diagnostics. We use Sloan Digital Sky Survey emission-line catalogs overlapping the volume- and mass-limited REsolved Spectroscopy Of a Local VolumE (RESOLVE) and Environmental COntex (ECO) surveys to determine the AGN percentage in SEL dwarfs. Our photoionization grids show that the [Oiii]/Hβversus [Sii]/Hαdiagram (Siiplot) and [Oiii]/Hβversus [Oi]/Hαdiagram (Oiplot) are less metallicity sensitive and more successful in identifying dwarf AGN than the popular [Oiii]/Hβversus [Nii]/Hαdiagnostic (Niiplot or “BPT diagram”). We identify a new category of “star-forming AGN” (SF-AGN) classified as star-forming by the Niiplot but as AGN by the Siiand/or Oiplots. Including SF-AGN, we find thez∼ 0 AGN percentage in dwarfs with SELs to be ∼3%–16%, far exceeding most previous optical estimates (∼1%). The large range in our dwarf AGN percentage reflects differences in spectral fitting methodologies between catalogs. The highly complete nature of RESOLVE and ECO allows us to normalize strong emission-line galaxy statistics to the full galaxy population, reducing the dwarf AGN percentage to ∼0.6%–3.0%. The newly identified SF-AGN are mostly gas-rich dwarfs with halo mass <1011.5M, where highly efficient cosmic gas accretion is expected. Almost all SF-AGN also have low metallicities (Z≲ 0.4Z), demonstrating the advantage of our method. 
    more » « less
  2. Abstract M82 is an archetypal starburst galaxy in the local Universe. The central burst of star formation, thought to be triggered by M82's interaction with other members in the M81 group, is driving a multiphase galaxy-scale wind away from the plane of the disk that has been studied across the electromagnetic spectrum. Here, we present new velocity-resolved observations of the [Cii] 158μm line in the central disk and the southern outflow of M82 using the upGREAT instrument on board SOFIA. We also report the first detections of velocity-resolved (ΔV= 10 km s−1) [Cii] emission in the outflow of M82 at projected distances of ≈1–2 kpc south of the galaxy center. We compare the [Cii] line profiles to observations of CO and Hiand find that likely the majority (>55%) of the [Cii] emission in the outflow is associated with the neutral atomic medium. We find that the fraction of [Cii] actually outflowing from M82 is small compared to the bulk gas outside the midplane (which may be in a halo or tidal streamers), which has important implications for observations of [Cii] outflows at higher redshift. Finally, by comparing the observed ratio of the [Cii] and CO intensities to models of photodissociation regions, we estimate that the far-ultraviolet (FUV) radiation field in the disk is ∼103.5G0, in agreement with previous estimates. In the outflow, however, the FUV radiation field is 2–3 orders of magnitudes lower, which may explain the high fraction of [Cii] arising from the neutral medium in the wind. 
    more » « less
  3. We present Atacama Large Millimeter/submillimeter Array (ALMA) sub-kiloparsec- to kiloparsec-scale resolution observations of the [C II], CO (9–8), and OH+(11–01) lines along with their dust continuum emission toward the far-infrared (FIR) luminous quasar SDSS J231038.88+185519.7 atz = 6.0031, to study the interstellar medium distribution, the gas kinematics, and the quasar-host system dynamics. We decompose the intensity maps of the [C II] and CO (9–8) lines and the dust continuum with two-dimensional elliptical Sérsic models. The [C II] brightness follows a flat distribution with a Sérsic index of 0.59. The CO (9–8) line and the dust continuum can be fit with an unresolved nuclear component and an extended Sérsic component with a Sérsic index of ∼1, which may correspond to the emission from an active galactic nucleus dusty molecular torus and a quasar host galaxy, respectively. The different [C II] spatial distribution may be due to the effect of the high dust opacity, which increases the FIR background radiation on the [C II] line, especially in the galaxy center, significantly suppressing the [C II] emission profile. The dust temperature drops with distance from the center. The effective radius of the dust continuum is smaller than that of the line emission and the dust mass surface density, but is consistent with that of the star formation rate surface density. This may indicate that the dust emission is a less robust tracer of the dust and gas distribution but is a decent tracer of the obscured star formation activity. The OH+(11–01) line shows a P-Cygni profile with an absorption at ∼–400 km s−1, which may indicate an outflow with a neutral gas mass of (6.2 ± 1.2)×108 Malong the line of sight. We employed a three-dimensional tilted ring model to fit the [C II] and CO (9–8) data cubes. The two lines are both rotation dominated and trace identical disk geometries and gas motions. This suggest that the [C II] and CO (9–8) gas are coplanar and corotating in this quasar host galaxy. The consistent circular velocities measured with [C II] and CO (9–8) lines indicate that these two lines trace a similar gravitational potential. We decompose the circular rotation curve measured from the kinematic model fit to the [C II] line into four matter components (black hole, stars, gas, and dark matter). The quasar-starburst system is dominated by baryonic matter inside the central few kiloparsecs. We constrain the black hole mass to be 2.97+0.51-0.77 × 109M; this is the first time that the dynamical mass of a black hole has been measured atz ∼ 6. This mass is consistent with that determined using the scaling relations from quasar emission lines. A massive stellar component (on the order of 109 M) may have already existed when the Universe was only ∼0.93 Gyr old. The relations between the black hole mass and the baryonic mass of this quasar indicate that the central supermassive black hole may have formed before its host galaxy. 
    more » « less
  4. The broad-line region (BLR) of active galactic nuclei (AGNs) traces gas close to the central supermassive black hole (BH). Recent reverberation mapping (RM) and interferometric spectro-astrometry data have enabled detailed investigations of the BLR structure and dynamics as well as estimates of the BH mass. These exciting developments have motivated comparative investigations of BLR structures using different broad emission lines. In this work, we have developed a method to simultaneously model multiple broad lines of the BLR from a single-epoch spectrum. We applied this method to the five strongest broad emission lines (Hα, Hβ, Hγ, Paβ, and He Iλ5876) in the UV-to-near-IR spectrum of NGC 3783, a nearby Type I AGN that has been well studied by RM and interferometric observations. Fixing the BH mass to the published value, we fit these line profiles simultaneously to constrain the BLR structure. We find that the differences between line profiles can be explained almost entirely as being due to different radial distributions of the line emission. We find that using multiple lines in this way also enables one to measure some important physical parameters, such as the inclination angle and virial factor of the BLR. The ratios of the derived BLR time lags are consistent with the expectation of theoretical model calculations and RM measurements. 
    more » « less
  5. Abstract We present new Atacama Large Millimeter/submillimeter Array observations that, for the first time, detect hydrogen and helium radio recombination lines from a protoplanetary disk. We imaged the Orion Nebula Cluster at 3.1 mm with a spectral setup that covered then= 42 → 41 transitions of hydrogen (H41α) and helium (He41α). The unprecedented sensitivity of these observations enables us to search for radio recombination lines toward the positions of ∼200 protoplanetary disks. We detect H41αfrom 17 disks, all of which are HST-identified “proplyds.” The detected H41αemission is spatially coincident with the locations of proplyd ionization fronts, indicating that proplyd H41αemission is produced by gas that has been photoevaporated off the disk and ionized by UV radiation from massive stars. We measure the fluxes and widths of the detected H41αlines and find line fluxes of ∼30–800 mJy km s−1and line widths of ∼30–90 km s−1. The derived line widths indicate that the broadening of proplyd H41αemission is dominated by outflowing gas motions associated with external photoevaporation. The derived line fluxes, when compared with measurements of 3.1 mm free–free flux, imply that the ionization fronts of H41α-detected proplyds have electron temperatures of ∼6000–11,000 K and electron densities of ∼106–107cm−3. Finally, we detect He41αtoward one H41α-detected source and find evidence that this system is helium-rich. Our study demonstrates that radio recombination lines are readily detectable in ionized photoevaporating disks, providing a new way to measure disk properties in clustered star-forming regions. 
    more » « less