---
title: Milliarcsecond X-ray Astrometry
url: https://www.emergentmind.com/topics/milliarcsecond-x-ray-astrometry
type: topic
---

# Milliarcsecond X-ray Astrometry

Searching arXiv for recent papers on milliarcsecond X-ray astrometry, including lensed AGN and Chandra astrometric methods.
Milliarcsecond X-ray astrometry denotes the localization of X-ray emission with positional precision at the milli-arcsecond scale, typically by combining high-resolution X-ray imaging with external astrometric references or with the geometric amplification supplied by strong gravitational lensing. In the contemporary literature, the method has been developed most explicitly for strongly lensed AGN and quasars, where parametric lens modelling and Bayesian or maximum-likelihood inference are used to infer source-plane X-ray positions far below the native angular resolution of Chandra. Reported applications include the localization of X-ray emission in CLASS B0712+472 and CLASS B1608+656 within 11 mas and 9 mas from the radio source, respectively [2203.04245], the measurement of a 1$\sigma$ quasi-elliptical uncertainty region of $0.5 \times 1.3$ milli-arcsecond for HE 0435-1223 [2505.08077], and a 99.7% confidence constraint of $0.\!\!''030 \times 0.\!\!''020$ for the radio-quiet quasar GraL J065904.1+162909 [2509.22797]. Related Chandra work on pulsars shows that Gaia-referenced frame alignment can also deliver absolute positions with typical astrometric precision $\sim 10$ mas and proper-motion statistical uncertainties down to $1.3$ mas yr$^{-1}$ [2602.18436].

## 1. Strong lensing as an X-ray astrometric amplifier

The gravitational-lensing formulation used in this area is the standard lens equation
$$
\beta = \theta - \nabla \psi(\theta),
$$
or equivalently, in one HE 0435-1223 analysis,
$$
\theta = \beta + \nabla \psi(\theta),
$$
with $\theta$ the image-plane coordinate, $\beta$ the source-plane coordinate, and $\psi(\theta)$ the scaled two-dimensional lens potential [2505.08077] [2509.22797]. In the J0659 analysis, the deflector is modeled with two mass components, including a Singular Isothermal Ellipsoid plus external shear for “deflector 1” and a singular isothermal sphere for “deflector 2” [2509.22797].

The operative principle is that strong lensing is achromatic but highly astigmatic: near the caustics, small separations in the source plane are stretched by factors of order $\mu^t$ along one direction and compressed along the orthogonal direction [2505.08077]. In mathematical terms, the local source-to-image mapping is described by the magnification matrix
$$
A = \partial \beta / \partial \theta = I - \partial^2 \psi / \partial \theta \partial \theta,
$$
whose eigenvalues are $(1/\mu_r, 1/\mu_t)$ [2505.08077]. This lens-mediated anisotropic amplification is the basis for obtaining source-plane constraints that are much smaller than the $\sim 0.5''$ Chandra PSF.

In the lensed-AGN literature, this framework is used to address a specific observational limitation: high-$z$ AGN are too distant to be spatially resolved with current or upcoming X-ray facilities [2505.08077]. A plausible implication is that lensing functions not merely as a flux amplifier but as a geometry-dependent astrometric transform, enabling relative localization of emission components that would otherwise remain unresolved.

## 2. Parametric lens modelling and external astrometric anchors

Milliarcsecond X-ray astrometry in lensed AGN depends on a lens model that reproduces the observed image configuration to very high precision. For GraL J065904.1+162909, Gaia DR3 and HST observations are used to infer a mass model for the deflector that reproduces the positions of the quasar lensed images to milliarcsecond precision [2509.22797]. The analysis employs LENSMODEL under a flat $\Lambda$CDM cosmology with $H_0=67.5$ and $\Omega_m=0.308$, minimizes
$$
\chi^2 = \sum_{i=1..4} \left[\left(\frac{RA_{{obs},i}-RA_{{mod},i}}{\sigma_{RA,i}}\right)^2 + \left(\frac{Dec_{{obs},i}-Dec_{{mod},i}}{\sigma_{Dec,i}}\right)^2\right],
$$
and reports a best-fit yielding $\chi^2 \approx 10^{-3}$, reproducing each lensed image to $\le 1$ mas [2509.22797].

For HE 0435-1223, optical astrometry comes from Gaia DR3, which detects all four lensed images with sub-mas precision, $\sigma_{\rm RA}, \sigma_{\rm Dec} \simeq 0.1$–$0.2$ mas [2505.08077]. The best Gaia-constrained lens model is an SIE+external shear, and the “zero-point” of co-registration is set by the four Gaia image positions mapped back to a single source position in the source plane [2505.08077]. In this construction, Gaia supplies the absolute relative geometry against which X-ray image centroids can be registered.

This reliance on external astrometric anchors also appears outside the lensing context. In the Chandra pulsar proper-motion study, unresolved X-ray flux from stars in the Gaia catalog is used in addition to X-ray bright point sources for alignment, improving uncertainties [2602.18436]. That work selects “classical” X-ray point sources and Gaia DR3 stars, assigns each source a PSF model using a precomputed MARX-based library, and fits simultaneously for frame translations, source positions, pulsar position, and proper motion [2602.18436]. Although its primary target class is different, it establishes the general importance of Gaia as an absolute astrometric reference frame for milliarcsecond-regime X-ray analyses.

## 3. Statistical inference, PSF forward modelling, and confidence regions

A central methodological feature of the lensed-quasar studies is pixel-level forward modelling of Chandra data. For HE 0435-1223, each ObsID is modeled by simulating high-fidelity Chandra PSFs of four point sources at the lens-predicted image positions using SAOTrace + MARX with the actual aspect solution and energy spectrum of each image [2505.08077]. The data are cut to a $6'' \times 6''$ region around the target and binned to $0.246''$ pixels. If $n_i$ denotes the observed counts in pixel $i$ and $\lambda_i(\Delta x,\Delta y; A_1 \ldots A_4)$ the model prediction, the fit maximizes the Poisson likelihood or equivalently minimizes the Cash statistic
$$
C(\Delta x,\Delta y,\{A_k\}) = -2 \ln L
= -2 \sum_i \left[n_i \ln \lambda_i - \lambda_i - \ln(n_i!)\right].
$$
At each trial $(\Delta x,\Delta y)$, the image flux normalizations are adjusted to minimize $C$ [2505.08077].

For J0659, the Bayesian formulation is explicit. A triangular grid of 200 trial source positions is selected near the inner caustic. For each source position, 1,000 MARX ray-traces are run, yielding predicted counts $N_{{mod},i,j}(k)$ in pixel $k$ of image $i$. The likelihood is
$$
L(d|s_j) = \prod_{i=1..4} \prod_k {\rm Poiss}\!\left(N_{{obs},i}(k)\mid \alpha_i N_{{mod},i,j}(k) + B_i(k)\right),
$$
with per-image normalizations $\alpha_i$, and with a flat prior the posterior satisfies
$$
P(s_j|d) \propto \exp[\ell(s_j)].
$$
Defining $\Delta \chi^2(s_j) = -2[\ell(s_j)-\ell_{\max}]$, the 2D contours $\Delta \chi^2=2.30, 6.18, 11.83$ correspond to 68.3%, 95.5%, and 99.7% confidence levels under 2 d.o.f. [2509.22797].

Confidence regions are then extracted either from Wilks’s theorem or from the Hessian of the objective function. In HE 0435-1223, the residual uncertainty in $(\Delta x,\Delta y)$ for each ObsID is derived from Wilks’s theorem, so that $\Delta C = C-C_{\min}$ is distributed as $\chi^2$ with 2 degrees of freedom [2505.08077]. In J0659, around the posterior peak the $2 \times 2$ covariance matrix $\Sigma$ is computed from the Hessian of $-\ell$, and in the rotated principal frame the 99.7% contour satisfies
$$
[\Delta x' \ \Delta y'] \Sigma^{-1}
\begin{bmatrix}
\Delta x'\\
\Delta y'
\end{bmatrix}
= 11.83.
$$
Its semi-axes are
$$
a = \sqrt{11.83\,\lambda_1}, \qquad b = \sqrt{11.83\,\lambda_2},
$$
with orientation
$$
\phi = \tfrac12 \atan\!\left[\frac{2\Sigma_{xy}}{\Sigma_{xx}-\Sigma_{yy}}\right]
$$
[2509.22797].

The same likelihood-centric design appears in Gaia-referenced Chandra astrometry for pulsars. There, events are modeled by
$$
L_{i,s,f} = f_s\,P_{s,f}\bigl(y_i-\hat y_{i,s,f}\bigr) + (1-f_s)\,\frac{1}{\pi R_s^2\,C_{{exp},f,s}},
$$
with
$$
\hat y_{i,s,f} = (x_s-\delta_f)+t_i(\mu_s-\dot\delta_f),
$$
and the total likelihood $L=\prod_{s,f,i}L_{i,s,f}$ is maximized, after which the Hessian of $-\ln L$ is inverted to yield a Gaussian covariance matrix $\Sigma$ for all fit parameters [2602.18436]. This suggests that milliarcsecond X-ray astrometry is now defined less by a single instrument capability than by a coupled framework of forward PSF modelling, external registration, and global likelihood inference.

## 4. Observational implementations and reported localization performance

The best-documented lensed-AGN results span several source classes and redshifts.

| System | Reported astrometric result | Principal implication |
|---|---|---|
| CLASS B0712+472 | X-ray source spatially located within 11 mas from the radio source | X-ray emission co-spatial with radio and optical [2203.04245] |
| CLASS B1608+656 | X-ray source spatially located within 9 mas from the radio source | X-ray co-spatial with radio, but displaced with respect to optical emission at $1\sigma$ level [2203.04245] |
| HE 0435-1223 | $1\sigma$ quasi-elliptical region of $0.5 \times 1.3$ mas | Evidence for a projected 3 mas optical–X-ray offset [2505.08077] |
| GraL J065904.1+162909 | 99.7% ellipse of $0.\!\!''030 \times 0.\!\!''020$ | Maximum possible X-ray–optical distance 44.7 mas [2509.22797] |

For CLASS B0712+472 and CLASS B1608+656, the 2022 study reports the localization of the X-ray emission from two strongly lensed AGN at $z=1.34$ and $z=1.394$, respectively, using a novel method that combines parametric lens modelling with a Bayesian analysis [2203.04245]. The paper states that this high astrometric precision improves on the limitations of existing X-ray instruments by two orders of magnitude and opens a path to search for offset and binary AGN at $z>1$ [2203.04245].

For HE 0435-1223, the X-ray imaging consists of 11 ACIS-S observations on Chandra, totaling 285 ks and $\sim 10\,000$ counts [2505.08077]. The X-ray source position is measured within a $1\sigma$ quasi-elliptical region of $0.5 \times 1.3$ milli-arcsecond, corresponding to about $150$ pc$^2$ at $z=1.689$ [2505.08077]. The authors state that, by referencing all X-ray images to the Gaia positions under the same lens model, they effectively achieve tens-of-$\mu$as relative astrometry in the source plane even though the raw Chandra PSF is $\sim 0.5''$ [2505.08077].

For J0659 at $z=3.083$, the inferred mass model reproduces the image positions to milliarcsecond precision, and archival Chandra observations are used to constrain the X-ray origin to be within a $0.\!\!''030 \times 0.\!\!''020$ ellipse at the 99.7% level [2509.22797]. The methodology is described as enhancing Chandra’s spatial resolution at high-$z$ by a factor of 6 [2509.22797].

Outside the lensing regime, the pulsar study reports absolute positions referenced to Gaia with typical astrometric precision $\sim 10$ mas and proper-motion statistical uncertainties down to $1.3$ mas yr$^{-1}$, which the abstract describes as the most precise X-ray PM achieved to date [2602.18436]. The paper’s scope is proper motion rather than source-plane deprojection, but its inclusion is relevant because it demonstrates that Gaia-tied Chandra astrometry can operate in the same milli-arcsecond scale, albeit with different observables and different error budgets.

## 5. Optical–X-ray offsets, co-spatiality, and astrophysical interpretation

The empirical motivation for milliarcsecond X-ray astrometry is the possibility that the X-ray and optical centroids of a distant AGN need not coincide. In CLASS B0712+472, the X-ray emission is reported to be co-spatial with the radio and optical emission [2203.04245]. In CLASS B1608+656, the X-ray emission is co-spatial with radio but displaced with respect to the optical emission at $1\sigma$ level, and the source is therefore positioned as an offset AGN candidate [2203.04245].

HE 0435-1223 provides the clearest quantified optical–X-ray displacement. The best X-ray source position is reported as $(\Delta x,\Delta y)=(-2.99,+0.16)$ mas relative to the Gaia optical centroid at $(0,0)$, yielding
$$
\Delta \theta = \sqrt{(\Delta RA)^2 + (\Delta Dec)^2} = 3.0 \pm 0.5 \ {\rm mas},
$$
or $\simeq 26 \pm 4$ pc [2505.08077]. To test coincidence, the difference in minimum Cash statistic between the global best fit and the best solution constrained to lie at the Gaia position is $\Delta C \simeq 23.7$; since $\Delta C \sim \chi^2(2)$, the corresponding $p$-value is $\simeq 3.9\times 10^{-4}$, i.e. a $3.36\sigma$ rejection of coincidence [2505.08077]. The offset is interpreted as most likely associated to a portion of the X-ray emission arising from an X-ray jet or outflow [2505.08077].

For J0659, the most probable X-ray source coincides with the optical source predicted by the lens model [2509.22797]. At 99.7% confidence, however, the maximum possible separation on the contour is
$$
\Delta_{\max} = \sqrt{a^2+b^2} = 0.0447'' \simeq 44.7 \ {\rm mas},
$$
which corresponds to $\lesssim 340$ pc at $z_s=3.083$ for $1''=7.52$ kpc [2509.22797]. The same pipeline can be applied to energy-filtered events, giving soft $(0.5$–$2$ keV$)$ emission localized to $0.040'' \times 0.020''$ and hard $(2$–$7$ keV$)$ emission localized to $0.030'' \times 0.015''$, both consistent with the optical position though with peak positions differing by $\sim 30$ mas [2509.22797].

These results bear directly on the interpretation of quasar inner structure. The authors of the lensed-quasar papers connect such offsets to X-ray jets, outflows, offset AGN candidates, and the possible presence of binary/offset AGN systems [2203.04245] [2505.08077]. A plausible implication is that positional non-coincidence can function as a morphological diagnostic when direct imaging of the underlying structure is unavailable.

## 6. X-ray varstrometry, dual-source hypotheses, and broader extensions

In HE 0435-1223, the astrometric framework is extended to “X-ray varstrometry,” a centroid-jitter test for unresolved multiple emitters. If two X-ray-emitting components are separated by $\theta$ and at least one varies in flux, the centroid jitter has rms amplitude
$$
\sigma_{astro} = \theta \cdot \left[\frac{q}{1+q}\right] \cdot \sqrt{\frac{\sigma_F^2}{\langle F\rangle^2}},
$$
where $q \le 1$ is the mean flux ratio, $\langle F\rangle$ the mean total flux, and $\sigma_F$ its rms variability [2505.08077]. For HE 0435-1223, $\sigma_F/\langle F\rangle \simeq 14\%$ is measured after excluding two epochs strongly affected by microlensing, and the measured perpendicular jitter is $\sigma_{noise} \lesssim 3.3$ mas [2505.08077]. Imposing $\sigma_{astro} \le 3.3$ mas and $q=1$ gives $\theta \le 47$ mas $(1\sigma)$, or $\lesssim 410$ pc at $z=1.689$ [2505.08077]. Thus X-ray varstrometry places upper limits on sub-kpc dual/offset AGN.

J0659 presents a complementary model-selection perspective. The paper states that, in future applications, single-source and two-source hypotheses can be compared through a Bayes factor
$$
B_{2/1} =
\frac{\iint L(d|s_1,s_2)\pi(s_1,s_2)\,ds_1 ds_2}
{\int L(d|s)\pi(s)\,ds},
$$
with the aim of detecting dual or binary AGN at projected separations $\lesssim 10$–$100$ pc at $z\sim 2$ [2509.22797]. This is presented as an extension rather than a demonstrated detection in that source.

The broader methodological extension is toward spectrally resolved astrometry. The J0659 work explicitly describes a novel approach that extends the methodology toward investigating the origin of the soft and hard X-ray emitting regions in lensed quasars [2509.22797]. It further elaborates on the potential of upcoming broadband and spectrally resolved X-ray astrometric studies to probe complex quasar morphology at high-$z$ and to identify dual and binary AGN candidates [2509.22797].

## 7. Limitations, precision floors, and research prospects

The current literature shows that milli-arcsecond X-ray astrometry is not a property of the detector alone but of the entire inference chain: external astrometry, lens or PSF model fidelity, count statistics, and registration strategy. In the pulsar study, a point source with $N_s$ counts and PSF rms width $r_s$ has positional uncertainty
$$
\sigma_s \simeq r_s/\sqrt{N_s},
$$
and the frame-registration precision in the background-free approximation satisfies
$$
\frac{1}{\sigma_{\rm frame}^2} = \sum_s \frac{1}{\sigma_s^2} = \sum_s \frac{N_s}{r_s^2},
$$
where the sum is over field point sources including Gaia stars [2602.18436]. This makes explicit that dense reference grids of faint sources can materially tighten the astrometric solution.

Systematics remain important. By comparison to the VLBI proper motion of PSR B2224+65, the pulsar analysis estimates a systematic floor $\sigma_{sys}\approx 2$ mas yr$^{-1}$, while other tests indicate systematic errors $\lesssim 2$–$3$ mas yr$^{-1}$ [2602.18436]. For lensed quasars, the reported precisions and contours depend on the adopted lens model, the fidelity of MARX or SAOTrace-based PSF simulations, and the validity of the co-registration assumptions across epochs [2505.08077] [2509.22797]. This suggests that the achieved milliarcsecond performance is an inferred source-plane or Gaia-referenced quantity, not a direct imaging resolution in the conventional sense.

The prospects articulated in the literature are correspondingly specific. The 2022 lensed-AGN study states that the demonstrated method opens a path to search for offset and binary AGN at $z>1$, and to directly test supermassive black hole formation models in a redshift range that has been mostly underconstrained to date [2203.04245]. The HE 0435-1223 paper states that the technique opens the door to sub-kiloparsec X-ray imaging at high redshift for hundreds to thousands of strongly lensed quasars [2505.08077]. The J0659 study points to upcoming broadband and spectrally resolved applications and notes that the methodology has been applied, so far, to five lensed quasars, including J0659 [2509.22797].

Taken together, these studies define milliarcsecond X-ray astrometry as an emerging high-precision regime in which Chandra data, Gaia astrometry, parametric lens models, and likelihood-based forward modelling are combined to test co-spatiality, measure offsets, and constrain unresolved multi-component X-ray structure in distant systems [2203.04245] [2505.08077] [2509.22797] [2602.18436].

Source: https://www.emergentmind.com/topics/milliarcsecond-x-ray-astrometry