Next Article in Journal
Recent Progress in Modeling the Macro- and Micro-Physics of Radio Jet Feedback in Galaxy Clusters
Next Article in Special Issue
Dense Molecular Environments of B[e] Supergiants and Yellow Hypergiants
Previous Article in Journal
Investigating Gravitationally Lensed Quasars Observable by Nancy Grace Roman Space Telescope
Previous Article in Special Issue
Synthetic Light Curve Design for Pulsating Binary Stars to Compare the Efficiency in the Detection of Periodicities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New Insight into the FS CMa System MWC 645 from Near-Infrared and Optical Spectroscopy

by
Andrea Fabiana Torres
1,2,*,
María Laura Arias
1,2,
Michaela Kraus
3,
Lorena Verónica Mercanti
1,2 and
Tõnis Eenmäe
4
1
Departamento de Espectroscopía, Facultad de Ciencias Astronómicas y Geofísicas, Universidad Nacional de La Plata, Paseo del Bosque S/N, La Plata B1900FWA, Argentina
2
Instituto de Astrofísica de La Plata (CCT La Plata-CONICET, UNLP), Paseo del Bosque S/N, La Plata B1900FWA, Argentina
3
Astronomical Institute, Czech Academy of Sciences, Fričova 298, 251 65 Ondřejov, Czech Republic
4
Tartu Observatory, University of Tartu, Observatooriumi 1, 61602 Tõravere, Estonia
*
Author to whom correspondence should be addressed.
Galaxies 2023, 11(3), 72; https://doi.org/10.3390/galaxies11030072
Submission received: 24 April 2023 / Revised: 31 May 2023 / Accepted: 1 June 2023 / Published: 10 June 2023
(This article belongs to the Special Issue Theory and Observation of Active B-type Stars)

Abstract

:
The B[e] phenomenon is manifested by a heterogeneous group of stars surrounded by gaseous and dusty circumstellar envelopes with similar physical conditions. Among these stars, the FS CMa-type objects are suspected to be binary systems, which could be experiencing or have undergone a mass-transfer process that could explain the large amount of material surrounding them. We aim to contribute to the knowledge of a recently confirmed binary, MWC 645, which could be undergoing an active mass-transfer process. We present near-infrared and optical spectra, identify atomic and molecular spectral features, and derive different quantitative properties of line profiles. Based on publicly available photometric data, we search for periodicity in the light curve and model the spectral energy distribution. We have detected molecular bands of CO in absorption at 1.62 μ m and 2.3 μ m for the first time. We derive an upper limit for the effective temperature of the cool binary component. We found a correlation between the enhancement of the H α emission and the decrease in optical brightness that could be associated with mass-ejection events or an increase in mass loss. We outline the global properties of the envelope, possibly responsible for brightness variations due to a variable extinction, and briefly speculate on different possible scenarios.

1. Introduction

In their evolution, some B-type stars undergo phases that are still puzzling for astrophysicists even after years of study since they develop certain peculiarities that are not yet well understood. The B[e] phenomenon displayed by several B-type stars in their optical spectra is an example of them. Its manifestation can be seen through the presence of permitted and forbidden low-excitation emission lines of neutral and low ionization metals arising from circumstellar (CS) gas and large infrared excess due to CS dust [1]. The phenomenon is associated with stars with different initial masses, isolated or in binary systems, transiting different evolutionary stages, such as supergiants, compact planetary nebulae, Herbig Ae/Be stars, and symbiotic systems [2]. Despite the cited differences among the stars, the physical conditions of their CS gaseous and dusty envelopes are similar, which is a crucial factor when seeking to understand the development of the phenomenon. Furthermore, since the CS envelopes veil the photospheric features of the central objects, it is difficult to determine their spectral types and evolutionary states. Therefore, stars without a proper classification comprise the category named “Unclassified B[e] stars” (UnclB[e]).
Nearly a decade after Lamers’ classification, Miroshnichenko proposed a new group called FS CMa stars [3]. The observational defining criteria are (1) the presence of a hot star (between O9 and A2 spectral types) continuum with emission lines of H I, Fe II, O I, [Fe II], [O I], Ca II; (2) an infrared spectral energy distribution (SED) that shows a large excess with a maximum at 10–30 μ m and a strong decrement beyond these wavelengths; and (3) a star location outside a region of star formation. According to these properties, almost all the UnclB[e] objects in Lamers et al.’s publication are in this group [4], which has approximately seventy members between confirmed and candidate ones [5]. They are suspected to be binaries at a post-mass-exchange evolutionary phase, with a secondary component fainter and cooler than the primary or degenerate [3,6,7]. Since the predictions of mass-loss rates of single-stars theory [8,9] cannot explain the existence of a large amount of CS matter, a mass-transfer process in a binary system could be a likely explanation. However, only a few objects of this class have been confirmed as binary systems, probably due to the scarcity of available observational data to discover them and the difficulties in detecting signs of binarity due to the presence of CS matter, the intrinsic stellar variability, and the low brightness of most members of the FS CMa group [10]. Recently, Miroshnichenko et al. [11] published a review of FS CMa objects, where they reported fifteen stars as binaries and six as binary system candidates.
MWC 645 (= V2211 Cyg, α = 21:53:27.49, δ = +52:59:58.01; V = 13.0, H-K = 1.53, J-H = 1.267) was originally included in the supplement of the Catalogue of Mount Wilson about A and B stars with bright H I spectral lines [12]. The presence of strong double emission lines of Fe II and [Fe II] (with a radial velocity difference between the red and blue peaks of 150 km s 1 ) and triple-peaked profiles of the H γ and H δ transitions were reported by Swings and Allen [13]. They also remarked striking spectral similarities between MWC 645 and η Car. Also, permitted and forbidden transitions of low excitation and ionization potential belonging to Ti II, Cr II, [O I], and [N II] were observed [1]. Photometric variations were found by Gottlieb and Liller [14] with an amplitude of 0.3 mag and a possible period of 23.6 years. A deep spectroscopic study was done by Jaschek et al. [15] that revealed no stellar absorption features. These authors concluded that possibly MWC 645 is a late B-type object based on the absence of He II lines and the weakness of the He I lines at λ 6678 Å and λ 7065 Å possibly detected once, each one at different years. They did not find spectral transitions from C, Ne, and Mg atoms or ions, but they found lines of K I and Cu II (typically seen in stellar types later than F) and Zr II (usually seen in stars later than A0-type). They highlighted the extreme spectroscopic variability of MWC 645 over the years. Lamers et al. [2] included it in the UnclB[e] stars group. MWC 645 has IRAS flux ratios that locate it in the region occupied by OH/IR stars [16]. Zickgraf [17] detected a characteristic asymmetric profile for the emission metal lines, with a steep red flank and a blue wing. He reported the splitting in the central emission of [O I] and [Fe II] lines and a peculiar emission profile of the H α line showing a broad blue and a narrow red component with a full width at half maximum (FWHM) of 5.0 Å and 1.3 Å, respectively. He proposed a latitude-dependent wind model with a large optical depth dust disk at an intermediate inclination to explain the asymmetric line profiles and their splitting. Marston and McCollum [18] obtained H α narrow band imaging and found no visible extended emission associated with the star.
Recently, Nodyarov et al. [19,20] studied high-resolution optical spectra of MWC 645 taken in two different years. He found absorption lines of neutral metals, such as Li I, Ca I, Fe I, Ti I, V I, and Ni I, typically present in cool stellar spectra, with a different average radial velocity in each spectrum, that revealed the binary nature of the object. However, they did not find any absorption line typical of a B-type object in any of their spectra. They disentangled the contribution of each stellar component and estimated their surface temperatures and luminosities (T e f f = 18,000 ± 2000 K and 4250 ± 250 K, log (L/L ) = 4.0 ± 0.5 and 3.1 ± 0.3 for the hot and cool components, respectively). Low-resolution near-IR spectra displayed emission lines of the H I Paschen and Brackett series, as well as of Fe II, O I, N I, and He I. Photometric monitoring in the optical and near-IR regions showed quasi-cyclic variations of both short and long periods (months and ∼4 years, respectively). The authors conclude that the star can be classified as an FS CMa-type object, where its intermediate-mass components (7 M and 2.8 M ) undergo an ongoing mass-transfer process. According to the shape displayed by the spectral energy distribution with weak emission peaks at about 10 μ m and 18 μ m, they inferred the presence of silicates in an optically thin dusty shell. MWC 645 is one of the eight FS CMa objects in which absorption lines of neutral metals typical of late-type secondaries have been detected. To contribute to the study of this intriguing object, we decided to observe it in the near-IR to search for signatures of both stars, mainly of the cool component, that help to characterize it. In addition, the acquisition of new optical spectra and the public availability of data (spectroscopic and photometric) that could shed some light on this complex system motivated us to analyze them. The paper is organized as follows: We present the infrared and optical observations used in this work in Section 2. In Section 3 and Section 4, we analyze the data. In Section 5, we discuss the results. Finally, Section 6 contains the main conclusions.

2. Observations

2.1. Near-Infrared Spectra

Near-infrared spectra were taken using the Gemini Near-Infrared Spectrograph (GNIRS, [21]) attached to the 8 m telescope at GEMINI-North (Hawaii) under the programs GN-2017A-Q-62, GN-2018A-Q-406, and GN-2022B-Q-225. On 6 June 2017, we obtained K-band spectra in long-slit mode centered at 2.35 μ m. The instrumental configuration used was a 110.5 l/mm grating, a 0.3 arcsec slit, and the short camera (0.15 arcsec/pix). We also acquired spectra with the same configuration but in cross-dispersed mode centered at 2.19 μ m and 2.36 μ m on 24 and 30 July 2018, respectively. The effective spectral coverage by these set-ups was 0.90 μ m–2.27 μ m and 0.85 μ m–2.45 μ m, respectively, with gaps between the orders (the interval 1.36 μ m to 1.46 μ m is unusable due to saturated telluric lines). The resulting mean spectral resolving power of the spectra was R ∼ 5500. On 24 August 2022, L-band spectra were obtained with a different long-slit configuration: a 31.7 l/mm grating, a 0.1 arcsec slit, and the long camera (0.05 arcsec/pix), with two different central wavelengths (3.48 and 4.00 μ m). This configuration resulted in R ∼ 5100. The spectra were taken in two ABBA nodding sequences along the slit. To account for telluric absorption, a late-B- or an early-A-type star close to the target in both time and position was observed. Stars of these spectral types are featureless in the observed wavelength range, except for hydrogen absorption lines that can be successfully removed in the reduction process by fitting theoretical line profiles. Flats were also acquired. The data were reduced with the Image Reduction and Analysis Facility (IRAF)/Gemini tasks. The sky contribution was removed by subtracting the AB pairs. The spectra were flat-fielded and telluric corrected. The wavelength calibration was performed using the telluric lines. The data were normalized to unity.
A and B positions were added to increase the signal-to-noise ratio (S/N). The final S/N ratio varies for the different spectral ranges, as it is affected by the quality of the telluric correction. Some regions are very polluted with telluric lines and it was impossible to make a complete cancellation, thus some residuals remain. In addition, for some spectral regions heavily crowded by emission lines, it becomes difficult to make accurate S/N ratio estimates. Table 1 summarizes the mean values of the S/N ratio for all our GNIRS near-IR observations.

2.2. Complementary Data

Optical observations were carried out at Ondřejov Observatory, Czech Republic, using the Coudé spectrograph [22] attached to the Perek 2 m telescope. We obtained spectra with a resolving power of R ∼ 12,000 covering a spectral range from 6262 Å to 6735 Å on 12 and 13 September 2018. We also acquired a spectrum centered at 8600 Å. We used a grating of 830.77 l/mm with a SITe 2030 × 800 CCD and a slit width of 0.7 arcsec. Additional observations were done at Tartu Observatory, Estonia, with the 1.5 m Cassegrain reflector AZT-12 on 1 November 2021, using the long-slit spectrograph ASP-32 with a 600 l/mm. The wavelength coverage extended from 5450 Å to 7480 Å. Data were processed using standard IRAF tasks. Spectra were bias and flat-field corrected, wavelength calibrated, heliocentric velocity corrected, and flux normalized.
We also searched for available optical spectra in the BeSS database [23]. We downloaded sixteen spectra taken between 2019 and 2022, with a resolving power R ∼ 14,000/16,000 in the spectral range 6500–6600 Å. In addition, we collected a lower resolution spectrum (R ∼ 5000) acquired on 30 August 2019 that covers the range 6150–7000 Å. The spectra were corrected by heliocentric velocity and normalized to the continuum using the standard IRAF tasks. We chose the same sample of continuum points to normalize all spectra. The telluric correction has not been applied.
In addition, we extracted from the public database ASAS-SN (All-Sky Automated Survey for Supernovae; Shappee et al. [24], Kochanek et al. [25]) (https://www.astronomy.ohio-state.edu/asassn/ (accessed on 2 February 2023)) survey photometric data of this star obtained over eight years. The collection is composed of V-band magnitudes from 16 December 2014 to 29 November 2018 and data in the g-band from 12 April 2018 up to 17 January 2023. Furthermore, we collected ground- and spaced-based multicolor photometry from Vizier service from 0.3 μ m to 140 μ m.

3. Analysis of the IR Data

Figure 1 shows the near-IR spectrum of MWC 645 from 8400 to 13,600 Å. It displays numerous emission lines, particularly of the H I Paschen series. The strongest lines, except those of H I, correspond to O I, Fe II, and the Ca II triplet. Many transitions of N I can be identified in emission along this spectral range. Forbidden lines of [Fe II] and [S II] are also present. The moderate-resolution data reveal several absorption lines of Fe I that could be associated with the cool stellar companion.
We carefully searched for He I lines in our spectra. If the lines are present, they are incipient and hidden in the noise. The most intense transitions in the interval λ 8400–24,500 Å (cited in the NIST database [26]) correspond to λ 10,830 Å and λ 20,587 Å. Nodyarov et al. [20] reported the presence of the He I λ 10,830 Å transition in emission from a low-resolution spectrum (R ∼ 700). We identified a group of Fe II lines with the three emission peaks at the interval 10,826–10,862 μ m. However, the bluest feature of this group (see Figure 1) is broad, and thus, the He I λ 10,830 Å line could be blended with the Fe II lines. Our data have a higher resolution than the spectrum of Nodyarov et al., but not sufficient to separate the Fe II lines from the He I line. These authors also identified the He I λ 20,587 Å line. Unfortunately, it lies outside our spectral coverage.
The H-band spectrum of MWC 645 (upper panel of Figure 2) is dominated by the H I Brackett series and several permitted and forbidden lines of Fe II. In the K-band, the most intense feature is the Br γ line (see lower panel of Figure 2), which stands out among several emission lines corresponding to Fe II, [Fe II], and presumably [Ni II]. The Mg II doublet at λ λ 21,374 Å and 21,437 Å is also in emission. The Pfund series extends from 2.3 microns longward. In addition, absorption features of neutral metals characteristic of late-type stars, such as Ca I, Mg I, and Na I, are present. For the first time, we have detected the presence of CO band heads in absorption around 2.3 μ m and around 1.6 μ m, which are typical photospheric features of late-type luminous stars.
Figure 3 shows the first obtained L-band spectrum of MWC 645 in two different spectral regions. The first interval between 33,310 Å and 36,410 Å is relatively featureless (see left panel), except for the presence of the permitted emission line of Fe II λ 35,423 Å, which is clearly seen above the continuum level, and some H I lines of the Humphreys series in emission, where the strongest is the one corresponding to the 20-6 transition. Lower-order members of the Humphreys series can be seen in the second spectral interval (right panel) that ranges from 38,520 Å to 41,730 Å, where the strong emission of the Br α line can also be observed. We searched for absorption bands of the first-overtone of silicon monoxide (SiO) around 4 μ m but found none.
Figure 4 plots three H I lines: Pa β , Br γ , and Br α . Their profiles are single-peaked but asymmetric. We measured the total equivalent width (EW) of each line using the ‘e’ function in the IRAF splot routine. The total measured EWs are 184 Å, 15 Å, and 48 Å, respectively. The percentage uncertainty of the EW measurements is 5%. Unfortunately, the wavelength calibration of our near-IR spectra is not accurate enough to determine reliable radial velocity measurements.

Co Absorption Bands

Figure 5 (upper panel) displays the second-overtone band heads of 12 CO in absorption from the 2018 H-band spectrum. The lower panel compares the K-band spectra taken in 2017 (in red) and 2018 (in black), where the variation of the first-overtone band heads of 12 CO is clearly seen. The positions of the 13 CO band heads are also marked and clearly detected in the spectrum from 2017. The 2018 spectrum is too weak to see these faint features.
The strength of the CO absorption bands in the near-IR spectra of classical late-type stars depends on the stellar effective temperature, T e f f , and surface gravity, log g [27,28]. The CO absorption becomes deeper when the effective temperature decreases and the luminosity increases. Thus, hot star spectra display no trace of CO features (T e f f 5800 K–6000 K, Ali et al. [29]), and dwarf stars present weaker CO absorption bands than supergiants. To characterize the cool companion of MWC 645, responsible for the CO absorption features, and estimate its fundamental parameters, we used the IRTF (NASA Infrared Telescope Facility) Spectral Library [30,31], which collects stellar spectra observed with the spectrograph SpeX at a resolving power of R ∼ 2000 and an S/N ratio of about 100 at λ < 4 μ m. We looked for late-type stars with spectral types between F and M and luminosity classes between I and V to compare their spectra with our spectrum from 2017 in the wavelength range from 2.26 to 2.44 μ m. Figure 6 shows this comparison. The MWC 645 spectrum (solid black line) was degraded to the resolution of the template spectrum (dashed red line), which corresponds to a G0 Ib-II star (HD 185018). The intensity of the first 12 CO band head of MWC 645 coincides reasonably well with that of the early G-type star; however, the rest of the band heads are less intense. The blue edge of the CO(2-0) band head might present an incipient emission. The absorption of the first 13 CO band head is more intense than that displayed by the library star. According to Wallace and Hinkle [32], the 13 CO isotope is prominent in the supergiants and giants but is not apparent in the dwarfs, although its strength also depends on the initial rotation velocity of the star and the mixing processes that can cause a surface enrichment in 13 C. Otherwise, the absorption lines of neutral metals are less intense than those in the template spectrum, indicating an earlier spectral type (F8-F9 subtypes). Furthermore, the lack of SiO band heads at 4 μ m, often observed in K0-type stars and later, also points towards an earlier type [33].
Winge et al. [34] presented a spectroscopic library of late spectral-type stellar templates in the K-band at a resolving power of R ∼ 5900. The authors plotted the equivalent width (EW) of the first CO overtone as a function of T e f f (see their Figure 3) for a stellar sample with T e f f in the range 3200–5200 K and different luminosity classes. They measured the EW from the blue edge of the (2-0) band head to the blue edge of the (3-1) band head, more precisely in the window 2.293–2.322 μ m. We measured the EW of the CO(2,0) band head from the spectrum of 2017 and obtained 2.55 ± 0.5 Å. A visual extrapolation of the relation seen in the figure between EW and T e f f in the hottest edge of the plot gives an estimation of T e f f around 5200 ± 100 K. From the spectrum obtained in 2018, we measured an EW of the CO(2,0) band head equal to 1.22 ± 0.1 and estimated a T e f f 5300 ± 100 K.

4. Analysis of the Optical Data

4.1. Photometric Light Curve

Figure 7 shows the light curve of MWC 645 taken from ASAS-SN. We applied the relationship derived by Nodyarov et al. [20] to convert the g-band magnitudes into V-band magnitudes. The optical photometry in the V-band acquired by the authors mentioned above is also included. The dates of the spectroscopic observations presented in this work are marked in the light curve as vertical lines.
The brightness fluctuations of the star up to July 2022 have been reported by Nodyarov et al. [20], who suggested that a new minimum in the light curve might take place in the second half of 2022. As can be seen in the plot, the star continued fading up to the end of October approximately, reaching a minimum of ∼0.1 mag brighter than the minimum that occurred in August/September 2018. Then, it began to strengthen in brightness again.
Nodyarov et al. [20] searched for periodicity, excluding visual magnitudes greater than 13.2 mag from their analysis. They derived 69, 145, and 295 days. They attributed the quasi-cyclic photometric variations to variable CS extinction. We applied the Lomb–Scargle method using the IRSA (https://irsa.ipac.caltech.edu/irsaviewer/timeseries (accessed on 30 November 2022)) time series tool to the V-band light curve shown in Figure 7. We discarded the magnitudes with errors greater than 0.03. The scan of periodic signals with values below ten days gave strong peaks at one day and harmonics of the sidereal day due to the observing cycle.
The periodogram for periods greater than one day is shown in Figure 8. The six peaks at or above a confidence level of 20 in the power spectrum correspond to periods of approximately 65, 112, 162, 298, 461, and 709 days. We dismissed the last period since the time coverage of the observations is not enough for its precise determination. The phase diagram for each period shows a large scatter of the magnitude points and a small amplitude in their modulation (∼0.2–0.3 mag). We should note that the 65- and 298-day periods were also found by Nodyarov et al. [20]. As our data spread over a more extended baseline than the one used by the authors mentioned above, this might be a possible explanation for the differences in the other identified periods.

4.2. Spectroscopic Data

The peculiar profile of the H α line of MWC 645, composed of a broad blue-shifted peak and a narrow red-shifted one, can be seen in Figure 9, where all the spectra are normalized to the continuum level. This plot shows not only the profile changes over different years (spectra from the same year are displayed in the same color) but also daily. We note that as the BeSS spectra are not corrected by telluric lines, some of the profiles present one or two absorption features superimposed on the blue-shifted emission peak corresponding to water vapor lines of the Earth’s atmosphere, which affect the shape of the profile. A variation in the emission strength of both peaks is seen. Using the ‘e’ task in the IRAF splot routine, we measured the intensity of the blue (V) and red (R) emission peaks and the total equivalent width of the H α profile, except for the lowest resolution spectrum. Table 2 presents these values and the calculated V/R ratios. We can see that the V/R ratio presents changes over four years, even doubling its value. We note that the ratio of V/R ∼ 0.3 corresponds to observations close in time (except for one). For this subset, the changes in EW might be mainly due to the continuum-level variations.
The average radial velocity of the blue and red emission components derived from the Ondřejov spectra are −225 ± 5 km s 1 and −31 ± 3 km s 1 , respectively, which are in agreement with the values reported by Zickgraf [17] of −218 km s 1 and −30 km s 1 and Nodyarov et al. [20] of −252 ± 9 km s 1 and −30 ± 2 km s 1 , respectively. Fitting a Gaussian profile to the narrow red component of the H α line, we obtained an average FWHM of 90 ± 1 km s 1 . To fit the broad emission component, we built a profile with a red wing symmetrical to the observed blue one and obtained an average FWHM of 318 ± 4 km s 1 .
Even though the BeSS material is not accurate enough to measure radial velocities, we have estimated them from the different spectra for both emission peaks fitting the components with Gaussian profiles. The average value is −229 km s 1 and −26 km s 1 for the blue and red emission peaks, respectively. In Figure 9, a variation in the central wavelength of the H α red emission peak (which is not distorted by telluric lines) can be observed from the different spectra; however, as the wavelength calibration of the BeSS spectra is not well suited for radial velocity determination, we cannot confirm if this change is real. The average FWHMs of the blue and red components are 256 ± 4 km s 1 and 80 ± 2 km s 1 , respectively. The H α broad component line profile from the BeSS spectra has a smaller average FWHM than the Ondřejov spectra. In the latter, the broad component presents a blue wing with a gentler slope that is also outlined for the red wing, giving a greater width.
Apart from the H α line, our medium-resolution spectroscopic observations and the 2019 low-resolution BeSS spectrum also show the lines of [O I] λ λ 6300, 6364 Å. These lines appear single-peaked, although asymmetric, as was previously mentioned by Zickgraf [17], which might suggest that the lines are composed of two blended components (see Figure 10). Several permitted Fe II lines and the forbidden lines of [N II] λ 6583 Å and [S II] λ λ 6716, 6731 Å are apparent. We calculated the average heliocentric radial velocity of the emission lines of the Ondřejov spectra and the standard error of the mean, obtaining −43 ± 2 km s 1 . Jaschek et al. [15] derived −76 ± 5 km s 1 , and Nodyarov et al. [20] found −61 ± 4.3 km s 1 , which indicates a variation in radial velocity. The He I λ 6678 Å transition is absent, as in the spectra studied by the last-mentioned authors and Zickgraf [17].

4.3. Global Properties of the Circumstellar Material

Figure 11 shows the spectral energy distribution (SED) of MWC 645, built from the photometry publicly available from the ultraviolet to the far-IR (0.3 μ m–140 μ m). The low-resolution spectrum taken in 2021 in the spectral region of the H α line is also included.
To derive the global physical properties of the CS material, we considered a simple model presented by Marchiano et al. [35] and Arias et al. [36]. The numerical code allows for obtaining the SED assembled by different envelope components. The model assumes the presence of a spherical envelope composed of gas close to the star (≤ 5 R * ) and (or) dust further away from it (≥ 100 R * ) [37,38]. The emergent flux is computed from the central star and the envelope (considering that the latter can be reduced to an equivalent shell), applying a plane-parallel solution for the transfer equation. The optical depth τ λ G and the source function characterize the gaseous shell, which can be described by adding as free parameters the electron temperature T G and the effective radius R G . The dusty region is treated using an analogous scheme with similar parameters describing the shell: an optical depth τ λ D , a temperature T D , and an effective radius R D . The model allows several dust shell components to be added. The interstellar extinction is also included by an optical depth τ λ I S M . The absorption A( λ ) is related to each optical depth through the expression τ = 0.4 ln ( 10 ) A( λ ). Using the law given by Cardelli et al. [39], it can be written as A( λ ) = [R V a(1/ λ ) + b ( 1 / λ ) ] E(B-V), where R V is the total to selective extinction and E(B-V) is the color excess. We took R V I S M = 3.1 for the interstellar dust and tried different values of R V D greater than 3.1 for the CS dust shell components. The temperature of the dust grains depends on the stellar radiation and the distance from the star center [40]. That is, T D (r) = T e f f W(r) [ 1 / ( 4 + p ) ] , where W(r) is the geometrical dilution factor. The parameter p depends on the nature of the dust, but it is usually on the order of one. Furthermore, as the equilibrium temperature should be lower than the dust condensation temperature (typically around 1500 K) to allow the formation of grains, it constrains the distance where condensation can occur, e.g., for a T e f f of 18,000 K, the condensation distance is about 249 R * .
The code gives the observed flux normalized to that at a reference wavelength λ r e f ; we chose λ r e f = 0.55 μ m. We assumed that the flux of the central object results from the contribution of both stellar fluxes, for which we considered the Kurucz [41] atmosphere models. Regarding what we know about the stars, we selected models between 17,000 and 20,000 K for the hot binary component and between 4500 and 6000 K for the cool one and explored the flux contribution of the hot and cool components to the total flux in the range of 70%–90% and 30%–10%, respectively (similar percentages were suggested by Nodyarov et al. [20]). Figure 11 displays our best fit of the theoretical SED (solid blue line) to the observed data. The best-fitting model was computed by taking the photospheric fluxes from a star with T e f f = 18,000 K, log g = 4.0, and R * = 3.73 R and a cool star with T e f f = 5000 K and log g = 1.5. The contribution of each stellar flux to the total flux is 80% and 20% for the hot and cold components, respectively. The resulting envelope has one gaseous shell at R G = 1.15 R with T G = 16780 K and τ V = 0.1. The dusty region comprises three different shells with the following parameters: R D 1 = 348 R * , T D 1 = 1310 K, R D 2 = 3750 R * , T D 2 = 507 K, R D 3 = 0.02 pc, T D 3 = 98 K. The computed total CS visual absorption is A V D = 0.097 mag. We obtained a color excess due to the interstellar medium of E(B-V) I S M = 0.98 ± 0.02 mag, which results in a total visual absorption A V = 3.13 ± 0.11 mag. This value agrees with the one derived by Nodyarov et al. [20]. A disagreement between the theoretical and observed SEDs in the region of log λ = 1.0–1.5 is observed. We should note that the employed code can model the thermal emission of the dust, but it cannot address the computation of silicate bands. Thus, the presence of silicate particles might be responsible for the observed difference between the SEDs at 10 μ m and 18 μ m. In fact, Nodyarov et al. [20] have already reported weak emission bumps at these wavelengths.
Figure 11. Spectral energy distribution of MWC 645. The open triangles represent the observed photometric data: optical bands (yellow), 2MASS (green) [42], WISE (violet) [43], MSX (light blue) [44], IRAS (red) [45], and AKARI (black) [46]. The error bars of the photometric data are included (in most photometric bands, they fall inside the symbols). The low-resolution spectrum acquired in 2021 over the H α region is also displayed (with a dashed black line). The solid red line shows the SED modeled considering the contribution of the photospheric fluxes from both stars, the thermal emission from a gaseous shell close to the system and the effect of the interstellar medium extinction. The solid blue line shows our best-fitting theoretical SED, obtained by adding to the SED plotted in red the contribution of three dusty shells surrounding the stellar system. The flux is normalized to that at λ r e f = 0.55 μ m and displayed on a logarithmic scale. The wavelengths are in microns.
Figure 11. Spectral energy distribution of MWC 645. The open triangles represent the observed photometric data: optical bands (yellow), 2MASS (green) [42], WISE (violet) [43], MSX (light blue) [44], IRAS (red) [45], and AKARI (black) [46]. The error bars of the photometric data are included (in most photometric bands, they fall inside the symbols). The low-resolution spectrum acquired in 2021 over the H α region is also displayed (with a dashed black line). The solid red line shows the SED modeled considering the contribution of the photospheric fluxes from both stars, the thermal emission from a gaseous shell close to the system and the effect of the interstellar medium extinction. The solid blue line shows our best-fitting theoretical SED, obtained by adding to the SED plotted in red the contribution of three dusty shells surrounding the stellar system. The flux is normalized to that at λ r e f = 0.55 μ m and displayed on a logarithmic scale. The wavelengths are in microns.
Galaxies 11 00072 g011

5. Discussion

The forest of emission features in the IR spectral ranges presented in this work account for the presence of large amounts of CS gas embedding and veiling MWC 645, which makes it challenging to characterize the components of the binary system. In addition, the existence of dust revealed by the strong IR excess of radiation above the photospheric fluxes is another element to consider when building a probable scenario.
In previous studies, the hot component of the system was considered to be a B-type star [13,15], until Nodyarov et al. [20] assigned it as an early-B subtype. This assignment was mainly constrained by the absence of He II lines and the possible presence of He I absorption lines suggested by Jaschek et al. [15] (see Section 1) in the optical spectral range and the identification of the He I lines in emission in the near-IR. We were not able to detect neutral helium lines in either emission or absorption in our data, at least with a signal intensity above the noise level. However, we cannot discard a possible blend of the He I line at 1.083 μ m (the transition of He I with the highest theoretical intensity that falls in our observed spectral ranges) with a group of Fe II emission lines precisely identified in the same spectral region. We detected Mg II lines in emission at 2.138 μ m and 2.144 μ m. According to the works of Clark and Steele [47] and Steele and Clark [48], who studied a representative sample of Be stars in the H- and K-bands, no evidence of He I features and the simultaneous presence of Br γ and Mg II lines in emission indicate a spectral type between B2 and B4. If He I lines are present, the spectral type is B3 or earlier. New high-resolution IR spectra could be valuable to clarify the presence of He I lines. Nevertheless, we cannot discard any variability in neutral helium lines.
The 12 CO absorption bands detected for the first time at 1.62 μ m and 2.3 μ m allow us to constrain the spectral type and effective temperature of the cool binary component. The spectral type derived by the best fit of the first 12 CO band head with that of a G0-type star does not agree with the strength of metallic lines, which are seen very weakly in the spectra of MWC 645. The effective temperature associated with this spectral type for supergiant/giant stars is around 5600 K [49,50]. Using the EW of the CO(2,0) band head, we determined a T e f f average value of about 5250 K. However, since the SED of the star in the near-IR has a significant contribution from the hot companion and from the ionized envelope (free-free emission), all spectral features from the cool companion are reduced in their intensity (due to the false continuum). Thus, the temperature determination should be interpreted as an upper limit to the effective temperature.
On the other hand, we might consider that instead of tracing a cool companion, the detected CO absorption might also uncover some mass-ejection episode, such as can be seen in the yellow hypergiant ρ Cas [51] or in the eruptive variable V838 Mon [52], when molecular absorption bands start to develop in the spectrum while the star’s brightness fades. In the case of ρ Cas, the CO bands turn from absorption into emission when the star reaches the next maximum brightness phase in which it is too hot for molecule condensation in its extended outer layers [53]. This finding agrees with the scenario that mass was ejected into the environment, radiating while it is expanding and cooling. Since we have only two observations, we cannot trace the evolution of the molecular absorption along the light curve. However, it is interesting to mention that the most intense CO absorption (seen in 2017) occurs when the brightness begins to increase after a seeming local minimum. This observation presents an emission peak blueward of the CO(2,0) absorption band head that is blue-shifted at about 308 km s 1 . If it is indeed due to CO, it might suggest a high-velocity molecular outflow. Furthermore, the considerably weaker CO absorption spectrum observed in 2018 might be interpreted as partially filled with circumstellar emission. Moreover, the absorption-line spectrum ascribed to the cool companion might have originated in an optically thick disk. Polster et al. [54] suggested that the absorption features seen in the spectrum of the FS CMa star, MWC 623, are formed in an equatorial disk viewed nearly edge-on, which acts as a pseudo-photosphere.
When we look at the light curve, we see global qualitative similarities between the range mostly recorded in V-band magnitudes up to the deepest minimum (which occurred in August/September 2018) and that traced by the g-band magnitudes from this minimum up to January 2023. In both time lapses of about four years, we can distinguish a well-outlined dip (at HJD-2450000 ∼ 7500 d and ∼ 9000 d) less intense than the main minimum. Although they present slightly different shapes and depths, and the light curve afterwards reaches a different maximum magnitude level, they are alike. This similar pattern suggests that the dominant source of variability is the same. Variable CS (or circumbinary) extinction along the line of sight due to dust clumps might be responsible for these photometric variations [55,56,57]. Our simple model to fit the SED allowed us to derive the global properties of the dusty envelope. Despite its spherical geometry, it traces different components of optically thin dust, the first located at a distance of ∼6 AU. This distance agrees with the innermost dusty disk radius of ∼5 AU derived for the star FS CMa through aperture-synthesis imaging in the L and N bands [58]. Material orbiting the system at this distance at Keplerian velocity would have a period of ∼5 years. A warped inner edge of the disk can also produce variable extinction [59].
Previous studies of MWC 645 have not reported variations in the V/R ratio of the H α line. Observations made over the last 30 years, although discontinuous in time, have not revealed an inversion in the intensity ratio of the peaks, always showing V/R < 1. We noticed that the V/R ratio derived from our spectra varies from 0.2 to 0.5 (see Table 2). We also calculated the V/R ratio from the intensities of both H α peaks included in the work of Nodyarov et al. [20]. Their observations are from 2004 to 2021 and not continuous. We found a variable V/R ratio with values between 0.3 and 0.9. The only two spectra from dates included in the light curve range, October 2016 (HJD-2450000 ∼ 7680 d) and November 2021 (HJD-2450000 ∼ 9545 d), present a V/R ratio of 0.7 and 0.3, respectively. Spectroscopic monitoring of the H α line could help scan and characterize this variability and search for any periodicity in the changes of the peak intensities that could be related to the rotation of a density perturbation in the disk [60] or an orbital motion [61]. Zickgraf computed line profiles assuming a latitude-dependent wind model with a dust disk and obtained similar line profiles to the observed H α profile with V/R < 1 for an intermediate inclination angle.
The dimming in optical brightness occurs over a long time, and, along this phase, the H α emission changes. The V/R ratio decreases, and the EW increases (with the smallest value of the V/R relation and the maximum EW at the minimum of the light curve). This fact seems to be associated with a change in the amount of the circumstellar material, and not only due to a natural increase in emission intensity during the light curve minimum. Also, the FWHM of the blue emission component is the largest at this point. This H α line strengthening might be attributed to enhanced mass loss (or a mass ejection episode). Similar variations in the emission of the spectral lines and brightness have been observed in the FS CMa star, MWC 728. The light curve of MWC 645 suggests a (possible) periodic behavior. If this were true, the H α enhancement might result from a mass transfer process during periastron passage in an eccentric binary system, with a period of the order of 4 years. A denser observational grid with high-quality spectra is needed to study the link between mass loss and brightness behavior.
A detailed calculation of the H α line profile based on a physically consistent model is a difficult task. However, it would be valuable to explore models with simplified assumptions to gain insight into the system geometry and the structure of the CS matter [9,54]. More complex scenarios considering non-conservative mass transfer between the binary components should be considered to draw a picture of the structures involved in the emission processes [62,63].

6. Conclusions

In this paper, we have studied the FS CMa-type object, MWC 645, a recently confirmed binary system. We have presented IR medium-resolution spectra covering the J-, H-, K-, and L-bands and identified the main spectral features. We have reported the presence of CO bands in absorption for the first time. We have searched for periodicity in the light curve and a possible correlation between its behavior and the spectroscopic optical data. We found that the photometric variations could be explained by variable extinction along the line of sight. In addition, we noted that the stellar brightness fading is accompanied by the enhancement of the H α line emission, which might be due to mass ejection events. Finally, a proper fitting to the observed SED was found, giving a global picture of the gaseous and dusty structures that could enshroud the binary.
Simultaneous optical and near-IR spectroscopy during the following brightness minimum would be very useful for tracing the onset and progress of the possible mass transfer. Such an understanding is utmost for deepening the comprehension of binary evolution in general and of the nature of this fascinating object in particular.

Author Contributions

Conceptualization, A.F.T., M.L.A. and M.K.; methodology, A.F.T., M.L.A. and M.K.; software, A.F.T., M.L.A. and L.V.M.; formal analysis, A.F.T. and M.L.A.; investigation, A.F.T., M.L.A., M.K., L.V.M. and T.E.; resources, A.F.T., M.L.A. and M.K.; data curation, A.F.T., M.L.A., M.K., L.V.M. and T.E.; writing—original draft preparation, review and editing, A.F.T., M.L.A., M.K., L.V.M. and T.E.; visualization, A.F.T., M.L.A. and M.K.; funding acquisition, A.F.T., M.L.A., M.K. and T.E. All authors have read and agreed to the published version of the manuscript.

Funding

A.F.T. and M.L.A. acknowledge financial support from the Universidad Nacional de La Plata (Programa de Incentivos 11/G160) and CONICET (PIP 1337), Argentina. M.K. acknowledges financial support from the Czech Science Foundation (GAČR, grant number 20-00150S). The Astronomical Institute of the Czech Academy of Sciences, Ondřejov, is supported by project RVO: 67985815. This project has received funding from the European Union’s Framework Programme for Research and Innovation Horizon 2020 (2014-2020) under the Marie Skłodowska-Curie Grant Agreement No. 823734. T.E. gratefully acknowledges financial support from the Estonian Ministry of Education and Research through the Estonian Research Council institutional research funding IUT40-1, from the European Union European Regional Development Fund project KOMEET 2014-2020.4.01.16-0029.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data involved in this research are available on request from the authors.

Acknowledgments

The authors are grateful to the referees, whose comments and suggestions helped to improve the paper. This work has made use of IRAF, which is distributed by the National Optical Astronomy Observatory, operated by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation; the BeSS database operated at LESIA, Observatoire de Meudon, France (http://basebe.obspm.fr (accessed on 2 February 2023)); the SIMBAD database and the VizieR catalog access tool, both operated at CDS, Strasbourg, France; the NASA Astrophysics Data System (ADS); the NASA IRTF (Infrared Telescope Facility) Spectral Library; and the NASA IRSA period search tool and the All-Sky Automated Survey for Supernovae (ASAS-SN). This paper is based on observations obtained at (i) Ondřejov Observatory (Czech Republic) with the Perek 2 m telescope; (ii) Tartu Observatory (Estonia); and (iii) the international Gemini Observatory, a program of NSF’s NOIRLab, which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation on behalf of the Gemini Observatory partnership: the National Science Foundation (United States), National Research Council (Canada), Agencia Nacional de Investigación y Desarrollo (Chile), Ministerio de Ciencia, Tecnología e Innovación (Argentina), Ministério da Ciência, Tecnologia, Inovações e Comunicações (Brazil), and Korea Astronomy and Space Science Institute (Republic of Korea) under program IDs GN-2017A-Q-62, GN-2018A-Q-406, and GN-2022B-Q-225. This work has made use of the ground-based research infrastructure of Tartu Observatory, funded through the projects TT8 (Estonian Research Council) and KosEST (EU Regional Development Fund). The authors thank L.S. Cidale for many fruitful discussions.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
CSCircumstellar
CCDCharged-coupled device
EWEquivalent width
VBlue peak intensity
RRed peak intensity
V/RBlue-to-red emission peak ratio
FWHMFull width at half maximum
HJDHeliocentric Julian date
IRInfrared
SEDSpectral energy distribution

References

  1. Allen, D.A.; Swings, J.P. The spectra of peculiar Be star with infrared excesses. Astron. Astrophys. 1976, 47, 293–302. [Google Scholar]
  2. Lamers, H.J.G.L.M.; Zickgraf, F.J.; de Winter, D.; Houziaux, L.; Zorec, J. An improved classification of B[e]-type stars. Astron. Astrophys. 1998, 340, 117–128. [Google Scholar]
  3. Miroshnichenko, A.S. Toward Understanding the B[e] Phenomenon. I. Definition of the Galactic FS CMa Stars. Astrophys. J. 2007, 667, 497–504. [Google Scholar] [CrossRef] [Green Version]
  4. Korčáková, D.; Sestito, F.; Manset, N.; Kroupa, P.; Votruba, V.; Šlechta, M.; Danford, S.; Dvořáková, N.; Raj, A.; Chojnowski, S.D.; et al. First detection of a magnetic field in low-luminosity B[e] stars. New scenarios for the nature and evolutionary stages of FS CMa stars. Astron. Astrophys. 2022, 659, A35. [Google Scholar] [CrossRef]
  5. Kuratova, A.K.; Miroshnichenko, A.S.; Zharikov, S.V.; Manset, N.; Khokhlov, S.A.; Raj, A.; Kusakin, A.V.; Reva, I.V.; Kokumbaeva, R.I.; Usenko, I.A.; et al. Spectroscopic Monitoring of the B[E] Objects Fscma and MOCAM. Odessa Astron. Publ. 2019, 32, 63. [Google Scholar] [CrossRef]
  6. Miroshnichenko, A.S.; Manset, N.; Kusakin, A.V.; Chentsov, E.L.; Klochkova, V.G.; Zharikov, S.V.; Gray, R.O.; Grankin, K.N.; Gandet, T.L.; Bjorkman, K.S.; et al. Toward Understanding the B[e] Phenomenon. II. New Galactic FS CMa Stars. Astrophys. J. 2007, 671, 828–841. [Google Scholar] [CrossRef]
  7. Miroshnichenko, A.S.; Chentsov, E.L.; Klochkova, V.G.; Zharikov, S.V.; Grankin, K.N.; Kusakin, A.V.; Gandet, T.L.; Klingenberg, G.; Kildahl, S.; Rudy, R.J.; et al. Toward Understanding the B[e] Phenomenon. III. Properties of the Optical Counterpart of IRAS 00470+6429. Astrophys. J. 2009, 700, 209–220. [Google Scholar] [CrossRef] [Green Version]
  8. Vink, J.S.; de Koter, A.; Lamers, H.J.G.L.M. Mass-loss predictions for O and B stars as a function of metallicity. Astron. Astrophys. 2001, 369, 574–588. [Google Scholar] [CrossRef]
  9. Carciofi, A.C.; Miroshnichenko, A.S.; Bjorkman, J.E. Toward Understanding the B[e] Phenomenon. IV. Modeling of IRAS 00470+6429. Astrophys. J. 2010, 721, 1079–1089. [Google Scholar] [CrossRef] [Green Version]
  10. Korčáková, D.; Miroshnichenko, A.S.; Zharikov, S.V.; Manset, N.; Danford, S.; Votruba, V.; Dvořáková, N. Mystery of lithium in FS CMa stars. Mem. Soc. Astron. Ital. 2020, 91, 118. [Google Scholar]
  11. Miroshnichenko, A.S.; Zharikov, S.V.; Manset, N.; Khokhlov, S.A.; Nodyarov, A.S.; Klochkova, V.G.; Danford, S.; Kuratova, A.K.; Mennickent, R.; Chojnowski, S.D.; et al. Recent Progress in Finding Binary Systems with the B[e] Phenomenon. Galaxies 2023, 11, 36. [Google Scholar] [CrossRef]
  12. Merrill, P.W.; Burwell, C.G. No. 682. Supplement to the Mount Wilson catalogue and bibliography of stars of classes B and A whose spectra have bright hydrogen lines. Contrib. Mt. Wilson Obs. Carnegie Inst. Wash. 1943, 682, 1–32. [Google Scholar]
  13. Swings, J.P.; Allen, D.A. MWC 645 and MWC 819: Two Stars Resembling Eta Carinae. Astrophys. Lett. 1973, 14, 65. [Google Scholar]
  14. Gottlieb, E.W.; Liller, W. Photometric histories of six infrared objects and three highly reddened blue supergiants. Astrophys. J. 1978, 225, 488–495. [Google Scholar] [CrossRef]
  15. Jaschek, M.; Andrillat, Y.; Jaschek, C. B[e] stars. III. MWC 645. Astron. Astrophys. Suppl. 1996, 120, 99–105. [Google Scholar] [CrossRef] [Green Version]
  16. Landaberry, S.J.C.; Pereira, C.B.; de Araújo, F.X. The eta Carinae spectrum of SS73 11 (=Ve 2–27). Astron. Astrophys. 2001, 376, 917–927. [Google Scholar] [CrossRef] [Green Version]
  17. Zickgraf, F.J. Kinematical structure of the circumstellar environments of galactic B[e] -type stars. Astron. Astrophys. 2003, 408, 257–285. [Google Scholar] [CrossRef]
  18. Marston, A.P.; McCollum, B. Extended shells around B[e] stars. Implications for B[e] star evolution. Astron. Astrophys. 2008, 477, 193–202. [Google Scholar] [CrossRef] [Green Version]
  19. Nodyarov, A.S.; Miroshnichenko, A.S.; Khokhlov, S.A.; Zharikov, S.V.; Manset, N.; Klochkova, V.G.; Usenko, I.A. High-Resolution Spectroscopy of the B[e] Star MWC 645. Odessa Astron. Publ. 2021, 34, 59. [Google Scholar] [CrossRef]
  20. Nodyarov, A.S.; Miroshnichenko, A.S.; Khokhlov, S.A.; Zharikov, S.V.; Manset, N.; Klochkova, V.G.; Grankin, K.N.; Arkharov, A.A.; Efimova, N.; Klimanov, S.; et al. Toward Understanding the B[e] Phenomenon. IX. Nature and Binarity of MWC645. Astrophys. J. 2022, 936, 129. [Google Scholar] [CrossRef]
  21. Elias, J.H.; Joyce, R.R.; Liang, M.; Muller, G.P.; Hileman, E.A.; George, J.R. Design of the Gemini near-infrared spectrograph. In Proceedings of the Ground-Based and Airborne Instrumentation for Astronomy, Orlando, FL, USA, 24–31 May 2006, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series; McLean, I.S., Iye, M., Eds.; SPIE, 2006; Volume 6269, p. 62694C. [Google Scholar] [CrossRef]
  22. Slechta, M.; Skoda, P. 2-meter telescope devices: Coudé slit spectrograph and HEROS. Publ. Astron. Inst. Czechoslov. Acad. Sci. 2002, 90, 1–4. [Google Scholar]
  23. Neiner, C.; de Batz, B.; Cochard, F.; Floquet, M.; Mekkas, A.; Desnoux, V. The Be Star Spectra (BeSS) Database. Astron. J. 2011, 142, 149. [Google Scholar] [CrossRef]
  24. Shappee, B.; Prieto, J.; Stanek, K.Z.; Kochanek, C.S.; Holoien, T.; Jencson, J.; Basu, U.; Beacom, J.F.; Szczygiel, D.; Pojmanski, G.; et al. All Sky Automated Survey for SuperNovae (ASAS-SN or “Assassin”). In Proceedings of the American Astronomical Society Meeting, Washington, DC, USA, 5–9 January 2014; p. 236.03. [Google Scholar]
  25. Kochanek, C.S.; Shappee, B.J.; Stanek, K.Z.; Holoien, T.W.S.; Thompson, T.A.; Prieto, J.L.; Dong, S.; Shields, J.V.; Will, D.; Britt, C.; et al. The All-Sky Automated Survey for Supernovae (ASAS-SN) Light Curve Server v1.0. Publ. ASP 2017, 129, 104502. [Google Scholar] [CrossRef]
  26. Kramida, A.; Ralchenko, Y.; Reader, J. ; NIST ASD Team. NIST Atomic Spectra Database (Ver. 5.10); National Institute of Standards and Technology: Gaithersburg, MD, USA, 2022. Available online: https://physics.nist.gov/asd (accessed on 20 May 2023).
  27. Kleinmann, S.G.; Hall, D.N.B. Spectra of Late-Type Standard Stars in the Region 2.0–2.5 Microns. Astrophys. J. Suppl. 1986, 62, 501. [Google Scholar] [CrossRef]
  28. Mármol-Queraltó, E.; Cardiel, N.; Cenarro, A.J.; Vazdekis, A.; Gorgas, J.; Pedraz, S.; Peletier, R.F.; Sánchez-Blázquez, P. A new stellar library in the region of the CO index at 2.3 μm. New index definition and empirical fitting functions. Astron. Astrophys. 2008, 489, 885–909. [Google Scholar] [CrossRef]
  29. Ali, B.; Carr, J.S.; Depoy, D.L.; Frogel, J.A.; Sellgren, K. Medium-Resolution Near-Infrared (2.15-2.35 micron) Spectroscopy of Late-Type Main-Sequence Stars. Astron. J. 1995, 110, 2415. [Google Scholar] [CrossRef] [Green Version]
  30. Cushing, M.C.; Rayner, J.T.; Vacca, W.D. An Infrared Spectroscopic Sequence of M, L, and T Dwarfs. Astrophys. J. 2005, 623, 1115–1140. [Google Scholar] [CrossRef] [Green Version]
  31. Rayner, J.T.; Cushing, M.C.; Vacca, W.D. The Infrared Telescope Facility (IRTF) Spectral Library: Cool Stars. Astrophys. J. Suppl. 2009, 185, 289–432. [Google Scholar] [CrossRef] [Green Version]
  32. Wallace, L.; Hinkle, K. Medium-Resolution Spectra of Normal Stars in the K Band. Astrophys. J. Suppl. 1997, 111, 445–458. [Google Scholar] [CrossRef] [Green Version]
  33. Heras, A.M.; Shipman, R.F.; Price, S.D.; de Graauw, T.; Walker, H.J.; Jourdain de Muizon, M.; Kessler, M.F.; Prusti, T.; Decin, L.; Vandenbussche, B.; et al. Infrared spectral classification of normal stars. Astron. Astrophys. 2002, 394, 539–552. [Google Scholar] [CrossRef]
  34. Winge, C.; Riffel, R.A.; Storchi-Bergmann, T. The Gemini Spectral Library of Near-IR Late-Type Stellar Templates and Its Application for Velocity Dispersion Measurements. Astrophys. J. Suppl. 2009, 185, 186–197. [Google Scholar] [CrossRef]
  35. Marchiano, P.E.; Cidale, L.S.; Brandi, E.; Muratore, M.F. Spectral Energy Distribution in the symbiotic system BI Cru. Bol. Asoc. Argent. Astron. Plata Argent. 2013, 56, 163–166. [Google Scholar]
  36. Arias, M.L.; Cidale, L.S.; Kraus, M.; Torres, A.F.; Aidelman, Y.; Zorec, J.; Granada, A. Near-infrared Spectra of a Sample of Galactic Unclassified B[e] Stars. Publ. ASP 2018, 130, 114201. [Google Scholar] [CrossRef] [Green Version]
  37. Zorec, J. Distances, Kinematics and Distribution of B[e] Stars in Our Galaxy. In Proceedings of the B[e] Stars, Astrophysics and Space Science Library; Hubert, A.M., Jaschek, C., Eds.; Springer: Dordrecht, The Netherlands, 1998; Volume 233, p. 27. [Google Scholar] [CrossRef]
  38. Muratore, M.F.; Kraus, M.; Liermann, A.; Schnurr, O.; Cidale, L.S.; Arias, M.L. Unveiling the evolutionary phase of B[e] supergiants. Bol. Asoc. Argent. Astron. Plata Argent. 2010, 53, 123–126. [Google Scholar]
  39. Cardelli, J.A.; Clayton, G.C.; Mathis, J.S. The Relationship between Infrared, Optical, and Ultraviolet Extinction. Astrophys. J. 1989, 345, 245. [Google Scholar] [CrossRef]
  40. Lamers, H.J.G.L.M.; Cassinelli, J.P. Introduction to Stellar Winds; Cambridge University Press: Cambridge, UK, 1999. [Google Scholar]
  41. Kurucz, R.L. Model atmospheres for G, F, A, B, and O stars. Astrophys. J. Suppl. 1979, 40, 1–340. [Google Scholar] [CrossRef]
  42. Skrutskie, M.F.; Cutri, R.M.; Stiening, R.; Weinberg, M.D.; Schneider, S.; Carpenter, J.M.; Beichman, C.; Capps, R.; Chester, T.; Elias, J.; et al. The Two Micron All Sky Survey (2MASS). Astron. J. 2006, 131, 1163–1183. [Google Scholar] [CrossRef]
  43. Wright, E.L.; Eisenhardt, P.R.M.; Mainzer, A.K.; Ressler, M.E.; Cutri, R.M.; Jarrett, T.; Kirkpatrick, J.D.; Padgett, D.; McMillan, R.S.; Skrutskie, M.; et al. The Wide-field Infrared Survey Explorer (WISE): Mission Description and Initial On-orbit Performance. Astron. J. 2010, 140, 1868–1881. [Google Scholar] [CrossRef]
  44. Egan, M.P.; Price, S.D.; Kraemer, K.E. The Midcourse Space Experiment Point Source Catalog Version 2.3. Proc. Am. Astron. Soc. Meet. Abstr. 2003, 203, 1301. [Google Scholar]
  45. Neugebauer, G.; Beichman, C.A.; Soifer, B.T.; Aumann, H.H.; Chester, T.J.; Gautier, T.N.; Gillett, F.C.; Hauser, M.G.; Houck, J.R.; Lonsdale, C.J.; et al. Early Results from the Infrared Astronomical Satellite. Science 1984, 224, 14–21. [Google Scholar] [CrossRef]
  46. Yamamura, I. The AKARI Far-Infrared Bright Source Catalogue. In Proceedings of the 38th COSPAR Scientific Assembly, Bremen, Germany, 18–25 July 2010; Volume 38, p. 2. [Google Scholar]
  47. Clark, J.S.; Steele, I.A. A representative sample of Be stars. II. K band spectroscopy. Astron. Astrophys. Suppl. 2000, 141, 65–77. [Google Scholar] [CrossRef] [Green Version]
  48. Steele, I.A.; Clark, J.S. A representative sample of Be stars III: H band spectroscopy. Astron. Astrophys. 2001, 371, 643–651. [Google Scholar] [CrossRef] [Green Version]
  49. Popper, D.M. Stellar masses. Annu. Rev. Astron Astrophys. 1980, 18, 115–164. [Google Scholar] [CrossRef]
  50. Kovtyukh, V.V. High-precision effective temperatures of 161 FGK supergiants from line-depth ratios. Mon. Not. RAS 2007, 378, 617–624. [Google Scholar] [CrossRef] [Green Version]
  51. Gorlova, N.; Lobel, A.; Burgasser, A.J.; Rieke, G.H.; Ilyin, I.; Stauffer, J.R. On the CO Near-Infrared Band and the Line-splitting Phenomenon in the Yellow Hypergiant ρ Cassiopeiae. Astrophys. J. 2006, 651, 1130–1150. [Google Scholar] [CrossRef] [Green Version]
  52. Geballe, T.R.; Rushton, M.T.; Eyres, S.P.S.; Evans, A.; van Loon, J.T.; Smalley, B. Infrared spectroscopy of carbon monoxide in V838 Monocerotis during 2002–2006. Astron. Astrophys. 2007, 467, 269–275. [Google Scholar] [CrossRef]
  53. Kraus, M.; Kolka, I.; Aret, A.; Nickeler, D.H.; Maravelias, G.; Eenmäe, T.; Lobel, A.; Klochkova, V.G. A new outburst of the yellow hypergiant star ρ Cas. Mon. Not. RAS 2019, 483, 3792–3809. [Google Scholar] [CrossRef]
  54. Polster, J.; Korčáková, D.; Manset, N. Time-dependent spectral-feature variations of stars displaying the B[e] phenomenon. IV. V2028 Cygni: Modelling of Hα bisector variability. Astron. Astrophys. 2018, 617, A79. [Google Scholar] [CrossRef] [Green Version]
  55. Miroshnichenko, A.; Ivezić, Ž.; Elitzur, M. On Protostellar Disks in Herbig Ae/Be Stars. Astrophys. J. Lett. 1997, 475, L41–L44. [Google Scholar] [CrossRef] [Green Version]
  56. Herbst, W.; Shevchenko, V.S. A Photometric Catalog of Herbig AE/BE Stars and Discussion of the Nature and Cause of the Variations of UX Orionis Stars. Astron. J. 1999, 118, 1043–1060. [Google Scholar] [CrossRef]
  57. Bouvier, J.; Grankin, K.N.; Alencar, S.H.P.; Dougados, C.; Fernández, M.; Basri, G.; Batalha, C.; Guenther, E.; Ibrahimov, M.A.; Magakian, T.Y.; et al. Eclipses by circumstellar material in the T Tauri star AA Tau. II. Evidence for non-stationary magnetospheric accretion. Astron. Astrophys. 2003, 409, 169–192. [Google Scholar] [CrossRef] [Green Version]
  58. Hofmann, K.H.; Bensberg, A.; Schertl, D.; Weigelt, G.; Wolf, S.; Meilland, A.; Millour, F.; Waters, L.B.F.M.; Kraus, S.; Ohnaka, K.; et al. VLTI-MATISSE L- and N-band aperture-synthesis imaging of the unclassified B[e] star FS Canis Majoris. Astron. Astrophys. 2022, 658, A81. [Google Scholar] [CrossRef]
  59. Bartlett, E.S.; Clark, J.S.; Coe, M.J.; Garcia, M.R.; Uttley, P. Timing and spectral analysis of the unusual X-ray transient XTE J0421+560/CI Camelopardalis. Mon. Not. RAS 2013, 429, 1213–1220. [Google Scholar] [CrossRef] [Green Version]
  60. Okazaki, A.T. Long-Term V/R Variations of Be Stars Due to Global One-Armed Oscillations of Equatorial Disks. Publ. ASJ 1991, 43, 75–94. [Google Scholar]
  61. Kriz, S.; Harmanec, P. A Hypothesis of the Binary Origin of Be Stars. Bull. Astron. Institutes Czechoslov. 1975, 26, 65. [Google Scholar]
  62. Deschamps, R.; Braun, K.; Jorissen, A.; Siess, L.; Baes, M.; Camps, P. Non-conservative evolution in Algols: Where is the matter? Astron. Astrophys. 2015, 577, A55. [Google Scholar] [CrossRef] [Green Version]
  63. Brož, M.; Mourard, D.; Budaj, J.; Harmanec, P.; Schmitt, H.; Tallon-Bosc, I.; Bonneau, D.; Božić, H.; Gies, D.; Šlechta, M. Optically thin circumstellar medium in the β Lyr A system. Astron. Astrophys. 2021, 645, A51. [Google Scholar] [CrossRef]
Figure 1. Normalized medium-resolution spectrum of MWC 645 taken with Gemini/GNIRS on July 2018 from 8500 Å to 13,600 Å. The normalized Ondřejov spectrum from 8400 Å to 8870 Å, acquired on September 2018, is shown in cyan. Main spectral lines are identified by colored markings. The spectral features of a given element (either permitted or forbidden and of different ionization states) are joined by a dashed line of the same color: hydrogen is indicated in red, oxygen in gray, iron in blue, calcium in pink, nitrogen in violet, sulfur in green, and helium in cyan. Wavelengths are given in angstroms.
Figure 1. Normalized medium-resolution spectrum of MWC 645 taken with Gemini/GNIRS on July 2018 from 8500 Å to 13,600 Å. The normalized Ondřejov spectrum from 8400 Å to 8870 Å, acquired on September 2018, is shown in cyan. Main spectral lines are identified by colored markings. The spectral features of a given element (either permitted or forbidden and of different ionization states) are joined by a dashed line of the same color: hydrogen is indicated in red, oxygen in gray, iron in blue, calcium in pink, nitrogen in violet, sulfur in green, and helium in cyan. Wavelengths are given in angstroms.
Galaxies 11 00072 g001
Figure 2. Normalized medium-resolution spectrum of MWC 645 taken with Gemini/GNIRS in 2018, covering the H-(upper panel) and K-bands (lower panel). Main spectral lines and molecular bands are identified by colored markings. The spectral features of a given element (either permitted or forbidden and of different ionization states) or molecule (of different isotopes) are joined by a dashed line of the same color: hydrogen is indicated in red, magnesium in gray, iron in blue, sodium in violet, nitrogen in green, and carbon monoxide in pink. Wavelengths are given in angstroms.
Figure 2. Normalized medium-resolution spectrum of MWC 645 taken with Gemini/GNIRS in 2018, covering the H-(upper panel) and K-bands (lower panel). Main spectral lines and molecular bands are identified by colored markings. The spectral features of a given element (either permitted or forbidden and of different ionization states) or molecule (of different isotopes) are joined by a dashed line of the same color: hydrogen is indicated in red, magnesium in gray, iron in blue, sodium in violet, nitrogen in green, and carbon monoxide in pink. Wavelengths are given in angstroms.
Galaxies 11 00072 g002
Figure 3. Normalized L-band spectrum of MWC 645 obtained in 2022. The emission lines of H I and Fe II are marked in red and blue, respectively. The “bump” longward of the Br α line is a remnant from telluric correction. Wavelengths are in angstroms.
Figure 3. Normalized L-band spectrum of MWC 645 obtained in 2022. The emission lines of H I and Fe II are marked in red and blue, respectively. The “bump” longward of the Br α line is a remnant from telluric correction. Wavelengths are in angstroms.
Galaxies 11 00072 g003
Figure 4. Strongest H I lines detected in our IR spectra of MWC 645. They display an asymmetric profile, where the red flank is steeper than the blue one. Wavelengths are in angstroms.
Figure 4. Strongest H I lines detected in our IR spectra of MWC 645. They display an asymmetric profile, where the red flank is steeper than the blue one. Wavelengths are in angstroms.
Galaxies 11 00072 g004
Figure 5. CO molecular bands of MWC 645. Some emission and absorption lines are also identified. The spectral features of a given element or molecule (of different isotopes) are indicated by colored markings and joined by a dashed line of the same color: hydrogen is indicated in red, magnesium in gray, iron in blue, sodium in violet, nitrogen in green, calcium in cyan, and carbon monoxide in pink. Upper panel: 12 CO second-overtone band heads seen in the H-band spectrum taken in 2018. Lower panel: 12 CO and 13 CO band heads in absorption detected in the K-band. The spectra obtained in 2017 (in red) and 2018 (in black) revealed the variability in the strength of the observed bands. Wavelengths are given in angstroms.
Figure 5. CO molecular bands of MWC 645. Some emission and absorption lines are also identified. The spectral features of a given element or molecule (of different isotopes) are indicated by colored markings and joined by a dashed line of the same color: hydrogen is indicated in red, magnesium in gray, iron in blue, sodium in violet, nitrogen in green, calcium in cyan, and carbon monoxide in pink. Upper panel: 12 CO second-overtone band heads seen in the H-band spectrum taken in 2018. Lower panel: 12 CO and 13 CO band heads in absorption detected in the K-band. The spectra obtained in 2017 (in red) and 2018 (in black) revealed the variability in the strength of the observed bands. Wavelengths are given in angstroms.
Galaxies 11 00072 g005
Figure 6. Comparison between the degraded spectrum of MWC 645 to R ∼ 2000 (solid black line) and a G0 Ib-II star, HD 185018 (dashed red line), where the CO(2-0) band head fits well. The other 12 CO absorption bands from MWC 645 are shallower than those of the template; perhaps they are filled by emission. Wavelengths are given in angstroms.
Figure 6. Comparison between the degraded spectrum of MWC 645 to R ∼ 2000 (solid black line) and a G0 Ib-II star, HD 185018 (dashed red line), where the CO(2-0) band head fits well. The other 12 CO absorption bands from MWC 645 are shallower than those of the template; perhaps they are filled by emission. Wavelengths are given in angstroms.
Galaxies 11 00072 g006
Figure 7. Light curve of MWC 645 from 16 December 2014 up to 11 November 2022 taken from ASAS-SN. The purple squares indicate the V-band magnitudes, and the green circles represent the g-band measurements converted to the V-band magnitudes. The conversion has been carried out with the relationship found by Nodyarov et al. [20]. Their optical photometric observations are also included (blue circles). Vertical solid blue lines mark the dates of our IR observations (2017, 2018, and 2022, respectively); dotted blue lines mark the dates of our optical spectra (2018 and 2021, respectively); and those dotted in green, gray, and red correspond to the spectra downloaded from the BeSS database taken in 2019, 2020, and 2022, respectively. Time is given in heliocentric Julian dates (HJD) minus 2.45 × 10 6 days.
Figure 7. Light curve of MWC 645 from 16 December 2014 up to 11 November 2022 taken from ASAS-SN. The purple squares indicate the V-band magnitudes, and the green circles represent the g-band measurements converted to the V-band magnitudes. The conversion has been carried out with the relationship found by Nodyarov et al. [20]. Their optical photometric observations are also included (blue circles). Vertical solid blue lines mark the dates of our IR observations (2017, 2018, and 2022, respectively); dotted blue lines mark the dates of our optical spectra (2018 and 2021, respectively); and those dotted in green, gray, and red correspond to the spectra downloaded from the BeSS database taken in 2019, 2020, and 2022, respectively. Time is given in heliocentric Julian dates (HJD) minus 2.45 × 10 6 days.
Galaxies 11 00072 g007
Figure 8. Fourier power spectrum of the ASAS-SN light curve of MWC 645. The periods in days are on a logarithmic scale. The green line shows the confidence level.
Figure 8. Fourier power spectrum of the ASAS-SN light curve of MWC 645. The periods in days are on a logarithmic scale. The green line shows the confidence level.
Galaxies 11 00072 g008
Figure 9. H α line variation of MWC 645. Spectra taken in 2018, 2019, 2020, and 2022 are displayed in blue, green, gray, and red, respectively. The spectrum with the lowest resolution is plotted with a dashed-line. The heliocentric radial velocity scale is shown in km s 1 .
Figure 9. H α line variation of MWC 645. Spectra taken in 2018, 2019, 2020, and 2022 are displayed in blue, green, gray, and red, respectively. The spectrum with the lowest resolution is plotted with a dashed-line. The heliocentric radial velocity scale is shown in km s 1 .
Galaxies 11 00072 g009
Figure 10. Example of the emission lines in the surroundings of the H α line of MWC 645. Ondřejov normalized spectra taken in 2018 on 12 September (in red line) and 13 September (in black line) and the 2019 low-resolution spectrum (in green) are shown with the main lines identified by colored markings. The spectral features of a given element (either permitted or forbidden and of different ionization states) are joined by a dashed line of the same color: hydrogen is indicated in red, oxygen in gray, iron in blue, nitrogen in violet, and sulfur in green. Wavelengths are given in angstroms.
Figure 10. Example of the emission lines in the surroundings of the H α line of MWC 645. Ondřejov normalized spectra taken in 2018 on 12 September (in red line) and 13 September (in black line) and the 2019 low-resolution spectrum (in green) are shown with the main lines identified by colored markings. The spectral features of a given element (either permitted or forbidden and of different ionization states) are joined by a dashed line of the same color: hydrogen is indicated in red, oxygen in gray, iron in blue, nitrogen in violet, and sulfur in green. Wavelengths are given in angstroms.
Galaxies 11 00072 g010
Table 1. Mean values of the S/N ratio for GNIRS near-IR observations in different spectral ranges.
Table 1. Mean values of the S/N ratio for GNIRS near-IR observations in different spectral ranges.
ObservationsSpectral RangeMean S/N Ratio
Program ID[Å]
GN-2017A-Q-6222,570–24,400200
GN-2018A-Q-40621,000–24,500200
GN-2018A-Q-40615,800–18,200100
GN-2018A-Q-40611,000–13,60090
GN-2018A-Q-4068400–11,00060
GN-2022B-Q-22533,311–36,446200
GN-2022B-Q-22538,500–41,786100
Table 2. H α line parameters of MWC 645. Column 1 indicates the observing date; column 2 the heliocentric Julian date (minus 2.45 × 10 6 d); column 3 the observatory/database where the spectrum was obtained; columns 4 and 5 the emission intensities of the blue and red peaks (V and R) in continuum units, respectively; column 6 the emission intensity ratio of both peaks (V/R); and column 7 the total equivalent width (EW) in Å. The measurement errors are on the order of 1% for the intensities and 10% for the EW values.
Table 2. H α line parameters of MWC 645. Column 1 indicates the observing date; column 2 the heliocentric Julian date (minus 2.45 × 10 6 d); column 3 the observatory/database where the spectrum was obtained; columns 4 and 5 the emission intensities of the blue and red peaks (V and R) in continuum units, respectively; column 6 the emission intensity ratio of both peaks (V/R); and column 7 the total equivalent width (EW) in Å. The measurement errors are on the order of 1% for the intensities and 10% for the EW values.
Obs. DateHJD-2450000ObservatoryVRV/REW
(yyyy-mm-dd) [Å]
2018-09-118373.3290Ondřejov118.4452.00.26−1980.3
2018-09-128374.2823Ondřejov132.8509.00.26−2122.8
2019-08-288724.3870BeSS 190.2268.90.33−1132.2
2019-08-308726.4279BeSS 275.1391.90.45−2130.7
2019-09-018728.4429BeSS95.7283.90.34−1237.2
2019-09-028729.4419BeSS105.7316.00.33−1291.7
2019-09-038730.4365BeSS101.5306.10.33−1353.7
2019-09-058732.4433BeSS100.1292.30.34−1147.3
2019-09-078734.3942BeSS109.7323.50.34−1372.1
2019-09-108737.3883BeSS99.8296.00.34−1292.5
2019-09-118738.3925BeSS132.3397.10.33−1642.3
2019-09-148741.4241BeSS110.5330.90.33−1348.1
2019-09-198746.3989BeSS100.5297.70.34−1129.5
2019-09-208747.3862BeSS103.1304.30.34−1264.6
2019-09-248751.3722BeSS109.1326.10.33−1283.5
2020-06-259025.5535BeSS94.1223.70.42−933.8
2020-07-129043.5126BeSS127.4264.80.48−1216.8
2021-11-019520.3161Tartu
2022-07-269787.5322BeSS47.3195.60.24−718.2
2022-12-099923.2775BeSS100.5321.20.31−1392.3
1 For details about the instruments and observers, please visit the web page BeSS database (http://basebe.obspm.fr/basebe/) (accessed on 2 February 2023). 2 Spectrum with the lowest resolution.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Torres, A.F.; Arias, M.L.; Kraus, M.; Mercanti, L.V.; Eenmäe, T. New Insight into the FS CMa System MWC 645 from Near-Infrared and Optical Spectroscopy. Galaxies 2023, 11, 72. https://doi.org/10.3390/galaxies11030072

AMA Style

Torres AF, Arias ML, Kraus M, Mercanti LV, Eenmäe T. New Insight into the FS CMa System MWC 645 from Near-Infrared and Optical Spectroscopy. Galaxies. 2023; 11(3):72. https://doi.org/10.3390/galaxies11030072

Chicago/Turabian Style

Torres, Andrea Fabiana, María Laura Arias, Michaela Kraus, Lorena Verónica Mercanti, and Tõnis Eenmäe. 2023. "New Insight into the FS CMa System MWC 645 from Near-Infrared and Optical Spectroscopy" Galaxies 11, no. 3: 72. https://doi.org/10.3390/galaxies11030072

APA Style

Torres, A. F., Arias, M. L., Kraus, M., Mercanti, L. V., & Eenmäe, T. (2023). New Insight into the FS CMa System MWC 645 from Near-Infrared and Optical Spectroscopy. Galaxies, 11(3), 72. https://doi.org/10.3390/galaxies11030072

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop