Revealing the Morpho-Kinematics of NGC 2371, a Planetary Nebula with a [WR] Central Star

1. Introduction

Planetary nebulae (PNe) mark one of the final stages in the evolution of low- and intermediate-mass stars ($1\lesssim M_{\mathrm{i}}/{\mathrm{M}}_{\odot }\lesssim 8$). After leaving the main sequence phase, these stars undergo intense mass-loss episodes, particularly as they enter the asymptotic giant branch (AGB) phase e.g., [1,2] and references therein. This process deposits most of the star’s mass in a circumstellar envelope through a slow and dense wind while shedding its outer layers. As a result, the exposed core increases its effective temperature ($T_{\mathrm{eff}}$) to values exceeding 20 kK. Subsequently, the post-AGB star develops a vigorous stellar wind and UV photon flux, which collectively compress, photoionise, and heat the circumstellar material, ultimately giving rise to the PN [3].

However, the scenario described above must be addressed despite its illustrative nature. There is growing evidence that most PNe host binary progenitor systems [4], and consequently, their evolution may have played a significant role in shaping the variety of morphologies observed in PNe see, e.g., [5]. One particularly promising scenario is the formation of PNe after the binary system experiences a common envelope evolution [6]. In this scenario, the companion star shares the envelope of an inflated AGB or red giant star, contributing to the shaping of the mass loss that produces toroidal structures and jet-like ejections, directly impacting the formation of the PN e.g., [7,8,9,10] and references therein.

One of the best tests of the common envelope scenario is to study the kinematics of proto-PNe and PNe see, e.g., [11]. This becomes particularly pertinent when considering complex structures and multiple ejection axes, which are evidenced by deep imaging and spectroscopy, bringing further complexity to the analysis see e.g., [12,13].

In this paper, we present a detailed characterisation of the bipolar PN NGC 2371. In Figure 1, we show a composite [O iii]$\lambda 5007$ image of this PN. A reverse greyscale image is saturated to enhance the faint outer structures of the nebula, while a standard greyscale image is superimposed to highlight the inner regions. The main visible structures that will be used later are labelled. Details of this observation are mentioned in the corresponding section. With respect to NGC 2371, various studies encompassing radio, infrared, optical, and UV spectra have been conducted to address the physical properties of this PN and its hydrogen-deficient [Wolf-Rayet]([WR])-type central star (CS) see, e.g., [14,15,16,17]. According to the analysis of Gaia data presented by [18], the distance to NGC 2371 is $d=1.66\pm 0.13$ kpc.

To our knowledge, Sabbadin et al. [19] conducted the initial effort to characterise the morpho-kinematics of this PN using various emission lines. These authors proposed that NGC 2371 has the shape of a three-axis ellipsoidal cavity with a toroidal central region and two external polar caps. However, it should be noted that the deep imaging of NGC 2371 has revealed that the so-called polar caps are part of the more extended bipolar structure, as illustrated in Figure 1.

More recently, Gómez-González et al. [20] presented one of the most comprehensive studies of this PN, in which they not only modelled the stellar atmosphere properties of the H-deficient [Wolf-Rayet]-type CS, but also obtained long-slit spectra from various morphological features to study its physical properties and abundances. These observations suggest that NGC 2371 exhibits a double structure that protrudes from a barrel-shaped central region [20].

Here, we present high-resolution spectra obtained at various epochs to meticulously study the kinematics of NGC 2371. Although our group has previously discussed preliminary results from some of these spectra, see [21,22], our current investigation is motivated by the recent work of Gómez-González et al. [20], prompting a detailed discussion of the formation scenario of NGC 2371 through a morpho-kinematic model. This paper is organised as follows. Section 2 outlines the details of our observations. Section 3 and Section 4 present our findings and their interpretation. Finally, we offer our concluding remarks in Section 5.

2. Observations

The [O iii]$\lambda 5007$ image in Figure 1 is used to illustrate the basic morphological structure of the nebula and to show the locations of the additional crossing slits (see Appendix A). This image was obtained on March 9, 2004, with the Arcadio Poveda 2.1 m telescope ($f/7.5$) at the Observatorio Astronómico Nacional in San Pedro Mártir, Mexico (OAN-SPM), using the “Rueda Italiana” filter wheel. A SITe CCD ($1024\times 1024$, 24 $\mu$m pix⁻¹, plate scale 0.3 arcsec) was used in unbinned mode. The interference filter is centred at $\lambda 5015$ Å with a bandwidth of $\Delta \lambda =51$ Å. The total exposure time was 600 s. Standard bias and flat-field correction, as well as cosmic-ray removal, were applied. The seeing during the observations was about 1.8 arcsec. This image is used only for morphological purposes, and no attempt was made to perform absolute flux calibration.

High-dispersion optical spectra were acquired using the Manchester Echelle Spectrometer MES [23] and the Arcadio Poveda Telescope (2.1 m; $f/7.5$) at the OAN-SPM during multiple observing runs (see Table 1). The average seeing during observations was approximately $1.7\pm 0.5$ arcsec. A CCD binning mode of $2\times 2$ was employed, utilising a SITe3 ($1024\times 1024$, 24 $\mu$m pix⁻¹, plate scale 0.62 arcsec) for runs A to E, and an E2V-4240 ($2048\times 2048$, 13.5 $\mu$m pix⁻¹, plate scale 0.35 arcsec) for runs F and G. The slit dimensions were set to 6.5 arcsec in length and 2 arcsec ($150\mu$m) in width. Due to the absence of cross-dispersion, a bandwidth filter of $\Delta \lambda =60$ Å was used to isolate the 114^th order for observing the [O iii] $\lambda$5007 emission line. Furthermore, a $\Delta \lambda =90$ Å bandwidth filter was used to isolate the 87^th order, covering the spectral range containing the H$\alpha$ and [N ii] $\lambda$6583 emission lines, but only for slits E1, E2, F2, and F3.

Wavelength calibration was conducted using a ThAr arc lamp, achieving an accuracy of $\pm 2$ km s⁻¹. The Full Width at Half Maximum (FWHM) of the arc lamp emission lines was estimated to be approximately $\approx 12$ km s⁻¹. The spectral scales for the SITe3 were 0.077 Å ([O iii]) and 0.100 Å (H$\alpha$+[N ii]), and for the E2V were 0.043 Å ([O iii]) and 0.057 Å (H$\alpha$+[N ii]).

Table 1 summarises the observations made by positioning the slit across various regions of the nebula. The data were calibrated using the NOIRLab 2.18 version of IRAF, which provides processing routines for long-slit spectroscopy [24,25,26]. Examples of slit paths of runs F and G are shown in Figures 2 and 4, while the others are shown in Appendix A.

3. Results

As can be seen in Figure 1, NGC 2371 displays an overall bipolar morphology. The central innermost region is dominated by a very bright, barrel-like structure. Embedded within this barrel are two compact knots of enhanced emission, located on opposite sides of the CS but not exactly collinear with it. Further out along this axis, we identify two fainter caps, and at the largest projected distances, two distorted outer lobes. The image has been processed to simultaneously enhance the bright core emission and the faint structures that outline the lobes. All these components, together with the CS, are labelled in the figure, and this nomenclature is adopted in the following sections.

With respect to kinematics, we focus on analysing the slits obtained in the F and G runs, although the resulting model was successfully tested on the other available slits (Appendix A).

In Figure 2, we show the position angle locations for run F. On the other hand, in Figure 3, the Position-Velocity (PV) diagrams (or maps) of the F1, F3, and F4 slits in the [O iii]$\lambda 5007$Å emission line (top row) are presented, alongside the F2 and F3 slits in the emission lines H$\alpha$ and [N ii]$\lambda 6583$Å (bottom row). Each PV diagram is labelled with the slit designation in the upper right corner, the corresponding filter name in the bottom left corner, and the slit orientation in the upper left corner.

We used the line splitting observed over the central star in spectra F1, F3, and F4 to estimate the systemic radial velocity, taking the midpoint between the blue- and red-shifted components in each case. From these measurements, we obtain mean systemic velocities of $v_{\mathrm{Hel}}^{\mathrm{sys}}=17\pm 2$ km s⁻¹ and $v_{\mathrm{LSR}}^{\mathrm{sys}}=8\pm 2$ km s⁻¹, in the heliocentric and Local Standard of Rest (LSR) frames, respectively.

The PV diagram for F1 reveals an expanding shell characterised by several internal structures, including notable asymmetry, although an elliptical structure can be fitted to the data. The angular size of this structure is $\simeq 50\pm 2$ arcsec, corresponding to the bright central region seen in Figure 2. The eastern segment of this structure exhibits a redshift, whereas the western segment shows a blueshift. Figure 2 also shows the slit intersecting the weaker regions of the lobe-like structures at the extremities, and the corresponding PV map reveals faint evidence of their expansion, with no significant inclination observed.

The more extreme velocities in this PV diagram are $-82$ and $+92\pm 2$ km s⁻¹, with respect to the systemic velocity. Furthermore, two prominent knots are discernible, each roughly $\sim 10$ arcsec in width and exhibiting a dispersion of $\sim 50$ km s⁻¹. Their centroids are symmetrically positioned at $\simeq 12$ arcsec from the central star and $\approx 30$ km s⁻¹ apart from each other. The orientation of these broad knots is opposite to that of the elliptical structure, with the eastern knot blueshifted and the western knot redshifted.

The PV map of slit F3, acquired at a location corresponding to the minor axis of the nebula, reveals a closed ellipsoidal structure. Despite the presence of internal features, it exhibits a higher degree of symmetry compared to F1. The angular extent of this region is approximately $\simeq 47$ arcsec, with the spectral line splitting at the position of the CS measuring $130\pm 2$ km s⁻¹. The NE and SW extremities of the ellipse are separated by a mere $\sim 10$ km s⁻¹, which are more pronounced on the PV map of F1.

The final slit of this run, observed in [O iii], is F4. The PV map corresponding to this slit reveals a particular point-symmetry. Starting from the position of the CS and progressing towards the NW, the centroid of the spectral line splitting shifts towards the red end of the spectrum. The intensity of the bluer component diminishes until it nearly vanishes at approximately 20 arcsec, whereas the redder component persists up to nearly 36 arcsec, stabilising its shift around 112 km s⁻¹. In contrast, in the SE direction, the centroid of the splitting appears to shift toward the blue, stabilising after $\approx 15$ arcsec, with the bluer expansion near $-100$ km s⁻¹, thus contributing to this apparent symmetry. Close to the ends of the spatial axis, there are two enhanced intensity emission regions, symmetrically located at $\simeq \pm 65$ arcsec from the CS, with corresponding velocities of $\simeq \pm 38$ km s⁻¹, respectively. These emission regions correspond to the zones in which the NW and SE caps are crossed by the F4 slit.

The bottom panels in Figure 3 present the PV diagrams for H$\alpha$ and [N ii]$\lambda 6583$ Å corresponding to slits F2 (left) and F3 (right). For slit F2, substantial emission in H$\alpha$, and even in He ii$\lambda 6560$, is detected within the brightest peak nebular emission. This observation suggests that both regions, spanning approximately 25 arcsec, exhibit very high excitation levels. However, within these regions, a pair of knots exhibits intense [N ii] brightness. The morphology of these regions in the direct images (Figure 1 and Figure 2) is radially elongated. Although the [N ii] emission and the morphology of these regions are indicative of collimated outflows (jets), it is noteworthy that these features are not perfectly aligned (separated $\approx 7$°). In [N ii], the maxima of these jet-like structures are located $\simeq 18\pm 2$ arcsec from the central star and exhibit velocities in opposite directions, as anticipated in collimated bipolar ejections. Despite their relatively low radial velocities, $+25\pm 2$ km s⁻¹(NE) and $-5\pm 2$ km s⁻¹(SW), it is important to consider that these values are projections in the plane of the sky. However, the notable aspect is the clear disparity in their magnitudes.

Conversely, the H$\alpha$ emission in the PV map corresponding to F3 closely resembles the [O iii] map. Similar to the case of F2, there is also a distinct emission of He ii$\lambda 6560$Å. The [N ii] emission is considerably more marginal, exhibiting a relatively broad region towards the southwest. Additionally, an arc-like structure is discernible in the blue region peaking at $-73\pm 2$ km s⁻¹ and extending from the central star to approximately 12 arcsec towards the northeast. Precisely at that location, a highly compact emission is detected at $\simeq 30$ km s⁻¹. Progressing towards the northeast from that point, extending nearly up to $\simeq 27\pm 2$ arcsec, a faintly observed splitting is evident, which tends to converge as it recedes from the central star, mirroring the shape and dimensions of the northeastern segment observed in [O iii]. This slit does not intersect the intense regions in [N ii] that define the luminous central zone of this nebula; however, as depicted in Figure 3 (left), numerous knots, smaller in size and intensity than the “jets”, are observed throughout that area, predominantly emitting in [N ii].

Figure 4. Colour-composite image of NGC 2371. The red, green, and blue channels correspond to the [N ii], H$\alpha$, and [O iii] optical images presented in Gómez-González et al. [20] (Figure 1). Slit positions are overlaid on the image for the run G. Refer to Table 1 for details.

Figure 4. Colour-composite image of NGC 2371. The red, green, and blue channels correspond to the [N ii], H$\alpha$, and [O iii] optical images presented in Gómez-González et al. [20] (Figure 1). Slit positions are overlaid on the image for the run G. Refer to Table 1 for details.

Figure 5 presents the PV maps corresponding to slits G1, G2, G3, and G4. These slits were observed at a $\mathrm{PA}=+37$°, crossing regions outside the main shell, which manifest as distinct lobes or large bubbles. The PV maps for G1 and G2 exhibit a predominance of the blue component, whereas for G3 and G4, the red component is more pronounced. Specifically, for G1, the emission spans approximately 66 arcsec, with expansion velocities ranging from 35 km s⁻¹ (NE) to 63 km s⁻¹ (SW), with the red velocity component being notably marginal.

The counterpart of G1 on the NW lobe is slit G4. A comparable structure is discerned, with a slight expansion in the SW region, manifesting a velocity of $v_{\exp }=34\pm 2$ km s⁻¹, which escalates towards the NE to approximately 58 km s⁻¹. It is imperative to interpret this latter value with caution, as the blue component tends to attenuate, eventually vanishing in this direction. This PV map also reveals a faint, stationary emission at +18 km s⁻¹in the NE extremity, extending over 10 arcsec. The total emission of this slit spans 82 arcsec.

As illustrated in Figure 4, slits G2 and G3 intersect the complex transition region between the barrel-like structure and the outer “lobes”. Both PV diagrams display a pronounced expansion pattern that subsequently decreases in amplitude. In the case of G2, the strongest expansion is directed toward the southwest (SW) along the slit and reaches $\sim 60$ km s⁻¹. The blue component decreases to 22 km s⁻¹ before appearing to increase once again. The red component vanishes when the expansion reaches its minimum. Additionally, there is a stationary emission centred at $-30$ km s⁻¹. The total emission spans 42 arcsec. The PV map for G3 reveals the largest expansion of 69 km s⁻¹, with another of 27 km s⁻¹ towards the SW along the slit. This slit also intersects a portion of the SW lobe, displaying a faint emission at the systemic velocity. Marginal emission is slightly observed towards the northeast (NE), but it is not possible to ascertain its expansion.

4. Discussion

4.1. The Morpho-Kinematic Model of NGC 2371

In order to further investigate the morphogenesis of NGC 2371, we produced a model using the 3D modelling tool for astrophysics, ShapeX Ver. 3.5.1 [27]. This software allows the user to define different structures with velocity patterns to produce synthetic images and PV diagrams that are then compared with those obtained from observations. A good model is evaluated by its ability to simultaneously reproduce both the observed PV diagrams and the nebular image. Our proposal is to obtain the fundamental structures, not in detail, given that in some cases there are insufficient elements to reach the microstructures. Guided by the velocity structure suggested by Gómez-González et al. [20] (see Figure 11 in that paper), we started by adopting axisymmetrical structures for NGC 2371.

To model the nebula, we use the ShapeX software in an interactive and iterative manner. Our procedure is described in detail in Appendix B. When applied to our data, we model the [O iii]$\lambda 5007$ Å emission in both the direct image and the PV maps, adopting an elliptical shell oriented SE–NW in the plane of the sky. Its major axis is inclined, with the NW side redshifted and the SE side blueshifted. We will refer to this ellipsoidal shell as the “barrel”, and we treat its polar regions as separate structures, hereafter referred to as the “NW cap” and the “SE cap”. In addition to these structures, we propose a system of bipolar lobes with its major axis slightly misaligned with the corresponding major axis of the barrel. Unlike most PNe with bipolar morphology, the ends of the lobes need to be spatially extended to match what we observe in the direct image, resulting in a wide waist. After the completion of the [O iii] emission model, we proceed to propose a structure for the internal knots emitting in [N ii]$\lambda 6583$ Å.

These low-excitation [N ii]-emitting knots were modelled as small, nearly aligned cylinders, corresponding to collimated ejections with slightly different distances with respect to the geometrical centre. The eastern knot was placed at $\mathrm{PA}=+70$°, with an inclination of $i=23$°, and a deprojected velocity of $V_{\mathrm{dep}}^{\mathrm{E}}=75$ km s⁻¹. In contrast, the brightest knot, located to the west, was modelled with $\mathrm{PA}=-120$°, an inclination of $i=-10$°, and a deprojected velocity of $V_{\mathrm{dep}}^{\mathrm{W}}=20$ km s⁻¹. The nature of these structures will be discussed in the following section.

Table 2 summarises the morpho-kinematic characteristics of all the proposed structures, including their spatial and kinematic variables after derojection. For a detailed explanation of how these parameters were derived, we refer the reader to Appendix B. In this table, the lobes correspond to the full bipolar structure. $R_{\mathrm{p}}$ and $R_{\mathrm{eq}}$ denote the semi-major and semi-minor axes, respectively, while $V_{\mathrm{p}}$ and $V_{\mathrm{eq}}$ are the polar and equatorial expansion velocities, relative to the systemic velocity. For the “barrel”, $V_{\mathrm{p}}$ is not applicable due to its truncated major axis, whereas $V_{\mathrm{eq}}$ does not apply to the caps and knots. $V/R$ represents the proportionality constant of the homologous velocity law adopted for each structure, while ${\tau }_{\mathrm{k}}$ corresponds to the kinematical age taken from variables measured on the model, using the formalism from [28] as well as the value of distance obtained from the Gaia data mentioned before [18].

4.2. The Brilliant Knots in [N ii]

The most brilliant knots in [N ii] are not fully aligned themselves. One possible interpretation of this configuration is that the two knots are uncorrelated, as they are not fully aligned and exhibit different velocities. Gómez-Gonález et al. [20] argue that the shape of the [N ii]-emitting clumps is more consistent with that of photo-evaporative flows instead of jet-like features. However, the high velocity of (at least) the Easter clump challenges this idea.

An alternative interpretation could be that these knots represent a pair of bipolar, high-velocity jets that were originally aligned (their misalignment on the sky is around 10°). As a hypothesis, we propose that this emission corresponds to a pair of jets, originally ejected in opposite directions at the same velocity. Given the intense emission in [N ii] from these knots, we assume that they result from shocks caused by the bipolar jet interacting with the inner shell. Consequently, their apparent distance from the central star to their current position corresponds to the projected distance between the star and the inner shell, which is determined from the images. We can then deproject this element to obtain the inclination angle of each knot with respect to the line of sight, allowing us to deproject each velocity vector. With this information, we estimated a 3D angular misalignment of both vectors of $\simeq 16$°. The results of using this hypothesis are presented in Table 2.

The East jet has a deprojected velocity of 75 km s⁻¹, while the West jet, which is the brightest, is also much slower ($\simeq 20$ km s⁻¹). These features may result from the interaction of a high-velocity jet with circumstellar matter. Given the electron densities of the West jet ($n_{\mathrm{e}}=1640$ cm⁻³) and its surroundings ($n_{\mathrm{e}}\simeq 650$ cm⁻³), obtained by Gómez-González et al. [20], we can reasonably estimate, using basic momentum conservation equations (see Appendix C), that a mass element could decelerate from 75 km s⁻¹ to 20 km s⁻¹. Thus, the slower velocity of 20 km s⁻¹ can be explained by a deviation from the original direction, which is consistent with the observed loss of alignment. We recognise that this is a rather simple explanation. However, the small angular separation of $\simeq 16$°, its strong intensity in [N ii], and the consistency of the density estimates justify a more detailed modelling of this feature, which lies beyond the scope of the present work. Our working hypothesis is that the jets were originally ejected along a symmetry axis passing through the central star, with an initial velocity of at least 75 km s⁻¹. For this reason, we do not regard the kinematical age of 7900 yr for the W jet as a robust value

The resulting ShapeX model, incorporating all proposed structures, is presented in Figure 6. It provides a satisfactory reproduction of the PV diagrams. In Figure 7, the synthetic spectra derived from the model are displayed superimposed on observed PV diagrams selected from Figure 3 and Figure 5, illustrating the consistency between the observations and the model.

4.3. The Origin of NGC 2371

The morpho-kinematic analysis presented in this work provides important clues to the origin and evolutionary history of NGC 2371. The different structural components show a wide range of kinematical ages: the inner barrel (${\tau }_k\simeq 2500$ yr), the caps ($\simeq 1600$–2600 yr), and the bipolar lobes ($\simeq 4400$ yr). This spread suggests that the nebula did not form in a single, instantaneous ejection but rather through multiple episodes of mass loss and shaping.

The bipolar lobes are very likely the oldest structure in NGC 2371, which might have formed as a result of the main ejection during the transition of the progenitor star from the AGB to the post-AGB phase. Interestingly, the morphology of the lobes suggests that this structure was decelerated during the outflow, possibly through interaction with the surrounding medium.

Within the estimated uncertainties, the inner barrel and its associated caps seem to have formed as a result of another mass-ejection event; notably, such a secondary ejection might be linked to the same process that produced the current H-deficient [WR] characteristics of the CS of NGC 2371. It is currently accepted that [WR] stars originate from the born-again scenario [29,30], in which the CS of a PN undergoes a late thermal pulse after leaving the thermally-pulsating AGB phase; during this episode, the surface H is converted into He and C and violently expelled into the old PN, naturally giving rise to a double-shell PN morphology [31]. Finally, the star evolves into a second post-AGB phase, developing once again a fast stellar wind and a high-ionising photon flux.

The caps projected on the tips of the barrel may trace subsequent, somewhat more collimated outflows that have influenced the outermost regions of the nebula. For the SE cap, no clear morphological or kinematical differences are found with respect to the barrel. The only obvious distinction is the lower surface brightness of both caps relative to the barrel. In contrast, the NW cap not only expands more rapidly but also appears slightly detached from the barrel. This asymmetry is consistent with episodic and not strictly symmetric bipolar ejections; however, local conditions in the circumstellar environment, such as density gradients or anisotropies in the surrounding medium, could likewise have affected the NW cap differently. These properties seem to be consistent with what is observed in other born-again PNe, where the inner shell is not completely bipolar but exhibits a collection of clumps and filaments inside the old PN [32].

Nevertheless, additional effects produced by a binary companion are also evident. The pronounced bipolar morphology of the lobes suggests that shaping by a companion played an important role prior to the formation of the barrel-like structure. Moreover, the slight misalignment between the lobes and the barrel likely reflects precession – possibly driven by binary interactions – or a change in the dominant mass-loss axis over time. For example, the bright [N ii] knots provide evidence of even later ejection events, or of jets interacting with denser circumstellar material. Their misalignment and asymmetric velocities strongly suggest episodic jet activity, possibly associated with intermittent accretion processes in a binary system. This interpretation is consistent with hydrodynamical simulations of binary-driven outflows, where collimated jets interact with pre-existing shells and produce localised distortions in morphology and kinematics e.g. [33]. An additional point to consider is that the derived kinematical ages of these knots may not directly reflect their true time of origin, but rather the significant alterations in trajectory and velocity that they have experienced. Such deviations could represent indirect evidence of their interaction with an inhomogeneous circumstellar medium or, alternatively, the imprint of variable jet-launching conditions linked to binarity.

Although there is currently no direct evidence confirming that the CS of NGC 2371 is a binary, the observational evidence presented here strongly supports this idea. Asteroseismological surveys using TESS and Kepler have demonstrated their potential to reveal pulsations and rotational modulations in PG 1159-type stars, with several objects already showing clear variability signatures e.g. [34,35]. These techniques could be applied to the nucleus of NGC 2371 to search for evidence of a companion. Moreover, binary CSs are increasingly recognised as key drivers in the formation and evolution of complex PNe e.g. [36], with confirmed examples illustrating the diverse shaping mechanisms related to binarity. Therefore, we propose that future time-domain photometry or time-resolved spectroscopy be pursued to test whether the CS of NGC 2371 shares a binary-driven evolutionary history.

Another important aspect to highlight is that the final model of NGC 2371 required the large-scale structures (the bipolar lobes and the barrel) to display surface irregularities in the form of small ridges and valleys. Interestingly, their orientations do not seem to be directly related to the flow of material expelled by the central star. This might indicate the action of additional energetic processes that locally modify both the morphology and the kinematics. Such processes could involve hydrodynamical instabilities, interaction with pre-existing circumstellar material, or even collimated outflows driven by binary interactions. The very shape of the bipolar lobes of NGC 2371 is unusual: while morphologically similar cases exist e.g. NGC 6369, [37], the differences suggest that more than one shaping mechanism may be at play. These comparisons point towards the possibility that the evolution of NGC 2371 involved additional processes beyond those accounted for in our current models, thus opening the way for future observational and theoretical studies to test these scenarios.

5. Conclusions

Our morpho-kinematic analysis of NGC 2371 shows that its present-day morphology is the outcome of multiple episodic mass-loss events rather than a single outflow. The combination of a barrel-shaped shell, extended bipolar lobes, younger polar caps, and misaligned [N ii] knots points to a complex evolutionary history. The derived kinematical ages, spanning from 1600 to 4400 years, indicate successive episodes of mass loss with different geometries and timescales.

The peculiar morphology of NGC 2371 – its nearly perpendicular bipolar lobes, bipolarity without a pronounced waist, and surface distortions that do not align with the central outflow – cannot be explained solely by standard axisymmetric wind interactions. Instead, our results indicate that a combination of factors played a decisive role in shaping this nebula: a late thermal pulse that produced the H-deficient [WR] CS, binary-driven interactions, and episodic jet activity. NGC 2371 therefore stands out as an exceptional and highly unusual PN, pointing to physical mechanisms that remain poorly understood or only marginally explored in current models of PN formation and evolution.