Double-peaked Lines Research Articles

Dark progenitors and massive descendants: A first ALMA perspective of radio-selected near-IR-dark galaxies in the COSMOS field

We present the first spectroscopic ALMA follow-up for a pilot sample of nine radio-selected near-IR-dark galaxies in the COSMOS field. These sources were initially selected as radio-detected sources (S3 GHz > 12.65 μJy) without an optical or near-IR (NIR) counterpart in the COSMOS2015 catalog (Ks ≳ 24.7 mag), and just three of them were subsequently detected in the deeper COSMOS2020. Several studies highlighted that this selection could provide a population of highly dust-obscured, massive, and star-bursting galaxies. With these new ALMA observations, we assess the spectroscopic redshifts of this pilot sample of sources and improve the quality of the physical properties estimated through SED-fitting. Moreover, we measure the quantity of molecular gas inside these galaxies and forecast their potential evolutionary path, finding that the radio-selected NIR-dark galaxies might likely represent a population of high-z progenitors of the massive and passive galaxies that were discovered at z ∼ 3. Finally, we present some initial constraints on the kinematics of the interstellar medium within the analyzed galaxies, reporting a high fraction (∼55%) of double-peaked lines that can be interpreted as the signature of a rotating structure in our targets or as the presence of major mergers in our sample. The results we present here show the scientific potential of (sub)mm observations for this elusive population of galaxies and highlight the potential contribution of these sources to the evolution of the massive and passive galaxies at high z.

Anatomy of an ionized bubble: NIRCam grism spectroscopy of the z=6.6 double-peaked Lyman-alpha emitter COLA1 and its environment

The increasingly neutral intergalactic gas at $z>6$ impacts the Lyman-alpha (Lyalpha ) flux observed from galaxies. One luminous galaxy, COLA1, stands out because of its unique double-peaked Lyalpha line at $z=6.6$, unseen in any simulation of reionization. Here, we present JWST/NIRCam wide-field slitless spectroscopy in a 21 arcmin2 field centered on COLA1. We find 141 galaxies spectroscopically selected through the doublet at $5.35<z<6.95$, with 40 of these sources showing Hbeta . For COLA1, we additionally detect $ as well as Hgamma . We measure a systemic redshift of $z=6.5917$ for COLA1, confirming the classical double-peak nature of the Lyalpha profile. This implies that it resides in a highly ionized bubble and that it is leaking ionizing photons with a high escape fraction of $f_ esc (LyC) =20$--$50$<!PCT!>, making it a prime laboratory to study Lyman continuum escape in the Epoch of Reionization. COLA1 shows all the signs of a prolific ionizer with a Lyalpha escape fraction of $81 Balmer decrement indicating no dust, a steep UV slope ($ UV 0.4$), and a star-formation surface density $ that of typical galaxies at similar redshift. We detect five galaxies in COLA1's close environment ($ z<0.02$). Exploiting the high spectroscopic completeness inherent to grism surveys, and using mock simulations that fully mimic the selection function, we show that the number of detected companions is very typical for a normal similarly UV-bright UV galaxy -- that is, the ionized bubble around COLA1 is unlikely to be due to an excessively large over-density. Instead, the measured ionizing properties suggest that COLA1 by itself might be powering the bubble required to explain its double-peaked Lyalpha profile ion 0.7$ pMpc), with only minor contributions from detected neighbors ($-19.5 UV

Doppler Tomography of the Circumstellar Disk of the Be Star κ Draconis

κ Draconis is a binary system with a classical Be star as the primary component. Its emission-line spectrum consists of hydrogen lines, notably the Hα line with peak intensity ratio (V/R) variations phase-locked with the orbital period P = 61.55 days. Among binaries demonstrating the Be phenomenon, κ Dra stands out as one of a few systems with a discernible mass of its secondary component. Based on more than 200 spectra obtained in 2014–2023, we verified the physical parameters and constructed the mass function. We used part of these data obtained in 2014–2021 to investigate regions in the circumstellar disk of the primary component that emit the Hα line using the Doppler tomography method. The results show that the disk has a non-uniform density distribution with a prominent enhancement at Vy ≈ 99 km s−1 and Vx≈−6 km s−1 that corresponds to a cloud-like source of the double-peaked Hα line profile. We argue that this enhancement’s motion is responsible for the periodic variations in the Hα V/R ratio, which is synchronised in orbital phase with the radial velocity (RV) of absorption lines from the atmosphere of the primary component.

Panic at the ISCO: Time-varying Double-peaked Broad Lines from Evolving Accretion Disks Are Common among Optically Variable AGNs

About 3%–10% of Type I active galactic nuclei (AGNs) have double-peaked broad Balmer lines in their optical spectra originating from the motion of gas in their accretion disk. Double-peaked profiles arise not only in AGNs, but occasionally appear during optical flares from tidal disruption events and changing-state AGNs. In this paper, we identify 250 double-peaked emitters (DPEs) among a parent sample of optically variable broad-line AGNs in the Zwicky Transient Facility (ZTF) survey, corresponding to a DPE fraction of 19%. We model spectra of the broad Hα emission-line regions and provide a catalog of the fitted accretion disk properties for the 250 DPEs. Analysis of power spectra derived from the 5 yr ZTF light curves finds that DPE light curves have similar amplitudes and power-law indices to other broad-line AGNs. Follow-up spectroscopy of 12 DPEs reveals that ∼50% display significant changes in the relative strengths of their red and blue peaks over long 10–20 yr timescales, indicating that broad-line profile changes arising from spiral arm or hotspot rotation are common among optically variable DPEs. Analysis of the accretion disk parameters derived from spectroscopic modeling provides evidence that DPEs are not in a special accretion state, but are simply normal broad-line AGNs viewed under the right conditions for the accretion disk to be easily visible. We include inspiraling supermassive black hole binary candidate SDSSJ1430+2303 in our analysis, and discuss how its photometric and spectroscopic variability is consistent with the disk-emitting AGN population in the ZTF survey.

A candidate for central tidal disruption event in the broad line AGN SDSS J1605 with double-peaked broad H β

ABSTRACT In this letter, motivated by double-peaked broad Balmer emission lines probably related to tidal disruption events (TDEs), a potential TDE candidate is reported in Sloan Digital Sky Survey (SDSS) J160536+134838 (=SDSS J1605) at z ∼ 0.44 having quasar-like spectrum but with double-peaked broad H β. The long-term Catalina Sky Survey (CSS) light curve can be naturally described by a main-sequence star of $2.82_{-0.19}^{+0.20}{\rm M_\odot }$ disrupted by the central black hole (BH) of $144_{-21}^{+26} \times 10^6{\rm M_\odot }$ in SDSS J1605. Meanwhile, the All-Sky Automated Survey for SuperNovae light curves afterwards show none apparent trend variability, indicating the bright CSS flare in SDSS J1605 unique and different enough from variability of normal active galactic nuclei (AGN). Furthermore, there is a consistency between the TDE model determined sizes of debris with the sizes of emission regions for the double-peaked broad H β described by the accretion disc model, supporting the disc-like broad emission line regions (BLRs) probably related to a central TDE in SDSS J1605. And the virial BH mass ∼7 times higher than the TDE model determined value can be naturally explained by R–L relation determined BLRs sizes very larger than the actual distance of emission regions related to TDEs debris in SDSS J1605. Although no clear conclusion on double-peaked broad lines absolutely related to TDEs, the results here provide clues to detect potential TDEs in AGN with double-peaked broad lines.

A Rotating Accretion Disk around MWC 297, a Young B1.5Ve Star

High-resolution spectra with iSHELL on IRTF in the K and M bands of the young, heavily accreting B1.5Ve star MWC 297 show numerous double-peaked CO lines. These CO lines originate in an inclined gaseous disk in Keplerian rotation. MWC 297 is the only early B star known to show a Keplerian disk in CO. Analysis of the spectra shows that 12CO 1–0 is optically thick for the low excitation lines. Even the 13CO 1–0 and 12CO 2–1 have somewhat optically thick lines at low J levels. We find that the CO emission in the disk can be fitted with CO being in a narrow ring at a radius of 12 au, with a temperature of 1500 K and a CO column density of 1.6 × 1018 cm−2. This model underestimates the line strength of high-J lines, indicating that they are excited by fluorescence. The CO overtone lines have a similar temperature. The 13CO lines are much brighter than expected from interstellar isotope ratios. The 13CO lines are wider than the 12CO ones, suggesting different excitation conditions. The same is true for 12CO 2–1. We see strong absorption in 12CO and 13CO 1–0 at low J levels, which is due to two cold foreground clouds. These clouds, one with a temperature of 8.3 K and a column density of 6.7 1017 cm−2 and the other one colder and with lower column density, can fully account for the observed extinction toward MWC 297.

Searching for Phase-Locked Variations of the Emission-Line Profiles in Binary Be Stars

There is growing evidence that many Be stars are parts of binary systems. As the B-type primaries are very fast rotators and their spectral lines may be distorted by the circumstellar material, it is not easy to measure their radial velocity directly from the spectral lines. It has been shown that some Be binaries exhibit peak intensity variations consisting of double-peaked Hα lines that are phase-locked with orbital periods. We searched for such variations in the spectra of 12 Be stars, including several known and suspected binaries. Our results include confirmation of the orbital periods in ν Geminorum, ϵ Capricorni, κ Draconis, 60 Cygni, and V2119 Cygni, its refinement in o Puppis, as well as suggesting hints for binarity in o Aquarii, BK Camelopardalis, and 10 Cassiopeae. Monitoring of the Hα line profile variations in β Canis Minoris for over the last 10 years gives further support to the existence of a 182.5-day period found earlier in a smaller set of data. A similar but still preliminary period (179.6 days) was found in the Hα line profile variations in ψ Persei. It is shown for the first time that ν Geminorum exhibits phase-locked variations in the Hα emission peak intensity ratio and, therefore, is a part of the inner binary in this triple system. Our results show that the mentioned phase-locked peak intensity variations are observed in more Be binary systems than previously known and can be used to search for binarity of Be stars when application of other methods is inconclusive.

Massive pre-main-sequence stars in M17

Context.The young massive-star-forming region M17 contains optically visible massive pre-main-sequence stars that are surrounded by circumstellar disks. Such disks are expected to disappear when these stars enter the main sequence. The physical and dynamical structure of these remnant disks are poorly constrained, especially the inner regions where accretion, photo-evaporation, and companion formation and migration may be ongoing.Aims.We aim to constrain the physical properties of the inner parts of the circumstellar disks of massive young stellar objects B243 (6M⊙) and B331 (12M⊙), two systems for which the central star has been detected and characterized previously despite strong dust extinction.Methods.Two-dimensional radiation thermo-chemical modelling with PRODIMOof double-peaked hydrogen lines of the Paschen and Brackett series observed with X-shooter was used to probe the properties of the inner disk of the target sources. The model was modified to treat these lines. Additionally, the dust structure was studied by fitting the optical and near-infrared spectral energy distribution.Results.B243 features a hot gaseous inner disk with dust at the sublimation radius at ~3 AU. The disk appears truncated at roughly 6.5 AU; a cool outer disk of gas and dust may be present, but it cannot be detected with our data. B331 also has a hot gaseous inner disk. A gap separates the inner disk from a colder dusty outer disk starting at up to ~100 AU. In both sources the inner disk extends to almost the stellar surface. Chemistry is essential for the ionization of hydrogen in these disks.Conclusions.The lack of a gap between the central objects and these disks suggests that they accrete through boundary-layer accretion. This would exclude the stars having a strong magnetic field. Their structures suggest that both disks are transitional in nature, that is to say they are in the process of being cleared, either through boundary-layer accretion, photo-evaporation, or through companion activity.

Evidence to disfavour dual-core systems leading to double-peaked narrow emission lines

ABSTRACTIn this paper, an interesting method is proposed to test dual-core systems for double-peaked narrow emission lines through a precious dual-core system with double-peaked narrow Balmer lines in the main galaxy but with single-peaked narrow Balmer lines in the companion galaxy. Under a dual-core system, considering narrow Balmer (Hα and Hβ) emissions (fe,α and fe,β) from a companion galaxy that are covered by the SDSS fiber for the main galaxy and narrow Balmer emissions (fc,α and fc,β) from the companion galaxy covered by the SDSS fiber for the companion galaxy, the same flux ratios fe,α/fc,α = fe,β/fc,β can be expected, owing to the totally similar physical conditions of each narrow Balmer emission region. Next, the precious dual-core system in SDSS J2219–0938 is discussed. After subtracting the pPXF code determined stellar lights, double-peaked narrow Balmer emission lines are confirmed in the main galaxy with a confidence level higher than 5σ, but single-peaked narrow Balmer emission lines are confirmed in the companion galaxy. Through measured fluxes of the emission components, fe,α/fc,α is about 0.82, which is different from fe,β/fc,β ∼ 0.52, which disfavours a dual-core system for the double-peaked narrow Balmer emission lines in SDSS J2219–0938.

Cyanopolyyne Chemistry in the L1544 Prestellar Core: New Insights from GBT Observations

We report a comprehensive study of the cyanopolyyne chemistry in the prototypical prestellar core L1544. Using the 100 m Robert C. Byrd Green Bank Telescope, we observe three emission lines of HC3N, nine lines of HC5N, five lines of HC7N, and nine lines of HC9N. HC9N is detected for the first time toward the source. The high spectral resolution (∼0.05 km s−1) reveals double-peak spectral line profiles with the redshifted peak a factor 3–5 brighter. Resolved maps of the core in other molecular tracers indicate that the southern region is redshifted. Therefore, the bulk of the cyanopolyyne emission is likely associated with the southern region of the core, where free carbon atoms are available to form long chains, thanks to the more efficient illumination of the interstellar field radiation. We perform a simultaneous modeling of the HC5N, HC7N, and HC9N lines to investigate the origin of the emission. To enable this analysis, we performed new calculation of the collisional coefficients. The simultaneous fitting indicates a gas kinetic temperature of 5–12 K, a source size of 80″, and a gas density larger than 100 cm−3. The HC5N:HC7N:HC9N abundance ratios measured in L1544 are about 1:6:4. We compare our observations with those toward the well-studied starless core TMC-1 and with the available measurements in different star-forming regions. The comparison suggests that a complex carbon chain chemistry is active in other sources and is related to the presence of free gaseous carbon. Finally, we discuss the possible formation and destruction routes in light of the new observations.

Chemistry and dynamics of the prestellar core L1544

We aim to quantify the effect of chemistry on the infall velocity in the prestellar core L1544. Previous observational studies have found evidence for double-peaked line profiles for the rotational transitions of several molecules, which cannot be accounted for with the models presently available for the physical structure of the source, without ad hoc up-scaling of the infall velocity. We ran one-dimensional hydrodynamical simulations of the collapse of a core with L1544-like properties (in terms of mass and outer radius), using a state-of-the-art chemical model with a very large chemical network combined with an extensive description of molecular line cooling, determined via radiative transfer simulations, with the aim of determining whether these expansions of the simulation setup (as compared to previous models) can lead to a higher infall velocity. After running a series of simulations where the simulation was sequentially simplified, we found that the infall velocity is almost independent of the size of the chemical network or the approach to line cooling. We conclude that chemical evolution does not have a large impact on the infall velocity, and that the higher infall velocities that are implied by observations may be the result of the core being more dynamically evolved than what is now thought, or alternatively the average density in the simulated core is too low. However, chemistry does have a large influence on the lifetime of the core, which varies by about a factor of two across the simulations and grows longer when the chemical network is simplified. Therefore, although the model is subject to several sources of uncertainties, the present results clearly indicate that the use of a small chemical network leads to an incorrect estimate of the core lifetime, which is naturally a critical parameter for the development of chemical complexity in the precollapse phase.

A chemical study of nine star-forming regions with evidence of infall motion

ABSTRACT The study of the physical and chemical properties of gas infall motion in the molecular clumps helps us understand the initial stages of star formation. We used the FTS wide-sideband mode of the IRAM 30-m telescope to observe nine infall sources with significant double-peaked blue line profile. The observation frequency ranges are 83.7–91.5 GHz and 99.4–107.2 GHz. We have obtained numbers of molecular line data. Using xclass, a total of 7–27 different molecules and isotopic transition lines have been identified in these nine sources, including carbon chain molecules, such as CCH, c-C3H2 and HC3N. According to the radiation transfer model, we estimated the rotation temperatures and column densities of these sources. Chemical simulations adopting a physical model of HMSFRs are used to fit the observed molecular abundances. The comparison shows that most sources are in the early HMPO stage, with the inner temperature around several 10 K.

Modelling the flare in NGC 1097 from 1991 to 2004 as a tidal disruption event

ABSTRACT In the letter, interesting evidence is reported to support a central tidal disruption event (TDE) in the known AGN NGC 1097. Considering the motivations of TDE as one probable origination of emission materials of double-peaked broad emission lines and also as one probable explanation to changing-look AGN, it is interesting to check whether there are clues to support a TDE in NGC 1097, not only a changing-look AGN but also an AGN with double-peaked broad emission lines. Under the assumption that the onset of broad Hα emission was due to a TDE, the 13 yr-long (1991–2004) variability of double-peaked broad Hα line flux in NGC 1097 can be well predicted by theoretical TDE model, with a (1 − 1.5) M⊙ main-sequence star tidally disrupted by the central BH with TDE model determined mass about (5 − 8) × 107 M⊙. The results provide interesting evidence to not only support TDE-related origin of double-peaked broad line emission materials but also support TDE as an accepted physical explanation to physical properties of changing-look AGN.

Dependence of Narrow-line Region Sizes on [O iii] Luminosity in Low-redshift Active Galactic Nuclei with Double-peaked Broad Balmer Emission Lines

In this paper, simple but interesting results are reported on the upper limits of narrow-line region (NLR) sizes of a small sample of 38 low-redshift (z < 0.1) active galactic nuclei (AGN) with double-peaked broad emission lines (double-peaked BLAGN), in order to check whether the NLR sizes in type-1 AGN (broad line) and type-2 AGN (narrow line) obey a similar empirical dependence on [O iii] luminosity. In order to correct the inclination effects on projected NLR sizes of type-1 AGN, the accretion disk origin is commonly applied to describe the double-peaked broad Hα line, leading to the determined inclination angles of central disk-like broad-line regions of 38 double-peaked BLAGN. Then, considering the fixed Sloan Digital Sky Survey (SDSS) fiber radius, the upper limits of the NLR sizes of the 38 double-peaked BLAGN can be estimated. Meanwhile, a strong linear correlation between continuum luminosity and [O iii] luminosity is applied to confirm that the [O iii] emissions of the 38 double-peaked BLAGN are totally covered in the SDSS fibers. Considering the reddening-corrected measured [O iii] luminosity, the upper limits of the NLR sizes of the 38 double-peaked BLAGN are within a 99.9999% confidence interval of the expected results from the empirical relation between NLR size and [O iii] luminosity in type-2 AGN. In the current understanding, there are no challenges to the unified model of AGN through the space properties of NLRs.

ASAS J071404+7004.3 – a close, bright nova-like cataclysmic variable with gusty winds

ABSTRACT Despite being bright (V ≃ 11.8) and nearby (d = 212 pc) ASAS J071404+7004.3 has only recently been identified as a nova-like cataclysmic variable. We present time-resolved optical spectroscopy obtained at the Isaac Newton and the Hiltner and McGraw-Hill Telescopes, together with Swift X-ray and ultraviolet observations. We combined these with TESS photometry and find a period of 3.28 h and a mass transfer rate of $4\!-\!9\times 10^{-9}\, {\mathrm{M_{\odot }\, yr}^{-1}}$. Historical photometry shows at least one low state establishing the system as a VY Scl star. Our high-cadence spectroscopy also revealed rapidly changing winds emanating from the accretion disc. We have modelled these using the Monte Carlo python code and shown that all the emission lines could emanate from the wind – which would explain the lack of double-peaked lines in such systems. In passing, we discuss the effect of variability on the position of cataclysmic variables in the Gaia Hertzsprung–Russell diagram.

Mapping the Morphology and Kinematics of a Lyα-selected Nebula at z = 3.15 with MUSE

Abstract Recent wide-field integral-field spectroscopy has revealed the detailed properties of high-redshift Lyα nebulae, most often targeted due to the presence of an active galactic nucleus (AGN). Here, we use VLT/MUSE to resolve the morphology and kinematics of a nebula initially identified due to strong Lyα emission at z ∼ 3.2 (LABn06). Our observations reveal a two-lobed Lyα nebula, at least ∼173 pkpc in diameter, with a light-weighted centroid near a mid-infrared source (within ≈17.2 pkpc) that appears to host an obscured AGN. The Lyα emission near the AGN is also coincident in velocity with the kinematic center of the nebula, suggesting that the nebula is both morphologically and kinematically centered on the AGN. Compared to AGN-selected Lyα nebulae, the surface-brightness profile of this nebula follows a typical exponential profile at large radii (&gt;25 pkpc), although at small radii, the profile shows an unusual dip at the location of the AGN. The kinematics and asymmetry are similar to, and the C iv and He ii upper limits are consistent with, other AGN-powered Lyα nebulae. Double-peaked and asymmetric line profiles suggest that Lyα resonant scattering may be important in this nebula. These results support the picture of the AGN being responsible for powering a Lyα nebula that is oriented roughly in the plane of the sky. Further observations will explore whether the central surface-brightness depression is indicative of either an unusual gas or dust distribution or variation in the ionizing output of the AGN over time.

In Search of Infall Motion in Molecular Clumps. III. HCO+ (1-0) and H13CO+ (1-0) Mapping Observations toward Confirmed Infall Sources

Abstract The study of infall motion helps us to understand the initial stages of star formation. In this paper, we use the IRAM 30 m telescope to make mapping observations of 24 infall sources confirmed in previous work. The lines we use to track gas infall motions are HCO+ (1-0) and H13CO+ (1-0). All 24 sources show HCO+ emissions, while 18 sources show H13CO+ emissions. The HCO+ integrated intensity maps of 17 sources show clear clumpy structures; for the H13CO+ line, 15 sources show clumpy structures. We estimated the column density of HCO+ and H13CO+ using the RADEX radiation transfer code, and the obtained [HCO+]/[H2] and [H13CO+]/[HCO+] of these sources are about 10−11–10−7 and 10−3–1, respectively. Based on the asymmetry of the line profile of the HCO+, we distinguish these sources: 19 sources show blue asymmetric profiles, and the other sources show red profiles or symmetric peak profiles. For eight sources that have double-peaked blue line profiles and signal-to-noise ratios greater than 10, the RATRAN model is used to fit their HCO+ (1-0) lines, and to estimate their infall parameters. The mean V in of these sources is 0.3–1.3 km s−1, and the M ̇ in is about 10−3–10−4 M ⊙ yr−1, which is consistent with the results of intermediate or massive star formation in previous studies. The V in estimated from the Myers model is 0.1–1.6 km s−1, and the M ̇ in is within 10−3–10−5 M ⊙ yr−1. In addition, some identified infall sources show other star-forming activities, such as outflows and maser emissions. Especially for those sources with a double-peaked blue asymmetric profile, most of them have both infall and outflow evidence.

Solid-Density Ion Temperature from Redshifted and Double-Peaked Stark Line Shapes.

Heβ spectral line shapes are important for diagnosing temperature and density in many dense plasmas. This work presents Heβ line shapes measured with high spectral resolution from solid-density plasmas with minimized gradients. The line shapes show hallmark features of Stark broadening, including quantifiable redshifts and double-peaked structure with a significant dip between the peaks; these features are compared to models through a Markov chain MonteCarlo framework. Line shape theory using the dipole approximation can fit the width and peak separation of measured line shapes, but it cannot resolve an ambiguity between electron density n_{e} and ion temperature T_{i}, since both parameters influence the strength of quasistatic ion microfields. Here a line shape model employing a full Coulomb interaction for the electron broadening computes self-consistent line widths and redshifts through the monopole term; redshifts have different dependence on plasma parameters and thus resolve the n_{e}-T_{i} ambiguity. The measured line shapes indicate densities that are 80-100% of solid, identifying a regime of highly ionized but well-tamped plasma. This analysis also provides the first strong evidence that dense ions and electrons are not in thermal equilibrium, despite equilibration times much shorter than the duration of x-ray emission; cooler ions may arise from nonclassical thermalization rates or anomalous energy transport. The experimental platform and diagnostic technique constitute a promising new approach for studying ion-electron equilibration in dense plasmas.

Structured velocity field in the inner envelope of B335: ALMA observations of rare CO isotopologues

Context. Studies of Class 0 objects allow to characterize the dynamical processes taking place at the onset of the star formation process and to determine the physical mechanisms responsible for the outcome of the collapse. Observations of dense gas tracers allow for the characterization of key kinematics of the gas that are directly involved in the star formation process, such as infall, outflow, and rotation. Aims. This work is aimed at investigating the molecular line velocity profiles of the Class 0 protostellar object B335 and attempts to place constraints on the infall motions happening in the circumstellar gas of the object. Methods. We present observations of C17O (1–0), C18O (1–0), and 12CO (2–1) transitions along with an analysis of spectral profiles at envelope radii between 100 and 860 au. Results. C17O emission presents a double-peaked line profile distributed in a complex velocity field. Both peaks present an offset of 0.2–1 km s−1 from the systemic velocity of the source in the probed area. The optical depth of the C17O emission has been estimated and found to be less than 1, suggesting that the two velocity peaks trace two distinct velocity components of the gas in the inner envelope. Conclusions. After discarding possible motions that could produce the complex velocity pattern, such as rotation and outflow, we conclude that infall motions are responsible for producing the velocity field. Because inside-out symmetric collapse cannot explain those observed profiles, it is suggested that these are produced by non-isotropic accretion from the envelope into the central source along the outflow cavity walls.

AGNs on the Move: A Search for Off-nuclear AGNs from Recoiling Supermassive Black Holes and Ongoing Galaxy Mergers with the Zwicky Transient Facility

Abstract A supermassive black hole (SMBH) ejected from the potential well of its host galaxy via gravitational wave recoil carries important information about the mass ratio and spin alignment of the pre-merger SMBH binary. Such a recoiling SMBH may be detectable as an active galactic nucleus (AGN) broad-line region offset by up to 10 kpc from a disturbed host galaxy. We describe a novel methodology using forward modeling with The Tractor to search for such offset AGNs in a sample of 5493 optically variable AGNs detected with the Zwicky Transient Facility (ZTF). We present the discovery of nine AGNs that may be spatially offset from their host galaxies and are candidates for recoiling SMBHs. Five of these offset AGNs exhibit double-peaked broad Balmer lines, which may have arisen from unobscured accretion disk emission, and four show radio emission indicative of a relativistic jet. The fraction of double-peaked emitters in our spatially offset AGN sample is significantly larger than the 16% double-peaked emitter fraction observed for ZTF AGNs overall. In our sample of variable AGNs we also identified 52 merging galaxies, including a new spectroscopically confirmed dual AGN. Finally, we detected the dramatic rebrightening of SDSS 1133, a previously discovered variable object and recoiling SMBH candidate, in ZTF. The flare was accompanied by the reemergence of strong P Cygni line features, indicating that SDSS 1133 may be an outbursting luminous blue variable star.