Flare Stars Research Articles

Unexpected stratification in the equivalent-duration distributions of flare stars

ABSTRACT In this study, the results obtained from the statistical models of flares observed with ground-based telescopes and TESS observations for the UV Ceti-type stars AD Leo and V1005 Ori are presented. The distribution of flare equivalent duration versus total flare duration for the flares of both stars, observed with different telescopes, has been modelled using the one-phase exponential association (OPEA) model. Although no bifurcation was observed in the models derived with the flare data of V1005 Ori, the OPEA model exhibited bifurcation in almost all observation data of AD Leo. While many similar studies have shown that such models derived for a star do not vary over time and exhibit a single distribution in the same observation season, the stratification in AD Leo’s distributions and consequently in its model is an unexpected result.

Separating states in astronomical sources using hidden Markov models: with a case study of flaring and quiescence on EV Lac

ABSTRACT We present a new method to distinguish between different states (e.g. high and low, quiescent and flaring) in astronomical sources with count data. The method models the underlying physical process as latent variables following a continuous-space Markov chain that determines the expected Poisson counts in observed light curves in multiple passbands. For the underlying state process, we consider several autoregressive processes, yielding continuous-space hidden Markov models of varying complexity. Under these models, we can infer the state that the object is in at any given time. The continuous state predictions from these models are then dichotomized with the help of a finite mixture model to produce state classifications. We apply these techniques to X-ray data from the active dMe flare star EV Lac, splitting the data into quiescent and flaring states. We find that a first-order vector autoregressive process efficiently separates flaring from quiescence: flaring occurs over 30 per cent–40 per cent of the observation durations, a well-defined persistent quiescent state can be identified, and the flaring state is characterized by higher plasma temperatures and emission measures.

Flaring Stars in a Nontargeted Millimeter-wave Survey with SPT-3G

We present a flare star catalog from 4 yr of nontargeted millimeter-wave survey data from the South Pole Telescope (SPT). The data were taken with the SPT-3G camera and cover a 1500 deg2 region of the sky from 20h40m0s to 3h20m0s in right ascension and from −42° to −70° in declination. This region was observed on a nearly daily cadence from 2019 to 2022 and chosen to avoid the plane of the galaxy. A short-duration transient search of this survey yields 111 flaring events from 66 stars, increasing the number of both flaring events and detected flare stars by an order of magnitude from the previous SPT-3G data release. We provide cross-matching to Gaia DR3, as well as matches to X-ray point sources found in the second ROSAT all-sky survey. We have detected flaring stars across the main sequence, from early-type A stars to M dwarfs, as well as a large population of evolved stars. These stars are mostly nearby, spanning 10–1000 pc in distance. Most of the flare spectral indices are constant or gently rising as a function of frequency at 95/150/220 GHz. The timescale of these events can range from minutes to hours, and the peak ν L ν luminosities range from 1027 to 1031 erg s−1 in the SPT-3G frequency bands.

The White-light Superflares from Cool Stars in GWAC Triggers

M-type stars are the ones that flare most frequently, but how big their maximum flare energy can reach is still unknown. We present 163 flares from 162 individual M2 through L1-type stars that triggered the GWAC, with flare energies ranging from 1032.2 to 1036.4 erg. The flare amplitudes range from △G = 0.84 to ∼10 mag. Flare energy increases with stellar surface temperature (T eff) but both △G and equivalent duration log10(ED) seem to be independent of T eff. Combining periods detected from light curves of TESS and K2, spectra from LAMOST, the Sloan Digital Sky Survey, the 2.16 m telescope, and the Gaia DR3 data, we found that these GWAC flare stars are young. For the stars that have spectra, we found that these stars are in or very near the saturation region, and log10(LHα/Lbol) is lower for M7–L1 stars than for M2–M6 stars. We also studied the relation between GWAC flare bolometric energy E bol and stellar hemispherical area S, and found that log10Ebol (in erg) increases with increasing S (in square centimeters), and the maximum flare energy log10Ebol,max≥log10S+14.25 . For M7–L1 stars, there seem to be other factors limiting their maximum flare energies in addition to the stellar hemispherical area.

Strengthening nuclear symmetry energy constraints using multiple resonant shattering flares of neutron stars with realistic mass uncertainties

ABSTRACT With current and planned gravitational-wave (GW) observing runs, coincident multimessenger timing of resonant shattering flares (RSFs) and GWs may soon allow for neutron star (NS) asteroseismology to be used to constrain the nuclear symmetry energy, an important property of fundamental nuclear physics that influences the composition and equation of state of NSs. In this work, we examine the effects of combining multiple RSF detections on these symmetry energy constraints, and consider how realistic uncertainties in the masses of the progenitor NSs may weaken them. We show that the detection of subsequent multimessenger events has the potential to substantially improve constraints beyond those obtained from the first, and that this improvement is insensitive to the mass of the NSs that produce the RSFs and its uncertainty. This sets these asteroseismic constraints apart from bulk NS properties such as radius, for which the NS mass is highly important, meaning that any multimessenger RSF and GW events can equally improve our knowledge of fundamental physics.

Detailed cool star flare morphology with CHEOPS and TESS

Context. White-light stellar flares are proxies for some of the most energetic types of flares, but their triggering mechanism is still poorly understood. As they are associated with strong X and ultraviolet emission, their study is particularly relevant to estimate the amount of high-energy irradiation onto the atmospheres of exoplanets, especially those in their stars’ habitable zone. Aims. We used the high-cadence, high-photometric capabilities of the CHEOPS and TESS space telescopes to study the detailed morphology of white-light flares occurring in a sample of 130 late-K and M stars, and compared our findings with results obtained at a lower cadence. Methods. We employed dedicated software for the reduction of 3 s cadence CHEOPS data, and adopted the 20 s cadence TESS data reduced by their official processing pipeline. We developed an algorithm to separate multi-peak flare profiles into their components, in order to contrast them to those of single-peak, classical flares. We also exploited this tool to estimate amplitudes and periodicities in a small sample of quasi-periodic pulsation (QPP) candidates. Results. Complex flares represent a significant percentage (≳30%) of the detected outburst events. Our findings suggest that high-impulse flares are more frequent than suspected from lower-cadence data, so that the most impactful flux levels that hit close-in exoplanets might be more time-limited than expected. We found significant differences in the duration distributions of single and complex flare components, but not in their peak luminosity. A statistical analysis of the flare parameter distributions provides marginal support for their description with a log-normal instead of a power-law function, leaving the door open to several flare formation scenarios. We tentatively confirmed previous results about QPPs in high-cadence photometry, report the possible detection of a pre-flare dip, and did not find hints of photometric variability due to an undetected flare background. Conclusions. The high-cadence study of stellar hosts might be crucial to evaluate the impact of their flares on close-in exoplanets, as their impulsive phase emission might otherwise be incorrectly estimated. Future telescopes such as PLATO and Ariel, thanks to their high-cadence capability, will help in this respect. As the details of flare profiles and of the shape of their parameter distributions are made more accessible by continuing to increase the instrument precision and time resolution, the models used to interpret them and their role in star-planet interactions might need to be updated constantly.

Searching for gamma-ray emission from stellar flares

Abstract Flares from magnetically active dwarf stars should produce relativistic particles capable of creating γ-rays. So far, the only isolated main sequence star besides the Sun to have been detected in γ-rays is TVLM 513-46546. Detecting γ-ray flares from more dwarf stars can improve our understanding of their magnetospheric properties, and could also indicate a diminished likelihood of their planets’ habitability. In this work, we stack data from the Fermi Gamma-ray Space Telescope during a large number of events identified from optical and X-ray flare surveys. We report an upper limit of γ-ray emission from the population of flare stars. Stacking results towards control positions are consistent with a non-detection. We compare these results to observed Solar γ-ray flares and against a model of emission from neutral pion decay. The upper limit is consistent with Solar flares when scaled to the flare energies and distances of the target stars. As with Solar flares, the neutral pion decay mechanism for γ-ray production is also consistent with these results.

Localizing Stellar Activity on Low-mass Stars with the Transiting Exoplanet Survey Satellite

Magnetic reconnection is the underlying cause of stellar flares which are linked to regions of high magnetic activity, like star spots. To understand trends in stellar activity, we need to study how it is heightened by rapid rotation and deep convection in young low-mass stars. We analyze light curves of such stars observed by the Transiting Exoplanet Survey Satellite to investigate trends in flare timing with starspot modulation signals. Trends in flare properties and spot modulation can provide a means to “localize” which face of a star flares more frequently and better understand their association with active regions. We present an analysis of light curves from M and K dwarfs with no companions from five nearby and young moving groups spanning ages ∼20–150 Myr. We discuss a technique to analyze the distribution of flares and star spots and describe our results, which reveal a tentative correlation.

Apache Point Observatory (APO)/SMARTS Flare Star Campaign Observations. I. Blue Wing Asymmetries in Chromospheric Lines during Mid-M-Dwarf Flares from Simultaneous Spectroscopic and Photometric Observation Data

We conducted the time-resolved simultaneous optical spectroscopic and photometric observations of mid-M-dwarf flare stars YZ CMi, EV Lac, and AD Leo. Spectroscopic observations were obtained using Apache Point Observatory 3.5 m and Small and Moderate Aperture Research Telescope System 1.5 m telescopes during 31 nights. Among the 41 detected flares, seven flares showed clear blue wing asymmetries in the Hα line, with various correspondences in flare properties. The duration of the blue wing asymmetries range from 20 minutes to 2.5 hr, including a flare showing the shift from blue to red wing asymmetry. Blue wing asymmetries can be observed during both white-light and candidate non-white-light flares. All of the seven flares showed blue wing asymmetries also in the Hβ line, but there are large varieties on which other chromospheric lines showed blue wing asymmetries. One among the 7 flares was also observed with soft X-ray spectroscopy, which enabled us to estimate the flare magnetic field and length of the flare loop. The line-of-sight velocities of the blueshifted components range from –73 to –122 km s−1. Assuming that the blueshifts were caused by prominence eruptions, the mass of upward-moving plasma was estimated to be 1015–1019 g, which are roughly on the relation between flare energy and erupting mass expected from solar coronal mass ejections (CMEs). Although further investigations are necessary for understanding the observed various properties, these possible prominence eruptions on M-dwarfs could evolve into CMEs, assuming the similar acceleration mechanism with solar eruptions.

XMM-Newton-discovered Fast X-ray Transients: host galaxies and limits on contemporaneous detections of optical counterparts

ABSTRACT Extragalactic fast X-ray transients (FXTs) are a class of soft (0.3–10 keV) X-ray transients lasting a few hundred seconds to several hours. Several progenitor mechanisms have been suggested to produce FXTs, including supernova shock breakouts, binary neutron star mergers, or tidal disruptions involving an intermediate-mass black hole and a white dwarf. We present detailed host studies, including spectroscopic observations of the host galaxies of seven XMM-Newton-discovered FXTs. The candidate hosts lie at redshifts 0.0928 <z < 0.645 implying peak X-ray luminosities of 1043 erg s−1<LX < 1045 erg s−1 and physical offsets of 1 kpc < rproj < 22 kpc. These observations increase the number of FXTs with a spectroscopic redshift measurement by a factor of 2, although we note that one event is re-identified as a Galactic flare star. We infer host star formation rates and stellar masses by fitting the combined spectroscopic and archival photometric data. We also report on a contemporaneous optical counterpart search to the FXTs in Pan-STARRS and ATLAS by performing forced photometry at the position of the FXTs. We do not find any counterpart in our search. Given our constraints, including peak X-ray luminosities, optical limits, and host properties, we find that XRT 110 621 is consistent with an supernova shock breakout (SN SBO) event. Spectroscopic redshifts of likely host galaxies for four events imply peak X-ray luminosities that are too high to be consistent with SN SBOs, but we are unable to discard either the binary neutron star or white dwarf–intermediate-mass black hole tidal disruption event scenarios for these FXTs.

Periodicity of Flare Frequency of Active Flare Stars Detected by Tess

Periodicity of Flare Frequency of Active Flare Stars Detected by Tess

Solar–Stellar Connection: X-Ray Flares to Energetic (E > 10 MeV) Particle Events

Energetic particle environments are an important factor for the viability of life on exoplanets surrounding flare stars. In the heliosphere, large gradual solar energetic (E > 10 MeV) particle (SEP) events are produced by shocks from fast coronal mass ejections (CMEs). Extensive observations of solar X-ray flares, CMEs, and SEP events can provide guidance for flare star models of stellar energetic particle (StEP) events, for which stellar flares, but only rarely the associated CMEs, are observed. Comparing an extensive list of peak fluxes, timescales, and peak temperatures of 585 ≥ M3.0 solar X-ray flares with the occurrence of associated SEP events of peak flux Ip > 1.4 proton flux units, enhanced with proxy decametric–hectometric type II radio bursts, we determine guidelines for StEP event outcomes, given only stellar X-ray flare inputs. Longer timescales and lower peak temperatures of X-ray flares with a given peak X-ray flux Fp are more favorable for occurrence of associated SEP/StEP events, which, however, are only a minority of all solar flare outcomes. Most solar flares do not result in SEP events, invalidating scaling laws between stellar flares, CMEs, and StEP events. We discuss recent observations and models of the flare–CME relationship and suggest that StEP intensities Ip may often be overestimated.

The statistical OPEA models indicate an anomaly in the solar X-ray flare energy levels

We have recognized an anomaly on the solar activity cycle, while we researching on the question that there can be two different flare mechanisms working on the stellar surface. We discuss an anomaly in the solar flare activity levels in the cycle minima and maxima from the 20th to 22nd Solar Activity Cycles, depending on the models and analyses of the X-ray data. In this study, we used the GOES satellite data accumulated from 1976 to 1989 and detected 670 solar flares. We have separately carried out two independent statistical analyses for the examination of the solar flare activity. Firstly, we derived the equivalent durations distribution via the total durations of flares. Contrary to our expectations based on flare stars, it seems that this distribution is modelled by four different fits, instead of just one model fit. We found that the model parameters, such as plateau, half−time, and span, vary via time in the different trends from the Solar Activity Cycle. Secondly, we have derived the cumulative frequency models for solar flares. Contrary to the studies in the literature, we found that these frequency models do not take shape depending on the minimum or maximum of a solar activity cycle. If these results are not an unexpected anomaly, this study indicates that the existence of a different trend can be revealed apart from the solar activity cycle for the flare energy level variation, when the solar flare behaviour is separately determined for each cycle and they are compared with each other.

Commensal transient searches in eight short gamma-ray burst fields

ABSTRACT A new generation of radio telescopes with excellent sensitivity, instantaneous uv coverage, and large fields of view, are providing unprecedented opportunities for performing commensal transient searches. Here, we present such a commensal search in deep observations of short gamma-ray burst fields carried out with the MeerKAT radio telescope in South Africa at 1.3 GHz. These four hour observations of eight different fields span survey lengths of weeks to months. We also carry out transient searches in time slices of the full observations, at time-scales of 15 min and 8 s. We find 122 variable sources on the long time-scales, of which 52 are likely active galactic nuclei, but there are likely also some radio flaring stars. While the variability is intrinsic in at least two cases, most of it is consistent with interstellar scintillation. In this study, we also place constraints on transient rates based on state-of-the-art transient simulations codes. We place an upper limit of 2 × 10−4 transients per day per square degree for transients with peak flux of 5 mJy, and an upper limit of 2.5 × 10−2 transients per day per square degree for transients with a fluence of 10 Jy ms, the minimum detectable fluence of our survey.

Fine Structures of Radio Bursts from Flare Star AD Leo with FAST Observations

Radio bursts from nearby active M-dwarfs have been frequently reported and extensively studied in solar or planetary paradigms. Whereas, their substructures or fine structures remain rarely explored despite their potential significance in diagnosing the plasma and magnetic field properties of the star. Such studies in the past have been limited by the sensitivity of radio telescopes. Here we report the inspiring results from the high time-resolution observations of a known flare star AD Leo with the Five-hundred-meter Aperture Spherical radio Telescope. We detected many radio bursts in the 2 days of observations with fine structures in the form of numerous millisecond-scale sub-bursts. Sub-bursts on the first day display stripe-like shapes with nearly uniform frequency drift rates, which are possibly stellar analogs to Jovian S-bursts. Sub-bursts on the second day, however, reveal a different blob-like shape with random occurrence patterns and are akin to solar radio spikes. The new observational results suggest that the intense emission from AD Leo is driven by electron cyclotron maser instability, which may be related to stellar flares or interactions with a planetary companion.

Three-dimensional SPH Simulations of FU Orionis Star Flares in the Clumpy Accretion Model

One of the early hypotheses about the origin of FU Orionis star (FUORs) outbursts explains them based on the fall of gas clumps from the remnants of protostellar clouds onto protoplanetary disks surrounding young stars. To calculate the consequences of such an event, we make 3D hydrodynamic simulations using the smoothed particle hydrodynamics method. It is shown that the fall of the clump on the disk in the vicinity of the star actually causes a burst of the star’s accretion activity, resembling in its characteristics the flares of known FUORs. In the region of incidence, an inhomogeneous gas ring is formed, which is inclined relative to the outer disk. During several revolutions around the star, this ring combines with the inner disk and forms a tilted disk. In the process of evolution, the inner disk expands, and its inclination relative to the outer disk decreases. After 100 revolutions, the angle of inclination is a few degrees. This result is of interest in connection with the discovery in recent years of protoplanetary disks, the inner region of which is inclined relative to the outer one. Such structures are usually associated with the existence in the vicinity of a star of a massive body (planet or brown dwarf), whose orbit is inclined relative to the plane of the disk. The results of our modeling indicate the possibility of an alternative explanation for this phenomenon.

Minute-cadence observations of the LAMOST fields with the TMTS – III. Statistical study of the flare stars from the first two years

ABSTRACTTsinghua University-Ma Huateng Telescopes for Survey (TMTS) aims to detect fast-evolving transients in the Universe, which has led to the discovery of thousands of short-period variables and eclipsing binaries since 2020. In this paper, we present the observed properties of 125 flare stars identified by TMTS within the first two years, with an attempt to constrain their eruption physics. As expected, most of these flares were recorded in late-type red stars with GBP − GRP >2.0 mag; however, the flares associated with bluer stars tend to be on average more energetic and have broader profiles. The peak flux (Fpeak) of the flare is found to depend strongly on the equivalent duration (ED) of the energy release, i.e. Fpeak∝ED0.72 ± 0.04, which is consistent with results derived from the Kepler and Evryscope samples. This relation is likely to be related to the magnetic loop emission, while, for the more popular non-thermal electron heating model, a specific time evolution may be required to generate this relation. We notice that flares produced by hotter stars have a flatter Fpeak - ED relation compared to that from cooler stars. This is related to the statistical discrepancy in light-curve shape of flare events with different colours. In spectra from LAMOST, we find that flare stars have apparently stronger H α emission than inactive stars, especially at the low-temperature end, suggesting that chromospheric activity plays an important role in producing flares. On the other hand, the subclass with frequent flares is found to show H α emission of similar strength in its spectra to that recorded with only a single flare but similar effective temperature, implying that chromospheric activity may not be the only trigger for eruptions.

Starspot Mapping with Adaptive Parallel Tempering. II. Application to TESS Data for M-dwarf Flare Stars AU Microscopii, YZ Canis Minoris, and EV Lacertae

Starspots and stellar flares are indicators of stellar magnetic activity. The magnetic energy stored around spots is thought to be the origin of flares, but the connection is not completely understood. To investigate the relation between spot locations deduced from light curves and the occurrence of flares therein, we perform starspot modeling for the TESS light curves of three M-dwarf flare stars, AU Mic, YZ CMi, and EV Lac, using the code implemented in Paper I. The code enables us to deduce multiple stellar/spot parameters by the adaptive parallel tempering algorithm efficiently. We find that flare occurrence frequency is not necessarily correlated with the rotation phases of the light curve for each star. The result of starspot modeling shows that any spot is always visible to the line of sight in all phases, and we suggest that this can be one of the reasons why there is no or low correlation between rotation phases and flare frequency. In addition, the amplitude and shape of the light curve for AU Mic and YZ CMi have varied in two years between different TESS cycles. The result of starspot modeling suggests that this can be explained by the variations of spot size and latitude.

Constraints on Stellar Flare Energy Ratios in the NUV and Optical from a Multiwavelength Study of GALEX and Kepler Flare Stars

We present a multiwavelength study of stellar flares on primarily G-type stars using overlapping time domain surveys in the near-ultraviolet (NUV) and optical regimes. The NUV (the Galaxy Evolution Explorer; GALEX) and optical (Kepler) wavelength domains are important for understanding energy fractionations in stellar flares, and for constraining the associated incident radiation on a planetary atmosphere. We follow up on the NUV flare detections presented in Brasseur et al., using coincident Kepler long (1557 flares) and short (two flares) cadence light curves. We find no evidence of optical flares at these times, and place limits on the flare energy ratio between the two wave bands. We find that the energy ratio is correlated with GALEX band energy, and extends over a range of about 3 orders of magnitude in the ratio of the upper limit of Kepler band flare energy to NUV flare energy at the same time for each flare. The two flares with Kepler short-cadence data indicate that the true Kepler band energy may be much lower than the long-cadence-based upper limit. A similar trend appears for the bulk flare energy properties of nonsimultaneously observed flares on the same stars. We provide updated models to describe the flare spectral energy distribution from the NUV through the optical including continua and emission lines to improve upon blackbody-only models. The spread of observed energy ratios is much larger than encompassed by these models and suggests new physics is at work. These results call for a better understanding of NUV flare physics and provide a cautionary tale about using only optical flare measurements to infer the UV irradiation of close-in planets.

A search for transient, monochromatic light from the Galactic plane

ABSTRACT The Galactic plane was searched for transient, monochromatic light at optical and near-IR wavelengths to detect pulses shorter than 1 s. An objective-prism Schmidt telescope of 0.28-m aperture and a CMOS camera were used to observe 973 square deg, with 8864 exposures of 1-s each, within a strip 2.1 deg wide along the Galactic plane, from Galactic longitude −4 deg to +248 deg. All exposures were analysed for transient, monochromatic sources using a ‘difference image’ algorithm that yielded 11 candidate sources. All 11 sources were found to be associated with either astrophysical emission-line objects or aircraft with subsecond blinking lights. Our survey ‘rediscovered’ many Wolf–Rayet stars, M dwarf flare stars, and planetary nebulae. It also identified an aircraft, of unknown type, that apparently had a nearly monochromatic lamp and a xenon lamp. This survey would have revealed optical and near-IR pulses having a power of ∼180 GW (wavelength dependent) if emitted by a 10-m aperture laser located 1 kiloparsec away. These non-detections of laser pulses from the Galactic plane, including a 10-deg region towards the Galactic centre, add to the non-detections from more than 5000 nearby stars. Indeed, all-sky surveys for emission-line objects (e.g. ionized gas, supernovae remnants, and active galactic nuclei) would have revealed lasers of a wide range of average brightness, wavelength, and cadence. The absence of beacons reveals more of a SETI desert, notably at the intensely surveyed optical and radio wavelengths.