
JetCurry I. Reconstructing Three-Dimensional Jet Geometry from
Two-Dimensional Images
Sailee M. Sawanta,1,Katie Kosaka,c,1,Kunyang Lia,d,1,Sayali S. Avachata,e,1,Eric S. Perlmanaand
Debasis Mitrab
aDepartment of Aerospace, Physics and Space Sciences, Florida Institute of Technology, 150 W. University Blvd., Melbourne, FL, 32901, USA
bDepartment of Computer Engineering and Sciences, Florida Institute of Technology, 150 W. University Blvd., Melbourne, FL, 32901, USA
cPhysics Department, University of Warwick, Coventry, CV4 7AL, UK
dCenter for Relativistic Astrophysics, School of Physics, Georgia Institute of Technology, 837 State Street, Atlanta, GA, 30332, USA
eInter-University Centre for Astronomy and Astrophysics, Pune, Maharashtra, 411007, India
ARTICLE INFO
Keywords:
galaxies: jets
galaxies: individual (M87)
methods: numerical
Applied computing: physical sciences
and engineering
Computing methodologies: modeling
and simulation
ABSTRACT
We present a three-dimensional (3-D) visualization of jet geometry using numerical methods based
on a Markov Chain Monte Carlo (MCMC) and limited memory Broyden–Fletcher–Goldfarb–Shanno
(BFGS) optimized algorithm. Our aim is to visualize the 3-D geometry of an active galactic nucleus
(AGN) jet using observations, which are inherently two-dimensional (2-D) images. Many AGN jets
display complex structures that include hotspots and bends. The structure of these bends in the jet’s
frame may appear quite different than what we see in the sky frame, where it is transformed by our
particular viewing geometry. The knowledge of the intrinsic structure will be helpful in understanding
the appearance of the magnetic field and hence emission and particle acceleration processes over the
length of the jet. We present the JetCurry algorithm to visualize the jet’s 3-D geometry from its 2-D
image. We discuss the underlying geometrical framework and outline the method used to decompose
the 2-D image. We report the results of our 3-D visualization of the jet of M87, using the test case
of the knot D region. Our 3-D visualization is broadly consistent with the expected double helical
magnetic field structure of knot D region of the jet. We also discuss the next steps in the development
of the JetCurry algorithm.
1. Introduction
Relativistic jets transport mass and energy from sub-
parsec central regions to Mpc-scale lobes, with a kinetic
power comparable to that of their host galaxies and active
galactic nuclei (AGNs). This profoundly influences the evo-
lution of the hosts, nearby galaxies, and the surrounding in-
terstellar and intracluster medium (Silk et al.,2012;Fabian,
2012). The generation of such flows is tied to the process
of accretion onto (likely) rotating black holes, where the
magneto-rotational instability can couple the black hole’s
spin and magnetic field to the disk or flow to produce high-
latitude outflows at speeds close to the speed of light (Meier
et al.,2001). While these jets have a dominant direction of
motion (i.e., outward from the black hole), they often have
bends and features (both stationary and moving) that are
either perpendicular or aligned relative to the jet at some
angle. Deciphering the true nature of these features and
their geometry, relation, and dynamical meaning within the
flow is a difficult problem, as any astronomical images we
acquire are of necessity two-dimensional (2-D) views of
three-dimensional (3-D) objects.
The problem of reconstructing 3-D information from 2-
D images is common to many fields, but it is particularly
critical in astronomy. In most other cases, e.g., medical
ssawant2011@my.fit.edu (S.M. Sawant)
ORCID(s): 0000-0002-7987-0310 (S.M. Sawant); 0000-0002-0867-8946
(K. Li); 0000-0002-5550-5693 (S.S. Avachat); 0000-0002-3099-1664 (E.S.
Perlman); 0000-0002-4351-1252 (D. Mitra)
1These authors contributed equally.
imaging, one may take images of a source from multiple
viewpoints to aid reconstruction. However, this is not pos-
sible in astronomy, so we must rely on other methods. For
instance, Steffen et al. (2011), Wenger et al. (2012), Wenger
et al. (2013), Cormier (2013), Sabatini et al. (2018), and
Lagattuta et al. (2019) used symmetries inherent in, respec-
tively, planetary nebulae and galaxies, plus 2-D images, to
infer and reconstruct 3-D visualizations of these objects.
This field is, in fact, rapidly growing in astronomy, as can
be seen by the vast number of subjects explored on the
3DAstrophysics blog2.
In astrophysical jets, the problem is rather different.
Unlike in galaxies or planetary nebulae, we cannot make
assumptions such as spherical, elliptical or disk symmetry,
or rotation. However, we can assume a dominant direction of
propagation. As an example of the typical knotted structure
of AGN jets, we show in Figure 1, a broad view of the M87
jet, one of the nearest of the class at 16.7 Mpc distance,
taken from Meyer et al. (2013). In every single image, the
M87 jet shows an amazing complexity of features, including
knots, helical undulations, shocks, and a variety of other
morphological structures, many of which are oriented at
some odd angles with respect to the overall jet direction. As
shown by Meyer et al. (2013), some of the features in the
M87 jet seem to move with apparent velocities up to about
6𝑐within the inner 12.
′′0of the jet, with a general decline
in apparent speed with increasing distance from the nucleus.
However, there are some nearly stationary components that
2https://3dastrophysics.wordpress.com/
Sawant, Kosak, Li, Avachat et al.: Preprint submitted to Elsevier. Accepted. Page 1 of 10
arXiv:2210.03033v1 [astro-ph.HE] 6 Oct 2022