Search

We report new radio observations of SDSS J090122.37+181432.3, a strongly lensed star-forming galaxy at $z=2.26$. We image 1.4 GHz (L-band) and 3 GHz (S-band) continuum using the VLA and 1.2 mm (band 6) continuum with ALMA, in addition to the CO(7-6) and CI(${\rm ^3P_2\rightarrow ^3\!P_1}$) lines, all at $\lesssim1.^{\prime\prime}7$ resolution. Based on the VLA integrated flux densities, we decompose the radio spectrum into its free-free (FF) and non-thermal components. The infrared-radio correlation (IRRC) parameter $q_{\rm TIR}=2.65_{-0.31}^{+0.24}$ is consistent with expectations for star forming galaxies. We obtain radio continuum-derived SFRs that are free of dust extinction, finding $\rm {620}_{-220}^{+280}\,M_\odot\,yr^{-1}$, $\rm {230}_{-160}^{+570}\,M_\odot\,yr^{-1}$, and $\rm {280}_{-120}^{+460}\,M_\odot\,yr^{-1}$ from the FF emission, non-thermal emission, and when accounting for both emission processes, respectively, in agreement with previous results. We estimate the gas mass from the CI(${\rm ^3P_2\rightarrow ^3\!P_1}$) line as $M_{\rm gas}=(1.2\pm0.2)\times10^{11}\,M_\odot$, which is consistent with prior CO(1-0)-derived gas masses. Using our new IR and radio continuum data to map the SFR, we assess the dependence of the Schmidt-Kennicutt relation on choices of SFR and gas tracer for $\sim{\rm kpc}$ scales. The different SFR tracers yield different slopes, with the IR being the steepest, potentially due to highly obscured star formation in J0901. The radio continuum maps have the lowest slopes and overall fidelity for mapping the SFR, despite producing consistent total SFRs. We also find that the Schmidt-Kennicutt relation slope is flattest when using CO(7-6) or CI(${\rm ^3P_2\rightarrow ^3\!P_1}$) to trace gas mass, suggesting that those transitions are not suitable for tracing the bulk molecular gas in galaxies like J0901., Comment: 28 pages, 13 figures, 9 tables

Strong gravitational lensing offers constraints on the Hubble constant that are independent of other methods. However, those constraints are subject to uncertainties in lens models. Previous studies suggest that using an elliptical power law + external shear (EPL+XS) for the lensing galaxy can yield results that are precise but inaccurate. We examine such models by generating and fitting mock lenses which produces multiple images of a background quasar-like point source. Despite using the same model for input and output, we find statistical bias in the Hubble constant on the order of 3% to 5%, depending on whether the elliptical lenses have noise or not. The phase space distribution has a `flared' shape that causes the mass power law slope to be underestimated and the Hubble constant to be overestimated. The bias varies with image configuration, which we quantify through annulus length between images with the first and second time delays ($\Delta r_{1,2}$). The statistical bias is worse for configurations that have narrow annuli (e.g., symmetric cross configurations). Assuming a source at redshift 2.0 and an EPL+XS lens at redshift 0.3, we find that the bias can be reduced, but not eliminated, if we limit the sample to systems with annulus lengths $\Delta r_{1,2} \gtrsim 0.3$ arcsec. As lens samples grow, it may be helpful to prioritize this range of image configurations for follow-up observation and analysis.

We consider gravitational lensing of a background source by a finite system of point-masses. The problem of determining the maximum possible number of lensed images has been completely resolved in the single-plane setting (where the point masses all reside in a single lens plane), but this problem remains open in the multiplane setting. We construct examples of $K$-plane point-mass gravitational lens ensembles that produce $\prod_{i=1}^K (5g_i-5)$ images of a single background source, where $g_i$ is the number of point masses in the $i^\text{th}$ plane. This gives asymptotically (for large $g_i$ with $K$ fixed) $5^K$ times the minimal number of lensed images. Our construction uses Rhie's single-plane examples and a structured parameter-rescaling algorithm to produce preliminary systems of equations with the desired number of solutions. Utilizing the stability principle from differential topology, we then show that the preliminary (nonphysical) examples can be perturbed to produce physically meaningful examples while preserving the number of solutions. We provide numerical simulations illustrating the result of our construction, including the positions of lensed images as well as the structure of the critical curves and caustics. We observe an interesting ``caustic of multiplicity'' phenomenon that occurs in the nonphysical case and has a noticeable effect on the caustic structure in the physically meaningful perturbative case., Comment: 20 pages, 8 figures. The paper will appear in the Journal of Mathematical Physics

The observation of the bending of light by mass, now known as gravitational lensing, was key in establishing general relativity as one of the pillars of modern physics. In the past couple of decades, there has been increasing interest in using gravitational lensing to test general relativity beyond the weak deflection limit. Black holes and neutron stars produce the strong gravitational fields needed for such tests. For a rotating compact object, the distinction between prograde and retrograde photon trajectories becomes important. In this paper, we explore subtleties that arise in interpreting the bending angle in this context and address the origin of seemingly contradictory results in the literature. We argue that analogies that cannot be precisely quantified present a source of confusion.

Integral field spectroscopy of high-redshift galaxies has become a powerful tool for understanding their dynamics and evolutionary states. However, in the case of gravitationally lensed systems, it has proved difficult to model both lensing and intrinsic kinematics in a way that takes full advantage of the information available in the spectral domain. In this paper, we introduce a new method for pixel-based source reconstruction that alters standard regularization schemes for two-dimensional data in a way that leverages kinematic information in a physically motivated but flexible fashion, and that is better suited to the three-dimensional nature of integral field data. To evaluate the performance of this method, we compare its results to those of a more traditional two-dimensional non-parametric approach using mock ALMA observations of a typical high-redshift dusty star-forming galaxy. We find that 3D regularization applied to an entire data cube reconstructs a source's intensity and velocity structure more accurately than 2D regularization applied to separate velocity channels. Cubes reconstructed with 3D regularization also have more uniform noise and resolution properties and are less sensitive to the signal-to-noise ratio of individual velocity channels than the results of 2D regularization. Our new approach to modeling integral field observations of lensed systems can be implemented without making restrictive a priori assumptions about intrinsic kinematics, and opens the door to new observing strategies that prioritize spectral resolution over spatial resolution (e.g., for multi-configuration arrays like ALMA)., Comment: Updated to match published version; 22 pages, 11 figures

We report X-ray observations of the most distant known gravitationally lensed quasar, J0439+1634 at $z=6.52$, which is also a broad absorption line (BAL) quasar, using the XMM-Newton Observatory. With a 130 ks exposure, the quasar is significantly detected as a point source at the optical position with a total of 358$^{+19}_{-19}$ net counts using the EPIC instrument. By fitting a power-law plus Galactic absorption model to the observed spectra, we obtain a spectral slope of $\Gamma=1.45^{+0.10}_{-0.09}$. The derived optical-to-X-ray spectral slope $\alpha_{\rm{ox}}$ is $-2.07^{+0.01}_{-0.01}$, suggesting that the X-ray emission of J0439+1634 is weaker by a factor of 18 than the expectation based on its 2500 Angstrom luminosity and the average $\alpha_{\rm{ox}}$ vs. luminosity relationship. This is the first time that an X-ray weak BAL quasar at $z>6$ has been observed spectroscopically. Its X-ray weakness is consistent with the properties of BAL quasars at lower redshift. By fitting a model including an intrinsic absorption component, we obtain intrinsic column densities of $N_{\rm{H}}=2.8^{+0.7}_{-0.6}\times10^{23}\,\rm{cm}^{-2}$ and $N_{\rm{H}}= 4.3^{+1.8}_{-1.5}\times10^{23}\,\rm{cm}^{-2}$, assuming a fixed $\Gamma$ of 1.9 and a free $\Gamma$, respectively. The intrinsic rest-frame 2--10 keV luminosity is derived as $(9.4-15.1)\times10^{43}\,\rm{erg\,s}^{-1}$, after correcting for lensing magnification ($\mu=51.3$). The absorbed power-law model fitting indicates that J0439+1634 is the highest redshift obscured quasar with a direct measurement of the absorbing column density. The intrinsic high column density absorption can reduce the X-ray luminosity by a factor of $3-7$, which also indicates that this quasar could be a candidate of intrinsically X-ray weak quasar., Comment: 8 pages, 3 figures, 1 table; Accepted for publication in ApJL