Spectral Energy Distribution Modeling Research Articles

Multiwavelength study of different flaring and low-activity states of blazar 4C+21.35

ABSTRACT Blazars, a class of active galactic nuclei, emit over the entire accessible electromagnetic spectrum and modelling of their broad-band spectral energy distribution (SED) is the key to constrain the underlying emission mechanisms. Here we report the results on the one-zone leptonic emission modelling carried out on the blazar 4C+21.35 using multiwavelength data spanning over the period 2008–2018. Broad-band SED modelling using γ-ray data from Fermi-Large Area Telescope, X-ray data from Swift-XRT and AstroSat, and UV–optical data from Swift-UVOT, AstroSat, and Catalina Real-Time Transient Survey was carried out at seven different epochs, including three γ-ray flaring episodes and four quiescent periods (three long-term averaged ones and one during AstroSat observing period). Our SED modelling suggests that two compact emission regions originating at a different time outside the broad-line region and moving away from the core with variation primarily in the jet electron spectra can explain the emission from the high-, moderate-, and low-activity periods. The emissions from high- and first low-activity states are likely to have originated in the first region. The moderate- and second low-activity states are likely due to the second emission region with fresh particle acceleration/injection at a later time.

Tracking the year-to-year variation in the spectral energy distribution of the narrow-line Seyfert 1 galaxy Mrk 335

ABSTRACT Multiwavelength monitoring of Mrk 335 with Swift between 2007 and 2019 are used to construct annual spectral energy distributions (SEDs) and track year-to-year changes. Non-contemporaneous archival data prior to 2007 are used to build a bright state SED. In this work, the changes are examined and quantified to build the foundation for future SED modelling. The yearly SEDs trace a downward trend on the average, with the X-ray portion varying significantly and acquiring further lower values in the past two years when compared to the optical/UV portion of SED. The bolometric Eddington ratios derived using optical/UV to X-ray SEDs and the calculated X-ray luminosities show a gradual decrease over the monitoring period. Changes in the parameters over time are examined. Principal component analysis suggests that the primary variability is in the X-ray properties of Mrk 335. When looking at the broader picture of Mrk 335 and its behaviour, the X-rays, accounting most of the variability in the 13-yr data, are possibly driven by physical processes related to the corona or absorption whereas the modest optical–UV variations suggest their origin within the accretion disc. These results are consistent with the previous interpretation of Mrk 335 using the timing analyses on the monitoring data and spectral modelling of deep observations.

The AGN contribution to the UV–FIR luminosities of interacting galaxies and its role in identifying the main sequence

ABSTRACT Emission from active galactic nuclei (AGNs) is known to play an important role in the evolution of many galaxies including luminous and ultraluminous systems (U/LIRGs), as well as merging systems. However, the extent, duration, and exact effects of its influence are still imperfectly understood. To assess the impact of AGNs on interacting systems, we present a spectral energy distribution (SED) analysis of a sample of 189 nearby galaxies. We gather and systematically re-reduce archival broad-band imaging mosaics from the ultraviolet to the far-infrared using data from GALEX, SDSS, 2MASS, IRAS, WISE, Spitzer, and Herschel. We use spectroscopy from Spitzer/IRS to obtain fluxes from fine-structure lines that trace star formation and AGN activity. Utilizing the SED modelling and fitting tool cigale, we derive the physical conditions of the interstellar medium, both in star-forming regions and in nuclear regions dominated by the AGN in these galaxies. We investigate how the star formation rates (SFRs) and the fractional AGN contributions (fAGN) depend on stellar mass, galaxy type, and merger stage. We find that luminous galaxies more massive than about $10^{10} \,\rm {M}_{*}$ are likely to deviate significantly from the conventional galaxy main-sequence relation. Interestingly, infrared AGN luminosity and stellar mass in this set of objects are much tighter than SFR and stellar mass. We find that buried AGNs may occupy a locus between bright starbursts and pure AGNs in the fAGN–[Ne v]/[Ne ii] plane. We identify a modest correlation between fAGN and mergers in their later stages.

Correlation between optical and γ-ray flux variations in bright flat spectrum radio quasars

ABSTRACT Blazars are known to show flux variations over a range of energies from low-energy radio to high-energy γ-rays. Cross-correlation analysis of the optical and γ-ray light curves in blazars shows that flux variations are generally correlated in both bands, however, there are exceptions. We explored this optical–GeV connection in four flat spectrum radio quasars by a systematic investigation of their long-term optical and γ-ray light curves. On analysis of the four sources, namely 3C 273, 3C 279, PKS 1510−089, and CTA 102, we noticed different behaviours between the optical and GeV flux variations. We found instances when (i) the optical and GeV flux variations are closely correlated, (ii) there are optical flares without γ-ray counterparts, and (iii) γ-ray flares without optical counterparts. To understand these diverse behaviours, we carried out broad-band spectral energy distribution (SED) modelling of the sources at different epochs using a one-zone leptonic emission model. The optical–UV emission is found to be dominated by emission from the accretion disc in the sources PKS 1510−089, CTA 102, and 3C 273, while in 3C 279, the synchrotron radiation from the jet dominates the optical–UV emission. Our SED analysis indicates that (i) correlated optical and γ-ray flux variations are caused by changes in the bulk Lorentz factor (Γ), (ii) γ-ray flares without optical counterparts are due to increase in Γ and/or the electron energy density, and (iii) an optical flare without γ-ray counterpart is due to increase in the magnetic field strength.

Resolving Decades of Periodic Spirals from the Wolf–Rayet Dust Factory WR 112

Abstract WR 112 is a dust-forming carbon-rich Wolf–Rayet (WC) binary with a dusty circumstellar nebula that exhibits a complex asymmetric morphology, which traces the orbital motion and dust formation in the colliding winds of the central binary. Unraveling the complicated circumstellar dust emission around WR 112 therefore provides an opportunity to understand the dust formation process in colliding-wind WC binaries. In this work, we present a multi-epoch analysis of the circumstellar dust around WR 112 using seven high spatial resolution (FWHM ∼ 0.″3–0.″4) N-band (λ ∼ 12 μm) imaging observations spanning almost 20 yr and that includes images obtained from Subaru/COMICS in 2019 October. In contrast to previous interpretations of a face-on spiral morphology, we observe clear evidence of proper motion of the circumstellar dust around WR 112 consistent with a nearly edge-on spiral with a θ s = 55° half-opening angle and a ∼20 yr period. The revised near edge-on geometry of WR 112 reconciles previous observations of highly variable nonthermal radio emission that was inconsistent with a face-on geometry. We estimate a revised distance to WR 112 of kpc based on the observed dust expansion rate and a spectroscopically derived WC terminal wind velocity of km s−1. With the newly derived WR 112 parameters, we fit optically thin dust spectral energy distribution models and determine a dust production rate of M ⊙ yr−1, which demonstrates that WR 112 is one of the most prolific dust-making WC systems known.

Tracing the Coevolution Path of Supermassive Black Holes and Spheroids with AKARI-selected Ultraluminous IR Galaxies at Intermediate Redshifts

Abstract We present the stellar population and ionized-gas outflow properties of ultraluminous IR galaxies (ULIRGs) at z = 0.1–1.0 that are selected from the AKARI far-IR all-sky survey. We construct a catalog of 1077 ULIRGs to examine feedback effects after major mergers. Of the 1077 ULIRGs, 202 are spectroscopically identified by SDSS and Subaru/FOCAS observations. Thanks to the deeper depth and higher resolution of AKARI compared to the previous Infrared Astronomical Satellite (IRAS) survey and reliable identification from the Wide-field Infrared Survey Explorer (WISE) mid-IR pointing, the sample is unique in identifying optically faint (i ∼ 20) IR-bright galaxies, which could be missed in previous surveys. A self-consistent spectrum and broadband spectral energy distribution (SED) decomposition method, which constrains stellar population properties in SED modeling based on spectral fitting results, has been employed for 149 ULIRGs whose optical continua are dominated by host galaxies. They are massive galaxies ( –1012 M ) associated with intense star formation activities (SFR ∼ 200–2000 M yr−1). The sample covers a range of active galactic nucleus (AGN) bolometric luminosity of 1010–1013 L , and the outflow velocity measured from the [O iii] 5007 Å line shows a correlation with AGN luminosity. Eight galaxies show extremely fast outflows with velocity up to 1500–2000 km s−1. However, the coexistence of vigorous starbursts and strong outflows suggests the star formation has not been quenched during the ULIRG phase. By deriving the stellar mass and mass fraction of the young stellar population, we find no significant discrepancies between stellar properties of ULIRGs with weak and powerful AGNs. The results are not consistent with the merger-induced evolutionary scenario, which predicts that star formation–dominated ULIRGs will show smaller stellar masses and younger stellar populations compared to AGN-dominated ULIRGs.

A disc reflection model for ultra-soft narrow-line Seyfert 1 galaxies

ABSTRACT We present a detailed analysis of the XMM–Newton observations of five narrow-line Seyfert 1 galaxies (NLS1s). They all show very soft continuum emission in the X-ray band with a photon index of Γ ≳ 2.5. Therefore, they are referred to as ‘ultra-soft’ NLS1s in this paper. By modelling their optical/UV–X-ray spectral energy distribution (SED) with a reflection-based model, we find indications that the disc surface in these ultra-soft NLS1s is in a higher ionization state than other typical Seyfert 1 AGN. Our best-fitting SED models suggest that these five ultra-soft NLS1s have an Eddington ratio of λEdd = 1–20 assuming available black hole mass measurements. In addition, our models infer that a significant fraction of the disc energy in these ultra-soft NLS1s is radiated away in the form of non-thermal emission instead of the thermal emission from the disc. Due to their extreme properties, X-ray observations of these sources in the iron band are particularly challenging. Future observations, e.g. from Athena, will enable us to have a clearer view of the spectral shape in the iron band and thus distinguish the reflection model from other interpretations of their broad-band spectra.

AT 2017gbl: a dust obscured TDE candidate in a luminous infrared galaxy

ABSTRACT We present the discovery with Keck of the extremely infrared (IR) luminous transient AT 2017gbl, coincident with the Northern nucleus of the luminous infrared galaxy (LIRG) IRAS 23436+5257. Our extensive multiwavelength follow-up spans ∼900 d, including photometry and spectroscopy in the optical and IR, and (very long baseline interferometry) radio and X-ray observations. Radiative transfer modelling of the host galaxy spectral energy distribution and long-term pre-outburst variability in the mid-IR indicate the presence of a hitherto undetected dust obscured active galactic nucleus (AGN). The optical and near-IR spectra show broad ∼2000 km s−1 hydrogen, He i, and O i emission features that decrease in flux over time. Radio imaging shows a fast evolving compact source of synchrotron emission spatially coincident with AT 2017gbl. We infer a lower limit for the radiated energy of 7.3 × 1050 erg from the IR photometry. An extremely energetic supernova would satisfy this budget, but is ruled out by the radio counterpart evolution. Instead, we propose AT 2017gbl is related to an accretion event by the central supermassive black hole, where the spectral signatures originate in the AGN broad line region and the IR photometry is consistent with re-radiation by polar dust. Given the fast evolution of AT 2017gbl, we deem a tidal disruption event (TDE) of a star a more plausible scenario than a dramatic change in the AGN accretion rate. This makes AT 2017gbl the third TDE candidate to be hosted by a LIRG, in contrast to the so far considered TDE population discovered at optical wavelengths and hosted preferably by post-starburst galaxies.

Molecular gas and star formation activity in luminous infrared galaxies in clusters at intermediate redshifts

We investigate the role of dense megaparsec-scale environments in processing molecular gas of cluster galaxies as they fall into the cluster cores. We selected a sample of ∼20 luminous infrared galaxies (LIRGs) belonging to intermediate-redshift clusters, mainly from the Herschel Lensing Survey and the Local Cluster Substructure Survey. These galaxies include MACS J0717.5+3745 at z = 0.546 and Abell 697, 963, 1763, and 2219 at z = 0.2 − 0.3. We performed spectral energy distribution modeling from the far-infrared to ultraviolet of the LIRGs, which span cluster-centric distances within r/r200 ≃ 0.2 − 1.6. We observed the LIRGs in CO(1→0) or CO(2→1) with the Plateau de Bure interferometer and its successor NOEMA, as part of five observational programs carried out between 2012 and 2017. We compared the molecular gas to stellar mass ratio M(H2)/M⋆, star formation rate (SFR), and depletion time (τdep) of the LIRGs with those of a compilation of cluster and field star-forming galaxies from the literature. The targeted LIRGs have SFR, M(H2)/M⋆, and τdep that are consistent with those of both main-sequence (MS) field galaxies and star-forming galaxies from the comparison sample. However we find that the depletion time, normalized to the MS value, tentatively increases with increasing r/r200, with a significance of 2.8σ, which is ultimately due to a deficit of cluster-core LIRGs with τdep ≳ τdep, MS. We suggest that a rapid exhaustion of the molecular gas reservoirs occurs in the cluster LIRGs and is indeed effective in suppressing their star formation and ultimately quenching them. This mechanism may explain the exponential decrease of the fraction of cluster LIRGs with cosmic time. The compression of the gas in LIRGs, possibly induced by intra-cluster medium shocks, may be responsible for the short timescales that are observed in a large fraction of cluster-core LIRGs. Some of our LIRGs may also belong to a population of infalling filament galaxies.

A Census of Star Formation in the Outer Galaxy. II. The GLIMPSE360 Field.

We have conducted a study of star formation in the outer Galaxy from 65°< l < 265°in the region observed by the GLIMPSE360 program. This Spitzer warm mission program mapped the plane of the outer Milky Way with IRAC at 3.6 and 4.5 μm. We combine the IRAC, Wide-field Infrared Survey Explorer (WISE), and Two Micron All Sky Survey catalogs and our previous results from another outer Galaxy survey and identify a total of 47,338 young stellar objects (YSOs) across the field spanning >180° in Galactic longitude. Using the DBSCAN method on the combined catalog, we identify 618 clusters or aggregations of YSOs having five or more members. We identify 10,476 class I, 29,604 class II, and 7325 anemic class II/class III YSOs. The ratio of YSOs identified as members of clusters was 25,528/47,338, or 54%. We found that 100 of the clusters identified have previously measured distances in the WISE H II survey. We used these distances in our spectral energy distribution (SED) fitting of the YSOs in these clusters, of which 96 had YSOs with <3σ fits. We used the derived masses from the SED model fits to estimate the initial mass function (IMF) in the inner and outer Galaxy clusters; dividing the clusters by galactocentric distances, the slopes were Γ = 1.87 ± 0.31 above 3 M ⵙ for R Gal < 11.5 kpc and Γ = 1.15 ± 0.24 above 3 M ⵙ for R Gal > 11.5 kpc. The slope of the combined IMF was found to be Γ = 1.92 ± 0.42 above 3 M ⵙ. These values are consistent with each other within the uncertainties and with literature values in the inner Galaxy high-mass star formation regions. The slopes are likely also consistent with a universal Salpeter IMF.

Revisiting the Impact of Dust Production from Carbon-rich Wolf–Rayet Binaries

Abstract We present a dust spectral energy distribution (SED) and binary stellar population analysis revisiting the dust production rates (DPRs) in the winds of carbon-rich Wolf–Rayet (WC) binaries and their impact on galactic dust budgets. DustEM SED models of 19 Galactic WC “dustars” reveal DPRs of M ⊙ yr−1 and carbon dust condensation fractions, χ C , between 0.002% and 40%. A large (0.1–1.0 μm) dust grain size composition is favored for efficient dustars where χ C ≳ 1%. Results for dustars with known orbital periods verify a power-law relation between χ C , orbital period, WC mass-loss rate, and wind velocity consistent with predictions from theoretical models of dust formation in colliding-wind binaries. We incorporated dust production into Binary Population and Spectral Synthesis (BPASS) models to analyze dust production rates from WC dustars, asymptotic giant branch stars (AGBs), red supergiants (RSGs), and core-collapse supernovae (SNe). BPASS models assuming constant star formation (SF) and a coeval 106 M ⊙ stellar population were performed at low, Large Magellanic Cloud (LMC)–like, and solar metallicities (Z = 0.001, 0.008, and 0.020). Both constant SF and coeval models indicate that SNe are net dust destroyers at all metallicities. Constant SF models at LMC-like metallicities show that AGB stars slightly outproduce WC binaries and RSGs by factors of 2–3, whereas at solar metallicities WC binaries are the dominant source of dust for ∼60 Myr until the onset of AGBs, which match the dust input of WC binaries. Coeval population models show that, for “bursty” SF, AGB stars dominate dust production at late times (t ≳ 70 Myr).

Constraining the Dimensionality of SN Ia Spectral Variation with Twins

Abstract SNe Ia continue to play a key role in cosmological measurements. Their interpretation over a range in redshift requires a rest-frame spectral energy distribution model. For practicality, these models are parameterized with a limited number of parameters and are trained using linear or nonlinear dimensionality reduction. This work focuses on the related problem of estimating the number of parameters underlying SN Ia spectral variation (the dimensionality). I present a technique for using the properties of high-dimensional space and the counting statistics of “twin” SNe Ia to estimate this dimensionality. Applying this method to the supernova pairings from Fakhouri et al. shows that a modest number of parameters (three to five, not including extinction) explain those data well. The analysis also finds that the intrinsic parameters are approximately Gaussian-distributed. The limited number of parameters hints that improved SED models are possible that may enable substantial reductions in SN cosmological uncertainties with current and near-term data sets.

Effective Opacity of the Intergalactic Medium from Galaxy Spectra Analysis

Abstract We measure the effective opacity ( ) of the intergalactic medium from the composite spectra of 281 Lyman-break galaxies in the redshift range 2 ≲ z ≲ 3. Our spectra are taken from the COSMOS Lyα Mapping And Tomographic Observations survey derived from the Low Resolution Imaging Spectrometer on the W.M. Keck I telescope. We generate composite spectra in two redshift intervals and fit them with spectral energy distribution (SED) models composed of simple stellar populations. Extrapolating these SED models into the Lyα forest, we measure the effective Lyα opacity ( ) in the 2.02 ≤ z ≤ 2.44 range. At z = 2.22, we estimate from a power-law fit to the data. These measurements are consistent with estimates from quasar analyses at z &lt; 2.5 indicating that the systematic errors associated with normalizing quasar continua are not substantial. We provide a Gaussian processes model of our results and previous measurements that describes the steep redshift evolution in from z = 1.5–4.

KELT-9 b’s Asymmetric TESS Transit Caused by Rapid Stellar Rotation and Spin–Orbit Misalignment

Abstract KELT-9 b is an ultra-hot Jupiter transiting a rapidly rotating, oblate early-A-type star in a polar orbit. We model the effect of rapid stellar rotation on KELT-9 b’s transit light curve using photometry from the Transiting Exoplanet Survey Satellite to constrain the planet’s true spin–orbit angle and to explore how KELT-9 b may be influenced by stellar gravity darkening. We constrain the host star’s equatorial radius to be 1.089 ± 0.017 times as large as its polar radius and its local surface brightness to vary by ∼38% between its hot poles and cooler equator. We model the stellar oblateness and surface brightness gradient and find that it causes the transit light curve to lack the usual symmetry around the time of minimum light. We take advantage of the light-curve asymmetry to constrain KELT-9 b’s true spin–orbit angle ( ), agreeing with Gaudi et al. that KELT-9 b is in a nearly polar orbit. We also apply a gravity-darkening correction to the spectral energy distribution model from Gaudi et al. and find that accounting for rapid rotation gives a better fit to available spectroscopy and yields a more reliable estimate for the star’s polar effective temperature.

Time-delay Measurement of Mg ii Broad-line Response for the Highly Accreting Quasar HE 0413-4031: Implications for the Mg ii–based Radius–Luminosity Relation

Abstract We present the monitoring of the active galactic nuclei continuum and Mg ii broad-line emission for the quasar HE 0413-4031 (z = 1.38) based on the six-year monitoring by the South African Large Telescope (SALT). We manage to estimate a time-delay of days in the rest frame of the source using seven different methods: interpolated cross-correlation function, discrete correlation function (DCF), z-transformed DCF, JAVELIN, two estimators of data regularity (Von Neumann, Bartels), and χ 2 method. This time-delay is below the value expected from the standard radius–luminosity relation. However, based on the monochromatic luminosity of the source and the spectral energy distribution modeling, we interpret this departure as the shortening of the time-delay due to the higher accretion rate of the source, with the inferred Eddington ratio of ∼0.4. The Mg ii line luminosity of HE 0413-4031 responds to the continuum variability as , which is consistent with the light-travel distance of the location of Mg ii emission at R out ∼ 1018 cm. Using the data of 10 other quasars, we confirm the radius–luminosity relation for the broad Mg ii line, which was previously determined for the broad Hβ line for lower-redshift sources. In addition, we detect a general departure of higher-accreting quasars from this relation in analogy to the Hβ sample. After the accretion-rate correction of the light-travel distance, the Mg ii–based radius–luminosity relation has a small scatter of only 0.10 dex.

No significant correlation between radial velocity planet presence and debris disc properties

ABSTRACT We investigate whether the tentative correlation between planets and debris discs which has been previously identified can be confirmed at high significance. We compile a sample of 201 stars with known planets and existing far-infrared observations. The sample is larger than those studied previously since we include targets from an unpublished Herschel survey of planet hosts. We use spectral energy distribution modelling to characterize Kuiper belt analogue debris discs within the sample, then compare the properties of the discs against a control sample of 294 stars without known planets. Survival analysis suggests that there is a significant (p ∼ 0.002) difference between the disc fractional luminosity distributions of the two samples. However, this is largely a result of the fact that the control sample contains a higher proportion of close binaries and of later-type stars; both of these factors are known to reduce disc detection rates. Considering only Sun-like stars without close binary companions in each sample greatly reduces the significance of the difference (p ∼ 0.3). We also find no evidence for a difference in the disc fractional luminosities of stars hosting planets more or less massive than Saturn (p ∼ 0.9). Finally, we find that the planet hosts have cooler discs than the control stars, but this is likely a detection bias, since the warmest discs in the control sample are also the faintest, and would thus be undetectable around the more distant planet hosts. Considering only discs in each sample that could have been detected around a typical planet host, we find p ∼ 0.07 for the temperatures.

An obscured AGN population hidden in the VIPERS galaxies: identification through spectral energy distribution decomposition

ABSTRACT The detection of X-ray emission constitutes a reliable and efficient tool for the selection of active galactic nuclei (AGNs), although it may be biased against the most heavily absorbed AGNs. Simple mid-infrared (IR) broad-band selection criteria identify a large number of luminous and absorbed AGNs, yet again host contamination could lead to non-uniform and incomplete samples. Spectral energy distribution (SED) decomposition is able to decouple the emission from the AGN versus that from star-forming regions, revealing weaker AGN components. We aim to identify the obscured AGN population in the VIMOS Public Extragalactic Redshift Survey in the Canada–France–Hawaii Telescope Legacy Survey W1 field through SED modelling. We construct SEDs for 6860 sources and identify 160 AGNs at a high confidence level using a Bayesian approach. Using optical spectroscopy, we confirm the nature of ∼85 per cent of the AGNs. Our AGN sample is highly complete (∼92 per cent) compared to mid-IR colour-selected AGNs, including a significant number of galaxy-dominated systems with lower luminosities. In addition to the lack of X-ray emission (80 per cent), the SED fitting results suggest that the majority of the sources are obscured. We use a number of diagnostic criteria in the optical, IR, and X-ray regimes to verify these results. Interestingly, only 35 per cent of the most luminous mid-IR-selected AGNs have X-ray counterparts suggesting strong absorption. Our work emphasizes the importance of using SED decomposition techniques to select a population of type II AGNs, which may remain undetected by either X-ray or IR colour surveys.

Gamma-Ray Flares in the Long-term Light Curve of 3C 454.3

Abstract 3C 454.3 is frequently observed in the flaring state. The long-term light curve of this source has been analyzed with 9 yr (2008 August–2017 July) of data from the Fermi-LAT detector. We have identified five flares and one quiescent state. The flares have substructures with many peaks during the flaring phase. We have estimated the rise and decay time of the flares and compared with flares of other similar sources. The modeling of gamma-ray spectral energy distributions shows in most cases that a log-parabola function gives the best fit to the data. We have done time-dependent leptonic modeling of two of the flares, for which simultaneous multiwavelength data are available. These two long-lasting flares, Flare-2A and Flare-2D, continued for 95 and 133 days, respectively. We have used the average values of Doppler factor, injected luminosity in electrons, size of the emission region, and the magnetic field in the emission region in modeling these flares. The emission region is assumed to be in the broad-line region in our single-zone model. The energy losses (synchrotron, synchrotron self-Compton, external Compton) and escape of electrons from the emission region have been included while doing the modeling. Although the total jet powers required to model these flares with the leptonic model are higher compared to other sources, they are always found to be lower than the Eddington luminosity of 3C 454.3. We also select some flaring peaks and show that the time variation of the Doppler factor or the injected luminosity in electrons over short timescales can explain their light curves.

SUGAR: An improved empirical model of Type Ia supernovae based on spectral features

Context.Type Ia supernovae (SNe Ia) are widely used to measure the expansion of the Universe. Improving distance measurements of SNe Ia is one technique to better constrain the acceleration of expansion and determine its physical nature.Aims.This document develops a new SNe Ia spectral energy distribution (SED) model, called the SUpernova Generator And Reconstructor (SUGAR), which improves the spectral description of SNe Ia, and consequently could improve the distance measurements.Methods.This model was constructed from SNe Ia spectral properties and spectrophotometric data from the Nearby Supernova Factory collaboration. In a first step, a principal component analysis-like method was used on spectral features measured at maximum light, which allowed us to extract the intrinsic properties of SNe Ia. Next, the intrinsic properties were used to extract the average extinction curve. Third, an interpolation using Gaussian processes facilitated using data taken at different epochs during the lifetime of an SN Ia and then projecting the data on a fixed time grid. Finally, the three steps were combined to build the SED model as a function of time and wavelength. This is the SUGAR model.Results.The main advancement in SUGAR is the addition of two additional parameters to characterize SNe Ia variability. The first is tied to the properties of SNe Ia ejecta velocity and the second correlates with their calcium lines. The addition of these parameters, as well as the high quality of the Nearby Supernova Factory data, makes SUGAR an accurate and efficient model for describing the spectra of normal SNe Ia as they brighten and fade.Conclusions.The performance of this model makes it an excellent SED model for experiments like the Zwicky Transient Facility, the Large Synoptic Survey Telescope, or the Wide Field Infrared Survey Telescope.

The path to Z And-type outbursts: The case of V426 Sagittae (HBHA 1704-05)

Context. The star V426 Sge (HBHA 1704-05), originally classified as an emission-line object and a semi-regular variable, brightened at the beginning of August 2018, showing signatures of a symbiotic star outburst. Aims. We aim to confirm the nature of V426 Sge as a classical symbiotic star, determine the photometric ephemeris of the light minima, and suggest the path from its 1968 symbiotic nova outburst to the following 2018 Z And-type outburst. Methods. We re-constructed an historical light curve (LC) of V426 Sge from approximately the year 1900, and used original low- (R ∼ 500–1500; 330–880 nm) and high-resolution (R ∼ 11 000–34 000; 360–760 nm) spectroscopy complemented with Swift-XRT and UVOT, optical UBVRCIC and near-infrared JHKL photometry obtained during the 2018 outburst and the following quiescence. Results. The historical LC reveals no symbiotic-like activity from ∼1900 to 1967. In 1968, V426 Sge experienced a symbiotic nova outburst that ceased around 1990. From approximately 1972, a wave-like orbitally related variation with a period of 493.4 ± 0.7 days developed in the LC. This was interrupted by a Z And-type outburst from the beginning of August 2018 to the middle of February 2019. At the maximum of the 2018 outburst, the burning white dwarf (WD) increased its temperature to ≳2 × 105 K, generated a luminosity of ∼7 × 1037 (d/3.3 kpc)2 erg s−1 and blew a wind at the rate of ∼3 × 10−6 M⊙ yr−1. Our spectral energy distribution models from the current quiescent phase reveal that the donor is a normal M4-5 III giant characterised with Teff ∼ 3400 K, RG ∼ 106 (d/3.3 kpc) R⊙ and LG ∼ 1350 (d/3.3 kpc)2 L⊙ and the accretor is a low-mass ∼0.5 M⊙ WD. Conclusions. During the transition from the symbiotic nova outburst to the quiescent phase, a pronounced sinusoidal variation along the orbit develops in the LC of most symbiotic novae. The following eventual outburst is of Z And-type, when the accretion by the WD temporarily exceeds the upper limit of the stable burning. At this point the system becomes a classical symbiotic star.