Light Curve Morphology Research Articles

Single-pulse Gamma-Ray Bursts Have Prevalent Hard-to-soft Spectral Evolution

We analyze the spectral evolution of 62 bright Fermi gamma-ray bursts with large enough signal-to-noise to allow for time-resolved spectral analysis. We develop a new algorithm to test for single-pulse morphology that is insensitive to the specific shape of pulses. Instead, it only checks whether or not there are multiple, isolated, or statistical significant peaks in the light curve. In addition, we carry out a citizen science test to assess light-curve morphology and spectral evolution. We find that, no matter the adopted assessment method, bursts characterized by single-peaked prompt emission light curves have a greater tendency to also have a consistently decaying peak energy or hard-to-soft spectral evolution. This contrasts with the behavior of multipeaked bursts, for which the tendency is to have a peak frequency that is not monotonically decreasing. We discuss this finding in the theoretical framework of internal/external shocks and find it to be consistent with at least some single-pulse bursts associated with particularly high-density environments.

Exocomet Models in Transit: Light Curve Morphology in the Optical—Near Infrared Wavelength Range

Following the widespread practice of exoplanetary transit simulations, various presumed components of an extrasolar system can be examined in numerically simulated transits, including exomoons, rings around planets, and the deformation of exoplanets. Template signals can then be used to efficiently search for light curve features that mark specific phenomena in the data, and they also provide a basis for feasibility studies of instruments and search programs. In this paper, we present a method for exocomet transit light curve calculations using arbitrary dust distributions in transit. The calculations, spanning four distinct materials (carbon, graphite, pyroxene, and olivine), and multiple dust grain sizes (100–300 nm, 300–1000 nm, and 1000–3000 nm) encompass light curves in VRJHKL bands. We also investigated the behavior of scattering colors. We show that multicolor photometric observations are highly effective tools in the detection and characterization of exocomet transits. They provide information on the dust distribution of the comet (encoded in the light curve shape), while the color information itself can reveal the particle size change and material composition of the transiting material, in relation to the surrounding environment. We also show that the typical cometary tail can result in the wavelength dependence of the transit timing. We demonstrate that multi-wavelength observations can yield compelling evidence for the presence of exocomets in real observations.

The JWST weather report from the nearest brown dwarfs I: multiperiod JWST NIRSpec + MIRI monitoring of the benchmark binary brown dwarf WISE 1049AB

ABSTRACT We report results from 8 h of JWST/MIRI low resolution spectroscopic (LRS) monitoring directly followed by 7 h of JWST/NIRSpec prism spectroscopic monitoring of the benchmark binary brown dwarf WISE 1049AB, the closest, brightest brown dwarfs known. We find water, methane, and CO absorption features in both components, including the 3.3 μm methane absorption feature and a tentative detection of small grain ($\lt $ 1μm) silicate absorption at $\gt $8.5 μm in WISE 1049A. Both components vary significantly ($\gt 1~{{\rm per\ cent}}$), with WISE 1049B displaying larger variations than WISE 1049A. Using K-means clustering, we find three main transition points in wavelength for both components of the binary: (1) change in behaviour at $\sim$2.3 μm coincident with a CO absorption bandhead, (2) change in behaviour at 4.2 μm, close to the CO fundamental band at $\lambda \gt $ 4.4 µm, and (3) change in behaviour at 8.3–8.5 µm, potentially corresponding to silicate absorption. We interpret the light curves observed with both NIRSpec and MIRI as likely stemming from (1) a deep pressure level driving the double-peaked variability seen in WISE 1049B at wavelengths $\lt $2.3 and $\gt $8.5 µm, (2) an intermediate pressure level shaping the light-curve morphology between 2.3 and 4.2 µm, and (3) a higher altitude pressure level producing single-peaked and plateaued light-curve behaviour between 4.2 and 8.5 µm.

The Hertzsprung progression of classical Cepheids in the Gaia era

ABSTRACT A new fine grid of non-linear convective pulsation models for the so-called ‘bump Cepheids’ is presented to investigate the Hertzprung progression (HP) phenomenon shown by their light and radial pulsation velocity curves. The period corresponding to the centre of the HP is investigated as a function of various model assumptions, such as the efficiency of superadiabatic convection, the mass–luminosity relation, and the metal and helium abundances. The assumed mass–luminosity relation is found to significantly affect the phenomenon but variations in the chemical composition as well as in the stellar mass (at fixed mass–luminosity relation) also play a key role in determining the value of the HP centre period. Finally, the predictive capability of the presented theoretical scenario is tested against observed light curves of bump Cepheids in the ESA Gaia data base, also considering the variation of the pulsation amplitudes and of the Fourier parameters R21 and Φ21 with the pulsation period. A qualitative agreement between theory and observations is found for what concerns the evolution of the light curve morphology as the period moves across the HP centre, as well for the pattern in period–amplitude, period–R21, and period–Φ21 planes. A larger sample of observed Cepheids with accurate light curves and metallicities is required in order to derive more quantitative conclusions.

A hidden population of massive white dwarfs: two spotted K+WD binaries

Abstract The identification and characterization of massive (≳ 0.8 M⊙) white dwarfs is challenging in part due to their low luminosity. Here we present two candidate single-lined spectroscopic binaries, Gaia DR3 4014708864481651840 and 5811237403155163520, with K-dwarf primaries and optically dark companions. Both have orbital periods of P ∼ 0.45 days and show rotational variability, ellipsoidal modulations, and high-amplitude radial velocity variations. Using light curves from the Transiting Exoplanet Survey Satellite (TESS), radial velocities from ground-based spectrographs, and spectral energy distributions, we characterize these binaries to describe the nature of the unseen companion. We find that both systems are consistent with a massive white dwarf companion. Unlike simple ellipsoidal variables, star spots cause the light curve morphology to change between TESS sectors. We attempt to constrain the orbital inclination using PHOEBE binary light curve models, but degeneracies in the light curves of spotted stars prevent a precise determination. Finally, we search for similar objects using Gaia DR3 and TESS, and comment on these systems in the context of recently claimed compact object binaries.

A Possible Mechanism for the “Late Phase” in Stellar White-light Flares

M dwarf flares observed by the Transiting Exoplanet Survey Satellite (TESS) sometimes exhibit a peak-bump light-curve morphology, characterized by a secondary, gradual peak well after the main, impulsive peak. A similar late phase is frequently detected in solar flares observed in the extreme ultraviolet from longer hot coronal loops distinct from the impulsive flare structures. White-light emission has also been observed in off-limb solar flare loops. Here, we perform a suite of one-dimensional hydrodynamic loop simulations for M dwarf flares inspired by these solar examples. Our results suggest that coronal plasma condensation following impulsive flare heating can yield high electron number density in the loop, allowing it to contribute significantly to the optical light curves via free-bound and free–free emission mechanisms. Our simulation results qualitatively agree with TESS observations: the longer evolutionary timescale of coronal loops produces a distinct, secondary emission peak; its intensity increases with the injected flare energy. We argue that coronal plasma condensation is a possible mechanism for the TESS late-phase flares.

Late-time HST UV and optical observations of AT 2018cow: extracting a cow from its background

ABSTRACT The bright, blue, rapidly evolving AT 2018cow is a well-studied peculiar extragalactic transient. Despite an abundance of multiwavelength data, there still is no consensus on the nature of the event. We present our analysis of three epochs of Hubble Space Telescope (HST) observations spanning the period from 713 to 1474 d post-burst, paying particular attention to uncertainties of the transient photometry introduced by the complex background in which AT 2018cow resides. Photometric measurements show evident fading in the UV and more subtle but significant fading in the optical. During the last HST observation, the transient’s optical/UV colours were still bluer than those of the substantial population of compact, young, star-forming regions in the host of AT 2018cow, suggesting some continued transient contribution to the light. However, a compact source underlying the transient would substantially modify the resulting spectral energy distribution, depending on its contribution in the various bands. In particular, in the optical filters, the complex, diffuse background poses a problem for precise photometry. An underlying cluster is expected for a supernova occurring within a young stellar environment or a tidal-disruption event (TDE) within a dense older one. While many recent works have focused on the supernova interpretation, we note the substantial similarity in UV light-curve morphology between AT 2018cow and several tidal disruption events around supermassive black holes. Assuming AT 2018cow arises from a TDE-like event, we fit the late-time emission with a disc model and find MBH = 103.2 ± 0.8 M⊙. Further observations are necessary to determine the late-time evolution of the transient and its immediate environment.

MOBSTER – VII. Using light curves to infer magnetic and rotational properties of stars with centrifugal magnetospheres

ABSTRACT Early-type B stars with strong magnetic fields and rapid rotation form centrifugal magnetospheres (CMs), as the relatively weak stellar wind becomes magnetically confined and centrifugally supported above the Kepler co-rotation radius. CM plasma is concentrated at and above the Kepler co-rotation radius at the intersection between the rotation and magnetic field axis. Stellar rotation can cause these clouds of material to intersect the viewer’s line of sight, leading to photometric eclipses. However, for stars with strong ($\sim 10\, {\rm kG}$) magnetic fields and rapid rotation, CMs can become optically thick enough for emission to occur via electron scattering. Using high-precision space photometry from a sample of stars with strong H α emission, we apply simulated light curves from the rigidly rotating magnetosphere model to directly infer magnetic and rotational properties of these stars. By comparing the values inferred from photometric modelling to those independently determined by spectropolarimetry, we find that magnetic obliquity angle β, viewer inclination i, and critical rotation fraction W can be approximately recovered for three of the four stars studied here. However, there are large discrepancies between the optical depth at the Kepler radius τK expected from magnetometry, and the values required to match the observations. We show that τK of order unity is needed to reasonably match the light-curve morphology of our sample stars.

Forward modelling of brightness variations in Sun-like stars

Context. The amplitude and morphology of light curves of Sun-like stars change substantially with increasing rotation rate: brightness variations are amplified and become more regular. This has not been explained so far. Aims. We develop a modelling approach for calculating brightness variations of stars with various rotation rates and use it to explain the observed trends in stellar photometric variability. Methods. We combined numerical simulations of magnetic flux emergence and transport with a model for stellar brightness variability to calculate synthetic light curves of stars as observed by the Kepler telescope. We computed the distribution of the magnetic flux on the stellar surface for various rotation rates and degrees of active-region nesting (i.e. the tendency of active regions to emerge in the vicinity of recently emerged regions). Using the resulting maps of the magnetic flux, we computed the rotational variability of our simulated stellar light curves as a function of rotation rate and nesting of magnetic features and compared our calculations to Kepler observations. Results. We show that both the rotation rate and the degree of nesting have a strong impact on the amplitude and morphology of stellar light curves. In order to explain the variability of most of the Kepler targets with known rotation rates, we need to increase the degree of nesting to values that are much higher than the values on the Sun. Conclusions. The suggested increase in nesting with the rotation rate can provide clues about the flux emergence process for high levels of stellar activity.

Examining the Rotation Period Distribution of the 40 Myr Tucana–Horologium Association with TESS

The Tucana–Horologium association (Tuc-Hor) is a 40 Myr old moving group in the southern sky. In this work, we measure the rotation periods of 313 Tuc-Hor objects with TESS light curves derived from TESS full-frame images and membership lists driven by Gaia EDR3 kinematics and known youth indicators. We recover a period for 81.4% of the sample and report 255 rotation periods for Tuc-Hor objects. From these objects we identify 11 candidate binaries based on multiple periodic signals or outlier Gaia DR2 and EDR3 renormalized unit weight error values. We also identify three new complex rotators (rapidly rotating M dwarf objects with intricate light-curve morphology) within our sample. Along with the six previously known complex rotators that belong to Tuc-Hor, we compare their light-curve morphology between TESS Cycle 1 and Cycle 3 and find that they change substantially. Furthermore, we provide context for the entire Tuc-Hor rotation sample by describing the rotation period distributions alongside other youth indicators such as Hα and Li equivalent width, as well as near-ultraviolet and X-ray flux. We find that measuring rotation periods with TESS is a fast and effective means to confirm members in young moving groups.

A survey for variable young stars with small telescopes: VI – Analysis of the outbursting Be stars NSW 284, gaia 19eyy, and VES 263

ABSTRACT This paper is one in a series reporting results from small telescope observations of variable young stars. Here, we study the repeating outbursts of three likely Be stars based on long-term optical, near-infrared, and mid-infrared photometry for all three objects, along with follow-up spectra for two of the three. The sources are characterized as rare, truly regularly outbursting Be stars. We interpret the photometric data within a framework for modelling light-curve morphology, and find that the models correctly predict the burst shapes, including their larger amplitudes and later peaks towards longer wavelengths. We are thus able to infer the start and end times of mass loading into the circumstellar discs of these stars. The disc sizes are typically 3 – 6 times the areas of the central star. The disc temperatures are ∼40 per cent, and the disc luminosities are ∼10 per cent of those of the central Be star, respectively. The available spectroscopy is consistent with inside-out evolution of the disc. Higher excitation lines have larger velocity widths in their double-horned shaped emission profiles. Our observations and analysis support the decretion disc model for outbursting Be stars.

What can Gaussian processes really tell us about supernova light curves? Consequences for Type II(b) morphologies and genealogies

ABSTRACT Machine learning has become widely used in astronomy. Gaussian process (GP) regression in particular has been employed a number of times to fit or resample supernova (SN) light curves, however by their nature typical GP models are not suited to fit SN photometric data and they will be prone to overfitting. Recently GP resampling was used in the context of studying the morphologies of Type II and IIb SNe and they were found to be clearly distinct with respect to four parameters: the rise time (trise), the magnitude difference between 40 and 30 d post-explosion (Δm40–30), the earliest maximum (post-peak) of the first derivative (dm1), and minimum of the second derivative (dm2). Here we take a close look at GP regression and its limitations in the context of SN light curves in general, and we also discuss the uncertainties on these specific parameters, finding that dm1 and dm2 cannot give reliable astrophysical information. We do reproduce the clustering in trise–Δm40–30 space, although it is not as clear cut as previously presented. The best strategy to accurately populate the trise–Δm40–30 space will be to use an expanded sample of high-quality light curves [such as those in the Asteroid Terrestrial-impact Last Alert System (ATLAS) transient survey] and analytical fitting methods. Finally, using the bpass fiducial models, we predict that future photometric studies will reveal clear clustering of the Type IIb and II light curve morphologies with a distinct continuum of transitional events.

Theoretical Stellar Pulsation Physics

AbstractPulsating stars play a crucial role in the calibration of the cosmic distance scale as well as in tracing the properties of the associated stellar populations. In the era of large observational surveys and precise astrometric missions, it is crucial to rely on accurate stellar pulsation models able to predict the observed behaviors for different physical assumptions. Indeed, the relations currently used in the literature to derive individual and mean distances of mainly radially pulsating stars such as Cepheids and RR Lyrae are well physically understood, but are also known to depend on a number of often unknown parameters. Recent extensive sets of stellar pulsation models developed by various authors show how variations in the physical assumptions can affect the theoretical prediction of the instability strip boundaries, the morphology and amplitude of light and radial velocity curves, and the consequent Period-Luminosity, Period-Luminosity-Color and Period-Wesenheit relations. These aspects are discussed in the framework of current open problems in the field of classical pulsating stars.

A Population of Dipper Stars from the Transiting Exoplanet Survey Satellite Mission

Dipper stars are a classification of young stellar objects that exhibit dimming variability in their light curves, dropping in brightness by 10%–50%, likely induced by occultations due to circumstellar disk material. This variability can be periodic, quasiperiodic, or aperiodic. Dipper stars have been discovered in young stellar associations via ground-based and space-based photometric surveys. We present the detection and characterization of the largest collection of dipper stars to date: 293 dipper stars, including 234 new dipper candidates. We have produced a catalog of these targets, which also includes young stellar variables that exhibit predominately burst-like variability and symmetric variability (equal parts bursting and dipping). The total number of catalog sources is 414. These variable sources were found in a visual survey of TESS light curves, where dip-like variability was observed. We found a typical age among our dipper sources of <5 Myr, with the age distribution peaking at ≈2 Myr, and a tail of the distribution extending to ages older than 20 Myr. Regardless of the age, our dipper candidates tend to exhibit infrared excess, which is indicative of the presence of disks. TESS is now observing the ecliptic plane, which is rich in young stellar associations, so we anticipate many more discoveries in the TESS data set. A larger sample of dipper stars would enhance the census statistics of light-curve morphologies and dipper ages.

Lightcurves and Rotations of Trans-Neptunian Objects in the 2:1 Mean Motion Resonance with Neptune

We report the rotational lightcurves of 21 trans-Neptunian objects (TNOs) in Neptune’s 2:1 mean motion resonance obtained with the 6.5 m Magellan-Baade telescope and the 4.3 m Lowell Discovery Telescope. The main survey’s goal is to find objects displaying a large lightcurve amplitude that is indicative of contact binaries or highly elongated objects. In our sample, two 2:1 resonant TNOs showed a significant short-term lightcurve amplitude: 2002 VD130 and (531074) 2012 DX98. The full lightcurve of 2012 DX98 infers a periodicity of 20.80 ± 0.06 hr and amplitude of 0.56 ± 0.03 mag, whereas 2002 VD130 rotates in 9.85 ± 0.07 hr with a 0.31 ± 0.04 mag lightcurve amplitude. Based on lightcurve morphology, we classify (531074) 2012 DX98 as a likely contact binary but 2002 VD130 as a likely single elongated object. Based on our sample and the lightcurves reported in the literature, we estimate the lower percentage of nearly equal-sized contact binaries at only 7%–14% in the 2:1 resonance, which is comparable to the low fraction reported for the dynamically cold classical TNOs. This low contact binary fraction in the 2:1 Neptune resonance is consistent with the lower estimate of the recent numerical modeling. We report the Sloan g′, r′, and i′ surface colors of 2002 VD130, which is an ultra-red TNO whereas 2012 DX98 is a very red object based on published surface colors.

Large-amplitude periodic outbursts and long-period variables in the VVV VIRAC2-β data base

ABSTRACT The VISTA Variables in the Via Lactea (VVV) survey obtained near-infrared photometry towards the Galactic bulge and the southern disc plane for a decade (2010–2019). We designed a modified Lomb–Scargle method to search for large-amplitude ($\Delta K_{s, 2-98{{\ \rm per\ cent}}}$ &amp;gt; 1.5 mag) mid to long-term periodic variables (P&amp;gt; 10 d) in the 2nd version of VVV Infrared Astrometric Catalogue (VIRAC2-β). In total, 1520 periodic sources were discovered, including 59 candidate periodic outbursting young stellar objects (YSOs), based on the unique morphology of the phase-folded light curves, proximity to Galactic H ii regions and mid-infrared colours. Five sources are spectroscopically confirmed as accreting YSOs. Both fast-rise/slow-decay and slow-rise/fast-decay periodic outbursts were found, but fast-rise/slow-decay outbursts predominate at the highest amplitudes. The multiwavelength colour variations are consistent with a variable mass accretion process, as opposed to variable extinction. The cycles are likely to be caused by dynamical perturbations from stellar or planetary companions within the circumstellar disc. An additional search for periodic variability amongst YSO candidates in published Spitzer-based catalogues yielded a further 71 candidate periodic accretors, mostly with lower amplitudes. These resemble cases of pulsed accretion but with unusually long periods and greater regularity. The majority of other long-period variables are pulsating dusty Miras with smooth and symmetric light curves. We find that some Miras have redder W3 − W4 colours than previously thought, most likely due to their surface chemical compositions.

New Time-resolved, Multi-band Flares in the GJ 65 System with gPhoton

Abstract Characterizing the distribution of flare properties and occurrence rates is important for understanding habitability of M-dwarf exoplanets. The Galaxy Evolution Explorer (GALEX) space telescope observed the GJ 65 system, composed of the active, flaring M stars BL Cet and UV Cet, for 15,900 s (∼4.4 hr) in two ultraviolet (UV) bands. The contrast in flux between flares and the photospheres of cool stars is maximized at UV wavelengths, and GJ 65 is the brightest and nearest flaring M-dwarf system with significant GALEX coverage. It therefore represents the best opportunity to measure low-energy flares with GALEX. We construct high-cadence lightcurves from calibrated photon events and find 13 new flare events with near-UV (NUV) energies ranging from 1028.5–1029.5 erg and recover one previously reported flare with an energy of 1031 erg. The newly reported flares are among the smallest M-dwarf flares observed in the UV with sufficient time resolution to discern lightcurve morphology. The estimated flare frequency at these low energies is consistent with extrapolation from the distributions of higher-energy flares on active M dwarfs measured by other surveys. The largest flare in our sample is bright enough to exceed the local nonlinearity threshold of the GALEX detectors, which precludes color analysis. However, we detect quasi-periodic pulsations during this flare in both the far-UV and NUV bands at a period of ∼50 s, which we interpret as a modulation of the flare’s chromospheric thermal emission through periodic triggering of reconnection by external MHD oscillations in the corona.

The puzzling story of flare inactive ultra fast rotating M dwarfs – I. Exploring their magnetic fields

ABSTRACT Stars which are rapidly rotating are expected to show high levels of activity according to the activity–rotation relation. However, previous TESS studies have found ultra fast rotating (UFR) M dwarfs with periods less than 1 d displaying low levels of flaring activity. As a result, in this study, we utilize VLT/FORS2 spectro-polarimetric data of 10 M dwarf UFR stars between spectral types ∼M2–M6 all with Prot &amp;lt; 1, to detect the presence of a magnetic field. We divide our sample into rotation period bins of equal size, with one star having many more flares in the TESS light curve than the other. We also provide an analysis of the long-term variability within our sample using TESS light curves taken during Cycles 1 and 3 (up to 3 yr apart). We identify 605 flares from our sample which have energies between 2.0 × 1031 and 5.4 × 1034 erg. Although we find no significant difference in the flare rate between the Cycles, two of our targets display changes in their light-curve morphology, potentially caused by a difference in the spot distribution. Overall, we find five stars (50 per cent) in our sample have a detectable magnetic field with strengths ∼1–2 kG. Of these five, four were the more flare active stars within the period bins with one being the less flare active star. It would appear the magnetic field strength may not be the answer to the lack of flaring activity and supersaturation or magnetic field configuration may play a role. However, it is clear the relationship between rotation and activity is more complex than a steady decrease over time.

GRB 191016A: a highly collimated gamma-ray burst jet with magnetized energy injection

ABSTRACT Long gamma-ray burst GRB 191016A was a bright and slow rising burst that was detected by the Swift satellite and followed up by ground based Liverpool Telescope (LT). LT follow up started 2411 s after the Swift Burst Alert Telescope (BAT) trigger using imager IO:O around the time of the late optical peak. From 3987–7687 s, we used the LT polarimeter RINGO3 to make polarimetric and photometric observations of the GRB simultaneously in the V, R, and I bands. The combined optical light curve shows an initial late peak followed by a decline until 6147 s, 6087 s, and 5247 s for I, R, and V filters respectively followed by a flattening phase. There is evidence of polarization at all phases including polarization ($P = 14.6 \pm 7.2 {{\ \rm per\ cent}}$) which is coincident with the start of the flattening phase. The combination of the light curve morphology and polarization measurement favours an energy injection scenario where slower magnetized ejecta from the central engine catches up with the decelerating blast wave. We calculate the minimum energy injection to be ΔE/E &amp;gt; 0.36. At a later time, combining the optical light curve from Burst Observer and Optical Transient Exploring System (BOOTES) (reported via GCN) and IO:O we see evidence of a jet break with jet opening angle 2°.

How Temporal Symmetry Defines Morphology in BATSE Gamma-Ray Burst Pulse Light Curves

Abstract We present compelling evidence that most gamma-ray burst (GRB) pulse light curves can be characterized by a smooth single-peaked component coupled with a more complex emission structure that is temporally symmetric around the time of the pulse peak. The model successfully fits 86% of Burst and Transient Source Experiment GRB pulses bright enough for structural properties to be measured. Surprisingly, a GRB pulse’s light-curve morphology can be accurately predicted by the pulse asymmetry and the stretching/compression needed to align the structural components preceding the temporal mirror with the time-reversed components following it. Such a prediction is only possible because GRB pulses exhibit temporal symmetry. Time-asymmetric pulses include fast rise exponential decays, rollercoaster pulses, and asymmetric u-pulses, while time-symmetric pulses include u-pulses and crowns. Each morphological type is characterized by specific asymmetries, stretching parameters, durations, and alignments between the smooth and structured components, and a delineation in the asymmetry/stretching distribution suggests that symmetric pulses and asymmetric pulses may belong to separate populations. Furthermore, pulses belonging to the short GRB class exhibit similar morphologies to the long GRB class, but appear to simply occur on shorter timescales.