Search

In spite of their key role in galaxy evolution and several decades of observational efforts, the census of supernova remnants (SNRs) in our Galaxy remains incomplete. Theoretical predictions based on the local supernova rate estimate the expected number of SNRs in the Galaxy to be $\gtrsim$ 1000. By contrast, the number of detected SNRs amounts to about 300. High-resolution, wide-area radio surveys at low frequencies are ideal tools with which to find missing SNRs, given the prominence of these sources at low radio frequencies. We aim to find missing SNRs using proprietary data from the LOFAR Two-Metre Sky Survey (LoTSS) at 144~MHz. We used LoTSS total intensity maps of two Galactic regions, one with $39^\mathrm{o} < l < 66^\mathrm{o}$ and $|b|< 2.5^\mathrm{o}$, and the other with $145^\mathrm{o} < l < 150^\mathrm{o}$ and $|b| < 3^\mathrm{o}$, in addition to mid-infrared (MIR) data from the Wide-field Infrared Survey Explorer (WISE) all-sky survey to search for SNR candidates. We report the discovery of 14 new SNR candidates selected on the basis of their morphology at 144 MHz and a lack of MIR emission. We also follow up on 24 previously reported SNR candidates, inferring their spectral index between the LoTSS frequency (144 MHz) and the frequency at which they were reported. The high resolution and sensitivity of LoTSS observations has resulted in the detection of 14 new SNR candidates. In order to unambiguously confirm the SNR nature of these candidates, follow-up X-ray observations are required with facilities such as eROSITA., Comment: Accepted in A&A

A new package for nonlinear least squares fitting is introduced in this paper. This package implements a recently developed algorithm that, for certain types of nonlinear curve fitting, reduces the number of nonlinear parameters to be fitted. One notable feature of this method is the absence of initialization which is typically necessary for nonlinear fitting gradient-based algorithms. Instead, just some bounds for the nonlinear parameters are required. Even though convergence for this method is guaranteed for exponential decay using the max-norm, the algorithm exhibits remarkable robustness, and its use has been extended to a wide range of functions using the Euclidean norm. Furthermore, this data-fitting package can also serve as a valuable resource for providing accurate initial parameters to other algorithms that rely on them.

We present the first deep (72 hours of observations) radio image of the Euclid Deep Field North (EDFN) obtained with the LOw-Frequency ARray (LOFAR) High Band Antenna (HBA) at 144 MHz. The EDFN is the latest addition to the LOFAR Two-Metre Sky Survey (LoTSS) Deep Fields and these observations represent the first data release for this field. The observations produced a 6" resolution image with a central r.m.s. noise of $32\,\mu$Jy\,beam$^{-1}$. A catalogue of $\sim 23,000$ radio sources above a signal-to-noise ratio (SNR) threshold of 5 is extracted from the inner circular 10 deg$^2$ region. We discuss the data analysis and we provide a detailed description of how we derived the catalogue of radio sources and on the issues related to direction-dependent calibration and their effects on the final products. Finally, we derive the radio source counts at 144 MHz in the EDFN using catalogues of mock radio sources to derive the completeness correction factors. The source counts in the EDFN are consistent with those obtained from the first data release of the other LoTSS Deep Fields (ELAIS-N1, Lockman Hole and Bootes), despite the different method adopted to construct the final catalogue and to assess its completeness., Comment: accepted for publication in A&A

A multivalued projection is an idempotent linear relation with invariant domain. We characterize multivalued projections that are operator ranges (called semiclosed) and provide several formulae of them. Moreover, we study the decomposability and continuity of multivalued projections, and describe nilpotent relations.

Context. Stars in transition phases, like those showing the B[e] phenomenon and luminous blue variables (LBVs), undergo strong, often irregular mass ejection events. The prediction of these phases in stellar evolution models is therefore extremely difficult if not impossible. As a result, their effective temperatures, their luminosities and even their true nature are not fully known. Aims. A suitable procedure to derive the stellar parameters of these types of objects is to use the BCD spectrophotometric classification system, based on the analysis of the Balmer discontinuity. The BCD parameters ({\lambda}_1, D) are independent of interstellar extinction and circumstellar contributions. Methods. We obtained low-resolution spectra for 14 stars with the B[e] phenomenon and LBVs. Using the BCD method, we derived the stellar and physical parameters. The study was complemented with the information provided by the JHK colour-colour diagram. Results. For each star, the BCD system gives a complete set of fundamental parameters and related quantities such as luminosity and distance. We confirmed HK Ori, HD 323771 and HD 52721 as pre-main sequence HAe/B[e], AS 202 and HD 85567 as FS CMa-type, and HD 62623 as sgB[e] stars. We classified Hen 3-847, CD-24 5721, and HD 53367 as young B[e] stars or FS CMa-type candidates, and HD 58647 as a slightly evolved B[e] star. In addition, Hen 3-1398 is an sgB[e] and MWC 877, CPD-59 2854 and LHA 120-S 65 are LBV candidates. The stellar parameters of the latter two LBVs are determined for the first time. Conclusions. Our results emphasise that the BCD system is a highly valuable tool to derive stellar parameters and physical properties of B-type stars in transition phases. This method can be combined with near-IR colour-colour diagrams to determine or confirm the evolutionary stage of emission-line stars with dust disks.

We report the discovery of a ring of low surface brightness radio emission around the Calvera pulsar, a high Galactic latitude, isolated neutron star, in the LOFAR Two-metre Sky Survey (LoTSS). It is centered at $\alpha=14\mathrm{h}11\mathrm{m}12.6\mathrm{s}$, $\delta=+79^\mathrm{o}23'15"$, has inner and outer radii of $14.2'$ and $28.4'$, and an integrated flux density at 144 MHz of $1.08\pm0.15$ Jy. The ring center is offset by $4.9'$ from the location of the Calvera pulsar. H$\alpha$ observations with the Isaac Newton Telescope show no coincident optical emission, but do show a small ($\sim20"$) optical structure internal to the ring. We consider three possible interpretations for the ring: that it is an H~II region, a supernova remnant (SNR), or an Odd Radio Circle (ORC). The positional coincidence of the ring, the pulsar, and an X-ray-emitting non-equilibrium ionisation plasma previously detected, lead us to prefer the SNR interpretation. If the source is indeed a SNR and its association with the Calvera pulsar is confirmed, then Calvera's SNR, or G118.4+37.0, will be one of few SNRs in the Galactic halo., Comment: Accepted in A&A

It is already accepted that Be stars are surrounded by circumstellar envelopes, which are mostly compatible with a disc geometry in Keplerian rotation. We aim to obtain a more complete characterisation of the properties of the circumstellar environment of Be stars that helps to constrain the theoretical models of the Be phenomenon. We present near-infrared, medium-resolution spectra of a sample of galactic Be stars with different spectral subtypes and luminosity classes. We measure different parameters of the hydrogen recombination lines from the Paschen, Brackett, Pfund, and Humphreys series, and use them to diagnose physical conditions in the circumstellar environment. We analysed the equivalent-width (EW) ratio between Br$\alpha$ and Br$\gamma$ lines and different diagrams of flux ratios. We also identify lines from He I, C I, N I, O I, Na I, Mg I, Mg II, Si I, Fe I, and Fe II. Analysing the EW measurements of particular He I, Mg II, Fe I, Fe II and O I lines, we find that for some lines they correlate with the spectral type of the star. Particularly, the emission of the O I~$\lambda\,1.3168\,\mu$m line decreases towards the later spectral types. We present an atlas of 22 Be stars, that covers a wide infrared (IR) spectral range with quasi-simultaneous observations. From a detailed analysis, we define new complementary criteria to Mennickent's classification of Be stars according to their disc opacity. Some objects in our sample present compact thick envelopes, while in others the envelope is extended and optically thin. The correlation between the full widths at half maximum (FWHM) and the peak separation ($\Delta \mathrm{V}$) versus $V\,\sin\,i$ for the Br10, Br$\delta$, and Hu14 lines reveals that the broadening mechanism is rotational. The Ly$\beta$ fluorescence is a key mechanism to explain the intensity of the emission of Mg II and O I lines., Comment: 27 pages, 14 figures. Accepted for publication in A&A

Core-collapse supernovae (CCSNe) produce large ($\gtrsim 0.1 \, {\rm M}_\odot$) masses of dust, and are potentially the primary source of dust in the Universe, but much of this dust may be destroyed before reaching the interstellar medium. Cassiopeia A (Cas A) is the only supernova remnant where an observational measurement of the dust destruction efficiency in the reverse shock is possible at present. We determine the pre- and post-shock dust masses in Cas A using a substantially improved dust emission model. In our preferred models, the unshocked ejecta contains $0.6-0.8 \, {\rm M}_\odot$ of $0.1 \, {\rm \mu m}$ silicate grains, while the post-shock ejecta has $0.02-0.09 \, {\rm M}_\odot$ of $5-10 {\, {\rm nm}}$ grains in dense clumps, and $2 \times 10^{-3} \, {\rm M}_\odot$ of $0.1 \, {\rm \mu m}$ grains in the diffuse X-ray emitting shocked ejecta. The implied dust destruction efficiency is $74-94 \%$ in the clumps and $92-98 \%$ overall, giving Cas A a final dust yield of $0.05-0.30 \, {\rm M}_\odot$. If the unshocked ejecta grains are larger than $0.1 \, {\rm \mu m}$, the dust masses are higher, the destruction efficiencies are lower, and the final yield may exceed $0.5 \, {\rm M}_\odot$. As Cas A has a dense circumstellar environment and thus a much stronger reverse shock than is typical, the average dust destruction efficiency across all CCSNe is likely to be lower, and the average dust yield higher. This supports a mostly-stellar origin for the cosmic dust budget., Comment: 10 pages, 9 figures, MNRAS accepted

Characterizations of the star, minus and diamond orders of operators are given in various contexts and the relationship between these orders is made more transparent. Moreover, we introduce a new partial order of operators which provides a unified scenario for studying the other three orders.

In the group of B stars with spectroscopic peculiarities, we can find the Be and the B[e] stars. The Be stars are early-type rapid rotators that present, as their principal characteristic, emission lines of hydrogen and singly ionized metals due to the presence of a gaseous envelope. The B[e] stars present in their spectra heterogeneous features that reveal the presence of regions with very different properties in a gaseous and dusty envelope. Our goal is to study the evolution of the disks around peculiar B stars through the variability of their physical properties and dynamical structure, as well as to set constraints on different models and disk forming mechanisms. Throughout the last decade, we have carried out temporal monitoring of a sample of objects in the near infrared using spectroscopic facilities at the Gemini and Las Campanas Observatories. In the present work, we focus on the classical Be star 12 Vul, for which also optical spectra have been collected quasi-simultaneously. We observed variability in the hydrogen line profiles of 12 Vul, attributed to dissipating and building-up processes of the circumstellar envelope. Also, we found that this Be star presented the $^{12}$CO band heads in emission in one observation. The emission of this molecule has not been previously reported in a Be star, while it is a common feature among B[e] stars. We obtained parameters to describe the $^{12}$CO emitting region and propose different scenarios to explain this intriguing emission., Comment: 6 pages, 3 figures. Accepted for publication in A&A

In this article we apply proper splittings of matrices to develop an iterative process to approximate solutions of matrix equations of the form TX = W. Moreover, by using the partial order induced by positive semidefinite matrices, we obtain equivalent conditions to the convergence of this process. We also include some speed comparison results of the convergence of this method. In addition, for all matrix T we propose a proper splitting based on the polar decomposition of T.

A significant number of Be stars show a second Balmer discontinuity (sBD) attributed to an extended circumstellar envelope (CE). The fast rotational velocity of Be stars undoubtedly plays a significant role in the formation of the CE. However, Bn stars, which are also B-type rapidly rotating stars, do not all present clear evidence of being surrounded by circumstellar material. We aim to characterize the populations of Be and Bn stars, and discuss the appearance of the sBD as a function of the stellar parameters. We expect to find new indices characterizing the properties of CEs in Be stars and properties relating Be and Bn stars. Correlations of the aspect and intensity of the sBD and the emission in the H$\alpha$ line with the stellar parameters and the $V\!\sin i$ are presented. Some Bn stars exhibit the sBD in absorption, which may indicate the presence of rather dense CEs. Six Bn stars show emission in the H$\alpha$ line, so they are reclassified as Be stars. The sBD in emission appears in Be stars with $V\!\sin i \lesssim 250$ km\,s$^{-1}$, and in absorption in both Be and Bn stars with \mbox{$V\!\sin i \gtrsim 50$ km\,s$^{-1}$}. Low-mass Be and Bn stars share the same region in the Hertzsprung-Russell diagram. The distributions of rotational to critical velocity ratios of Be and Bn stars corresponding to the current stellar evolutionary stage are similar, while distributions inferred for the zero-age main sequence have different skewness. We found emission in the H$\alpha$ line and signs of a CE in some Bn stars, which motivated us to think that Bn and Be stars probably belong to the same population. It should be noted that some of the most massive Bn stars could display the Be phenomenon at any time. The similarities found among Be and Bn stars deserve to be more deeply pursued., Comment: 30 pages, 9 tables, 19 figures. Accepted for publication in A&A

Ion acceleration in the MeV range can be routinely achieved with table-top laser technology. One of the current challenges is to improve the energy coupling from the laser to the proton beam without increasing the laser peak power. Introducing nanostructures at the front target surface was shown to be beneficial for an efficient transfer of energy to the electrons. In this manuscript, we study by using full-scale three-dimensional particle-in-cell simulations and finite laser pulses, the process when a proposed optimal target with triangular nanostructure (previously found to allow 97% laser energy absorption) is used . We demonstrate that the absorbed laser energy does not depend on the dimensionality in the range of parameters presented. We also present an analytical model for laser absorption that includes deviations from the ideal conditions. This is supported by a numerical parameter study that establishes the tolerance with respect to the nanostructure size, use of different ion species, existence of preplasma, etc. We found that altering the target thickness or using different ions does not affect the absorption, but it does affect the energy redistribution among the different plasma species. The optimal configuration ($h = 1~\lambda,~ w = 0.7~ \lambda$) is robust with respect to the target fabrication errors. However, high contrast laser pulses are required, because a pre-plasma layer with a thickness on the order of 0.5 lambda is enough to lower the laser absorption by more than a 10% in a non-optimal scenario.

Mixed-morphology supernova remnants (MM SNRs) are a mysterious class of objects that display thermal X-ray emission within their radio shell. They are an older class of SNRs, and as such are profoundly affected by the environment into which they evolve. VRO 42.05.01 is a MM SNR of puzzling morphology in the direction of the Galactic anticentre. Low-frequency radio observations of supernova remnants are sensitive to synchrotron electrons accelerated in the shock front. We aim to compare the low-frequency emission to higher frequency observations to understand the environmental and shock acceleration conditions that have given rise to the observed properties of this source. We present a LOFAR High Band Antenna map centred at 143 MHz of the region of the Galactic plane centred at $l = 166^\rm{o}, b = 3.5^\rm{o}$ at 143 MHz, with a resolution of 148'' and an rms noise of 4.4 mJy bm$^{-1}$ . Our map is sensitive to scales as large as $6^\rm{o}$. We compared the LOw Frequency ARay (LOFAR) observations to archival higher frequency radio, infrared, and optical data to study the emission properties of the source in different spectral regimes. We did this both for the SNR and for OA 184, an H II region within our field of view. We find that the radio spectral index of VRO 42.05.01 increases at low radio frequencies; i.e. the LOFAR flux is higher than expected from the measured spectral index value at higher radio frequencies. This observed curvature in the low-frequency end of the radio spectrum occurs primarily in the brightest regions of the source, while the fainter regions present a roughly constant power-law behaviour between 143 MHz and 2695 MHz. We favour an explanation for this steepening whereby radiative shocks have high compression ratios and electrons of different energies probe different length scales across the shocks, therefore sampling regions of different compression ratios., Comment: Accepted in A&A

The advanced stages of several high-mass stars are characterized by episodic mass loss shed during phases of instability. Key for assigning these stars a proper evolutionary state is to assess the composition and geometry of their ejecta alongside the stellar properties. We selected five hot LBV candidates in M33 to refine their classification, investigate their circumstellar environments and explore their evolutionary properties. Being accessible targets in the near-infrared, we conducted medium-resolution spectroscopy with GNIRS/GEMINI in the $K-$band to investigate their molecular circumstellar environments. Two stars were found to display CO emission, which was modeled to emerge from a circumstellar or circumbinary Keplerian disk/ring. The identification of the carbon isotope $^{13}$C and, for one of the two stars, a significantly low $^{12}$CO/$^{13}$CO ratio, implies an evolved stellar state. As both CO emission stars are highly luminous and hence do not undergo a red supergiant phase, we suggest that stripping processes and equatorial high-density ejecta due to fast rotation are responsible for the enrichment of the stellar surface with processed material from the core. A candidate B[e]SG displays an absorption CO profile, which may be attributed to a jet or stellar pulsations. The featureless infrared spectra of two stars suggest a low-density gas shell or dissipation of the molecule due to the ionizing temperature of the star. We propose spectroscopic monitoring of our targets to evaluate the stability of the CO molecule and assess the time-dependent dynamics of the circumstellar gas structures., Comment: 13 pages, 8 figures, 4 tables. MNRAS in press

Cassiopeia A is one of the best-studied supernova remnants. Its shocked ejecta emits brightly in radio and X-rays. Its unshocked ejecta can be studied through infrared emission, the radio-active decay of $^{44}$Ti, and low frequency free-free absorption due to cold gas internal to the shell. Free-free absorption is affected by the mass, geometry, temperature, and ionisation conditions in the absorbing gas. Observations at the lowest radio frequencies constrain a combination of these properties. We use LOFAR LBA observations at 30-77 MHz and L-band VLA observations to compare $u-v$-matched images with a common resolution of 17". We simultaneously fit, per pixel, for the emission measure and the ratio of the emission from the unabsorbed front of the shell versus the absorbed back of the shell. We explore the effects that low temperatures and a high degree of clumping can have on the derived physical properties, such as mass and density. We also compile published radio flux measurements, fit for the absorption processes that occur in the radio band, and consider how they affect the secular decline of the source. We find a mass in the unshocked ejecta of $M = 2.95 \pm {0.48} \,M_{\odot}$ for an assumed gas temperature of $T=100$ K. This estimate is reduced for colder gas temperatures and if the ejecta are clumped. We measure the reverse shock to have a radius of $114$" $\pm $6". We also find that a decrease in the amount of mass in the unshocked ejecta (as more and more material meets the reverse shock and heats up) cannot account for the observed low frequency behaviour of the secular decline rate. To reconcile our low frequency absorption measurements with models that predict little mass in the unshocked ejecta we need the ejecta to be very clumped, or the temperature in the cold gas to be low ($\sim10$ K). Both conditions can jointly contribute to the high absorption., Comment: Accepted for publication in A&A v2: including the DOI, language edits

The interaction of high-intensity laser pulses and solid targets provides a promising way to create compact, tunable and bright XUV attosecond sources that can become a unique tool for a variety of applications. However, it is important to control the polarization state of this XUV radiation, and to do so in the most efficient regime of generation. Using the relativistic electronic spring (RES) model and particle-in-cell (PIC) simulations, we show that the polarization state of the generated attosecond pulses can be tuned in a wide range of parameters by adjusting the polarization and angle of incidence of the laser radiation. In particular, we demonstrate the possibility of producing circularly polarized attosecond pulses in a wide variety of setups., Comment: 6 pages, 3 figures

We assume that stars may undergo surface differential rotation to study its impact on the interpretation of $V\!\sin i$ and on the observed distribution $\Phi(u)$ of ratios of true rotational velocities $u=V/V_\rm c$ ($V_\rm c$ is the equatorial critical velocity). We discuss some phenomena affecting the formation of spectral lines and their broadening, which can obliterate the information carried by $V\!\sin i$ concerning the actual stellar rotation. We studied the line broadening produced by several differential rotational laws, but adopted Maunder's expression $\Omega(\theta)=\Omega_o(1+\alpha\cos^2\theta)$ as an attempt to account for all of these laws with the lowest possible number of free parameters. We studied the effect of the differential rotation parameter $\alpha$ on the measured $V\!\sin i$ parameter and on the distribution $\Phi(u)$ of ratios $u=V/V_\rm c$. We conclude that the inferred $V\!\sin i$ is smaller than implied by the actual equatorial linear rotation velocity $V_\rm eq$ if the stars rotate with $\alpha<0$, but is larger if the stars have $\alpha>0$. For a given $|\alpha|$ the deviations of $V\!\sin i$ are larger when $\alpha<0$. If the studied Be stars have on average $\alpha<0$, the number of rotators with $V_\rm eq\simeq0.9V_\rm c$ is larger than expected from the observed distribution $\Phi(u)$; if these stars have on average $\alpha>0$, this number is lower than expected. We discuss seven phenomena that contribute either to narrow or broaden spectral lines, which blur the information on the rotation carried by $V\!\sin i$ and, in particular, to decide whether the Be phenomenon mostly rely on the critical rotation. We show that two-dimensional radiation transfer calculations are needed in rapid rotators to diagnose the stellar rotation more reliably., Comment: To appear in A&A