The EDGE-CALIFA collaboration has obtained resolved IFU and CO spectroscopy for 126 nearby galaxies, selected from the CALIFA main sample. We have assembled the spatially resolved products of the survey into Astropy-compatible pixel tables that reduce the oversampling in the original FITS images and facilitate comparison of pixels across different images. By joining these pixel tables to lower dimensional tables that provide profiles, integrated spectra, or global properties, it is possible to investigate the dependence of star formation rate on both local and global conditions. In this short video talk I provide an overview of EDGE, with examples of the use of the database to conduct analysis and generate plots.

We use the AREPO numerical code to model the structure of a Milky Way like galaxy (MW) via a suite of simulations composed of a stellar disc and bulge, a dark matter halo, and a gaseous disc under isothermal conditions. For each model, we produce longitude velocity (l-v) maps of the gas surface densities to extract the skeletons of the main features (arms, bar), and the contours defining the terminal velocities of the gas. We compare these with observations via a number of diagnostic tools, and select the model that best reproduces the main observed features of the Milky Way.

We employ the Feedback In Realistic Environments (FIRE-2) physics model to study how the properties of giant molecular clouds (GMCs) evolve during galaxy mergers. Due to the rarity of mergers in the local Universe, samples of nearby merging galaxies suitable for studies of individual GMCs are limited. Idealized simulations provide us with a new window to study GMC evolution during a merger, and assist in interpreting observations. We conduct a pixel-by-pixel analysis of the simulated molecular gas properties in both undisturbed control galaxies and galaxy mergers. The simulated GMC-pixels follow a similar trend in a diagram of velocity dispersion (σv) versus gas surface density (Σmol) as observed in normal spiral galaxies in the Physics at High Angular resolution in Nearby GalaxieS (PHANGS) survey. For simulated mergers, we see a significant increase in both the Σmol and σv for GMC-pixels by a factor of 5 – 10, which put these pixels to be above the trend of PHANGS galaxies in the σv vs Σmol diagram. This deviation indicates that GMCs in the simulated merger are more gravitationally unbound and have higher virial parameter (αvir) of 10 – 100, which is much larger than that of simulated control galaxies. Furthermore, we find that the increase in αvir generally happens at the same time as the increase in global star formation rate (SFR), which suggests feedback is playing a role in dispersing the gas. The correspondence between high αvir and SFR also suggests some other physical mechanisms besides self-gravity are helping the GMCs in starburst mergers to collapse and form stars.

We present new neutral hydrogen (Hi) observations of the nearby galaxy NGC 2403 to determine the nature of a low-column density cloud that was detected earlier by the Green Bank Telescope.

We find that this cloud is the tip of a complex of filaments of extraplanar gas that is coincident with the main disk. The total Hi mass of the complex is 2 × 107 M⊙ or 0.6% of the total Hi mass of the galaxy. The main structure, previously referred to as the 8-kpc filament, is now seen to be even more extended, along a 20 kpc stream.

Feedback from supernovae (SNe) is an essential mechanism that self-regulates the growth of galaxies. We build an SN feedback model based on high-resolution simulations of superbubble and SN-driven outflows for the physical understanding of the galaxy–CGM connection. Using an Eulerian hydrodynamic code Athena++, we find universal scaling relations for the time evolution of superbubble momentum, when the momentum and time are scaled by those at the shell-formation time. We then develop an SN feedback model utilizing Voronoi tessellation, and implement it into the GADGET3-Osaka smoothed particle hydrodynamic code. We show that our stochastic thermal feedback model produces galactic outflow that carries the metals high above the galactic plane but with weak suppression of star formation. Additional mechanical feedback further suppresses star formation. Therefore, we argue that both thermal and mechanical feedback is necessary for the SN feedback model of galaxy evolution when an individual SN bubble is unresolved.

The winds observed around asymptotic giant branch (AGB) stars are generally attributed to radiation pressure on dust, which is formed in the extended dynamical atmospheres of these pulsating, strongly convective stars. Current radiation-hydrodynamical models can explain many of the observed features, and they are on the brink of delivering a predictive theory of mass loss. This review summarizes recent results and ongoing work on winds of AGB stars, discussing critical ingredients of the driving mechanism, and first results of global 3D RHD star-and-wind-in-a-box simulations. With such models it becomes possible to follow the flow of matter, in full 3D geometry, all the way from the turbulent, pulsating interior of an AGB star, through its atmosphere and dust formation zone into the region where the wind is accelerated by radiation pressure on dust. Advanced instruments, which can resolve the stellar atmospheres, where the winds originate, provide essential data for testing the models.

The unidentified infrared (UIR) bands, whose carriers are thought to be organics, have been widely observed in various astrophysical environments. However, our knowledge of the detailed chemical composition and formation process of the carriers is still limited. We have synthesized laboratory organics named Quenched Nitrogen-included Carbonaceous Composite (QNCC) by quenching plasma produced from nitrogen gas and hydrocarbon solids. Infrared and X-ray analyses of QNCC showed that infrared properties of QNCC well reproduce the UIR bands observed in novae and amine structures contained in QNCC play an important role in the origin of the broad 8 m feature, which characterizes the UIR bands in novae. QNCC is at present the best laboratory analog of organic dust formed around dusty classical novae, which carries the UIR bands in novae via thermal emission process [Endo et al.(2021)].

Most of a star’s mass is bound in a hydrostatic equilibrium in which pressure balances gravity. But if at some near-surface layer additional outward forces overcome gravity, this can transition to a supersonic, outflowing wind, with the sonic point, where the outward force cancels gravity, marking the division between hydrostatic atmosphere and wind outflow. This talk will review general issues with such transonic initiation of a stellar wind outflow, and how this helps set the wind mass loss rate. The main discussion contrasts the flow initiation in four prominent classes of steady-state winds: (1) the pressure-driven coronal wind of the sun and other cool stars; (2) line-driven winds from OB stars; (3) a two-stage initiation model for the much denser winds from Wolf-Rayet (WR) stars; and (4) the slow “overflow” mass loss from highly evolved giant stars. A follow on discussion briefly reviews eruptive mass loss, with particular focus on the giant eruption of η Carinae.

Classical Wolf-Rayet (WR) stars mark an important stage in the late evolution of massive stars. As hydrogen-poor massive stars, these objects have lost their outer layers, while still losing further mass through strong winds indicated by their prominent emission line spectra. Wolf-Rayet stars have been detected in a variety of different galaxies. Their strong winds are a major ingredient of stellar evolution and population synthesis models. Yet, a coherent theoretical picture of their strong mass-loss is only starting to emerge. In particular, the occurrence of WR stars as a function of metallicity (Z) is still far from being understood.

To uncover the nature of the complex and dense winds of Wolf-Rayet stars, we employ a new generation of model atmospheres including a consistent solution of the wind hydrodynamics in an expanding non-LTE situation. With this technique, we can dissect the ingredients driving the wind and predict the resulting mass-loss for hydrogen-depleted massive stars. Our modelling efforts reveal a complex picture with strong, non-linear dependencies on the luminosity-to-mass ratio and Z with a steep, but not totally abrupt onset for WR-type winds in helium stars. With our findings, we provide a theoretical motivation for a population of helium stars at low Z, which cannot be detected via WR-type spectral features. Our study of massive He-star atmosphere models yields the very first mass-loss recipe derived from first principles in this regime. Implementing our first findings in stellar evolution models, we demonstrate how traditional approaches tend to overpredict WR-type mass loss in the young Universe.

To look at propagating winds from evolved stars into the interstellar medium is to look at how they are sustained. To understand their origins, we must look to the circumstances that create them in the first instance. In this article, I examine the physical conditions under which pulsation-enhanced, dust-driven winds are first generated. These initial conditions can help constrain the late evolutionary stages of these stars and provide insight into the mechanisms that cause the mass loss itself.

Progenitors of Type Ib and Ic supernovae (SNe) are stripped envelope stars and provide important clues on the mass-loss history of massive stars. Direct observations of the progenitors before the supernova explosion would provide strong constraints on the exact nature of SN Ib/Ic progenitors. Given that stripped envelope massive stars can have an optically thick wind as in the case of Wolf-Rayet stars, the influence of the wind on the observational properties needs to be properly considered to correctly infer progenitor properties from pre-SN observations. Non-LTE stellar atmosphere models indicate that the optical brightness could be greatly enhanced with an optically thick wind because of lifting-up of the photosphere from the stellar surface to the wind matter, and line and free-free emissions. So far, only a limited number of SN Ib/Ic progenitor candidates have been reported, including iPTF13bvn, SN 2017ein and SN 2019yvr. We argue that these three candidates are a biased sample, being unusually bright in the optical compared to what is expected from typical SN Ib/Ic progenitors, and that mass-loss enhancement during the final evolutionary stage can explain their optical properties.

Water is a ubiquitous molecule in circumstellar envelopes (CSEs). Its emission has been detected at a wide range of distances from the central oxygen-rich evolved star. In particular, the water maser transition at 22 GHz, typically extending from about 5–20 stellar radii to as far as several hundred stellar radii from the star, has been commonly used to probe the structure and dynamics of the intermediate regions of the CSE where dust is condensing and the inner wind is being accelerated. The advent of ALMA has opened the door to high-angular resolution mapping of much higher excitation transitions of water, probing the inner regions of the CSEs, some of which are anticipated to exhibit maser action. The ALMA ATOMIUM large program observed many such transitions towards a sample of AGB stars & red supergiants. The preliminary results show that while some transitions depart only slightly from LTE, others clearly show signs of maser action. The Gaussian fitting of the non-diffuse/compact part of some of the (quasi) thermal & maser transitions reveal interesting velocity gradients, signatures of outflowing and infalling motions hence providing important constraints for stellar wind models.

The cosmic origin of fluorine is still under debate. Asymptotic giant branch (AGB) stars are among the few suggested candidates to efficiently synthesis F in our Galaxy, however their relative contribution is not clear. In this paper, we briefly review the theoretical studies from stellar yield models of the F synthesis and chemical equilibrium models of the F-containing molecules in the outflow around AGB stars. Previous detections of the F-bearing species towards AGB and post-AGB stars are also highlighted. We suggest that high-resolution ALMA observations of the AlF, one of the two main carriers of F in the outflow of AGB stars, can provide a reliable tracer of the F-budget in AGB stars. This will be helpful to quantify the role of AGB stars in the Galactic F budget.

Astrochemical models treat dust surfaces as ice covered. We investigate the effects of implementing increased bare dust binding energies of CO and S-bearing species on the chemistry in the outflows of asymptotic giant branch (AGB) stars. We demonstrate the potential for improving agreement with observations in the outflow of IK Tau.

Increasing the binding energies to measured and computationally derived values in high mass-loss AGB outflows increased the production of daughter species. Switching from a high binding energy on bare dust to weaker binding to ice, the gas phase abundance increased at a radius in agreement with observations of IK Tau, suggesting that displacement of bound species could contribute to this observational puzzle. Using a strong binding to bare dust, a gas phase increase was not observed, however parent species concentrations had to be increased by around a factor of four to explain observed concentrations.

We present interferometric continuum and molecular line emission maps obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) of OH231.8+4.2, a well studied bipolar nebula around an asymptotic giant branch (AGB) star that is key to understand the remarkable changes in nebular morphology and kinematics during the short transition from the AGB to the Planetary Nebula (PN) phase. The excellent angular resolution of our maps (∼20 mas ∼30 AU) allows us to scrutinize the central nebular regions of OH231.8+4.2, which hold the clues to understanding how this iconic object assembled its complex nebular architecture. We report, for the first time in this object and others of its kind (i.e. pre-PNe with massive bipolar outflows), the discovery of a rotating circumbinary disk of radius ∼30 AU traced by NaCl, KCl, and H2O emission lines. The disk lies at the base of a young bipolar wind with signs of rotation as well. A compact spatially resolved dust disk is found perpendicular to the bipolar outflow. We also identify a point-like continuum source, which likely represents the central Mira star enshrouded by a ∼3 R* shell or disk of hot (∼1400 K) freshly formed dust. The point source is slightly off-centre from the disk centroid, enabling us for the first time to place constraints to the orbital separation of the central binary system.

Binary interaction with a stellar or planetary companion has been proposed to be the driving mechanism behind large-scale asymmetries, such as spirals and disks, observed within AGB outflows. We developed the first chemical kinetics model that takes the effect of a stellar companions’s UV radiation into account. The presence of a stellar companion can initiate a rich photochemistry in the inner wind. Its impact is determined by the intensity of the UV radiation and the extinction the radiation experiences. The outcome of the inner wind photochemistry depends on the balance between two-body reactions and photoreactions. If photoreactions dominate, the outflow can appear molecule-poor. If two-body reactions dominate, chemical complexity within the outflow can increase, yielding daughter species with a large inner wind abundance. A comprehensive view on the molecular content of the outflow, especially combined with abundance profiles, can point towards the presence of a stellar companion.

One of the big challenges for 21st century stellar astrophysics is the impact of binary interactions on stellar evolution. Such interactions are believed to play a key role in the death throes of 1-8 M⊙ stars, as they evolve from the AGB stars into Planetary Nebulae. X-ray surveys of UV-emitting AGB stars show that ∼40% of objects with FUV emission and GALEX FUV/NUV flux ratios ≳0.2 have variable X-ray emission characterized by very high temperatures (Tx∼35-160 MK) and luminosities (Lx∼0.002-0.2L⊙). We hypothesize that such AGB stars have accretion and (accretion-powered) outflows associated with a close binary companion. UV spectroscopy with HST/STIS of our brightest object (Y Gem) shows the presence of infalling and outflowing gas, providing direct kinematic confirmation of this hypothesis. However, the UV-emitting AGB star population is dominated by objects with little or no FUV emission, and we do not know whether the UV emission from these is intrinsic to the AGB star or extrinsic (i.e., due to binarity). Here we present the first results from a large grid of simple chromospheric models to help discriminate between the intrinsic and extrinsic mechanisms of UV emission for AGB stars.

Thanks to the long-term collaborations between nuclear and astrophysics, we have good understanding on the origin of elements in the universe, except for the elements around Ti and some neutron-capture elements. From the comparison between observations of nearby stars and Galactic chemical evolution models, a rapid neutron-capture process associated with core-collapse supernovae is required. The production of C, N, F and some minor isotopes depends on the rotation of massive stars, and the observations of distant galaxies with ALMA indicate rapid cosmic enrichment. It might be hard to find very metal-poor or Population III (and dust-free) galaxies at very high redshifts even with JWST.

We run numerical simulations of massive colliding wind binaries, and quantify the accretion onto the secondary under different conditions. We set 3D simulation of a LBV–WR system and vary the LBV mass loss rate to obtain different values of wind momentum ratio η. We show that the mean accretion rate for stationary systems fits a power law Macc∝ η–1.6 for a wide range of η, until for extremely small η saturation in the accretion is reached. We find that the stronger the primary wind, the smaller the opening angle of the colliding wind structure (CWS), and compare it with previous analytical estimates. We demonstrate the efficiency of clumpy wind in penetrating the CWS and inducing smaller scale clumps that can be accreted. We propose that simulations of colliding winds can reveal more relations as the ones we found, and can be used to constrain stellar parameters.

Maser properties can be measured with milli-arcsec precision over multiple epochs using ALMA, cm- and mm-wave VLBI and e-MERLIN. This allows: (i) Tracing SiO maser proper motions in the pulsation-dominated zone; (ii) Quantifying clumpiness, variability and asymmetry of the wind traced by masers; (iii) Contrasting behaviour from OH masers even at similar distances from the star; (iv) Measuring magnetic fields. Mass lost from the star, traced by SiO masers, is likely to take decades to reach ∼5 stellar radii. At 5–50 stellar radii, once dust is well formed, 22-GHz H2O masers show the wind accelerating through the escape velocity; its overall direction is away from the star but the velocity field is complex. In a few cases (so far), highly-directed, localised ejecta are seen. Magnetic fields appear to be stellar-centred and strong enough to influence wind kinematics. Recent ALMA and other observations have shown that otherwise inconspicuous companions shape a majority of evolved star winds, whilst advanced models demonstrate how, for some situations, this is compatible with masers showing negligible rotation proper motions. The long-term monitoring achievable at radio frequencies complements the multi-transition maser studies and analysis of thermal lines and dust at shorter wavelengths.