SEARCH · Search NASA

Results for “triangularity”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,225 records · Page 68

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Toroidal modeling of runaway avalanche in DIII-D discharges

Abstract A toroidal modeling tool is developed to study the runaway electron (RE) avalanche production process in tokamak plasmas, by coupling the Rosenbluth–Putvinski avalanche model (Rosenbluth and Putvinski 1997 Nucl. Fusion 37 1355) with an n = 0 magneto-hydrodynamic (MHD) solver. Initial value numerical simulations are carried out for two DIII-D discharges with different plasma shapes (one near circular, and the other with high elongation). It is found that, assuming the same level of about 1% seed current level, the Rosenbluth–Putvinski model somewhat underestimates the RE plateau current for the circular-shaped plasma, as compared with that measured in DIII-D experiments. For an elongated, higher current plasma, simulations find strong runaway current avalanche production despite the lack of measured plateau RE current in experiments. A possible reason for this discrepancy is a lack of additional RE dissipation physics in the present two-dimensional model. Systematic scans of the plasma boundary shape, at fixed pre-disruption plasma current, find that the plasma elongation helps to reduce the RE avalanche production, confirming recent results obtained with an analytic model (Fülöp et al 2020 J. Plasma Phys. 86 474860101). The effect of the plasma triangularity (either positive or negative), on the other hand, has a minor effect. On the physics side, the avalanche process involves two competing mechanisms associated with the electric field. On the one hand, a stronger electric field produces a higher instantaneous avalanche growth rate. On the other hand, a fast growing RE current quickly reduces the fraction of the conduction current together with the electric field, and hence a faster dissipation of the poloidal flux. As a final result of these two dynamic processes, the runaway plateau current is not always the largest with the strongest initial electric field. These results lay the foundation for future self-consistent inclusion of the MHD dynamics and the RE amplification process.

Physics↗

DOE OSTI2021

Exploration of RMP ELM control on ITER-similar shape (ISS) in KSTAR

Recent control advances at KSTAR enabled us to not only establish but also to sustain the ITER-similar shape (ISS) in a stationary manner, matching the ITER-like dimensionless shape parameters in lower single null configuration with vastly contrasting upper/lower triangularities. Here, the kinetic parameters of the ISS are different from the typical KSTAR configurations seen in ELM-RMP suppression experiments. We experimentally demonstrate the ISS-compatible resonant magnetic perturbation (RMP) control in KSTAR for the first time, using the n = 2, +90 degrees phasing RMP.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

DIII-D research advancing the physics basis for optimizing the tokamak approach to fusion energy

DIII-D physics research addresses critical challenges for the operation of ITER and the next generation of fusion energy devices. This is done through a focus on innovations to provide solutions for high performance long pulse operation, coupled with fundamental plasma physics understanding and model validation, to drive scenario development by integrating high performance core and boundary plasmas. Substantial increases in off-axis current drive efficiency from an innovative top launch system for EC power, and in pressure broadening for Alfven eigenmode control from a co-/counter-I p steerable off-axis neutral beam, all improve the prospects for optimization of future long pulse/steady state high performance tokamak operation. Fundamental studies into the modes that drive the evolution of the pedestal pressure profile and electron vs ion heat flux validate predictive models of pedestal recovery after ELMs. Understanding the physics mechanisms of ELM control and density pumpout by 3D magnetic perturbation fields leads to confident predictions for ITER and future devices. Validated modeling of high-Z shattered pellet injection for disruption mitigation, runaway electron dissipation, and techniques for disruption prediction and avoidance including machine learning, give confidence in handling disruptivity for future devices. For the non-nuclear phase of ITER, two actuators are identified to lower the L–H threshold power in hydrogen plasmas. With this physics understanding and suite of capabilities, a high poloidal beta optimized-core scenario with an internal transport barrier that projects nearly to Q = 10 in ITER at ~8 MA was coupled to a detached divertor, and a near super H-mode optimized-pedestal scenario with co-I p beam injection was coupled to a radiative divertor. The hybrid core scenario was achieved directly, without the need for anomalous current diffusion, using off-axis current drive actuators. Also, a controller to assess proximity to stability limits and regulate β N in the ITER baseline scenario, based on plasma response to probing 3D fields, was demonstrated. Finally, innovative tokamak operation using a negative triangularity shape showed many attractive features for future pilot plant operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

The dependence of confinement on the isotope mass in the core and the edge of AUG and JET-ILW H-mode plasmas

Experiments in ASDEX Upgrade (AUG) and JET with the ITER-like wall (JET-ILW) are performed to separate the pedestal and core contributions to confinement in H-modes with different main ion masses. A strong isotope mass dependence in the pedestal is found which is enhanced at high gas puffing. This is because the ELM type changes when going from D to H for matched engineering parameters, which is likely due to differences in the inter ELM transport with isotope mass. The pedestal can be matched in H and D plasmas by varying only the triangularity and keeping the engineering parameters relevant for core transport the same. With matched pedestals Astra/TGLF (Sat1geo) core transport simulations predict the experimental profiles equally well for H and D. These core transport simulations show a negligible mass dependence and no gyro-Bohm scaling is observed. However, to match the experimental observations at medium β it is required to take the fast-ion dilution and rotation into account. This is not enough for high β plasmas where for the first time a profile match between H and D plasmas was achieved experimentally. Under these conditions quasilinear modelling with TGLF over predicts the transport in the core of H and D plasmas alike.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Experimental study on the role of the target electron temperature as a key parameter linking recycling to plasma performance in JET-ILW*

Changes in global and edge plasma parameters (H98(y,2), dimensionless collisionality ν*, core density peaking, separatrix density ne,sep) with variations in the D2 fueling rate and divertor configuration are unified into a single trend when mapped to $\langle$Te,ot$\rangle$, the spatially averaged spectroscopically derived outer target electron temperature. Dedicated JET with the ITER-like wall (JET-ILW) experiments in combination with an extended JET-ILW database of unseeded low-triangularity H-mode plasmas spanning a wide range of D2 fueling rates, Ip, Bt and heating power have demonstrated the importance of $\langle$Te,ot$\rangle$ as a key physics parameter linking the recycling particle source and detachment with plasma performance. The remarkably robust H98(y,2) trend with $\langle$Te,ot$\rangle$ is connected to a strong inverse correlation between $\langle$Te,ot$\rangle$, ne,sep and ν*, thus directly linking changes in the divertor recycling moderated by $\langle$Te,ot$\rangle$ with the previously established relationship between ν*, core density peaking and core pressure resulting in a degradation in core plasma performance with decreasing $\langle$Te,ot$\rangle$ (increasing ν*). Additionally, a strong inverse correlation between the separatrix to pedestal density ratio, ne,sep/ne,ped, and $\langle$Te,ot$\rangle$ is also established, with the rise in ne,sep/ne,ped saturating at $\langle$Te,ot$\rangle$ > 10 eV. A strong reduction in H98(y,2) is observed as $\langle$Te,ot$\rangle$ is driven from 30 to 10 eV via additional D2 gas fueling, while the divertor remains attached. Consequently, the pronounced performance degradation in attached divertor conditions has implications for impurity seeding radiative divertor scenarios, in which H98(y,2) is already low (~0.7) before impurities are injected into the plasma since moderate gas fueling rates are required to promote high divertor neutral pressure. A favorable pedestal pressure, pe,ped, dependence on Ip has also been observed, with an overall increase in pe,ped at Ip = 3.4 MA as $\langle$Te,ot$\rangle$ is driven down from attached to high-recycling divertor conditions. In contrast, pe,ped is reduced with decreasing $\langle$Te,ot$\rangle$ in the lower Ip branches. Further work is needed to (i) clarify the potential role of edge opacity on the observed favorable pedestal pressure Ip scaling; as well as to (ii) project the global and edge plasma performance trends with $\langle$Te,ot$\rangle$ to reactor-scale devices to improve predictive capability of the coupling between recycling and confined plasma fueling in what are foreseen to be more opaque edge plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Optimizing the HSX stellarator for microinstability by coil-current adjustments

Abstract The optimization of helically symmetric experiment (HSX) for reduced microinstability has been achieved by examining a large set of configurations within a neighborhood of the standard operating configuration. This entailed generating a database of more than 10 6 magnetic-field configurations for HSX by varying the currents in external coils. Using a set of volume-averaged metrics and gyrokinetic simulations, this database has helped to identify a set of configurations that can be used to regulate trapped-electron-mode stability in HSX. This set of configurations is also found to correlate flux-surface elongation and triangularity with an increase in magnetic-well depth, an increase in rotational transform, and low neoclassical heat-flux relative to the standard quasi-helically-symmetric configuration. These results demonstrate sensitivity of plasma behavior in response to changes in a 3D magnetic field to both neoclassical and gyrokinetic models, and the experimental potential in HSX to explore turbulence optimization. This perturbative optimization approach is not unique to HSX, and can readily be deployed on existing fusion devices to identify novel magnetic-fields to be used in turbulence-optimization experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Stability and transport of gyrokinetic critical pedestals

Abstract A gyrokinetic threshold model for pedestal width–height scaling prediction is applied to multiple devices. A shaping and aspect ratio scan is performed on National Spherical Torus Experiment (NSTX) equilibria, finding Δ ped = 0.92 A 1.04 κ − 1.24 0.38 δ β θ , ped 1.05 for the wide-pedestal branch with pedestal width Δ ped , aspect ratio A , elongation κ , triangularity δ , and normalized pedestal height β θ , ped . The width–transport scaling is found to vary significantly if the pedestal height is varied either with a fixed density or fixed temperature, showing how fueling and heating sources affect the pedestal density and temperature profiles for the kinetic-ballooning-mode (KBM) limited profiles. For an NSTX equilibrium, at fixed density, the wide branch is Δ ped = 0.028 ( q e / Γ e − 1.7 ) 1.5 ∼ η e 1.5 and at fixed temperature Δ ped = 0.31 ( q e / Γ e − 4.7 ) 0.85 ∼ η e 0.85 , where q e and Γ e are turbulent electron heat and particle fluxes and η e = ∇ ln ⁡ T e / ∇ ln ⁡ n e for an electron temperature T e and density n e . Pedestals close to the KBM limit are shown to have modified turbulent transport coefficients compared to the strongly driven KBMs. The role of flow shear is studied as a width–height scaling constraint and pedestal saturation mechanism for a standard and lithiated wide pedestal discharge. Finally, the stability, transport, and flow shear constraints are combined and examined for an NSTX experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Predicted thresholds for RMP ELM suppression access in double-null configurations

Modeling of DIII-D plasmas spanning shapes from single to double-null (DN) reveals new insights into the nature of resonant magnetic perturbation (RMP) conditions necessary for edge-localized mode (ELM) suppression. The suppression of ELMs with RMPs has proven difficult in DN configurations, where no device has thus far reported any hints of suppression. Modeling using the GPEC code finds a reduced high-field side response closer to DN shaping. The resulting synthetic diagnostic measurements are consistent with what has been observed in experiments on DIII-D, validating the plasma response model in this regime. While common metrics for suppression do not illustrate a clear distinction between single-null (SN) and DN cases, the pedestal top resonant field does show a ∼20% decrease at DN shaping in modeling. Field penetration is assessed using linear tearing theory, which demonstrates a lack of sufficient pedestal top resonant flux in the DN shape, requiring at least 1.5 × greater RMP coil currents than what was used in experiment. Analysis from drift kinetic simulations further indicate up to 2 × larger critical island widths are required at the pedestal top for tearing mode growth compared to SN cases. Effective island widths inferred from 3D ideal MHD are also analyzed, where maximum widths in lab coordinates indicate a threshold of ∼18–24× the ion gyroradius for sufficient profile flattening for ELM suppression. These results suggest that ELM suppression may be possible in DN with sufficiently large RMP coil amplitude. Future prospects of achieving RMP ELM suppression in the DN configuration may involve going to lower triangularity, which is also highlighted in this work.

RMP ELM suppression↗

DOE OSTI2026

Overview of the ASDEX Upgrade results

After a 26-month vent ASDEX Upgrade (AUG) went back in operation with a newly designed upper W-divertor suitable for alternative divertor configurations (featuring in-vessel coils and cryo-pump). Parameter scans and an extensive set of measurements were obtained and their interpretation is ongoing. Prompted by the ITER wall change, dedicated experiments on non-boronized plasma startup were contrasted to that employing asymmetric and more symmetric boronizations. The asymmetric boronization proved to be as beneficial as the more symmetric one, which is in contrast to previous model calculations assuming perfect sticking of boron (measurements suggest sticking ≈ 0.3). In the startup phase also the impurity influxes at the outboard limiters were investigated contrasting the unboronized case featuring cold edges (low- Z radiation) to the boronized case, in which the lifetime of the boron layers could be estimated. Pedestal stability investigations revealed that the quasi-continuous exhaust (QCE) regime is obtained when ballooning modes are active in the vicinity of the separatrix and the global peeling-ballooning stability is high enough. Thus, at high enough shaping and high gas flux both can be achieved and QCE is a consequence. The closely related enhanced D-alpha (EDA) mode is not clearly distinguishable from QCE, e.g. the quasi-coherent mode characteristic for EDA also shows up in QCE. In QCE the impurity transport is behaving benign as could be measured for Ne with a novel analysis method making use of a comprehensive set of CXRS measurements. For high radiative fractions the regime of the X-point radiator (XPR) is accessible at AUG and the understanding of its access conditions and behaviour is further developed. Due to the localized radiative cooling at the X-point the XPR can be well diagnosed and thus controlled. For negative triangularity shapes, further experiments at increased shaping resulted in strongly heated L-mode plasmas avoiding ELMs. Two integrated modelling approaches towards ITER suggested that core W-accumulation will be no issue for ITER and that the fusion yield in ITER may be Q = 12 (i.e. ITPA20-IL scaling is too pessimistic). Further, investigations of the ITER ramp-down in AUG provide insights into maintaining position control. Various aspects of shattered pellet injection were investigated in AUG and one of the results show that with increasing Ne fraction the radiation during the current quench increases and the current decay becomes faster.

ASDEX Upgrade↗

DOE OSTI2026

Overview of the MAST Upgrade physics programme: testing novel concepts at low aspect ratio to inform future devices

The research programme performed on the Mega Amp Spherical Tokamak (MAST) Upgrade experiment has made significant advances in developing the physics understanding of low aspect ratio tokamaks in support of the operation of ITER and design of fusion powerplants. High performance plasma scenarios have been developed to facilitate a broad programme of experiments, in which confinement is constrained by the presence of m/n = 2/1 modes that cause substantial losses of fast ions. The onset of these modes coincides with the q = 2 surface residing in a local minimum in the toroidal current density profile. The maximum electron temperature at the pedestal top, T e,ped is limited with gas fuelling to ∼350 eV to maintain regular ELMs; higher T e,ped results in a transition to a non-stationary ELM-free regime. The operational space of spherical tokamaks has been expanded into small and ELM-free regimes. Strong shaping of the last closed flux surface can induce a transition from large to small ELMs, and ELM suppression with resonant magnetic perturbations has been observed for the first time in a low aspect ratio tokamak. Negative triangularity shaping has induced a transition from ELMy H-mode to a high-performance L-mode regime for the first time in a low aspect ratio tokamak. In studies of fast ion confinement, losses of fast particles due to Global Alfvén Eigenmodes have been identified. Interactions between fast ions generated by off-axis neutral beam injection and thermal neutrals can result in significant losses of fast ions. Experiments with on- and off-axis neutral beam injection exhibit a flux pumping mechanism, where the central safety factor is held to ∼1 in the absence of sawteeth. In studies of pedestal physics, it has been found that elevated main chamber neutral pressures result in an increase in the electron density and reduction in the temperature at the pedestal top. Advances in understanding plasma exhaust include the integration of a high-performance plasma core with detached outer divertors in the X-point target configuration. A newly commissioned lower divertor cryopump reduces the lower divertor neutral pressure by up to 50%, with minimal effect on the main chamber or upper divertor. New measurements and SOLPS-ITER simulations emphasise the importance of plasma–neutral interactions on divertor detachment in the conditions accessible in experiments. Real-time control of the ionisation front location in both divertor chambers independently has been demonstrated in double null experiments, enabled by the tightly baffled divertor chambers.

MAST Upgrade↗

DOE OSTI2026

Progress and innovations in the TCV tokamak research programme

Research on the Tokamak à Configuration Variable addresses a wide range of key questions relevant to ITER and future fusion power plants. Over the past two years, highly productive experimental campaigns have led to major advances across several areas: the ITER baseline scenario; pedestal properties in low-collisionality, peeling-limited conditions; and the development of high- β N , non-inductive regimes. Alternative high-confinement scenarios have likewise received significant attention, with remarkable progress in quasi-continuous exhaust operation, X-point radiator plasmas, and negative triangularity configurations. Substantial achievements were also made in the mitigation or benign termination of runaway electron beams, in elucidating fast-ion loss mechanisms, and in improving exhaust behaviour in both conventional and alternative divertor geometries. These experimental results have been strongly supported by advances in modelling and their direct application to the experiment, ranging from gyrokinetic simulations of core and pedestal turbulence to fluid-based studies of scrape-off layer and divertor physics in diverse geometries. Plasma control has taken on an increasingly important role, with model-based and data-driven approaches now closely intertwined with physics studies. This article provides a overview of these recent activities, together with a brief outlook on forthcoming upgrades and next steps.

EPFL↗

DOE OSTI2022

Solving quasi-free and quadratic Lindblad master equations for open fermionic and bosonic systems

The dynamics of Markovian open quantum systems are described by Lindblad master equations. For fermionic and bosonic systems that are quasi-free, i.e. with Hamiltonians that are quadratic in the ladder operators and Lindblad operators that are linear in the ladder operators, we derive the equation of motion for the covariance matrix. This determines the evolution of Gaussian initial states and the steady states, which are also Gaussian. Using ladder super-operators (a.k.a. third quantization), here we show how the Liouvillian can be transformed to a many-body Jordan normal form which also reveals the full many-body spectrum. Extending previous work by Prosen and Seligman, we treat fermionic and bosonic systems on equal footing with Majorana operators, shorten and complete some derivations, also address the odd-parity sector for fermions, give a criterion for the existence of bosonic steady states, cover non-diagonalizable Liouvillians also for bosons, and include quadratic systems. In extension of the quasi-free open systems, quadratic open systems comprise additional Hermitian Lindblad operators that are quadratic in the ladder operators. While Gaussian states may then evolve into non-Gaussian states, the Liouvillian can still be transformed to a useful block-triangular form, and the equations of motion for k-point Green’s functions form a closed hierarchy. Based on this formalism, results on criticality and dissipative phase transitions in such models are discussed in a companion paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2019

Modeling the Kondo effect of a magnetic atom adsorbed on graphene

The Kondo effect due to a magnetic atom adsorbed on graphene is investigated using two theoretical approaches, including the non-equilibrium Green's function method within the slave-boson mean-field approximation and a newly developed tight-binding (TB) model with an effectively infinite Hubbard U term. Both methods reveal the presence of a Kondo peak in the local density of state (LDOS) of graphene near the Fermi level. A sharp Kondo peak is only observed in the odd-neighboring C sites of the C atom directly below the magnetic atom placed in top position. The peak intensity is found to decay quickly with respect to the distance between the C site and the magnetic atom. A theoretical model of scanning tunneling microscopy (STM) of graphene is presented as a means to identify the manifestation of the Kondo effect via direct topographic STM measurements. STM simulations show that in the proximity of the magnetic atom, only the C atoms of one sublattice are visible at low bias potential and thus a triangular lattice can be seen, in agreement with recent experiments. Further away from the magnetic atom, the Kondo effect dies down and eventually both sublattices are visible, leading to the recovery of the hexagonal lattice. Furthermore, this work not only provides a new TB method to study the Kondo effect in graphene, but also presents both scanning tunneling spectroscopy and microscopy fingerprints of the Kondo effect to facilitate its experimental verification.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Scanning tunneling microscopy study of epitaxial Fe 3 GeTe 2 monolayers on Bi 2 Te 3

Abstract Introducing magnetism to the surface state of topological insulators, such as Bi 2 Te 3 , can lead to a variety of interesting phenomena. We use scanning tunneling microscopy (STM) to study a single quintuple layer (QL) of the van der Waals magnet Fe 3 GeTe 2 (FGT) that is grown on Bi 2 Te 3 via molecular beam epitaxy. STM topographic images show that the FGT grows as free-standing islands on Bi 2 Te 3 and outwards from Bi 2 Te 3 steps. Atomic resolution imaging shows triangular lattices of 390 ± 10 pm for FGT and 430 ± 10 pm for Bi 2 Te 3 , consistent with the respective bulk crystals. A moiré pattern is observed on FGT regions with a periodicity of 4.3 ± 0.4 nm that can be attributed solely to this lattice mismatch and thus indicates zero rotational misalignment. While most of the surface is covered by a single QL of the FGT, there are small double QL regions, as well as regions with distinct chemical terminations due to an incomplete QL. The most common partial QL surface termination is the FeGe layer, in which the top two atomic layers are missing. This termination has a distinctive electronic structure and a 3 x 3 R 30 ∘ reconstruction overlaid on the moiré pattern in STM images. Magnetic circular dichroism measurements confirm these thin FGT films are ferromagnetic withT C ∼190 K.

Materials Science↗

DOE OSTI2020

Quantum electronics and optics at the interface of solid neon and superfluid helium

Here, we predict a new quantum electronic structure at the interface between two condensed phases of noble-gas elements: solid neon and superfluid helium. An excess electron injected onto this interface self-confines its wavefunction into a nanometric dome structure. Its size varies with pressure and optical transitions cover a broad mid-infrared spectrum. A collection of such electrons can form a classical Wigner crystal resembling a quantum-dot array on a triangular lattice. Under ultrafast laser illumination, this Wigner crystal can exhibit the quantum optical phenomenon of superradiance in the picosecond time scale. The ultralong spin-coherence time and micron-scale deterministic configurability allow the electrons in this system to serve as quantum information carriers. Their spin states can be controlled and readout by on-chip single-electron devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

An open-endcap blade trap for radial-2D ion crystals

We present the design and experimental demonstration of an open-endcap radio frequency trap to confine ion crystals in the radial-two dimensional (2D) structural phase. The central axis of the trap is kept free of obstructions to allow for site-resolved imaging of ions in the 2D crystal plane, and the confining potentials are provided by four segmented blade electrodes. We discuss the design challenges, fabrication techniques, and voltage requirements for implementing this open-endcap trap. Finally, we validate its operation by confining up to 29 ions in a 2D triangular lattice, oriented such that both in-plane principal axes of the 2D crystal lie in the radial direction.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2023

An finite element analysis surrogate model with boundary oriented graph embedding approach for rapid design

Abstract In this work, we present a boundary oriented graph embedding (BOGE) approach for the graph neural network to assist in rapid design and digital prototyping. The cantilever beam problem has been solved as an example to validate its potential of providing physical field results and optimized designs using only 10 ms. Providing shortcuts for both boundary elements and local neighbor elements, the BOGE approach can embed unstructured mesh elements into the graph and performs an efficient regression on large-scale triangular-mesh-based finite element analysis (FEA) results, which cannot be realized by other machine-learning-based surrogate methods. It has the potential to serve as a surrogate model for other boundary value problems. Focusing on the cantilever beam problem, the BOGE approach with 3-layer DeepGCN model achieves the regression with mean square error (MSE) of 0.011 706 (2.41% mean absolute percentage error) for stress field prediction and 0.002 735 MSE (with 1.58% elements having error larger than 0.01) for topological optimization. The overall concept of the BOGE approach paves the way for a general and efficient deep-learning-based FEA simulator that will benefit both industry and Computer Aided Design (CAD) design-related areas.

42 ENGINEERING↗

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)↗