Post-starburst Galaxies Research Articles

DESI Massive Poststarburst Galaxies at z ∼ 1.2 Have Compact Structures and Dense Cores

Poststarburst galaxies (PSBs) are young quiescent galaxies that have recently experienced a rapid decrease in star formation, allowing us to probe the fast-quenching period of galaxy evolution. In this work, we obtained Hubble Space Telescope (HST)/WFC3 F110W imaging to measure the sizes of 171 massive ( log(M*/M⊙)∼11) spectroscopically identified PSBs at 1 < z 1.3 selected from the DESI Survey Validation luminous red galaxy sample. This statistical sample constitutes an order of magnitude increase from the ∼20 PSBs with space-based imaging and deep spectroscopy. We perform structural fitting of the target galaxies with pysersic and compare them to quiescent and star-forming galaxies in the 3D-HST survey. We find that these PSBs are more compact than the general population of quiescent galaxies, lying systematically ∼0.1 dex below the established size–mass relation. However, their central surface mass densities are similar to those of their quiescent counterparts ( log(Σ1kpc/(M⊙kpc−2))∼10.1 ). These findings are easily reconciled by later ex situ growth via minor mergers or a slight progenitor bias. These PSBs are round in projection (b/a median ∼ 0.8), suggesting that they are primarily spheroids, not disks, in 3D. We find no correlation between the time since quenching and light-weighted PSB sizes or central densities. This disfavors apparent structural growth due to the fading of centralized starbursts in this galaxy population. Instead, we posit that the fast quenching of massive galaxies at this epoch occurs preferentially in galaxies with preexisting compact structures.

Pulling Back the Curtain on Shocks and Star Formation in NGC 1266 with Gemini-NIFS

We present Gemini near-infrared integral field spectrograph K-band observations of the central 400 pc of NGC 1266, a nearby (D ≈ 30 Mpc) post-starburst galaxy with a powerful multiphase outflow and a shocked interstellar medium. We detect seven H2 rovibrational emission lines excited thermally to T ∼ 2000 K, and weak Brγ emission, consistent with a fast continuous shock (or C-shock). With these bright H2 lines, we observe the spatial structure of the shock with an unambiguous tracer for the first time. The Brγ emission is concentrated in the central ≲100 pc, indicating that any remaining star formation in NGC 1266 is in the nucleus, while the surrounding cold molecular gas has little ongoing star formation. Though it is unclear what fraction of this Brγ emission is from star formation or the active galactic nuclei (AGN), assuming it is entirely due to star formation we measure an instantaneous star formation rate of 0.7 M ⊙ yr−1, though the star formation rate may be significantly higher in the presence of additional extinction. NGC 1266 provides a unique laboratory to study the complex interactions between AGN, outflows, shocks, and star formation, all of which are necessary to unravel the evolution of the post-starburst phase.

To High Redshift and Low Mass: Exploring the Emergence of Quenched Galaxies and Their Environments at 3 < z < 6 in the Ultra-deep JADES MIRI F770W Parallel

We present the robust selection of high-redshift quiescent galaxies (QG) and poststarburst (PSB) galaxies using ultra-deep NIRCam and MIRI imaging from the JWST Advanced Deep Extragalactic Survey (JADES). At 3 < z < 6, MIRI 7.7 μm imaging provides rest-frame J band, which is commonly used to break the degeneracy between old stellar populations and dust attenuation at lower redshifts. We identify 23 passively evolving galaxies in UVJ color space in a mass-limited (log M ⋆/M ⊙ ≥ 8.5) sample over 8.8 arcmin2. An evaluation of the contribution of the 7.7 μm shows that JADES-like NIRCam coverage (9+ photometric bands) can compensate for lacking the J band at these redshifts; however, more limited three-band selections perform better with MIRI. Our sample is characterized by rapid quenching timescales (∼100–600 Myr) with formation redshifts z f ≲ 9 and includes a potential record-holding massive QG at zphot=5.33−0.17+0.16 and two QGs with evidence for significant residual dust content (A V ∼ 1–2). In addition, we present a large sample of 12 log M ⋆/M ⊙ = 8.5–9.5 PSBs, demonstrating that UVJ selection can be extended to low mass. An analysis of the environment of our sample reveals that the group known as the Cosmic Rose contains a massive QG and a dust-obscured star-forming galaxy (a so-called Jekyll and Hyde pair) plus three additional QGs within ∼20 kpc. Moreover, the Cosmic Rose is part of a larger overdensity at z ∼ 3.7, which contains 7/12 of our low-mass PSBs. Another four low-mass PSBs are members of an overdensity at z ∼ 3.4; this result strongly indicates low-mass PSBs are preferentially associated with overdense environments at z > 3.

A fast-rotator post-starburst galaxy quenched by supermassive black-hole feedback at z = 3

AbstractThe most massive galaxies in the Universe stopped forming stars due to the time-integrated feedback from central supermassive black holes (SMBHs). However, the exact quenching mechanism is not yet understood, because local massive galaxies were quenched billions of years ago. Here we present JWST/NIRSpec integral-field spectroscopy observations of GS-10578, a massive, quiescent galaxy at redshift z = 3.064 ± 0.002. From its spectrum, we measure a stellar mass M⋆ = 1.6 ± 0.2 × 1011 M⊙ and a dynamical mass Mdyn = 2.0 ± 0.5 × 1011 M⊙. Half of its stellar mass formed at z = 3.7–4.6, and the system is now quiescent, with a current star-formation rate of less than 19 M⊙ yr−1. We detect ionized- and neutral-gas outflows traced by [O iii] emission and Na i absorption, with mass outflow rates 0.14–2.9 and 30–100 M⊙ yr−1, respectively. Outflow velocities reach vout ≈ 1,000 km s−1, comparable to the galaxy escape velocity. GS-10578 hosts an active galactic nucleus, evidence that these outflows are due to SMBH feedback. The neutral outflow rate is higher than the star-formation rate. Hence, this is direct evidence for ejective SMBH feedback, with a mass loading capable of interrupting star formation by rapidly removing its fuel. Stellar kinematics show ordered rotation, with spin parameter $${\lambda }_{{{{R}}}_{{\rm{e}}}}=0.62\pm 0.07$$ λ R e = 0.62 ± 0.07 , meaning GS-10578 is rotation-supported. This study presents direct evidence for ejective active galactic nucleus feedback in a massive, recently quenched galaxy, thus helping to clarify how SMBHs quench their hosts. The high value of $${\lambda }_{{{{R}}}_{{\rm{e}}}}$$ λ R e implies that quenching can occur without destroying the stellar disk.

Efficient NIRCam Selection of Quiescent Galaxies at 3 < z < 6 in CEERS

Substantial populations of massive quiescent galaxies at z ≥ 3 challenge our understanding of rapid galaxy growth and quenching over short timescales. In order to piece together this evolutionary puzzle, more statistical samples of these objects are required. Established techniques for identifying massive quiescent galaxies are increasingly inefficient and unconstrained at z > 3. As a result, studies report that as much as 70% of quiescent galaxies at z > 3 may be missed from existing surveys. In this work, we propose a new empirical color selection technique designed to select massive quiescent galaxies at 3 ≲ z ≲ 6 using JWST NIRCam imaging data. We use empirically constrained galaxy spectral energy distribution (SED) templates to define a region in the F277W − F444W versus F150W − F277W color plane that captures quiescent galaxies at z > 3. We apply these color selection criteria to the Cosmic Evolution Early Release Science (CEERS) Survey and use SED fitting on sources in the region to identify 44 candidate z ≳ 3 quiescent galaxies. Over half of these sources are newly discovered and, on average, exhibit specific star formation rates of poststarburst galaxies. Most of these sources would not be discovered using canonical UVJ diagrams. We derive volume density estimates of n ∼ 1–4 × 10−5 Mpc−3 at 3 < z < 5, finding excellent agreement with existing reports on similar populations in the CEERS field. Thanks to NIRCam’s wavelength coverage and sensitivity, this technique provides an efficient tool to search for large samples of these rare galaxies.

Extended Emission-line Regions in Poststarburst Galaxies Hosting Tidal Disruption Events

We report the discovery of an extended emission-line region (EELR) in MUSE observations of Markarian 950, a nearby (z = 0.01628) poststarburst (PSB) galaxy that hosted the tidal disruption event (TDE) iPTF 16fnl. The EELR requires a nonstellar ionizing continuum with a luminosity Lion,min≳1043 erg s−1, inconsistent with the current weak state (L IR,AGN < 2.5 × 1042 erg s−1) of the galactic nucleus. The ionized gas has low velocity (∼–50 km s−1) and low turbulence (σ gas ≲ 50 km s−1) and is kinematically decoupled from the stellar motions, indicating that the gas kinematics is not active galactic nucleus (AGN) driven. Markarian 950 is the third PSB galaxy to host a weak nuclear ionizing source as well as an EELR and a TDE. The overall properties of these three galaxies, including the kinematics and accretion history, are unusual but strikingly similar. We estimate that the incidence of EELRs in PSB-TDE hosts is a factor of ∼10 × higher than in other PSB galaxies. This suggests that a gas-rich postmerger environment is a key ingredient in driving elevated TDE rates. Based on the current observations, we cannot rule out that the EELRs may be powered through an elevated TDE rate in these galaxies. If the EELRs are not TDE powered, the presence of intermittent AGN activity, and in particular the fading of the AGN, may be associated with an increased TDE rate and/or an increased rate of detecting TDEs.

Not So Windy After All: MUSE Disentangles AGN-driven Winds from Merger-induced Flows in Galaxies along the Starburst Sequence

Poststarburst galaxies are believed to be in a rapid transition between major merger starbursts and quiescent ellipticals, where active galactic nucleus (AGN) feedback is suggested as one of the processes responsible for the quenching. To study the role of AGN feedback, we constructed a sample of poststarburst candidates with AGN and indications of ionized outflows in optical. We use MUSE/VLT observations to spatially resolve the properties of the stars and multiphase gas in five of them. All galaxies show signatures of interaction/merger in their stellar or gas properties, with some at an early stage of interaction with companions ∼50 kpc, suggesting that optical poststarburst signatures may be present well before the final starburst and coalescence. We detect narrow and broad kinematic components in multiple transitions in all the galaxies. Our detailed analysis of their kinematics and morphology suggests that, contrary to our expectation, the properties of the broad kinematic components are inconsistent with AGN-driven winds in three out of five galaxies. The two exceptions are also the only galaxies in which spatially resolved NaID P-Cygni profiles are detected. In some cases, the observations are more consistent with interaction-induced galactic-scale flows, an often overlooked process. These observations raise the question of how to interpret broad kinematic components in interacting and perhaps also in active galaxies, in particular when spatially resolved observations are not available or cannot rule out merger-induced galactic-scale motions. We suggest that NaID P-Cygni profiles are more effective outflow tracers, and use them to estimate the energy that is carried by the outflow.

Correction to: Chemical evolution of local post-starburst galaxies: Implications for the mass-metallicity relation

Correction to: Chemical evolution of local post-starburst galaxies: Implications for the mass-metallicity relation

Can Supernovae Quench Star Formation in High-z Galaxies?

JWST is providing a unique opportunity to directly study the feedback processes regulating star formation (SF) in early galaxies. The two z > 5 quiescent systems (JADES-GS-z7-01-QU and MACS0417-z5BBG) detected so far show a recent starburst after which SF is suppressed. To clarify whether such quenching is due to supernova (SN) feedback, we have developed a minimal physical model. We derive a condition on the minimum star formation rate, SFRmin , lasting for a time interval Δt b , required to quench SF in a galaxy at redshift z, with gas metallicity Z, and hosted by a halo of mass M h . We find that lower (z, Z, M h ) systems are more easily quenched. We then apply the condition to JADES-GS-z7-01-QU (z = 7.3, M ⋆ = 108.6 M ⊙) and MACS0417-z5BBG (z = 5.2, M ⋆ = 107.6 M ⊙) and find that SN feedback largely fails to reproduce the observed quenched SF history. Alternatively, we suggest that SF is rapidly suppressed by radiation-driven dusty outflows sustained by the high specific star formation rates (43 and 25 Gyr−1, respectively) of the two galaxies. Our model provides a simple tool to interpret the SF histories of post-starburst galaxies and unravel quenching mechanisms from incoming JWST data.

Evolution of gas flows along the starburst to post-starburst to quiescent galaxy sequence

ABSTRACT We measure velocity offsets in the $\rm{Na {}\rm {\small I}}$ λλ5890, 5896 ($\rm{Na {}\rm {\small D}}$) interstellar medium absorption lines to track how neutral galactic winds change as their host galaxies evolve. Our sample of ∼80 000 SDSS spectra at 0.010 &amp;lt; z &amp;lt; 0.325 includes starburst, post-starburst, and quiescent galaxies, forming an evolutionary sequence of declining star formation rate (SFR). We detect bulk flows across this sequence, mostly at higher host stellar masses (log(M*/M⊙)) &amp;gt; 10). Along this sequence, the fraction of outflows decreases (76 ± 2 per cent to 65 ± 4 per cent to a 3σ upper limit of 34 per cent), and the mean velocity offset changes from outflowing to inflowing (−84.6 ± 5.9 to −71.6 ± 11.4 to $76.6\pm 2.3\, \rm km s^{-1}$). Even within the post-starburst sample, wind speed decreases with time elapsed since the starburst ended. These results reveal that outflows diminish as galaxies age. For post-starbursts, there is evidence for an AGN contribution, especially to the speediest outflows: (1) SFR declines faster in time than outflow velocity, a decoupling arguing against massive stellar feedback; (2) of the few outflows strong enough to escape the interstellar medium (9/105), three of the four hosts with measured emission lines are Seyfert galaxies. For discy starburst galaxies, however, the trends suggest flows out of the stellar disc plane (with outflow 1/2-opening angle &amp;gt; 45°) instead of from the nucleus: the wind velocity decreases as the disc becomes more edge-on, and the outflow fraction, constant at $\sim 90~{{\ \rm per\ cent}}$ for disc inclinations i &amp;lt; 45°, steadily decreases from $\sim 90~{{\ \rm per\ cent}}$ to 20 per cent for i &amp;gt; 45°.

Post-starburst Galaxies in SDSS-IV MaNGA: Two Broad Categories of Evolutionary Pathways

We study the size–mass relation (SMR) and recent star formation history (SFH) of post-starburst (PSB) galaxies in the local Universe using spatially resolved spectroscopy from the final data release of MaNGA. Our sample includes 489 PSB galaxies: 94 cPSB galaxies with central PSB regions, 85 rPSB galaxies with ringlike PSB regions, and 310 iPSB galaxies with irregular PSB regions. When compared to control galaxies of similar star formation rate, redshift, and mass, a similar SMR is found for all types of PSB samples except the cPSB galaxies, which have smaller sizes at intermediate masses ( 9.5≲log10(M*/M⊙)≲10.5 ). The iPSB galaxies in the star-forming sequence (iPSB-SF) show no/weak gradients in D n (4000), EW(Hδ A ), and EW(Hα), consistent with the global star-forming status of this type of galaxy, while the quiescent iPSB (iPSB-Q) sample shows negative gradients in D n (4000) and positive gradients in EW(Hδ A ), indicating older stellar populations in the inner regions. Both the cPSB and rPSB samples show positive gradients in D n (4000) and negative gradients in EW(Hδ A ), indicating younger stellar populations in the inner regions. These results imply that the four types of PSB galaxies can be broadly divided into two distinct categories in terms of evolutionary pathway: (1) iPSB-SF and iPSB-Q, which have SMRs and SFHs similar to control galaxies, preferring an inside-out quenching process, and (2) rPSB and cPSB, which appear to be different stages of the same event and likely to follow the outside-in quenching process driven by disruption events such as mergers that result in a more compact structure as quenching proceeds.

Chemical evolution of local post-starburst galaxies: implications for the mass–metallicity relation

ABSTRACT We use the stellar fossil record to constrain the stellar metallicity evolution and star-formation histories of the post-starburst (PSB) regions within 45 local PSB galaxies from the MaNGA survey. The direct measurement of the regions’ stellar metallicity evolution is achieved by a new two-step metallicity model that allows for stellar metallicity to change at the peak of the starburst. We also employ a Gaussian process noise model that accounts for correlated errors introduced by the observational data reduction or inaccuracies in the models. We find that a majority of PSB regions (69 per cent at &amp;gt;1σ significance) increased in stellar metallicity during the recent starburst, with an average increase of 0.8 dex and a standard deviation of 0.4 dex. A much smaller fraction of PSBs are found to have remained constant (22 per cent) or declined in metallicity (9 per cent, average decrease 0.4 dex, standard deviation 0.3 dex). The pre-burst metallicities of the PSB galaxies are in good agreement with the mass–metallicity (MZ) relation of local star-forming galaxies. These results are consistent with hydrodynamic simulations, which suggest that mergers between gas-rich galaxies are the primary formation mechanism of local PSBs, and rapid metal recycling during the starburst outweighs the impact of dilution by any gas inflows. The final mass-weighted metallicities of the PSB galaxies are consistent with the MZ relation of local passive galaxies. Our results suggest that rapid quenching following a merger-driven starburst is entirely consistent with the observed gap between the stellar mass–metallicity relations of local star-forming and passive galaxies.

An Explanation for the Overrepresentation of Tidal Disruption Events in Post-starburst Galaxies

Tidal disruption events (TDEs) provide a valuable probe for studying the dynamics of stars in the nuclear environments of galaxies. Recent observations show that TDEs are strongly overrepresented in post-starburst or “green valley” galaxies, although the underlying physical mechanism remains unclear. Considering the possible interaction between stars and active galactic nucleus (AGN) disks, the TDE rates can be greatly changed compared to those in quiescent galactic nuclei. In this work, we revisit TDE rates by incorporating an evolving AGN disk within the framework of “loss cone” theory. We numerically evolve the Fokker–Planck equations by considering star–disk interactions, in situ star formation in the unstable region of the outer AGN disk, and evolution of the accretion process of supermassive black holes. We find that the TDE rates are enhanced by about 2 orders of magnitude shortly after the AGN transitions into an inactive stage. During this phase, the accumulated stars rapidly scatter into the loss cone due to the disappearance of the inner standard thin disk. Our results provide an explanation for the overrepresentation of TDEs in post-starburst galaxies.

A Multiwavelength View of Black Holes and Outflows in Post-starburst Galaxies

Abstract. Post-starburst galaxies (PSBs) have quenched (significant decline in star formation rate) both recently and rapidly (≲Gyr). They are thus promising in providing insights into activities that are happening at the early stage of quenching. While studies have suggested that black hole feedback in the form of active galactic nuclei (AGN) and outflows play important roles in quenching, the details of how they impact the host galaxies and their interplay with other quenching mechanisms are still not fully understood. We find that PSBs commonly show signatures of AGN activity but they appear to be weak and/or heavily obscured. These AGN might be able to drive outflows but they are likely not strong enough to remove gas from the host galaxy. Direct evidence of AGN quenching the star formation of the host galaxy is still missing and AGN likely quench by disturbing rather than expelling the gas.

Needle in a Haystack: Finding Supermassive Black Hole-related Flares in the Zwicky Transient Facility Public Survey

Transient accretion events onto supermassive black holes (SMBHs), such as tidal disruption events (TDEs), Bowen Fluorescence Flares (BFFs), and active galactic nuclei (AGNs), which are accompanied by sudden increases of activity, offer a new window onto the SMBH population, accretion physics, and stellar dynamics in galaxy centers. However, such transients are rare and finding them in wide-field transient surveys is challenging. Here we present the results of a systematic real-time search for SMBH-related transients in Zwicky Transient Facility (ZTF) public alerts, using various search queries. We examined 345 rising events coincident with a galaxy nucleus, with no history of previous activity, of which 223 were spectroscopically classified. Of those, five (2.2%) were TDEs, one (0.5%) was a BFF, and two (0.9%) were AGN flares. Limiting the search to blue events, the fraction of TDEs nearly doubles to 4.1%, and no TDEs are missed. Limiting the search further to candidate post-starburst galaxies increases the relative number of TDEs to 16.7%, but the absolute numbers in such a search are small. The main contamination source is supernovae (95.1% of classified events), of which the majority (82.2% of supernovae) are of Type Ia. In a comparison set of 39 events with limited photometric history, the AGN contamination increases to ∼30%. Host galaxy offset is not a significant discriminant of TDEs in current ZTF data, but might be useful in higher-resolution data. Our results can be used to quantify the efficiency of various SMBH-related transient search strategies in optical surveys such as ZTF and the Legacy Survey of Space and Time.

Stars, Gas, and Star Formation of Distant Post-starburst Galaxies

We present a comprehensive multiwavelength study of five post-starburst galaxies with M * > 1011 M ⊙ at z ∼ 0.7, examining their stars, gas, and current and past star formation activities. Using optical images from the Subaru Telescope and Hubble Space Telescope, we observe a high incidence of companion galaxies and low surface brightness tidal features, indicating that quenching is closely related to interactions between galaxies. From optical spectra provided by the LEGA-C survey, we model the stellar continuum to derive the star formation histories and show the stellar masses of progenitors ranging from 2 × 109 M ⊙ to 1011 M ⊙, undergoing a burst of star formation several hundred million years prior to observation, with a decay timescale of ∼100 Myr. Our Atacama Large Millimeter/submillimeter Array observations detect CO(2-1) emission in four galaxies, with the molecular gas spreading over up to >1″, or ∼10 kpc, with a mass of up to ∼2 × 1010 M ⊙. However, star-forming regions are unresolved by either the slit spectra or 3 GHz continuum observed by the Very Large Array. Comparisons between the star formation rates (SFRs) and gas masses, and the sizes of CO emission and star-forming regions suggest a low star-forming efficiency. We show that the SFRs derived from IR and radio luminosities with commonly used calibrations tend to overestimate the true values because of the prodigious amount of radiation from old stars and the contribution from active galactic nuclei (AGNs), as the optical spectra reveal weak AGN-driven outflows.

Modeling Star Formation Histories of Post Starburst Galaxies with BAGPIPES

Post-starburst galaxies provide insights into quenching mechanisms. While these galaxies show signs of quiescence, obtaining accurate star formation rates (SFRs) is difficult. This research uses BAGPIPES, a Bayesian spectrum-fitting code, to obtain model spectra of 13 galaxies. From these models, we calculate SFRs over 10 and 100 Myr durations. We compare these SFRs to published Hα and total infrared (TIR) SFRs to examine biases in these measurements. We find that the BAGPIPES-inferred 10, 100 Myr SFRs are comparable with Hα and TIR traced SFRs, respectively, with longer-duration tracers showing systematically higher SFRs.

LOFAR Detection of Post-starburst Galaxies

We report detection of 104 radio counterparts in a sample of 1204 Sloan Digital Sky Survey-detected post-starburst (PSB) galaxies using the 144 MHz LOw-Frequency ARray (LOFAR) Two-meter Sky Survey Data Release 2 (LoTSS2). The detection rate of optically selected PSB galaxies in the 144 MHz images is ∼8.6% while that in the 1.4 GHz was estimated to be ∼1.8% using the Very Large Array-Faint Images of the Radio Sky at Twenty centimeter images. Although the 144 MHz radio luminosity shows no significant correlation with luminosity at other wave bands from optical to far-infrared, the fraction of radio-detected PSB galaxies shows a clear trend with absolute magnitude and colors in optical and Wide-field Infrared Survey Explorer infrared bands. A preliminary analysis suggests that the radio emission is dominated by Active Galactic Nuclei.

The Correlations of Stellar Tidal Disruption Rates with Properties of Massive Black Holes and Their Host Galaxies

Stars can be either disrupted as tidal disruption events (TDEs) or swallowed whole by massive black holes (MBHs) at galactic centers when they approach sufficiently close to these MBHs. In this work, we investigate the correlations of such stellar consumption rates with both the MBH mass M BH and the inner slope of the host-galaxy mass density distribution α. We introduce a simplified analytical power-law model with a power-law stellar-mass density distribution surrounding MBHs and separate the contributions of two-body relaxation and stellar orbital precession for the stellar orbital angular momentum evolution in nonspherical galaxy potentials. The stellar consumption rates derived from this simplified model can be well consistent with the numerical results obtained with a more realistic treatment of stellar distributions and dynamics around MBHs, providing an efficient way to estimate TDE rates. The origin of the correlations of stellar consumption rates with M BH and α is explained by the dependence of this analytical model on those MBH/host-galaxy properties and by the separation of the stellar angular momentum evolution mechanisms. We propose that the strong positive correlation between the rates of stellar consumption due to two-body relaxation and α provides one interpretation for the overrepresentation of TDEs found in some rare E+A/poststarburst galaxies. We find high TDE rates for giant stars, up to those for solar-type stars. Understanding the origin of the correlations of the stellar consumption rates will be necessary for obtaining the demographics of MBHs and their host galaxies via TDEs.

Star formation and molecular gas properties of post-starburst galaxies

ABSTRACT Post-starburst galaxies are believed to be in a rapid transition between major merger starbursts and quiescent ellipticals. Their optical spectrum is dominated by A-type stars, suggesting a starburst that was quenched recently. While optical observations suggest little ongoing star formation, some have been shown to host significant molecular gas reservoirs. This led to the suggestion that gas depletion is not required to end the starburst, and that star formation is suppressed by other processes. We present NOEMA CO(1−0) observations of 15 post-starburst galaxies with emission lines consistent with active galactic nucleus (AGN) photoionization. We collect post-starburst candidates with molecular gas measurements from the literature, with some classified as classical E + A, while others with line ratios consistent with AGN and/or shock ionization. Using far-infrared observations, we show that systems that were reported to host exceptionally large molecular gas reservoirs host in fact obscured star formation, with some systems showing star formation rates comparable to ULIRGs. Among E + A galaxies with molecular gas measurements, 7 out of 26 (26 per cent) host obscured starbursts. Using far-infrared observations, post-starburst candidates show similar SFR–$M_{\mathrm{H_2}}$ and Kennicutt–Schmidt relations to those observed in star-forming and starburst galaxies. In particular, there is no need to hypothesize star formation quenching by processes other than the consumption of molecular gas by star formation. The combination of optical, far-infrared, and CO observations indicates that some regions within these galaxies have been recently quenched, while others are still forming stars in highly obscured regions. All this calls into question the traditional interpretation of such galaxies.