Simulations Of Galaxy Evolution Research Articles

Powering galactic superwinds with small-scale AGN winds

ABSTRACT We present a new implementation for active galactic nucleus (AGN) feedback through small-scale, ultrafast winds in the moving-mesh hydrodynamic code arepo. The wind is injected by prescribing mass, momentum, and energy fluxes across a spherical boundary centred on a supermassive black hole according to available constraints for accretion disc winds. After sweeping-up a mass equal to their own, small-scale winds thermalize, powering energy-driven outflows with dynamics, structure, and cooling properties in excellent agreement with those of analytic wind solutions. Momentum-driven solutions do not easily occur, because the Compton cooling radius is usually much smaller than the free-expansion radius of the small-scale winds. Through various convergence tests, we demonstrate that our implementation yields wind solutions, which are well converged down to the typical resolution achieved in cosmological simulations. We test our model in hydrodynamic simulations of isolated Milky Way – mass galaxies. Above a critical AGN luminosity, initially spherical, small-scale winds power bipolar, energy-driven superwinds that break out of the galactic nucleus, flowing at speeds $\gt 1000 \rm \, km \, s^{-1}$ out to $\sim 10 \, \rm kpc$. These energy-driven outflows result in moderate, but long-term, reduction in star formation, which becomes more pronounced for higher AGN luminosities and faster small-scale winds. Suppression of star formation proceeds through a rapid mode that involves the removal of the highest density, nuclear gas, and through a slower mode that effectively halts halo gas accretion. Our new implementation makes it possible to model AGN-driven winds in a physically meaningful and validated way in simulations of galaxy evolution, the interstellar medium and black hole accretion flows.

An Accreting, Anomalously Low-mass Black Hole at the Center of Low-mass Galaxy IC 750

Abstract We present a multiwavelength study of the active galactic nucleus in the nearby (D = 14.1 Mpc) low-mass galaxy IC 750, which has circumnuclear 22 GHz water maser emission. The masers trace a nearly edge-on, warped disk ∼0.2 pc in diameter, coincident with the compact nuclear X-ray source that lies at the base of the ∼kiloparsec-scale extended X-ray emission. The position–velocity structure of the maser emission indicates that the central black hole (BH) has a mass less than 1.4 × 105 M ⊙. Keplerian rotation curves fitted to these data yield enclosed masses between 4.1 × 104 M ⊙ and 1.4 × 105 M ⊙, with a mode of 7.2 × 104 M ⊙. Fitting the optical spectrum, we measure a nuclear stellar velocity dispersion km s−1. From near-infrared photometry, we fit a bulge mass of (7.3 ± 2.7) × 108 M ⊙ and a stellar mass of 1.4 × 1010 M ⊙. The mass upper limit of the intermediate-mass BH in IC 750 falls roughly two orders of magnitude below the M BH–σ * relation and roughly one order of magnitude below the M BH–M Bulge and M BH–M * relations—larger than the relations’ intrinsic scatters of 0.58 ± 0.09 dex, 0.69 dex, and 0.65 ± 0.09 dex, respectively. These offsets could be due to larger scatter at the low-mass end of these relations. Alternatively, BH growth is intrinsically inefficient in galaxies with low bulge and/or stellar masses, which causes the BHs to be undermassive relative to their hosts, as predicted by some galaxy evolution simulations.

Quantifying inhomogeneities in the HI distributions of simulated galaxies

The NIHAO cosmological simulations form a collection of a hundred high-resolution galaxies. We used these simulations to test the impact of stellar feedback on the morphology of the HI distribution in galaxies. We ran a subsample of twenty of the galaxies with different parameterizations of stellar feedback, looking for differences in the HI spatial distribution and morphology. We found that different feedback models do leave a signature in HI, and can potentially be compared with current and future observations. These findings can help inform future modeling efforts in the parameterization of stellar feedback in cosmological simulations of galaxy formation and evolution.

Metallicity gradients in small and nearby spiral galaxies

ABSTRACT Spectra of H ii regions obtained with Gemini/GMOS are used to derive the radial metallicity gradients of four small, low-mass spiral galaxies. The analysis of the outer disc of one of them, NGC 1058, uncovers the characteristic flattening found in similar extended disc galaxies. After combining these data with published long-slit observations of nearby spiral galaxies, no evidence for a dependence of the disc scale length-normalized metallicity gradients with stellar mass is found, down to log (M⋆/M⊙) ∼ 8.5. The abundance gradients derived from these observations are compared to predictions from recent cosmological simulations of galaxy evolution, finding that in several cases the simulations fail to reproduce the mean steepening of the gradients, expressed in dex kpc−1, with decreasing stellar mass for present-day galaxies, or do not extend to sufficiently small stellar masses for a meaningful comparison. The mean steepening of the abundance gradients (in dex kpc−1) with decreasing disc scale length is in qualitative agreement with predictions from the inside-out model of Boissier & Prantzos, although the predicted slopes are systematically steeper than observed. This indicates the necessity of including processes such as outflows and radial mixing in similar models of galactic chemical evolution. Published spatially resolved metallicity and photometric data of dwarf irregular galaxies suggest that significant, but transitory, metallicity gradients can develop for systems that have experienced recent (t < 100 Myr) enhanced star formation in their inner discs.

Population Estimates for Electromagnetically Distinguishable Supermassive Binary Black Holes

Abstract Distinguishing the photon output of an accreting supermassive black hole binary system from that of a single supermassive black hole accreting at the same rate is intrinsically difficult because the majority of the light emerges from near the innermost stable orbits of the black holes. However, there are two possible signals that can distinctively mark binaries, both arising from the gap formed in circumbinary accretion flows inside approximately twice the binary separation. One of these is a “notch” cut into the thermal spectra of these systems in the IR/optical/UV, the other a periodically varying excess hard X-ray luminosity whose period is of order the binary orbital period. Using data from detailed galaxy evolution simulations, we estimate the distribution function in mass, mass ratio, and accretion rate for accreting supermassive binary black holes (SMBBHs) as a function of redshift and then transform this distribution function into predicted source counts for these two potential signals. At flux levels ≳10−13 erg cm−2 s−1, there may be ∼O(102) such systems in the sky, mostly in the redshift range 0.5 ≲ z ≲ 1. Roughly 10% should have periods short enough (≲5 yr) to detect the X-ray modulation; this is also the period range accessible to Pulsar Timing Array observations.

Supernovae feedback propagation: the role of turbulence

Modelling the propagation of supernova (SN) bubbles, in terms of energy, momentum and spatial extent, is critical for simulations of galaxy evolution which do not capture these scales. To date, small scale models of SN feedback predict that the evolution of above-mentioned quantities can be solely parameterised by average quantities of the surrounding gas, such as density. However, most of these studies neglect the turbulent motions of this medium. In this paper, we study the propagation and evolution of SNe in turbulent environments. We confirm that the time evolution of injected energy and momentum can be characterised by the average density. However, the details of the density structure of the interstellar medium play a crucial role in the spatial extent of the bubble, even at a given average density. We demonstrate that spherically symmetric models of SN bubbles do not model well their spatial extent, and therefore cannot not be used to design sub-grid models of SNe feedback at galactic and cosmological scales.

Figuring Out Gas & Galaxies in Enzo (FOGGIE). I. Resolving Simulated Circumgalactic Absorption at 2 ≤ z ≤ 2.5

Abstract We present simulations from the new “Figuring Out Gas & Galaxies in Enzo” (FOGGIE) project. In contrast to most extant simulations of galaxy formation, which concentrate computational resources on galactic disks and spheroids with fluid and particle elements of fixed mass, the FOGGIE simulations focus on extreme spatial and mass resolution in the circumgalactic medium (CGM) surrounding galaxies. Using the Enzo code and a new refinement scheme, FOGGIE reaches spatial resolutions of 381 comoving h −1 pc and resolves extremely low masses (≲1–100 ) out to 100 comoving h −1 kpc from the central halo. At these resolutions, cloud and filament-like structures giving rise to simulated absorption are smaller, and better resolved, than the same structures simulated with standard density-dependent refinement. Most of the simulated absorption arises in identifiable and well-resolved structures with masses ≲104 , well below the mass resolution of typical zoom simulations. However, integrated quantities such as mass surface density and ionic covering fractions change at only the ≲30% level as resolution is varied. These relatively small changes in projected quantities—even when the sizes and distribution of absorbing clouds change dramatically—indicate that commonly used observables provide only weak constraints on the physical structure of the underlying gas. Comparing the simulated absorption features to the KODIAQ (Keck Observatory Database of Ionized Absorption toward Quasars) survey of z ∼ 2–3.5 Lyman limit systems, we show that high-resolution FOGGIE runs better resolve the internal kinematic structure of detected absorption and better match the observed distribution of absorber properties. These results indicate that circumgalactic medium resolution is key in properly testing simulations of galaxy evolution with circumgalactic observations.

The impact of black hole seeding in cosmological simulations

Most cosmological simulations of galaxy evolution include active galactic nucleus (AGN) feedback, typically seeding black holes with masses of $\geq 10^5\, h^{-1}\, \rm{M}_{\odot}$ when the dark matter halo exceeds a given threshold mass. Taylor & Kobayashi (2014) introduced a new model, which seeds black holes at $10^3\, h^{-1}\, \rm{M}_{\odot}$ based on gas properties alone, and motivated by the channel of black hole formation due to the collapse of the most massive first stars in the Universe. We compare the black hole mass when the dark matter halo mass is $10^{10}\, h^{-1}\, \rm{M}_{\odot}$ between the different seeding methods. We find that seeding based upon gas properties gives a distribution of black hole masses with $\langle \log M_{\rm{BH}} {/ \rm{M}_{\odot}} \rangle = (5.18 \pm 0.54)$ when dark matter halo mass is $10^{10}\, h^{-1}\, \rm{M}_{\odot}$, consistent with the {seeding criteria} used in other simulations. However, the evolution of individual galaxies can be strongly affected by the different seeding mechanisms. We also find that the mean value of the distribution of black hole masses at a given halo mass evolves over time, with higher masses at higher redshifts, indicative of downsizing. Our results can inform more physically motivated black hole and AGN feedback models in cosmological simulations and semi-analytic models.

Dynamic localized turbulent diffusion and its impact on the galactic ecosystem

Modelling the turbulent diffusion of thermal energy, momentum, and metals is required in all galaxy evolution simulations due to the ubiquity of turbulence in galactic environments. The most commonly employed diffusion model, the Smagorinsky model, is known to be over-diffusive due to its strong dependence on the fluid velocity shear. We present a method for dynamically calculating a more accurate, locally appropriate, turbulent diffusivity: the dynamic localised Smagorinsky model. We investigate a set of standard astrophysically-relevant hydrodynamical tests, and demonstrate that the dynamic model curbs over-diffusion in non-turbulent shear flows and improves the density contrast in our driven turbulence experiments. In galactic discs, we find that the dynamic model maintains the stability of the disc by preventing excessive angular momentum transport, and increases the metal-mixing timescale in the interstellar medium. In both our isolated Milky Way-like galaxies and cosmological simulations, we find that the interstellar and circumgalactic media are particularly sensitive to the treatment of turbulent diffusion. We also examined the global gas enrichment fractions in our cosmological simulations, to gauge the potential effect on the formation sites and population statistics of Population III stars and supermassive black holes, since they are theorised to be sensitive to the metallicity of the gas out of which they form. The dynamic model is, however, not for galaxy evolution studies only. It can be applied to all astrophysical hydrodynamics simulations, including those modelling stellar interiors, planetary formation, and star formation.

Stellar Radiation Is Critical for Regulating Star Formation and Driving Outflows in Low-mass Dwarf Galaxies

Abstract Effective stellar feedback is used in models of galaxy formation to drive realistic galaxy evolution. Models typically include energy injection from supernovae (SNe) as the dominant form of stellar feedback, often in some form of sub-grid recipe. However, it has been recently suggested that pre-SN feedback (stellar winds or radiation) is necessary in high-resolution simulations of galaxy evolution to properly regulate star formation and properties of the interstellar medium (ISM). Following these processes is computationally challenging, so many prescriptions model this feedback approximately, accounting for the local destruction of dense gas clouds around newly formed stars in lieu of a full radiative transfer calculation. In this Letter we examine high-resolution simulations (1.8 pc) of an isolated dwarf galaxy with detailed stellar feedback tracked on a star-by-star basis. By following stellar ionizing radiation with an adaptive ray-tracing radiative transfer method, we test its importance in regulating star formation and driving outflows in this galaxy. We find that including ionizing radiation reduces the star formation rate (SFR) by over a factor of 5, and is necessary to produce the ISM conditions needed for SNe to drive significant outflows. We find that a localized approximation for radiation feedback is sufficient to regulate the SFR on short timescales, but does not allow significant outflows. Short- and long-range radiation effects are both important in driving the evolution of our low-metallicity, low-mass dwarf galaxy. Generalizing these results to more massive galaxies would be a valuable avenue of future research.

Star formation at low rates: impact of lacking massive stars on the evolution of dwarf galaxies

AbstractIn recent years dedicated observations have uncovered star formation at extremely low rates in dwarf galaxies, tidal tails, ram-pressure stripped gas clouds, and the outskirts of galactic disks. At the same time, numerical simulations of galaxy evolution have advanced to higher spatial and mass resolutions, but have yet to account for the underfilling of the uppermost mass bins of stellar initial mass function (IMF) at low star-formation rates. In such situations, simulations may simply scale down the IMF, without realizing that this unrealistically results in fractions of massive stars, along with fractions of massive star feedback energy (e.g., radiation and SNII explosions). Not properly accounting for such parameters has consequences for the self-regulation of star formation, the energetics of galaxies, as well as for the evolution of chemical abundances. Here we present numerical simulations of dwarf galaxies with low star-formation rates allowing for two extreme cases of the IMF: a “filled” case with fractional massive stars vs. a truncated IMF, at which the IMF is built bottom-up until the gas reservoir allows the formation of a last single star at an uppermost mass. The aim of the study is to demonstrate the different effects on galaxy evolution with respect to self-regulation, feedback, and chemistry. The case of a stochastic sampled IMF is situated somewhere in between these extremes.

Understanding mechanical feedback from HERGs and LERGs

AbstractThe properties of ∼ 1000 high-excitation and low-excitation radio galaxies (HERGs and LERGs) selected from the Heywood et al. (2016) 1 – 2 GHz VLA survey of Stripe 82 are investigated. The HERGs in this sample are generally found in host galaxies with younger stellar populations than LERGs, consistent with other work. The HERGs tend to accrete at a faster rate than the LERGs, but there is more overlap in the accretion rates of the two classes than has been found previously. We find evidence that mechanical feedback may be significantly underestimated in hydrodynamical simulations of galaxy evolution; 84 % of this sample release more than 10 % of their energy in mechanical form. Mechanical feedback is significant for many of the HERGs in this sample as well as the LERGs; nearly 50 % of the HERGs release more than 10 % of their energy in their radio jets.

Recovering the Physical Properties of Molecular Gas in Galaxies from CO SLED Modeling

Abstract Modeling of the spectral line energy distribution (SLED) of the CO molecule can reveal the physical conditions (temperature and density) of molecular gas in Galactic clouds and other galaxies. Recently, the Herschel Space Observatory and ALMA have offered, for the first time, a comprehensive view of the rotational J = 4−3 through J = 13−12 lines, which arise from a complex, diverse range of physical conditions that must be simplified to one, two, or three components when modeled. Here we investigate the recoverability of physical conditions from SLEDs produced by galaxy evolution simulations containing a large dynamical range in physical properties. These simulated SLEDs were generally fit well by one component of gas whose properties largely resemble or slightly underestimate the luminosity-weighted properties of the simulations when clumping due to nonthermal velocity dispersion is taken into account. If only modeling the first three rotational lines, the median values of the marginalized parameter distributions better represent the luminosity-weighted properties of the simulations, but the uncertainties in the fitted parameters are nearly an order of magnitude, compared to approximately 0.2 dex in the “best-case” scenario of a fully sampled SLED through J = 10−9. This study demonstrates that while common CO SLED modeling techniques cannot reveal the underlying complexities of the molecular gas, they can distinguish bulk luminosity-weighted properties that vary with star formation surface densities and galaxy evolution, if a sufficient number of lines are detected and modeled.

The Many Routes to AGN Feedback

The energy released by Active Galactic Nuclei (AGN) in the form of radiation, winds, or radio plasma jets, is known to impact on the surrounding interstellar medium. The result of these processes, known as AGN (negative) feedback, is suggested to prevent gas, in and around galaxies, from cooling, and to remove, or at least redistribute, gas by driving massive and fast outflows, hence playing a key role in galaxy evolution. Given its importance, a large effort is devoted by the astronomical community to trace the effects of AGN on the surrounding gaseous medium and to quantify their impact for different types of AGN. This review briefly summarizes some of the recent observational results obtained in different wavebands, tracing different phases of the gas. I also summarize the new insights they have brought, and the constraints they provide to numerical simulations of galaxy formation and evolution. The recent addition of deep observations of cold gas and, in particular, of cold molecular gas, has brought some interesting surprises and has expanded our understanding of AGN and AGN feedback.

Quantifying the AGN-driven outflows in ULIRGs (QUADROS) – I: VLT/Xshooter observations of nine nearby objects

Although now routinely incorporated into hydrodynamic simulations of galaxy evolution, the true importance of the feedback effect of the outflows driven by active galactic nuclei (AGN) remains uncertain from an observational perspective. This is due to a lack of accurate information on the densities, radial scales and level of dust extinction of the outflow regions. Here we use the unique capabilities of VLT/Xshooter to investigate the warm outflows in a representative sample of 9 local (0.06 $<$ z $<$ 0.15) ULIRGs with AGN nuclei and, for the first time, accurately quantify the key outflow properties. We find that the outflows are compact (0.05 $<$ R$_{[OIII]}$ $<$ 1.2 kpc), significantly reddened (median E(B-V)$\sim$0.5 magnitudes), and have relatively high electron densities (3.4 $<$ log$_{10}$ n$_e$ (cm$^{-3}$) $<$ 4.8). It is notable that the latter densities -- obtained using trans-auroral [SII] and [OII] emission-line ratios -- exceed those typically assumed for the warm, emission-line outflows in active galaxies, but are similar to those estimated for broad and narrow absorption line outflow systems detected in some type 1 AGN. Even if we make the most optimistic assumptions about the true (deprojected) outflow velocities, we find relatively modest mass outflow rates ($0.07 < \dot{M} < 11$ M$_{sol}$ yr$^{-1}$) and kinetic powers measured as a fraction of the AGN bolometric luminosities ($4\times10^{-4} < \dot{E}/L_{BOL} <1$%). Therefore, although warm, AGN-driven outflows have the potential to strongly affect the star formation histories in the inner bulge regions ($r \sim$ 1kpc) of nearby ULIRGs, we lack evidence that they have a significant impact on the evolution of these rapidly evolving systems on larger scales.

CANDELS: Elevated Black Hole Growth in the Progenitors of Compact Quiescent Galaxies at z ∼ 2

Abstract We examine the fraction of massive ( ) compact star-forming galaxies (cSFGs) that host an active galactic nucleus (AGN) at . These cSFGs are likely the direct progenitors of the compact quiescent galaxies observed at this epoch, which are the first population of passive galaxies to appear in large numbers in the early Universe. We identify cSFGs that host an AGN using a combination of Hubble WFC3 imaging and Chandra X-ray observations in four fields: the Chandra Deep Fields, the Extended Groth Strip, and the UKIDSS Ultra Deep Survey field. We find that (65/166) of cSFGs at host an X-ray detected AGN. This fraction is 3.2 times higher than the incidence of AGN in extended star-forming galaxies with similar masses at these redshifts. This difference is significant at the level. Our results are consistent with models in which cSFGs are formed through a dissipative contraction that triggers a compact starburst and concurrent growth of the central black hole. We also discuss our findings in the context of cosmological galaxy evolution simulations that require feedback energy to rapidly quench cSFGs. We show that the AGN fraction peaks precisely where energy injection is needed to reproduce the decline in the number density of cSFGs with redshift. Our results suggest that the first abundant population of massive quenched galaxies emerged directly following a phase of elevated supermassive black hole growth and further hints at a possible connection between AGN and the rapid quenching of star formation in these galaxies.

Kiloparsec-scale Simulations of Star Formation in Disk Galaxies. IV. Regulation of Galactic Star Formation Rates by Stellar Feedback

Abstract Star formation from the interstellar medium of galactic disks is a basic process controlling the evolution of galaxies. Understanding the star formation rate (SFR) in a local patch of a disk with a given gas mass is thus an important challenge for theoretical models. Here we simulate a kiloparsec region of a disk, following the evolution of self-gravitating molecular clouds down to subparsec scales, as they form stars that then inject feedback energy by dissociating and ionizing UV photons and supernova explosions. We assess the relative importance of each feedback mechanism. We find that H2-dissociating feedback results in the largest absolute reduction in star formation compared to the run with no feedback. Subsequently adding photoionization feedback produces a more modest reduction. Our fiducial models that combine all three feedback mechanisms yield, without fine-tuning, SFRs that are in excellent agreement with observations, with H2-dissociating photons playing a crucial role. Models that only include supernova feedback—a common method in galaxy evolution simulations—settle to similar SFRs, but with very different temperatures and chemical states of the gas, and with very different spatial distributions of young stars.

A deconvolution technique to correct deep images of galaxies from instrumental scattered light

Deep imaging of the diffuse light that is emitted by stellar fine structures and outer halos around galaxies is often now used to probe their past mass assembly. Because the extended halos survive longer than the relatively fragile tidal features, they trace more ancient mergers. We use images that reach surface brightness limits as low as 28.5−29 mag arcsec-2 (g-band) to obtain light and color profiles up to 5–10 effective radii of a sample of nearby early-type galaxies. These were acquired with MegaCam as part of the CFHT MATLAS large programme. These profiles may be compared to those produced using simulations of galaxy formation and evolution, once corrected for instrumental effects. Indeed they can be heavily contaminated by the scattered light caused by internal reflections within the instrument. In particular, the nucleus of galaxies generates artificial flux in the outer halo, which has to be precisely subtracted. We present a deconvolution technique to remove the artificial halos that makes use of very large kernels. The technique, which is based on PyOperators, is more time efficient than the model-convolution methods that are also used for that purpose. This is especially the case for galaxies with complex structures that are hard to model. Having a good knowledge of the point spread function (PSF), including its outer wings, is critical for the method. A database of MegaCam PSF models corresponding to different seeing conditions and bands was generated directly from the deep images. We show that the difference in the PSFs in different bands causes artificial changes in the color profiles, in particular a reddening of the outskirts of galaxies having a bright nucleus. The method is validated with a set of simulated images and applied to three representative test cases: NGC 3599, NGC 3489, and NGC 4274, which exhibits a prominent ghost halo for two of them. This method successfully removes this.

Converging on the Initial Mass Function of Stars

Understanding the origin of stellar masses—the initial mass function (IMF)— remains one of the most challenging problems in astrophysics. The IMF is a key ingredient for simulations of galaxy formation and evolution, and is used to calibrate star formation relations in extra-galactic observations. Modeling the IMF directly in hydrodynamical simulations has been attempted in several previous studies, but the most important processes that control the IMF remain poorly understood. This is because predicting the IMF from direct hydrodynamical simulations involves complex physics such as turbulence, magnetic fields, radiation feedback and mechanical feedback, all of which are difficult to model and the methods used have limitations in terms of accuracy and computational efficiency. Moreover, a physical interpretation of the simulated IMFs requires a numerically converged solution at high resolution, which has so far not been convincingly demonstrated. Here we present a resolution study of star cluster formation aimed at producing a converged IMF. We compare a set of magnetohydrodynamical (MHD) adaptive-mesh-refinement simulations with three different implementations of the thermodynamics of the gas: 1) with an isothermal equation of state (EOS), 2) with a polytropic EOS, and 3) with a simple stellar heating feedback model. We show that in the simulations with an isothermal or polytropic EOS, the number of stars and their mass distributions depend on the numerical resolution. By contrast, the simulations that employ the simple radiative feedback module demonstrate convergence in the number of stars formed and in their IMFs.

Delay-time distribution of core-collapse supernovae with late events resulting from binary interaction

Most massive stars, the progenitors of core-collapse supernovae, are in close binary systems and may interact with their companion through mass transfer or merging. We undertake a population synthesis study to compute the delay-time distribution of core-collapse supernovae, that is, the supernova rate versus time following a starburst, taking into account binary interactions. We test the systematic robustness of our results by running various simulations to account for the uncertainties in our standard assumptions. We find that a significant fraction, $15^{+9}_{-8}$%, of core-collapse supernovae are `late', that is, they occur 50-200 Myrs after birth, when all massive single stars have already exploded. These late events originate predominantly from binary systems with at least one, or, in most cases, with both stars initially being of intermediate mass ($4-8M_{\odot}$). The main evolutionary channels that contribute often involve either the merging of the initially more massive primary star with its companion or the engulfment of the remaining core of the primary by the expanding secondary that has accreted mass at an earlier evolutionary stage. Also, the total number of core-collapse supernovae increases by $14^{+15}_{-14}$% because of binarity for the same initial stellar mass. The high rate implies that we should have already observed such late core-collapse supernovae, but have not recognized them as such. We argue that $\phi$ Persei is a likely progenitor and that eccentric neutron star - white dwarf systems are likely descendants. Late events can help explain the discrepancy in the delay-time distributions derived from supernova remnants in the Magellanic Clouds and extragalactic type Ia events, lowering the contribution of prompt Ia events. We discuss ways to test these predictions and speculate on the implications for supernova feedback in simulations of galaxy evolution.