Follow-Up of Extended Shells around B[e] Stars

B[e] stars are massive B type emission line stars in different evolutionary stages ranging from pre-main sequence to post-main sequence. Due to their mass loss and ejection events these objects deposit huge amounts of mass and energy into their environment and enrich it with chemically processed material, contributing significantly to the chemical and dynamical evolution of their host galaxies. However, the large-scale environments of these enigmatic objects have not attracted much attention. The first and so far only catalog reporting the detection of extended shells around a sample of B[e] stars was an Ha imaging survey carried out in the year 2001, and was limited to bright targets in the northern hemisphere. We have recently started a follow-up of those targets to detect possible evolution of their nebulae in the plane of the sky over a baseline of two decades. Furthermore, we extend our survey to southern targets and fainter northern ones to complement and complete our knowledge on large-scale ejecta surrounding B[e] stars. Besides imaging in Ha and selected nebular lines, we utilize long-slit and 3D spectral observations across the nebulae to derive their physical properties. We discovered pronounced nebula structures around 15 more objects, resulting in a total of 27 B[e] stars with a large-scale nebula. Here we present our (preliminary) results for three selected objects: the two massive supergiants MWC137 and MWC 314, and the unclassified B[e] star MWC 819.


Introduction
B[e] stars are massive B-type emission-line stars. Their spectra exhibit intense Balmer emission, numerous forbidden lines, as well as permitted lines of low ionized metals. These stars can be found in different evolutionary stages, from pre-to post-main sequence and spreading over a large mass range (from sub-solar up to 70 M for the currently known members 1 ). In addition, their spectral energy distributions display a strong infrared excess due to the hot, 500 to 1000 K, circumstellar dust (e.g., [1]). This dust is suggested to be located in a ring or disk-like structure based on polarimetric observations. The confinement of the dust into a disk has been confirmed by interferometric observations of the closest and brightest objects [2][3][4]. Furthermore, detailed analysis of interferometric and high-resolution spectroscopic observations revealed that the stars are surrounded by multiple rings of gas that revolve the central object on (quasi-)Keplerian orbits (see, e.g., [3][4][5][6][7] ). The location of the dust can (but must not) coincide with some individual gas rings. The distribution of rings around each object is found to be unique and independent on whether the star is presumably single or part of a binary system [6]. The composition and dynamics of these small-scale discs (several AU), in particular around the B[e] supergiants, have been studied intensively to unveil their possible formation mechanisms for which critical rotation, binary interaction, as well as slow-wind solutions and pulsation triggered ejections are among the suggested scenarios (for details see reviews [8,9]).
In addition to the small-scale dusty discs, B[e] stars are surrounded by large-scale (up to several pc) ionized gas, which has not been studied much. A pioneering work [10] presented an Hα imaging survey for 25 northern B[e] stars, and report the detection of extended shells around 50% of the targets with various shapes: spherical and ring nebulae, spiral-arm like features, bipolar and unipolar-lobe structures. This diversity in shapes indicates that the formation mechanism is most likely not a unique process. Furthermore, the physical mechanism creating such mass ejection is not well known. Therefore, to have a statistically meaningful sample of B[e] targets with extended nebulae, and to investigate details of their mass-loss history we started an observational campaign with two goals. Firstly, to perform a follow-up of the targets with already detected nebulae. Secondly, to enlarge the sample of B[e] stars with confirmed nebulae by adding fainter, yet unexplored objects located in the north and by extending the survey to the southern sky. In this proceedings we present results for selected objects of our project.

Observations and Data Reduction
Our follow-up Hα imaging survey is mostly carried out with the Nordic Optical Telescope 2 (NOT) at La Palma and with the Danish Telescope in La Silla, for the northern and southern targets, respectively. Average exposure times are 15-30 min. For deeper images of selected targets GMOS [11,12] at Gemini South was used while for very large nebulae we utilized the Isaac Newton Telescope (INT). Complementing long-slit spectroscopic observations have been secured already for a small sample of objects. These observations extend over the wavelength range from 6350 to 6850 Å with a spectral resolution of R ∼ 10 000 and have been carried out at the NOT as well. This wavelength range covers nebular emission lines-Hα 6563 Å, [NII] λλ6548,6583, [SII] λλ 6716,6731. In addition, we have exploited the Scanning Fabry-Perot Interferometer (FPI) mounted at the Russian 6m Big Telescope Alt-azimuth (BTA). More observational details of our survey will be presented in a forthcoming paper (Liimets et al. in prep). In this proceedings we will restrict to details about observations and data reduction for the targets under discussion. These are the two supergiant stars MWC 137, and MWC 314, and the object MWC 819, which currently has an unknown evolutionary state.
The Hα images of the B[e] supergiant MWC 137 were obtained with the NOT and the Danish Telescope in 2016 and 2019, respectively. For the spectral observations the long-slit possiblities at the NOT as well as the scanning FPI mode at the BTA were exploited. Further observational details of the data used to investigate MWC 137 can be found from [13,14].
The imaging observations of the large-scale bipolar nebula of MWC 314 have been performed with the Wide Field Camera at the INT, which provides a Field of View (FOV) 33 × 33 and a pixel scale 0 .33 pix −1 . On 2019 October 3 five dithered frames with an exposure time of 10 min each were acquired totalling 50 min. The narrowband Hα filter No. 196 was used, which includes also the [NII] doublet. Standard routines in IRAF 3 were used for the data reduction. Spectral observations of MWC 314 were performed at the BTA with the multi-mode focal reducer SCORPIO-2 [15], which provides FOV 6 × 6 with a 0 .4 pixel scale in the 2 × 2 binning mode. To cover the entire nebula, it had to be split into four fields. The 3D spectroscopic cubes obtained in the scanning FPI mode have a resolution of R∼16000 in the narrow spectral range containing the Hα emission line. The observed four fields were merged into a mosaic cube with an angular size 18 × 12 . Details of the FPI data reduction can be found in [16].
A new deep Hα image of MWC 819 was acquired on 2020 November 15 with GMOS attached to Gemini South (program ID GS-2020B-Q-210). The FOV was 5 .5 × 5 .5 with a pixel scale 0 .16 pix −1 in the utilized 2x2 binning mode. The observations were taken with filter G0336, which covers, besides Hα, also the nebular lines from the [NII] λλ6548,6583 doublet. To remove the gaps between the CCDs the total exposure time of 10 min was divided into 5 subframes of 2 min each. For the data reduction the official Gemini Observatory software DRAGONS [17] was used.
Finally, we have at our hands the raw Hα images from 2001 observed with the Mt. Palomar 60 inch telescope. Details about the observations are given in [10]. We have reduced the images using standard IRAF routines.

Results and Discussion
Our Hα imaging survey has been completed. To the already published Hα survey of 25 B[e] stars [10], we added 32 objects, 9 in the north and 23 in the southern sky. Similarly to the first survey, we report detecting large-scale nebular features around 50% of the observed stars, that is 4 in the north and 11 in the south. The full description of our survey will be published separately (Liimets et [13,14,18]. In the following, we summarize our work on MWC 137 and present preliminary results for two additional objects: MWC 314 and MWC 819.

MWC 137
The object MWC 137 has been classified as a supergiant star that has just left the main sequence. Its mass and age have been estimated to be 37 +9 −5 M and 4.7 ± 0.8 Myr, respectively [14]. Our detailed investigations of the various data sets presented in [13,14] have led to the following picture of the environs of the B[e] supergiant MWC 137.
The central star is surrounded by hot CO gas forming a subarcsecond size disc displaying Keplerian rotation. Further out from the center a prominent arcminute scale ionized nebula is visible ( Figure 1). In the Figure, the emission from the Hα+[NII] doublet is depicted in blue, while with red we present the warm and cool dust observed with the Spitzer Space Telescope in the 3.6 µm band. The dust visible in the infrared has a slightly larger size than the optical nebula, appart from the north-west (NW) direction, where the warm and cool dust seems to be absent (see Figure 1). The ionized gas has an elliptical shape with varying intensity areas, presenting several brigther filaments. We compared our newly obtained Hα image from 2019 with the image from 2001 [10]. Our baseline of 18.1 years did not reveal any morphological changes. Neither could we yet detect any expansion in the plane of the sky.
We further investigated the ionzied gas in the optical range using our various spectral observations. In Figure 2 left we present our radial velocity (RV) measurements of the combined [SII] λλ 6716,6731 doublet obtained with the scanning FPI which is providing spectral coverage over the whole FOV. The RV pattern is complex, which is also supported with our long-slit spectra obtained with the NOT+ALFOSC 4 at various position angles across the nebula (see Figure 4 in [13] and Figure 5 in [14]). However, the northern part is mostly blue shifted and the southern part red shifted. In addition, from our long-slit spectra we could estimate electron densities utilizing the same [SII] doublet and assuming a nebular temperature of 10,000 K. Similarly to the RV, the electron densities show a varying pattern across the nebula, from 0 to 800 cm −3 (Figure 2 right). The higher densities are at the regions of intense emission and at the central area.  Moving to even larger distances from the center of MWC 137 we find that the optical nebula is surrounded by a vast amount of cold molecular gas (see Figure 12 in [13]) detected with the Atacama Pathfinder EXperiment (APEX). The various molecular clouds resolved in the radio regime have masses from 15 to 245 M and seem to embrace the ionized nebula in the west, south, and east, clearly implying an interstellar origin rather than resulting from the evolution of the central star.

MWC 314
MWC 314 is one of the most luminous stars in the Galaxy [19]. The debate over its nature is ongoing-a Luminous Blue Variable (LBV, e.g., [20]) or a B[e] supergiant (e.g., [21]). It seems to be established though that the central object is a contact binary. But the masses of the system and of the individual stars are also controversial. Suggested are a total mass of the system of 66 ± 9 M [22] and of ∼20 M [23]. This discrepancy renders it difficult to assign the object a proper evolutionary state. MWC 314 is surrounded by a huge bipolar nebula extending 18' across ( [10], Figure 3), which considering the distance of 3 kpc [19] translates into a physical size larger than 13 pc.
We started the analysis with our Hα image from 2019, which we compared to the one from 2001 to find any morphological changes and/or expansion in the plane of the sky. Comparing the bipolar lobes in the east-west direction, no morphological changes were detected. For the expansion in the plane of the sky, we exploited the so-called magnification method (see e.g., [24,25]), which is able to detect expansion down to 0.1 pixels on the residual image that has been created from the images taken at two different dates. To properly apply the magnification method, further processing of the images was needed, which included pixel by pixel, seeing, and flux matching. Further details of this processing can be found from [14,25]. Unfortunately, our baseline of 18 years did not reveal any expansion in the plane of the sky. Next, from inspection of our deep image, a possible third lobe NW of the central object might be identified, which is barely seen on the image from 2001 [10]. To further investigate this lobe, we included it into our spectroscopic observations with the FPI. It is almost fully covered by the pointing P1 (see Figure 3). The other pointings were chosen such that they cover most of the north-east (NE) lobe (P2) and the lobe in the west (W, P3 and P4). The proper motion of the central star, according to Gaia DR2 [26], is −2.3 mas yr −1 and −4.8 mas yr −1 in Right Ascension and Declination, respectively. A blue arrow, indicating the direction of this movement, is drawn in Figure 3, which fits well with a plausible swept-back nature of the bipolar nebula similar to the swept-back nebula detected around the classical nova GK Persei 1901 (see Figure 2 in [27]). This could also mean that the western lobe, which has been previously considered in the literature as the plausible counterpart to the NE lobe, must not necessarily be connected to MWC 314. In fact, the NE lobe appears to be rather extended towards the west and south, but these features have no clear counterparts in NE direction.
To test whether the NW structure is the actual counterpart to the NE lobe, velocity measurements over the entire nebula are required. Unfortunately, the 3D spectroscopic observation (see Figure 4) had no signal in that area. However, as pointed out by Bruce Balick during the conference, our proposed new lobe in the NW seems to be replenished with gas, while the two previously considered lobes are not. This is, so far, the only argument speaking against the scenario of the new lobe to be related to the MWC 314. Without further observations, in particular a repetition of the velocity measurements, we cannot confirm yet the possible connection of the NW lobe with MWC 314. Thirdly, we analysed the velocity field of the lobes of MWC 314 which is presented in Figure 4. As was already mentioned, the signal-to-noise was too low to get any measurements at the position of the NW lobe. In addition, due to problems with subtraction of the bright airglow emission lines, the RV measurements suffer from artefacts, which appear in Figure 4 as concentric rings with artificially higher radial velocities. Further analysis using new data reduction algorithms is required to remove these spurious signals. Nevertheless, from the RV figure we can tentatively conclude that the northern parts are predominantly blue shifted while the southern parts are redder, broadly from −15 km s −1 to +40 km s −1 , respectively.

MWC 819
MWC 819 has long been listed as an unclassified B[e] star, but recently, [28] concluded from their analysis of new optical and infrared data that it might be a compact protoplanetary nebula. Consequently, the central object should currently be evolving into a white dwarf. On the left side panel of Figure 5 we present the image from 2001 [10] and on the right side panel we show our Hα image from 2020. The FOV of each image is 3 .6 × 3 .1 and MWC 819 is the bright source in the center of each image.
From the 2001 data MWC 819 was identified as having a unipolar nebula [10], visible as the faint roundish nebular feature in the NE. Our new and considerably deeper image resolves the previously detected features in more detail and reveals an additional possibly related nebular ejecta. In particular, the "unipolar lobe" (marked with a big black circle in our 2020 image) for which a zoom-in is presented with an orange square (FOV 38 × 36 ), resembles an individual nebula around a faint, possibly background object that was not detected in the 2001 image. To confirm this hypothesis, additional observations and analyses are needed. We would also like to draw the readers attention to the smaller nebular features in the south-west (SW) marked with two smaller black circles. Those features were not mentioned in [10]. However, after being clearly resolved in our deeper image, faint counterparts can be detected also in the 2001 image. Furthermore, with red we mark previously not detected features. One in the NE and the other in the NW direction. The latter is especially interesting, because due to its position perpendicular to the main nebular features in the NE to SW direction, it might be possible that a jet is emanating from MWC 819. In fact, jets have been discovered in a few other B[e] stars (e.g., MWC 137 [13,29], MWC 922 [30], Z CMa [31]).

Conclusions
We have performed a new Hα imaging survey of 32 B[e] stars. Of these, 9 targets are located in the northern hemisphere, and 23 in the southern one. Moreover, we re-observed the 25 B[e] stars, whose environment had been imaged about two decades ago [10]. Our study has two main goals: (i) to search for indications of large-scale ejecta around the yet unstudied objects, and (ii) to investigate the morphology of the nebulae and their expansion dynamics. Furthermore, all objects with detected circumstellar nebulae are observed spectroscopically in a systematic way to derive the physical parameters such as temperature and density distribution across the nebulae together with their ionization structure and chemical composition. The combined information will be used to address the issue of the yet unknown mass-loss history of these enigmatic objects. Moreover, we expect that the detailed information about the nebulae gained during our project will help to settle the yet unknown evolutionary state of many of the objects.
Our survey of the newly imaged objects has revealed that about 50% (15 objects out of 32) of the B[e] stars have detectable nebular features. This detection rate is similar to the previously reported one (12 out of 25 objects, [10]). The shapes of nebulae display a high diversity, and a full catalog with all the images is in preparation.
While for most objects the analysis is still ongoing and, in particular, a wealth of spectroscopic data still needs to be acquired, we have presented our (preliminary) results for three objects. Two of them (MWC 137 and MWC 314) had been previously imaged, but our current baseline of about 18 years provided neither evidence for changes in nebular morphology nor indications for nebular expansions. However, both objects have rather large distances (about 5 and 3 kpc, respectively), so that our current spatial resolution is still too coarse to expect detectable changes, given the moderate radial velocities measured in both nebulae. Our deep image of the presumably bipolar nebula around MWC 314 pointed to a third, faint lobe filled with emission, which would fit better to the direction of the proper motion of the star. Whether it can indeed be associated with MWC 314 requires better radial velocity measurements of the entire nebular structure than the one we could secure so far.
For the third object presented, MWC 819, we have also an old set of recorded data at hand. Our new, significantly deeper image suggests that one of the previously detected nebular features has been mistakenly identified as uni-polar lobe connected to MWC 819. It seems that this feature could be associated instead to a faint back-ground source, which was not visible on the old image. However, our new, high-quality image displays other, yet unreported signs for ejected matter, maybe in relation with a jet. Confirmation is still pending and requires high-quality spectra, which are currently under way, to determine the velocity structure of the detected features.

Institutional Review Board Statement: Not applicable
Informed Consent Statement: Not applicable Data Availability Statement: Most of the data presented in these proceedings can be acquired using the following public archives: NOT FITS Header Archive 5 for NOT data, Isaac Newton Group Archive 6 for INT data, Gemini Observatory Archive 7 for GMOS data, ESO Archive 8 for APEX data, and The General Observation Archive 9 for BTA data. Observations obtained with the Danish and Mt. Palomar 60 inch telescope can be obtained from TL on a reasonable request.
Newton Group of Telescopes in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias. We thank Bruce Balick for a valuable discussions on the nebula of MWC 314.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: We would like to emphasize that stars with the B[e] phenomenon are not predicted by any of the currently known stellar evolution theories. This might be related to the fact that these models are purely 1D and do not account for aspherical mass loss, which is typically seen in B[e] stars. Hence this limit of 70 M might not be real and more massive B[e] stars might exist, which just have not been discovered yet.