Spatial field reconstruction with INLA Smole, Majda; Rino-Silvestre, João; González-Gaitán, Santiago ...
Astronomy and astrophysics (Berlin),
01/2023, Volume:
669
Journal Article
Peer reviewed
Open access
Aims.
Monte Carlo radiative transfer (MCRT) simulations are a powerful tool for understanding the role of dust in astrophysical systems and its influence on observations. However, due to the strong ...coupling of the radiation field and medium across the whole computational domain, the problem is non-local and non-linear, and such simulations are computationally expensive in the case of realistic 3D inhomogeneous dust distributions. We explore a novel technique for post-processing MCRT output to reduce the total computational run time by enhancing the output of computationally less expensive simulations of lower-quality.
Methods.
We combined principal component analysis (PCA) and non-negative matrix factorisation (NMF) as dimensionality reduction techniques together with Gaussian Markov random fields and the integrated nested Laplace approximation (INLA), an approximate method for Bayesian inference, to detect and reconstruct the non-random spatial structure in the images of lower signal-to-noise ratios or with missing data.
Results.
We tested our methodology using synthetic observations of a galaxy from the SKIRT Auriga project - a suite of high-resolution magnetohydrodynamic Milky Way-sized galaxies simulated in cosmological environment using a ‘zoom-in' technique. With this approach, we are able to reproduce high-photon-number reference images ~5 times faster with median residuals below ~20%.
ABSTRACT
We present polarimetric maps of the Circinus galaxy nucleus in the BVRI bands, obtained with VLT/FORS2. Circinus is the closest Seyfert 2 galaxy and harbours an archetypal obscured active ...galactic nucleus (AGN). Recent high angular resolution imaging revealed that a major fraction of its mid-infrared (MIR) emission is coming from the polar region. Previously, we demonstrated that these observations are consistent with a model of a compact dusty disc and a hyperboloid shell, resembling a hollow cone on larger scales. Here, we focus on the AGN core, up to 40 pc from the central engine, and compare the observations to the radiative transfer models. Polarization maps reveal a conical structure, coinciding with the ionization cone. The wavelength-dependence of the polarization degree indicates that scattering on dust grains is producing polarization. The observed polarization degree ($\approx 1-3{{\ \rm per\ cent}}$) is lower than predicted by the models; however, this is only a lower limit, since stellar emission dominates the total flux in the optical. The observed polarization angle (≈30°) is reproduced by the model of a dusty disc with a hollow cone that is illuminated by a tilted anisotropic central source. An accretion disc aligned with the ionization cone axis, and alternative dust geometries, such as a paraboloid shell, or a torus enveloped by ambient dust, are inconsistent with the data. We conclude that the optical polarimetric imaging supports earlier evidence for the presence of dust in the polar region, tentatively associated with dusty outflows.
We present an analysis of the diversity of V-band light-curves of hydrogen-rich type II supernovae. Analyzing a sample of 116 supernovae, several magnitude measurements are defined, together with ...decline rates at different epochs, and time durations of different phases. It is found that magnitudes measured at maximum light correlate more strongly with decline rates than those measured at other epochs: brighter supernovae at maximum generally have faster declining light-curves at all epochs. We find a relation between the decline rate during the "plateau" phase and peak magnitudes, which has a dispersion of 0.56 mag, offering the prospect of using type II supernovae as purely photometric distance indicators. Our analysis suggests that the type II population spans a continuum from low-luminosity events which have flat light-curves during the "plateau" stage, through to the brightest events which decline much faster. A large range in optically thick phase durations is observed, implying a range in progenitor envelope masses at the epoch of explosion. During the radioactive tails, we find many supernovae with faster declining light-curves than expected from full trapping of radioactive emission, implying low mass ejecta. It is suggested that the main driver of light-curve diversity is the extent of hydrogen envelopes retained before explosion. Finally, a new classification scheme is introduced where hydrogen-rich events are typed as simply "SN II" with an "s2" value giving the decline rate during the "plateau" phase, indicating its morphological type.
Abstract
We present the largest and most homogeneous collection of near-infrared (NIR) spectra of Type Ia supernovae (SNe Ia): 339 spectra of 98 individual SNe obtained as part of the Carnegie ...Supernova Project-II. These spectra, obtained with the FIRE spectrograph on the 6.5 m Magellan Baade telescope, have a spectral range of 0.8–2.5
μ
m. Using this sample, we explore the NIR spectral diversity of SNe Ia and construct a template of spectral time series as a function of the light-curve-shape parameter, color stretch
s
BV
. Principal component analysis is applied to characterize the diversity of the spectral features and reduce data dimensionality to a smaller subspace. Gaussian process regression is then used to model the subspace dependence on phase and light-curve shape and the associated uncertainty. Our template is able to predict spectral variations that are correlated with
s
BV
, such as the hallmark NIR features: Mg
ii
at early times and the
H
-band break after peak. Using this template reduces the systematic uncertainties in
K
-corrections by ∼90% compared to those from the Hsiao template. These uncertainties, defined as the mean
K
-correction differences computed with the color-matched template and observed spectra, are on the level of 4 × 10
−4
mag on average. This template can serve as the baseline spectral energy distribution for light-curve fitters and can identify peculiar spectral features that might point to compelling physics. The results presented here will substantially improve future SN Ia cosmological experiments, for both nearby and distant samples.
It is well known that massive stars (M > 8 M☉) evolve up to the collapse of the stellar core, resulting in most cases in a supernova (SN) explosion. Their heterogeneity is related mainly to different ...configurations of the progenitor star at the moment of the explosion and to their immediate environments. We present photometry and spectroscopy of SN 2010bt, which was classified as a Type IIn SN from a spectrum obtained soon after discovery and was observed extensively for about 2 months. After the seasonal interruption owing to its proximity to the Sun, the SN was below the detection threshold, indicative of a rapid luminosity decline. We can identify the likely progenitor with a very luminous star (log L/L☉ 7) through comparison of Hubble Space Telescope images of the host galaxy prior to explosion with those of the SN obtained after maximum light. Such a luminosity is not expected for a quiescent star, but rather for a massive star in an active phase. This progenitor candidate was later confirmed via images taken in 2015 (∼5 yr post-discovery), in which no bright point source was detected at the SN position. Given these results and the SN behavior, we conclude that SN 2010bt was likely a Type IIn SN and that its progenitor was a massive star that experienced an outburst shortly before the final explosion, leading to a dense H-rich circumstellar environment around the SN progenitor.
Studying galaxies at different cosmic epochs entails several observational effects that need to be taken into account to compare populations across a large time-span in a consistent manner. We use a ...sample of 166 nearby galaxies that hosted type Ia supernovae (SNe Ia) and have been observed with the integral field spectrograph MUSE as part of the AMUSING survey. Here, we present a study of the systematic errors and bias on the host stellar mass with increasing redshift, which are generally overlooked in SNe Ia cosmological analyses. We simulate observations at different redshifts (0.1 <
z
< 2.0) using four photometric bands (
griz
, similar to the Dark Energy Survey-SN program) to then estimate the host galaxy properties across cosmic time. We find that stellar masses are systematically underestimated as we move towards higher redshifts, due mostly to different rest-frame wavelength coverage, with differences reaching 0.3 dex at
z
∼ 1. We used the newly derived corrections as a function of redshift to correct the stellar masses of a known sample of SN Ia hosts and derive cosmological parameters. We show that these corrections have a small impact on the derived cosmological parameters. The most affected is the value of the mass step Δ
M
, which is reduced by ∼0.004 (6% lower). The dark energy equation of state parameter
w
changes by Δ
w
∼ 0.006 (0.6% higher) and the value of Ω
m
increases at most by 0.001 (∼0.3%), all within the derived uncertainties of the model. While the systematic error found in the estimate of the host stellar mass does not significantly affect the derived cosmological parameters, it is an important source of systematic error that needs to be corrected for as we enter a new era of precision cosmology.
We report results from the Supernova Photometric Classification Challenge (SNPhotCC), a publicly released mix of simulated supernovae (SNe), with types (Ia, Ibc, and II) selected in proportion to ...their expected rates. The simulation was realized in the
griz
g
r
i
z
filters of the Dark Energy Survey (DES) with realistic observing conditions (sky noise, point-spread function, and atmospheric transparency) based on years of recorded conditions at the DES site. Simulations of non–Ia-type SNe are based on spectroscopically confirmed light curves that includeunpublishednon-Ia samples donated from the Carnegie Supernova Project (CSP), the Supernova Legacy Survey (SNLS), and the Sloan Digital Sky Survey-II (SDSS-II). A spectroscopically confirmed subset was provided for training. We challenged scientists to run their classification algorithms and report a type and photo-
z
z
for each SN. Participants from 10 groups contributed 13 entries for the sample that included a host-galaxy photo-
z
z
for each SN and nine entries for the sample that had no redshift information. Several different classification strategies resulted in similar performance, and for all entries the performance was significantly better for the training subset than for the unconfirmed sample. For the spectroscopically unconfirmed subset, the entry with the highest average figure of merit for classifying SNe Ia has an efficiency of 0.96 and an SN Ia purity of 0.79. As a public resource for the future development of photometric SN classification and photo-
z
z
estimators, we have released updated simulations with improvements based on our experience from the SNPhotCC, added samples corresponding to the Large Synoptic Survey Telescope (LSST) and the SDSS-II, and provided the answer keys so that developers can evaluate their own analysis.
Abstract
Some supernovae, such as pair-instability supernovae, are predicted to have a duration of more than a year in the observer frame. To constrain the rates of supernovae lasting for more than a ...year, we conducted a long-term deep transient survey using Hyper Suprime-Cam (HSC) on the 8.2 m Subaru telescope. HSC is a wide-field (a 1.75 deg
2
field-of-view) camera and it can efficiently conduct transient surveys. We observed the same 1.75 deg
2
field repeatedly using the
g-
,
r-
,
i-
, and
z-
band filters with the typical depth of 26 mag for four seasons (from late 2016 to early 2020). Using these data, we searched for transients lasting for more than a year. Two supernovae were detected in two continuous seasons, one supernova was detected in three continuous seasons, but no transients lasted for all four seasons searched. The discovery rate of supernovae lasting for more than a year with the typical limiting magnitudes of 26 mag is constrained to be
. All the long-lasting supernovae we found are likely Type IIn supernovae and our results indicate that about 40% of Type IIn supernovae have long-lasting light curves. No plausible pair-instability supernova candidates lasting for more than a year are discovered. By comparing the survey results and survey simulations, we constrain the luminous pair-instability supernova rate up to
z
≃ 3 is of the order of 100 Gpc
−3
yr
−1
at most, which is 0.01–0.1% of the core-collapse supernova rate.
We report the discovery of distant RR Lyrae stars, including the most distant known in the Milky Way, using data taken in the g-band with the Dark Energy Camera as part of the High cadence Transient ...Survey (HiTS; 2014 campaign). We detect a total of 173 RR Lyrae stars over a ∼120 deg2 area, including both known RR Lyrae and new detections. The heliocentric distances dH of the full sample range from 9 to >200 kpc, with 18 of them beyond 90 kpc. We identify three sub-groups of RR Lyrae as members of known systems: the Sextans dwarf spheroidal galaxy, for which we report 46 new discoveries, and the ultra-faint dwarf galaxies Leo IV and Leo V. Following an MCMC methodology, we fit spherical and ellipsoidal profiles of the form (R) ∼ Rn to the radial density distribution of RR Lyrae in the Galactic halo. The best fit corresponds to the spherical case, for which we obtain a simple power-law index of , consistent with recent studies made with samples covering shorter distances. The pulsational properties of the outermost RR Lyrae in the sample (dH > 90 kpc) differ from the ones in the halo population at closer distances. The distribution of the stars in a period-amplitude diagram suggest they belong to Oosterhoff-intermediate or Oosterhoff II groups, similar to what is found in the ultra-faint dwarf satellites around the Milky Way. The new distant stars discovered represent an important addition to the few existing tracers of the Milky Way potential in the outer halo.
Supernova (SN) 2017cbv in NGC 5643 is one of a handful of Type Ia supernovae (SNe Ia) reported to have excess blue emission at early times. This paper presents extensive BVRIYJHKs-band light curves ...of SN 2017cbv, covering the phase from −16 to +125 days relative to B-band maximum light. The SN 2017cbv reached a B-band maximum of 11.710 0.006 mag, with a postmaximum magnitude decline of Δm15(B) = 0.990 0.013 mag. The SN suffered no host reddening based on Phillips intrinsic color, the Lira-Phillips relation, and the CMAGIC diagram. By employing the CMAGIC distance modulus = 30.58 0.05 mag and assuming H0 = 72 km s−1 Mpc−1, we found that 0.73 M 56Ni was synthesized during the explosion of SN 2017cbv, which is consistent with estimates using reddening- and distance-free methods via the phases of the secondary maximum of the near-IR- (NIR-) band light curves. We also present 14 NIR spectra from −18 to +49 days relative to the B-band maximum light, providing constraints on the amount of swept-up hydrogen from the companion star in the context of the single degenerate progenitor scenario. No Paβ emission feature was detected from our postmaximum NIR spectra, placing a hydrogen mass upper limit of 0.1 M . The overall optical/NIR photometric and NIR spectral evolution of SN 2017cbv is similar to that of a normal SN Ia, even though its early evolution is marked by a flux excess not seen in most other well-observed normal SNe Ia. We also compare the exquisite light curves of SN 2017cbv with some Mch delayed detonation models and sub-Mch double detonation models.