Open Access
Issue
A&A
Volume 688, August 2024
Article Number A161
Number of page(s) 16
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/202348225
Published online 15 August 2024

© The Authors 2024

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1. Introduction

With the advent of wide-field surveys over the last two decades, new classes or subgroups of transients have been emerging, necessitating a revision of our understanding of the final stages of stellar evolution. One such grouping occupies a range of luminosities between those typical of supernovae (SNe) and classical novae, and are usually referred to as intermediate-luminosity optical transients (Berger et al. 2009) or gap transients (e.g. Kasliwal 2012). Here, we focus on these and the distinction between luminous red novae (LRNe) and intermediate-luminosity red transients (ILRTs).

Luminous red novae have a number of distinctive features. Ranging in luminosity from −4 to −15 mag, they display double-peaked light curves (e.g. Blagorodnova et al. 2017; Pastorello & Fraser 2019; Pastorello et al. 2019a). The initial peak is thought to result from the cooling emission of the common envelope after its ejection in a stellar-merger scenario (Metzger & Pejcha 2017; Blagorodnova et al. 2021; Matsumoto & Metzger 2022). Around the time of the first peak, the spectra are rather blue and display strong Balmer and Fe II emission lines (Cai et al. 2019, 2022a). Later, there is a second peak that is often attributed to matter ejection or hydrogen recombination (Pastorello & Fraser 2019; Matsumoto & Metzger 2022). By the time of the second peak, the Balmer lines are weaker and a forest of metal lines (mostly in absorption) appears (Pastorello et al. 2019a). At late times, the spectra resemble those of an M-type star, displaying molecular absorption bands such as TiO and VO (Martini et al. 1999; Lynch et al. 2004; Pastorello & Fraser 2019; Pastorello et al. 2021). The strongest evidence for the stellar-merger origin of LRNe comes from V1309 Sco, which showed a decrease in period consistent with an inspiralling binary (Tylenda et al. 2011).

On the other hand, ILRTs usually rise to peak over about two weeks (though this is based on a limited sample of about ten objects), and their peak absolute magnitudes tend to fall between −11.5 and −14.5 (e.g. Cai et al. 2018, 2021; Stritzinger et al. 2020a). They display a post-peak plateau, somewhat reminiscent of type IIP SNe. ILRTs show calcium emission features at all stages of their evolution: in particular, the [Ca II] λλ7291,7323 doublet, characteristic of ILRTs, and the λλλ8498,8542,8662 near-infrared (NIR) Ca II triplet. There are several proposed mechanisms for the production of ILRTs, including electron-capture SNe from super-asymptotic giant branch star progenitors (Thompson et al. 2009; Pumo et al. 2009), massive stars in dusty cocoons experiencing outbursts (Bond et al. 2009; Humphreys et al. 2011), and accretion onto a compact object (Kashi et al. 2010).

Luminous red novae and ILRTs can be said to belong to the same broad family of red transient phenomena, and as such it can be difficult to disentangle the two observationally. For example, the [Ca II] λλ7291,7323 doublet, characteristic of ILRTs, is also seen in some LRNe (Cai et al. 2019). Furthermore, the forest of metal lines characteristic of LRNe can also be seen in ILRTs (e.g. NGC 300-2008-OT1; Humphreys et al. 2011). The ILRTs SN 2008S and NGC 300-2008-OT1 each had a detected progenitor that likely did not survive the explosive event (Botticella et al. 2009; Adams et al. 2016; Valerin et al., in prep.), supporting the idea that such events are terminal.

AT 2022fnm ( α = 12 h 25 m 33 . s 713 $ \alpha=12^{\mathrm{h}}25^{\mathrm{m}}33{{{{\overset{\text{s}}{.}}}}}713 $, δ = + 45 ° 41 12 . 45 $ \delta =+45\circ41{{\prime}}12{{{\overset{\prime\prime}{.}}}}45 $, J2000; see Fig. 1) was first discovered on 31 March 2022 (UTC dates are used throughout this paper) by the Asteroid Terrestrial-impact Last Alert System (ATLAS) survey (Tonry et al. 2018; Smith et al. 2020) and was given the internal survey designation ATLAS22kez (Tonry et al. 2022). It was classified on 2 April 2022 as a luminous blue variable by Bostroem et al. (2022) as part of the Global Supernova Project (Howell & Global Supernova Project 2017; Howell 2019).

thumbnail Fig. 1.

B-band image taken at the NOT (ALFOSC) at +3 d showing the location of AT 2022fnm ( α = 12 h 25 m 33 . s 713 $ \alpha=12^{\mathrm{h}}25^{\mathrm{m}}33{{{{\overset{\text{s}}{.}}}}}713 $, δ = + 45 ° 41 12 . 45 $ \delta =+45\circ41{{\prime}}12{{{\overset{\prime\prime}{.}}}}45 $, J2000) in its host galaxy, NGC 4389.

The region of the sky containing AT 2022fnm was well sampled by the ATLAS survey, which allowed us to constrain the epoch of the onset of the outburst to within a 3 d window. We set the outburst epoch as MJD 59 667.98 ± 1.5 based on the midpoint between the first ATLAS detection (MJD 59 669.47) and the last non-detection (MJD 59 666.50, with a limiting magnitude of 20.17 in the ATLAS c band).

AT 2022fnm is associated with the galaxy NGC 4389, for which there is a spread of distances reported in the literature. We adopted a luminosity distance of 10.3 Mpc, as determined from the measured redshift of 0.00239 for a Hubble constant of 69.6 km s−1 Mpc−1 (Bennett et al. 2014). This results in a distance modulus of μ = 30.1 ± 0.6 mag. This distance falls roughly between the recent Tully-Fisher distance from the Cosmicflows4 catalogue (μ = 29.855 ± 0.40 mag; Tully et al. 2013) and the kinematic distance after correcting for local velocity fields, μ = 30.70 ± 0.15 mag, from the NASA Extragalactic Database1(NED). We take AT 2022fnm to have a foreground extinction of AV = 0.040 mag from Schlafly & Finkbeiner (2011, from NED), which assumes an RV of 3.1 (Schultz & Wiemer 1975); motivated by the blue spectra and unremarkable colour evolution, we assume negligible host-galaxy extinction.

This paper is organised as follows: in Section 2 we discuss the acquisition and reduction of the optical and infrared data. In Section 3 we present the photometric evolution of the object and compare it to other gap transients; a full list can be found in Table A.1. In Section 4 we address the spectroscopic evolution. Section 5 contains our overall discussion. The data logs, both photometric and spectroscopic, along with the details of the objects used for comparison can be found in Appendix A.

2. Observational data

2.1. Optical and infrared imaging

We obtained optical imaging from a number of telescope and instrument configurations: Nordic Optical Telescope (NOT) + ALFOSC (as part of the NUTS2 programme2), Liverpool Telescope (LT) + IO:O, Gran Telescopio Canarias (GTC) + OSIRIS, Asiago Schmidt 67/92 + Moravian, and the Las Cumbres Observatory Global Telescope (LCOGT) network of 1m telescopes. A detailed log is given in Table A.2. The images were reduced in a standard manner, comprising trimming, bias subtraction and flat-field correction. In the case of the NOT+ALFOSC and Asiago data, we performed these reductions with the FOSCGUI pipeline3 developed by E. Cappellaro. The IO:O data were reduced using the IO:O pipeline4. The Las Cumbres Observatory data were automatically reduced by the BANZAI pipeline5. Imaging was also obtained from the ATLAS and Zwicky Transient Facility (ZTF; Bellm et al. 2019) surveys. These data are reduced automatically using bespoke pipelines (Smith et al. 2020; Masci et al. 2019; Magnier et al. 2020). All data were taken in Sloan ugriz, Johnson-Cousins BV, and ATLAS co filter systems; the ATLAS o band is roughly equivalent to r+i and ATLAS c to g+r.

We also obtained five epochs of NIR imaging using NOT+NOTCam; the log is given in Table A.3. The images were reduced using the NOTCam QUICKLOOK reduction package6. A master flat was created using bright and dim sky flats and bad pixels were masked. The wide-field camera of NOTCam suffers significant optical distortion for which we accounted using a distortion model. Sky images were created from the dithered on-source exposures, which allowed us to subtract the dark current along with the sky background. For both optical and NIR imaging, we performed point-spread-function-fitting photometry using the SuperNOva PhotometrY (SNOoPY) package7.

2.2. Optical spectroscopy

We obtained ten epochs of long-slit optical spectra spanning +2 d to +90 d from a number of different telescopes; the observation log can be seen in Table A.4.

We reduced the NOT+ALFOSC spectra using the FOSCGUI data reduction pipeline, whilst the Keck+LRIS spectrum was reduced using the LPIPE data-reduction pipeline (Perley 2019). The Telescopio Nazionale Galileo (TNG)+DOLORES and GTC+OSIRIS spectra were reduced manually using standard IRAF8 tasks. The reduction process was similar for each instrument, with the two-dimensional (2D) spectra trimmed, bias-subtracted, and flat-fielded before being extracted to one dimension. Arc-lamp spectra were used to calibrate the wavelengths of the spectra and a master sensitivity function was created from a standard-star spectrum that was taken with an identical instrumental setup, before a telluric correction was applied.

3. Photometric evolution

3.1. Light curve and colour

The multi-band light curve of AT 2022fnm is shown in Fig. 2. It exhibits a relatively smooth evolution in all bands, with a plateau between approximately +15 and +50 d in the redder bands.

thumbnail Fig. 2.

Multi-band light curve of AT 2022fnm. The inverted open triangles represent upper limits. The vertical line marks the epoch of the onset of the outburst. ATLAS data beyond +40 d are also shown with open symbols as they are noisier.

As the o-band limit is the closest to the discovery epoch, we used it to calculate the rise time by fitting it with a second-order polynomial; the resulting rise time is 5.3 ± 1.5 d, which is relatively well constrained compared to the limited sample of ILRTs with secure rise times, and we find a peak apparent magnitude of 17.36 mag in o, corresponding to an absolute magnitude of −12.7 ± 0.1. Overall, the light curve of AT 2022fnm is also reminiscent of LRNe, though it only displays a relatively short plateau rather than a second peak.

We compare the r-band light curve of AT 2022fnm with ILRTs and LRNe in Fig. 3 (references in A.1). For the first 100 d, the light-curve evolution of SN 2008S does not differ markedly from that of AT 2022fnm, but the rate of decline is slower and it does not display the same clear peak and plateau that is visible in the AT 2022fnm light curve. AT 2018aes (Cai et al. 2021) has a similar peak to AT 2022fnm, but its behaviour over the first 40 d after peak does not resemble AT 2022fnm at all, instead showing much slower evolution, but after +40 d their decline rates are quite similar. PTF 10fqs, albeit fainter than AT 2022fnm, shows a similar evolution, though it has a clear plateau that lasts until about +70 d.

thumbnail Fig. 3.

r-band comparison of AT 2022fnm with other ILRTs and LRNe from the literature. All objects have been corrected for reddening with the values listed in Table A.1 and respective references. The reference epoch for each object has been set to its first peak.

In Fig. 4 we show the colour evolution of AT 2022fnm in B − V, g − r, and V − Ks; all display a clear redward progression with time following a brief blueward excursion in the first ten days for B − V and g − r (we do not have V − Ks data covering this period). AT 2021afy, a LRN, shows qualitatively similar behaviour to AT 2022fnm in B − V, with a particularly close match in evolution from about +40 d. On the other hand, SN 2008S, the prototypical ILRT, also shows similar evolution in this colour. The limited data for AT 2022fnm in V − Ks make comparison more difficult, but in general, we see a redward movement in both ILRTs and LRNe, though the rate varies from object to object. SN 2008S, despite its slightly less steep redward movement, is qualitatively the best match to AT 2022fnm.

thumbnail Fig. 4.

Colour evolution of AT 2022fnm and comparison objects. We calculated the colours with observations taken within 0.3 d of each other except in the case of V − Ks, for which we allowed a ±3.5 d window, owing to the modest cadence of the NIR observations.

In Fig. 5 we show both the rise time of AT 2022fnm and our comparison objects and their decline rates over the first 15 d and over the period between +50 d and +100 d (where applicable). The purpose of this plot was to determine where AT 2022fnm fell in the parameter space and whether or not it showed a particular kinship with either LRNe or ILRTs. At least in the case of our limited sample, it seems that LRNe occupy a larger portion of the parameter space both in terms of rise times and decline rates, whilst the ILRTs are clustered together. This is consistent with the findings of Cai et al. (2022b). AT 2022fnm occupies a region consistent with both LRNe and ILRTs.

thumbnail Fig. 5.

Rise time and decline rates plotted against the (first) r-band peak for AT 2022fnm and comparison objects. Rise times plotted with filled markers are calculated against explosion or outburst epoch, whilst those with unfilled markers are calculated against the detection epoch. Rise times are calculated to the peak observed magnitude, which could be more sensitive to photometric uncertainties for objects whose light curve flattens around the peak.

3.2. Pseudo-bolometric light curve and spectral energy distribution evolution

In Fig. 6 we show the pseudo-bolometric light curve of AT 2022fnm along with evolution of the corresponding black-body temperature and radius. The light curve was computed based on our uBgcVroizJHKs data using the SUPERBOL package9 (Nicholl 2018). Owing to a gap in the data coverage, we only show the first 80 d (cf. Fig. 2). After this point the interpolation between epochs and extrapolation across different filters become unreliable, and the pseudo-bolometric light curve is dominated by numerical artefacts. We see that the black-body radius first decreases (around +5 d), before increasing again at around +15 d. Hydrogen recombination could increase the transparency, but then another energy source would be required to stop and reverse this behaviour. As such, a source of ionising photons is required to push the photosphere back out; this will be addressed in Section 5. The photospheric velocities evident in the bottom panel of Fig. 6 are 4000 km s−1 and higher in the first week, whilst the spectra do not suggest velocities this high. However, the spectral lines in emission form farther out in an optically thin region of the circumstellar material (CSM), whereas the velocities derived from the black-body fitting reflect an optically thick continuum that lies at a smaller radius.

thumbnail Fig. 6.

Pseudo-bolometric properties. Top panel: pseudo-bolometric light curve of AT 2022fnm. Second panel: temperature versus time of a black body that has been fitted to the pseudo-bolometric light curve. Third panel: log of the radius versus time of a black body that has been fitted to the pseudo-bolometric light curve. Bottom panel: expansion velocity (radius over time) versus time of a black body that has been fitted to the pseudo-bolometric light curve. Note that the radius is simply being divided by the elapsed time at each point to generate this crude measure. The data have been corrected for redshift and foreground extinction. Points corresponding to the SED fits from Fig. 8 have been overplotted for comparison in the form of red triangles. The figure has been truncated at +80 d as the light curve is dominated by numerical artefacts beyond this epoch.

We attempted to reproduce the pseudo-bolometric light curve with a simple analytical diffusion model that accounts for the luminosity of the cooling ejecta and the 56Ni decay process. Details will be presented in Valerin et al. (in prep.), but the fundamental design is based on Chatzopoulos et al. (2012) and Arnett (1980, 1982). The resulting model light curve for AT 2022fnm is shown in Fig. 7. To mimic both the sharp first peak as well as the subsequent slow decline, we divided the ejecta into two regions. The outer layers (which can be identified with the stellar envelope) are faster and less dense, and therefore release their internal energy over a short timescale. On the other hand, the inner parts of the ejecta contain the bulk of the mass and are slower, resulting in a longer photon diffusion time. Hydrogen recombination, though not directly treated, is accounted for in a rough manner by means of the so-called diffusion time in the bulk material. The light curve is well reproduced by the model, with both the peak and the decline captured.

thumbnail Fig. 7.

Attempt to reproduce the pseudo-bolometric light curve with a simple analytical diffusion model. All three components of the model are shown along with their total, which fits the observed luminosity well. The model is somewhat degenerate, and the details of this are discussed in the text, as are the parameters. The best-fitting explosion epoch from our modelling is MJD 59 669.38, which is consistent with our adopted outburst epoch (MJD 59 667.98 ± 1.5).

We note that the mass and the velocity of the ejecta are coupled, which means that an equivalent model could be contrived with a higher mass and velocity or, equally, with a lower mass and velocity, so we cannot use this model to constrain the energy of the outburst reliably. Nevertheless, we were able to reproduce the pseudo-bolometric light curve with a plausible set of parameters. The model considers the integrated radiated energy of the bulk material and the envelope separately, but combined they give a total integrated radiated energy of Etot = 4.24 × 1048 erg. The best-fitting explosion epoch from our modelling is MJD 59 669.38, which agrees within the uncertainties with our adopted outburst epoch (MJD 59 667.98 ± 1.5). Overall, the model is consistent with a weak eruption or explosion of a star of relatively low mass. The full set of model parameters can be found in Table 1. More sophisticated modelling may be used in the future to provide more robust explosion parameters.

Table 1.

Analytical diffusion model parameters.

In Fig. 8 we show a sequence of spectral energy distribution (SED) fits for AT 2022fnm. It can be seen that at late times the SED is best fit by two black bodies, owing to the development of a NIR excess. The NIR photometric data are sparse, though it does seem that the parameters of the second black body stay fairly constant between +104 d and +120 d.

thumbnail Fig. 8.

Sequence of SED fits for AT 2022fnm. At late times a second black-body component (red) is required, and the sum of the two components is shown in green.

4. Spectral evolution

The spectral sequence for AT 2022fnm is shown in Fig. 9. It initially displays a strong blue continuum that weakens with time, and in the later epochs the continuum clearly becomes red. There is a large decrease in relative luminosity bluewards of Hβ between the +13 d and +23 d spectra. The development of such a red continuum is common to both ILRTs and LRNe.

thumbnail Fig. 9.

Optical spectral sequence of AT 2022fnm. The spectra have been calibrated against the photometry. The offsets applied to each spectrum are shown in the plot, and the +90 d spectrum has been multiplied by a factor of 2 for clarity. Negative fluxes have been set to zero. The +2 d spectrum is plotted in grey with a version binned to 10 Å overplotted in black to compensate for the noise of the spectrum. The absorption features seen at approximately 6725 Å and 8200 Å are artefacts.

Balmer lines dominate the spectra throughout the evolution of AT 2022fnm. Balmer lines are seen in both ILRTs and LRNe, though at late times the lines typically decrease in prominence significantly in LRNe. The complete evolution of Hα and Hβ is plotted in Fig. 10 with the spectra having been binned to a common sampling.

thumbnail Fig. 10.

Evolution of Hα and Hβ in AT 2022fnm over time. All spectra have been binned to a common pixel size and normalised against the peak of Hα in the left panel and Hβ in the right. Only spectra with both Hα and Hβ coverage are shown. Vertical lines marking ±1000, 2000, and 3000 km s−1 have been added to guide the eye.

It can be seen in Fig. 10 that the Hα emission feature narrows with time and that Hβ shows similar behaviour. An absorption line is present bluewards of Hα at 6485 Å in the +33 d spectrum, which we associate with Hα. A similar absorption is present bluewards of Hβ at 4802 Å. Both have velocities of approximately −3500 km s−1. In the +33 d spectrum, two absorption features can be seen at 6574 Å and 6592 Å, just redwards of Hα, and if they are associated with it, then they have velocities of approximately +500 km s−1 and +1300 km s−1, respectively.

There is a noticeable Hα absorption trough in the +2 d spectrum at approximately 6559 Å, corresponding to a velocity of around −150 km s−1 and with a full width at half maximum of approximately 240 km s−1. This feature is not present in the +3 d or later spectra. If this feature is real, then it must result from a thin shell of material in the immediate vicinity of the outburst that was quickly overrun.

To assess the evolution of the Balmer lines, we fitted a Lorentzian profile at each epoch to adequately capture the electron-scattering wings. We find that over the first 15 d after the onset of the outburst, the full width at half maximum velocity stays at approximately 625–750 km s−1 but drops monotonically to around 400 km s−1 at +90 d. In Fig. 11 we compare one early (+8 d) and one late spectrum (+90 d) of AT 2022fnm to spectra of other objects. In the upper-left inset of each panel, we show a zoom-in of the Hα region for each object. In the upper-left inset of panel a, the Hα profile of AT 2022fnm resembles that of AT 2017be and SN 2008S in both shape and velocity, though the peak of the latter appears to be slightly redshifted. In the upper-left inset of panel b, Hα appears to be slightly blueshifted at late times in the case of AT 2022fnm and AT 2017be, but slightly redshifted in the case of SN 2008S and perhaps AT 2018hso, though the signal-to-noise ratio of the line is not particularly good for that object. At early times all of the objects (AT 2022fnm, SN 2008S, AT 2017be, AT 2018aes, and AT 2018hso) are dominated by the Balmer lines, as is typical in ILRTs and LRNe. The overall spectral evolution of AT 2022fnm resembles that of the ILRTs AT 2017be and AT 2018aes. The LRN AT 2018hso resembles AT 2022fnm at early times, but at late times it displays a red continuum and there is little resemblance between the two objects in terms of spectral lines, firmly placing AT 2022fnm in ILRT territory, as far as its spectra are concerned.

thumbnail Fig. 11.

Comparison of an early spectrum and late spectrum of AT 2022fnm with other objects. The spectra have been normalised against the sum of the flux densities except in the case of the insets, which have been normalised against the peak of Hα. The phases are calculated relative to the epoch of the onset of the outburst. The upper-left inset shows Hα. The +83 d AT 2018hso spectrum has been trimmed to 4500–9000 Å and binned to 10 Å, owing to noise.

Moving on from the Balmer lines to consider oxygen, we note that there is the narrow absorption at around 7770 Å, which is visible in the +33 d spectrum as well as the preceding spectra. If this line is O Iλ7774, then it is unusual to see it in absorption; it does not appear to be visible in galactic emission in the 2D spectra, so it is unclear whether this is a result of subtraction. Typically, the O Iλ7774 line is seen in emission in airglow spectra. O I emission at 8446 Å is apparent by +57 d, growing to much greater prominence by +90 d. This O I emission is not pronounced in SN 2008S or AT 2017be at early times, but is clear at later times in these objects, as can be seen in Fig. 11.

It is difficult to confirm the presence of Fe II above the noise in the +2 d spectrum, but it is clearly present from the +3 d spectrum all the way up through the +90 d spectrum. In particular, there are emission lines at 4924 Å, 5018 Å, and 5169 Å, as well as a forest of emission lines between approximately 4490 Å and 4620 Å.

Turning to calcium, initially weak, we note that the Ca IIλ8498,8542,8662 NIR triplet becomes apparent within the first two weeks and is extremely pronounced at late times. The triplet is also weak at early times in AT 2017be (Fig. 11), but, as in AT 2022fnm, it is very prominent at later times. The [Ca II] λλ7291,7323 doublet, characteristic of ILRTs, though present from the +13 d spectrum onwards, is unusually weak for an ILRT. Like AT 2022fnm, the late-time spectrum of the ILRT AT 2017be does not show particularly prominent emission in the [Ca II] λλ7291,7323 doublet (Fig. 11). The doublet, a key signature of ILRTs, is not typically seen in LRNe, although it was observed in AT 2018hso as a very weak emission feature (Cai et al. 2019). The Ca H λ3969 and Ca K λ3934 lines are initially very prominent in absorption, but by +57 d they are instead in emission. The change from Ca H and K absorption to emission can be seen in Fig. 11 where Ca H and K emission is also seen in the ILRTs SN 2008S and AT 2018aes as well as in the LRN AT 2018hso.

The late-time spectra are reminiscent of interacting SNe, although we note that many CSM-dominated transients appear similar at these phases. We see no evidence for spectral features associated with nucleosynthesis (e.g. [O II] λλ6300,6364), although these can be suppressed in interacting core-collapse SNe due to higher densities. O Iλ8446 is present, and it is possible that this is due to Bowen fluorescence (Bowen 1947). In this case, the Lyman-β photons required presumably arise from CSM interaction. Given that the Ca IIλ8498,8542,8662 NIR triplet and that lines due to Fe II are present, whilst the [Ca II] λλ7291,7323 doublet is very weak (Fig. 9) and [Fe II] lines are absent, suggests that the density remained high (Kozma & Fransson 1998; Dessart & Hillier 2011).

5. Discussion

The similarity of the light curve of the ILRT AT 2022fnm with light curves of LRNe reflects that in both cases a comparable amount of energy is injected into a comparable amount of hydrogen-rich material. The object displays a clear NIR excess (apparent by approximately +100 d), typical of ILRTs (Cai et al. 2021) and also seen in a number of LRNe (Pastorello et al. 2019b; Stritzinger et al. 2020b). On the other hand, the object does not display the double-peaked light curve associated with LRNe, though this feature is not particularly prominent in some LRNe (e.g. AT 2018hso; Cai et al. 2019). Nor does AT 2022fnm display the overwhelming forest of metal lines in absorption along with the weakening of the Balmer lines typically seen at intermediate times when LRNe have their second peak, or the molecular lines associated with LRNe at late times. In certain respects, AT 2022fnm resembles AT 2018hso, though AT 2018hso is a LRN with a double-peaked light curve.

Unlike LRNe, which display a second peak in photospheric radius after the maximum, the photospheric radii of ILRTs typically undergo a monotonic decrease (Botticella et al. 2009; Cai et al. 2019, 2021). AT 2022fnm, on the other hand, displays aberrant behaviour in this regard, which precludes the use of this metric as a tidy method of separating LRNe and ILRTs. In Fig. 6 over the first 20 d we see the photosphere moving inwards and then outwards again. The decrease in radius could potentially be explained by increasing transparency due to hydrogen recombination, which is indeed seen at around 6000 K in SN IIP models, for example (Kasen & Woosley 2009).

After the decrease in radius, a source of ionising photons would then be required to push the photosphere outwards. The spectral evolution provides some clues: the continuum towards the blue end weakens between +13 d and +23 d, which is potentially indicative of increased opacity as the bluer wavelengths experience greater scattering. Furthermore, the beginning of the increase in radius at around +15 d also corresponds to the beginning of the plateau seen in the light curve. Taken together, these observations could naïvely be seen to suggest increased interaction, with the outer layers of ejecta hitting the CSM and leading to the creation of photons and hence ionisation, increasing the opacity. Alternatively, we may be seeing a scenario similar to the increase in opacity routinely observed in SNe IIP, where the shock ionises the outer envelope, leading to a decrease in transparency. The radius of AT 2022fnm shows an overall variation of approximately a factor of 1.6, whilst the low-luminosity classical IIP SNe SN 2005cs and SN 2009md display variation of approximately a factor of 2.5 (Jäger et al. 2020). However, in contrast to AT 2022fnm, the early evolution of SN 2005cs shows a monotonic increase in radius over time (unfortunately, no early radius measurements are available for SN 2009md; Jäger et al. 2020).

The increase in radius in AT 2022fnm from around +15 d is puzzling given that no temperature rise corresponding to the proposed ionisation is seen, though the decline rate of the luminosity does slow at this epoch. Additionally, the equivalent width of Hα (plotted in Fig. 12) is decreasing at this point. Similar, though less marked, behaviour was seen in the case of the ILRT SNhunt120 (Stritzinger et al. 2020a), whose photospheric radius also dipped early before climbing again.

thumbnail Fig. 12.

Evolution of the spectral flux density and the equivalent width of Hα. Top panel: evolution of the spectral flux density of Hα over time. Middle panel: evolution of the equivalent width of Hα over time. Bottom panel: evolution of the ratio of Hα to Hβ spectral flux density over time.

It is unclear whether the requirement of a second black-body component after around +100 d is due to the presence of multiple ejecta components or the presence of pre-existing dust. If there is pre-existing dust, it must not exist in the immediate vicinity of AT 2023fnm or it would be quickly overrun by the ejecta. The second black-body component has a temperature of 2200 ± 350 K, consistent with the sublimation temperature for amorphous carbon dust of approximately 1900 K (Landt et al. 2023). Unfortunately, the point at which a second black-body component becomes necessary is approximately coincident with the epoch of our final spectrum, and the lack of later spectra makes its difficult to make any further determination.

Spectroscopically, AT 2022fnm more closely resembles an ILRT than a LRN. Whilst the development of a red continuum is common to both LRNe and ILRTs, the presence of the [Ca II] λλ7291,7323 doublet is a key signature of ILRTs, though, as noted, it has also been seen in the LRN AT 2018hso. Ca H and K emission is visible in the +57 d and +90 d spectra, and in earlier spectra it is visible in absorption. However, these lines are also seen in emission in some LRNe and ILRTs (e.g. AT 2018hso and SN 2008S). The Ca IIλ8498,8542,8662 NIR triplet is also prominent in both ILRTs and LRNe, as seen in Fig. 11. Nevertheless, AT 2022fnm does not develop the forest of metal lines typically seen in absorption in LRNe, nor does it develop the strong molecular lines associated with LRNe at late times. This highlights the need for late-time spectra to disentangle these two classes of objects.

The absorption lines bluewards of Hα and Hβ (at 6485 Å and 4802 Å, respectively), which we identify as components of these lines, have a velocity of about −3500 km s−1. This is likely associated with the ejecta as the photosphere recedes in mass coordinates, revealing further material. Indeed, this recession in mass coordinates is apparent in Fig. 6 from around +30 d, when the photospheric radius stays roughly constant. There are also further absorption features redwards of Hα at 6574 Å and 6592 Å, with velocities of 500 km s−1 and 1300 km s−1, respectively.

AT 2022fnm can be said to be an ambiguous object in that it shows features consistent with LRNe and ILRTs, straddling the phase space covered by both categories of objects. AT 2022fnm is not alone in this respect, with certain other objects, such as AT 2018hso, sharing features of both the LRN and ILRT classifications (Cai et al. 2019). The simple analytical modelling we performed shows that the object is consistent with a weak eruption or explosion of a star of relatively low mass, though, as mentioned, the parameters of the model are highly degenerate. The presence of pre-existing dust, as evidenced by the NIR excess, is consistent with an ILRT; however, both pre-existing and newly formed dust can also be found in core-collapse SNe, such as SN 1987A (Roche et al. 1989; Lucy et al. 1989) and the type IIP SNe SN 2002hh (Meikle et al. 2006) and SN 2004et (Kotak et al. 2009), as well as in stellar-merger scenarios (Bond et al. 2003; Pastorello et al. 2019b), so it cannot be used to discriminate between these possibilities. Nevertheless, the preponderance of evidence supports the classification of AT 2022fnm as an ILRT, and the characteristics of the object are consistent with a weak stellar explosion or eruption. We expect the next generation of surveys to yield hundreds of LRNe and ILRTs, and suggest that a combination of relatively high-cadence optical and NIR photometry and spectroscopy will be needed if we wish to attribute various processes associated with the evolution of single and massive binary systems.


7

SNOoPY is a package for SN photometry using point-spread-function fitting and/or template subtraction developed by E. Cappellaro. A package description can be found at http://sngroup.oapd.inaf.it/ecsnoopy.html

8

IRAF is a software suite for the reduction and analysis of astronomical data. A description can be found at https://iraf-community.github.io/

Acknowledgments

We thank the anonymous referee for carefully reading the manuscript. S. Moran acknowledges support from the Magnus Ehrnrooth Foundation and the Vilho, Yrjö and Kalle Väisälä Foundation. RK acknowledges support from the Research Council of Finland (340613). MF is supported by a Royal Society – Science Foundation Ireland University Research Fellowship. AP, EC and NER and AR acknowledge support of the PRIN-INAF 2022 project “Shedding light on the nature of gap transients: from the observations to the models”. Y.-Z. Cai is supported by the National Natural Science Foundation of China (NSFC, Grant 12303054), the Yunnan Fundamental Research Projects (Grant 202401AU070063) and the International Centre of Supernovae, Yunnan Key Laboratory (Grant 202302AN360001). S. Mattila acknowledges support from the Academy of Finland project 350458. CPG acknowledges financial support from the Secretary of Universities and Research (Government of Catalonia) and by the Horizon 2020 Research and Innovation Programme of the European Union under the Marie Skłodowska-Curie and the Beatriu de Pinós 2021 BP 00168 programme, from the Spanish Ministerio de Ciencia e Innovación (MCIN) and the Agencia Estatal de Investigación (AEI) 10.13039/501100011033 under the PID2020-115253GA-I00 HOSTFLOWS project, and the programme Unidad de Excelencia María de Maeztu CEX2020-001058-M. NER acknowledges partial support from MIUR, PRIN 2017 (grant 20179ZF5KS) and from the Spanish MICINN grant PID2019-108709GB-I00 and FEDER funds. PL wishes to thank the Swedish Research Council for support. AVF’s group at UC Berkeley has received financial assistance from the Christopher R. Redlich Fund, Gary and Cynthia Bengier, Clark and Sharon Winslow, Alan Eustace (W.Z. is a Bengier-Winslow-Eustace Specialist in Astronomy), William Draper, Timothy and Melissa Draper, Briggs and Kathleen Wood, Sanford Robertson (T.G.B. is a Draper-Wood-Robertson Specialist in Astronomy), and numerous other donors; additional support was provide by NASA/HST grant AR-14295 from the Space Telescope Science Institute (STScI), which is operated by the Association of Universities for Research in Astronomy (AURA), Inc., under NASA contract NAS5-26555. This work is supported by the National Natural Science Foundation of China (NSFC grants 12033003 and 11633002), the Scholar Programme of Beijing Academy of Science and Technology (DZ:BS202002) and the Tencent Xplorer Prize. Based in part on observations made with the Nordic Optical Telescope, owned in collaboration by the University of Turku and Aarhus University, and operated jointly by Aarhus University, the University of Turku and the University of Oslo, representing Denmark, Finland and Norway, the University of Iceland and Stockholm University at the Observatorio del Roque de los Muchachos, La Palma, Spain, of the Instituto de Astrofisica de Canarias. The data presented here were obtained in part with ALFOSC, which is provided by the Instituto de Astrofisica de Andalucia (IAA) under a joint agreement with the University of Copenhagen and NOT. Based in part on observations made with the Gran Telescopio Canarias (GTC), installed in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias, in the island of La Palma. The W. M. Keck Observatory is operated as a scientific partnership among the California Institute of Technology, the University of California and NASA; the observatory was made possible by the generous financial support of the W. M. Keck Foundation. We thank WeiKang Zheng and Yi Yang for assistance with the Keck observations. This work has made use of data from the Asteroid Terrestrial-impact Last Alert System (ATLAS) project. ATLAS is primarily funded to search for near-Earth objects through NASA grants NN12AR55G, 80NSSC18K0284, and 80NSSC18K1575; byproducts of the NEO search include images and catalogues from the survey area. The ATLAS science products have been made possible through the contributions of the University of Hawaii Institute for Astronomy, the Queen’s University Belfast, the Space Telescope Science Institute, and the South African Astronomical Observatory. This work makes use of observations from the LCOGT network. We acknowledge ESA Gaia, DPAC, and the Photometric Science Alerts Team.

References

  1. Adams, S. M., Kochanek, C. S., Prieto, J. L., et al. 2016, MNRAS, 460, 1645 [NASA ADS] [CrossRef] [Google Scholar]
  2. Arnett, W. D. 1980, ApJ, 237, 541 [Google Scholar]
  3. Arnett, W. D. 1982, ApJ, 253, 785 [Google Scholar]
  4. Bellm, E. C., Kulkarni, S. R., Graham, M. J., et al. 2019, PASP, 131, 018002 [Google Scholar]
  5. Bennett, C. L., Larson, D., Weiland, J. L., & Hinshaw, G. 2014, ApJ, 794, 135 [Google Scholar]
  6. Berger, E., Soderberg, A. M., Chevalier, R. A., et al. 2009, ApJ, 699, 1850 [NASA ADS] [CrossRef] [Google Scholar]
  7. Blagorodnova, N., Kotak, R., Polshaw, J., et al. 2017, ApJ, 834, 107 [NASA ADS] [CrossRef] [Google Scholar]
  8. Blagorodnova, N., Klencki, J., Pejcha, O., et al. 2021, A&A, 653, A134 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Bond, H. E., Henden, A., Levay, Z. G., et al. 2003, Nature, 422, 405 [CrossRef] [Google Scholar]
  10. Bond, H. E., Bedin, L. R., Bonanos, A. Z., et al. 2009, ApJ, 695, L154 [Google Scholar]
  11. Bostroem, K. A., Valenti, S., Sand, D. J., et al. 2022, TNS Classif. Rep., 2022-857, 1 [Google Scholar]
  12. Botticella, M. T., Pastorello, A., Smartt, S. J., et al. 2009, MNRAS, 398, 1041 [NASA ADS] [CrossRef] [Google Scholar]
  13. Bowen, I. S. 1947, PASP, 59, 196 [NASA ADS] [CrossRef] [Google Scholar]
  14. Cai, Y. Z., Pastorello, A., Fraser, M., et al. 2018, MNRAS, 480, 3424 [NASA ADS] [Google Scholar]
  15. Cai, Y. Z., Pastorello, A., Fraser, M., et al. 2019, A&A, 632, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Cai, Y. Z., Pastorello, A., Fraser, M., et al. 2021, A&A, 654, A157 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Cai, Y. Z., Pastorello, A., Fraser, M., et al. 2022a, A&A, 667, A4 [Google Scholar]
  18. Cai, Y., Reguitti, A., Valerin, G., & Wang, X. 2022b, Universe, 8, 493 [NASA ADS] [CrossRef] [Google Scholar]
  19. Chatzopoulos, E., Wheeler, J. C., & Vinko, J. 2012, ApJ, 746, 121 [Google Scholar]
  20. Dessart, L., & Hillier, D. J. 2011, MNRAS, 410, 1739 [NASA ADS] [Google Scholar]
  21. Howell, D. 2019, Am. Astron. Soc. Meet. Abstr., 233, 258.16 [NASA ADS] [Google Scholar]
  22. Howell, D. A., & Global Supernova Project 2017, Am. Astron. Soc. Meet. Abstr., 230, 318.03 [NASA ADS] [Google Scholar]
  23. Humphreys, R. M., Bond, H. E., Bedin, L. R., et al. 2011, ApJ, 743, 118 [NASA ADS] [CrossRef] [Google Scholar]
  24. Jäger, Z., Jr., Vinkó, J., Bíró, B. I., et al. 2020, MNRAS, 496, 3725 [CrossRef] [Google Scholar]
  25. Kasen, D., & Woosley, S. E. 2009, ApJ, 703, 2205 [NASA ADS] [CrossRef] [Google Scholar]
  26. Kashi, A., Frankowski, A., & Soker, N. 2010, ApJ, 709, L11 [NASA ADS] [CrossRef] [Google Scholar]
  27. Kasliwal, M. M. 2012, PASA, 29, 482 [NASA ADS] [CrossRef] [Google Scholar]
  28. Kasliwal, M. M., Kulkarni, S. R., Arcavi, I., et al. 2011, ApJ, 730, 134 [NASA ADS] [CrossRef] [Google Scholar]
  29. Kotak, R., Meikle, W. P. S., Farrah, D., et al. 2009, ApJ, 704, 306 [NASA ADS] [CrossRef] [Google Scholar]
  30. Kozma, C., & Fransson, C. 1998, ApJ, 497, 431 [NASA ADS] [CrossRef] [Google Scholar]
  31. Landt, H., Mitchell, J. A. J., Ward, M. J., et al. 2023, ApJ, 945, 62 [NASA ADS] [CrossRef] [Google Scholar]
  32. Lucy, L. B., Danziger, I. J., Gouiffes, C., & Bouchet, P. 1989, in IAU Colloq. 120: Structure and Dynamics of the Interstellar Medium, eds. G. Tenorio-Tagle, M. Moles, & J. Melnick, 350, 164 [NASA ADS] [CrossRef] [Google Scholar]
  33. Lynch, D. K., Rudy, R. J., Russell, R. W., et al. 2004, ApJ, 607, 460 [CrossRef] [Google Scholar]
  34. Magnier, E. A., Schlafly, E. F., Finkbeiner, D. P., et al. 2020, ApJS, 251, 32 [Google Scholar]
  35. Martini, P., Wagner, R. M., Tomaney, A., et al. 1999, AJ, 118, 1034 [NASA ADS] [CrossRef] [Google Scholar]
  36. Masci, F. J., Laher, R. R., Rusholme, B., et al. 2019, PASP, 131, 018003 [Google Scholar]
  37. Matsumoto, T., & Metzger, B. D. 2022, ApJ, 938, 5 [NASA ADS] [CrossRef] [Google Scholar]
  38. Meikle, W. P. S., Mattila, S., Gerardy, C. L., et al. 2006, ApJ, 649, 332 [CrossRef] [Google Scholar]
  39. Metzger, B. D., & Pejcha, O. 2017, MNRAS, 471, 3200 [NASA ADS] [CrossRef] [Google Scholar]
  40. Nicholl, M. 2018, Res. Notes Am. Astron. Soc., 2, 230 [NASA ADS] [Google Scholar]
  41. Pastorello, A., & Fraser, M. 2019, Nat. Astron., 3, 676 [Google Scholar]
  42. Pastorello, A., Mason, E., Taubenberger, S., et al. 2019a, A&A, 630, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Pastorello, A., Chen, T. W., Cai, Y. Z., et al. 2019b, A&A, 625, L8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Pastorello, A., Fraser, M., Valerin, G., et al. 2021, A&A, 646, A119 [EDP Sciences] [Google Scholar]
  45. Pastorello, A., Valerin, G., Fraser, M., et al. 2023, A&A, 671, A158 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Perley, D. A. 2019, PASP, 131, 084503 [Google Scholar]
  47. Pumo, M. L., Turatto, M., Botticella, M. T., et al. 2009, ApJ, 705, L138 [NASA ADS] [CrossRef] [Google Scholar]
  48. Roche, P. F., Aitken, D. K., Smith, C. H., & James, S. D. 1989, Nature, 337, 533 [NASA ADS] [CrossRef] [Google Scholar]
  49. Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103 [Google Scholar]
  50. Schultz, G. V., & Wiemer, W. 1975, A&A, 43, 133 [NASA ADS] [Google Scholar]
  51. Smith, K. W., Smartt, S. J., Young, D. R., et al. 2020, PASP, 132, 085002 [Google Scholar]
  52. Stritzinger, M. D., Taddia, F., Fraser, M., et al. 2020a, A&A, 639, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Stritzinger, M. D., Taddia, F., Fraser, M., et al. 2020b, A&A, 639, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Thompson, T. A., Prieto, J. L., Stanek, K. Z., et al. 2009, ApJ, 705, 1364 [NASA ADS] [CrossRef] [Google Scholar]
  55. Tonry, J. L., Denneau, L., Heinze, A. N., et al. 2018, PASP, 130, 064505 [Google Scholar]
  56. Tonry, J., Denneau, L., Weiland, H., et al. 2022, TNSTR, 2022-836, 1 [NASA ADS] [Google Scholar]
  57. Tully, R. B., Courtois, H. M., Dolphin, A. E., et al. 2013, AJ, 146, 86 [NASA ADS] [CrossRef] [Google Scholar]
  58. Tylenda, R., Hajduk, M., Kamiński, T., et al. 2011, A&A, 528, A114 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

Appendix A: Tables

Table A.1.

Comparison objects.

Table A.2.

Optical photometry. B-band and V-band photometry are given in Vega magnitudes, the rest in AB magnitudes.

Table A.3.

NIR photometry given in Vega magnitudes.

Table A.4.

Optical spectroscopy. Instrumental resolutions were determined based on measurement of the [O I] λ5577.34 skyline except in the case of those marked with an asterisk, where reference values were used.

All Tables

Table 1.

Analytical diffusion model parameters.

Table A.1.

Comparison objects.

Table A.2.

Optical photometry. B-band and V-band photometry are given in Vega magnitudes, the rest in AB magnitudes.

Table A.3.

NIR photometry given in Vega magnitudes.

Table A.4.

Optical spectroscopy. Instrumental resolutions were determined based on measurement of the [O I] λ5577.34 skyline except in the case of those marked with an asterisk, where reference values were used.

All Figures

thumbnail Fig. 1.

B-band image taken at the NOT (ALFOSC) at +3 d showing the location of AT 2022fnm ( α = 12 h 25 m 33 . s 713 $ \alpha=12^{\mathrm{h}}25^{\mathrm{m}}33{{{{\overset{\text{s}}{.}}}}}713 $, δ = + 45 ° 41 12 . 45 $ \delta =+45\circ41{{\prime}}12{{{\overset{\prime\prime}{.}}}}45 $, J2000) in its host galaxy, NGC 4389.

In the text
thumbnail Fig. 2.

Multi-band light curve of AT 2022fnm. The inverted open triangles represent upper limits. The vertical line marks the epoch of the onset of the outburst. ATLAS data beyond +40 d are also shown with open symbols as they are noisier.

In the text
thumbnail Fig. 3.

r-band comparison of AT 2022fnm with other ILRTs and LRNe from the literature. All objects have been corrected for reddening with the values listed in Table A.1 and respective references. The reference epoch for each object has been set to its first peak.

In the text
thumbnail Fig. 4.

Colour evolution of AT 2022fnm and comparison objects. We calculated the colours with observations taken within 0.3 d of each other except in the case of V − Ks, for which we allowed a ±3.5 d window, owing to the modest cadence of the NIR observations.

In the text
thumbnail Fig. 5.

Rise time and decline rates plotted against the (first) r-band peak for AT 2022fnm and comparison objects. Rise times plotted with filled markers are calculated against explosion or outburst epoch, whilst those with unfilled markers are calculated against the detection epoch. Rise times are calculated to the peak observed magnitude, which could be more sensitive to photometric uncertainties for objects whose light curve flattens around the peak.

In the text
thumbnail Fig. 6.

Pseudo-bolometric properties. Top panel: pseudo-bolometric light curve of AT 2022fnm. Second panel: temperature versus time of a black body that has been fitted to the pseudo-bolometric light curve. Third panel: log of the radius versus time of a black body that has been fitted to the pseudo-bolometric light curve. Bottom panel: expansion velocity (radius over time) versus time of a black body that has been fitted to the pseudo-bolometric light curve. Note that the radius is simply being divided by the elapsed time at each point to generate this crude measure. The data have been corrected for redshift and foreground extinction. Points corresponding to the SED fits from Fig. 8 have been overplotted for comparison in the form of red triangles. The figure has been truncated at +80 d as the light curve is dominated by numerical artefacts beyond this epoch.

In the text
thumbnail Fig. 7.

Attempt to reproduce the pseudo-bolometric light curve with a simple analytical diffusion model. All three components of the model are shown along with their total, which fits the observed luminosity well. The model is somewhat degenerate, and the details of this are discussed in the text, as are the parameters. The best-fitting explosion epoch from our modelling is MJD 59 669.38, which is consistent with our adopted outburst epoch (MJD 59 667.98 ± 1.5).

In the text
thumbnail Fig. 8.

Sequence of SED fits for AT 2022fnm. At late times a second black-body component (red) is required, and the sum of the two components is shown in green.

In the text
thumbnail Fig. 9.

Optical spectral sequence of AT 2022fnm. The spectra have been calibrated against the photometry. The offsets applied to each spectrum are shown in the plot, and the +90 d spectrum has been multiplied by a factor of 2 for clarity. Negative fluxes have been set to zero. The +2 d spectrum is plotted in grey with a version binned to 10 Å overplotted in black to compensate for the noise of the spectrum. The absorption features seen at approximately 6725 Å and 8200 Å are artefacts.

In the text
thumbnail Fig. 10.

Evolution of Hα and Hβ in AT 2022fnm over time. All spectra have been binned to a common pixel size and normalised against the peak of Hα in the left panel and Hβ in the right. Only spectra with both Hα and Hβ coverage are shown. Vertical lines marking ±1000, 2000, and 3000 km s−1 have been added to guide the eye.

In the text
thumbnail Fig. 11.

Comparison of an early spectrum and late spectrum of AT 2022fnm with other objects. The spectra have been normalised against the sum of the flux densities except in the case of the insets, which have been normalised against the peak of Hα. The phases are calculated relative to the epoch of the onset of the outburst. The upper-left inset shows Hα. The +83 d AT 2018hso spectrum has been trimmed to 4500–9000 Å and binned to 10 Å, owing to noise.

In the text
thumbnail Fig. 12.

Evolution of the spectral flux density and the equivalent width of Hα. Top panel: evolution of the spectral flux density of Hα over time. Middle panel: evolution of the equivalent width of Hα over time. Bottom panel: evolution of the ratio of Hα to Hβ spectral flux density over time.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.