ABSTRACT
Heating from active galactic nuclei (AGNs) is thought to stabilize cool-core clusters, limiting star formation and cooling flows. We employ radiative magnetohydrodynamic (MHD) simulations to ...model light AGN jet feedback with different accretion modes (Bondi–Hoyle–Lyttleton and cold accretion) in an idealized Perseus-like cluster. Independent of the probed accretion model, accretion efficiency, jet density and resolution, the cluster self-regulates with central entropies and cooling times consistent with observed cool-core clusters in this non-cosmological setting. We find that increased jet efficiencies lead to more intermittent jet powers and enhanced star formation rates. Our fiducial low-density jets can easily be deflected by orbiting cold gaseous filaments, which redistributes angular momentum and leads to more extended cold gas distributions and isotropic bubble distributions. In comparison to our fiducial low momentum-density jets, high momentum-density jets heat less efficiently and enable the formation of a persistent cold gas disc perpendicular to the jets that is centrally confined. Cavity luminosities measured from our simulations generally reflect the cooling luminosities of the intracluster medium and correspond to averaged jet powers that are relatively insensitive to short periods of low-luminosity jet injection. Cold gas structures in our MHD simulations with low momentum-density jets generally show a variety of morphologies ranging from discy to very extended filamentary structures. In particular, magnetic fields are crucial to inhibit the formation of unrealistically massive cold gas discs by redistributing angular momentum between the hot and cold phases and by fostering the formation of elongated cold filaments that are supported by magnetic pressure.
Abstract
New observational facilities are probing astrophysical transients such as stellar explosions and gravitational-wave sources at ever-increasing redshifts, while also revealing new features in ...source property distributions. To interpret these observations, we need to compare them to predictions from stellar population models. Such models require the metallicity-dependent cosmic star formation history (
(
Z
,
z
)
) as an input. Large uncertainties remain in the shape and evolution of this function. In this work, we propose a simple analytical function for
(
Z
,
z
)
. Variations of this function can be easily interpreted because the parameters link to its shape in an intuitive way. We fit our analytical function to the star-forming gas of the cosmological TNG100 simulation and find that it is able to capture the main behavior well. As an example application, we investigate the effect of systematic variations in the
(
Z
,
z
)
parameters on the predicted mass distribution of locally merging binary black holes. Our main findings are that (i) the locations of features are remarkably robust against variations in the metallicity-dependent cosmic star formation history, and (ii) the low-mass end is least affected by these variations. This is promising as it increases our chances of constraining the physics that govern the formation of these objects (
https://github.com/LiekeVanSon/SFRD_fit/tree/7348a1ad0d2ed6b78c70d5100fb3cd2515493f02/
).
The pattern speed with which galactic bars rotate is intimately linked to the amount of dark matter in the inner regions of their host galaxies. In particular, dark matter haloes act to slow down ...bars via torques exerted through dynamical friction. Observational studies of barred galaxies tend to find that bars rotate fast, while hydrodynamical cosmological simulations of galaxy formation and evolution in the Lambda cold dark matter (ΛCDM) framework have previously found that bars slow down excessively. This has led to a growing tension between fast bars and the ΛCDM cosmological paradigm. In this study we revisit this issue, using the Auriga suite of high-resolution, magneto-hydrodynamical cosmological zoom-in simulations of galaxy formation and evolution in the ΛCDM framework, finding that bars remain fast down to
z
= 0. In Auriga, bars form in galaxies that have higher stellar-to-dark matter ratios and are more baryon-dominated than in previous cosmological simulations; this suggests that in order for bars to remain fast, massive spiral galaxies must lie above the commonly used abundance matching relation. While this reduces the aforementioned tension between the rotation speed of bars and ΛCDM, it accentuates the recently reported discrepancy between the dynamically inferred stellar-to-dark matter ratios of massive spirals and those inferred from abundance matching. Our results highlight the potential of using bar dynamics to constrain models of galaxy formation and evolution.
2002cx-like supernovae are a sub-class of sub-luminous Type Ia supernovae (SNe). Their light curves and spectra are characterized by distinct features that indicate strong mixing of the explosion ...ejecta. Pure turbulent deflagrations have been shown to produce such mixed ejecta. Here, we present hydrodynamics, nucleosynthesis and radiative-transfer calculations for a 3D full-star deflagration of a Chandrasekhar-mass white dwarf. Our model is able to reproduce the characteristic observational features of SN 2005hk (a prototypical 2002cx-like supernova), not only in the optical, but also in the near-infrared. For that purpose we present, for the first time, five near-infrared spectra of SN 2005hk from −0.2 to 26.6 d with respect to B-band maximum. Since our model burns only small parts of the initial white dwarf, it fails to completely unbind the white dwarf and leaves behind a bound remnant of ∼1.03 M - consisting mainly of unburned carbon and oxygen, but also enriched by some amount of intermediate-mass and iron-group elements from the explosion products that fall back on the remnant. We discuss possibilities for detecting this bound remnant and how it might influence the late-time observables of 2002cx-like SNe.
We present a machine-learning (ML) approach for estimating galaxy cluster masses from Chandra mock images. We utilize a Convolutional Neural Network (CNN), a deep ML tool commonly used in image ...recognition tasks. The CNN is trained and tested on our sample of 7896 Chandra X-ray mock observations, which are based on 329 massive clusters from the simulation. Our CNN learns from a low resolution spatial distribution of photon counts and does not use spectral information. Despite our simplifying assumption to neglect spectral information, the resulting mass values estimated by the CNN exhibit small bias in comparison to the true masses of the simulated clusters (−0.02 dex) and reproduce the cluster masses with low intrinsic scatter, 8% in our best fold and 12% averaging over all. In contrast, a more standard core-excised luminosity method achieves 15%-18% scatter. We interpret the results with an approach inspired by Google DeepDream and find that the CNN ignores the central regions of clusters, which are known to have high scatter with mass.
In the version of the single-degenerate scenario of Type Ia supernovae (SNe Ia) studied here, a carbon-oxygen white dwarf explodes close to the Chandrasekhar limit after accreting material from a ...non-degenerate helium (He) companion star. In the present study, we employ the Stellar GADGET code to perform three-dimensional hydrodynamical simulations of the interaction of the SN Ia ejecta with the He companion star taking into account its orbital motion and spin. It is found that only 2%-5% of the initial companion mass is stripped off from the outer layers of He companion stars due to the supernova (SN) impact. The dependence of the unbound mass (or the kick velocity) on the orbital separation can be fitted to a good approximation by a power law for a given companion model. After the SN impact, the outer layers of a He donor star are significantly enriched with heavy elements from the low-expansion-velocity tail of SN Ia ejecta. The total mass of accumulated SN-ejecta material on the companion surface reaches about greater than or equal to 10 super(-3) M. for different companion models. This enrichment with heavy elements provides a potential way to observationally identify the surviving companion star in SN remnants. Finally, by artificially adjusting the explosion energy of the W7 explosion model, we find that the total accumulation of SN ejecta on the companion surface is also dependent on the explosion energy with a power-law relation to a good approximation.
SN 2014J in M82 is the closest detected Type Ia supernova (SN Ia) in at least 28 yr and perhaps in 410 yr. Despite its small distance of 3.3 Mpc, SN 2014J is surprisingly faint, peaking at ...V = 10.6 mag, and assuming a typical SN Ia luminosity, we infer an observed visual extinction of A
V
= 2.0 ± 0.1 mag. But this picture, with R
V
= 1.6 ± 0.2, is too simple to account for all observations. We combine 10 epochs (spanning a month) of HST/Space Telescope Imaging Spectrograph (STIS) ultraviolet through near-infrared spectroscopy with HST/Wide Field Camera 3 (WFC3), Katzman Automatic Imaging Telescope, and FanCam photometry from the optical to the infrared and nine epochs of high-resolution TRES (Tillinghast Reflection Echelle Spectrograph) spectroscopy to investigate the sources of extinction and reddening for SN 2014J. We argue that the wide range of observed properties for SN 2014J is caused by a combination of dust reddening, likely originating in the interstellar medium of M82, and scattering off circumstellar material. For this model, roughly half of the extinction is caused by reddening from typical dust (E(B − V) = 0.45 mag and R
V
= 2.6) and roughly half by scattering off Large Magellanic Cloud-like dust in the circumstellar environment of SN 2014J.
In a companion paper, Seitenzahl et al. have presented a set of three-dimensional delayed detonation models for thermonuclear explosions of near-Chandrasekhar-mass white dwarfs (WDs). Here, we ...present multidimensional radiative transfer simulations that provide synthetic light curves and spectra for those models. The model sequence explores both changes in the strength of the deflagration phase (which is controlled by the ignition configuration in our models) and the WD central density. In agreement with previous studies, we find that the strength of the deflagration significantly affects the explosion and the observables. Variations in the central density also have an influence on both brightness and colour, but overall it is a secondary parameter in our set of models. In many respects, the models yield a good match to the observed properties of normal Type Ia supernovae (SNe Ia): peak brightness, rise/decline time-scales and synthetic spectra are all in reasonable agreement. There are, however, several differences. In particular, the models are systematically too red around maximum light, manifest spectral line velocities that are a little too high and yield I-band light curves that do not match observations. Although some of these discrepancies may simply relate to approximations made in the modelling, some pose real challenges to the models. If viewed as a complete sequence, our models do not reproduce the observed light-curve width-luminosity relation (WLR) of SNe Ia: all our models show rather similar B-band decline rates, irrespective of peak brightness. This suggests that simple variations in the strength of the deflagration phase in Chandrasekhar-mass deflagration-to-detonation models do not readily explain the observed diversity of normal SNe Ia. This may imply that some other parameter within the Chandrasekhar-mass paradigm is key to the WLR, or that a substantial fraction of normal SNe Ia arise from an alternative explosion scenario.
We investigate the brightness distribution expected for thermonuclear explosions that might result from the ignition of a detonation during the violent merger of white dwarf (WD) binaries. Violent WD ...mergers are a subclass of the canonical double degenerate scenario where two carbon-oxygen (CO) WDs merge when the larger WD fills its Roche lobe. Determining their brightness distribution is critical for evaluating whether such an explosion model could be responsible for a significant fraction of the observed population of Type Ia supernovae (SNe Ia). We argue that the brightness of an explosion realized via the violent merger model is mainly determined by the mass of 56Ni produced in the detonation of the primary CO WD. To quantify this link, we use a set of sub-Chandrasekhar mass WD detonation models to derive a relationship between primary WD mass (m
WD) and expected peak bolometric brightness (M
bol). We use this m
WD-M
bol relationship to convert the masses of merging primary WDs from binary population models to a predicted distribution of explosion brightness. We also investigate the sensitivity of our results to assumptions about the conditions required to realize a detonation during violent mergers of WDs. We find a striking similarity between the shape of our theoretical peak-magnitude distribution and that observed for SNe Ia: our model produces a M
bol distribution that roughly covers the range and matches the shape of the one observed for SNe Ia. However, this agreement hinges on a particular phase of mass accretion during binary evolution: the primary WD gains ∼0.15-0.35 M from a slightly evolved helium star companion. In our standard binary evolution model, such an accretion phase is predicted to occur for about 43 per cent of all binary systems that ultimately give rise to binary CO WD mergers. We also find that with high probability, violent WD mergers involving the most massive primaries ( 1.3 M, which should produce bright SNe) have delay times 500 Myr.
Calculations of synthetic spectropolarimetry are one means to test multidimensional explosion models for Type Ia supernovae. In a recent paper, we demonstrated that the violent merger of a 1.1 and ...0.9 M⊙ white dwarf binary system is too asymmetric to explain the low polarization levels commonly observed in normal Type Ia supernovae. Here, we present polarization simulations for two alternative scenarios: the sub-Chandrasekhar mass double-detonation and the Chandrasekhar mass delayed-detonation model. Specifically, we study a 2D double-detonation model and a 3D delayed-detonation model, and calculate polarization spectra for multiple observer orientations in both cases. We find modest polarization levels (<1 per cent) for both explosion models. Polarization in the continuum peaks at ∼0.1–0.3 per cent and decreases after maximum light, in excellent agreement with spectropolarimetric data of normal Type Ia supernovae. Higher degrees of polarization are found across individual spectral lines. In particular, the synthetic Si ii λ6355 profiles are polarized at levels that match remarkably well the values observed in normal Type Ia supernovae, while the low degrees of polarization predicted across the O i λ7774 region are consistent with the non-detection of this feature in current data. We conclude that our models can reproduce many of the characteristics of both flux and polarization spectra for well-studied Type Ia supernovae, such as SN 2001el and SN 2012fr. However, the two models considered here cannot account for the unusually high level of polarization observed in extreme cases such as SN 2004dt.