Neutral Gas-phase Research Articles

Multiphase Black Hole Feedback and a Bright [C ii] Halo in a LoBAL Quasar at z ∼ 6.6

Although the mass growth of supermassive black holes during the epoch of reionization is expected to play a role in shaping the concurrent growth of their host galaxies, observational evidence of feedback at z ≳ 6 is still sparse. We perform the first multiscale and multiphase characterization of black-hole-driven outflows in the z ∼ 6.6 quasar J0923+0402 and assess how these winds impact the cold gas reservoir. We employ the SimBAL spectral synthesis to fit broad absorption line features and find a powerful ionized outflow on a ≲210 pc scale, with a kinetic power ∼2%–100% of the quasar luminosity. Atacama Large Millimeter/submillimeter Array (ALMA) observations of [C ii] emission allow us to study the morphology and kinematics of the cold gas. We detect high-velocity [C ii] emission, likely associated with a cold neutral outflow at ∼0.5−2 kpc scale in the host galaxy, and a bright extended [C ii] halo with a size of ∼15 kpc. For the first time at such an early epoch, we accurately constrain the outflow energetics in both the ionized and the atomic neutral gas phases. We find such energetics to be consistent with expectations for an efficient feedback mechanism, and both ejective and preventative feedback modes are likely at play. The scales and energetics of the ionized and atomic outflows suggest that they might be associated with different quasar accretion episodes. The results of this work indicate that strong black hole feedback is occurring in quasars at z ≳ 6 and is likely responsible for shaping the properties of the cold gas reservoir up to circumgalactic scales.

Extragalactic Magnetism with SOFIA (SALSA Legacy Program). VII. A Tomographic View of Far-infrared and Radio Polarimetric Observations through MHD Simulations of Galaxies

The structure of magnetic fields in galaxies remains poorly constrained, despite the importance of magnetism in the evolution of galaxies. Radio synchrotron and far-infrared (FIR) polarization and polarimetric observations are the best methods to measure galactic scale properties of magnetic fields in galaxies beyond the Milky Way. We use synthetic polarimetric observations of a simulated galaxy to identify and quantify the regions, scales, and interstellar medium (ISM) phases probed at FIR and radio wavelengths. Our studied suite of magnetohydrodynamical cosmological zoom-in simulations features high-resolutions (10 pc full-cell size) and multiple magnetization models. Our synthetic observations have a striking resemblance to those of observed galaxies. We find that the total and polarized radio emission extends to approximately double the altitude above the galactic disk (half-intensity disk thickness of h I radio ∼ h PI radio = 0.23 ± 0.03 kpc) relative to the total FIR and polarized emission that are concentrated in the disk midplane (h I FIR ∼ h PI FIR = 0.11 ± 0.01 kpc). Radio emission traces magnetic fields at scales of ≳300 pc, whereas FIR emission probes magnetic fields at the smallest scales of our simulations. These scales are comparable to our spatial resolution and well below the spatial resolution (<300 pc) of existing FIR polarimetric measurements. Finally, we confirm that synchrotron emission traces a combination of the warm neutral and cold neutral gas phases, whereas FIR emission follows the densest gas in the cold neutral phase in the simulation. These results are independent of the ISM magnetic field strength. The complementarity we measure between radio and FIR wavelengths motivates future multiwavelength polarimetric observations to advance our knowledge of extragalactic magnetism.

The impact of cosmic rays on the interstellar medium and galactic outflows of Milky Way analogues

Abstract During the last decade, cosmological simulations have managed to reproduce realistic and morphologically diverse galaxies, spanning the Hubble sequence. Central to this success was a phenomenological calibration of the few included feedback processes, whilst glossing over higher complexity baryonic physics. This approach diminishes the predictive power of such simulations, preventing to further our understanding of galaxy formation. To tackle this fundamental issue, we investigate the impact of cosmic rays (CRs) and magnetic fields on the interstellar medium (ISM) and the launching of outflows in a cosmological zoom-in simulation of a Milky Way-like galaxy. We find that including CRs decreases the stellar mass of the galaxy by a factor of 10 at high redshift and ∼4 at cosmic noon, leading to a stellar mass to halo mass ratio in good agreement with abundance matching models. Such decrease is caused by two effects: i) a reduction of cold, high-density, star-forming gas, and ii) a larger fraction of SN events exploding at lower densities, where they have a higher impact. SN-injected CRs produce enhanced, multi-phase galactic outflows, which are accelerated by CR pressure gradients in the circumgalactic medium of the galaxy. While the mass budget of these outflows is dominated by the warm ionised gas, warm neutral and cold gas phases contribute significantly at high redshifts. Importantly, our work shows that future JWST observations of galaxies and their multi-phase outflows across cosmic time have the ability to constrain the role of CRs in regulating star formation.

ALMA reveals a compact and massive molecular outflow driven by the young AGN in a nearby ULIRG

ABSTRACT The ultraluminous infrared galaxy F13451+1232 is an excellent example of a galaxy merger in the early stages of active galactic nucleus (AGN) activity, a phase in which AGN-driven outflows are expected to be particularly important. However, previous observations have determined that the mass outflow rates of the warm ionized and neutral gas phases in F13451+1232 are relatively modest, and there has been no robust detection of molecular outflows. Using high-spatial resolution Atacama Large Millimeter/submillimeter Array CO(1–0) observations, we detect a kiloparsec-scale circumnuclear disc, as well as extended (r ∼ 440 pc), intermediate-velocity (300 &amp;lt; |v| &amp;lt; 400 km s−1) cold molecular gas emission that cannot be explained by rotational disc motions. If interpreted as AGN-driven outflows, the mass outflow rates associated with this intermediate-velocity gas are relatively modest ($\dot{M}_\mathrm{out}=22$–27 M⊙ yr−1); however, we also detect a compact (rout &amp;lt; 120 pc), high-velocity (400 &amp;lt; v &amp;lt; 680 km s−1) cold molecular outflow near the primary nucleus of F13451+1232, which carries an order of magnitude more mass ($\dot{M}_\mathrm{out}$ ∼ 230 M⊙ yr−1) than (and several times the kinetic power of) the previously detected warmer phases. Moreover, the similar spatial scales of this compact outflow and the radio structure indicate that it is likely accelerated by the small-scale (r ∼ 130 pc) AGN jet in the primary nucleus of F13451+1232. Considering the compactness of the nuclear outflow and intermediate-velocity non-rotating gas that we detect, we argue that high-spatial resolution observations are necessary to properly quantify the properties of AGN-driven outflows and their impacts on host galaxies.

Infrared bands of neutral gas-phase carbon clusters in a broad spectral range.

The identification of species in the interstellar medium requires precise and molecule-specific spectroscopic information in the laboratory framework, in broad spectral ranges and under conditions relevant to interstellar environments. In this work, we measure the gas-phase infrared spectra of neutral carbon clusters, CN (N = 6-11), in a molecular beam. The CN distribution is formed by photofragmentation of C60 molecules, concurrently showing a top-down formation mechanism. A broad spectral range in the infrared between 500-3200 cm-1 (20-3.125 μm) is investigated. We observe strong bands between 5 and 6 μm, in conjunction with novel features in the 3 μm region. Density functional theory calculations reveal that these short wavelength modes correspond to combination bands with significant infrared intensity. Moreover, we identify the N ≤ 10 clusters as linear, while C11 adopts a ring configuration, placing the linear-to-ring transition at N = 11 under our molecular beam conditions. The linearity of C10 is discussed based on the formation pathway from larger clusters in energetic conditions. Given the vast and very precise infrared information already been released from the James Webb Space Telescope mission, this infrared spectroscopic data set in conjunction with information on formation mechanisms is of major relevance for identifying neutral carbon clusters in astronomical environments.

Detection and prebiotic chemistry of possible glycine precursor molecule methylenimine towards the hot molecular core G10.47+0.03

Abstract Amino acids are essential for the synthesis of protein. Amino acids contain both amine (R–NH2) and carboxylic acid (R–COOH) functional groups, which help to understand the possible formation mechanism of life in the universe. Among the 20 types of amino acids, glycine (NH2CH2COOH) is known as the simplest non-essential amino acid. In the last 40 years, all surveys of NH2CH2COOH in the interstellar medium, especially in the star-formation regions, have failed at the millimetre and sub-millimetre wavelengths. We aimed to identify the possible precursors of NH2CH2COOH, because it is highly challenging to identify NH2CH2COOH in the interstellar medium. Many laboratory experiments have suggested that methylenimine (CH2NH) plays a key role as a possible precursor of NH2CH2COOH in the star-formation regions via the Strecker synthesis reaction. After spectral analysis using the local thermodynamic equilibrium (LTE) model, we successfully identified the rotational emission lines of CH2NH towards the hot molecular core G10.47+0.03 using the Atacama Compact Array (ACA). The estimated column density of CH2NH towards G10.47+0.03 is (3.40 ± 0.2) × 1015 cm−2 with a rotational temperature of 218.70 ± 20 K, which is estimated from the rotational diagram. The fractional abundance of CH2NH with respect to H2 towards G10.47+0.03 is 2.61 × 10−8. We found that the derived abundance of CH2NH agree fairly well with the existing two-phase warm-up chemical modelling abundance value of CH2NH. We discuss the possible formation pathways of CH2NH within the context of hot molecular cores, and we find that CH2NH is likely mainly formed via neutral–neutral gas-phase reactions of CH3 and NH radicals towards G10.47+0.03.

Computational approaches to identify sites of phosphorylation.

Due to their oftentimes ambiguous nature, phosphopeptide positional isomers can present challenges in bottom-up mass spectrometry-based workflows as search engine scores alone are often not enough to confidently distinguish them. Additional scoring algorithms can remedy this by providing confidence metrics in addition to these search results, reducing ambiguity. Here we describe challenges to interpreting phosphoproteomics data and review several different approaches to determine sites of phosphorylation for both data-dependent and data-independent acquisition-based workflows. Finally, we discuss open questions regarding neutral losses, gas-phase rearrangement, and false localization rate estimation experienced by both types of acquisition workflows and best practices for managing ambiguity in phosphosite determination.

Gas-Phase Characterization of Hypervalent Carbon Compounds Bearing 7-6-7-Ring Skeleton: Penta- versus Tetra-Coordinate Isomers.

In this study, we afford explicit characterizations of the electronic and geometrical structures of recently reported hypervalent penta-coordinate carbon compounds by using gas-phase characterization techniques: photodissociation spectroscopy (PDS) and ion mobility-mass spectrometry (IM-MS). In particular for a compound with moderately electron-donating ligands, bearing p-methylthiophenyl substituents, the coexistence of tetra- and penta-coordinate isomers is confirmed, consistent with solution characterizations. It is in sharp contrast to the exclusive tetra-coordinate form (with normal valence of the central carbon atom) in the single crystal. This suggests that a non-polar environment makes the penta-coordinate structure thermodynamically most stable. This delicate difference between the tetra- and penta-coordinate structures, which depends on the environment, is a close reflection of the lower activation barrier of the SN 2 reaction found in neutral solvent or gas-phase reactions.

Coupling plasma physics and chemistry in the PIC model of electric propulsion: Application to an air-breathing, low-power Hall thruster

This work represents a first attempt to include the complex variety of electron-molecule processes in a full kinetic particle-in-cell/test particle Monte Carlo model for the plasma and neutral gas phase in a Hall thruster. Particular emphasis has been placed on Earth’s atmosphere species for the air-breathing concept. The coupling between the plasma and the gas phase is self-consistently captured by assuming the cold gas approximation and considering gas-wall and gas recycling from the walls due to ion neutralization. The results showed that, with air molecular propellants, all the most relevant thruster performance figures degraded relative to the nominal case using Xe propellant. The main reasons can be ascribed to a reduced ionization cross-section, a larger gas ionization mean free path due to lighter mass air species, and additional electron collisional power losses. While vibrational excitations power losses are negligible, dissociation and electronic excitations compete with the ionization channel. In addition, for molecular oxygen, the large dissociation leads to even faster atoms, further reducing their transit time inside the discharge channel. Future studies are needed to investigate the role of non-equilibrium vibrational kinetics and metastable states for stepwise ionization.

A Multiwavelength View of IC 860: What Is in Action inside Quenching Galaxies * * Herschel is an ESA space observatory with science instruments provided by European-led principal investigator consortia and with important participation from NASA.

We present a multiwavelength study of IC 860, a nearby post-starburst galaxy at the early stage of transitioning from blue and star forming to red and quiescent. Optical images reveal a galaxy-wide, dusty outflow originating from a compact core. We find evidence for a multiphase outflow in the molecular and neutral gas phase from the CO position–velocity diagram and NaD absorption features. We constrain the neutral mass outflow rate to be ∼0.5 M ⊙ yr−1, and the total hydrogen mass outflow rate to be ∼12 M ⊙ yr−1. Neither outflow component seems able to escape the galaxy. We also find evidence for a recent merger in the optical images, CO spatial distribution, and kinematics, and evidence for a buried active galactic nucleus in the optical emission line ratios, mid-IR properties, and radio spectral shape. The depletion time of the molecular gas reservoir under the current star formation rate is ∼7 Gyr, indicating that the galaxy could stay at the intermediate stage between the blue and red sequence for a long time. Thus the timescales for a significant decline in star formation rate (quenching) and gas depletion are not necessarily the same. Our analysis supports the quenching picture where outflows help suppress star formation by disturbing rather than expelling the gas and shed light on possible ongoing activities in similar quenching galaxies.

METAL: The Metal Evolution, Transport, and Abundance in the Large Magellanic Cloud Hubble Program. III. Interstellar Depletions, Dust-to-Metal, and Dust-to-Gas Ratios versus Metallicity

Abstract The metallicity and gas density dependence of interstellar depletions, the dust-to-gas (D/G), and dust-to-metal (D/M) ratios have important implications for how accurately we can trace the chemical enrichment of the universe, either by using FIR dust emission as a tracer of the ISM or by using spectroscopy of damped Lyα systems to measure chemical abundances over a wide range of redshifts. We collect and compare large samples of depletion measurements in the Milky Way (MW), Large Magellanic Cloud (LMC) (Z = 0.5 Z ⊙), and Small Magellanic Cloud (SMC) (Z = 0.2 Z ⊙). The relations between the depletions of different elements do not strongly vary between the three galaxies, implying that abundance ratios should trace depletions accurately down to 20% solar metallicity. From the depletions, we derive D/G and D/M. The D/G increases with density, consistent with the more efficient accretion of gas-phase metals onto dust grains in the denser ISM. For log N(H) &gt; 21 cm−2, the depletion of metallicity tracers (S, Zn) exceeds −0.5 dex, even at 20% solar metallicity. The gas fraction of metals increases from the MW to the LMC (factor 3) and SMC (factor 6), compensating for the reduction in total heavy element abundances and resulting in those three galaxies having the same neutral gas-phase metallicities. The D/G derived from depletions are respective factors of 2 (LMC) and 5 (SMC) higher than the D/G derived from FIR, 21 cm, and CO emission, likely due to the combined uncertainties on the dust FIR opacity and on the depletion of carbon and oxygen.

Cool outflows in MaNGA: a systematic study and comparison to the warm phase

ABSTRACT This paper investigates the neutral gas phase of galactic winds via the Na i Dλλ5890, 5895Å feature within z ∼ 0.04 MaNGA galaxies and directly compares their incidence and strength to the ionized winds detected within the same parent sample. We find evidence for neutral outflows in 127 galaxies (∼5 per cent of the analysed line-emitting sample). ${\rm Na\, \small {\rm I}\ D}$ winds are preferentially seen in galaxies with dustier central regions and both wind phases are more often found in systems with elevated star formation rate (SFR) surface densities, especially when there has been a recent upturn in the star formation activity according to the $\rm SFR_{5 \, Myr}/SFR_{800 \, Myr}$ parameter. We find the ionized outflow kinematics to be in line with what we measure in the neutral phase. This demonstrates that, despite their small contributions to the total outflow mass budget, there is value to collecting empirical measurements of the ionized wind phase to provide information on bulk motion in the outflow. Depending on dust corrections applied to the ionized gas diagnostics, the neutral phase has ∼1.2−1.8 dex higher mass outflow rates ($\dot{M}_{\rm out}$), on average, compared to the ionized phase. We quantify scaling relations between $\dot{M}_{\rm out}$ and the strengths of the physical wind drivers (SFR, LAGN). Using a radial–azimuthal stacking method, and by considering inclination dependencies, we find results consistent with biconical outflows orthogonal to the disc plane. Our work complements other multiphase outflow studies in the literature that consider smaller samples, more extreme objects, or proceed via stacking of larger samples.

Infrared Spectroscopy of Jet-cooled “GrandPAHs” in the 3–100 μm Region

Abstract Although large polycyclic aromatic hydrocarbons (PAHs) are likely to be responsible for IR emission of gaseous and dusty regions, their neutral experimental high-resolution gas-phase IR spectra—needed to construct accurate astronomical models—have so far remained out of reach because of their nonvolatility. Applying laser desorption to overcome this problem, we report here the first IR spectra of the jet-cooled large PAHs coronene (C24H12), peropyrene (C26H14), ovalene (C32H14), and hexa(peri)benzocoronene (C42H18) in the 3–100 μm region. Apart from providing experimental spectra that can be compared directly to astronomical data, such IR spectra are crucial for assessing the accuracy of theoretically predicted spectra used to interpret interstellar IR emission. Here we use the experimental spectra to evaluate the performance of conventional calculations using the harmonic approximation, as well as calculations with an anharmonic (GVPT2) treatment. The harmonic prediction agrees well with the experiment between 100 and 1000 cm−1 (100 and 10 μm) but shows significant shortcomings in the combination band (1600–2000 cm−1, 6.25–5 μm) and CH-stretch (2950–3150 cm−1, 3.4–3.17 μm) regions. Especially the CH-stretch region is known to be dominated by the effects of anharmonicity, and we find that large PAHs are no exception. However, for the CH out-of-plane region (667–1000 cm−1, 15–10 μm) the anharmonic treatment that significantly improves the predicted spectra for small PAHs leads to large and unrealistic frequency shifts, and intensity changes for large PAHs, thereby rendering the default results unreliable. A detailed analysis of the results of the anharmonic treatment suggests a possible route for improvement, although the underlying cause for the large deviations remains a challenge for theory.

Oscillation of Gas Density in the Gas Filament Remained by a Streamer Discharge in Water

When a streamer discharge occurs in water, several luminous plasma filaments will be created in the water during the discharge. After the discharge, these plasma filaments turn into neutral gas phase and remain in water. The gas filament remained in water is a good object for studying the basic processes involved in the streamer propagation. We investigated the evolution of the gas filaments remained in water after a streamer discharge at different experimental conditions. We recorded eight successive images during one discharge pulse. The density of gas in the gas filament and the radius of the gas filament were measured from the obtained images. We found that the radius of the gas filament and the density of gas in the gas filament are almost not influenced by the impulse voltage within the range studied. While the conductivity of water has strong effect on the radius of the gas filament and the density of gas in the gas filament. The radius of the gas filament becomes thicker and expands faster as the conductivity of water becomes larger. The density of gas in the gas filament remained in water oscillates between 400 to 800 kg/m3 with an duration of ~10 μs during the expansion period of 4–39 μs after the HV pulse starts. Both the impulse voltage and the conductivity of water do not affect the oscillation duration of the density of gas in the gas filament.

A simple method to calculate solution-phase free energies of charged species in computational electrocatalysis

Determining the adsorption potential of adsorbed ions in the field of computational electrocatalysis is of great interest to study their interaction with the electrode material and the solvent, and to map out surface phase diagrams and reaction pathways. Calculating the adsorption potentials of ions with density functional theory and comparing across various ions requires an accurate reference energy of the ion in solution and electrons at the same electrochemical scale. Here we highlight a previously used method for determining the reference free energy of solution phase ions using a simple electrochemical thermodynamic cycle, which allows this free energy to be calculated from that of a neutral gas-phase or solid species and an experimentally measured equilibrium potential, avoiding the need to model solvent around the solution phase ion in the electronic structure calculations. While this method is not new, we describe its use and utility in detail and show that this same method can be used to find the free energy of any ion from any reaction, as long as the half-cell equilibrium potential is known, even for reactions that do not transfer the same number of protons and electrons. To illustrate its usability, we compare the adsorption potentials obtained with DFT of I*, Br*, Cl*, and SO4 * on Pt(111) and Au(111) and OH* and Ag* on Pt(111) with those measured experimentally and find that this simple and computationally affordable method reproduces the experimental trends.

Diagnosing the interstellar medium of galaxies with far-infrared emission lines

Context. Atomic fine structure lines have been detected in the local Universe and at high redshifts over the past decades. The [C II] emission line at 158 μm is an important observable as it provides constraints on the interstellar medium (ISM) cooling processes. Aims. We develop a physically motivated framework to simulate the production of far-infrared line emission from galaxies in a cosmological context. This first paper sets out our methodology and describes its first application: simulating the [C II] 158 μm line emission in the local Universe. Methods. We combine the output from EAGLE cosmological hydrodynamical simulations with a multi-phase model of the ISM. Gas particles are divided into three phases: dense molecular gas, neutral atomic gas, and diffuse ionised gas (DIG). We estimate the [C II] line emission from the three phases using a set of CLOUDY cooling tables. Results. Our results agree with previous findings regarding the contribution of these three ISM phases to the [C II] emission. Our model shows good agreement with the observed L[C II–star formation rate (SFR) relation in the local Universe within 0.4 dex scatter. Conclusions. The fractional contribution to the [C II] line from different ISM phases depends on the total SFR and metallicity. The neutral gas phase dominates the [C II] emission in galaxies with SFR ∼ 0.01–1 M⊙ yr−1, but the ionised phase dominates at lower SFRs. Galaxies above solar metallicity exhibit lower L[C II/SFR ratios for the neutral phase. In comparison, the L[C II/SFR ratio in the DIG is stable when metallicity varies. We suggest that the reduced size of the neutral clouds, caused by increased SFRs, is the likely cause for the L[C II deficit at high infrared luminosities, although EAGLE simulations do not reach these luminosities at z = 0.

Thermal desorption effects on fragment ion production from multi-photon ionized uridine and selected analogues.

Experiments on neutral gas-phase nucleosides are often complicated by thermal lability. Previous mass spectrometry studies of nucleosides have identified enhanced relative production of nucleobase ions (e.g. uracil+ from uridine) as a function of desorption temperature to be the critical indicator of thermal decomposition. On this basis, the present multi-photon ionization (MPI) experiments demonstrate that laser-based thermal desorption is effective for producing uridine, 5-methyluridine, and 2′-deoxyuridine targets without thermal decomposition. Our experiments also revealed one notable thermal dependence: the relative production of the sugar ion C5H9O4+ from intact uridine increased substantially with the desorption laser power and this only occurred at MPI wavelengths below 250 nm (full range studied 222–265 nm). We argue that this effect can only be rationalized plausibly in terms of changing populations of different isomers, tautomers, or conformers in the target as a function of the thermal desorption conditions. Furthermore, the wavelength threshold behavior of this thermally-sensitive MPI channel indicates a critical dependence on neutral excited state dynamics between the absorption of the first and second photons. The experimental results are complemented by density functional theory (DFT) optimizations of the lowest-energy structure of uridine and two further conformers distinguished by different orientations of the hydroxymethyl group on the sugar part of the molecule. The energies of the transitions states between these three conformers are low compared with the energy required for decomposition.

The Quest for a Plausible Formation Route of Formyl Cyanide in the Interstellar Medium: a State-of-the-art Quantum-chemical and Kinetic Approach

Abstract Interstellar complex organic molecules are assumed to be mainly formed on dust–grain surfaces. However, neutral gas-phase reactions in the interstellar medium can play an important role. In this paper, by investigating the reaction between aldehydes and the cyano radical, we show that both formaldehyde (CH2O) and acetaldehyde (CH3CHO) can lead to the formation of formyl cyanide (HCOCN). Owing to accurate quantum-chemical computations followed by rate constant evaluations, we have been able to suggest and validate an effective mechanism for the formation of HCOCN, one of the molecules observed in the ISM. Quite interestingly, the mechanism starting from CH2O is very effective at a low temperature, while that involving CH3CHO becomes more efficient at temperatures above 200 K.

A twist on the reaction of the CN radical with methylamine in the interstellar medium: new hints from a state-of-the-art quantum-chemical study

ABSTRACTDespite the fact that the majority of current models assume that interstellar complex organic molecules (iCOMs) are formed on dust–grain surfaces, there is some evidence that neutral gas-phase reactions play an important role. In this paper, we investigate the reaction occurring in the gas phase between methylamine (CH3NH2) and the cyano (CN) radical, for which only fragmentary and/or inaccurate results have been reported to date. This case study allows us to point out the pivotal importance of employing quantum-chemical calculations at the state of the art. Since the two major products of the CH3NH2 + CN reaction, namely the CH3NH and CH2NH2 radicals, have not been spectroscopically characterized yet, some effort has been made for filling this gap.

Multiphase outflows in post-starburst E+A galaxies - II. A direct connection between the neutral and ionized outflow phases

ABSTRACT Post-starburst E+A galaxies are systems that hosted a powerful starburst that was quenched abruptly. Simulations suggest that these systems provide the missing link between major merger ULIRGs and red and dead ellipticals, where AGN feedback is responsible for the expulsion or destruction of the molecular gas. However, many details remain unresolved and little is known about AGN-driven winds in this short-lived phase. We present spatially resolved IFU spectroscopy with MUSE/VLT of SDSS J124754.95-033738.6, a post-starburst E+A galaxy with a recent starburst that started 70 Myr ago and ended 30 Myr ago, with a peak SFR of $\sim 150\, \mathrm{M_{\odot }\,yr^{ -1}}$. We detect disturbed gas throughout the entire field of view, suggesting triggering by a minor merger. We detect fast-moving multiphased gas clouds, embedded in a double-cone face-on outflow, which are traced by ionized emission lines and neutral NaID emission and absorption lines. We find remarkable similarities between the kinematics, spatial extents, and line luminosities of the ionized and neutral gas phases, and propose a model in which they are part of the same outflowing clouds, which are exposed to both stellar and AGN radiation. Our photoionization model provides consistent ionized line ratios, NaID absorption optical depths and EWs, and dust reddening. Using the model, we estimate, for the first time, the neutral-to-ionized gas mass ratio (about 20), the sodium neutral fraction, and the size of the outflowing clouds. This is one of the best ever observed direct connections between the neutral and ionized outflow phases in AGN.