Molecular Gas Kinematics and Star Formation Properties of the   Strongly-Lensed Quasar Host Galaxy RXS J1131-1231

T. K. Daisy Leung; Dominik A. Riechers; Riccardo Pavesi (Cornell)

arXiv:1701.07830·astro-ph.GA·February 20, 2017

Molecular Gas Kinematics and Star Formation Properties of the Strongly-Lensed Quasar Host Galaxy RXS J1131-1231

T. K. Daisy Leung, Dominik A. Riechers, Riccardo Pavesi (Cornell)

PDF

TL;DR

This study uses high-resolution observations of molecular gas in a strongly-lensed quasar host galaxy to reveal its rotating disk structure, star formation activity, and black hole growth, providing insights into galaxy evolution at intermediate redshift.

Contribution

It presents detailed lens modeling and molecular gas analysis of RXS J1131-1231, demonstrating the galaxy's rotating disk and star formation properties, and compares black hole and bulge growth at z~0.65.

Findings

01

The host galaxy has an extended rotating disk with a CO size of ~6 kpc.

02

Star formation rate is ~120 solar masses per year, similar to z~1 disk galaxies.

03

No evidence of star formation quenching; black hole growth outpaces bulge assembly.

Abstract

We report observations of CO(J=2-1) and CO(J=3-2) line emission towards the quadruply-lensed quasar RXS J1131-1231 at z = 0.654 obtained using the Plateau de Bure Interferometer (PdBI) and the Combined Array for Research in Millimeter-wave Astronomy (CARMA). Our lens modeling shows that the asymmetry in the double-horned CO(J = 2-1) line profile is mainly a result of differential lensing, where the magnification factor varies from ~3 to ~9 across different kinematic components. The intrinsically symmetric line profile and a smooth source-plane velocity gradient suggest that the host galaxy is an extended rotating disk, with a CO size of R_CO~6 kpc and a dynamical mass of M_dyn~8x10^10 Msun. We also find a secondary CO-emitting source near RXS J1131-1231 whose location is consistent with the optically-faint companion reported in previous studies. The lensing-corrected molecular gas…

Figures15

Click any figure to enlarge with its caption.

Tables6

Table 1. Table 1 Observed Properties of RXJ1131 and its companion

Parameter	Unit	Value
$z_{CO (2 - 1)}$		0.6541 $\pm$ 0.0002
$I_{CO (2 - 1)}$	Jy km s^-1	20.3 $\pm$ 0.6
$S_{CO (2 - 1)}$ ^a^aPeak flux density in the intensity map.	Jy km s^-1 beam^-1	8.12 $\pm$ 0.30
FWHM_CO(2-1) ^b^bFrom fitting a double Gaussian to the observed CO( $J$ = 2 $\to$ 1) spectrum (Figure 1).	km s^-1	179 $\pm$ 9, 255 $\pm$ 28
FWHM_CO(2-1) ^c^cFrom fitting a double Gaussian with a common FWHM to the de-lensed CO( $J$ = 2 $\to$ 1) spectrum (Figure 5).	km s^-1	220 $\pm$ 72
$I_{CO (3 - 2)}$	Jy km s^-1	35.7 $\pm$ 6.9

Table 2. Table 2 Photometry data

Wavelength	Frequency	Flux Density	Instrument
( $µm$ )	(GHz)	(mJy)
Combined/Unresolved
1.25	239834	1.009 $\pm$ 0.090	CTIO/J-Band
1.65	181692	1.448 $\pm$ 0.121	CTIO/H-Band
2.17	138153	2.064 $\pm$ 0.160	CTIO/Ks-Band
3.4	88174.2	7.027 $\pm$ 0.142	WISE/W1
3.6	83275.7	5.618 $\pm$ 0.002	Spitzer/IRAC
4.5	66620.5	7.803 $\pm$ 0.002	Spitzer/IRAC
4.6	65172.3	8.872 $\pm$ 0.163	WISE/W2
5.8	51688.4	10.720 $\pm$ 0.005	Spitzer/IRAC
8.0	37474.1	14.470 $\pm$ 0.004	Spitzer/IRAC
12	24982.7	21.960 $\pm$ 0.425	WISE/W3
12	24982.7	$<$ 400	IRAS
22	13626.9	55.110 $\pm$ 1.878	WISE/W4
24	12491.4	70.204 $\pm$ 0.026	Spitzer/MIPS
25	11991.7	$<$ 500	IRAS
60	4996.54	$<$ 600	IRAS
100	2997.92	$<$ 1000	IRAS
250	1199.17	289.4 $\pm$ 9.6	Herschel/SPIRE
350	856.55	168.2 $\pm$ 8.6	Herschel/SPIRE
500	599.585	56.8 $\pm$ 8.8	Herschel/SPIRE
1387.93	216	$<$ 2.49	CARMA
2152.82	139.256	1.23 $\pm$ 0.22	PdBI
Foreground Lensing Galaxy (deblended bands)
0.555	540167	0.056 $\pm$ 0.006	HST-ACS/V-Band
0.814	368295	0.238 $\pm$ 0.013	HST-ACS/I-Band
1.6	187370	0.539 $\pm$ 0.041	HST-NICMOS(NIC2)/H-Band
3.6	83275.7	0.585 $\pm$ 0.003 ^a^aFlux obtained using aperture photometry after subtracting the emission of RXJ1131 from the total emission.	Spitzer/IRAC
4.5	66620.5	1.794 $\pm$ 0.003 ^a^aFlux obtained using aperture photometry after subtracting the emission of RXJ1131 from the total emission.	Spitzer/IRAC
5.8	51688.4	3.163 $\pm$ 0.006 ^a^aFlux obtained using aperture photometry after subtracting the emission of RXJ1131 from the total emission.	Spitzer/IRAC
8.0	37474.1	4.589 $\pm$ 0.006 ^a^aFlux obtained using aperture photometry after subtracting the emission of RXJ1131 from the total emission.	Spitzer/IRAC
2152.82	139.256	0.799 $\pm$ 0.082	PdBI
61414	4.8815	0.866 $\pm$ 0.027	VLA
Background Galaxy RXJ1131 (deblended bands)
0.555	540167	0.009 $\pm$ 0.004 ^b^bA contribution from the quasar has been removed (see C06), and thus the flux density corresponds to the host galaxy only.	HST-ACS/V-Band
0.814	368295	0.041 $\pm$ 0.005 ^b^bA contribution from the quasar has been removed (see C06), and thus the flux density corresponds to the host galaxy only.	HST-ACS/I-Band
1.6	187370	0.133 $\pm$ 0.004 ^b^bA contribution from the quasar has been removed (see C06), and thus the flux density corresponds to the host galaxy only.	HST-NICMOS(NIC2)/H-Band
3.6	83275.7	5.034 $\pm$ 0.002	Spitzer/IRAC
4.5	66620.5	6.009 $\pm$ 0.002	Spitzer/IRAC
5.8	51688.4	7.557 $\pm$ 0.003	Spitzer/IRAC
8.0	37474.1	9.881 $\pm$ 0.004	Spitzer/IRAC
2152.82	139.256	0.39 $\pm$ 0.12 ^c^cFlux extracted from the residual map after subtracting a point-source model. For SED modeling, we use $S_{2 m m}$ = 0.27 $\pm$ 0.08 mJy to exclude synchrotron emission (see §3.3).	PdBI
61414	4.8815	1.273 $\pm$ 0.042	VLA

Table 3. Table 3 Lens parameters constrained by models of seven velocity channels

Parameters		Median values
Offset in RA	(″)	0.004 $\pm$ 0.027
Offset in Dec	(″)	0.003 $\pm$ 0.027
Axial Ratio		0.56 $\pm$ 0.16
Position Angle	(deg)	103 $\pm$ 22
Einstein Radius ^a^aThis corresponds to mass of $M (θ$ $<$ $θ_{E})$ = (7.42 $\pm$ 0.02) $\times 10^{11}$ $M_{⊙}$ within the Einstein radius.	(″)	1.833 $\pm$ 0.002

Table 4. Table 4 Magnification factors of various kinematic components in CO( J 𝐽 J = 2 → → \rightarrow 1)

Velocity Range (km s^-1)	Source 1 $μ_{L}$	Source 2 $μ_{L}$
346 $-$ 260	6.7 $\pm$ 2.5	7.2 $\pm$ 5.6
238 $-$ 153	7.6 $\pm$ 1.6
131 $-$ 45	8.7 $\pm$ 2.0
24 $-$ -62	4.1 $\pm$ 0.9
-84 $-$ -170	4.2 $\pm$ 0.6
-191 $-$ -277	4.3 $\pm$ 2.4
-300 $-$ -385	3.1 $\pm$ 0.9
weighted average	4.4
median	5.5

Table 5. Table 5 SED fitting results

Parameters		With 24µm	Without 24µm
$T_{d}$	(K)	54 ${}_{- 10}^{+ 8}$	55 ${}_{- 20}^{+ 21}$
$β$		1.6 ${}_{- 0.4}^{+ 0.5}$	2.2 ${}_{- 0.3}^{+ 0.3}$
$α$		1.6 ${}_{- 0.6}^{+ 0.5}$	8.5 ${}_{- 6.2}^{+ 7.0}$
$λ_{0}$ ^a^aObserved-frame wavelength where $τ_{ν}$ = 1	( $µm$ )	559 ${}_{- 324}^{+ 278}$	365 ${}_{- 120}^{+ 111}$
$λ_{peak}$ ^b^bObserved-frame wavelength of the SED peak	( $µm$ )	159 ${}_{- 40}^{+ 19}$	155 ${}_{- 43}^{+ 38}$
$f_{norm, 500 µm}$ ^c^cObserved-frame flux density at 500 $µm$	(mJy)	55 ${}_{- 13}^{+ 13}$	59 ${}_{- 6}^{+ 6}$
$L_{FIR}$ ^d^dRest-frame 42.5 $-$ 122.5 $µm$ luminosity	(10¹² $L_{⊙}$ )	3.81 ${}_{- 1.97}^{+ 1.92}$	4.24 ${}_{- 2.00}^{+ 2.17}$
$M_{d}$ ^e^eDerived assuming an absorption mass coefficient of $κ$ = 2.64 m² kg^-1 at $λ$ = 125.0 $µm$ (Dunne et al., 2003)	(10⁸ $M_{⊙}$ )	16 ${}_{- 12}^{+ 5}$	14 ${}_{- 7}^{+ 5}$

Table 6. Table 6 Physical properties of RXJ1131 and its companion

Parameter	Unit	Value
$r_{32}$		0.78 $\pm$ 0.37
FWHM_{CO(2-1),RXJ1131} ^a^aFrom fitting a double Gaussian with a common FWHM to the de-lensed spectrum.	km s^-1	220 $\pm$ 72
FWHM_{CO(2-1),RXJ1131} ^b^bFrom fitting a single Gaussian to the de-lensed spectrum.	km s^-1	600 $\pm$ 160
FWHM_{CO(2-1),companion} ^b^bFrom fitting a single Gaussian to the de-lensed spectrum.	km s^-1	73 $\pm$ 43
$M_{gas, RXJ1131}$	$10^{10}$ $M_{⊙}$	1.38 $\pm$ 0.33
$M_{gas, companion}$	$10^{9}$ $M_{⊙}$	1.92 $\pm$ 0.09
$R_{CO, RXJ1131}$	kpc	6.2 $\pm$ 3.0
$R_{CO, companion}$	kpc	4.2 $\pm$ 2.8
$M_{dyn, RXJ1131}$	$10^{10}$ $M_{⊙}$	8.3 $\pm$ 1.9 ^c^cExcluding systematic uncertainties.
$M_{dyn, companion}$	$10^{9}$ $M_{⊙}$	3.5 $\pm$ 2.3 ^c^cExcluding systematic uncertainties.
$f_{gas}$	%	18 $\pm$ 4 ^d^dExcluding uncertainties in the dynamical masses.
$f_{mol}$	%	34 $\pm$ 16
$L_{IR}$	$10^{12}$ $L_{⊙}$	$\sim$ 1.5
$L_{FIR}$	$10^{11}$ $L_{⊙}$	6.9 $\pm$ 3.6
SFR_FIR	$M_{⊙}$ yr^-1	120 $\pm$ 63
$M_{dust}$	$10^{8}$ $M_{⊙}$	$\sim$ 3
GDR		54 $\pm$ 13
$τ_{depl}$	Myr	102 $\pm$ 25
$M_{*}$	$10^{10}$ $M_{⊙}$	3.0 $\pm$ 1.0
$M_{BH}$ ^e^eSluse et al. (2012).	$10^{7}$ $M_{⊙}$	$\sim$ 8
$M_{BH} / M_{bulge}$	%	$>$ 0.27 ${}_{- 0.08}^{+ 0.11}$

Equations2

V = V_{0} + \frac{2}{π} V_{a} arctan (\frac{R}{R _{t}}),

V = V_{0} + \frac{2}{π} V_{a} arctan (\frac{R}{R _{t}}),

Peer Reviews

No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.

Videos

No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.

Full text

Molecular Gas Kinematics and Star Formation Properties of the Strongly-Lensed Quasar Host Galaxy RXS J1131 $-$ 1231

T. K. Daisy Leung, Dominik A. Riechers, and Riccardo Pavesi

Department of Astronomy, Space Sciences Building, Cornell University, Ithaca, NY 14853, USA; [email protected]

Abstract

We report observations of CO( $J$ = 2 $\rightarrow$ 1) and CO( $J$ = 3 $\rightarrow$ 2) line emission towards the quadruply-lensed quasar RXS J1131 $-$ 1231 at $z$ = 0.654 obtained using the Plateau de Bure Interferometer (PdBI) and the Combined Array for Research in Millimeter-wave Astronomy (CARMA). Our lens modeling shows that the asymmetry in the double-horned CO( $J$ = 2 $\rightarrow$ 1) line profile is mainly a result of differential lensing, where the magnification factor varies from $\sim$ 3 to $\sim$ 9 across different kinematic components. The intrinsically symmetric line profile and a smooth source-plane velocity gradient suggest that the host galaxy is an extended rotating disk, with a CO size of $R_{\rm CO}$$\sim$ 6 kpc and a dynamical mass of $M_{\rm dyn}$ $\sim$ 8 $\times 10^{10}$ $M_{\odot}$ . We also find a secondary CO-emitting source near RXS J1131 $-$ 1231 whose location is consistent with the optically-faint companion reported in previous studies. The lensing-corrected molecular gas masses are $M_{\rm gas}$ = (1.4 $\pm$ 0.3) $\times 10^{10}$ $M_{\odot}$ and (2.0 $\pm$ 0.1) $\times 10^{9}$ $M_{\odot}$ for RXS J1131 $-$ 1231 and the companion, respectively. We find a lensing-corrected stellar mass of $M_{*}$ = (3 $\pm$ 1) $\times 10^{10}$ $M_{\odot}$ and a star formation rate of SFRFIR = (120 $\pm$ 63) $M_{\odot}$ yr*-1*, corresponding to a specific SFR and star formation efficiency comparable to $z$ $\sim$ 1 disk galaxies not hosting quasars. The implied gas mass fraction of $\sim$ 18 $\pm$ 4% is consistent with the previously-observed cosmic decline since $z$ $\sim$ 2. We thus find no evidence for quenching of star formation in RXS J1131 $-$ 1231. This agrees with our finding of an elevated $M_{\rm BH}$$/$$M_{\rm bulge}$ ratio of $>$ 0.27 ${}^{+0.11}_{-0.08}$ % compared to the local value, suggesting that the bulk of its black hole mass is largely in place while its stellar bulge is still assembling.

Subject headings:

infrared: galaxies – galaxies: high-redshift – galaxies: ISM – galaxies: evolution – quasars: individual (RXS J1131 $-$ 1231) – radio lines: ISM

††slugcomment: Accepted to the ApJ

1. Introduction

Many recent studies of galaxy evolution have been focused on investigating the interplay between star formation and active galactic nucleus (AGN) activity across cosmic epochs (e.g., Di Matteo et al., 2005; Alexander et al., 2005; Hopkins et al., 2006; Coppin et al., 2008; Page et al., 2012; Simpson et al., 2012; Lamastra et al., 2013). It is currently not well-understood when and how the supermassive black holes (SMBHs) and stellar populations of present-day massive galaxies were assembled, but it is clear that the co-moving star formation rate and the black hole accretion rate densities both increased substantially since $z$ $>\,3$ and reached their climax at $z$ $\sim$ $2$ , followed by a rapid decline toward $z$ $\sim$ [math] (e.g., Hopkins & Beacom, 2006; Madau & Dickinson, 2014). A leading explanation for this decline is the decrease in molecular gas content and star formation efficiency (e.g., Erb et al., 2006; Carilli & Walter, 2013; Walter et al., 2014), but direct molecular gas measurements at intermediate redshift $(0.2\,<\,z\,<\,1)$ that could confirm this explanation remain largely limited to spatially unresolved CO observations of a modest sample of $\sim$$30$ ultra-luminous infrared galaxies (ULIRGs; Combes et al., 2011, 2013).

Meanwhile, empirical scaling relations such as the $M_{\rm BH}$$-$$M_{\rm bulge}$ relation (e.g., Magorrian et al., 1998; Häring & Rix, 2004) have been established locally, suggesting a co-eval growth between local SMBHs and their host galaxies. Attempts to extend this relation out to higher redshifts, beyond the peak epoch of star formation and AGN activity, have been made in recent years. These studies find that high- $z$ AGN host galaxies do not appear to follow the same $M_{\rm BH}$$-$$M_{\rm bulge}$ relation as nearby spheroidal galaxies. (e.g., Walter et al., 2004; Borys et al., 2005; McLure et al., 2006; Peng et al., 2006; Riechers et al., 2008; Coppin et al., 2008; Alexander et al., 2008). Yet, the $M_{\rm BH}$$-$$M_{\rm bulge}$ relation remains poorly-constrained at intermediate redshifts due to the difficulty in separating the stellar component contributing to the optical emission from that of the bright AGN. This stems from the limited resolving power, which restricts the dynamic range that can be achieved at positions near the AGN. Strong gravitational lensing provides the magnification necessary to spatially separate the AGN emission from the host galaxy stellar emission, significantly improving the dynamic range that can be achieved in studies of their host galaxies with current instruments.

The quasar RXS J113151.62 $-$ 123158 (hereafter RXJ1131) at $z_{\textrm{s, QSO}}$ = $0.658$ (Sluse et al., 2003, hereafter S03) is a particularly well-suited source for studying the evolution of molecular gas properties in quasar host galaxies and the connection between SMBHs and their host galaxies at intermediate redshift given its unique lensing configuration. The stellar emission in the host galaxy of RXJ1131 is lensed into an Einstein ring of $1$$\farcs$$83$ in radius that is clearly separated from the quadruply imaged quasar emission (Claeskens et al., 2006, hereafter C06). The foreground lens is an elliptical galaxy at $z_{\textrm{L}}$ = $0.295$ (S03). Reis et al. (2014) report a high spin parameter of $a$ $\sim$ $0.9$ for the moderate-mass black hole in RXJ1131 ( $M_{\rm BH}$ = $8$$\times 10^{7}$$M_{\odot}$ ; Sluse et al., 2012), with an intrinsic bolometric luminosity of $L_{\rm bol,X}$ = $1.3$$\times 10^{45}$ ergs s*-1* (Pooley et al., 2007).

In this paper, we present CO( $J$ = 2 $\rightarrow$ 1) and CO( $J$ = 3 $\rightarrow$ 2) line observations and broadband photometry spanning rest-frame UV to radio wavelengths towards RXJ1131 to study the properties of its molecular gas, dust and stellar populations. In §2, we outline details of the observations and of our data reduction procedures. In §3, we report results for the CO( $J$ = 2 $\rightarrow$ 1) and CO( $J$ = 3 $\rightarrow$ 2) emission and broadband photometry in RXJ1131. In §4, we present lens modeling and dynamical modeling of the CO( $J$ = 2 $\rightarrow$ 1) data and spectral energy distribution (SED) modeling of the photometric data. In §5, we discuss the ISM properties of the host galaxy of RXJ1131 and compare them to other galaxy populations at low and high redshift. Finally, we summarize the main results and present our conclusions in §6. We use a concordance $\Lambda$ CDM cosmology throughout this paper, with parameters from the WMAP9 results: $H_{0}$ = $69.32$ km s*-1* Mpc*-1*, $\Omega_{\rm M}$ = $0.29$ , and $\Omega_{\Lambda}$ = $0.71$ (Hinshaw et al., 2013).

2. Observations

2.1. CARMA CO( $J$ = 3 $\rightarrow$ 2)

Observations of the CO( $J$ = 3 $\rightarrow$ 2) rotational line ( $\nu_{\rm rest}$ = 345.79599 GHz) towards RXJ1131 at $z_{\rm s,QSO}$ = 0.658 were carried out with the Combined Array for Research in Millimeter-wave Astronomy (CARMA; Program ID: cf0098; PI: D. Riechers) in the D array configuration. The line frequency is redshifted to $\nu_{\rm obs}$ = 209.10443 GHz at the quasar redshift. Observations were carried out on 2014 February 02 under poor 1.5 mm weather conditions and on 2014 February 17 under good 1.5 mm weather conditions. This resulted in a total on-source time of 2.94 hours after flagging poor visibility data. The correlator setup provides a bandwidth of 3.75 GHz in each sideband and a spectral resolution of 12.5 MHz ( $\sim$ 17.9 km s*-1*). The line was placed in the lower sideband with the local oscillator tuned to $\nu_{\rm LO}\sim$ 216 GHz. The radio quasars J1127 $-$ 189 (first track) and 3C273 (second track) were observed every 15 minutes for pointing, amplitude, and phase calibration. Mars was observed as the primary absolute flux calibrator and 3C279 was observed as the bandpass calibrator for both tracks.

Given that the phase calibrator used for the first track was faint and was observed under poor weather conditions and that the phase calibrator used for the second track was far from our target source, the phase calibration is subpar, with an rms scatter $\sim$ 50 $\arcdeg$ over a baseline length of $\sim$ 135 m. We thus conservatively estimate a calibration accuracy of $\sim$ 40% based on the flux scale uncertainties, the gain variations over time, and the phase scatter on the calibrated data. We therefore treat the derived CO( $J$ = 3 $\rightarrow$ 2) line intensity with caution and ensure that our physical interpretation of this system and the conclusion of this paper do not rely on this quantity.

The miriad package was used to calibrate the visibility data. The calibrated visibility data were imaged and deconvolved using the CLEAN algorithm with “natural” weighting. This yields a synthesized clean beam size of 3 $\farcs$ 2 $\times$ 1 $\farcs$ 9 at a position angle (PA) of 8°for the lower sideband image cube. The final rms noise is $\sigma$ = 13.3 mJy beam*-1* over a channel width of 25 MHz. An rms noise of $\sigma$ = 0.83 mJy beam*-1* is reached by averaging over the line-free channels in both sidebands.

2.2. PdBI CO( $J$ = 2 $\rightarrow$ 1)

Observations of the CO( $J$ = 2 $\rightarrow$ 1) rotational line ( $\nu_{\rm rest}$ = $230.53800$ GHz) towards RXJ1131 were carried out using the IRAM Plateau de Bure Interferometer (PdBI; Program ID: S14BX; PI: D. Riechers). Based on the CARMA CO( $J$ = 3 $\rightarrow$ 2) line redshift of $z_{\rm CO(3-2)}$ = 0.655, the CO( $J$ = 2 $\rightarrow$ 1) line is redshifted to $\nu_{\rm obs}$ = 139.256 GHz. Two observing runs were carried out on 2014 December 06 and 2015 February 05 under good weather conditions in the C and D array configurations, respectively. This resulted in 3.75 hours of cumulative six antenna-equivalent on-source time after discarding unusable visibility data. The 2 mm receivers were used to cover the redshifted CO( $J$ = 2 $\rightarrow$ 1) line and the underlying continuum emission, employing a correlator setup that provides an effective bandwidth of 3.6 GHz (dual polarization) and a native spectral resolution of 1.95 MHz ( $\sim$ $4.2$ km s*-1*). The nearby quasars B1127 $-$ 145 and B1124 $-$ 186 were observed every 22 minutes for pointing, secondary amplitude, and phase calibration, and B1055 $+$ 018 was observed as the bandpass calibrator for both tracks. MWC349 and 3C279 were observed as primary flux calibrators for the C and D array observations, respectively, yielding calibration accuracy better than 15%.

The gildas package was used to calibrate and analyze the visibility data. The calibrated visibility data were imaged and deconvolved using the CLEAN algorithm with “natural” weighting. This yields a synthesized clean beam size of 4 $\farcs$ 44 $\times$ 1 $\farcs$ 95 (PA = 13°). The final rms noise is $\sigma$ = 1.45 mJy beam*-1* over 10 MHz (21.5 km s*-1*). The continuum image at $\nu_{\rm cont}\sim$ 139 GHz is produced by averaging over 3.16 GHz of line-free bandwidth. This yields an rms noise of 0.082 mJy beam*-1*.

2.3. Karl G. Jansky Very Large Array (Archival)

Our analysis also uses archival data of the 4.885 GHz radio continuum obtained with the Karl G. Jansky Very Large Array (VLA; Program ID: AW741; PI: Wucknitz). Observations were carried out on 2008 December 29 under excellent weather conditions in the A array configuration for a total of $\sim$ 7 hours on-source time. The C-band receivers were used with a continuum mode setup, providing a bandwidth of 50 MHz for the two IF bands with full polarization. The nearby radio quasar J1130 $-$ 149 was observed every 10 minutes for pointing, amplitude, and phase calibration. J1331 $+$ 305 was observed as the primary flux calibrator, and J0319 $+$ 415 was observed as the bandpass calibrator, yielding $\sim$ 10% calibration accuracy. We used aips to calibrate the visibility data. The calibrated visibility data were imaged and deconvolved using the CLEAN algorithm with robust = 0, which was chosen to obtain a higher resolution image given high SNR. This yields a synthesized clean beam size of 0 $\farcs$ 49 $\times$ 0 $\farcs$ 35 (PA = 0.18 $\arcdeg$ ) and a final rms noise of $\sigma$ = 13 $\mu$ Jy beam*-1*.

2.4. HST (Archival)

We obtained an HST image taken with the ACS using the F555W filter ( $V$ -band) from the Hubble Legacy Archive. The details of the observations can be found in C06. We apply an astrometric correction to the optical image by adopting the VLA 5 GHz map as the reference coordinate frame. We shift the latter to the east by 0 $\farcs$ 5963 in R.A. and $+$ 0 $\farcs$ 8372 in Dec., which is consistent with the typical astrometric precision (1 ${}^{\prime\prime}-$ 2*′′*) of images from the Hubble Legacy Archive111http://hla.stsci.edu/hla_faq. html. This astrometric correction is critical to avoid artificial spatial offsets between different emitting regions and to carry out our lens modeling, in which the absolute position of the foreground lensing galaxy is based on its coordinates in the high-resolution optical image. The VLA image is calibrated using a well-monitored phase calibrator, with absolute positional accuracy of $\sim$ 2 mas. For this reason, the absolute alignment between the VLA image and other interferometric images reported in this paper are expected to have an astrometric precision better than 0 $\farcs$ 1, modulo uncertainties related to the SNR and phase instability.

3. Results

3.1. CO( $J$ = 2 $\rightarrow$ 1) Emission

We detect CO( $J$ = 2 $\rightarrow$ 1) line emission towards the background source RXJ1131 in the PdBI data at $\gtrsim$ 27 $\sigma$ significance. Based on this measurement, we refine the redshift of RXJ1131 to $z_{\rm CO}$ = 0.6541 $\pm$ 0.0002222This redshift is derived by fitting a double-Gaussian to the de-lensed spectrum (Figure 5) instead of the observed spectrum (Figure 1) to avoid biases in our redshift determination due to differential lensing (see §4.1.2). . The emission is spatially and kinematically resolved with a highly asymmetric double-horned line profile as shown in Figure 1. Fitting a double Gaussian results in peak flux densities of 75.3 $\pm$ 2.6 mJy and 24.0 $\pm$ 2.0 mJy, and a FWHM of 179 $\pm$ 9 km s*-1* and 255 $\pm$ 28 km s*-1* for the two components, respectively. The peaks are separated by $\Delta v_{\rm sep}$ = 400 $\pm$ 12 km s*-1*. The total integrated line flux is 20.3 $\pm$ 0.6 Jy km s*-1*.

We construct a zeroth order moment map, red/blue channel maps, and first and second moment maps, as shown in Figure 2, using the $uv$ -continuum subtracted data cube over a velocity range of $\Delta v$ $\sim$ 750 km s*-1*. The higher-order moment maps are produced using unbinned channel maps with 3 $\sigma$ clipping. The peak flux density is 8.12 $\pm$ 0.30 Jy km s*-1* beam*-1* in the intensity map. Observed properties of the CO( $J$ = 2 $\rightarrow$ 1) emission line are summarized in Table 1.

The deconvolved source size FWHM obtained from fitting a single two-dimensional Gaussian to the integrated line emission in the image plane yields 5 $\farcs$ 1 $\pm$ 0 $\farcs$ 7 $\times$ 3 $\farcs$ 7 $\pm$ 0 $\farcs$ 7, which is consistent with that obtained by visibility-plane fitting within the uncertainties. Since the spatial distribution of the observed CO emission is unlikely to be fully described by a simple Gaussian and appears to be a superposition of at least two components (top left panel of Figure 2), we also fit two Gaussians to the intensity map. This yields deconvolved source sizes of 3 $\farcs$ 8 $\pm$ 0 $\farcs$ 4 $\times$ 1 $\farcs$ 9 $\pm$ 0 $\farcs$ 4 and 3 $\farcs$ 6 $\pm$ 0 $\farcs$ 3 $\times$ 1 $\farcs$ 5 $\pm$ 0 $\farcs$ 3, separated by $\sim$ 2 $\farcs$ 2 in RA and $\sim$ 1 $\farcs$ 7 in Dec. The deconvolved source sizes of both models suggest that the gravitationally lensed CO emission is more extended than the optical “Einstein ring”, which has a diameter of $\sim$ 3 $\farcs$ 6 (i.e., the “Einstein ring” formed by CO emission is likely to have a larger diameter compared to the optical one). This is consistent with the centroid position of the redshifted emission, which is along the quasar arc seen in the optical image, and the blueshifted emission, which is offset further to the SE (top right of Figure 2). Therefore, the CO-emitting region in RXJ1131 is likely to be more extended than its stellar and quasar emission.

We also place an upper limit on HNC( $J$ = 2 $\rightarrow$ 1) line emission in the foreground galaxy at $z\sim$ 0.295. Assuming a typical line width of 300 km s*-1*, this corresponds to a 3 $\sigma$ limit of 0.35 Jy km s*-1* beam*-1*.

3.2. CO( $J$ = 3 $\rightarrow$ 2) Emission

We detect CO( $J$ = 3 $\rightarrow$ 2) line emission towards RXJ1131 in the CARMA data at $\gtrsim$ 5 $\sigma$ significance. The CO( $J$ = 3 $\rightarrow$ 2) spectrum appears to be consistent with a double-peaked profile, as shown in Figure 3, where we over-plot spectra of the CO( $J$ = 2 $\rightarrow$ 1) and CO( $J$ = 3 $\rightarrow$ 2) lines. We extract the peak fluxes and their corresponding uncertainties for the blue and red wing independently. We find a peak line flux of 5.13 $\pm$ 1.43 Jy km s*-1* beam*-1* for the blue wing, indicating a $\gtrsim$ 3 $\sigma$ detection for this component alone, and a peak line flux of 11.45 $\pm$ 1.99 Jy km s*-1* beam*-1* for the red wing, indicating a $\sim$ 6 $\sigma$ detection. We measure a line intensity of 35.7 $\pm$ 6.9 Jy km s*-1* (Table 1) by summing up fluxes over the FWZI linewidth used to infer the CO( $J$ = 2 $\rightarrow$ 1) line intensity.

Assuming that the spatial extent of CO( $J$ = 2 $\rightarrow$ 1) and CO( $J$ = 3 $\rightarrow$ 2) is similar and therefore the emission is magnified by the same amount, the measured line intensities correspond to a brightness temperature ratio of $r_{\rm 32}$ = $T_{\small\mbox{CO($ J $\,=\,3\,$ \rightarrow $\,2)}}$$/$$T_{\small\mbox{CO($ J $\,=\,2\,$ \rightarrow $\,1)}}$ = 0.78 $\pm$ 0.37. The quoted error bar is derived by adding the uncertainties associated with the CO line intensities and those from absolute flux calibrations in quadrature. This brightness temperature ratio is consistent with thermalized excitation within the uncertainties, as commonly observed in nuclear regions of nearby ULIRGs and high- $z$ quasars (e.g., Weiß et al., 2007; Riechers et al., 2011b; Carilli & Walter, 2013), but also with the lower excitation seen in normal star-forming disks (e.g., Dannerbauer et al., 2009; Carilli & Walter, 2013; Daddi et al., 2015).

3.3. Continuum Emission

No 1.4 mm continuum emission is detected at the position of CO( $J$ = 3 $\rightarrow$ 2) down to a 3 $\sigma$ limit of 2.49 mJy beam*-1*. This is consistent with the spectrum shown in Figure 3.

We detect 2.2 mm continuum emission at an integrated flux density of 1.2 $\pm$ 0.2 mJy, with a peak flux of $S_{\nu}$ = 799 $\pm$ 88 $\mu$ Jy beam*-1* centered on the lensing galaxy (Figure 4). Slightly extended emission along the lensing arc is also detected. This suggests that we detect emission in both the foreground and the background galaxy and that the emission is marginally resolved along its major axis. We subtract a point source model in the visibility-plane to remove the unresolved part of the emission, which we here assume to be dominated by the foreground galaxy. The emission in the residual map coincides spatially with the lensing arc. We measure a flux density of $S_{\nu}$ = 0.39 $\pm$ 0.08 mJy for this residual component. This flux density is consistent with the difference between the integrated and the peak flux density measured in the original continuum map ( $\sim$ 0.4 mJy). We therefore adopt $S_{\nu}$ = 0.39 $\pm$ 0.12 mJy as the best estimate for the 2 mm continuum flux of the background galaxy (RXJ1131). We here quote a conservative error bar, which is derived by adding the uncertainty associated with the flux density of the point-source model ( $\delta S_{\nu}$ = 0.088 mJy) with that of the peak flux in the residual map (0.08 mJy) in quadrature. We caution that this does not account for the systematic uncertainties of the de-blending procedure, where we have assigned 100% of the point source flux to the foreground galaxy. We report the peak flux in the original map ( $S_{\nu}$ = 799 $\pm$ 88 $\mu$ Jy beam*-1*) for the foreground galaxy, which is the best estimate possible at the resolution of our observations, but we acknowledge that a non-negligible contribution from the background source to the peak flux cannot be ruled out.

The VLA C-band continuum image in Figure 4 shows resolved emission from the jets and the core of the foreground elliptical galaxy as well as emission toward the background quasar. Multiple peaks are seen along the arc with their centroids coincident with the optical emission from the quasar. We extract the flux densities for the lensing arc and the radio core in Table 3.4. We find a spectral index of $\alpha^{\rm 2mm}_{\rm 6cm}$ = $-$ 0.02 $\pm$ 0.07 for the foreground galaxy and $\alpha^{\rm 2mm}_{\rm 6cm}$ = $-$ 0.35 $\pm$ 0.21 for the background galaxy by fitting power-laws ( $S_{\nu}\propto\nu^{\alpha}$ ) to their continuum fluxes at 5 GHz and 2 mm. The spectral slope derived for the background source is flatter than the typical slope of pure synchrotron emission ( $\alpha$ $\sim$ $-$ 0.7; e.g., Andreani et al., 1993). This likely suggests that at least a fraction of the observed 2 mm emission arises from thermal dust emission. This spectral slope would be even shallower if the background source contributes to the unresolved fraction of the 2 mm flux. In this case, the 2 mm flux of the foreground galaxy would be lower than the value reported here and lead to a slope steeper than $\alpha^{\rm 2mm}_{\rm 6cm}$ $\sim$ $-$ 0.02, which is flatter than that typical of elliptical galaxies. Assuming a spectral slope of $\alpha$ $\sim$ $-$ 0.7 to account for synchrotron radiation in RXJ1131, we expect a flux density of $S_{\rm 2mm}$ = 0.122 $\pm$ 0.004 mJy at 2 mm. The flux excess of $S_{\rm 2mm}$ = 0.27 $\pm$ 0.08 mJy therefore likely arises due to thermal dust emission.

3.4. Photometry

We compile mid-IR (MIR) to far-IR broadband photometry from various catalogs available on the NASA/IPAC Infrared Science Archive (IRSA) in Table 3.4 with aperture corrections when warranted. These data were obtained using the Cerro Tololo Inter-American Observatory (CTIO) for the Two Micron All Sky Survey (2MASS; Skrutskie et al., 2006), the Wide-field Infrared Survey Explorer (WISE; Wright et al., 2010), the Infrared Astronomical Satellite (IRAS; Neugebauer et al., 1984), and the Multiband Imaging Photometer (MIPS; Rieke et al., 2004) and Mid-infrared Infrared Array Camera (IRAC; Fazio et al., 2004) on the Spitzer Space Telescope. We retrieve PBCD (level 2) Spitzer/IRAC images from the Spitzer Heritage Archive and perform aperture photometry on the channel 1 image to extract the flux density at 3.6 $\mu$ m since it is not available from the IRSA archive.

The emission in the IRAC images is slightly extended. We thus use an HST image ( $\sim$ 0 $\farcs$ 07 resolution) to determine the origin of their centroids, all of which are found to be centered at the position corresponding to the lensed emission from the background galaxy. To recover the diffuse background emission, we subtract a point source model centered on the lensing galaxy, using the average FWHM found by fitting a Gaussian profile to several field stars with the imexam routine of IRAF. We perform aperture photometry on the residual image to obtain decomposed flux measurements of the background galaxy. The photometry for the foreground galaxy is then obtained by subtracting the background emission from the observed total flux. The resulting photometry in Table 3.4 is obtained after performing an aperture correction described in the IRAC Instrument Handbook333http://irsa.ipac.caltech.edu/data/SPITZER/docs/irac/iracinstrumenthandbook/ to correct for the fact that the imaging was calibrated using a 12*′′* aperture, which is larger than the aperture (5 $\farcs$ 8) we used to perform aperture photometry.

We fit a power-law spectrum to the decomposed IRAC photometry to disentangle the background and foreground emission from the total flux observed in the MIPS 24 $\micron$ band. The spectral indices corresponding to the best-fitting curves are $\alpha$ = $-1.8$ and $\alpha$ = $-0.85$ for the lensing galaxy and RXJ1131, respectively. The latter is consistent with the mean 3.6 $-$ 8 $\micron$ spectral slope of $\alpha$ = $-$ 1.07 $\pm$ 0.53 found for unobscured AGN (Stern et al., 2005). An extrapolation of the fit to 24 $\micron$ yields 33.96 $\pm$ 0.01 mJy and 25.19 $\pm$ 0.03 mJy for the foreground galaxy and RXJ1131, respectively. The uncertainties are the standard deviations of the extrapolated fluxes obtained from two independent Monte Carlo simulations, each of 500 iterations. We incorporate the decomposed 24 $\micron$ data in our SED fitting to provide some constraints on the Wien tail beyond the dust peak of the SED of RXJ1131. Details of the SED modeling are presented in §4.5.

Extraction of the Herschel/SPIRE photometry at 250, 350, and 500 $\micron$ was carried out using sussextractor within the Herschel Interactive Processing Environment (HIPE; Ott, 2010) on Level 2 maps obtained from the Herschel Science Archive. These maps were processed by the SPIRE pipeline version 13.0 within HIPE. The sussextractor task estimates the flux density from an image convolved with a kernel derived from the SPIRE beam. The flux densities measured by sussextractor were confirmed by using the Timeline Fitter, which performs photometry by fitting a 2D elliptical Gaussian to the Level 1 data at the source position given by the output of sussextractor. The fluxes obtained from both methods are consistent within the uncertainties.

4. Analysis

4.1. Lens Modeling

At the angular resolution of the CO( $J$ = 2 $\rightarrow$ 1) data, the source is resolved $\gtrsim$ 2 resolution elements. Given the extent of the lensed emission (see Figure 2), this implies that we do not resolve structures (e.g. knots and arcs) of the lensed emission in our CO( $J$ = 2 $\rightarrow$ 1) data. Nevertheless, the high spectral resolution of these data provides kinematic information on spatial scales smaller than the beam (see Figure 2). Hence, we reconstruct the intrinsic line profile and source-plane velocity structure by carrying out a parametric lens modeling over different channel slices of the interferometric data using our lensing code uvmcmcfit (Bussmann et al. 2015a; see Bussmann et al. 2015b for details of the code). Our approach follows a similar strategy as Riechers et al. (2008), who reconstruct a source-plane velocity gradient and constrain the gas dynamics in the $z$ $>$ 4 quasar host galaxy of PSS J2322 $+$ 1944, which is also lensed into an Einstein ring configuration. To ensure adequate SNRs for lens modeling, we bin the frequency channels by a factor of five to produce seven independent $\Delta v$ $\sim$ 105 km s*-1* channels (dashed line in Figure 5) that cover the full linewidth of $\sim$ 750 km s*-1*.

We model the lens mass distribution using a singular isothermal ellipsoid (SIE) profile, which is described by five free parameters: the positional offset in R.A. and Dec. relative to an arbitrary chosen fixed coordinate in the image, the Einstein radius, the axial ratio, and the position angle. Positional offset between the foreground galaxy and the pre-defined coordinate is initialized using the VLA radio continuum map. We impose a uniform prior of $\pm$ 0 $\farcs$ 05 in both $\Delta$ R.A. and $\Delta$ Dec., motivated by the astrometry uncertainties in the VLA image as well as the uncertainties provided by previous SIE lens model (C06). We initialize the Einstein radius based on the model parameters reported by C06 and impose a uniform prior using $\pm$ 3 $\sigma$ of their uncertainties. The sources are modeled using elliptical Gaussian profiles, which are parameterized by six free parameters: the positional offset in R.A. and Dec. relative to the lens, the intrinsic flux density, the effective radius, the axial ratio, and the position angle. The position of each source is allowed to vary between $\pm$ 1 $\farcs$ 5 (i.e., within the Einstein radius) and the effective radius is allowed to vary from 0 $\farcs$ 01 $-$ 2*′′*.

Our code uses an Markov Chain Monte Carlo (MCMC) approach to sample the posterior probability distribution function (PDF) of the model parameters. In each model, we require a target acceptance rate of $\sim$ 0.25 $-$ 0.5 and check for chain convergence by inspecting trace plots and by requiring that the samples are obtained beyond at least an autocorrelation time. We thus employ $\sim$ 50,000 samples as the initial “burn-in” phase to stabilize the Markov chains (which we then discard) and use the final $\sim$ 5,000 steps, sampled by 128 walkers, to identify the posterior. Here, we identify the best-fit model and the quoted uncertainties using the median and the 68% confidence intervals in the marginal PDFs.

We first obtain a preliminary lens model for each channel slice independently, where their lens parameters are allowed to vary and are initialized according to the aforementioned way. We obtain the final model by repeating the modeling over each slice but fixing their lens parameters to the overall median in the preliminary models, as listed in Table 3. This ensures that all models share the same lens profile. The magnification factors in Table 4 are determined by taking the ratio between the image plane flux and the source plane flux of each model.

Our model parameters in Table 3, describing the mass distribution of the lensing galaxy, are consistent (within the uncertainties) with that of the SIE model presented by C06. We find a mass of $M(\theta$ $<$ $\theta_{\textrm{E}})$ = (7.47 $\pm$ 0.02) $\times$ 1011 $M_{\odot}$ within the Einstein radius.

4.1.1 Interpretation of the Source-plane Morphology

The reconstructed source locations, as represented by magenta ellipses in Figure 6, demonstrate an intrinsic velocity gradient across the source plane, which is consistent with a kinematically-ordered disk-like galaxy. Additional support to the disk conjecture can be found in the double-horned line profile (Figure 1) and the observed (image plane) velocity field (Figure 2). Furthermore, C06 also find that the reconstructed source plane emission in optical-NIR is best-reproduced using a $n$ = 1 Sersic profile. We thus interpret RXJ1131 as a disk galaxy.

A better fit is found for the lens model of the red-most channel if we add a second source component (see top left panel in Figure 6). This is consistent with previous results reported by Brewer & Lewis (2008, hereafter B08), who find an optically faint companion (component F in their paper) $\sim$ 2.4 kpc in projection from the AGN host galaxy in $V$ -band, and with C06, who find evidence for an interacting galaxy near RXJ1131. Spatially, the red velocity component of the CO emission also consistent with this component F. It is therefore likely that we detect CO( $J$ = 2 $\rightarrow$ 1) emission in a companion galaxy.

We decompose the total line flux into two components: one from RXJ1131 and the other from its companion. Since the companion is only detected in the red-most channel, we derive its intrinsic gas mass using the best-fit flux densities and magnification factors obtained from the models of this channel. Assuming a brightness temperature ratio of $r_{\rm 21}$ = 1 between CO( $J$ = 2 $\rightarrow$ 1) and CO( $J$ = 1 $\rightarrow$ 0) lines and a CO luminosity-to-H2 mass conversion factor of $\alpha_{\rm CO}$ = 0.8 $M_{\odot}~{}($ K km pc ${}^{2})^{-1}$ , we find a molecular gas mass of $M_{\rm gas}$ = $($ 1.92 $\pm$ 0.09 $)$ $\times 10^{9}$ $M_{\odot}$ . For the molecular gas mass in RXJ1131, we derive its intrinsic line flux over the FWZI linewidth using the respective magnification factors listed in Table 4, which to first order takes into account the effect of differential lensing. This yields $I_{\small\mbox{CO($ J $\,=\,2\,$ \rightarrow $\,1)}}$ = 2.93 $\pm$ 0.70 Jy km s*-1*, where the uncertainty includes those on the magnification factors. Adopting the same brightness temperature ratio and $\alpha_{\rm CO}$ as used for the companion, this corresponds to a gas mass of $M_{\rm gas}$ = $($ 1.38 $\pm$ 0.33 $)$ $\times 10^{10}$ $M_{\odot}$ , which implies a gas mass ratio of $\sim$ 7:1 between RXJ1131 and its companion.

The spatial resolution of the data in hand is a few arcsec, which implies that despite the high SNR and spectral resolution, constraints on the intrinsic sizes of the lensed galaxies are modest, and thus the magnification factors may be under-predicted (see e.g., Bussmann et al., 2015b; Dye et al., 2015; Rybak et al., 2015).

4.1.2 Spatial Extent and Differential Lensing

In the image-plane integrated line map shown in Figure 2, the redshifted component is cospatial with the Einstein ring that is seen in the optical image, with most of its apparent flux originating along the lensing arc, whereas centroid of the blueshifted emission is offset to the SE of the lensing arc. This suggests that the CO-emitting region in RXJ1131 is extended. To further illustrate this, we show the channel maps of 21.5 km s*-1* width and a spatial spectra map of 1 $\farcs$ 5 resolution in Figure 7 and Figure 8, respectively. These figures show that redshifted emission is present to the west, peaking toward the lensing arc (black crosses in Figure 7), and shifts to the east with decreasing velocity (blue wing). This is consistent with the source plane positions in our models and is suggestive of an extended CO emitting region.

Previous studies of RXJ1131 find evidence for differential lensing across the HST $V$ -, $I$ -, and $H$ -bands, where the magnification factor varies from 10.9 to 7.8 (C06). This indicates that the emission from different stellar populations within the host galaxy have various spatial extents and positions with respect to the caustic. The best-fit lens models obtained here for different CO channels show that differential lensing also plays an important role in the observed CO( $J$ = 2 $\rightarrow$ 1) emission, with a magnification factor $(\mu_{\rm L})$ that varies from 8.7 to 3.1 across different kinematic components (Table 4). The asymmetry in the line profile (Figure 1 and Figure 5) is therefore predominantly a result of the redshifted CO-emitting gas being more strongly-magnified than the blueshifted component. A secondary reason is likely due to the inclusion of the emission of the companion in the most redshifted velocity channels. The variation in $\mu_{\rm L}$ found across channels is consistent with the source plane positions relative to the caustics in Figure 6, where the red wing emission mainly originates near the cusp of the caustic and the blue wing emission is located beyond the caustics. In fact, the intrinsic line flux of the redshifted and blueshifted emission in RXJ1131 (after subtracting a contribution from the companion) are $I_{\small\mbox{CO($ J $\,=\,2\,$ \rightarrow $\,1)}}$ = 1.26 $\pm$ 0.23 Jy km s*-1* and 1.25 $\pm$ 0.23 Jy km s*-1*, respectively, implying an intrinsically symmetric line profile (Figure 5). This is consistent with the symmetric source-plane velocity gradient in our lens model (Figure 6 and Figure 9).

4.2. CO( $J$ = 2 $\rightarrow$ 1) Kinematics

Fitting two Gaussians with a common FWHM to the “intrinsic” CO( $J$ = 2 $\rightarrow$ 1) line profile of RXJ1131 (after correcting for lensing using the magnification factors for various channels and separating the emission from RXJ1131 and its companion), we find a roughly symmetric double-horned profile with a flux ratio of 1.2 $\pm$ 0.4 between the peaks, which are separated by $\Delta v_{\rm sep}$ = 387 $\pm$ 45 km s*-1*, and each with a FWHM of 220 $\pm$ 72 km s*-1*. The peak separation obtained from this “intrinsic” line profile is slightly lower than that obtained from the observed spectrum (i.e., without lensing corrections). This discrepancy is likely a result of differential lensing, which causes the line peak of the red wing to shift towards higher velocity channels, and thereby biasing the centroid of one Gaussian to higher velocity than otherwise. To facilitate a comparison (§5.1.1) with previous works, which were observed at lower spectral resolution, we also fit a single-Gaussian to the intrinsic line profiles. This yields FWHMs of 600 $\pm$ 160 km s*-1* for RXJ1131 and 73 $\pm$ 43 km s*-1* for the companion galaxy.

A clear velocity gradient and a high velocity dispersion ( $\gtrsim$ 400 km s*-1*) near the central region are seen in Figure 2. While beam smearing is inevitably the dominant factor in the observed velocity dispersion at the spatial resolution of these data, the exceedingly high velocity dispersion may hint at potential perturbations from the AGN, or internal turbulence due to interactions with the companion, and/or instability due to the large gas content. Therefore, in this scenario, RXJ1131 is consistent with a disrupted disk galaxy hosting an optically bright quasar and in the process of merging.

4.3. CO( $J$ = 2 $\rightarrow$ 1) Dynamical Modeling

As discussed in §4.1.1, we interpret RXJ1131 as a disk galaxy as it displays a kinematically-ordered velocity gradient in the source-plane velocity map of the CO emission, a symmetric double-horned line profile (Figure 5, Figure 6 and, left panel of Figure 9), and a disk-like morphology in the source-plane reconstruction of the optical-NIR emission (C06). We extract a one dimensional position-velocity (PV) profile by assuming that the source-plane centroids of different velocity components obtained from dynamical lens modeling are dominated by the tangential component of the true velocity vector of a rotating disk, i.e., each velocity component would be seen as lying along the major axis of a rotating disk if observed with sufficiently high angular resolution (see right panel of Figure 9). In this process, the positions for each velocity component (plotted as data points in the right panel of Figure 9) are extracted along the best-fitted major axis, which is along a PA of 121°.

We then attempt to characterize the molecular gas kinematics using an empirically-motivated disk model (e.g., Courteau, 1997; Puech et al., 2008; Miller et al., 2011):

[TABLE]

where $V$ is the observed velocity, $V_{0}$ is the velocity at dynamical center, $V_{a}$ is the asymptotic velocity, and $R_{t}$ is the “turnover” radius at which the rising part of the curve begins to flatten. We perform non-linear least squares fitting using an orthogonal distance regression to find the best-fit parameters, taking into account the uncertainties in both velocity (channel width) and distance offset. We also place an upper limit on $R_{t}$ $<$ 15 kpc to keep this parameter physical (e.g., Puech et al., 2008; Miller et al., 2011). The parameter uncertainties are inferred based on a Monte Carlo simulation of 500 iterations, where the input parameters are perturbed according to random Gaussian distributions with standard deviations corresponding to their uncertainties. Using this model, we find $V_{a}$ = 988 $\pm$ 618 km s*-1*, $R_{t}$ = 10.9 $\pm$ 7.8 kpc, and $V_{0}$ = 0 $\pm$ 9 km s*-1*. However, since the emission is not resolved along the flat regime of the rotation curve, the asymptotic velocity and the “turnover” radius are poorly constrained. In particular, $V_{a}$ and $R_{t}$ are highly correlated with a Pearson coefficient $R$ = 0.998, and $-$ 0.400 between $V_{a}$ and $V_{0}$ .

The asymptotic velocity ( $V_{a}$ ) — an extrapolation of the model out to radius beyond the disk scale-length and half-light radius — is not equivalent to the maximum observed velocity ( $V_{\rm max}$ ), which is commonly used in literature to parameterize disk rotation. The arctangent model is most commonly used in studies of the Tully-Fisher relation, where an extrapolation to V2.2 (velocity at 2.2 disk scale-length or $\sim$ 1.375 half-light radius, or $\sim$ 0.7 $R_{\rm opt}$ 444Radius enclosing 83% of the light distribution.) is typically adopted as the rotation velocity ( $V_{\rm max}$ in their terminology), since this corresponds to the radius at which the velocity of a pure exponential disk peaks (Courteau & Rix, 1997). Here, we adopt the maximum observed velocity $V_{\rm rot}$ = 303 $\pm$ 55 km s*-1* at 6 $\pm$ 3 kpc from the dynamical center as a proxy to the rotation velocity. This radius corresponds to $\sim$ 0.6 $R_{e}$ , where $R_{e}$ is the half-light radius $\sim$ 10.3 kpc inferred from the HST $I$ -band lens model (C06; converted to our cosmology). We note that the source plane half-light radius varies substantially with wavelength. In particular, the half-light radius is found to be $\sim$ 4 kpc and $\sim$ 7 kpc in $V$ -band (B08) and $H$ -band (C06), respectively. The CO gas is thus of similar spatial extent as in the $H$ and $I$ -bands.

In the rest-frame, emission in the observed-frame $H$ -band corresponds to NIR emission $(\sim$ 1 $\micron)$ , tracing radiation from the accretion disk surrounding the central AGN and also from old and evolved stellar populations; $I$ -band corresponds to roughly the optical $V$ -band, tracing stellar radiation from existing, less massive (i.e., longer-lasting) stars; $V$ -band corresponds to roughly $U$ -band, tracing radiation from massive young stars in the host galaxy. Hence, the relative compactness observed in the $V$ -band may be explained in part due to the fact that the emission in this band is more susceptible to dust extinction than in other bands and/or dominated by a central starburst caused by higher concentrations of star-forming gas towards the central regions — owing to gravitational perturbations induced from interactions with the companion (e.g., Di Matteo et al., 2005). This would be consistent with the picture that old stars form first and constitute the bulge component of a spiral galaxy, and that nuclear starbursts (in the inner few kpc) can be triggered at a later time as the progenitor disk galaxy interacts with other galaxies to form a larger bulge.

4.4. Dynamical Mass

Assuming the gas to be virialized, the dynamical mass can be approximated by $M_{\rm dyn}$ $\sim$ $\sigma^{2}R/G$ , where $\sigma$ is the velocity dispersion, or the rotational velocity in the case of a rotating disk model $($ i.e., $\sigma$ = $V_{\rm rot}\,\sin\,i)$ . Using a rotational velocity $V_{\rm rot}\,\sin\,i$ = 303 km s*-1* (see §4.3), we find a dynamical mass of $M_{\rm dyn}$ $\sin^{2}i$ $(<$ 6 kpc $)$ = 1.3 $\times 10^{11}$ $M_{\odot}$ enclosed within the CO-emitting region in RXJ1131. If we instead consider the CO( $J$ = 2 $\rightarrow$ 1) line peak separation $(\Delta v_{\rm sep}/2\sim$ 200 km s*-1* $)$ as the rotation velocity, we find $M_{\rm dyn}$ $\sin^{2}i$ $(<$ 6 kpc $)$ = 5.8 $\times 10^{10}$ $M_{\odot}$ . We derive an inclination angle of 56.4 $\arcdeg$ from the morphological axial ratio of $a/b\sim$ 1 $\farcs$ 8 $/$ 3 $\farcs$ 25, which we estimate from the source-plane image reconstructed by C06 (Figure 3 in their paper). This corresponds to an inclination-corrected dynamical mass of 8.3 $\times 10^{10}$$M_{\odot}$ $<$ $M_{\rm dyn}$ $<$ 25 $\times 10^{10}$$M_{\odot}$ . Our estimate should be considered at best an upper limit since the gas in RXJ1131 is unlikely to be virialized. In the following sections, we use the lower limit (8.3 $\pm$ 1.9) $\times 10^{10}$ $M_{\odot}$ as the dynamical mass as it is derived in a manner similar to what is commonly used in literature (e.g., Solomon et al. 1997, hereafter S97; Downes & Solomon 1998, hereafter DS98; Greve et al. 2005, hereafter G05).

Using the velocity dispersion ( $\sigma$ = 30 km s*-1*) obtained by fitting a single Gaussian to the de-lensed line profile of the companion and a half-light radius of $R_{\rm CO}$ = 4.2 $\pm$ 2.8 kpc from the best-fit lens model, we find a dynamical mass of $M_{\rm dyn}$ = (3.5 $\pm$ 2.3) $\times 10^{9}$ $M_{\odot}$ for the companion, assuming an inclination angle of $i$ = 30 $\arcdeg$ . The uncertainty here only includes that of the CO source size. On the other hand, we find $M_{\rm dyn}$ $\sin^{2}30\arcdeg$ = 5.8 $\times 10^{8}$ $M_{\odot}$ if we adopt the better-constrained $V$ -band source size of $\sim$ 700 pc (B08). Since the $V$ -band based dynamical mass measurement is substantially lower than the gas mass, the $V$ -band emitting region may appear to be much smaller than its true extent due to dust obscuration.

The CO-based dynamical mass estimates correspond to a mass ratio of $\sim$ 24:1 between RXJ1131 and the companion, with a gas mass ratio of $\sim$ 7:1 derived in §4.1.1. We thus classify the system as a gas-rich,“wet” minor merger.

4.5. SED Modeling

We fit dust SED models to the 24 µm $-$ 2.2 mm photometry using a modified-blackbody (MBB) function with a power-law attached to the Wien side to account for the MIR excess due to emission of warm and small dust grains. The IRAS 60 $\micron$ and 100 $\micron$ upper limits are included to constrain the dust peak. Here, we use a flux density of $S_{\rm 2mm}$ = 0.27 $\pm$ 0.08 mJy derived in §3.3 instead of the deblended flux listed in Table 3.4 to exclude a potential contribution due to synchrotron emission (see §3.3) in the dust SED modeling. An uncertainty from absolute flux calibration of $\sim$ 15% is added in quadrature to the PdBI 2 mm continuum photometry in our fitting procedure.

The fit is performed using the code mbb_emcee (e.g., Riechers et al., 2013; Dowell et al., 2014), which samples the posterior distributions using an MCMC approach and uses instrumental response curves to perform color correction. The model is described by five free parameters: the rest-frame characteristic dust temperature ( $T_{d}$ ), the emissivity index ( $\beta$ ), the power-law index ( $\alpha$ ), the flux normalization at 500 $\micron$ ( $f_{\rm norm}$ ), and the observed-frame wavelength at which the emission becomes optically thick ( $\lambda_{0}$ ). We impose a uniform prior with an upper limit of 100 K on $T_{d}$ (see e.g., Sajina et al., 2012), a Gaussian prior centered around 1.9 with a standard deviation of 0.3 on $\beta$ , and a uniform prior with an upper limit of 1000 $\micron$ on $\lambda_{0}$ . We check for chain convergence by requiring that the autocorrelation length of each parameter is less than the number of steps taken for the burn-in phase (which are then discarded). Here we report the statistical means and the 68% confidence intervals in the marginal PDFs as the best-fit parameters, as listed in Table 5. The best-fit models are shown in Figure 10 along with the broadband photometry that is listed in Table 3.4.

In the first model, we attempt to constrain the power-law index by including the 24 $\micron$ data. Based on the resulting posterior PDFs, we find an apparent IR luminosity (rest-frame 8 $-$ 1000 $\micron$ ) of 8.22 ${}^{+2.75}_{-2.98}$$\times 10^{12}$ $L_{\odot}$ , a far-IR luminosity (rest-frame 42.5 $-$ 122.5 $\micron$ ) of 3.81 ${}^{+1.92}_{-1.97}$$\times 10^{12}$ $L_{\odot}$ , and a dust mass of 16 ${}^{+5}_{-12}$$\times 10^{8}$ $M_{\odot}$ , none of which are corrected for lensing magnification. For the mass absorption coefficient, we adopt $\kappa$ = 2.64 m2kg*-1* at rest-frame 125.0 $\micron$ (Dunne et al., 2003). The dust mass uncertainty does not include that of the absorption coefficient.

A fit including the MIR 24 $\micron$ photometry is likely an upper limit on the far-IR luminosity due solely to star formation in the AGN host galaxy. If we instead fit for a model excluding this constraint, two major consequences are immediately apparent. First, the power-law index is poorly-constrained (see Table 5). Second, the steep power-law implies only a small contribution from the power-law regime to the total IR luminosity as compared to the graybody component. Thus, the far-IR luminosity in this model should, in principle, correspond to a lower limit on the cold dust emission. Using the best-fit parameters for this model, we find a total IR luminosity $L_{\rm IR}$ (rest-frame 8 $-$ 1000 µm) of 8.67 ${}^{+5.27}_{-5.27}$$\times 10^{12}$ $L_{\odot}$ , a far-IR luminosity $L_{\rm FIR}$ of 4.24 ${}^{+2.17}_{-2.00}$$\times 10^{12}$ $L_{\odot}$ and a dust mass $M_{\rm dust}$ of 14 ${}^{+5}_{-7}$$\times 10^{8}$ $M_{\odot}$ , all of which are not lensing-corrected. Taken at face value, this implies an FIR-to-IR luminosity ratio of $\sim$ 49 $\pm$ 38%.

The dust temperature from both models is similar to that of ULIRGs at 0.6 $<$ $z$ $<$ 1.0 (54 $\pm$ 5 K; Combes et al. 2013, hereafter C13). The far-IR luminosity is comparable in both models, which is not surprising given the lack of constraints in the MIR. For the subsequent analysis, we adopt the physical quantities from the first model (i.e., with constraints at 24 µm). The choice of SED model does not affect the derived star formation rate (SFR) given the similar far-IR luminosity, and their dust masses are consistent within the uncertainties. We correct for lensing using the median magnification factor $(\mu_{\rm L}\,=\,5.5)$ from the CO lens models. This yields a $L_{\rm FIR}$ of $($ 6.9 $\pm$ 3.6 $)$$\times 10^{11}$ $L_{\odot}$ and an intrinsic total IR luminosity of $\sim$ 1.5 $\times 10^{12}$ $(5.5/\mu_{\rm L})$ $L_{\odot}$ , implying that RXJ1131 can be classified as an ULIRG. Assuming a Salpeter (1955) initial mass function (IMF), we find a SFRFIR of 120 $\pm$ 63 $M_{\odot}$ yr*-1* using a standard conversion (Kennicutt, 1998).

We derive the stellar mass of RXJ1131 by fitting SED models to the rest-frame UV-to-mm photometry using the high- $z$ version of the magphys code (da Cunha et al., 2008, 2015). Two sets of stellar templates modeled using either the Bruzual & Charlot (2003) or the unpublished Charlot & Bruzual 2007 stellar population synthesis code are provided in the magphys package. We adopt the former set. To minimize contaminations from the quasar, we only fit to the HST, Herschel, and PdBI data, where both the HST and the PdBI 2 mm photometry are de-blended from the AGN (see bottom section of Table 3.4). The input photometry are corrected for lensing using their respective magnification factors to account for differential lensing (light blue circles in Figure 10). We thus find a stellar mass of $M_{*}$ = 2.95 ${}^{+1.32}_{-0.86}$$\times 10^{10}$ $M_{\odot}$ , which is the median value of the posterior probability distribution and the uncertainties are derived from the 16th and 84th percentiles. We note that the models are over-fitted with a best-fit $\chi^{2}$ = 0.41 which is unsurprising due to sparse sampling of the SED compare to the number of free parameters. The resulting dust mass and IR luminosity are consistent with those obtained from the MBB $+$ power-law models within the uncertainties, albeit some differences in the assumptions behind the two methods. The consistency may be attributed to the large uncertainties arising from the lack of photometric constraints on the models and the fact that the best-fit parameters from the MBB method are similar to those of the magphys method.

5. Discussion

5.1. ISM Properties

In this section, we derive the gas properties of the merging system RXJ1131 based on CO( $J$ = 2 $\rightarrow$ 1) and compare them with those reported by C13555The far-IR luminosity in C13 is derived based on 60 $\micron$ and 100 $\micron$ IRAS fluxes, and using a different definition of $L_{\rm FIR}$ : rest-frame 40 $-$ 500 µm. Following this convention, we find a far-IR luminosity of $L_{\rm FIR}$ = $($ 8.8 $\pm$ 0.4 $)$$\times 10^{11}$$(\mu_{\rm L}$ /5.5 $)\mbox{$ {}^{-1} $}$ $L_{\odot}$ and a SFR of $($ 150 $\pm$ 70 $)$ $M_{\odot}$ yr*-1* for RXJ1131. — the largest sample of CO-detected ULIRGs at similar redshift (0.6 $<$ $z$ $<$ 1.0). Their results are based on spatially unresolved CO( $J$ = 2 $\rightarrow$ 1) and CO( $J$ = 4 $\rightarrow$ 3) line observations with the IRAM 30-m single-dish telescope.

5.1.1 Linewidths and Sizes

The FWHM linewidth of $\Delta v$ $\sim$ 600 $\pm$ 160 km s*-1* found for RXJ1131 by fitting a single Gaussian is considerably larger than the statistical average in the C13 sample (370 km s*-1*) and local ULIRGs (median: 300 $\pm$ 85 km s*-1*, with the largest being 480 km s*-1*; S97, ). Linewidths exceeding 500 km s*-1* are also commonly observed in high- $z$ starburst galaxies (e.g., G05, ) and high- $z$ quasar host galaxies (e.g., Coppin et al., 2008), which are believed to originate from mergers. The wider CO linewidth observed in RXJ1131 also supports a merger picture.

The CO gas in RXJ1131 is $\sim$ 6 $\pm$ 3 kpc in radius (in the source plane), which is more extended than the average of 3.5 $\pm$ 2.3 kpc in a sample of disk-like U/LIRGs studied by Ueda et al. (2014), but consistent with their range of 1.1 $-$ 9.3 kpc. Our CO size is also consistent with that of high- $z$ ( $z$ $>$ 1) galaxies ( $R\sim$ 4 $-$ 20 kpc; G05; Daddi et al., 2010; Riechers et al., 2011a; Ivison et al., 2011) and local U/LIRGs in the Gao & Solomon (1999) sample ( $R\lesssim$ 10 kpc).

5.1.2 Gas Mass Fractions and Gas-to-dust Ratio

We find a dynamical gas mass fraction of $f_{\rm gas}$ = $M_{\rm gas}$ / $M_{\rm dyn}$ = 18 $\pm$ 4% and a baryonic gas mass fraction of $f_{\textrm{mol}}$ = $M_{\rm gas}$ /( $M_{\rm gas}+M_{*}$ ) = 34 $\pm$ 16% for the merger system (i.e., RXJ1131 and companion). Recent studies find that the baryonic gas fraction of starburst galaxies has decreased from $f_{\rm mol}$ $\sim$ 40% to $\lesssim$ 10% between $z$ $\sim$ 2 and $z$ $\sim$ 0 (S97; Gao & Solomon, 2004; Tacconi et al., 2006), and from $f_{\textrm{mol}}$ $\sim$ 50% to $\sim$ 5% between the same redshift range for “normal star-forming” galaxies (Geach et al., 2011; Saintonge et al., 2011; Tacconi et al., 2013)666These authors use the “Galactic” value of $\alpha_{\rm CO}$ = 4.6 $M_{\odot}~{}($ K km pc ${}^{2})^{-1}$ to compute the molecular gas mass.. Both the dynamical and baryonic gas mass fractions of RXJ1131+companion are thus consistent with the trend of decreasing molecular gas content since $z$ $\sim$ 2 which has been suggested as the cause for the decline in sSFR and cosmic star formation history towards $z$ $\sim$ 0 (e.g., Tacconi et al., 2013; Carilli & Walter, 2013; Genzel et al., 2015).

Using the lensing-corrected dust mass, we find a galactic-scale gas-to-dust ratio (GDR) of 54 $\pm$ 13. This would be higher by a factor of two if we were to adopt a dust mass from the other SED fit that is unconstrained at 24 $\micron$ . This GDR is lower than the statistical average of 206 in the C13 sample but is well within the broad range of values measured over their entire sample ( $\sim$ 1 $-$ 770). Our ratio is also consistent with high- $z$ SMGs (Bothwell et al., 2013) and local ULIRGs (Wilson et al., 2008), but lower than that of the Milky Way by $\sim 7\sigma$ (ignoring systematic uncertainties; Li & Draine, 2001; Zubko et al., 2004; Draine et al., 2007).

There are a number of systematic uncertainties associated with the derived gas-to-dust ratio, in particular the mass opacity coefficient $\kappa$ , the $\alpha_{\rm CO}$ conversion factor, and the brightness temperature ratio $r_{\rm 21}$ . If we instead use the “Galactic” $\alpha_{\rm CO}$ value, which may be more appropriate for some ULIRGs (e.g., Papadopoulos et al., 2012) and minor mergers (Narayanan et al., 2012), the gas mass (and thus gas-to-dust ratio) would be $\sim$ 6 times higher. We note that this gas mass is physically possible based on the dynamical mass constraints derived in §4.4. On the other hand, we would also obtain a higher gas mass if we were to assume sub-thermal excitation between CO( $J$ = 2 $\rightarrow$ 1) and CO( $J$ = 1 $\rightarrow$ 0) emission. We also note that the gas-to-dust ratio derived for RXJ1131 may be biased low as the gas is likely to be more extended than the optically thick dust. Consequently, the overall magnification factor for the CO gas may be lower than the optically thick dust, which dominates the far-IR luminosity. This would lead to an overestimation of the dust mass by adopting the CO magnification factor for the dust.

5.1.3 Star Formation Efficiency and specific SFR

To first order, the star formation efficiency $($ SFE = $L_{\rm FIR}$ $/$$M_{\rm gas})$ indicates the star formation rate per unit solar mass of molecular gas available in a galaxy. Using a wavelength range of 40 $-$ 500 $\micron$ defined in C13 for the far-IR luminosity, we find an SFE of 58 $\pm$ 10 $L_{\odot}$ $M_{\odot}^{-1}$ , which is on the low end among other U/LIRGs at $z$ $<$ 0.6 (S97; Combes et al., 2011) but consistent with those of low- $z$ spiral galaxies ( $z$ $<$ 0.1; Solomon & Vanden Bout, 2005) and high- $z$ disk-like galaxies, which are also IR luminous galaxies with $L_{\rm IR}$ $\sim$ 1012 $M_{\odot}$ (Daddi et al., 2008, 2010). This suggests that the merger system is converting gas into stars at an efficiency similar to those of “normal” star-forming disk-like galaxies rather than starburst galaxies (Tacconi et al., 2008; Riechers et al., 2011a, C13). This is in agreement with its disk-like kinematic signatures and its extended molecular gas distribution. Assuming the star formation continues at the current rate without gas replenishment, the SFE corresponds to a gas depletion time of $\tau$ = 102 $\pm$ 25 Myr.

The specific star formation rate (sSFR = SFR/ $M_{*}$ ) of 4 ${}^{+2.6}_{-2.4}$ Gyr*-1* derived for RXJ1131 is $\lesssim$ 1.5 $\sigma$ above the main sequence according to the redshift-dependent “main sequence” relation in Tacconi et al. (2013, and references therein). Given that RXJ1131 shares similar star formation rate, star formation efficiency, and CO disk size as other “main sequence” disk galaxies, the small elevation in sSFR over the main sequence at $z$ $\sim$ 0.7 suggests that the star formation activity in RXJ1131 may be enhanced by interactions with the companion.

The host galaxy of RXJ1131 is an extended disk with low star formation efficiency in a minor merger system. Therefore, removal of angular momentum of the gas via gravitational torque is likely inefficient to convert the entire gas disk into a massive stellar bulge. In this case, the disk component may be retained upon merging with the companion. This scenario is consistent with the results from recent simulations, which suggest that bulge formation maybe suppressed in gas-rich mergers, thereby allowing the formation of large disk galaxies with low bulge-to-disk ratios (Springel & Hernquist, 2005; Robertson et al., 2006; Hopkins et al., 2009). This also supports the idea that not all mergers will transform into elliptical galaxies, as in the classical picture (Toomre & Toomre, 1972).

5.2. Systemic Redshift and Velocity Offset

Sluse et al. (2007) report two sets of AGN lines observed in RXJ1131. The first set of lines is at $z$ $\sim$ 0.654, including the narrow component of the Balmer lines, the $[$ OIII $]$ 4959, 5007Å lines, and the MgII 2798Å absorption line; the second set is at $z_{\rm s,QSO}$ $\sim$ 0.658, including the broad component of the Balmer lines and the MgII 2798Å emission line. Using the CO line center redshift as the systemic redshift, we find that the redshift of the first set is fully consistent with the systemic redshift. This supports previous claims that [OIII] lines, tracing the narrow line region (NLR), are good proxies to the true systemic redshift (e.g., Vrtilek, 1985; Nelson, 2000). On the other hand, the second set of lines is redshifted by $\sim$ 715 km s*-1*.

Velocity offsets between broad line region (BLR) and NLR lines have been reported in literature. Richards et al. (2002) find a median offset of $\sim$ 100 $\pm$ 270 km s*-1* between [MgII] and [OIII] lines in a sample of $>$ 3800 quasars, and Bonning et al. (2007) report a mean offset of $\sim$ 100 $\pm$ 210 km s*-1* between the broad component of H $\beta$ and [OIII] lines in a sample of $\sim$ 2600 quasars at 0.1 $<$$z$$<$ 0.8, where only $\lesssim$ 20 of them (i.e., $<$ 1%) are found to have offsets $>$ 800 km s*-1* and $\sim$ 1% are found to have offsets $>$ 500 km s*-1* 777Bonning et al. (2007) report the fraction of objects with offset velocities greater than 500, 800, 1000, 1500, 2000, and 2500 km s*-1*. We therefore quote the two fractions corresponding to offset velocities closest to that of RXJ1131 ( $\sim$ 715 km s*-1*) in this discussion. Thus, large velocity offsets between BLR and NLR lines comparable to that of RXJ1131 are uncommon but have been observed in some cases.

The observed velocity offset between the BLR and NLR lines may be explained by a recoiling black hole (BH), where the BLR is moving at high velocity relative to the bulk of its host galaxy (Madau & Quataert, 2004; Bonning et al., 2007; Loeb, 2007). Depending on the initial conditions of the black hole pair (e.g., black hole mass ratio, spin-orbit orientation, spin magnitude), numerical relativity simulations have shown that recoil velocities can reach up to $v_{\rm kick}$ $\sim$ 4000 km s*-1* for spinning BHs, with typical recoil velocities of $v_{\rm kick}$ $\sim$ 100 $-$ 500 km s*-1* (e.g., Libeskind et al., 2006; Campanelli et al., 2007). Several sources have been proposed as recoiling BH candidates (Komossa et al., 2008; Civano et al., 2010; Steinhardt et al., 2012). However, Decarli et al. (2014) have recently refuted such scenario for one of the candidates — SDSS J0927+2943 — by finding that the redshift of its BLR lines is indeed consistent with its CO systemic redshift. This is in contrast with RXJ1131, where our CO observations confirm the redshifted BLR lines compared to the CO systemic redshift. Since this scenario requires a coalesced BH, it would imply that RXJ1131 is a product of a previous merger, which is not implausible and might also explain the highly spinning BH in RX1131 ( $a$ $\sim$ 0.9; Reis et al., 2014).

Alternative scenarios e.g., outflow/inflow of gas in the BLR, viewing angle towards the accretion disk, and obscuration in the clumpy accretion disk are more commonly invoked to explain velocity offsets between BLR and systemic redshift. Since the BLR lines of RXJ1131 show positive velocity offsets with respect to its systemic redshift, it may imply that the observed BLR line emission is dominated by the gas that is flowing into the central BH, or by the receding component of the accretion disk, owing to the viewing angle or the obscuration in the accretion disk. Sluse et al. (2007) report a covering factor of 20% for the accretion disk in RXJ1131 based on its broad MgII 2798Å absorption line at $z$ = 0.654. Additionally, the centroids of the BLR lines in RXJ1131 may be biased towards longer wavelengths due to microlensing (e.g., Sluse et al., 2007, 2012), which may have magnified the redshifted component of the compact BLR more strongly than its blueshifted component.

5.3. The $M_{\rm BH}$$-$$M_{\rm bulge}$ Relation

We find a $M_{\rm BH}$$/$$M_{\rm bulge}$ ratio of $>$ 0.27 ${}^{+0.11}_{-0.08}$ % using the black hole mass of $M_{\rm BH}$ $\sim$ 8 $\times 10^{7}$ $M_{\odot}$ (Sluse et al., 2012) and the stellar mass derived in §4.5 as an upper limit to the bulge mass. This ratio is consistent with those of other intermediate- $z$ radio-loud AGNs (McLure et al., 2006) but is higher than those of nearby AGNs (Häring & Rix, 2004). Our results therefore support the emerging picture that quasars grow faster and/or earlier than their host galaxies at higher redshifts (e.g., Walter et al., 2004; Peng et al., 2006; McLure et al., 2006; Riechers et al., 2008). The elevated $M_{\rm BH}$$/$$M_{\rm bulge}$ ratio of RXJ1131 compared to local AGNs suggests that the bulk of the black hole mass of RXJ1131 is largely in place while its stellar bulge is still assembling.

6. Summary and Conclusions

We present PdBI CO( $J$ = 2 $\rightarrow$ 1) and CARMA CO( $J$ = 3 $\rightarrow$ 2) observations towards the quadruply-imaged quasar RXJ1131 at $z_{\rm CO}$ $\sim$ $0.654$ , making this the first resolved CO study at intermediate redshift. Using the CO line intensities, we find a brightness temperature ratio of $r_{32}$ = $0.78\,\pm\,0.37$ between the CO( $J$ = 2 $\rightarrow$ 1) and CO( $J$ = 3 $\rightarrow$ 2) lines, consistent with thermalized excitation but also with the lower excitation seen in normal star-forming disks. We also detect marginally resolved 2 mm continuum emission underlying the CO( $J$ = 2 $\rightarrow$ 1) line and resolved radio continuum emission at 5 GHz in archival VLA data in both the foreground lensing galaxy and RXJ1131.

Based on our lens modeling analysis of different CO( $J$ = 2 $\rightarrow$ 1) velocity channels, we find a secondary CO-emitting source near RXJ1131 whose spatial position is consistent with those of an optically faint companion reported in previous optical studies (\al@Claeskens06a,Brewer08a; \al@Claeskens06a,Brewer08a). The magnification factor inferred for the CO emission in RXJ1131 is found to vary from $\mu_{\rm L}$ $\sim$ 3 to $\sim$ 9 across channels. This is indicative of an extended molecular gas distribution in the host galaxy of RXJ1131, where the different kinematic components of the gas are magnified inhomogeneously, similar to what was found for the $z$$>$ 4 quasar PSS J2322 $+$ 1944 (Riechers et al., 2008). Upon correcting for lensing magnification and subtracting a contribution from the companion, we find an intrinsically symmetric double-horned CO( $J$ = 2 $\rightarrow$ 1) line profile for RXJ1131. This together with a symmetric source-plane velocity gradient argues for a rotating disk in RXJ1131, in good agreement with previous findings (C06). Physical quantities derived for RXJ1131 and the companion throughout this paper are summarized in Table 6.

Based on the lensing-corrected CO( $J$ = 2 $\rightarrow$ 1) line intensities, we find an intrinsic gas mass of $M_{\rm gas}$ = $($ 1.38 $\pm$ 0.33 $)$ $\times 10^{10}$ $M_{\odot}$ for RXJ1131 and $($ 1.92 $\pm$ 0.09 $)$ $\times 10^{9}$ $M_{\odot}$ for the companion, corresponding to a gas mass ratio of $\sim$ 7:1. Using the source-plane CO( $J$ = 2 $\rightarrow$ 1) size of $R$ $\sim$ 6 kpc, we find a dynamical mass of $M_{\rm dyn}$ $\sim$ 8 $\times 10^{10}$ $M_{\odot}$ for RXJ1131. The dynamical gas mass fraction of $f_{\rm gas}$ = $M_{\rm gas}$ / $M_{\rm dyn}$ $\sim$ 18% and baryonic gas mass fraction of $f_{\textrm{mol}}$ = $M_{\rm gas}$ /( $M_{\rm gas}+M_{*}$ ) $\sim$ 34% are consistent with the trend of decreasing molecular gas content since $z$ $\sim$ 2 (e.g., Lagos et al., 2011; Tacconi et al., 2013, ;C13) which has been suggested as the cause for the decline in sSFR and cosmic star formation history towards $z$ $\sim$ 0 (e.g., Tacconi et al., 2013; Carilli & Walter, 2013; Genzel et al., 2015). The CO-based dynamical mass ratio of $\sim$ 24:1 between RXJ1131 and the companion, and a gas mass ratio of $\sim$ 7:1 suggest that the system is a gas-rich, “wet” minor merger.

Fitting dust SED models to the IR-to-mm photometry, we derive a lensing-corrected dust mass of $M_{\rm dust}$ $\sim$ 3 $\times 10^{8}$ $M_{\odot}$ , an infrared luminosity of $L_{\rm IR}$ $\sim$ 1.5 $\times 10^{12}$ $(5.5/\mu_{\rm L})$ $L_{\odot}$ , and a far-IR luminosity that corresponds to a SFRFIR $\sim$ 120 $M_{\odot}$ yr*-1*. These physical properties suggest that the merger system is dusty in nature with on-going star formation activity occuring at a rate comparable to local ULIRGs/mergers and high- $z$ massive disk galaxies (da Cunha et al., 2010; Daddi et al., 2010). We also derive a stellar mass of $M_{*}$ $\sim$ 3 $\times 10^{10}$ $M_{\odot}$ by fitting SED models to the rest-frame UV-to-mm photometry, which have been corrected for their respective magnification factors before performing the fit to account for differential lensing effect.

The source-plane distribution of the gas and stellar populations of different ages indicates that the CO gas is of similar spatial extent as the old and long-lasting stellar populations, whereas regions of recent star formation may be embedded within the molecular gas reservoir as a result of gas accumulation driven by interactions with the companion. Based on dynamical mass constraints, we cannot rule out the possibility that the compact star formation in the host galaxy may be heavily dust-obscured. Hence, the true extent of recent star formation may be as extended as the molecular gas reservoir.

While properties such as CO linewidth, SFR, and gas mass found in RXJ1131 are consistent with those of local ULIRGs and high- $z$ starburst galaxies, its SFE is comparable to those of nearby and high- $z$ disk galaxies rather than starburst systems. This is in good agreement with its disk-like kinematic signatures and its extended molecular gas distribution. We find a specific star formation rate (sSFR $\sim$ 4 Gyr*-1*) that is $\lesssim$ 1.5 $\sigma$ higher than those of “main sequence” galaxies. The slight elevation in sSFR over the main sequence suggests that the on-going star formation activity in RXJ1131 could be enhanced by interactions with the companion. Recent simulations have illustrated that the disk component of a gas-rich progenitor galaxy with low SFE can be retained upon merging since the efficiency at removing angular momentum of the gas via gravitational torques provided by stellar components is reduced in such a system (Springel & Hernquist, 2005; Robertson et al., 2006; Hopkins et al., 2009). As such, the extended gas disk of RXJ1131 together with its low SFE may indicate that the star formation in RXJ1131 could form a larger stellar bulge in the remnant disk galaxy upon coalescing. This picture is in agreement with the one based on the $M_{\rm BH}$$-$$M_{\rm bulge}$ relation, where we find an elevated $M_{\rm BH}$$/$$M_{\rm bulge}$ ratio of $>$ 0.27 ${}^{+0.11}_{-0.08}$ % for RXJ1131 compared to the local value. This suggests that the stellar bulge of RXJ1131 is still assembling in order to evolve onto the local relation.

We find that the redshift inferred from the NLR lines reported in previous studies are consistent with the systemic redshift as measured from the CO line, but that the BLR lines are redshifted by $\sim$ 715 km s*-1*. We raise several plausible scenarios that may explain the observed velocity offset, e.g., outflow/inflow of gas in the BLR, kinematics of the accretion disk, geometric effects, microlensing, and a recoiling black hole from merger event. The latter scenario might also explain the high black hole spin parameter of $a$ = 0.87 ${}^{+0.08}_{-0.15}$ reported by Reis et al. (2014), but further evidence is needed to confirm or rule out this scenario.

Theoretical studies have suggested that negative feedback from an AGN may remove a large fraction of the molecular gas from its host galaxy, thereby quenching its star formation (e.g., Di Matteo et al., 2005). In this study, we find that the star formation efficiency and specific SFR of RXJ1131 are comparable to those of $z$ $\sim$ 1 $-$ 1.5 disk galaxies, which are not known to host quasars, and that its molecular gas mass fraction is consistent with the observed cosmic decline for star-forming galaxies since $z$ $\sim$ 2 $-$ 3. Hence, we find no evidence of negative AGN feedback on the cold molecular gas fraction and on the star formation activity in RXJ1131. Future observations at higher resolution will allow us to better constrain the molecular gas kinematics and dynamics of RXJ1131 to investigate any potential interplay with the quasar on smaller scales. More broadly, systematic studies of the correlations between the molecular gas fraction, stellar mass, and AGN luminosity at different redshifts will enable us to better understand the relative importance of AGN feedback and of the evolution in the molecular gas mass fraction on the decline of star formation history and black hole accretion history.

We thank the referee for providing detailed and constructive comments that helped to improve the clarity of this manuscript. DR and RP acknowledge support from the National Science Foundation under grant number AST-1614213 to Cornell University. RP acknowledges support through award SOSPA3-008 from the NRAO. DR acknowledges the hospitality at the Aspen Center for Physics and the Kavli Institute for Theoretical Physics during part of the writing of this manuscript. This work is based on observations carried out under project number S14BX with the IRAM NOEMA Interferometer. IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain). Support for CARMA construction was derived from the Gordon and Betty Moore Foundation, the Kenneth T. and Eileen L. Norris Foundation, the James S. McDonnell Foundation, the Associates of the California Institute of Technology, the University of Chicago, the states of Illinois, California, and Maryland, and the National Science Foundation. Ongoing CARMA development and operations are supported by the National Science Foundation under a cooperative agreement and by the CARMA consortium universities. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. This research made use of data obtained with Herschel, an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA. This research has made use of NASA’s Astrophysics Data System Bibliographic Services. This work is based in part on observations made with the NASA/ESA Hubble Space Telescope, and obtained from the Hubble Legacy Archive, which is a collaboration between the Space Telescope Science Institute (STScI/NASA), the Space Telescope European Coordinating Facility (ST-ECF/ESA) and the Canadian Astronomy Data Centre (CADC/NRC/CSA). This work is based in part on observations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. This publication made use of data products from the Wide-field Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration. This publication made use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. This research made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. This research made use of Astropy, a community-developed core Python package for Astronomy (Astropy Collaboration et al., 2013). This research made use of APLpy, an open-source plotting package for Python hosted at http://aplpy.github.com. Facilities: IRAM PdBI, CARMA, VLA, Herschel(SPIRE), WISE, IRAS, 2MASS, Spitzer(IRAC, MIPS), HST(ACS, NICMOS)

Bibliography102

The reference list from the paper itself. Each links out to its DOI / PubMed record.

1Alexander et al. (2005) Alexander, D. M., Smail, I., Bauer, F. E., et al. 2005, Nature, 434, 738 · doi ↗
2Alexander et al. (2008) Alexander, D. M., Brandt, W. N., Smail, I., et al. 2008, AJ, 135, 1968 · doi ↗
3Andreani et al. (1993) Andreani, P., La Franca, F., & Cristiani, S. 1993, MNRAS, 261, L 35 · doi ↗
4Astropy Collaboration et al. (2013) Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A 33 · doi ↗
5Bonning et al. (2007) Bonning, E. W., Shields, G. A., & Salviander, S. 2007, Ap J, 666, L 13 · doi ↗
6Borys et al. (2005) Borys, C., Smail, I., Chapman, S. C., et al. 2005, Ap J, 635, 853 · doi ↗
7Bothwell et al. (2013) Bothwell, M. S., Smail, I., Chapman, S. C., et al. 2013, MNRAS, 429, 3047 · doi ↗
8Brewer & Lewis (2008) Brewer, B. J., & Lewis, G. F. 2008, MNRAS, 390, 39 (B 08) · doi ↗

TL;DR

Contribution

Findings

Abstract

Peer Reviews

Videos

Molecular Gas Kinematics and Star Formation Properties of the Strongly-Lensed Quasar Host Galaxy RXS J1131−-−1231

Abstract

Subject headings:

1. Introduction

2. Observations

2.1. CARMA CO(JJJ = 3 →\rightarrow→ 2)

2.2. PdBI CO(JJJ = 2 →\rightarrow→ 1)

2.3. Karl G. Jansky Very Large Array (Archival)

2.4. HST (Archival)

3. Results

3.1. CO(JJJ = 2 →\rightarrow→ 1) Emission

3.2. CO(JJJ = 3 →\rightarrow→ 2) Emission

3.3. Continuum Emission

3.4. Photometry

4. Analysis

4.1. Lens Modeling

4.1.1 Interpretation of the Source-plane Morphology

4.1.2 Spatial Extent and Differential Lensing

4.2. CO(JJJ = 2 →\rightarrow→ 1) Kinematics

4.3. CO(JJJ = 2 →\rightarrow→ 1) Dynamical Modeling

4.4. Dynamical Mass

4.5. SED Modeling

5. Discussion

5.1. ISM Properties

5.1.1 Linewidths and Sizes

5.1.2 Gas Mass Fractions and Gas-to-dust Ratio

5.1.3 Star Formation Efficiency and specific SFR

5.2. Systemic Redshift and Velocity Offset

5.3. The M_{\rm BH}$$-$$M_{\rm bulge} Relation

6. Summary and Conclusions

Molecular Gas Kinematics and Star Formation Properties of the Strongly-Lensed Quasar Host Galaxy RXS J1131 $-$ 1231

2.1. CARMA CO( $J$ = 3 $\rightarrow$ 2)

2.2. PdBI CO( $J$ = 2 $\rightarrow$ 1)

3.1. CO( $J$ = 2 $\rightarrow$ 1) Emission

3.2. CO( $J$ = 3 $\rightarrow$ 2) Emission

4.2. CO( $J$ = 2 $\rightarrow$ 1) Kinematics

4.3. CO( $J$ = 2 $\rightarrow$ 1) Dynamical Modeling

5.3. The $M_{\rm BH}$$-$$M_{\rm bulge}$ Relation