A new generation DIFMAP for Modelfitting Interferometric Data and Estimating Variances, Biases and Correlations

Agniva Roychowdhury,Eileen T. Meyer

2023-12-20

Abstract:We present the program `Catalogue of proper motions in extragalactic jets from Active galactic Nuclei with Very large Array Studies' or CAgNVAS, with the objective of using archival and new VLA observations to measure proper motions of jet components beyond hundred parsecs. This objective requires extremely high accuracy in component localization. Interferometric datasets are noisy and often lack optimal coverage of the visibility plane, making interpretation of subtleties in deconvolved imaging inaccurate. Fitting models to complex visibilities, rather than working in the imaging plane, is generally preferred as a solution when one needs the most accurate description of the true source structure. In this paper, we present a new generation version of $\texttt{DIFMAP}$ (\texttt{ngDIFMAP}) to model and fit interferometric closure quantities developed for the CAgNVAS program. \texttt{ngDIFMAP} uses a global optimization algorithm based on simulated annealing, which results in more accurate parameter estimation especially when the number of parameters is high. Using this package we demonstrate the ramifications of amplitude and phase errors, as well as loss of $u-v$ coverage, on parameters estimated from visibility data. The package can be used to accurately predict variance, bias, and correlations between parameters. Our results demonstrate the limits on information recovery from noisy interferometric data, with a particular focus on the accurate reporting of errors on measured quantities.

Instrumentation and Methods for Astrophysics

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to improve the accuracy of modeling and fitting radio interferometric data, especially when dealing with large - scale (more than 100 parsecs) observational data from active galactic nucleus (AGN) jets. Specifically, the author has developed a new version of DIFMAP (ngDIFMAP), aiming to estimate parameters more accurately through global optimization algorithms and closure quantity measurements, and to evaluate the variance, bias and correlation of these parameters.

### Main Problems

1. **Improve the accuracy of model fitting**:
- The traditional DIFMAP tool performs poorly under high - precision requirements, especially when the number of parameters is large.
- The new version ngDIFMAP introduces a global optimization algorithm based on simulated annealing to improve the accuracy of parameter estimation.

2. **Handle closure quantity measurements**:
- Closure quantity measurements (such as closure phase and closure amplitude) can partially eliminate the influence of antenna gain errors, thus providing a more robust estimate of the source structure.
- ngDIFMAP can fit these closure quantity measurements, further improving the reliability of the results.

3. **Evaluate uncertainties**:
- The influence of amplitude and phase errors and insufficient \(u - v\) coverage on parameter estimation has been studied.
- The Monte Carlo method and covariance matrix analysis are used to evaluate the variance, bias and correlation of parameters.

4. **Research on large - scale jets**:
- Through VLA archive data and new observational data, construct a large - scale active galactic nucleus jet dynamic catalogue (CAgNVAS) to study the physical properties of large - scale jets and their influence on the interstellar environment.

### Mathematical Formulas

- **Complex visibility**:
\[
V(u,v,w)=\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}A(l,m)I(l,m)e^{-i[ul + vm+w(\sqrt{1 - l^{2}-m^{2}}- 1)]}dl dm\sqrt{1 - l^{2}-m^{2}}
\]

- **Observed visibility**:
\[
\tilde{V}_{ij}=H_{ij}G_{i}G_{j}^{*}[V_{ij}+\epsilon_{ij}]
\]

- **Observed visibility with time - dependent error**:
\[
\tilde{V}_{ij}=[1 + Y_{ij}(t)]e^{j\delta\phi_{ij}}[V_{ij}+\epsilon_{ij}]
\]

- **Closure phase**:
\[
\tilde{\Psi}_{ijk}(t)=[\tilde{\Psi}_{ij}+\tilde{\Psi}_{jk}+\tilde{\Psi}_{ki}]=[\Psi_{ij}+\Psi_{jk}+\Psi_{ki}]+[\epsilon'_{ijk}]
\]

- **Closure phase variance**:
\[
\sigma_{clp}=\sqrt{\frac{\sigma_{ij}^{2}}{|\tilde{V}_{ij}|^{2}}+\frac{\sigma_{jk}^{2}}{|\tilde{V}_{jk}|^{2}}+\frac{\sigma_{ki}^{2}}{|\tilde{V}_{ki}|^{2}}}
\]

- **Closure amplitude**:
\[
A_{ijkl}=\frac{|\tilde{V}_{ij}||\tilde{V}_{kl}|}{|\tilde{V}_{ik}||\tilde{V}_{jl}|}=\frac{|V_{ij}||V_{kl}|}{|V_{ik}||V_{jl}|}
\]

- **Closure amplitude error**:
\[
\sigma_{cl - amp}=A_{ijkl}\sqrt{\frac{\sigma_{ij}^{2}}{|\til