Introducing the Condor Array Telescope: III. The expansion and age of the shell of the dwarf nova Z Camelopardalis, and detection of a second, larger shell

The existence of a vast nova shell surrounding the prototypical dwarf nova Z Camelopardalis (Z Cam) proves that some old novae undergo metamorphosis to appear as dwarf novae thousands of years after a nova eruption. The expansion rates of ancient nova shells offer a way to constrain both the time between nova eruptions and the time for post-nova mass transfer rates to decrease significantly, simultaneously testing nova thermonuclear runaway models and hibernation theory. Previous limits on the expansion rate of part of the Z Cam shell constrain the inter-eruption time between Z Cam nova events to be $>$ 1300 years. Deeper narrow-band imaging of the ejecta of Z Cam with the Condor Array Telescope now reveals very low surface brightness areas of the remainder of the shell. A second, even fainter shell is also detected, concentric with and nearly three times the size of the"inner"shell. This is the first observational support of the prediction that concentric shells must surround the frequently-erupting novae of relatively massive white dwarfs. The Condor images extend our Z Cam imaging baseline to 15 years, yielding the inner shell's expansion rate as $v = 83 \pm 37$ km s$^{-1}$ at 23 degrees South of West, in excellent agreement with our 2012 prediction. This velocity corresponds to an approximate age of $t = 2672^{-817}_{+2102}$ yr. While consistent with the suggestion that the most recent nova eruption of Z Cam was the transient recorded by Chinese Imperial astrologers in the year 77 BCE, the age uncertainty is still too large to support or disprove a connection with Z Cam.


INTRODUCTION
Cataclysmic variables (CVs) all comprise a white dwarf (WD) accreting matter from its binary companion (Warner 1995).While the donor stars can be helium WDs, red giants or subgiants, the vast majority of known CVs contain approximately solar composition main sequence star or brown dwarf secondaries.The recent accretion rate / determines the subtype of CV that we observe, while the history of / over multiple Gyr drives the binary's long-term evolution.
The continued accretion of hydrogen-rich matter onto a WD results in an increasingly dense and massive envelope supported by degenerate electron pressure.Nuclear reactions and their concomitant heating near the accreted envelope base do not, at first, lead to envelope expansion because degenerate electron pressure is temperature independent.Once the nuclear heating timescale becomes ★ E-mail: mshara@amnh.org(MMS) shorter than the timescale required to carry away heat (Shara 1981) a thermonuclear runaway ensues.The resulting classical nova (Starrfield et al. 1972;Prialnik et al. 1978) fuses several percent of its envelope's hydrogen to helium in a few minutes, generating enough energy to eject the accreted envelope and achieving peak luminosities up to 10 6  ⊙ .Nova eruptions cease when all or most of the accreted envelope has been ejected (Prialnik et al. 1979).
Nova shells are expected to remain bright enough to be detected for multiple centuries (Tappert et al. 2020) after eruption.The likely sixcentury-old shell of Nova Sco 1437 CE (Shara et al. 2017b) is easily imaged with 1-m-class telescopes.A recent H  imaging survey of 47 CVs focused on (high mass transfer rate) nova-like variables, but did not detect any new nova shells (Sahman & Dhillon 2022).We will return to this important negative result in Section 3, where we demonstrate that significantly lower surface brightness imaging than has hitherto been accomplished is essential to detecting the oldest nova ejecta.
The envelope masses required to initiate a nova thermonuclear runaway on a WD range from ∼ 10 −7  ⊙ for WDs close to 1.4 ⊙ , to as large as ∼ 10 −3  ⊙ for CV WDs in the 0.6 ⊙ range (Yaron et al. 2005).The time required to accrete a critical envelope mass (i.e. the time between successive nova eruptions) on ∼ 1.4 ⊙ WDs can be as short as 1 yr (Darnley et al. 2014), to longer than 1 Myr on low mass WDs (Yaron et al. 2005).Because the critical envelope masses are all much smaller than the masses of CV donors, all CVs, including those never having been observed to undergo a nova eruption, will eventually undergo (many) nova eruptions (see also Ford 1978).Every CV not currently undergoing a nova eruption is detached (Hillman et al. 2020) or re-accumulating an envelope onto its WD in preparation for the next nova event.
Large / (typically > 10 −8  ⊙ yr −1 ) are seen in CVs with orbital periods  ∼ 5 − 10 hr; these are the so-called nova-like binaries.In the century before an eruption, and for decades afterwards, most novae are indistinguishable from nova-like binaries, exhibiting high / (Collazzi et al. 2009).The accretion disks of CVs exhibiting / slightly less than those of nova-like variables are unstable, and undergo transitions from a hot to a cold state (Lasota 2001;Dubus et al. 2018).This leads to the rapid, episodic dumping of much of the accretion disk onto the WD-a dwarf nova (DN) outburst (Osaki 1974).The liberation of gravitational potential energy by the sudden accretion of most of the mass in an accretion disk (∼ 10 −10  ⊙ ) can liberate a few  ⊙ for about a week.
The subclass of CVs with / close to the nova-like-DN boundary are named after the prototype Z Camelopardalis (Z Cam).When Z Cam's / from the secondary star into the accretion disk surrounding the WD is too large to produce dwarf nova outbursts, the system resembles a nova-like variable.Z Cam returns to undergoing DN outbursts when / declines below the critical level (Osaki 1974;Buat-Ménard et al. 2001;Simonsen et al. 2014).
The hibernation scenario of CVs (Shara et al. 1986;Prialnik & Shara 1986;Livio & Shara 1987;Kovetz et al. 1988;Hillman et al. 2020) predicts that the mass lost during a nova eruption should lead to an increased separation of the WD from its donor, and a concomitant decrease in / several centuries after an eruption, once irradiation from the post-nova, cooling WD fades.If this prediction is correct then some old novae, initially observed to be nova-like variables, should begin to undergo dwarf nova eruptions as / declines.The best-known example is GK Persei (nova Per 1901) which began DN outbursts 65 years after its nova outburst.An even more extreme case is V1017 Cen, which underwent a DN event that may have triggered WD surface nuclear burning just 14 years after its classical nova outburst in 2005 (Aydi et al. 2022).Kato (2022) lists 10 old novae that are currently undergoing DN eruptions.
The hibernation scenario also predicts that some DN should still be surrounded by the ejecta of their last nova eruption.The nova shell of the prototypical dwarf nova Z Cam was the first discovered (Shara et al. 2007), and it remains the largest and oldest such shell known (Shara et al. 2012).As of 2023, three dwarf novae are known to be surrounded by the ejecta of classical novae (Miszalski et al. 2016;Shara et al. 2017a,b).
A comparison of narrow-band H  + [N ii] images of the Z Cam shell, taken with the Mayall 4-m telecope in 2007 and 2010, constrained but did not detect expansion of the nova shell (Shara et al. 2012), with an upper limit of 0.17 arcsec yr −1 .Including the important effect of deceleration as the ejecta sweeps up interstellar matter in its snowplow phase, the observed upper limit on shell expansion constrains the nova event to have occurred at least 1300 years ago.This was the first strong test of the prediction of nova thermonuclear runaway theory that the inter-outburst times of most classical novae are longer than 1000 years.
At its Gaia-determined distance of just 214 pc (Bailer-Jones et al. 2021), Z Cam must achieve a peak visual magnitude in the range ∼ +1 to −3 when it undergoes a nova eruption, and it is brighter than 2nd magnitude for 2 to 3 weeks.The intriguing suggestions that Z Cam was a bright nova recorded by Chinese imperial astrologers in October and November 77 BCE (Johansson 2007) or more recent transients in Ursa Major (Warner 2015) is consistent with our expansion upper limit (see Hoffmann (2019) for a critique of the latter suggestion).We predicted (Shara et al. 2012) that if Z Cam is indeed the nova of 77 BCE, its ejecta should currently be expanding at ∼ 85 km s −1 , corresponding to 110 mas yr −1 .Since 15 years have passed since our 2007 epoch of observations, and a shell expansion of 15 yr × 110 mas yr −1 = 1.65 arcsec should be straightforward to measure, we undertook a program of narrow-band imaging of Z Cam as part of the science verification of the new Condor Array Telescope (Lanzetta et al. 2023), hereafter Paper I.
In Section 2 we briefly describe the Condor telescope, emphasising its features optimised for extremely low surface brightness imaging.We also describe the two epochs' images that we used to search for shell expansion.In Section 3 we present the Condor images to demonstrate the telescope's remarkable low surface brightness sensitivity.The methodology used to compare the two epochs' images is described in Section 4, and the measured angular and spatial expansion velocities are presented.Our results are briefly summarized in Section 5.

OBSERVATIONS
Because this is one of the first papers to present scientific results based on the Condor Array Telescope, we briefly describe Condor and its capabilities before describing the datasets used in our analyses.

The Condor Array Telescope
Condor is an "array telescope" that consists of six apochromatic refracting telescopes of objective diameter 180 mm, each equipped with a large-format (9576 × 6388 pixels), very low read-noise (1.2 e −1 ), very rapid read-time (< 1 s) CMOS camera.It is of particular relevance here that Condor is equipped with a set of He ii 468.6, [O iii] 500.7,He i 587.5, H  656.3, [N ii] 658.4,and [S ii] 671.6 narrow-band filters (one per telescope).Condor has no interferometric capability; its six co-aligned refractors enable the acquisition of images of a single field through up to six different filters simultaneously.Condor is located at the Dark Skies New Mexico observatory new Animas, New Mexico.Details of the telescope are described in Lanzetta et al. (2023).

The Data
We observed Z Cam using the Mosaic CCD camera on the Kitt Peak National Observatory (KPNO) Mayall 4-m telescope in 2007 and using Condor in 2021.
The Kitt Peak National Observatory (KPNO) Mayall 4-m telescope with the Mosaic CCD camera at prime focus was used to observe Z Cam on the nights of 2007 January 23 and 24.The Mosaic CCD camera incorporates eight 2048 × 4096 pixel 2 SITe thinned CCDs, which were used with 2 × 2 on chip binning to yield an effective 8192 × 8192 pixel 2 format.At prime focus, the camera obtains an image scale of 0.26 arcsec per binned pixel, which results in a 36 × 36 arcmin 2 field of view.We obtained 106 images of total exposure time 45,300 s through the KPNO k1009 narrow-band filter.In the f/3.1 beam of the Mayall telescope the filter's wavelength centroid is 656.3 nm and its FWHM 8.1 nm, which is sensitive to the sum of H  656.3 nm and [N ii] 658.4 nm.The telescope was dithered between exposures.
The images were bias subtracted and flat fielded using standard procedures, and the individual images were aligned using daophot (Stetson 1987).The images were then stitched together using mon-tage2, which is a mosaicking program within daophot.See (Shara et al. 2012) for further details, and the images of Z Cam.
We used Condor with its complement of narrow-band filters to observe Z Cam in November 2021.The Condor CMOS cameras obtain an image scale of 0.85 arcsec per pixel, which results in the same 2.3 × 1.5 deg 2 (134 × 91 arcmin 2 ) field of view for each of the six telescopes.Each exposure was 600 s in length, and we obtained a total exposure time of 63,600 s through each filter with one of the six telescopes.The array was dithered by a random offset of 15 arcmin between exposures.Images of the dusk and dawn twilight sky were obtained every night, and bias observations were obtained at the end of every night.
The observations were processed through the Condor data pipeline (Papers I and II), with steps involving bias subtraction, flat fielding and background subtraction, astrometric calibration, and photometric calibration.It is of particular relevance here that the astrometric calibrations exhibit systematic differences between the transformed pixel and celestial coordinates of < 0.1 arcsec.

EXTENDED NEBULOSITY
The six narrow-band images of Z Cam are shown as a mosaic in Fig. 1.The [S ii], co-added H  and [N ii], and [O iii] images were colorized red, green, and blue, respectively, summed as a single image, and displayed in Fig. 2. Despite the fact that the light-collecting area of Condor is ∼ 80 times smaller than that of the Mayall 4-m telescope, Condor clearly images fainter features (see Fig. 2) than the 4-m telescope (see figs. 1 and 3 of Shara et al. 2012) in an exposure time just 1.4 times longer, through filters of comparable transmission efficiency and FWHM.Both the KPNO and Condor images are sky noise limited, so read noise cannot account for the greater Condor sensitivity.It is the combination of Condor's exceptionally clean point spread function and careful image flattening and sky subtraction in the Condor data reduction pipeline that are responsible for the increased sensitivity (Lanzetta et al. 2023).
In addition to the bright arc to the SW of Z Cam, linear features to the NE, SE, and SW of Z Cam seen in Shara et al. (2012) are clearly visible in the H  and [N ii] Condor images (see Fig. 1).Notably absent in the Mayall 4-m telescope images is any nebulosity North or NW of Z Cam.A faint arc, brightest in [N ii] and H  (green in Fig. 2), is now clearly seen in the Condor images to the NE and North of Z Cam.The surface brightness of this arc is ∼ 26.2 mag arcsec −2 .This is at least 1 mag arcsec −2 more sensitive than the recent search of the environs of 47 CVs (Sahman & Dhillon 2022), which did not locate any new nova shells.
Together with the previously detected nebulosities, this faint arc forms a circular shell ∼ 40 arcmin in diameter, though notably not centered on Z Cam itself (circled in red in Fig. 2).Radial velocities of the gas throughout the circular shell will be needed to establish with certainty that the northern arc is related to the brighter Z Cam nebulosities, and whether the entire shell is a single entity.
An even fainter arc is seen in the H  and [N ii] images a few arcmin north and east of the newly detected NE arc (see Fig. 3).If the mass of the WD in Z Cam is ∼ 1 M ⊙ (Shafter 1983), and its average accretion rate is ∼ 10 −9  ⊙ yr −1 (Hartley et al. 2005), then we expect nova eruptions separated by ∼16-46 kyr (Yaron et al. 2005).The observed 21.1 mas yr −1 proper motion of Z Cam (Bailer-Jones et al. 2021) demonstrates that the time between successive nova eruptions is significantly less than the time required (∼ 100 kyr) for Z Cam to cross its own shell.It is thus tempting to ascribe the very faint outer arc to a previous eruption of Z Cam.
In Fig. 4 we test this hypothesis by showing the entire Condor fieldof-view (2.5 × 1.5 degrees) of the co-added H  + [N ii] images, after 10×10 pixel block averaging.A nearly circular, shell-like structure is seen to encircle the now fully-visible inner shell.No trace of this outer shell is seen in the (equally deep) He ii, He i, [O iii] or [S ii] images, or their sum, after 10 × 10 pixel binning.This argues strongly against Galactic cirrus as the source of these very low surface brightness features.
In Fig. 5 we compare the co-added H  + [N ii] image with the Galex FUV and NUV archival images.The NE arc is clearly seen in the NUV image, but not the FUV image.The outer shell is not seen in either of the NUV or FUV images.
In Fig. 6 we show azimuthally-averaged radial flux profiles for each of our six narrow-band images, normalized by their median values.The profiles are given by the weighted mean flux per 100 pixel-thick radial bin taken from the approximate center of the outer, nearly spherical shell.The standard error on the mean is shown by the shaded regions, which are only slightly thicker than the lines themselves.The inner shell, with radius ∼ 15 arcmin is clearly seen in H , [N ii] and [S ii].The outer shell, with radius ∼40 arcmin, is strongest in H .
Have we detected the first instance of concentric nova shells due to successive nova eruptions?Nova theory predicts that such nested shells must exist, and be easiest to detect surrounding massive white dwarfs which erupt on (relatively short, for novae) timescales of order just 1-10 kyr.Though many other such multiple shells should exist, no other shell has ever been imaged to the depth of our images of Z Cam, explaining why it is now the first candidate multiple-shell nova.Deeper imaging and/or spectroscopic characterization of the outer shell will be needed to unambiguously determine its nature.

MEASURING THE SHELL EXPANSION
The brightest and sharpest-edged feature in the complex of nebulosity surrounding Z Cam is the arc ∼ 15 arcmin SW of the star (Shara et al. 2012).The 3-year imaging baseline of KPNO images only allowed us to place an upper limit on the expansion of Z Cam's ejecta.The SE quasi-linear feature is too faint in the KPNO and Condor images to yield any expansion estimate.
To quantify the expansion of the SW arc, we employed two methods of determining the spatial shift of nebulosity between the two epochs (2007 and 2021): sum of squared residuals (SSR) minimization and cross-correlation peak measurement.Their means yield our estimate of the shell's velocity, time since eruption, and uncertainties in both these quantities.

Image alignment, reprojection, and calibration
We retrieved a WCS solution for the 2007 Kitt Peak H  + [N ii] image from astrometry.net (Lang et al. 2010).Using the adaptive resampling algorithm of reproject (Robitaille et al. 2020), we projected the Kitt Peak image onto the pixel grid of the sum of the Condor H  + [N ii] images, matching the latter's 0.85 arcsec/pixel resolution.To ensure an accurate alignment between the images, we matched sources in the Condor and reprojected Kitt Peak images.With the matched source lists, we used astroalign (Beroiz et al. 2020) to determine the affine transformation that matched the Kitt Peak centroids to the Condor centroids, which was then applied to the Kitt Peak image.As noted above, this astrometric matching is accurate to better than 0.1 arcsec across the entire Condor image.
The SW arc (see Fig. 7) was analyzed using a cropped segment of the aligned field.Within this subfield, we subtracted background values for both images and, for SSR minimization, we additionally multiplied the Condor image by a scaling factor determined by the flux of the arc's nebulosity.

Sum of squared residuals minimization
It is helpful to refer to Fig. 7 in the following.To determine the position of the SW arc, we produced vertical flux profiles spanning ±15 pixels of the brightest section of the leading edge in both the Kitt Peak and Condor images.The leading edge was defined via the following procedure.In the Kitt Peak image, we drew by eye The image has been downsampled by a factor of 32 with averaged pixel blocks and is displayed inverted with linear scaling.The now-clear "inner" shell encircling Z Cam ("inner shell", solid yellow outline) is surrounded by an outer, roughly circular nebulosity (though flattened on the southern rim) which is approximately 1.5 degrees in diameter ("outer shell", between dashed red rings).This larger feature may be the remnant of the nova eruption previous to the one that produced the now fully-seen, inner shell.The nebulosity to the north and east of the inner shell is also visible ("NE arc", solid light blue outline).
MNRAS 000, 1-13 ( 2023) This gave a set of points that followed the arc with vertical scatter of ∼ 2 pixels.We fitted these points with a 10th-degree polynomial and took the result to define the leading edge of the arc.We limited the horizontal range of the edge (i.e. in the Right Ascension direction) to the brightest continuous region for which the process just described gave good polynomial fits.This gave a SW arc edge of length 533 pixels, or 7.64 arcmin (see Fig. 7).
For each pixel along the above-defined leading edge of the arc, we took vertical cuts (i.e. in the Declination direction) spanning ±15 pixels above and below that point.Measuring the flux values at these points for a given image produced 533 vertical profiles each 30 pixels in length, which were then averaged to retrieve a single "net" vertical flux profile for the entire arc (see subfigures on the bottom and right of Fig. 8).
The SSR, calculated as the sum of the the squared differences between the two images' net vertical flux profiles, quantifies the alignment of the feature between the epochs.By shifting the Kitt Peak image in 0.1 pixel increments up to maximum southward and westward shifts of 2.0 pixels and 3.0 pixels and remeasuring the net vertical flux profiles for each shift, we determined the shift for which the SSR is minimized.This coincides with the shift that best aligns the arc in the Kitt Peak image to the arc in the Condor image.This shift gives the net motion of the arc over the nearly 15 year period between the two images and thus the expansion velocity.Image shifting was computed using ndimage's (Gommers et al. 2022) shift function with fifth order spline interpolation.The SSR for the selected shifts is coded by color in the central subfigure of Fig. 8.
We determined the velocity of the arc by applying this method with bootstrapping (Simpson & Mayer-Hasselwander 1986), randomly selecting 533 points with replacement along the leading edge of length 533 pixels, recording the shift that minimized the SSR, and repeating this process 10,000 times.For some bootstrapped samples, the shift that produced the minimum SSR was found to lie outside of the tested shifts, giving a minimum point on the edge of the search range.We discarded these tests.The mean and standard deviation of the resulting distribution of shifts were taken to represent the most likely shift and a conservative estimate of the uncertainty for this method.

Cross-correlation peak measurement
Displacement between two images that share a feature can additionally be quantified via cross-correlation.The product of the Fourier transform of an image with the complex conjugate of the Fourier transform of the displaced image yields the Fourier transform of the cross-correlation.Thus the inverse Fourier transform of this product gives the cross-correlation plane, across which coordinates correspond to shifts of one image relative to the other.The position of peak cross-correlation amplitude is representative of the translation that best aligns the features in the two images.Thus we take this point to be the maximum-likelihood shift of the arc between the two epochs.To reduce contaminating the cross-correlation function with sources other than the relevant nebulosity, we masked out regions farther than approximately 10 pixels away from the arc as well as any nearby or overlapping extraneous sources.We calculated the crosscorrelation plane using the scikit-image (van der Walt et al. 2014) masked normalized cross-correlation function (Padfield 2012).To measure cross-correlation on a subpixel scale, we rescaled the images before cross-correlation by a factor of 5.The cross-correlation plane calculated from the images upscaled by a factor of 10 to a 0.1 pixel resolution is coded by color in Fig. 9.
We employed bootstrapping to determine the distribution of likely cross-correlation maxima.This method is advantageous for retrieving measurement statistics on the cross-correlation of distributions with unknown errors without added assumptions (Sciacchitano 2019).This approach has seen astrophysical use in, for example, multiwavelength blazar light curve studies (e.g.Max-Moerbeck et al. 2014).As the flux uncertainties of the Kitt Peak image are not known, we limit our implementation of the method (described in Peterson et al. 1998) to the "random subset selection" procedure.For each of 10,000 iterations, we randomly selected with replacement a number of pixels equal to the total size of the image and removed unselected pixels from the cross-correlation calculation by masking them.Unlike SSR minimisation, in which multiply-selected vertical profiles can proportionally count towards the resulting minimum, multiply-selected pixels in images cannot easily be weighted in their contribution to the cross-correlation.Therefore we instead counted these pixels only once.As this reduces the size of each bootstrapped sample by a factor of approximately 1/, we expect the uncertainties to be somewhat larger than a hypothetical fully-sampled bootstrapping method (Peterson et al. 1998).As with the SSR minimisation method, the mean and standard deviation of the distribution of crosscorrelation maxima were interpreted as the shift and a conservative estimate of its uncertainty.

Determining the motion of the SW arc
Of the 10,000 bootstrapped samples taken for the SSR minimisation for the SW arc, 7,653 had valid minima within the range of tested shifts.The covariance of the resulting points produced the error ellipse shown in Fig. 8. Z Cam is the Gaia DR3 source 1123169888190445568, located at a distance of 214 pc (Bailer-Jones et al. 2021).These data suggest the SW arc is expanding to the west and south at 65 ± 38 km s −1 and 47 ± 19 km s −1 , respectively.This corresponds to a net motion of 80 ± 42 km s −1 at 36 • south of west.
The mean cross-correlation peak of the SW arc appears at a west and south translation corresponding to 88 ± 29 km s −1 and 18 ± 11 km s −1 , respectively, or a net motion of 90 ± 32 km s −1 at 12 • south of west (see Fig. 9).
The combined result of these two methods, a shift and formal uncertainty as the means of the two resulting shifts and uncertainties, gives the expansion of the SW arc to be a net motion of 83±37 km s −1 at 23 • south of west.
Taking the time since the eruption of Z Cam given by ejecta in the snowplow phase to be  = /4 (Oort 1951), where the distance from Z Cam to the SW arc is  = 0.9 pc, we find  = 2772 years via SSR minimisation, and  = 2473 years via the cross-correlation peak.The modest difference between these two ages (just 299 years, or 11% of the age we adopt below) is encouraging, and suggests that the formal errors associated with each method above may be significant overestimates.We nonetheless adopt the mean value and mean uncertainty for our adopted age of the Z Cam shell: 2672 −817 +2102 yr.

CONCLUSIONS
Deep narrow-band images of the ejecta of Z Cam, taken with the Condor Array Telescope, have revealed the nearly complete, circular shell, ∼ 30 arcmin in diameter, of this old nova/dwarf nova.An ∼80 arcmin diameter, fainter concentric shell is also detected.If the faint outer shell is confirmed, Z Cam would be the first instance of a pair of faint, concentric shells surrounding an old nova, representing the two most recent eruptions of the system's white dwarf.
We have measured the expansion velocity of the brightest ejecta to be 83 ± 37 km s −1 , and its age to be 2672 −817 +2102 yr.The suggestion that the "new" star first seen by Chinese Imperial astrologers in 77 BCE (Ho 1962) was a nova eruption of Z Cam (Johansson 2007)

Figure 1 .
Figure 1.Clockwise from the top right, Condor images of the area surrounding Z Cam in [O iii], H , [S ii], [N ii], He i, and He ii.The total exposure time in each of the six narrow-band filters was 63,600 sec.Each image was taken with one of the six 180 mm apochromatic refracting telescopes that constitute the Condor Array (Lanzetta et al. 2023).Images have been smoothed via a Gaussian kernel ( = 2.5 pixels) and are displayed with linear scaling.

Figure 2 .
Figure 2. Top: Area surrounding Z Cam in Condor [S ii], H  + [N ii], and [O iii] in red, green, and blue, respectively.Z Cam is circled in red.The images have been smoothed via a Gaussian kernel ( = 2.5 pixels) and are displayed with asinh scaling.The current direction of motion (from Gaia, Bailer-Jones et al. 2021) of Z Cam is shown with a dashed cyan arrow, and the solid yellow arrow points to the location of Z Cam 20,000 years ago.Bottom: The area of the region containing the SW arc, SE linear feature, and NE linear feature (denoted by the white box) is shown enlarged.

Figure 3 .
Figure 3. Area surrounding Z Cam in co-added Condor [N ii] and H . The image has been smoothed via a Gaussian kernel ( = 2.5 pixels) and is displayed with linear scaling.Note the faint nebulosity several arcmin north and east of the NE arc.

Figure 4 .
Figure 4.The entire Condor field in co-added [N ii] and H  (top) and the same with annotations of notable features overlayed (bottom).The image has been downsampled by a factor of 32 with averaged pixel blocks and is displayed inverted with linear scaling.The now-clear "inner" shell encircling Z Cam ("inner shell", solid yellow outline) is surrounded by an outer, roughly circular nebulosity (though flattened on the southern rim) which is approximately 1.5 degrees in diameter ("outer shell", between dashed red rings).This larger feature may be the remnant of the nova eruption previous to the one that produced the now fully-seen, inner shell.The nebulosity to the north and east of the inner shell is also visible ("NE arc", solid light blue outline).MNRAS 000, 1-13 (2023)

Figure 5 .
Figure 5.Comparison of GALEX FUV, GALEX NUV, and Condor H  + [N ii] images in the region surrounding Z Cam.The Condor image has been downsampled by a factor of 16 with averaged pixel blocks and is displayed inverted with linear scaling.

Figure 6 .
Figure 6.Azimuthally-averaged radial flux profiles for each band, normalized by their median values.The profiles are given by the weighted mean flux per 100 pixel-thick radial bin taken from the approximate center of the outer arc.Standard error on the mean is shown by the shaded regions, which are only slightly thicker than the lines themselves.

Figure 7 .
Figure 7. H  + [N ii] Kitt Peak image in the region of the SW arc.The leading edge determined by 10th-degree polynomial fitting is shown in red.
is consistent with our deduced expansion age of Z Cam's ejecta.(We also cannot rule out all the suggestions of Warner 2015 of more recent transients, but see also Hoffmann 2019; Hoffmann & Vogt 2020 and Neuhäuser & Neuhäuser 2021).The age uncertainty of the last eruption of Z Cam remains large, so the identity of the underlying star of the transient of 77 BCE is far from settled.ACKNOWLEDGEMENTS MMS, KML and JTG gratefully acknowledge the support of NSF award 2108234.This material is based upon work supported by the National Science Foundation under Grants 1910001, 2107954, and 2108234.Based on observations at Kitt Peak National Observatory, NSF's National Optical Astronomy Observatories (NOAO Prop.2006B-0417; PI:Shara), which is operated by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation.This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia),processed by the Gaia Data Processing and Analysis Consortium (DPAC, https: //www.cosmos.esa.int/web/gaia/dpac/consortium).Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement.MMS acknowledges the longtime interest of Göran Johansson, who first suggested that Z Cam was the nova of 77 BCE, and correspondence with S. Hoffmann and R. Neühauser about historical novae.Facilities: KPNO, Condor Software: astroalign (Beroiz et al. 2020), astrometry.net(Lang et al. 2010), astropy (Astropy Collaboration et al. 2013, 2018, 2022), daophot (Stetson 1987), ndimage (Gommers et al.

Figure 8 .
Figure 8. Sum of squared residuals between the Condor and Kitt Peak vertical flux profiles for different pixel shifts of the SW arc, with color data corresponding to the fully sampled leading edge.The point X corresponds to the mean minimum point from 10,000 bootstrapped samples of the arc.The direction of the velocity vector is shown by the black arrow.Vertical profiles for the mean point as well as the vertices of the corresponding error ellipse (A, B, C, and D) are shown below and to the right.