Galactic Disk Research Articles

Scaling relations between mass components and rotational velocity in disk galaxies

Abstract In this study, we re-evaluate the baryonic Tully–Fisher relation (BTFR) by analyzing the correlations between maximum rotational velocity and various mass components, including stellar mass, atomic hydrogen (H i) mass, baryonic mass, and dark matter mass in a sample of 141 disk galaxies from the SPARC (Spitzer Photometry and Accurate Rotation Curves) database, augmented by recent data on stellar and dark matter masses. We apply multiple statistical methods, including Monte Carlo orthogonal distance regression (MCODR), Monte Carlo least-squares (MCLS), and traditional least-squares (LS), to assess the impact of different fitting techniques on the derived scaling relations between the mass components (stellar, H i, and dark matter) and maximum rotational velocities of these galaxies. We find that the selection of statistical methods significantly influences the derived slopes and intercepts the relation between maximum rotational velocity and mass components. The MCODR method that accounts for errors in both variables consistently produces steeper slopes, suggesting a stronger correlation between stellar mass and rotational dynamics compared to other methods. In contrast, the MCLS method tends to yield flatter slopes, highlighting the sensitivity of this approach to outliers. Despite the variations in slope and intercept across different methods, the fundamental relation between baryonic mass and rotational velocity remains robust. We have also compared dark matter mass derived from different halo models [NFW (Navarro–Frenk–White) versus combined NFW + Dekel–Zhao profiles] and noted that the slope from the NFW profile is slightly steeper than that from the combined profile, highlighting the sensitivity of scaling relations to the selection of halo model. Overall, this study reinforces the robustness of the BTFR across different mass components in disk galaxies while emphasizing the critical role of statistical methods and dark matter profiles in analyzing galactic dynamics.

Disc distortion revisited

Abstract We revisit the dynamics of razor-thin, stone-cold, self-gravitating discs. By recasting the equations into standard cylindrical coordinates, the linearised vertical dynamics of an exponential disc can be followed for several gigayears on a laptop in a few minutes. An initially warped disc rapidly evolves into a flat inner region and an outward-propagating spiral corrugation wave that rapidly winds up and would quickly thicken a disc with non-zero radial velocity dispersion. The Sgr dwarf galaxy generates a similar warp in the Galactic disc as it passes through pericentre, and the warp generated by the dwarf’s last pericentre ∼35 Myr ago is remarkably similar to the warp traced by the Galaxy’s HI disc. The resemblance to the observed warp is fleeting but its timing is perfect. For the adopted parameters the amplitude of the model warp is a factor 3 too small, but there are several reasons for this being so. The marked flaring of our Galaxy’s low-α disc just outside the solar circle can be explained as a legacy of earlier pericentres.

When LAMOST meets Gaia DR3. Exploring the metallicity of open clusters

Open clusters (OCs) are excellent probes as their age and abundance can be tightly constrained, allowing us to explore the distribution of metallicity and composition across the disk of the Milky Way. By conducting a comprehensive analysis of the metallicity of OCs, we can obtain valuable information about the history of their chemical enrichment. Moreover, by observing stars in different regions of the Milky Way, we can identify significant spatial structures in their chemical composition and abundance. This enables us to understand stellar birth radii through chemical tagging. Nevertheless, it remains challenging to infer the original positions of OCs using current data alone. The aim of this study is to investigate the distribution of metallicity in the solar neighborhood using a large dataset from Gaia DR3 combined with LAMOST spectra. With accurate ages and metallicity measurements, we can determine birth radii for the stars and attempt to understand their migration pattern. We chose a total of 1131 OCs within 3 Kpc of the Sun from the Gaia DR3 and LAMOST DR8 low-resolution spectral database (R=1800). We used an artificial neural network to correct the LAMOST data by incorporating high-resolution spectral data from GALAH DR3 (R=28000). The average metallicity of the OCs was determined based on the reliable Fe/H values for their members. We then examined the distribution of metallicity across different regions within the Galaxy and inferred birth radii of the OCs from their age and metallicity. The correction method presented here can partially eliminate the systematic offset for LAMOST data. We discuss the metallicity trend as a function of Galactocentric distance and the guiding radii. We also compare these observational results with those from chemo-dynamic simulations. Values derived from observational metallicity data are slightly lower than predicted values when the uncertainties are not considered. However, the metallicity gradients are consistent with previous calculations. Finally, we investigated the birthplace of OCs and find hints that the majority of OCs near the Sun have migrated from the outer Galactic disk.

Ultraviolet Radiation Fields in Star-forming Disk Galaxies: Numerical Simulations with TIGRESS-NCR

With numerical simulations that employ adaptive ray-tracing (ART) for radiative transfer at the same time as evolving gas magnetohydrodynamics, thermodynamics, and photochemistry, it is possible to obtain a high-resolution view of ultraviolet (UV) fields and their effects in realistic models of the multiphase interstellar medium. Here, we analyze results from TIGRESS-NCR simulations, which follow both far-UV (FUV) wavelengths, important for photoelectric heating and polycyclic aromatic hydrocarbon excitation, and the Lyman continuum (LyC), which photoionizes hydrogen. Considering two models, representing solar neighborhood and inner-galaxy conditions, we characterize the spatial distribution and time variation of UV radiation fields, and quantify their correlations with gas. We compare four approximate models for the FUV to simulated values to evaluate alternatives when full ART is infeasible. By convolving FUV radiation with density, we produce mock maps of dust emission. We introduce a method to calibrate mid-IR observations, for example from JWST, to obtain high-resolution gas surface density maps. We then consider the LyC radiation field, finding most of the gas exposed to this radiation to be in ionization–recombination equilibrium and to have a low neutral fraction. Additionally, we characterize the ionization parameter as a function of the environment. Using a simplified model of the LyC radiation field, we produce synthetic maps of emission measure (EM). We show that the simplified model can be used to extract an estimate of the neutral fraction of the photoionized gas and mean free path of ionizing radiation from observed EM maps in galaxies.

Photometry and kinematics of dwarf galaxies from the Apertif HI survey

Understanding the dwarf galaxy population in low density environments (in the field) is crucial for testing the current Lambda Cold Dark Matter cosmological model. The increase in diversity toward low-mass galaxies is seen as an increase in the scatter of scaling relations, such as the stellar mass--size and the baryonic Tully-Fisher relation (BTFR), and is also demonstrated by recent in-depth studies of an extreme sub-class of dwarf galaxies with low surface brightnesses but large physical sizes called ultra-diffuse galaxies (UDGs). We aim to select dwarf galaxies independent of their stellar content and to make a detailed study of their gas and stellar properties. We selected galaxies from the APERture Tile In Focus (Apertif) survey and applied a constraint on their $i$-band absolute magnitude in order to exclude high-mass systems. The sample consists of 24 galaxies, 22 of which are resolved in by at least three beams, and they span mass ranges of 8.6 $ lesssim 9.7 and a stellar mass range of 8.0 $ lesssim 9.7 (with only three galaxies having log ( msun)>9). We determined the geometrical parameters of the and stellar disks, built kinematic models from the data using and extracted surface brightness profiles in the g-, r- and i- bands from the Pan-STARRS 1 photometric survey. We used these measurements to place our galaxies on the stellar mass--size relation and the BTFR, and we compared them with other samples from the literature. We find that at a fixed stellar mass, our dwarfs have larger optical effective radii than isolated optically selected dwarfs from the literature, and we found misalignments between the optical and morphologies for some of our sample. For most of our galaxies, we used the morphology to determine their kinematics, and we stress that deep optical observations are needed to trace the underlying stellar disks. Standard dwarfs in our sample follow the same BTFR of high-mass galaxies, whereas UDGs are slightly offset toward lower rotational velocities, in qualitative agreement with results from previous studies. Finally, our sample features a fraction (25<!PCT!>) of dwarf galaxies in pairs that is significantly larger with respect to previous estimates based on optical spectroscopic data.

Bidimensional Exploration of the warm-Temperature Ionised gaS (BETIS)

The extraplanar diffuse ionised gas (eDIG) is a key component for understanding the feedback processes that connect galactic discs and their halos. In this paper, we present the second study of the Bidimensional Exploration of the warm-Temperature Ionised Gas (BETIS) project, the aim of which is to explore the possible ionisation mechanisms and characteristics of the eDIG. We use a sample of eight edge-on galaxies observed with the Multi-Unit Spectroscopic Explorer (MUSE) integral field spectrograph (IFS) and apply the methodology developed in the first paper of the BETIS project for obtaining binned emission line maps. We find that the vertical and radial profiles of the [N II]/Hα, [S II]/Hα, [O III]/Hβ, and [O I]/Hα ratios depict a complex ionisation structure within galactic halos – which is influenced by the spatial distribution of H II regions across the galactic plane as observed from our line of sight–, with Lyman continuum photon leakage from OB associations constituting the main ionisation source. Moreover, the electron temperature and S+/S ionisation ratio also exhibit a dependency on the distribution of H II regions within the galactic discs. Our analysis excludes low-mass, hot, and evolved stars (HOLMES) as viable candidates for secondary ionisation sources to elucidate the unusual behaviour of the line ratios at greater distances from the galactic midplane. In contrast, we ascertain that shocks induced in the interstellar medium by star formation(SF)-related feedback mechanisms represent a promising secondary ionisation source of the eDIG. We present a suite of models integrating ionisation mechanisms arising from fast shocks and photoionisation associated with star formation. When applied to the classical Baldwin–Phillips–Terlevich (BPT) diagrams, these models reveal that the ionisation budget of the eDIG ranges from 20% to 50% across our sample, with local variations of up to 20% within individual galaxy halos. This contribution correlates with the presence of filaments and other structural components observed within galaxy halos. The presence of shocks is additionally supported by the observation of a high density of high [O I]/Hα ratios, which is characteristic of shock-compressed ionised gas, and is likely induced by feedback from regions of intense SF within the galactic disc. These results demonstrate consistency across all galaxies analysed in this sample.

Origin of the spectral features observed in the cosmic-ray spectrum

Recent measurements reveal the presence of several features in the cosmic-ray (CR) spectrum. In particular, the proton and helium spectra exhibit a spectral hardening at $ GV $ and a spectral steepening at $ TV $, followed by the well-known knee-like feature at $ PV $. The spectra of heavier nuclei also harden at $ GV $, while no claim can be currently made about the presence of the $ TV $ softening, due to low statistics. In addition, the B/C ratio also exhibits a hardening at $ GeV/n $ and seems to be rather shallow at $ TeV/n We propose a possible explanation of the observed spectral features in the framework of a composite diffusion scenario and considering different classes of sources. The proposed scenario is based on two assumptions. First, in the Galactic disk, where magnetic field lines are mainly oriented along the Galactic plane, particle scattering is assumed to be very inefficient. Therefore, the transport of CRs from the disk to the halo is set by the magnetic field line random walk induced by large-scale turbulence. Second, we propose that the spectral steepening at $ TV $ is related to the typical maximum rigidity reached in the acceleration of CRs by the majority of supernova remnants, while we assume that only a fraction of sources, contributing to $ 10-20<!PCT!> $ of the CR population, can accelerate particles up to sim PV rigidities. Within this framework we show that it is possible to reproduce the proton and helium spectra from GV to multi-PV; the $p$/He ratio; the spectra of CRs from lithium to iron; the $ p $ flux and the $ p /p$ ratio; and the abundance ratios B/C, B/O, C/O, Be/C, Be/O, and Be/B. We also discuss the Be Be ratio in view of the recent AMS02 preliminary measurements.

Recovering chemical bimodalities in observed edge-on stellar disks: Insights from AURIGA simulations

The well-known bimodal distribution of Milky Way disk stars in the [α/Fe]–metallicity plane is often used to define thick and thin disks. In external edge-on galaxies, there have been attempts to identify this type of bimodality using integral field spectroscopy (IFS) data. However, for unresolved stellar populations, observations only contain integrated information, making these studies challenging. We assessed the ability to recover chemical bimodalities in IFS observations of edge-on galaxies, using 24 Milky Way-mass galaxies from the AURIGA zoom-in cosmological simulations. We first analyzed the distribution of single stellar particles in the [Mg/Fe]–[Fe/H] plane, finding that bimodality is frequent but not ubiquitous and often unclear. Then we produced mock IFS [Mg/Fe] and [Fe/H] maps of galaxies seen edge on, and considered integrated stellar-population properties (projected and spatially binned). We investigated how the distribution of stars in the [Mg/Fe]–[Fe/H] plane is affected by edge-on projection and spatial binning. Bimodality is preserved, while distributions change their shapes. Naturally, broad distributions of individual star particles are narrowed into smaller [Mg/Fe] and [Fe/H] ranges for spatial bins. We observe continuous distributions from high [Mg/Fe] and low [Fe/H], to lower [Mg/Fe] values and higher [Fe/H]. Despite being continuous, these distributions are bimodal in most cases. The overlap in [Fe/H] is small, and different [Mg/Fe] components show up as peaks instead of sequences (even when the latter are present for individual particles). The larger the spatial bins, the narrower the [Mg/Fe]–[Fe/H] distribution. This narrowing helps amplify the density of different [Mg/Fe] peaks, often leading to a clearer bimodality in mock IFS observations than for original star particles. We also assessed the correspondence of chemical bimodalities with the distinction between geometric thick and thin disks. Their individual particles have different distributions, but mostly overlap in [Mg/Fe] and [Fe/H]. However, integrated properties of geometric thick and thin disks in mock maps do mostly segregate into different regions of the [Mg/Fe]–[Fe/H] plane. In bimodal distributions, they correspond to the two distinct peaks. Our results show that this approach can be used for bimodality studies in future IFS observations of edge-on external galaxies.

Learning the Universe: GalactISM Simulations of Resolved Star Formation and Galactic Outflows across Main-sequence and Quenched Galactic Environments

We present a suite of six high-resolution chemodynamical simulations of isolated galaxies, spanning observed disk-dominated environments on the star-forming main sequence, as well as quenched, bulge-dominated environments. We compare and contrast the physics driving star formation and stellar feedback among the galaxies, with a view to modeling these processes in cosmological simulations. We find that the mass loading of galactic outflows is coupled to the clustering of supernova explosions, which varies strongly with the rate of galactic rotation Ω = v circ/R via the Toomre length, leading to smoother gas disks in the bulge-dominated galaxies. This sets an equation of state in the star-forming gas that also varies strongly with Ω, so that the bulge-dominated galaxies have higher midplane densities, lower velocity dispersions, and higher molecular gas fractions than their main-sequence counterparts. The star formation rate in five out of six galaxies is independent of Ω and is consistent with regulation by the midplane gas pressure alone. In the sixth galaxy, which has the most centrally concentrated bulge and thus the highest Ω, we reproduce dynamical suppression of the star formation efficiency in agreement with observations. This produces a transition away from pressure-regulated star formation.

Contemporaneous high-angular-resolution imaging of the AGB star W Hya in vibrationally excited H2O lines and visible polarized light with ALMA and VLT/SPHERE-ZIMPOL

We present contemporaneous high-angular-resolution millimeter imaging and visible polarimetric imaging of the nearby asymptotic giant branch (AGB) star W Hya to better understand the dynamics and dust formation within a few stellar radii. The star W Hya \ was observed in two vibrationally excited H$_2$O lines at 268 and 251 GHz with Atacama Large Millimeter/submillimeter Array (ALMA) at a spatial resolution of 16times 20 mas and at 748 and 820 nm at a resolution of 26times 27 mas with the Very Large Telescope (VLT)/Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE)-Zurich Imaging Polarimeter (ZIMPOL). ALMA's high spatial resolution allowed us to image strong emission of the vibrationally excited H$_2$O \ line at 268 GHz ($ K_a,K_c $) over the stellar surface instead of absorption against the continuum, which is expected for thermal excitation. Strong, spotty emission was also detected along and just outside the stellar disk limb at an angular distance of sim 40 mas (sim 1.9 $R_ star ), extending to sim 60 mas (sim 2.9 $R_ star ). Another H$_2$O \ line ($ K_a,K_c $) at 251 GHz with a similar upper-level energy was tentatively identified, which shows absorption over the stellar surface. This suggests that the emission over the surface seen in the 268 GHz H$_2$O \ line is suprathermal or even maser emission. The estimated gas temperature and H$_2$O \ density are consistent with the radiatively pumped masers. The 268 GHz H$_2$O \ line reveals global infall at up to sim 15 km s$^ within 2--3 $R_ star but outflows at up to sim 8 km s$^ \ are also present. The polarized intensity maps obtained in the visible reveal clumpy dust clouds forming within sim 40 mas (sim 1.9 $R_ star ) with a particularly prominent cloud in the SW quadrant and a weaker cloud in the east. The 268 GHz H$_2$O \ emission overlaps very well with the visible polarized intensity maps, which suggests that both the nonthermal and likely maser H$_2$O \ emission and the dust originate from dense, cool pockets in the inhomogeneous atmosphere within sim 2--3 $R_ star

Biases in Exoplanet Transmission Spectra Introduced by Limb-darkening Parametrization

One of the main endeavors of the field of exoplanetary sciences is the characterization of exoplanet atmospheres on a population level. The current method of choice to accomplish this task is transmission spectroscopy, where the apparent radius of a transiting exoplanet is measured at multiple wavelengths in search of atomic and molecular absorption features produced by the upper atmosphere constituents. To extract the planetary radius from a transit light curve, it is necessary to account for the decrease in luminosity away from the center of the projected stellar disk, known as limb darkening. Physically motivated parametrizations of limb darkening, in particular of the quadratic form, are commonly used in exoplanet transit light-curve fitting. Here, we show that such parametrizations can introduce significant wavelength-dependent biases in the transmission spectra currently obtained with all instrument modes of the JWST, and thus have the potential to affect atmospheric inferences. To avoid such biases, we recommend the use of standard limb-darkening parametrizations with wide uninformative priors that allow for nonphysical stellar intensity profiles in the transit fits, and thus for a complete and symmetrical exploration of the parameter space. We further find that fitting the light curves at the native resolution results in errors on the measured transit depths that are significantly smaller compared to light curves that are binned in wavelength before fitting, thus potentially maximizing the amount of information that can be extracted from the data.

Determining the Extents, Geometries, and Kinematics of Narrow-line Region Outflows in Nearby Seyfert Galaxies

Outflowing gas from supermassive black holes in the centers of active galaxies has been postulated as a major contributor to galactic evolution. To explore the interaction between narrow-line region (NLR) outflows and their host galaxies, we use Hubble Space Telescope (HST) Space Telescope Imaging Spectrograph (STIS) spectra and Wide Field Camera 3 (WFC3) images of 15 nearby (z < 0.02) active galactic nuclei (AGN) to determine the extents and geometries of their NLRs. We combine new HST WFC3 continuum and [O iii] λ5007 images of 11 AGN with four archival AGN to match existing spectra from HST STIS. For the six AGN with suitable long-slit coverage of their NLRs, we use isophotal fitting of ground-based images, continuum-subtracted [O iii] images, and the STIS spectra, to resolve, measure, and de-project the gas kinematics to the plane of the host galaxy disk and distinguish NLR outflows from galaxy rotation and/or kinematically disturbed gas. We find an average [O iii] extent of ∼680 pc with a correlation between gas extent and [O iii] luminosity of R [O III] ∝ L[OIII]0.39 . The measured extents depend strongly on the depth of the [O iii] images, highlighting the importance of adopting uniform thresholds when analyzing scaling relationships. The outflows reach from 39% to 88% of the full NLR extents, and we find that all six of the AGN with STIS coverage of their entire NLRs show strong kinematic evidence for outflows, despite previous uncertainty for these AGN. This suggests that NLR outflows are ubiquitous in moderate-luminosity AGN and that standard criteria for kinematic modeling are essential for identifying outflows.

Collisions Between Dark Matter Confined High Velocity Clouds and Magnetized Galactic Disks: The Smith Cloud. II.

We present high-resolution simulations of high velocity clouds (HVCs) colliding with the outer part of the Galactic disk. All of the simulations include a 3 × 108 M ⊙ dark matter subhalo. Three simulations model a dark matter subhalo without a gaseous component, while eight simulations model a dark matter subhalo accompanied by a gaseous cloud of mass 2–8 × 106 M ⊙. Half of the simulations include the coherent component of the Galaxy's magnetic field. Each simulation spans ∼40 million years before the collision and ∼40 million years after the collision. The collisions between the gas cloud and disk splash gas into the halo, punch half-kiloparsec-size holes in the disk, and form long-lived, multi-kiloparsec-size shells. Each shell encloses a bubble of relatively cool gas. Holes and shells of these scales would be observable, and some have been observed in the past. We determine the fate of the HVC gas, temperature, composition, and ionization state of the bubble and shell gas, size and longevity of the holes, and effects of cloud density. Simulations show that the clouds do not survive the chaos of passage through the disk, but instead become part of the splash, bubble, and shell. Some dark matter clouds may appear to carry material with them long after the collision, but this material is shell gas that was captured by the dark matter subhalo. These results have ramifications for the Smith Cloud and other clouds hypothesized to have hit the Galactic disk.

ISMGCC: Finding Gas Structures in Molecular Interstellar Medium Using Gaussian Decomposition and Graph Theory

Molecular line emissions are commonly used to trace the distribution and properties of molecular Interstellar Medium. However, the emissions are heavily blended on the Galactic disk toward the inner Galaxy because of the relatively large line widths and the velocity overlaps of spiral arms. Structure identification methods based on voxel connectivity in Position-Position-Velocity (PPV) data cubes often produce unrealistically large structures, which is the “over-linking” problem. Therefore, identifying molecular cloud structures in these directions is not trivial. We propose a new method based on Gaussian decomposition and graph theory to solve the over-linking problem, named InterStellar Medium Gaussian Component Clustering (ISMGCC). Using the Milky Way Imaging Scroll Painting (MWISP) 13CO(1–0) data in the range of 13.°5 ≤ l ≤ 14.°5, ∣b∣ ≤ 0.°5, and −100 ≤ V lsr ≤ +200 km s−1, our method identified three hundred molecular gas structures with at least 16 pixels. These structures contain 92% of the total flux in the raw data cube and show single-peaked line profiles on more than 93% of their pixels. The ISMGCC method could distinguish gas structures in crowded regions and retain most of the flux without global data clipping or assumptions on the structure geometry, meanwhile, allowing multiple Gaussian components for complicated line profiles.

A Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE)

We examine the prevalence of truncated star-forming discs in the Virgo cluster down to M* ≃ 107 M⊙. This work makes use of deep, high-resolution imaging in the Hα+[N II] narrow-band from the Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE) and optical imaging from the Next Generation Virgo Survey (NGVS). To aid in the understanding of the effects of the cluster environment on star formation in Virgo galaxies, we take a physically motivated approach to define the edge of the star-forming disc via a drop-off in the radial specific star formation rate profile. A comparison with the expected sizes of normal galactic discs provides a measure of how truncated star-forming discs are in the cluster. We find that truncated star-forming discs are nearly ubiquitous across all regions of the Virgo cluster, including beyond the virial radius (0.974 Mpc). The majority of truncated discs at large cluster-centric radii are of galaxies likely on their first infall. As the intra-cluster medium density is low in this region, it is difficult to explain this population with solely ram-pressure stripping. A plausible explanation is that these galaxies are undergoing starvation of their gas supply before ram-pressure stripping becomes the dominant quenching mechanism. A simple model of starvation shows that this mechanism can produce moderate disc truncations within 1−2 Gyr. This model is consistent with ‘slow-then-rapid’ or ‘delayed-then-rapid’ quenching, whereby the early starvation mode drives disc truncations without significant change to the integrated star formation rate, and the later ram-pressure stripping mode rapidly quenches the galaxy. The origin of starvation may be in the group structures that exist around the main Virgo cluster, which indicates the importance of understanding pre-processing of galaxies beyond the cluster virial radius.

A Possible Formation Scenario of the Gaia BH1: Inner Binary Merger in Triple Systems

Based on astrometric measurements and spectral analysis from Gaia DR3, two quiescent black hole (BH) binaries, Gaia BH1 and BH2, have been identified. Their origins remain controversial, particularly for Gaia BH1. By considering a rapidly rotating (ω/ω crit = 0.8) and strongly magnetized (B 0 = 5000 G) merger product, we find that, at typical Galactic metallicity, the merger product can undergo efficient chemically homogeneous evolution. This results in the merger product having a significantly smaller radius during its evolution compared to that of a normally evolving massive star. Under the condition that the initial triple stability is satisfied, we use the Multiple Stellar Evolution code and the MESA code to identify an initial hierarchical triple that can evolve into Gaia BH1. It initially consists of three stars with masses of 9.03 M ⊙, 3.12 M ⊙, and 1 M ⊙, with inner and outer orbital periods of 2.21 days and 121.92 days, and inner and outer eccentricities of 0.41 and 0.45, respectively. This triple initially experiences triple evolution dynamics instability (TEDI) followed by Roche lobe overflow (RLOF). During RLOF, the inner orbit shrinks, and tidal effects gradually suppress the TEDI. Eventually, the inner binary undergoes a merger through contact (or collision). Finally, using models of rapidly rotating and strongly magnetic stars, along with standard core-collapse supernova (SN) or failed supernova (FSN) models, we find that a postmerger binary (PMB) consisting of an 12.11 M ⊙ merger product and a 1 M ⊙ companion star (originally an outer tertiary) can avoid RLOF. After an SN or FSN with a low ejected mass of ∼0.22 M ⊙ and a low kick velocity ( 46−33+25kms−1 or 9−8+16kms−1 ), the PMB can form Gaia BH1 in the Galactic disk.

CLOUDY modeling suggests a diversity of ionization mechanisms for diffuse extraplanar gas

Context. The ionization of diffuse gas located far above the energetic midplane OB stars poses a challenge to the commonly accepted notion that radiation from OB stars is the primary ionization source for gas in galaxies. Aims. We investigated the sources of ionizing radiation, specifically leaking midplane H II regions and/or in situ hot low-mass evolved stars (HOLMES), in extraplanar diffuse ionized gas (eDIG) in a sample of eight nearby (17−52 Mpc) edge-on disk galaxies observed with the Multi Unit Spectroscopic Explorer (MUSE). Methods. We constructed a model for the photoionization of eDIG clouds and the propagation of ionizing radiation through the eDIG using subsequent runs of CLOUDY photoionization code. Our model includes radiation originating both from midplane OB stars and in situ evolved stars and its dilution and processing as it propagates in the eDIG. Results. We fit the model to the data using the vertical line ratio profiles of our sample galaxies, and find that while the ionization by in situ evolved stars is insignificant for most of the galaxies in our sample, it may be able to explain the enhanced high-ionization lines in the eDIG of the green valley galaxy ESO 544−27. Conclusions. Our results show that while leaking radiation from midplane H II regions is the primary ionization source for eDIG, in situ evolved stars can play a significant part in ionizing extraplanar gas in galaxies with low star forming rates.

An Edge-on Regular Disk Galaxy at z = 5.289

While rotation-supported gas disks are known to exist as early as at z ≈ 7, it is still a general belief that stellar disks form late in the Universe. This picture is now being challenged by the observations from the James Webb Space Telescope (JWST), which have revealed a large number of disk-like galaxies that could be at z > 3, with some being candidates at z > 7. As an early formation of stellar disks will greatly impact our theory of galaxy formation and evolution, it is important to determine when such systems first emerged. Here we present D-CEERS-RUBIES-z5289 at z = 5.289 ± 0.001, the second confirmed stellar disk at z > 5, discovered using the archival JWST NIRCam imaging and NIRSpec spectroscopic data. This galaxy has a highly regular edge-on disk morphology, extends to ∼6.2 kpc along its major axis, and has an effective radius of ∼1.3–1.4 kpc. Such a large stellar disk is yet to be produced in numerical simulations. By analyzing its 10-band spectral energy distribution using four different tools, we find that it has a high stellar mass of 109.5–10.0 M ⊙. Its age is in the range of 330–510 Myr, and it has a mild star formation rate of 10–30 M ⊙ yr−1. While the current spectroscopic data do not allow the derivation of its rotation curve, the width of its Hα line from the partial slit coverage on one side of the disk reaches ∼345 km s−1, which suggests that it could have a significant contribution from rotation.

On the Origin of the 107 K Hot Emitting Gas in the Circumgalactic Medium of the Milky Way

The presence of the ≈106 K gas in the circumgalactic medium of the Milky Way (MW) has been well established. However, the location and the origin of the newly discovered hot gas at “supervirial (SV)” temperatures of ≈107 K have been puzzling. This hot gas has been detected in both absorption and emission; here, we focus on the emitting gas only. We show that both the “virial” and the SV temperature gas, as observed in emission, occupy disk-like extraplanar regions, in addition to the diffuse virial temperature gas filling the halo of the MW. We perform idealized hydrodynamical simulations to show that the ≈107 K emitting gas is likely to be produced by stellar feedback in and around the Galactic disk. We further show that the emitting gas at both SV and virial temperatures in the extraplanar regions is metal enriched and is not in hydrostatic equilibrium with the halo but is continuously evolving.

Evaluating the gravitational wave detectability of globular clusters and the Magellanic Clouds for LISA

We used the stellar evolution code bpass and the gravitational wave (GW) simulation code legwork to simulate populations of compact binaries that may be detected by the future space-based GW detector LISA. Specifically, we simulate the Magellanic Clouds and binary populations mimicking several globular clusters, neglecting dynamical effects. We find that a handful of sources should be detectable in each of the Magellanic Clouds, but for globular clusters the amount of detectable sources will likely be less than one each. We compared our results to earlier research and find that our predicted numbers are several dozen times lower than both the results from calculations that used the stellar evolution code bse and take dynamical effects into account, and results from calculations that used the stellar evolution code SeBa for the Magellanic Clouds. Earlier research that compared bpass models for GW sources in the Galactic disk with bse models found a similarly sized discrepancy. We determine that this discrepancy is caused by differences between the stellar evolution codes, particularly in the treatment of mass transfer and common-envelope events in binaries: in bpass mass transfer is more likely to be stable and tends to lead to less orbital shrinkage in the common-envelope phase than in other codes. This difference results in fewer compact binaries with periods short enough to be detected by LISA existing in the bpass population. For globular clusters, we conclude that the impact of dynamical effects is uncertain based on the literature, but the differences in stellar evolution have an effect of a factor of 20 to 40 on the number of detectable binaries.