Showing 1–50 of 149 results for author: Zakamska, N

Cataclysmic variables from Sloan Digital Sky Survey -- V (2020-2023) identified using machine learning

Authors: Keith Inight, Boris T. Gänsicke, Axel Schwope, Scott F. Anderson, Elmé Breedt, Joel R. Brownstein, Sebastian Demasi, Susanne Friedrich, J. J. Hermes, Knox S. Long, Timothy Mulvany, Gautham A. Pallathadka, Mara Salvato, Simone Scaringi, Matthias R. Schreiber, Guy S. Stringfellow, John R. Thorstensen, Nadia L. Zakamska

Abstract: SDSS-V is carrying out a dedicated survey for white dwarfs, single and in binaries, and we report the analysis of the spectroscopy of 505 cataclysmic variables (CVs) and CV candidates obtained during the first 34 months of observations of SDSS-V. We developed a convolutional neural network (CNN) to aid with the identification of CV candidates among the over 2 million SDSS-V spectra obtained with t… ▽ More SDSS-V is carrying out a dedicated survey for white dwarfs, single and in binaries, and we report the analysis of the spectroscopy of 505 cataclysmic variables (CVs) and CV candidates obtained during the first 34 months of observations of SDSS-V. We developed a convolutional neural network (CNN) to aid with the identification of CV candidates among the over 2 million SDSS-V spectra obtained with the BOSS spectrograph. The CNN reduced the number of spectra that required visual inspection to $\simeq2$ per cent of the total. We identified 779 CV spectra among the CNN-selected candidates, plus an additional 37 CV spectra that the CNN misclassified, but that were found serendipitously by human inspection of the data. Analysing the SDSS-V spectroscopy and ancillary data of the 505 CVs in our sample, we report 62 new CVs, spectroscopically confirm 243 and refute 13 published CV candidates, and we report 68 new or improved orbital periods. We discuss the completeness and possible selection biases of the machine learning methodology, as well as the effectiveness of targeting CV candidates within SDSS-V. Finally, we re-assess the space density of CVs, and find $1.3\times 10^{-5}\,\mathrm{pc^{-3}}$. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Comments: Submitted to MNRAS. 23 pages, 25 figures. Supplementary data contains details of 505 CVs including spectra and light curves

JWST ERS Program Q3D: The pitfalls of virial BH mass constraints shown in a z = 3 quasar with an ultramassive host

Authors: Caroline Bertemes, Dominika Wylezalek, David S. N. Rupke, Nadia L. Zakamska, Sylvain Veilleux, Benjamin Beckmann, Andrey Vayner, Swetha Sankar, Yuzo Ishikawa, Nadiia Diachenko, Weizhe Liu, Yu-Ching Chen, Jerome Seebeck, Dieter Lutz, Guilin Liu

Abstract: We present JWST MIRI/NIRSpec observations of the extremely red quasar SDSS J165202.64+172852.3 at z~3, one of the most luminous quasars known to date, driving powerful outflows and hosting a clumpy starburst, amidst several interacting companions. We estimate the black hole (BH) mass of the system based on the broad H$α$ and H$β$ lines, as well as the Pa$β$ emission in the IR and MgII in the UV. W… ▽ More We present JWST MIRI/NIRSpec observations of the extremely red quasar SDSS J165202.64+172852.3 at z~3, one of the most luminous quasars known to date, driving powerful outflows and hosting a clumpy starburst, amidst several interacting companions. We estimate the black hole (BH) mass of the system based on the broad H$α$ and H$β$ lines, as well as the Pa$β$ emission in the IR and MgII in the UV. We recover a very broad range of mass estimates, with constraints ranging between log $M_{\rm BH}$=9 and 10.1, which is exacerbated if imposing a uniform BLR geometry at all wavelengths. Several factors may contribute to the large spread: measurement uncertainties (insufficient sensitivity to detect the broadest component of the faint Pa$β$ line, spectral blending, ambiguities in the broad/narrow component distinction), lack of virial equilibrium (in a system characterised by powerful outflows and rapid accretion), and uncertainties on the luminosity-inferred size of the broad line region, a.o. given central dust obscuration. We constrain the stellar mass via SED fitting, suggesting the host to be extremely massive at $10^{12.8\pm 0.5} M_\odot$ - ~2 dex above the characteristic mass of the Schechter fit to the z=3 stellar mass function. Notably, J1652's central BH might be interpreted as being either undermassive, overmassive, or in line with the BH mass-stellar mass relation, depending on the choice of assumptions. The recovered Eddington ratio varies accordingly, but exceeds 10% in any case. We put our results into context by providing an extensive overview and discussion of recent literature results and their associated assumptions. Our findings provide an important demonstration of the uncertainties inherent in virial BH mass estimates, which are of particular relevance in the JWST era given the growing number of studies on rapidly accreting quasars at high redshift. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Comments: 15 pages, 7 figures, submitted to A&A, comments welcome

VODKA-JWST: Synchronized growth of two SMBHs in a massive gas disk? A 3.8 kpc separation dual quasar at cosmic noon with JWST NIRSpec IFU

Authors: Yuzo Ishikawa, Nadia L. Zakamska, Yue Shen, Xin Liu, Yu-Ching Chen, Hsiang-Chih Hwang, Andrey Vayner, Sylvain Veilleux, David S. N. Rupke, Dominika Wylezalek, Arran C. Gross, Swetha Sankar, Nadiia Diachenko

Abstract: The search for dual supermassive black holes (SMBHs) is of immense interest in modern astrophysics. Galaxy mergers may be an important route to fuel and to produce SMBH pairs. Actively accreting SMBH pairs can be observed as a dual quasar, which are vital probes of SMBH growth. Gaia observations have enabled a novel technique to systematically search for such dual quasars at previously unreachable… ▽ More The search for dual supermassive black holes (SMBHs) is of immense interest in modern astrophysics. Galaxy mergers may be an important route to fuel and to produce SMBH pairs. Actively accreting SMBH pairs can be observed as a dual quasar, which are vital probes of SMBH growth. Gaia observations have enabled a novel technique to systematically search for such dual quasars at previously unreachable sub-kpc scales, based on the small jitters of the light centroid as the two quasars vary stochastically. Here we present the first detailed study of a 0.46'', 3.8 kpc separation, VODKA-selected dual quasar, J0749+2255, at $z=2.17$ using JWST/NIRSpec integral field unit spectroscopy. This is one of the most distant, small separation dual quasars identified today. Dual quasars at cosmic noon are not well characterized. We detect the faint ionized gas of the host galaxy, best traced by the narrow \ha\ emission. Line ratio diagnostics show a mix of ionization from the two quasars and intense star formation. The spatially-resolved spectra of the two quasars suggest that they have very similar black hole properties (two $M_{BH}\sim 10^9\ \textrm{M}_{\odot}$ with large Eddington ratio reaching $L/L_{Edd}\sim0.2$) hinting at the possible synchronized growth and accretion from the same gas supply. Surprisingly, the ionized gas kinematics suggest an extended, rotating disk rather than a disturbed system that would be expected in a major gas-rich galaxy merger. While it is unclear if J0749+2255 is representative of the dual quasar evolution, the observations with JWST revealed a major puzzle. It would be interesting to see what observations of other dual quasars will show. △ Less

Submitted 12 March, 2024; originally announced March 2024.

Comments: 21 pages, 9 figures, submitted

Journal ref: ApJ 2024

VODKA-JWST: A 3.8 kpc dual quasar at cosmic noon in a powerful starburst galaxy with JWST/MIRI IFU

Authors: Yu-Ching Chen, Yuzo Ishikawa, Nadia L. Zakamska, Xin Liu, Yue Shen, Hsiang-Chih Hwang, David Rupke, Andrey Vayner, Arran C. Gross, Weizhe Liu, Dominika Wylezalek, Sylvain Veilleux, Caroline Bertemes, Nadiia Diachenko, Swetha Sankar

Abstract: Dual quasars, two active supermassive black holes at galactic scales, represent crucial objects for studying the impact of galaxy mergers and quasar activity on the star formation rate (SFR) within their host galaxies, particularly at cosmic noon when SFR peaks. We present JWST/MIRI mid-infrared integral field spectroscopy of J074922.96+225511.7, a dual quasar with a projected separation of 3.8 ki… ▽ More Dual quasars, two active supermassive black holes at galactic scales, represent crucial objects for studying the impact of galaxy mergers and quasar activity on the star formation rate (SFR) within their host galaxies, particularly at cosmic noon when SFR peaks. We present JWST/MIRI mid-infrared integral field spectroscopy of J074922.96+225511.7, a dual quasar with a projected separation of 3.8 kilo-parsec at a redshift $z$ of 2.17. We detect spatially extended [Fe II] 5.34$\rm μ$m and polycyclic aromatic hydrocarbon (PAH) 3.3$μ$m emissions from the star formation activity in its host galaxy. We derive the SFR of 10$^{3.0\pm0.2}$ M$_{\odot}$ yr$^{-1}$ using PAH 3.3$μ$m, which is five times higher than that derived from the cutoff luminosity of the infrared luminosity function for galaxies at $z\sim2$. While the SFR of J0749+2255 agrees with that of star-forming galaxies of comparable stellar mass at the same redshifts, its molecular gas content falls short of expectations based on the molecular Kennicutt-Schmidt law. This discrepancy may result from molecular gas depletion due to the longer elevated stage of star formation, even after the molecular gas reservoir is depleted. We do not observe any quasar-driven outflow that impacts PAH and [Fe II] in the host galaxy based on the spatially resolved maps. From the expected flux in PAH-based star formation, the [Fe II] line likely originates from the star-forming regions in the host galaxy. Our study highlights the stardust nature of J0749+2255, indicating a potential connection between the dual quasar phase and intense star formation activities. △ Less

Submitted 9 March, 2024; v1 submitted 6 March, 2024; originally announced March 2024.

Comments: 12 pages, 7 figures, 1 table, submitted to ApJ, comments are welcome

Evidence for Intrinsic X-ray Weakness Among Red Quasars at Cosmic Noon

Authors: Yilun Ma, Andy Goulding, Jenny E. Greene, Nadia L. Zakamska, Dominika Wylezalek, Yan-Fei Jiang

Abstract: Quasar feedback is a key ingredient in shaping galaxy evolution. A rare population of extremely red quasars (ERQs) at $z=2-3$ are often associated with high-velocity [OIII]$\lambda5008$ outflows and may represent sites of strong feedback. In this paper, we present an X-ray study of 50 ERQs to investigate the link between the X-ray and outflow properties of these intriguing objects. Using hardness… ▽ More Quasar feedback is a key ingredient in shaping galaxy evolution. A rare population of extremely red quasars (ERQs) at $z=2-3$ are often associated with high-velocity [OIII]$\lambda5008$ outflows and may represent sites of strong feedback. In this paper, we present an X-ray study of 50 ERQs to investigate the link between the X-ray and outflow properties of these intriguing objects. Using hardness ratio analysis, we confirm that the ERQs are heavily obscured systems with gas column density reaching $N_\mathrm{H}=10^{23-24}\,\mathrm{cm^{-2}}$. We identify 17 X-ray-non-detected ERQs at high mid-infrared luminosities of $νL_\mathrm{ν,6μm}\gtrsim3\times10^{46}\,\mathrm{erg\,s^{-1}}$. By stacking the X-ray observations, we find that the non-detected ERQs are on average underluminous in X-rays by a factor of $\sim5$ for their bolometric luminosities. We consider such X-ray weakness to be due to both gas absorption and intrinsic factors. Moreover, we find that the X-ray-weak sources also display higher-velocity outflows. One option to explain this trend is that weaker X-rays facilitate more vigorous line-driven winds, which then accelerate the [OIII]-emitting gas to kpc-scales. Alternatively, super-Eddington accretion could also lead to intrinsic X-ray weakness and more powerful continuum-driven outflow. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: 22 pages, 10 figures, submitted to ApJ

JWST discovers an AGN ionization cone but only weak radiative-driven feedback in a powerful $z$$\approx$3.5 radio-loud AGN

Authors: Wuji Wang, Dominika Wylezalek, Carlos De Breuck, Joël Vernet, David S. N. Rupke, Nadia L. Zakamska, Andrey Vayner, Matthew D. Lehnert, Nicole P. H. Nesvadba, Daniel Stern

Abstract: We present the first results from a JWST program studying the role played by powerful radio jets in the evolution of the most massive galaxies at the onset of Cosmic Noon. Using NIRSpec integral field spectroscopy, we detect 24 rest-frame optical emission lines from the $z=3.5892$ radio galaxy 4C+19.71. 4C+19.71 contains one of the most energetic radio jets known, making it perfect for testing rad… ▽ More We present the first results from a JWST program studying the role played by powerful radio jets in the evolution of the most massive galaxies at the onset of Cosmic Noon. Using NIRSpec integral field spectroscopy, we detect 24 rest-frame optical emission lines from the $z=3.5892$ radio galaxy 4C+19.71. 4C+19.71 contains one of the most energetic radio jets known, making it perfect for testing radio-mode feedback on the interstellar medium (ISM) of a $M_{\star}\sim10^{11}\,\rm M_{\odot}$ galaxy. The rich spectrum enables line ratio diagnostics showing that the radiation from the active galactic nucleus (AGN) dominates the ionization of the entire ISM out to at least $25\,$kpc, the edge of the detection. Sub-kpc resolution reveals filamentary structures and emission blobs in the warm ionized ISM distributed on scales of $\sim5$ to $\sim20\,$kpc. A large fraction of the extended gaseous nebula is located near the systemic velocity. This nebula may thus be the patchy ISM which is illuminated by the AGN after the passage of the jet. A radiatively-driven outflow is observed within $\sim5\,$kpc from the nucleus. The inefficient coupling ($\lesssim 10^{-4}$) between this outflow and the quasar and the lack of extreme gas motions on galactic scales are inconsistent with other high-$z$ powerful quasars. Combining our data with ground-based studies, we conclude that only a minor fraction of the feedback processes is happening on $<25\,$kpc scales. △ Less

Submitted 15 January, 2024; v1 submitted 4 January, 2024; originally announced January 2024.

Comments: Accept for publication in Astronomy & Astrophysics, 10 pages and 5 figures in main text

The Mira Distance to M101 and a 4% Measurement of H0

Authors: Caroline D. Huang, Wenlong Yuan, Adam G. Riess, Warren Hack, Patricia A. Whitelock, Nadia L. Zakamska, Stefano Casertano, Lucas M. Macri, Massimo Marengo, John W. Menzies, Randall K. Smith

Abstract: The giant spiral galaxy M101 is host to the nearest recent Type Ia Supernova (SN 2011fe) and thus has been extensively monitored in the near-infrared to study the late-time lightcurve of the supernova. Leveraging this existing baseline of observations, we derive the first Mira-based distance to M101 by discovering and classifying a sample of 211 Miras with periods ranging from 240 to 400 days in t… ▽ More The giant spiral galaxy M101 is host to the nearest recent Type Ia Supernova (SN 2011fe) and thus has been extensively monitored in the near-infrared to study the late-time lightcurve of the supernova. Leveraging this existing baseline of observations, we derive the first Mira-based distance to M101 by discovering and classifying a sample of 211 Miras with periods ranging from 240 to 400 days in the supernova field. Combined with new HST WFC3/IR channel observations, our dataset totals 11 epochs of F110W (HST $YJ$) and 13 epochs of F160W (HST $H$) data spanning $\sim$2900 days. We adopt absolute calibrations of the Mira Period-Luminosity Relation based on geometric distances to the Large Magellanic Cloud and the water megamaser host galaxy NGC 4258, and find $μ_{\rm M101} = $ 29.10 $\pm$ 0.06 mag. This distance is in 1$σ$ agreement with most other recent Cepheid and Tip of the Red Giant Branch distance measurements to M101. Including the previous Mira-SNIa host, NGC 1559 and SN 2005df, we determine the fiducial SN Ia peak luminosity, $M^0_B = -19.27 \pm 0.09$ mag. With the Hubble diagram of SNe Ia, we derive $H_0 = 72.37 \pm 2.97 $ km s$^{-1}$Mpc$^{-1}$, a $4.1\%$ measurement of $H_0$ using Miras. We find excellent agreement with recent Cepheid distance ladder measurements of $H_0$ and confirm previous indications that the local universe value of $H_0$ is higher than the early-universe value at $\sim$ $95\%$ confidence. Currently, the Mira-based $H_0$ measurement is still dominated by the statistical uncertainty in the SN Ia peak magnitude. △ Less

Submitted 31 January, 2024; v1 submitted 13 December, 2023; originally announced December 2023.

Comments: 22 pages, 11 figures, accepted to ApJ

[OIII] 5007 emissions in extremely red quasars (ERQs) are compact

Authors: Marie Wingyee Lau, Serena Perrotta, Fred Hamann, Jarred Gillette, David S. N. Rupke, Andrey Vayner, Nadia L. Zakamska, Dominika Wylezalek

Abstract: ``Extremely red quasars'' (ERQs) are a non-radio-selected, intrinsically luminous population of quasars at cosmic noon selected by their extremely red colour from rest-frame UV to mid-IR. ERQs are uniquely associated with exceptionally broad and blueshifted [OIII] 5007 emission reaching speeds >6000 km s^-1. We obtained adaptive optics integral-field spectroscopic observations using Keck/OSIRIS an… ▽ More ``Extremely red quasars'' (ERQs) are a non-radio-selected, intrinsically luminous population of quasars at cosmic noon selected by their extremely red colour from rest-frame UV to mid-IR. ERQs are uniquely associated with exceptionally broad and blueshifted [OIII] 5007 emission reaching speeds >6000 km s^-1. We obtained adaptive optics integral-field spectroscopic observations using Keck/OSIRIS and Gemini/NIFS of a sample of 10 ERQs with bolometric luminosities (10^47.0-10^47.9) erg s^-1 at z ~(2.3-3.0). The goal is to measure the sizes and spatially-resolved kinematics of the [OIII]-emitting regions. We study the surface brightness maps and aperture-extracted spectra and model the point-spread functions. We identify signs of merger activities in the continuum emissions. We identify physically distinct [OIII] kinematic components that are bimodal and respectively trace ERQ-driven outflows of velocity dispersion >250 km s^-1 and dynamically quiescent interstellar media. We find that the ERQ-driven ionized outflows are typically at ~1 kpc scales whereas the quiescent ionized gas extend to a few kpc. Compared to normal quasars the extremely fast ERQ-driven [OIII] outflows tend to be more compact, supporting the notion that ERQs are in a young stage of quasar/galaxy evolution and represent systems with unique physical conditions beyond orientation differences with normal quasar populations. The kinematically quiescent [OIII] emissions in ERQs tend to be spatially-resolved but less extended than in normal quasars, which can be explained by global and patchy dust obscuration. The hint of ionization cones suggests some of the obscuration can be partially explained by a patchy torus. △ Less

Submitted 27 June, 2024; v1 submitted 6 December, 2023; originally announced December 2023.

Comments: Accepted by MNRAS

Measuring The Mass-Radius Relation of White Dwarfs Using Wide Binaries

Authors: Stefan Arseneau, Vedant Chandra, Hsiang-Chih Hwang, Nadia L. Zakamska, Gautham Adamane Pallathadka, Nicole R. Crumpler, J. J. Hermes, Kareem El-Badry, Hans-Walter Rix, Keivan G. Stassun, Boris T. Gaensicke, Joel R. Brownstein, Sean Morrison

Abstract: Measuring the mass-radius relation of individual white dwarfs is an empirically challenging task that has been performed for only a few dozen stars. We measure the white dwarf mass-radius relation using gravitational redshifts and radii of 137 white dwarfs in wide binaries with main sequence companions. We obtain the space velocities to these systems using the main sequence companion, and subtract… ▽ More Measuring the mass-radius relation of individual white dwarfs is an empirically challenging task that has been performed for only a few dozen stars. We measure the white dwarf mass-radius relation using gravitational redshifts and radii of 137 white dwarfs in wide binaries with main sequence companions. We obtain the space velocities to these systems using the main sequence companion, and subtract these Doppler redshifts from the white dwarfs' apparent motions, isolating their gravitational redshifts. We use Gaia data to calculate the surface temperatures and radii of these white dwarfs, thereby deriving an empirical gravitational redshift-radius relation. This work demonstrates the utility of low-resolution Galactic surveys to measure the white dwarf equation of state. Our results are consistent with theoretical models, and represent the largest sample of individual white dwarf gravitational redshift measurements to date. △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: 12 pages, 9 figures

Discovery of a proto-white dwarf with a massive unseen companion

Authors: Gautham Adamane Pallathadka, Vedant Chandra, Nadia L. Zakamska, Hsiang-Chih Hwang, Yossef Zenati, J. J. Hermes, Kareem El-Badry, Boris T. Gaensicke, Sean Morrison, Nicole R. Crumpler, Stefan Arseneau

Abstract: We report the discovery of SDSS~J022932.28+713002.7, a nascent extremely low-mass (ELM) white dwarf (WD) orbiting a massive ($> 1\,M_\odot$ at 2$σ$ confidence) companion with a period of 36 hours. We use a combination of spectroscopy, including data from the ongoing SDSS-V survey, and photometry to measure the stellar parameters for the primary pre-ELM white dwarf. The lightcurve of the primary WD… ▽ More We report the discovery of SDSS~J022932.28+713002.7, a nascent extremely low-mass (ELM) white dwarf (WD) orbiting a massive ($> 1\,M_\odot$ at 2$σ$ confidence) companion with a period of 36 hours. We use a combination of spectroscopy, including data from the ongoing SDSS-V survey, and photometry to measure the stellar parameters for the primary pre-ELM white dwarf. The lightcurve of the primary WD exhibits ellipsoidal variation, which we combine with radial velocity data and $\tt{PHOEBE}$ binary simulations to estimate the mass of the invisible companion. We find that the primary WD has mass $M_1$ = $0.18^{+0.02}_{-0.02}$ M$_\odot$ and the unseen secondary has mass $M_2$ = $1.19^{+0.21}_{-0.14}$ M$_\odot$. The mass of the companion suggests that it is most likely a near-Chandrasekhar mass white dwarf or a neutron star. It is likely that the system recently went through a Roche lobe overflow from the visible primary onto the invisible secondary. The dynamical configuration of the binary is consistent with the theoretical evolutionary tracks for such objects, and the primary is currently in its contraction phase. The measured orbital period puts this system on a stable evolutionary path which, within a few Gyrs, will lead to a contracted ELM white dwarf orbiting a massive compact companion. △ Less

Submitted 24 October, 2023; originally announced October 2023.

Comments: 21 Pages, 8 Figures

First results from the JWST Early Release Science Program Q3D: Powerful quasar-driven galactic scale outflow at $z=3$

Authors: Andrey Vayner, Nadia L. Zakamska, Yuzo Ishikawa, Swetha Sankar, Dominika Wylezalek, David S. N. Rupke, Sylvain Veilleux, Caroline Bertemes, Jorge K. Barrera-Ballesteros, Hsiao-Wen Chen, Nadiia Diachenko, Andy D. Goulding, Jenny E. Greene, Kevin N. Hainline, Fred Hamann, Timothy Heckman, Sean D. Johnson, Hui Xian Grace Lim, Weizhe Liu, Dieter Lutz, Nora Lutzgendorf, Vincenzo Mainieri, Ryan McCrory, Grey Murphree, Nicole P. H. Nesvadba , et al. (3 additional authors not shown)

Abstract: Quasar-driven galactic outflows are a major driver of the evolution of massive galaxies. We report observations of a powerful galactic-scale outflow in a $z=3$ extremely red, intrinsically luminous ($L_{\rm bol}\simeq 5\times 10^{47}$erg s$^{-1}$) quasar SDSSJ1652+1728 with the Near Infrared Spectrograph (NIRSpec) on board JWST. We analyze the kinematics of rest-frame optical emission lines and id… ▽ More Quasar-driven galactic outflows are a major driver of the evolution of massive galaxies. We report observations of a powerful galactic-scale outflow in a $z=3$ extremely red, intrinsically luminous ($L_{\rm bol}\simeq 5\times 10^{47}$erg s$^{-1}$) quasar SDSSJ1652+1728 with the Near Infrared Spectrograph (NIRSpec) on board JWST. We analyze the kinematics of rest-frame optical emission lines and identify the quasar-driven outflow extending out to $\sim 10$ kpc from the quasar with a velocity offset of ($v_{r}=\pm 500$ km s$^{-1}$) and high velocity dispersion (FWHM$=700-2400$ km s$^{-1}$). Due to JWST's unprecedented surface brightness sensitivity in the near-infrared -- we unambiguously show that the powerful high velocity outflow in an extremely red quasar (ERQ) encompasses a large swath of the host galaxy's interstellar medium (ISM). Using the kinematics and dynamics of optical emission lines, we estimate the mass outflow rate -- in the warm ionized phase alone -- to be at least $2300\pm1400$ $M_{\odot}$ yr$^{-1}$. We measure a momentum flux ratio between the outflow and the quasar accretion disk of $\sim$1 on kpc scale, indicating that the outflow was likely driven in a relatively high ($>10^{23}$cm$^{-2}$) column density environment through radiation pressure on dust grains. We find a coupling efficiency between the bolometric luminosity of the quasar and the outflow of 0.1$\%$, matching the theoretical prediction of the minimum coupling efficiency necessary for negative quasar feedback. The outflow has sufficient energetics to drive the observed turbulence seen in shocked regions of the quasar host galaxy, likely directly responsible for prolonging the time it takes for gas to cool efficiently. △ Less

Submitted 25 July, 2023; originally announced July 2023.

Comments: 14 pages, 7 figures, 1 table, submitted for publication in ApJ

Roman Early-Definition Astrophysics Survey Opportunity: Galactic Roman Infrared Plane Survey (GRIPS)

Authors: Roberta Paladini, Catherine Zucker, Robert Benjamin, David Nataf, Dante Minniti, Gail Zasowski, Joshua Peek, Sean Carey, Lori Allen, Javier Alonso-Garcia, Joao Alves, Friederich Anders, Evangelie Athanassoula, Timothy C. Beers, Jonathan Bird, Joss Bland-Hwathorn, Anthony Brown, Sven Buder, Luca Casagrande, Andrew Casey, Santi Cassisi, Marcio Catelan, Ranga-Ram Chary, Andre-Nicolas Chene, David Ciardi , et al. (45 additional authors not shown)

Abstract: A wide-field near-infrared survey of the Galactic disk and bulge/bar(s) is supported by a large representation of the community of Galactic astronomers. The combination of sensitivity, angular resolution and large field of view make Roman uniquely able to study the crowded and highly extincted lines of sight in the Galactic plane. A ~1000 deg2 survey of the bulge and inner Galactic disk would yiel… ▽ More A wide-field near-infrared survey of the Galactic disk and bulge/bar(s) is supported by a large representation of the community of Galactic astronomers. The combination of sensitivity, angular resolution and large field of view make Roman uniquely able to study the crowded and highly extincted lines of sight in the Galactic plane. A ~1000 deg2 survey of the bulge and inner Galactic disk would yield an impressive dataset of ~120 billion sources and map the structure of our Galaxy. The effort would foster subsequent expansions in numerous dimensions (spatial, depth, wavelengths, epochs). Importantly, the survey would benefit from early defintion by the community, namely because the Galactic disk is a complex environment, and different science goals will require trade offs. △ Less

Submitted 14 July, 2023; originally announced July 2023.

Comments: Submitted to the Roman Project on October 22 2021 in response to a call for white papers on early-definition Astrophysics opportunity

Varstrometry for Off-nucleus and Dual Sub-kpc AGN (VODKA): Very Long Baseline Array Searches for Dual or Off-nucleus Quasars and Small-scale Jets

Authors: Yu-Ching Chen, Xin Liu, Joseph Lazio, Peter Breiding, Sarah Burke-Spolaor, Hsiang-Chih Hwang, Yue Shen, Nadia L. Zakamska

Abstract: Dual and off-nucleus active supermassive black holes are expected to be common in the hierarchical structure formation paradigm, but their identification at parsec scales remains a challenge due to strict angular resolution requirements. We conduct a systematic study using the Very Long Baseline Array (VLBA) to examine 23 radio-bright candidate dual and off-nucleus quasars. The targets are selecte… ▽ More Dual and off-nucleus active supermassive black holes are expected to be common in the hierarchical structure formation paradigm, but their identification at parsec scales remains a challenge due to strict angular resolution requirements. We conduct a systematic study using the Very Long Baseline Array (VLBA) to examine 23 radio-bright candidate dual and off-nucleus quasars. The targets are selected by a novel astrometric technique ("varstrometry") from Gaia, aiming to identify dual or off-nucleus quasars at (sub)kilo-parsec scales. Among these quasars, 8 exhibit either multiple radio components or significant (>3$σ$) positional offsets between the VLBA and Gaia positions. The radio emission from the three candidates which exhibit multiple radio components is likely to originate from small-scale jets based on their morphology. Among the remaining five candidates with significant VLBA-Gaia offsets, three are identified as potential dual quasars at parsec scales, one is likely attributed to small-scale jets, and the origin of the last candidate remains unclear. We explore alternative explanations for the observed VLBA-Gaia offsets. We find no evidence for optical jets at kilo-parsec scales, nor any contamination to Gaia astrometric noise from the host galaxy; misaligned coordinate systems are unlikely to account for our offsets. Our study highlights the promise of the varstrometry technique in discovering candidate dual or off-nucleus quasars and emphasizes the need for further confirmation and investigation to validate and understand these intriguing candidates. △ Less

Submitted 10 October, 2023; v1 submitted 12 July, 2023; originally announced July 2023.

Comments: 17 pages, 8 figures, 4 tables, accepted to ApJ

Discovery of spectacular quasar-driven superbubbles in red quasars

Authors: Lu Shen, Guilin Liu, Zhicheng He, Nadia L. Zakamska, Eilat Glikman, Jenny E. Greene, Weida Hu, Guobin Mou, Dominika Wylezalek, David S. N. Rupke

Abstract: Quasar-driven outflows on galactic scales are a routinely invoked ingredient for galaxy formation models. We report the discovery of ionized gas nebulae as traced by [O III] $λ$5007 AA emission surrounding three luminous red quasars at $z \sim 0.4$ from Gemini Integral Field Unit (IFU) observations. All these nebulae feature unprecedented pairs of "superbubbles" extending $\sim$20 kpc in diameter,… ▽ More Quasar-driven outflows on galactic scales are a routinely invoked ingredient for galaxy formation models. We report the discovery of ionized gas nebulae as traced by [O III] $λ$5007 AA emission surrounding three luminous red quasars at $z \sim 0.4$ from Gemini Integral Field Unit (IFU) observations. All these nebulae feature unprecedented pairs of "superbubbles" extending $\sim$20 kpc in diameter, and the line-of-sight velocity difference between the red- and blue-shifted bubbles reaches up to $\sim$1200 km s$^{-1}$. Their spectacular dual-bubble morphology (in analogy to the Galactic "Fermi bubbles") and their kinematics provide unambiguous evidence for galaxy-wide quasar-driven outflows, in parallel with the quasi-spherical outflows similar in size from luminous Type-1 and -2 quasars at concordant redshift. These bubble pairs manifest themselves as a signpost of the short-lived superbubble ``break-out'' phase, when the quasar wind drives the bubbles to escape the confinement from the dense environment and plunge into the galactic halo with a high-velocity expansion. △ Less

Submitted 12 July, 2023; originally announced July 2023.

Comments: 47 pages, 15 figures, Published in Science Advances on Jul 12 2023

First results from the JWST Early Release Science Program Q3D: Benchmark Comparison of Optical and Mid-IR Tracers of a Dusty, Ionized Red Quasar Wind at z=0.435

Authors: D. S. N. Rupke, D. Wylezalek, N. L. Zakamska, S. Veilleux, C. Bertemes, Y. Ishikawa, W. Liu, S. Sankar, A. Vayner, H. X. G. Lim, R. McCrory, G. Murphree, L. Whitesell, L. Shen, G. Liu, J. K. Barrera-Ballesteros, H. -W. Chen, N. Diachenko, A. D. Goulding, J. E. Greene, K. N. Hainline, F. Hamann, T. Heckman, S. D. Johnson, D. Lutz , et al. (5 additional authors not shown)

Abstract: The [OIII] 5007 A emission line is the most common tracer of warm, ionized outflows in active galactic nuclei across cosmic time. JWST newly allows us to use mid-infrared spectral features at both high spatial and spectral resolution to probe these same winds. Here we present a comparison of ground-based, seeing-limited [OIII] and space-based, diffraction-limited [SIV] 10.51 micron maps of the pow… ▽ More The [OIII] 5007 A emission line is the most common tracer of warm, ionized outflows in active galactic nuclei across cosmic time. JWST newly allows us to use mid-infrared spectral features at both high spatial and spectral resolution to probe these same winds. Here we present a comparison of ground-based, seeing-limited [OIII] and space-based, diffraction-limited [SIV] 10.51 micron maps of the powerful, kiloparsec-scale outflow in the Type 1 red quasar SDSS J110648.32+480712.3. The JWST data are from the Mid-InfraRed Instrument (MIRI). There is a close match in resolution between the datasets (0."6), in ionization potential of the O$^{+2}$ and S$^{+3}$ ions (35 eV), and in line sensitivity (1e-17 to 2e-17 erg/s/cm$^2$/arcsec$^2$). The [OIII] and [SIV] line shapes match in velocity and linewidth over much of the 20 kpc outflowing nebula, and [SIV] is the brightest line in the rest-frame 3.5-19.5 micron range, demonstrating its usefulness as a mid-IR probe of quasar outflows. [OIII] is nevertheless intriniscally brighter and provides better contrast with the point-source continuum, which is strong in the mid-IR. There is a strong anticorrelation of [OIII]/[SIV] with average velocity, which is consistent with a scenario of differential obscuration between the approaching (blueshifted) and receding (redshifted) sides of the flow. The dust in the wind may also obscure the central quasar, consistent with models that attribute red quasar extinction to dusty winds. △ Less

Submitted 11 December, 2023; v1 submitted 21 June, 2023; originally announced June 2023.

Comments: resubmitted to match published version; minor changes

Journal ref: 2023ApJ...953L..26R

Varstrometry for Off-nucleus and Dual sub-Kpc AGN (VODKA). SDSS J1608+2716: A Sub-arcsec Quadruply Lensed Quasar at z=2.575

Authors: Junyao Li, Xin Liu, Yue Shen, Masamune Oguri, Arran C. Gross, Nadia L. Zakamska, Yu-Ching Chen, Hsiang-Chih Hwang

Abstract: We report Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) deep IR (F160W) imaging of SDSS J1608+2716. This system, located at a redshift of z=2.575, was recently reported as a triple quasar candidate with subarcsecond separations ($\sim0.25''$) based on selection from Gaia astrometry and follow-up Keck adaptive optics-assisted integral field unit spectroscopy. Our new HST deep IR imaging r… ▽ More We report Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) deep IR (F160W) imaging of SDSS J1608+2716. This system, located at a redshift of z=2.575, was recently reported as a triple quasar candidate with subarcsecond separations ($\sim0.25''$) based on selection from Gaia astrometry and follow-up Keck adaptive optics-assisted integral field unit spectroscopy. Our new HST deep IR imaging reveals the presence of a fourth point-like component located $\sim0.9''$ away from the triple system. Additionally, we detect an edge-on disk galaxy located in between the four point sources. The entire system exhibits a characteristic cusp structure in the context of strong gravitational lensing, and the observed image configuration can be successfully reproduced using a lens model based on a singular isothermal ellipsoid mass profile. These findings indicate that this system is a quadruply lensed quasar. Our results highlight the challenges associated with identifying dual/multiple quasars on $\sim$kpc scales at high redshifts, and emphasize the crucial role of deep, high-resolution IR imaging in robustly confirming such systems. △ Less

Submitted 9 September, 2023; v1 submitted 21 June, 2023; originally announced June 2023.

Comments: 10 pages, 3 figures, accepted by ApJL

Torus skin outflow in a near-Eddington quasar revealed by spectropolarimetry

Authors: Nadia L. Zakamska, Rachael M. Alexandroff

Abstract: Even when the direct view toward the active nucleus is obscured, nuclear emission propagating along other directions can scatter off surrounding material, become polarized and reach the observer. Spectropolarimetry can thus be an important tool in investigating the circumnuclear geometry and kinematics of quasars on scales that cannot yet be probed via direct observations. Here we discuss an intri… ▽ More Even when the direct view toward the active nucleus is obscured, nuclear emission propagating along other directions can scatter off surrounding material, become polarized and reach the observer. Spectropolarimetry can thus be an important tool in investigating the circumnuclear geometry and kinematics of quasars on scales that cannot yet be probed via direct observations. Here we discuss an intriguing class of quasars where the polarization position angle swings by large amounts (90 deg) within an emission line. We investigate a kinematic model in which the scattering dust or electrons are in an axisymmetric outflow. We propagate Stokes parameters in a variety of geometries of emitter, scatterer and observer. We use these models to predict polarization fraction, line profiles and polarization position angles and compare them to observations. We demonstrate that the swinging polarization angle can be a result of the geometry of the outflow and the orientation of the observer. Polarization properties of a near-Eddington extremely red quasar SDSS J1652 can be successfully explained by a model in which the quasar is surrounded by a geometrically thick disk, whose `skin' is outflowing at 1000 km/s and acts as the scatterer on scales of a few tens of pc. The line of sight to the observer in this source is within or close to the skin of the torus, in agreement with multi-wavelength data. Spectropolarimetric data and models presented here strongly support the thick-disk geometry of circumnuclear material suggested by recent numerical simulations of high-rate accretion flows onto black holes. △ Less

Submitted 2 September, 2023; v1 submitted 9 June, 2023; originally announced June 2023.

Comments: 17 pages, accepted to MNRAS

Varstrometry for Off-nucleus and Dual sub-Kpc AGN (VODKA): Investigating the Nature of J0823+2418 at $z=1.81$: a Likely Lensed Quasar

Authors: Arran C. Gross, Yu-Ching Chen, Adi Foord, Xin Liu, Yue Shen, Masamune Oguri, Andy Goulding, Hsiang-Chih Hwang, Nadia Zakamska, Yilun Ma, Liam Nolan

Abstract: Dual quasars at small physical separations are an important precursor phase of galaxy mergers, ultimately leading to the coalescence of the two supermassive black holes. Starting from a sample of dual/lensed quasar candidates discovered using astrometric jitter in Gaia data, we present a pilot case study of one of the most promising yet puzzling candidate dual quasars at cosmic noon (z$\sim$1.8).… ▽ More Dual quasars at small physical separations are an important precursor phase of galaxy mergers, ultimately leading to the coalescence of the two supermassive black holes. Starting from a sample of dual/lensed quasar candidates discovered using astrometric jitter in Gaia data, we present a pilot case study of one of the most promising yet puzzling candidate dual quasars at cosmic noon (z$\sim$1.8). Using multi-wavelength imaging and spectroscopy from X-ray to radio, we test whether the J0823+2418 system is two individual quasars in a bound pair at separation$\sim$0.64'', or instead a single quasar being gravitationally lensed by a foreground galaxy. We find consistent flux ratios ($\sim$1.25-1.45) between the two sources in optical, NIR, UV, and radio, and thus similar spectral energy distributions, suggesting a strong lensing scenario. However, differences in the radio spectral index, as well as changing X-ray fluxes, hint at either a dual quasar with otherwise nearly identical properties, or perhaps lensing-based time lag of $\sim$3 days paired with intrinsic variability. We find with lens mass modeling that the relative NIR positions and magnitudes of the two quasars and a marginally detected central galaxy are consistent with strong lensing. Archival SDSS spectra likewise suggest a foreground absorber via Mg II absorption lines. We conclude that J0823+2418 is likely a lensed quasar, and therefore that the VODKA sample contains a population of these lensed systems (perhaps as high as 50%) as well as dual quasars. △ Less

Submitted 6 June, 2023; originally announced June 2023.

Comments: 25 pages, 14 figures, submitted to ApJ

Cataclysmic Variables from Sloan Digital Sky Survey V -- the search for period bouncers continues

Authors: K. Inight, Boris T. Gänsicke, A. Schwope, S. F. Anderson, C. Badenes, E. Breedt, V. Chandra, B. D. R. Davies, N. P. Gentile Fusillo, M. J. Green, J. J. Hermes, I. Achaica Huamani, H. Hwang, K. Knauff, J. Kurpas, K. S. Long, V. Malanushenko, S. Morrison, I. J. Quiroz C., G. N. Aichele Ramos, A. Roman-Lopes, M. R. Schreiber, A. Standke, L. Stütz, J. R. Thorstensen , et al. (3 additional authors not shown)

Abstract: SDSS-V is carrying out a dedicated survey for white dwarfs, single and in binaries, and we report the analysis of the spectroscopy of cataclysmic variables (CVs) and CV candidates obtained during the final plug plate observations of SDSS. We identify eight new CVs, spectroscopically confirm 53 and refute eleven published CV candidates, and we report 21 new or improved orbital periods. Combined wit… ▽ More SDSS-V is carrying out a dedicated survey for white dwarfs, single and in binaries, and we report the analysis of the spectroscopy of cataclysmic variables (CVs) and CV candidates obtained during the final plug plate observations of SDSS. We identify eight new CVs, spectroscopically confirm 53 and refute eleven published CV candidates, and we report 21 new or improved orbital periods. Combined with previously published data, the orbital period distribution of the SDSS-V CVs does not clearly exhibit a period gap. This is consistent with previous findings that spectroscopically identified CVs have a larger proportion of short-period systems compared to samples identified from photometric variability. Remarkably, despite a systematic search, we find very few period bouncers. We estimate the space density of period bouncers to be $\simeq0.2\times10^{-6}\,\mathrm{pc}^{-3}$, i.e. they represent only a few per cent of the total CV population. This suggests that during their final phase of evolution, CVs either destroy the donor, e.g. via a merger, or that they become detached and cease mass transfer. △ Less

Submitted 11 September, 2023; v1 submitted 22 May, 2023; originally announced May 2023.

Comments: Accepted by MNRAS. Includes machine readable list of CVs

Journal ref: 2023MNRAS.525.3597I

Infrared spectroscopic confirmation of z~2 photometrically-selected obscured quasars

Authors: Yuzo Ishikawa, Ben Wang, Nadia L. Zakamska, Gordon T. Richards, Joseph F. Hennawi, Angelica B. Rivera

Abstract: The census of obscured quasar populations is incomplete, and remains a major unsolved problem, especially at higher redshifts, where we expect a greater density of galaxy formation and quasar activity. We present Gemini GNIRS near-infrared spectroscopy of 24 luminous obscured quasar candidates from the Sloan Digital Sky Survey's Stripe 82 region. The targets were photometrically selected using a W… ▽ More The census of obscured quasar populations is incomplete, and remains a major unsolved problem, especially at higher redshifts, where we expect a greater density of galaxy formation and quasar activity. We present Gemini GNIRS near-infrared spectroscopy of 24 luminous obscured quasar candidates from the Sloan Digital Sky Survey's Stripe 82 region. The targets were photometrically selected using a WISE/W4 selection technique that is optimized to identify IR-bright and heavily-reddened/optically-obscured targets at $z>1$. We detect emission lines of ${\rm Hα}$, ${\rm Hβ}$, and/or ${\rm[ O~III]}$ in 23 sources allowing us to measure spectroscopic redshifts in the range $1<z<3$ with bolometric luminosities spanning $L=10^{46.3}-10^{47.3}$ erg s$^{-1}$. We observe broad $10^3-10^4$ km s$^{-1}$ Balmer emissions with large ${\rm Hα}/{\rm Hβ}$ ratios, and we directly observe a heavily reddened rest-frame optical continuum in several sources, suggesting high extinction ($A_V\sim7-20$ mag). Our observations demonstrate that such optical/infrared photometric selection successfully recovers high-redshift obscured quasars. The successful identification of previously undetected red, obscured high-redshift quasar candidates suggests that there are more obscured quasars yet to be discovered. △ Less

Submitted 4 April, 2023; originally announced April 2023.

Comments: Accepted by MNRAS, 13 pages, 8 figures, 2 tables

Compact and Quiescent Circumgalactic Medium and Ly$α$ Halos around Extremely Red Quasars (ERQs)

Authors: Jarred Gillette, Marie Wingyee Lau, Fred Hamann, Serena Perrotta, David S. N. Rupke, Dominika Wylezalek, Nadia L. Zakamska, Andrey Vayner

Abstract: Red quasars may represent a young stage of galaxy evolution that provide important feedback to their host galaxies. We are studying a population of extremely red quasars (ERQs) with exceptionally fast and powerful outflows, at median redshift $z$ = 2.6. We present Keck/KCWI integral field spectra of 11 ERQs, which have a median color $i-W3$ = 5.9~mag, median… ▽ More Red quasars may represent a young stage of galaxy evolution that provide important feedback to their host galaxies. We are studying a population of extremely red quasars (ERQs) with exceptionally fast and powerful outflows, at median redshift $z$ = 2.6. We present Keck/KCWI integral field spectra of 11 ERQs, which have a median color $i-W3$ = 5.9~mag, median $\left\langle L_{\text{bol}} \right\rangle$ $\approx$ 5 $\times$ $10^{47}$ erg s$^{-1}$, Ly$α$ halo luminosity $\left\langle L_{\text{halo}} \right\rangle$ $=$ 5 $\times$ $10^{43}$ erg s$^{-1}$, and maximum linear size $>128$ kpc. The ERQ halos are generally similar to those of blue quasars, following known trends with $L_{\text{bol}}$ in halo properties. ERQs have halo symmetries similar to Type-I blue quasars, suggesting Type-I spatial orientations. ERQ $\left\langle L_{\text{halo}} \right\rangle$ is $\sim$2 dex below blue quasars, which is marginal due to scatter, but consistent with obscuration lowering photon escape fractions. ERQ halos tend to have more compact and circularly symmetric inner regions than blue quasars, with median exponential scale lengths of $\sim$9 kpc, compared to $\sim$16 kpc for blue quasars. When we include the central regions not available in blue quasar studies (due to PSF problems), the true median ERQ halo scale length is just $\sim$6 kpc. ERQ halos are also kinematically quiet, with median velocity dispersion 293 km s$^{-1}$, consistent with expected virial speeds. Overall we find no evidence for feedback on circumgalactic scales, and the current episode of quasar activity, perhaps due to long outflow travel times, has not been around long enough to affect the circumgalactic medium. We confirm the narrow Ly$α$ emission spikes found in ERQ aperture spectra are halo features, and are useful for systemic redshifts and measuring outflow speeds in other features. △ Less

Submitted 22 March, 2023; originally announced March 2023.

Comments: 19 pages, 9 figures, submitted to Monthly Notices of the Royal Astronomical Society

First results from the JWST Early Release Science Program Q3D: The Warm Ionized Gas Outflow in z ~ 1.6 Quasar XID 2028 and its Impact on the Host Galaxy

Authors: Sylvain Veilleux, Weizhe Liu, Andrey Vayner, Dominika Wylezalek, David S. N. Rupke, Nadia L. Zakamska, Yuzo Ishikawa, Caroline Bertemes, Jorge K. Barrera-Ballesteros, Hsiao-Wen Chen, Nadiia Diachenko, Andy D. Goulding, Jenny E. Greene, Kevin N. Hainline, Fred Hamann, Timothy Heckman, Sean D. Johnson, Hui Xian Grace Lim, Dieter Lutz, Nora Lutzgendorf, Vincenzo Mainieri, Roberto Maiolino, Ryan McCrory, Grey Murphree, Nicole P. H. Nesvadba , et al. (4 additional authors not shown)

Abstract: Quasar feedback may regulate the growth of supermassive black holes, quench coeval star formation, and impact galaxy morphology and the circumgalactic medium. However, direct evidence for quasar feedback in action at the epoch of peak black hole accretion at z ~ 2 remains elusive. A good case in point is the z = 1.6 quasar WISEA J100211.29+013706.7 (XID 2028) where past analyses of the same ground… ▽ More Quasar feedback may regulate the growth of supermassive black holes, quench coeval star formation, and impact galaxy morphology and the circumgalactic medium. However, direct evidence for quasar feedback in action at the epoch of peak black hole accretion at z ~ 2 remains elusive. A good case in point is the z = 1.6 quasar WISEA J100211.29+013706.7 (XID 2028) where past analyses of the same ground-based data have come to different conclusions. Here we revisit this object with the integral field unit of the Near Infrared Spectrograph (NIRSpec) on board the James Webb Space Telescope (JWST) as part of Early Release Science program Q3D. The excellent angular resolution and sensitivity of the JWST data reveal new morphological and kinematic sub-structures in the outflowing gas plume. An analysis of the emission line ratios indicates that photoionization by the central quasar dominates the ionization state of the gas with no obvious sign for a major contribution from hot young stars anywhere in the host galaxy. Rest-frame near-ultraviolet emission aligned along the wide-angle cone of outflowing gas is interpreted as a scattering cone. The outflow has cleared a channel in the dusty host galaxy through which some of the quasar ionizing radiation is able to escape and heat the surrounding interstellar and circumgalactic media. The warm ionized outflow is not powerful enough to impact the host galaxy via mechanical feedback, but radiative feedback by the AGN, aided by the outflow, may help explain the unusually small molecular gas mass fraction in the galaxy host. △ Less

Submitted 22 June, 2023; v1 submitted 15 March, 2023; originally announced March 2023.

Comments: 17 pages, 9 figures, accepted for publication in The Astrophysical Journal

Spatially resolved observations of outflows in the radio loud AGN of UGC 8782

Authors: Rogemar A. Riffel, Rogério Riffel, Marina Bianchin, Thaisa Storchi-Bergmann, Gabriel Luan Souza de Oliveira, Nadia L. Zakamska

Abstract: We use optical Integral Field Spectroscopy (IFU) to study the gas emission structure and kinematics in the inner 3.4$\times$4.9 kpc$^2$ region of the galaxy UGC 8782 (3C 293), host of a radio loud Active Galactic Nucleus (AGN). The observations were performed with the GMOS-IFU on the Gemini North telescope, resulting in a spatial resolution of $\sim725$ pc at the distance of the galaxy. While the… ▽ More We use optical Integral Field Spectroscopy (IFU) to study the gas emission structure and kinematics in the inner 3.4$\times$4.9 kpc$^2$ region of the galaxy UGC 8782 (3C 293), host of a radio loud Active Galactic Nucleus (AGN). The observations were performed with the GMOS-IFU on the Gemini North telescope, resulting in a spatial resolution of $\sim725$ pc at the distance of the galaxy. While the stars present ordered rotation following the orientation of the large scale disc, the gas shows a disturbed kinematics. The emission-line profiles present two kinematic components: a narrow ($σ\lesssim200$ km s$^{-1}$) component associated with the gas in the disc of the galaxy and a broad ($σ\gtrsim200$ km s$^{-1}$) component produced by gas outflows. Emission-line ratio diagrams indicate that the gas in the disc is excited by the AGN radiation field, while the emission of the outflow includes additional contribution of shock excitation due to the interaction of the radio jet with the environment gas. Deviations from pure rotation, of up to 30 km s$^{-1}$, are observed in the disc component and likely produced by a previous merger event. The broad component is blueshifted by $\sim150-500$ km s$^{-1}$ relative to the systemic velocity of the galaxy in all locations. We construct radial profiles of the mass outflow rate and kinetic power of the ionized gas outflows, which have the maximum values at $\sim1$ kpc from the nucleus with peak values of $\dot{M}_{\rm out,ΔR}=0.5\pm0.1$ M$_\odot$ yr$^{-1}$ and $\dot{K}_{\rm out,ΔR} =$(6.8$\pm$1.1)$\times$10$^{41}$ erg s$^{-1}$. The kinetic coupling efficiency of these outflows is in the range of 1$-$3 per cent, indicating that they could be powerful enough to affect the star formation in the host galaxy as predicted by theoretical simulations. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: Accepted for publication by MNRAS

First results from the JWST Early Release Science Program Q3D: Ionization cone, clumpy star formation and shocks in a $z=3$ extremely red quasar host

Authors: Andrey Vayner, Nadia L. Zakamska, Yuzo Ishikawa, Swetha Sankar, Dominika Wylezalek, David S. N. Rupke, Sylvain Veilleux, Caroline Bertemes, Jorge K. Barrera-Ballesteros, Hsiao-Wen Chen, Nadiia Diachenko, Andy D. Goulding, Jenny E. Greene, Kevin N. Hainline, Fred Hamann, Timothy Heckman, Sean D. Johnson, Hui Xian Grace Lim, Weizhe Liu, Dieter Lutz, Nora Lutzgendorf, Vincenzo Mainieri, Ryan McCrory, Grey Murphree, Nicole P. H. Nesvadba , et al. (3 additional authors not shown)

Abstract: Massive galaxies formed most actively at redshifts $z=1-3$ during the period known as `cosmic noon.' Here we present an emission-line study of an extremely red quasar SDSSJ165202.64+172852.3 host galaxy at $z=2.94$, based on observations with the Near Infrared Spectrograph (NIRSpec) integral field unit (IFU) on board JWST. We use standard emission-line diagnostic ratios to map the sources of gas i… ▽ More Massive galaxies formed most actively at redshifts $z=1-3$ during the period known as `cosmic noon.' Here we present an emission-line study of an extremely red quasar SDSSJ165202.64+172852.3 host galaxy at $z=2.94$, based on observations with the Near Infrared Spectrograph (NIRSpec) integral field unit (IFU) on board JWST. We use standard emission-line diagnostic ratios to map the sources of gas ionization across the host and a swarm of companion galaxies. The quasar dominates the photoionization, but we also discover shock-excited regions orthogonal to the ionization cone and the quasar-driven outflow. These shocks could be merger-induced or -- more likely, given the presence of a powerful galactic-scale quasar outflow -- these are signatures of wide-angle outflows that can reach parts of the galaxy that are not directly illuminated by the quasar. Finally, the kinematically narrow emission associated with the host galaxy presents as a collection of 1 kpc-scale clumps forming stars at a rate of at least 200 $M_{\odot}$ yr$^{-1}$. The ISM within these clumps shows high electron densities, reaching up to 3,000 cm$^{-3}$ with metallicities ranging from half to a third solar with a positive metallicity gradient and V band extinctions up to 3 magnitudes. The star formation conditions are far more extreme in these regions than in local star-forming galaxies but consistent with that of massive galaxies at cosmic noon. JWST observations reveal an archetypical rapidly forming massive galaxy undergoing a merger, a clumpy starburst, an episode of obscured near-Eddington quasar activity, and an extremely powerful quasar outflow simultaneously. △ Less

Submitted 25 July, 2023; v1 submitted 13 March, 2023; originally announced March 2023.

Comments: 19 pages, 8 figures. Accepted for publication in ApJ

The AGNIFS survey: spatially resolved observations of hot molecular and ionised outflows in nearby active galaxies

Authors: R. A. Riffel, T. Storchi-Bergmann, R. Riffel, M. Bianchin, N. L. Zakamska, D. Ruschel-Dutra, M. C. Bentz, L. Burtscher, D. M. Crenshaw, L. G. Dahmer-Hahn, N. Z. Dametto, R. I. Davies, M. R. Diniz, T. C. Fischer, C. M. Harrison, V. Mainieri, M. Revalski, A. Rodriguez-Ardila, D. J. Rosario, A. J. Schonell

Abstract: We present the hot molecular and warm ionised gas kinematics for 33 nearby ($0.001\lesssim z\lesssim0.056$) X-ray selected active galaxies using the H$_2 2.1218 μ$m and Br$γ$ emission lines observed in the K-band with the Gemini Near-Infrared Field Spectrograph (NIFS). The observations cover the inner 0.04$-$2 kpc of each AGN at spatial resolutions of 4$-$250 pc with a velocity resolution of… ▽ More We present the hot molecular and warm ionised gas kinematics for 33 nearby ($0.001\lesssim z\lesssim0.056$) X-ray selected active galaxies using the H$_2 2.1218 μ$m and Br$γ$ emission lines observed in the K-band with the Gemini Near-Infrared Field Spectrograph (NIFS). The observations cover the inner 0.04$-$2 kpc of each AGN at spatial resolutions of 4$-$250 pc with a velocity resolution of $σ_{\rm inst}\approx$20 ${\rm km s^{-1}}$. We find that 31 objects (94 per cent) present a kinematically disturbed region (KDR) seen in ionised gas, while such regions are observed in hot molecular gas for 25 galaxies (76 per cent). We interpret the KDR as being due to outflows with masses of 10$^2-$10$^7$ M$_\odot$ and 10$^0-$10$^4$ M$_\odot$ for the ionised and hot molecular gas, respectively. The ranges of mass-outflow rates ($\dot{M}_{\rm out}$) and kinetic power ($\dot{E}_{\rm K}$) of the outflows are 10$^{-3}-$10$^{1}$ M$_\odot$yr$^{-1}$ and $\sim$10$^{37}$$-$10$^{43}$ erg s$^{-1}$ for the ionised gas outflows, and 10$^{-5}$$-$10$^{-2}$ M$_\odot$ yr$^{-1}$ and 10$^{35}$$-$10$^{39}$ erg s$^{-1}$ for the hot molecular gas outflows. The median coupling efficiency in our sample is $\dot{E}_{K}/L_{\rm bol}\approx1.8\times10^{-3}$ and the estimated momentum fluxes of the outflows suggest they are produced by radiation-pressure in low-density environment, with possible contribution from shocks. △ Less

Submitted 14 March, 2023; v1 submitted 22 February, 2023; originally announced February 2023.

Comments: 37 pages, published in MNRAS - A few typos in the text and in the label of Fg 1 were corrected in this version

Journal ref: MNRAS, 521, 1832 (2023)

Cold Mode Gas Accretion on Two Galaxy Groups at z$\sim$2

Authors: Andrey Vayner, Nadia L. Zakamska, Sanchit Sabhlok, Shelley A. Wright, Lee Armus, Norman Murray, Gregory Walth, Yuzo Ishikawa

Abstract: We present Keck Cosmic Web Imager (KCWI) integral field spectroscopy (IFS) observations of rest-frame UV emission lines $\rm Lyα$, C IV $λλ$ 1548 Å, 1550Åand He II 1640 Åobserved in the circumgalactic medium (CGM) of two $z=2$ radio-loud quasar host galaxies. We detect extended emission on 80-90 kpc scale in $\rm Lyα$ in both systems with C IV, and He II emission also detected out to 30-50 kpc. Al… ▽ More We present Keck Cosmic Web Imager (KCWI) integral field spectroscopy (IFS) observations of rest-frame UV emission lines $\rm Lyα$, C IV $λλ$ 1548 Å, 1550Åand He II 1640 Åobserved in the circumgalactic medium (CGM) of two $z=2$ radio-loud quasar host galaxies. We detect extended emission on 80-90 kpc scale in $\rm Lyα$ in both systems with C IV, and He II emission also detected out to 30-50 kpc. All emission lines show kinematics with a blue and redshifted gradient pattern consistent with velocities seen in massive dark matter halos and similar to kinematic patterns of inflowing gas seen in hydrodynamical simulations. Using the kinematics of both resolved $\rm Lyα$ emission and absorption, we can confirm that both kinematic structures are associated with accretion. Combining the KCWI data with molecular gas observations with Atacama Large Millimeter/submillimeter Array (ALMA) and high spatial resolution of ionized gas with Keck OSIRIS, we find that both quasar host galaxies reside in proto-group environments at $z=2$. We estimate $1-6\times10^{10}$M$_\odot$ of warm-ionized gas within 30-50 kpc from the quasar that is likely accreting onto the galaxy group. We estimate inflow rates of 60-200 M$_\odot$yr$^{-1}$, within an order of magnitude of the outflow rates in these systems. In the 4C 09.17 system, we detect narrow gas streams associated with satellite galaxies, potentially reminiscent of ram-pressure stripping seen in local galaxy groups and clusters. We find that the quasar host galaxies reside in dynamically complex environments, with ongoing mergers, gas accretion, ISM stripping, and outflows likely playing an important role in shaping the assembly and evolution of massive galaxies at cosmic noon. △ Less

Submitted 30 November, 2022; originally announced December 2022.

Comments: 24 pages, 11 figures, 6 tabes. Accepted for publication in MNRAS

CSS1603+19: a low-mass polar near the cataclysmic variable period minimum

Authors: Yiqi Liu, Hsiang-Chih Hwang, Nadia L. Zakamska, John R. Thorstensen

Abstract: CSS1603+19 is a cataclysmic variable (CV) with an orbital period of 81.96 min, near the minimal period of cataclysmic variables. It is unusual in having a strong mid-infrared excess inconsistent with thermal emission from a brown dwarf companion. Here we present time-resolved multi-wavelength observations of this system. WISE photometry indicates that the mid-infrared excess displays a one-magnitu… ▽ More CSS1603+19 is a cataclysmic variable (CV) with an orbital period of 81.96 min, near the minimal period of cataclysmic variables. It is unusual in having a strong mid-infrared excess inconsistent with thermal emission from a brown dwarf companion. Here we present time-resolved multi-wavelength observations of this system. WISE photometry indicates that the mid-infrared excess displays a one-magnitude eclipsing-like variability during the orbit. We obtained near-infrared and optical spectroscopy using Gemini, MDM and APO telescopes. Near-infrared spectra show possible cyclotron features indicating that the white dwarf has a magnetic field of about 5MG. Optical and near-infrared spectra display double-peaked emission lines, with both components showing strong radial velocity variations during the orbital period and with the broad component leading the narrow component stably by about 0.2 of the orbital phase. We construct a physical model informed by existing observations of the system and determine that one component likely originates from the accretion column onto the magnetized white dwarf in synchronous rotation with the orbital motion and the other from the Roche overflow point. This allows us to constrain the masses of the binary components to be $M_1>0.24 M_{\odot}$ for the white dwarf accretor and $M_2=0.0644\pm0.0074 M_\odot$ for the donor. We classify the system as an AM Herculis star, or a polar. It has likely completed its stint on the period gap, but has not yet gone through the period bounce. △ Less

Submitted 17 April, 2023; v1 submitted 27 November, 2022; originally announced November 2022.

First results from the JWST Early Release Science Program Q3D: Turbulent times in the life of a $z \sim 3$ extremely red quasar revealed by NIRSpec IFU

Authors: Dominika Wylezalek, Andrey Vayner, David S. N. Rupke, Nadia L. Zakamska, Sylvain Veilleux, Yuzo Ishikawa, Caroline Bertemes, Weizhe Liu, Jorge K. Barrera-Ballesteros, Hsiao-Wen Chen, Andy D. Goulding, Jenny E. Greene, Kevin N. Hainline, Nora Lützgendorf, Fred Hamann, Timothy Heckman, Sean D. Johnson, Dieter Lutz, Vincenzo Mainieri, Roberto Maiolino, Nicole P. H. Nesvadba, Patrick Ogle, Eckhard Sturm

Abstract: Extremely red quasars, with bolometric luminosities exceeding $10^{47}$ erg s$^{-1}$, are a fascinating high-redshift population that is absent in the local universe. They are the best candidates for supermassive black holes accreting at rates at or above the Eddington limit, and they are associated with the most rapid and powerful outflows of ionized gas known to date. They are also hosted by mas… ▽ More Extremely red quasars, with bolometric luminosities exceeding $10^{47}$ erg s$^{-1}$, are a fascinating high-redshift population that is absent in the local universe. They are the best candidates for supermassive black holes accreting at rates at or above the Eddington limit, and they are associated with the most rapid and powerful outflows of ionized gas known to date. They are also hosted by massive galaxies. Here we present the first integral field unit (IFU) observations of a high-redshift quasar obtained by the Near Infrared Spectrograph (NIRSpec) on board the James Webb Space Telescope (JWST), which targeted SDSSJ165202.64+172852.3, an extremely red quasar at $z=2.94$. JWST observations reveal extended ionized gas - as traced by [OIII]$λ$5007Å- in the host galaxy of the quasar, its outflow, and the circumgalactic medium. The complex morphology and kinematics imply that the quasar resides in a very dense environment with several interacting companion galaxies within projected distances of 10-15 kpc. The high density of the environment and the large velocities of the companion galaxies suggest that this system may represent the core of a forming cluster of galaxies. The system is a good candidate for a merger of two or more dark matter halos, each with a mass of a few $10^{13}$ M$_\odot$ and traces potentially one of the densest knots at $z\sim3$. △ Less

Submitted 18 October, 2022; originally announced October 2022.

Comments: ApJL in press, 11 pages, please also see associated ESA / STScI press release from Oct 20th 2022

A close quasar pair in a disk-disk galaxy merger at z = 2.17

Authors: Yu-Ching Chen, Xin Liu, Adi Foord, Yue Shen, Masamune Oguri, Nianyi Chen, Tiziana Di Matteo, Miguel Holgado, Hsiang-Chih Hwang, Nadia Zakamska

Abstract: Galaxy mergers produce pairs of supermassive black holes (SMBHs), which may be witnessed as dual quasars if both SMBHs are rapidly accreting. The kiloparsec (kpc)-scale separation represents a physical regime sufficiently close for merger-induced effects to be important yet wide enough to be directly resolvable with the facilities currently available. Whereas many kpc-scale dual active galactic nu… ▽ More Galaxy mergers produce pairs of supermassive black holes (SMBHs), which may be witnessed as dual quasars if both SMBHs are rapidly accreting. The kiloparsec (kpc)-scale separation represents a physical regime sufficiently close for merger-induced effects to be important yet wide enough to be directly resolvable with the facilities currently available. Whereas many kpc-scale dual active galactic nuclei--the low-luminosity counterparts of quasars--have been observed in low-redshift mergers, no unambiguous dual quasar is known at cosmic noon ($z\approx2$), the peak of global star formation and quasar activity. Here we report multiwavelength observations of SDSS J0749+2255 as a kpc-scale, dual-quasar system hosted by a galaxy merger at cosmic noon ($z=2.17$). We discover extended host galaxies associated with the much brighter compact quasar nuclei (separated by 0.46" or 3.8 kpc) and low-surface-brightness tidal features as evidence for galactic interactions. Unlike its low-redshift and low-luminosity counterparts, SDSS J0749+2255 is hosted by massive compact disk-dominated galaxies. The apparent lack of stellar bulges and the fact that SDSS J0749+2255 already follows the local SMBH mass--host stellar mass relation, suggest that at least some SMBHs may have formed before their host stellar bulges. While still at kpc-scale separations where the host-galaxy gravitational potential dominates, the two SMBHs may evolve into a gravitationally bound binary system in around 0.22 Gyr. △ Less

Submitted 5 April, 2023; v1 submitted 22 September, 2022; originally announced September 2022.

Comments: 88 pages, 18 figures, 8 tables; published

Journal ref: Nature, 616, 45 (2023)

Statistics of Galactic-Scale Quasar Pairs at Cosmic Noon

Authors: Yue Shen, Hsiang-Chih Hwang, Masamune Oguri, Nianyi Chen, Tiziana Di Matteo, Yueying Ni, Simeon Bird, Nadia Zakamska, Xin Liu, Yu-Ching Chen, Kaitlin M. Kratter

Abstract: The statistics of galactic-scale quasar pairs can elucidate our understanding of the dynamical evolution of supermassive black hole (SMBH) pairs, the duty cycles of quasar activity in mergers, or even the nature of dark matter, but have been challenging to measure at cosmic noon, the prime epoch of massive galaxy and SMBH formation. Here we measure a double quasar fraction of… ▽ More The statistics of galactic-scale quasar pairs can elucidate our understanding of the dynamical evolution of supermassive black hole (SMBH) pairs, the duty cycles of quasar activity in mergers, or even the nature of dark matter, but have been challenging to measure at cosmic noon, the prime epoch of massive galaxy and SMBH formation. Here we measure a double quasar fraction of $\sim 6.2\pm0.5\times 10^{-4}$ integrated over $\sim 0.3-3$ arcsec separations (projected physical separations of $\sim 3-30\,{\rm kpc}$ at $z\sim 2$) in luminous ($L_{\rm bol}>10^{45.8}\,{\rm erg\,s^{-1}}$) unobscured quasars at $1.5<z<3.5$, using Gaia EDR3-resolved pairs around SDSS DR16 quasars. The measurement was based on a sample of 60 Gaia-resolved double quasars (out of 487 Gaia pairs dominated by quasar+star superpositions) at these separations, corrected for pair completeness in Gaia, which we quantify as functions of pair separation, magnitude of the primary, and magnitude contrast. The double quasar fraction increases towards smaller separations by a factor of $\sim 5$ over these scales. The division between physical quasar pairs and lensed quasars in our sample is currently unknown, requiring dedicated follow-up observations (in particular, deep, sub-arcsec-resolution IR imaging for the closest pairs). Intriguingly, at this point the observed pair statistics are in rough agreement with theoretical predictions both for the lensed quasar population in mock catalogs and for dual quasars in cosmological hydrodynamic simulations. Upcoming wide-field imaging/spectroscopic space missions such as Euclid, CSST and Roman, combined with targeted follow-up observations, will conclusively measure the abundances and host galaxy properties of galactic-scale quasar pairs, offset AGNs, and sub-arcsec lensed quasars across cosmic time. △ Less

Submitted 9 August, 2022; originally announced August 2022.

Comments: 19 pages, 9 figures; submitted to ApJ

Journal ref: ApJ, 943, 38 (2023)

Wide twin binaries are extremely eccentric: evidence of twin binary formation in circumbinary disks

Authors: Hsiang-Chih Hwang, Kareem El-Badry, Hans-Walter Rix, Chris Hamilton, Yuan-Sen Ting, Nadia Zakamska

Abstract: The Gaia mission recently revealed an excess population of equal-mass "twin" wide binaries, with mass ratio $q\gtrsim 0.95$, extending to separations of at least 1000 AU. The origin of this population is an enigma: twin binaries are thought to form via correlated accretion in circumbinary disks, but the typical observed protostellar disks have radii of $\sim100$ AU, far smaller than the separation… ▽ More The Gaia mission recently revealed an excess population of equal-mass "twin" wide binaries, with mass ratio $q\gtrsim 0.95$, extending to separations of at least 1000 AU. The origin of this population is an enigma: twin binaries are thought to form via correlated accretion in circumbinary disks, but the typical observed protostellar disks have radii of $\sim100$ AU, far smaller than the separations of the widest twins. Here, we infer the eccentricity distribution of wide twins from the distribution of their $v$-$r$ angles, i.e., the angle between the components' separation and relative velocity vectors. We find that wide twins must be on extremely eccentric orbits. For the excess-twin population at 400-1000 AU, we infer a near-delta function excess of high-eccentricity system, with eccentricity $0.95 \lesssim e \leq 1$. These high eccentricities for wide twins imply pericenter distances of order $10$ AU and suggest that their orbits were scattered via dynamical interactions in their birth environments, consistent with a scenario in which twins are born in circumbinary disks and subsequently widened. These results further establish twin wide binaries as a distinct population and imply that wide twins can be used as a probe of the dynamical history of stellar populations. △ Less

Submitted 7 July, 2022; v1 submitted 11 May, 2022; originally announced May 2022.

Comments: accepted for publication in The Astrophysical Journal Letter

Probing the Inner Circumgalactic Medium and Quasar Illumination around the Reddest `Extremely Red Quasar' (ERQ)

Authors: Marie Wingyee Lau, Fred Hamann, Jarred Gillette, Serena Perrotta, David S. N. Rupke, Dominika Wylezalek, Nadia L. Zakamska

Abstract: Dusty quasars might be in a young stage of galaxy evolution with prominent quasar feedback. A recently discovered population of luminous, extremely red quasars at $z\sim$~2--4 has extreme spectral properties related to exceptionally powerful quasar-driven outflows. We present Keck/KCWI observations of the reddest known ERQ, at $z=$\,2.3184, with extremely fast [\ion{O}{III}]~$λ$5007 outflow at… ▽ More Dusty quasars might be in a young stage of galaxy evolution with prominent quasar feedback. A recently discovered population of luminous, extremely red quasars at $z\sim$~2--4 has extreme spectral properties related to exceptionally powerful quasar-driven outflows. We present Keck/KCWI observations of the reddest known ERQ, at $z=$\,2.3184, with extremely fast [\ion{O}{III}]~$λ$5007 outflow at $\sim$6000~km~s$^{-1}$. The Ly$α$ halo spans $\sim$100~kpc. The halo is kinematically quiet, with velocity dispersion $\sim$300~km~s$^{-1}$ and no broadening above the dark matter circular velocity down to the spatial resolution $\sim$6~kpc from the quasar. We detect spatially-resolved \ion{He}{II}~$λ$1640 and \ion{C}{IV}~$λ$1549 emissions with kinematics similar to the Ly$α$ halo and a narrow component in the [\ion{O}{III}]~$λ$5007. Quasar reddening acts as a coronagraph allowing views of the innermost halo. A narrow Ly$α$ spike in the quasar spectrum is inner halo emission, confirming the broad \ion{C}{IV}~$λ$1549 in the unresolved quasar is blueshifted by $2240$~km~s$^{-1}$ relative to the halo frame. We propose the inner halo is dominated by moderate-speed outflow driven in the past and the outer halo dominated by inflow. The high central concentration of the halo and the symmetric morphology of the inner region are consistent with the ERQ being in earlier evolutionary stage than blue quasars. The \ion{He}{II}~$λ$1640/Ly$α$ ratio of the inner halo and the asymmetry level of the overall halo are dissimilar to Type~II quasars, suggesting unique physical conditions for this ERQ that are beyond orientation differences from other quasar populations. We find no evidence of mechanical quasar feedback in the Ly$α$-emitting halo. △ Less

Submitted 2 August, 2022; v1 submitted 11 March, 2022; originally announced March 2022.

Comments: 20 pages, 18 figures, published in MNRAS

Enhancement of Double-Close-Binary Quadruples

Authors: Gavin B. Fezenko, Hsiang-Chih Hwang, Nadia L. Zakamska

Abstract: Double-close-binary quadruples (2+2 systems) are hierarchical systems of four stars where two short-period binary systems move around their common center of mass on a wider orbit. Using Gaia Early Data Release 3, we search for comoving pairs where both components are eclipsing binaries. We present eight 2+2 quadruple systems with inner orbital periods of $<$ 0.4 days and with outer separations of… ▽ More Double-close-binary quadruples (2+2 systems) are hierarchical systems of four stars where two short-period binary systems move around their common center of mass on a wider orbit. Using Gaia Early Data Release 3, we search for comoving pairs where both components are eclipsing binaries. We present eight 2+2 quadruple systems with inner orbital periods of $<$ 0.4 days and with outer separations of $\gtrsim 1000$ AU. All but one system are newly discovered by this work, and we catalog their orbital information measured from their light curves. We find that the occurrence rate of 2+2 quadruples is 7.3 $\pm$ 2.6 times higher than what is expected from random pairings of field stars. At most a factor of $\sim$2 enhancement may be explained by the age and metallicity dependence of the eclipsing binary fraction in the field stellar population. The remaining factor of $\sim$3 represents a genuine enhancement of the production of short-period binaries in wide-separation ($>10^3$ AU) pairs, suggesting a close binary formation channel that may be enhanced by the presence of wide companions. △ Less

Submitted 31 January, 2022; originally announced February 2022.

Comments: 19 pages, 6 figures, accepted for publication by MNRAS

The eccentricity distribution of wide binaries and their individual measurements

Authors: Hsiang-Chih Hwang, Yuan-Sen Ting, Nadia L. Zakamska

Abstract: Eccentricity of wide binaries is difficult to measure due to their long orbital periods. With Gaia's high-precision astrometric measurements, eccentricity of a wide binary can be constrained by the angle between the separation vector and the relative velocity vector (the $v$-$r$ angle). In this paper, by using the $v$-$r$ angles of wide binaries in Gaia Early Data Release 3, we develop a Bayesian… ▽ More Eccentricity of wide binaries is difficult to measure due to their long orbital periods. With Gaia's high-precision astrometric measurements, eccentricity of a wide binary can be constrained by the angle between the separation vector and the relative velocity vector (the $v$-$r$ angle). In this paper, by using the $v$-$r$ angles of wide binaries in Gaia Early Data Release 3, we develop a Bayesian approach to measure the eccentricity distribution as a function of binary separations. Furthermore, we infer the eccentricities of individual wide binaries and make them publicly available. Our results show that the eccentricity distribution of wide binaries at $10^2$ AU is close to uniform and becomes superthermal at $>10^{3}$ AU, suggesting two formation mechanisms dominating at different separation regimes. The close binary formation, most likely disk fragmentation, results in a uniform eccentricity distribution at $<10^{2}$ AU. The wide binary formation that leads to highly eccentric wide binaries at $>10^{3}$ AU may be turbulent fragmentation and/or the dynamical unfolding of compact triples. With Gaia, measuring eccentricities is now possible for a large number of wide binaries, opening a new window to understanding binary formation and evolution. △ Less

Submitted 8 March, 2022; v1 submitted 2 November, 2021; originally announced November 2021.

Comments: Accepted for publication in MNRAS. Fig. 7 is the key result. The online catalog of individual eccentricity measurements is available at: https://drive.google.com/file/d/11D4FeLxQjOaNfP8tw5DAKANt2A7zLbpn/view?usp=sharing. The codes and pre-computed grid are available at: https://github.com/HC-Hwang/Eccentricity-of-wide-binaries

Wide binaries from the H3 survey: the thick disk and halo have similar wide binary fractions

Authors: Hsiang-Chih Hwang, Yuan-Sen Ting, Charlie Conroy, Nadia L. Zakamska, Kareem El-Badry, Phillip Cargile, Dennis Zaritsky, Vedant Chandra, Jiwon Jesse Han, Joshua S. Speagle, Ana Bonaca

Abstract: Due to the different environments in the Milky Way's disk and halo, comparing wide binaries in the disk and halo is key to understanding wide binary formation and evolution. By using Gaia Early Data Release 3, we search for resolved wide binary companions in the H3 survey, a spectroscopic survey that has compiled $\sim$150,000 spectra for thick-disk and halo stars to date. We identify 800 high-con… ▽ More Due to the different environments in the Milky Way's disk and halo, comparing wide binaries in the disk and halo is key to understanding wide binary formation and evolution. By using Gaia Early Data Release 3, we search for resolved wide binary companions in the H3 survey, a spectroscopic survey that has compiled $\sim$150,000 spectra for thick-disk and halo stars to date. We identify 800 high-confidence (a contamination rate of 4%) wide binaries and two resolved triples, with binary separations mostly between $10^3$-$10^5$ AU and a lowest [Fe/H] of $-2.7$. Based on their Galactic kinematics, 33 of them are halo wide binaries, and most of those are associated with the accreted Gaia-Sausage-Enceladus galaxy. The wide binary fraction in the thick disk decreases toward the low metallicity end, consistent with the previous findings for the thin disk. Our key finding is that the halo wide binary fraction is consistent with the thick-disk stars at a fixed [Fe/H]. There is no significant dependence of the wide binary fraction on the $α$-captured abundance. Therefore, the wide binary fraction is mainly determined by the iron abundance, not their disk or halo origin nor the $α$-captured abundance. Our results suggest that the formation environments play a major role for the wide binary fraction, instead of other processes like radial migration that only apply to disk stars. △ Less

Submitted 7 March, 2022; v1 submitted 2 November, 2021; originally announced November 2021.

Comments: Accepted for publication in MNRAS. Fig. 7 is the key plot. Online data are available at: https://drive.google.com/file/d/1TXKEwKBHcP4j9tMeJpa2hWDFHsEc1Qd3/view?usp=sharing. Comments are welcome!

The SNIa Runaway LP 398-9: Detection of Circumstellar Material and Surface Rotation

Authors: Vedant Chandra, Hsiang-Chih Hwang, Nadia L. Zakamska, Simon Blouin, Andrew Swan, Thomas R. Marsh, Ken J. Shen, Boris T. Gänsicke, J. J. Hermes, Odelia Putterman, Evan B. Bauer, Evan Petrosky, Vikram S. Dhillon, Stuart P. Littlefair, Richard P. Ashley

Abstract: A promising progenitor scenario for Type Ia supernovae (SNeIa) is the thermonuclear detonation of a white dwarf in a close binary system with another white dwarf. After the primary star explodes, the surviving donor can be spontaneously released as a hypervelocity runaway. One such runaway donor candidate is LP 398-9, whose orbital trajectory traces back $\approx 10^5$ years to a known supernova r… ▽ More A promising progenitor scenario for Type Ia supernovae (SNeIa) is the thermonuclear detonation of a white dwarf in a close binary system with another white dwarf. After the primary star explodes, the surviving donor can be spontaneously released as a hypervelocity runaway. One such runaway donor candidate is LP 398-9, whose orbital trajectory traces back $\approx 10^5$ years to a known supernova remnant. Here we report the discovery of carbon-rich circumstellar material around LP 398-9, revealed by a strong infrared excess and analyzed with follow-up spectroscopy. The circumstellar material is most plausibly composed of inflated layers from the star itself, mechanically and radioactively heated by the past companion's supernova. We also detect a 15.4 hr periodic signal in the UV and optical light curves of LP 398-9, which we interpret as surface rotation. The rotation rate is consistent with theoretical predictions from this supernova mechanism, and the brightness variations could originate from surface inhomogeneity deposited by the supernova itself. Our observations strengthen the case for this double-degenerate SNIa progenitor channel, and motivate the search for more runaway SNIa donors. △ Less

Submitted 13 April, 2022; v1 submitted 13 October, 2021; originally announced October 2021.

Comments: Accepted to MNRAS

Multi-phase outflows in high redshift quasar host galaxies

Authors: Andrey Vayner, Nadia Zakamska, Shelley A. Wright, Lee Armus, Norman Murray, Gregory Walth

Abstract: We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of six radio-loud quasar host galaxies at $z=1.4-2.3$. We combine the kpc-scale resolution ALMA observations with high spatial-resolution adaptive optics integral field spectrograph data of the ionized gas. We detect molecular gas emission in five quasar host galaxies and resolve the molecular interstellar medium using the… ▽ More We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of six radio-loud quasar host galaxies at $z=1.4-2.3$. We combine the kpc-scale resolution ALMA observations with high spatial-resolution adaptive optics integral field spectrograph data of the ionized gas. We detect molecular gas emission in five quasar host galaxies and resolve the molecular interstellar medium using the CO (3-2) or CO (4-3) rotational transitions. Clumpy molecular outflows are detected in four quasar host galaxies and in a merger system 21 kpc away from one quasar. Between the ionized and cold-molecular gas phases, the majority of the outflowing mass is in a molecular phase, while for three out of four detected multi-phase gas outflows, the majority of the kinetic luminosity and momentum flux is in the ionized phase. Combining the energetics of the multi-phase outflows, we find that their driving mechanism is consistent with energy-conserving shocks produced by the impact of the quasar jets with the gas in the galaxy. By assessing the molecular gas mass to the dynamics of the outflows, we estimate a molecular gas depletion time scale of a few Myr. The gas outflow rates exceed the star formation rates, suggesting that quasar feedback is a major mechanism of gas depletion at the present time. The coupling efficiency between the kinetic luminosity of the outflows and the bolometric luminosity of the quasar of 0.1-1% is consistent with theoretical predictions. Studying multi-phase gas outflows at high redshift is important for quantifying the impact of negative feedback in shaping the evolution of massive galaxies. △ Less

Submitted 30 September, 2021; originally announced October 2021.

Comments: 23 pages, 8 figures, 6 tabes. Accepted for publication in ApJ

A 99-minute Double-lined White Dwarf Binary from SDSS-V

Authors: Vedant Chandra, Hsiang-Chih Hwang, Nadia L. Zakamska, Boris T. Gaensicke, J. J. Hermes, Axel Schwope, Carles Badenes, Gagik Tovmassian, Evan B. Bauer, Dan Maoz, Matthias R. Schreiber, Odette F. Toloza, Keith P. Inight, Hans-Walter Rix, Warren R. Brown

Abstract: We report the discovery of SDSS J133725.26+395237.7 (hereafter SDSS J1337+3952), a double-lined white dwarf (WD+WD) binary identified in early data from the fifth generation Sloan Digital Sky Survey (SDSS-V). The double-lined nature of the system enables us to fully determine its orbital and stellar parameters with follow-up Gemini spectroscopy and Swift UVOT ultraviolet fluxes. The system is near… ▽ More We report the discovery of SDSS J133725.26+395237.7 (hereafter SDSS J1337+3952), a double-lined white dwarf (WD+WD) binary identified in early data from the fifth generation Sloan Digital Sky Survey (SDSS-V). The double-lined nature of the system enables us to fully determine its orbital and stellar parameters with follow-up Gemini spectroscopy and Swift UVOT ultraviolet fluxes. The system is nearby ($d = 113$ pc), and consists of a $0.51\, M_\odot$ primary and a $0.32\, M_\odot$ secondary. SDSS J1337+3952 is a powerful source of gravitational waves in the millihertz regime, and will be detectable by future space-based interferometers. Due to this gravitational wave emission, the binary orbit will shrink down to the point of interaction in $\approx 220$ Myr. The inferred stellar masses indicate that SDSS J1337+3952 will likely not explode as a Type Ia supernova (SN Ia). Instead, the system will probably merge and evolve into a rapidly rotating helium star, and could produce an under-luminous thermonuclear supernova along the way. The continuing search for similar systems in SDSS-V will grow the statistical sample of double-degenerate binaries across parameter space, constraining models of binary evolution and SNe Ia. △ Less

Submitted 26 August, 2021; originally announced August 2021.

Comments: Accepted to ApJ

Varstrometry for Off-nucleus and Dual sub-Kpc AGN (VODKA): Hubble Space Telescope Discovers Double Quasars

Authors: Yu-Ching Chen, Hsiang-Chih Hwang, Yue Shen, Xin Liu, Nadia L. Zakamska, Qian Yang, Jennifer I. Li

Abstract: Dual supermassive black holes (SMBHs) at $\sim$kpc scales are the progenitor population of SMBH mergers and play an important role in understanding the pairing and dynamical evolution of massive black holes in galaxy mergers. Because of the stringent resolution requirement and the apparent rareness of these small-separation pairs, there are scarce observational constraints on this population, with… ▽ More Dual supermassive black holes (SMBHs) at $\sim$kpc scales are the progenitor population of SMBH mergers and play an important role in understanding the pairing and dynamical evolution of massive black holes in galaxy mergers. Because of the stringent resolution requirement and the apparent rareness of these small-separation pairs, there are scarce observational constraints on this population, with few confirmed dual SMBHs at $<10$kpc separations at $z>1$. Here we present results from a pilot search for kpc-scale dual quasars selected with Gaia Data release 2 (DR2) astrometry and followed up with Hubble Space Telescope (HST) Wide Field Camera 3 dual-band (F475W and F814W) snapshot imaging. Our targets are quasars primarily selected with the varstrometry technique, i.e., light centroid jitter caused by asynchronous variability from both members in an unresolved quasar pair, supplemented by sub-arcsec pairs already resolved by Gaia DR2. We find an overall high fraction of HST-resolved pairs among the varstrometry-selected quasars (unresolved in Gaia DR2), $\sim 30-50\%$, increasing toward high redshift ($\sim 60-80\%$ at $z>1.5$). We discuss the nature of the 43 resolved sub-arcsec pairs based on HST and supplementary data. A substantial fraction ($\sim 40\%$) of these pairs are likely physical quasar pairs or gravitationally lensed quasars. We also discover a triple quasar candidate and a quadruply lensed quasar, which is among the smallest-separation quadruple lenses. These results provide important guidelines to improve varstrometry selection and follow-up confirmation of $\sim$kpc-scale dual SMBHs at high redshift. △ Less

Submitted 15 December, 2021; v1 submitted 3 August, 2021; originally announced August 2021.

Comments: 22 pages, 10 figures, 3 tables, accepted to ApJ

Journal ref: ApJ 925 162 (2022)

The impact of low luminosity AGN on their host galaxies: A radio and optical investigation of the kpc-scale outflow in MaNGA 1-166919

Authors: Aisha Al Yazeedi, Ivan Yu. Katkov, Joseph D. Gelfand, Dominika Wylezalek, Nadia L. Zakamska, Weizhe Liu

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

Submitted 15 May, 2021; originally announced May 2021.

Comments: 29 pages, 28 figures, 11 tables, accepted by The Astrophysical Journal for publication

A hidden population of high-redshift double quasars unveiled by astrometry

Authors: Yue Shen, Yu-Ching Chen, Hsiang-Chih Hwang, Xin Liu, Nadia Zakamska, Masamune Oguri, Jennifer I-Hsiu Li, Joseph Lazio, Peter Breiding

Abstract: Galaxy mergers occur frequently in the early universe and bring multiple supermassive black holes (SMBHs) into the nucleus, where they may eventually coalesce. Identifying post-merger-scale (i.e., <~a few kpc) dual SMBHs is a critical pathway to understanding their dynamical evolution and successive mergers. While serendipitously discovering kpc-scale dual SMBHs at z<1 is possible, such systems ar… ▽ More Galaxy mergers occur frequently in the early universe and bring multiple supermassive black holes (SMBHs) into the nucleus, where they may eventually coalesce. Identifying post-merger-scale (i.e., <~a few kpc) dual SMBHs is a critical pathway to understanding their dynamical evolution and successive mergers. While serendipitously discovering kpc-scale dual SMBHs at z<1 is possible, such systems are elusive at z>2, but critical to constraining the progenitors of SMBH mergers. The redshift z~2 also marks the epoch of peak activity of luminous quasars, hence probing this spatial regime at high redshift is of particular significance in understanding the evolution of quasars. However, given stringent resolution requirements, there is currently no confirmed <10 kpc physical SMBH pair at z>2. Here we report two sub-arcsec double quasars at z>2 discovered from a targeted search with a novel astrometric technique, demonstrating a high success rate (~50%) in this systematic approach. These high-redshift double quasars could be the long-sought kpc-scale dual SMBHs, or sub-arcsec gravitationally-lensed quasar images. One of these double quasars (at z=2.95) was spatially resolved with optical spectroscopy, and slightly favors the scenario of a physical quasar pair with a projected separation of 3.5 kpc (0.46"). Follow-up observations of double quasars discovered by this targeted approach will be able to provide the first observational constraints on kpc-scale dual SMBHs at z>2. △ Less

Submitted 7 May, 2021; originally announced May 2021.

Comments: Published in Nature Astronomy (April 2021)

Journal ref: Nature Astronomy, 5, 569 (2021)

The AGNIFS survey: distribution and excitation of the hot molecular and ionised gas in the inner kpc of nearby AGN hosts

Authors: R. A. Riffel, T. Storchi-Bergmann, R. Riffel, M. Bianchin, N. L. Zakamska, D. Ruschel-Dutra, A. J. Schonell, D. J. Rosario, A. Rodriguez-Ardila, T. C. Fischer, R. I. Davies, N. Z. Dametto, L. G. Dahmer-Hahn, D. M. Crenshaw, L. Burtscher, M. C. Bentz

Abstract: We use the Gemini NIFS instrument to map the H$_2 2.1218μ$m and Br$γ$ flux distributions in the inner 0.04-2 kpc of a sample of 36 nearby active galaxies ($0.001\lesssim z\lesssim0.056$) at spatial resolutions from 4 to 250 pc. We find extended emission in 34 galaxies. In $\sim$55% of them, the emission in both lines is most extended along the galaxy major axis, while in the other 45% the extent f… ▽ More We use the Gemini NIFS instrument to map the H$_2 2.1218μ$m and Br$γ$ flux distributions in the inner 0.04-2 kpc of a sample of 36 nearby active galaxies ($0.001\lesssim z\lesssim0.056$) at spatial resolutions from 4 to 250 pc. We find extended emission in 34 galaxies. In $\sim$55% of them, the emission in both lines is most extended along the galaxy major axis, while in the other 45% the extent follows a distinct orientation. The emission of H$_2$ is less concentrated than that of Br$γ$, presenting a radius that contains half of the flux 60% greater, on average. The H$_2$ emission is driven by thermal processes - X-ray heating and shocks - at most locations for all galaxies, where $0.4<H_2/Brγ<6$. For regions where H$_2$/Br$γ>6$ (seen in 40% of the galaxies), shocks are the main H$_2$ excitation mechanism, while in regions with H$_2$/Br$γ<0.4$ (25% of the sample) the H$_2$ emission is produced by fluorescence. The only difference we found between type 1 and type 2 AGN was in the nuclear emission-line equivalent widths, that are smaller in type 1 than in type 2 due to a larger contribution to the continuum from the hot dusty torus in the former. The gas masses in the inner 125 pc radius are in the range $10^1-10^4$ M$_\odot$ for the hot H$_2$ and $10^3-10^6$ M$_\odot$ for the ionised gas and would be enough to power the AGN in our sample for $10^5-10^8$ yr at their current accretion rates. △ Less

Submitted 7 April, 2021; originally announced April 2021.

Comments: 21 pages, accepted for publication in MNRAS

X-ray analysis of SDSS J165202.60+172852.4, an obscured quasar with outflows at peak galaxy formation epoch

Authors: Yuzo Ishikawa, Andy D. Goulding, Nadia L. Zakamska, Fred Hamann, Andrey Vayner, Sylvain Veilleux, Dominika Wylezalek

Abstract: We report on deep XMM-Newton and NuSTAR observations of the high redshift, z=2.94, extremely red quasar (ERQ), SDSS J165202.60+172852.4, with known galactic ionized outflows detected via spatially-resolved [OIII] emission lines. X-ray observations allow us to directly probe the accretion disk luminosity and the geometry and scale of the circumnuclear obscuration. We fit the spectra from the XMM-Ne… ▽ More We report on deep XMM-Newton and NuSTAR observations of the high redshift, z=2.94, extremely red quasar (ERQ), SDSS J165202.60+172852.4, with known galactic ionized outflows detected via spatially-resolved [OIII] emission lines. X-ray observations allow us to directly probe the accretion disk luminosity and the geometry and scale of the circumnuclear obscuration. We fit the spectra from the XMM-Newton/EPIC and NuSTAR detectors with a physically motivated torus model and constrain the source to exhibit a near Compton-thick column density of $N_H=(1.02^{+0.76}_{-0.41})\times10^{24}\textrm{ cm}^{-2}$, a near edge-on geometry with the line-of-sight inclination angle of $θ_i=85^{\circ}$, and a scattering fraction of $f_{sc}\sim 3$ %. The absorption-corrected, intrinsic 2-10 keV X-ray luminosity of $L_{\textrm{2-10}}=(1.4^{+1}_{-1}) \times10^{45} \textrm{ erg s}^{-1}$ reveals a powerful quasar that is not intrinsically X-ray weak, consistent with observed trends in other ERQs. We also estimate the physical properties of the obscuration, although highly uncertain: the warm ionized scattering density of $n_e \sim 7.5\times(10^2-10^3)\textrm{ cm}^{-3}$ and the obscuration mass of $M_{obsc} \sim 1.7\times(10^4-10^6) M_{\odot}$. As previously suggested with shallower X-ray observations, optical and infrared selection of ERQ has proved effective in finding obscured quasars with powerful outflow signatures. Our observations provide an in-depth view into the X-ray properties of ERQs and support the conclusions of severely photon-limited studies of obscured quasar populations at high redshifts. △ Less

Submitted 17 January, 2021; originally announced January 2021.

Comments: Accepted for publication in MNRAS. 12 pages, 6 Figures, 4 Tables

Powerful winds in high-redshift obscured and red quasars

Authors: Andrey Vayner, Nadia L. Zakamska, Rogemar A. Riffel, Rachael Alexandroff, Maren Cosens, Fred Hamann, Serena Perrotta, David S. N. Rupke, Thaisa Storchi Bergmann, Sylvain Veilleux, Greg Walth, Shelley Wright, Dominika Wylezalek

Abstract: Quasar-driven outflows must have made their most significant impact on galaxy formation during the epoch when massive galaxies were forming most rapidly. To study the impact of quasar feedback we conducted rest-frame optical integral field spectrograph (IFS) observations of three extremely red quasars (ERQs) and one type-2 quasar at $z=2-3$, obtained with the NIFS and OSIRIS instruments at the Gem… ▽ More Quasar-driven outflows must have made their most significant impact on galaxy formation during the epoch when massive galaxies were forming most rapidly. To study the impact of quasar feedback we conducted rest-frame optical integral field spectrograph (IFS) observations of three extremely red quasars (ERQs) and one type-2 quasar at $z=2-3$, obtained with the NIFS and OSIRIS instruments at the Gemini North and W. M. Keck Observatory with the assistance of laser-guided adaptive optics. We use the kinematics and morphologies of the [OIII] 5007Åand H$α$ 6563Åemission lines redshifted into the near-infrared to gauge the extents, kinetic energies and momentum fluxes of the ionized outflows in the quasars host galaxies. For the ERQs, the galactic-scale outflows are likely driven by radiation pressure in a high column density environment or due to an adiabatic shock. For the type-2 quasar, the outflow is driven by radiation pressure in a low column density environment or due to a radiative shock. The outflows in the ERQs carry a significant amount of energy ranging from 0.05-5$\%$ of the quasar's bolometric luminosity, powerful enough to have a significant impact on the quasar host galaxies. However, the outflows are likely only impacting the inner few kpc of each host galaxy. The observed outflow sizes are generally smaller than other ionized outflows observed at high redshift. The high ratio between the momentum flux of the ionized outflow and the photon momentum flux from the quasar accretion disk and high nuclear obscuration makes these ERQs great candidates for transitional objects where the outflows are likely responsible for clearing material in the inner regions of each galaxy, unveiling the quasar accretion disk at optical wavelengths. △ Less

Submitted 8 June, 2021; v1 submitted 12 January, 2021; originally announced January 2021.

Comments: 18 Pages, 12 figures, 2 tabes, published in MNRAS

Buckling bars in nearly face-on galaxies observed with MaNGA

Authors: K. M. Xiang, D. M. Nataf, E. Athanassoula, N. L. Zakamska, K. Rowlands, K. Masters, A. Fraser-McKelvie, N. Drory, K. Kraljic

Abstract: Over half of disk galaxies are barred, yet the mechanisms for bar formation and the life-time of bar buckling remain poorly understood. In simulations, a thin bar undergoes a rapid (<1 Gyr) event called "buckling," during which the inner part of the bar is asymmetrically bent out of the galaxy plane and eventually thickens, developing a peanut/X-shaped profile when viewed side-on. Through analyzin… ▽ More Over half of disk galaxies are barred, yet the mechanisms for bar formation and the life-time of bar buckling remain poorly understood. In simulations, a thin bar undergoes a rapid (<1 Gyr) event called "buckling," during which the inner part of the bar is asymmetrically bent out of the galaxy plane and eventually thickens, developing a peanut/X-shaped profile when viewed side-on. Through analyzing stellar kinematics of N-body model snapshots of a galaxy before, during, and after the buckling phase, we confirm a distinct quadrupolar pattern of out-of-plane stellar velocities in nearly face-on galaxies. This kinematic signature of buckling allows us to identify five candidates of currently buckling bars among 434 barred galaxies in the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) Survey, an integral field unit (IFU) spectroscopic survey that measures the composition and kinematic structure of nearby galaxies. The frequency of buckling events detected is consistent with the 0.5-1 Gyr timescale predicted by simulations. The five candidates we present more than double the total number of candidate buckling bars, and are the only ones found using the kinematic signature. △ Less

Submitted 11 January, 2021; originally announced January 2021.

Comments: 22 pages, 12 figures. Accepted for publication in ApJ

Variability, periodicity and contact binaries in WISE

Authors: Evan Petrosky, Hsiang-Chih Hwang, Nadia L. Zakamska, Vedant Chandra, Matthew J. Hill

Abstract: The time-series component of WISE is a valuable resource for the study of variable objects. We present an analysis of an all-sky sample of ~450,000 AllWISE+NEOWISE infrared light curves of likely variables identified in AllWISE. By computing periodograms of all these sources, we identify ~56,000 periodic variables. Of these, ~42,000 are short-period (P<1 day), near-contact or contact eclipsing bin… ▽ More The time-series component of WISE is a valuable resource for the study of variable objects. We present an analysis of an all-sky sample of ~450,000 AllWISE+NEOWISE infrared light curves of likely variables identified in AllWISE. By computing periodograms of all these sources, we identify ~56,000 periodic variables. Of these, ~42,000 are short-period (P<1 day), near-contact or contact eclipsing binaries, many of which are on the main sequence. We use the periodic and aperiodic variables to test computationally inexpensive methods of periodic variable classification and identification, utilizing various measures of the probability distribution function of fluxes and of timescales of variability. The combination of variability measures from our periodogram and non-parametric analyses with infrared colors from WISE and absolute magnitudes, colors and variability amplitude from Gaia is useful for the identification and classification of periodic variables. Furthermore, we show that the effectiveness of non-parametric methods for the identification of periodic variables is comparable to that of the periodogram but at a much lower computational cost. Future surveys can utilize these methods to accelerate more traditional time-series analyses and to identify evolving sources missed by periodogram-based selections. △ Less

Submitted 25 February, 2021; v1 submitted 8 December, 2020; originally announced December 2020.

Comments: 23 pages, 15 figures. Accepted for Publication in MNRAS. The full catalog of WISE variables, periodic variables, and binaries is available at https://zakamska.johnshopkins.edu/data.htm

Chemical abundances in Seyfert galaxies -- V. The discovery of shocked emission outside the AGN ionization axis

Authors: R. A. Riffel, O. L. Dors, M. Armah, T. Storchi-Bergmann, A. Feltre, G. F. Hagele, M. V. Cardaci, D. Ruschel-Dutra, A. C Krabbe, E. Perez-Montero, N. L. Zakamska, I. C. Freitas

Abstract: We present maps for the electron temperature in the inner kpc of three luminous Seyfert galaxies: Mrk 79, Mrk 348, and Mrk 607 obtained from Gemini GMOS-IFU observations at spatial resolutions of 110-280 pc. We study the distributions of electron temperature in active galaxies and find temperatures varying in the range from 8000 to >30000 K. Shocks due to gas outflows play an important role in the… ▽ More We present maps for the electron temperature in the inner kpc of three luminous Seyfert galaxies: Mrk 79, Mrk 348, and Mrk 607 obtained from Gemini GMOS-IFU observations at spatial resolutions of 110-280 pc. We study the distributions of electron temperature in active galaxies and find temperatures varying in the range from 8000 to >30000 K. Shocks due to gas outflows play an important role in the observed temperature distributions of Mrk 79 and Mrk 348, while standard photoionization models reproduce the derived temperature values for Mrk 607. In Mrk 79 and Mrk 348, we find direct evidence for shock-ionization with overall orientation orthogonal to the ionization axis, where shocks can be easily observed as the AGN radiation field is shielded by the nuclear dusty torus. This also indicates that even when the ionization cones are narrow, the shocks can be much wider-angle. △ Less

Submitted 3 December, 2020; v1 submitted 3 December, 2020; originally announced December 2020.

Comments: Accepted by MNRAS Letters

The non-monotonic, strong metallicity dependence of the wide-binary fraction

Authors: Hsiang-Chih Hwang, Yuan-Sen Ting, Kevin C. Schlaufman, Nadia L. Zakamska, Rosemary F. G. Wyse

Abstract: The metallicity dependence of the wide-binary fraction in stellar populations plays a critical role in resolving the open question of wide binary formation. In this paper, we investigate the metallicity ([Fe/H]) and age dependence of the wide-binary fraction (binary separations between $10^3$ and $10^4$ AU) for field F and G dwarfs within 500 pc by combining their metallicity and radial velocity m… ▽ More The metallicity dependence of the wide-binary fraction in stellar populations plays a critical role in resolving the open question of wide binary formation. In this paper, we investigate the metallicity ([Fe/H]) and age dependence of the wide-binary fraction (binary separations between $10^3$ and $10^4$ AU) for field F and G dwarfs within 500 pc by combining their metallicity and radial velocity measurements from LAMOST DR5 with the astrometric information from Gaia DR2. We show that the wide-binary fraction strongly depends on the metallicity: as metallicity increases, the wide-binary fraction first increases, peaks at [Fe/H]$\simeq 0$, and then decreases at the high metallicity end. The wide-binary fraction at [Fe/H]$=0$ is about two times larger than that at [Fe/H]$=-1$ and [Fe/H]$=+0.5$. This metallicity dependence is dominated by the thin-disk stars. Using stellar kinematics as a proxy of stellar age, we show that younger stars have a higher wide-binary fraction at fixed metallicity close to solar. We propose that multiple formation channels are responsible for the metallicity and age dependence. In particular, the positive metallicity correlation at [Fe/H]$<0$ and the age dependence may be due to the denser formation environments and higher-mass clusters at earlier times. The negative metallicity correlation at [Fe/H]$>0$ can be inherited from the similar metallicity dependence of close binaries, and radial migration may play a role in enhancing the wide-binary fraction around the solar metallicity. △ Less

Submitted 14 December, 2020; v1 submitted 6 October, 2020; originally announced October 2020.

Comments: Accepted by MNRAS. Fig 3 and 7 are the key results

A Gravitational Redshift Measurement of the White Dwarf Mass-Radius Relation

Authors: Vedant Chandra, Hsiang-Chih Hwang, Nadia L. Zakamska, Sihao Cheng

Abstract: The mass-radius relation of white dwarfs is largely determined by the equation of state of degenerate electrons, which causes the stellar radius to decrease as mass increases. Here we observationally measure this relation using the gravitational redshift effect, a prediction of general relativity that depends on the ratio between stellar mass and radius. Using observations of over three thousand w… ▽ More The mass-radius relation of white dwarfs is largely determined by the equation of state of degenerate electrons, which causes the stellar radius to decrease as mass increases. Here we observationally measure this relation using the gravitational redshift effect, a prediction of general relativity that depends on the ratio between stellar mass and radius. Using observations of over three thousand white dwarfs from the Sloan Digital Sky Survey and the Gaia space observatory, we derive apparent radial velocities from absorption lines, stellar radii from photometry and parallaxes, and surface gravities by fitting atmospheric models to spectra. By averaging the apparent radial velocities of white dwarfs with similar radii and, independently, surface gravities, we cancel out random Doppler shifts and measure the underlying gravitational redshift. Using these results, we empirically measure the white dwarf mass-radius relation across a wide range of stellar masses. Our results are consistent with leading theoretical models, and our methods could be used with future observations to empirically constrain white dwarf core composition and evolution. △ Less

Submitted 25 August, 2020; v1 submitted 28 July, 2020; originally announced July 2020.

Comments: 11 pages, 5 figures. Accepted to ApJ; corrected typos, added a reference and DOI

Computational Tools for the Spectroscopic Analysis of White Dwarfs

Authors: Vedant Chandra, Hsiang-Chih Hwang, Nadia L. Zakamska, Tamás Budavári

Abstract: The spectroscopic features of white dwarfs are formed in the thin upper layer of their stellar photosphere. These features carry information about the white dwarf's surface temperature, surface gravity, and chemical composition (hereafter 'labels'). Existing methods to determine these labels rely on complex ab-initio theoretical models which are not always publicly available. Here we present two t… ▽ More The spectroscopic features of white dwarfs are formed in the thin upper layer of their stellar photosphere. These features carry information about the white dwarf's surface temperature, surface gravity, and chemical composition (hereafter 'labels'). Existing methods to determine these labels rely on complex ab-initio theoretical models which are not always publicly available. Here we present two techniques to determine atmospheric labels from white dwarf spectra: a generative fitting pipeline that interpolates theoretical spectra with artificial neural networks, and a random forest regression model using parameters derived from absorption line features. We test and compare our methods using a large catalog of white dwarfs from the Sloan Digital Sky Survey (SDSS), achieving the same accuracy and negligible bias compared to previous studies. We package our techniques into an open-source Python module 'wdtools' that provides a computationally inexpensive way to determine stellar labels from white dwarf spectra observed from any facility. We will actively develop and update our tool as more theoretical models become publicly available. We discuss applications of our tool in its present form to identify interesting outlier white dwarf systems including those with magnetic fields, helium-rich atmospheres, and double-degenerate binaries. △ Less

Submitted 8 September, 2020; v1 submitted 22 July, 2020; originally announced July 2020.

Comments: 11 pages, 7 figures. Accepted for publication in MNRAS; updated references