Search NASA⌕ Search

SEARCH · Search NASA

Results for “stellar evolution”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Simulations of Interacting Binary Systems—Pathways to Radio-bright Gamma-Ray Burst Progenitors

Although the association of gamma-ray bursts (GRBs) with massive stellar death is on firm footing, the nature of the progenitor system and the key ingredients required for a massive star to produce a GRB remain open questions. Here, we investigate the evolution of a 15–25M ⊙ massive star with a 10–15M ⊙ black hole (BH) using Modules for Experiments in Stellar Astrophysics. We quantify companion-influenced angular-momentum evolution over a stellar lifetime for orbital periods where tides are significant, varying stellar and BH masses, initial stellar spin, and accretion and dynamo prescriptions while tracking mass loss and angular momentum. The final spin is set by tidal torques versus stellar winds. For binaries that initially avoid Roche-lobe overflow, tides can spin up the star, but a late-stage expansion can drive tidal stripping; associated mass and angular-momentum loss can suppress spin-up. We find that massive-star BH binaries at comparable mass ratios may be potential GRB progenitors for short orbital periods (∼20–5 × 10 2 days) and long orbital periods (∼2 × 10 3 –4 × 10 3 days), where our suite of lifetime simulations reveals a favored parameter space with negligible mass loss and enough spin angular momentum to power a GRB jet. For initially nonrotating stars, this provides a lower limit on the final spin above a threshold estimate consistent with forming a postcollapse BH mass of 5–10M ⊙ with a spin parameter ≥ 0.5. For initially rapidly rotating stars, tidal interactions may sustain high spin when mass loss is negligible because the binary is not tidally synchronized.

Hernandez, Angel [Univ. of Colorado, Boulder, CO (↗
DOE OSTI2025

Improving galaxy cluster selection with the outskirt stellar mass of galaxies

The number density and redshift evolution of optically selected galaxy clusters offer an independent measurement of the amplitude of matter fluctuations, 𝑆 8 . However, recent results have shown that clusters chosen by the redMaPPer algorithm show richness-dependent biases that affect the weak lensing signals and number densities of clusters, increasing uncertainty in the cluster mass calibration and reducing their constraining power. Here, in this work, we evaluate an alternative cluster proxy, outskirt stellar mass, 𝑀 out , defined as the total stellar mass within a [50, 100] kpc envelope centered on a massive galaxy. This proxy exhibits scatter comparable to redMaPPer richness, 𝜆, but is less likely to be subject to projection effects. We compare the Dark Energy Survey Year 3 redMaPPer cluster catalog with a 𝑀 out selected cluster sample from the Hyper-Suprime Camera survey. We use weak lensing measurements to quantify and compare the scatter of 𝑀 out and 𝜆 with halo mass. Our results show 𝑀 out has a scatter consistent with 𝜆, with a similar halo mass dependence, and that both proxies contain unique information about the underlying halo mass. We find 𝜆-selected samples introduce features into the measured Δ⁢Σ signal that are not well fit by a log-normal scatter only model, absent in 𝑀 out selected samples. Our findings suggest that 𝑀 out offers an alternative for cluster selection with more easily calibrated selection biases, at least at the generally lower richnesses probed here. Combining both proxies may yield a mass proxy with a lower scatter and more tractable selection biases, enabling the use of lower mass clusters in cosmology. Finally, we find the scatter and slope in the 𝜆 −𝑀 out scaling relation to be 0.49 ±0.02 and 0.38 ±0.09.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

On the Importance of the Convective Urca Process in 3D Simulations of a Simmering White Dwarf

Type Ia supernovae are bright thermonuclear explosions that are important to numerous areas of astronomy. However, the origins of these events are poorly understood. One proposed setting is that of a near Chandrasekhar mass white dwarf that undergoes runaway carbon burning in the core. During the thousand years leading up to the explosion, the white dwarf undergoes a simmering phase where slow carbon burning heats the core and drives convection. A poorly understood aspect of this phase is the convective Urca process, which links convection with weak nuclear reactions. We use the low Mach number code MAESTROeX to perform full 3D simulations as is required to accurately capture the turbulent convection. We present simulations with and without the A=23 convective Urca process, which have relaxed to a steady state. We characterize the effects of the convective Urca process on the neutrino losses, the nuclear energy generation, and the convective boundary. We find that the size of the convection zone is substantially reduced by the convective Urca process, though convection still extends past the Urca shell. Our findings on the structure of the convective zone and the compositional changes can be used to inform 1D stellar models that track the longer-timescale evolution.

FOS: Physical sciences↗
DOE OSTI2026

Proton-rich Production of Lanthanides: The vi Process

The astrophysical origin of the lanthanides is an open question in nuclear astrophysics. Besides the widely studied s, i, and r processes in moderately to strongly neutron-rich environments, an intriguing alternative site for lanthanide production could in fact be robustly proton-rich matter outflows from core-collapse supernovae under specific conditions—in particular, high-entropy winds with enhanced neutrino luminosity and fast dynamical timescales. In this environment, excess protons present after charged-particle reactions have ceased can continue to be converted to neutrons by (anti)neutrino interactions, producing a neutron-capture reaction flow up to A ∼ 200. This scenario, christened the νi process in a recent paper, has previously been discussed as a possibility. Here, we examine the prospects for the νi process through the lenses of stellar abundance patterns, bolometric light curves, and galactic chemical evolution models, with a particular focus on hypernovae as candidate sites. We identify specific lanthanide signatures for which the νi process can provide a credible supplement to the r/i processes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Fast explicit solutions for neutrino-electron scattering: Explicit asymptotic methods

Here, we present results of explicit asymptotic approximations applied to neutrino-electron scattering in a representative model of neutrino population evolution under conditions characteristic of core-collapse supernova explosions or binary neutron star mergers. It is shown that this approach provides stable solutions of these stiff systems of equations, with accuracy and time stepping comparable to that for standard implicit treatments such as backward Euler, fixed point iteration, and Anderson-accelerated fixed point iteration. Because each time step can be computed more rapidly with the explicit asymptotic approximation than with implicit methods, this suggests that algebraically stabilized explicit integration methods could be used to compute neutrino evolution coupled to hydrodynamics more efficiently in stellar explosions and mergers than the methods currently in use.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

16 O ⁡(𝑒,𝑒′⁢𝛼)⁢ 12 C measurements and the 12 C ⁡(𝛼,𝛾)⁢ 16 O astrophysical reaction rate

The 12 C ⁡(𝛼,𝛾)⁢ 16 O reaction, an important component of stellar helium burning, has a key role in nuclear astrophysics. It has significant impact on the evolution and final state of heavy to low mass stars, maximum mass of stellar formed black holes and also shapes the elemental abundances resulting from nucleosynthesis in such stars. Providing a reliable estimate for the energy dependence of this reaction at stellar helium burning temperatures has been a longstanding and important goal. In this work, we study the role of potential new 𝐸⁢1 and 𝐸⁢2 measurements of the 16 O ⁡(𝑒,𝑒′⁢𝛼)⁢ 12 C reaction in reducing the overall uncertainty in the astrophysical 𝑆 factor for the 12 C ⁡(𝛼,𝛾)⁢ 16 O 𝐸⁢1 and 𝐸⁢2 ground state capture extrapolated to a stellar energy of 300 keV. A multilevel 𝑅-matrix analysis is used to make extrapolations of the 𝑆 𝐸⁢1 (300 keV) and 𝑆 𝐸⁢2 (300 keV) factors for the 12 C ⁡(𝛼,𝛾)⁢ 16 O reaction from existing ground state capture data. Bayesian analysis is used to quantify the uncertainties in the extrapolations for both the existing data alone and also when possible new experimental data are included. In particular, we consider a new experiment that would make use of a high-intensity low-energy electron beam that impinges on a windowless oxygen gas target as a means to determine the total 𝐸⁢1 and 𝐸⁢2 ground state cross sections for this reaction. Here we find that the new data could significantly reduce the 𝑆⁡(300) uncertainties. Splitting the new data into high and low energy regions shows that both low and high energy data are effective in reducing the uncertainty.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Accurately simulating core-collapse self-interacting dark matter halos

The properties of satellite halos provide a promising probe for dark matter (DM) physics. Observations have motivated current efforts to explain surprisingly compact DM halos. If DM is not collisionless, but has strong self-interactions, halos can undergo gravothermal collapse, leading to higher densities in the central region of the halo. However, it is challenging to model this collapse phase from first principles. To improve on this, we sought to better understand the numerical challenges and convergence properties of self-interacting dark matter (SIDM) N-body simulations in the collapse phase. Especially, our aim was to better understand the evolution of satellite halos. To do so, we ran SIDM N-body simulations of a low-mass halo in isolation and within an external gravitational potential. The simulation set-up was motivated by the perturber of the stellar stream GD-1. We find that the halo evolution is very sensitive to energy conservation errors, and a SIDM kernel size that is too large can artificially speed up the collapse. Moreover, we demonstrate that the King model can describe the density profile at small radii for the late stages that we have simulated. Furthermore, for our most highly resolved simulation (N = 5 × 10 7 ) we have made the data public. It can serve as a benchmark. Overall, we find that the current numerical methods do not suffer from convergence problems in the late collapse phase and provide guidance on how to choose numerical parameters, for example that the energy conservation error is better kept well below 1%. This allows simulations to be run of halos that become concentrated enough to explain observations of GD-1-like stellar streams or strong gravitational lensing systems.

dark matter↗
DOE OSTI2025

Effects of Dark Matter Self-interactions on Sagittarius and its Stream

This work explores how assumptions regarding the particle-physics nature of dark matter can alter the evolution of the Sagittarius (Sgr) dwarf spheroidal galaxy and its expansive stellar stream. We run a large suite of N-body simulations to model the infall of a spherically symmetric Sgr-like dwarf, exploring how the presence of dark matter self-interactions impacts its evolution. For a scattering cross section of σ/m χ = 30 cm 2 g −1 (at orbital velocity scales), these interactions result in significantly less stellar mass and little to no dark matter bound to the progenitor at the present day. To isolate the cause of this mass loss, we introduce a novel technique for controlling which pairs of dark matter simulation particles can interact. This enables us to identify ram-pressure evaporation—the scattering of satellite and host dark matter particles—as the primary source of the enhanced mass loss. The rapid disintegration of the Sgr progenitor when self-interactions are allowed alters some key properties of the resulting stellar stream, most dramatically suppressing the presence of a “spur” on the apocenter of the trailing stream arm that correlates with the mass of the satellite at last pericenter. We demonstrate how the effects on the Sgr system scale with the particular choice of self-interaction cross section, which affects the degree of ram-pressure evaporation. These findings generalize beyond the Sgr system, underscoring that dwarf stellar streams and dwarf galaxies with close passages may serve as sensitive probes for dark matter self-interactions.

Milky Way dark matter halo↗
DOE OSTI2025

3D Convective Urca Process in a Simmering White Dwarf

Abstract A proposed setting for thermonuclear (Type Ia) supernovae is a white dwarf that has gained mass from a companion to the point of carbon ignition in the core. In the early stages of carbon burning, called the simmering phase, energy released by the reactions in the core drive the formation and growth of a core convection zone. One aspect of this phase is the convective Urca process, a linking of weak nuclear reactions to convection, which may alter the composition and structure of the white dwarf. The convective Urca process is not well understood and requires 3D fluid simulations to properly model the turbulent convection, an inherently 3D process. Because the neutron excess of the fluid both sets and is set by the extent of the convection zone, the realistic steady state can only be determined in simulations with real 3D mixing processes. Additionally, the convection is relatively slow (Mach number less than 0.005) and thus a low Mach number method is needed to model the flow over many convective turnovers. Using the MAESTROeX low Mach number hydrodynamic software, we present the first full-star 3D simulations of the A = 23 convective Urca process, spanning hundreds of convective turnover times. Our findings on the extent of mixing across the Urca shell, the characteristic velocities of the flow, the energy-loss rates due to neutrino emission, and the structure of the convective boundary can be used to inform 1D stellar models that track the longer-timescale evolution.

Boyd, Brendan (ORCID:0000000254199751)↗
DOE OSTI2025

Modeling resistive-inductive evolution of currents in Wendelstein 7-X

This research investigates the temporal evolution of the toroidal plasma current in the Wendelstein 7-X (W7-X) stellarator under different heating, fueling, and current drive scenarios. The THRIFT code has been modernized and its predictions of the evolution of the toroidal current have been compared against experimentally measured currents in W7-X. Good agreement is found with respect to the characteristic timescale between experimentally measured and simulated toroidal currents. The total bootstrap current is under-predicted owing to the applicability of the BOOTSJ model for the plasma collisionalities in question. Edge plasma resistivity is found to play an important role in the asymptotic behavior of the evolution of the current, indicating a possible limitation of the minimum plasma temperature when applying this model. Simulations of ECCD and heating power steps show THRIFT is capable of capturing the dynamical evolution of the current in response to changes in current sources. Future uses of THRIFT include validating and benchmarking other non-inductive current models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Formation and Redshift Evolution of Dark Matter Spikes

Dark matter density spikes forming around adiabatically growing black holes can dramatically enhance indirect and direct detection signals. Canonical predictions, however, assume a zero-mass seed in a purely dark matter environment and do not track the long-term dynamical impact of surrounding stars. We present a semi-analytic framework that first generalizes adiabatic spike formation to include finite seed masses, stellar cusps, and non-circular orbits, and then studies the subsequent cosmic evolution by solving coupled Fokker-Planck equations for the dark matter and stellar phase-space distributions, with a heating rate modulated by the cosmic star formation rate. Starting conservatively from canonical Gondolo-Silk spikes and marginalizing over astrophysical uncertainties, we find that stellar gravitational heating drives the inner slope towards $γ_χ\simeq 1.5$ within a few Gyrs (e.g by $z \lesssim 2$ for spikes formed at $z\simeq 10$), yielding overdensities two to four orders of magnitude below canonical expectations but still well above an NFW-like cusp. We provide redshift-dependent benchmarks for the column density and $J$-factor relevant to scattering, decay and annihilation signatures. Any robust interpretation of indirect dark matter signals from galactic nuclei must account for this evolution.

Herrera, Gonzalo [MIT, LNS; MIT, MKI; Harvard U.]↗
DOE OSTI2024

Investigating nonlinear magnetohydrodynamics in an optimized, reactor-scale quasi-axisymmetric stellarator

We use the new simulation capabilities of the extended-magnetohydrodynamic (MHD) code, M3D-C1, to investigate the nonlinear MHD properties of a reactor-scale quasisymmetric stellarator equilibrium. Our model captures the self-consistent evolution of the magnetic field, temperature, density, and flow profiles without imposing restrictions on the structure of the first. We include the effects of resistivity using a realistic temperature-dependent Spitzer model, along with a model for heat transport that captures the key physical characteristic, namely, strongly anisotropic diffusion in directions perpendicular and parallel to the magnetic field. We consider a quasi-axisymmetric, finite-pressure equilibrium that was optimized for self-consistent bootstrap current, quasi-symmetry, and energetic particle confinement. Our assessment finds that the equilibrium is highly unstable to interchange-like pressure-driven instabilities near the plasma edge. The initially unstable modes rapidly destabilize other modes in the direction of the N-fold rotational symmetry (toroidal, in this case). For this equilibrium, N = 2, meaning destabilization of a large number of even-numbered toroidal Fourier modes. Thus, field-periodicity is likely to be an important factor in the nonlinear MHD stability characteristics of optimized stellarators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Beyond Mass and Multiscale Environments: What Shapes Low Surface Brightness Galaxies? Evidence from MaNGA

The origin of low surface brightness (LSB) galaxies remains a key open question in galaxy formation, reflecting the balance between internal mechanisms and environmental influence. Using MaNGA integral-field spectroscopy, we investigate whether LSB and high surface brightness (HSB) galaxies of comparable stellar mass (9 < logM * < 10) occupy distinct environments or differ primarily through internal evolution. Our late-type sample comprises 113 central and 29 satellite LSB galaxies, and 374 central and 142 satellite HSB galaxies. We characterize environments on scales from ∼100 kpc to 10 Mpc, analyzing radial profiles of stellar mass surface density (Σ * ), star formation activity, and gas-phase metallicity. Central LSB and HSB galaxies inhabit similarly low-density large-scale (>200 kpc) environments, but LSB galaxies are more isolated on small scales (∼100 kpc). Even after matching in stellar mass and environment, LSB galaxies show systematically lower Σ * , Σ SFR , and metallicities, often hosting diffuse, weakly star-forming bulges embedded in extended disks. These results indicate that LSB structure and star formation are not primarily governed by a large-scale environment or halo mass. While secondary halo properties such as spin, concentration, or gas accretion history are often invoked, their environmental dependence appears weak. Instead, LSB–HSB differences for centrals likely reflect divergent assembly or interaction histories and internal processes—such as angular momentum-driven disk evolution or inefficient gas conversion—largely decoupled from a large-scale environment. Nonetheless, the environment still influences the observed star formation and chemical differences between central and satellite LSB galaxies.

Shen, Mengting [Xiamen University (China); SDSS Co↗
DOE OSTI2024

The Evolution of Galaxies and Clusters at High Spatial Resolution with Advanced X-ray Imaging Satellite (AXIS)

Stellar and black hole feedback heat and disperse surrounding cold gas clouds, launching gas flows off circumnuclear and galactic disks, producing a dynamic interstellar medium. On large scales bordering the cosmic web, feedback drives enriched gas out of galaxies and groups, seeding the intergalactic medium with heavy elements. In this way, feedback shapes galaxy evolution by shutting down star formation and ultimately curtailing the growth of structure after the peak at redshift 2–3. To understand the complex interplay between gravity and feedback, we must resolve both the key physics within galaxies and map the impact of these processes over large scales, out into the cosmic web. The Advanced X-ray Imaging Satellite (AXIS) is a proposed X-ray probe mission for the 2030s with arcsecond spatial resolution, large effective area, and low background. AXIS will untangle the interactions of winds, radiation, jets, and supernovae with the surrounding interstellar medium across the wide range of mass scales and large volumes driving galaxy evolution and trace the establishment of feedback back to the main event at cosmic noon. This white paper is part of a series commissioned for the AXIS Probe mission concept; additional AXIS white papers can be found at the AXIS website.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Magnetosphere Evolution and Precursor-driven Electromagnetic Signals in Merging Binary Neutron Stars

We detail new force-free simulations to investigate magnetosphere evolution and precursor electromagnetic (EM) signals from binary neutron stars. Our simulations fully follow a representative inspiral motion, capturing the intricate magnetospheric dynamics and their impact on EM outflows. We explore a range of stellar magnetic moment orientations and relative strengths, finding that the magnetospheres and Poynting flux evolution are strongly configuration dependent. The Poynting flux exhibits pulsations at twice the orbital frequency, 2Ω, and is highly anisotropic, following a power-law dependence on orbital frequency. The index ranges from 1 to 6, shaped by the intricate dynamics of the magnetospheres. Furthermore, we present the first computation of (1) the EM forces acting on the star surfaces, revealing the presence of torques that, for highly magnetized stars, could influence the orbital dynamics or break the crust; (2) the high-energy emission signals from these systems by adopting the established isolated pulsar theory. Assuming curvature radiation in the radiation-reaction limit, we find that photons could reach TeV–PeV energies in the last ∼ms for magnetic field strengths 10 10 –10 15 G. However, our analysis of single photon magnetic pair production suggests that these photons are unlikely to escape, with the MeV band emerging as a promising observational window for precursor high-energy emission. In this framework, we construct high-energy emission skymaps and light curves, exploring observational implications. Finally, we propose potential precursor radio emission and delayed afterglow echoes from magnetized outflows, which may contribute to late-time rebrightening in short gamma-ray bursts or to orphan afterglows.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The DESI One-Percent Survey: Modelling the clustering and halo occupation of all four DESI tracers with U CHUU

We present results from a set of mock lightcones for the DESI One-Percent Survey, created from the UCHUU simulation. This 8 h −3 Gpc 3 N-body simulation comprises 2.1 trillion particles and provides high-resolution dark matter (sub)haloes in the framework of the Planck-based ΛCDM cosmology. Employing the subhalo abundance matching (SHAM) technique, we populated the UCHUU (sub)haloes with all four DESI tracers – Bright Galaxy Survey (BGS), luminous red galaxies (LRGs), emission line galaxies (ELGs), and quasars (QSOs) – to z = 2.1. Our method accounts for redshift evolution as well as the clustering dependence on luminosity and stellar mass. The two-point clustering statistics of the DESI One-Percent Survey generally agree with predictions from UCHUU across scales ranging from 0.3 h −1 Mpc to 100 h −1 Mpc for the BGS and across scales ranging from 5 h −1 Mpc to 100 h −1 Mpc for the other tracers. We observed some differences in clustering statistics that can be attributed to incompleteness of the massive end of the stellar mass function of LRGs, our use of a simplified galaxy-halo connection model for ELGs and QSOs, and cosmic variance. We find that at the high precision of UCHUU, the shape of the halo occupation distribution (HOD) of the BGS and LRG samples is smaller bias values, likely due to cosmic variance. The bias dependence on absolute magnitude, stellar mass, and redshift aligns with that of previous surveys. These results provide DESI with tools to generate high-fidelity lightcones for the remainder of the survey and enhance our understanding of the galaxy-halo connection.

cosmology↗
DOE OSTI2026

Neural operator transformers capture bifurcating drift-wave turbulence in fusion plasma simulations

Self-consistent modeling of turbulence-driven transport is critical for optimizing confinement in magnetically confined fusion plasmas, such as tokamaks and stellarators. In particular, capturing the long-term co-evolution of turbulence, flow, and background plasma profiles remains computationally challenging. Direct numerical simulation of these multiscale, highly nonlinear processes is often demanding and impractical for real-time control or design optimization. To address this bottleneck, we investigate transformer-based neural operator partial differential equation surrogates for emulating the dynamics of drift-wave turbulence bifurcation mediated by zonal flows, using the modified Hasegawa–Wakatani (MHW) model as a prototypical system. We find that the finetuned neural operator model has excellent performance in capturing the multi-spatiotemporal-scales of MHW turbulence bifurcation and is robust to testing on rare and out-of-distribution dynamics. Specifically, we demonstrate that a single unified model accurately predicts both quasi-steady-state turbulence and a wide range of dynamical transition processes, such as nonlinear saturation, spontaneous suppression of turbulence, and the emergence of macroscopic zonal flows, over time horizons vastly exceeding the local turbulence correlation time. This computationally efficient approach establishes a strong foundation for fast, AI-based modeling of complex, multiscale phenomena in magnetized fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Repainting the colour–mass diagrams by unearthing the green mountain: dust-rich S0 galaxies in the colour–(galaxy stellar mass) diagram, and the colour–(black hole mass) relations for dust-poor versus dust-rich galaxies

ABSTRACT Lenticular galaxies are notoriously misclassified as elliptical galaxies and, as such, a (disc inclination)-dependent correction for dust is often not applied to the magnitudes of dusty lenticular galaxies. This results in overly red galaxy colours, impacting their distribution in the colour–magnitude diagram. It is revealed how this has led to an underpopulation of the ‘green valley’ by hiding a ‘green mountain’ of massive dust-rich lenticular galaxies – known to be built from gas-rich major mergers – within the ‘red sequence’ of colour–(stellar mass) diagrams. Correcting for dust, a ‘green mountain’ appears at M*,gal ∼ 1011 M⊙, along with signs of an extension to lower masses producing a ‘green range’ or ‘green ridge’ on the green side of the ‘red sequence’ and ‘blue cloud.’ The ‘red sequence’ is shown to be comprised of two components: a red plateau defined by elliptical galaxies with a near-constant colour and by lower-mass dust-poor lenticular galaxies, which are mostly a primordial population but may include faded/transformed spiral galaxies. The presence of the quasi-triangular-shaped galaxy evolution sequence, previously called the ‘Triangal’, is revealed in the galaxy colour–(stellar mass) diagram. It tracks the speciation of galaxies and their associated migration through the diagram. The connection of the ‘Triangal’ to previous galaxy morphology sequences (Fork, Trident, and Comb) is also shown herein. Finally, the colour–(black hole mass) diagram is revisited, revealing how the dust correction generates a blue–green sequence for the spiral and dust-rich lenticular galaxies that is offset from a green–red sequence defined by the dust-poor lenticular and elliptical galaxies.

Graham, Alister W. (ORCID:0000000264969414)↗