Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “centrality measures”

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 757 records · Page 42

Crustal and uppermost mantle structure near the Gloria Fault, North Atlantic, from ocean bottom seismometer surface wave observations

SUMMARY In this work, we present both 1-D and 3-D shear wave velocity (Vs) models of the oceanic crust and uppermost mantle below the Deep OCean Test ARray area, located ∼ 70 km north of the central section of the Gloria Fault, in the eastern North-Atlantic Ocean. The velocity models are inferred from the dispersion of surface waves recorded on ocean bottom seismometers. Dispersion measurements are obtained from the analysis of ambient seismic noise at short periods (< 14 s) and teleseismic surface waves at long periods (> 14 s) using the two-station method. The 1-D Vs model is inferred from the joint inversion of Rayleigh wave phase velocities and Love wave group and phase velocities. The 3-D tomographic model is obtained by inversion of 2-D Love wave group velocity maps as a function of depth, further constrained by the average of Love wave phase velocities obtained from ambient noise (4–9 s) and the average Rayleigh and Love wave phase velocities calculated from teleseismic data (14–44 s). The 1-D Vs model shows a sediment layer with a low velocity of 1.05 km s−1, similar to previous studies in the region. Below the sediments, we find an oceanic crust with velocities ranging from 3.3 to 4.5 km s−1. The model reaches an unusually high velocity of 4.9 km s−1 in a 20 km thick layer at depths between 16 and 36 km. We interpret this fast velocity layer as indicative of the presence of harzburgite, a residue of enhanced melting that might have been formed by the proximity between the Mid-Atlantic Ridge and the Azores mantle plume. At greater depths the velocity decreases, forming a low-velocity zone that reaches a minimum at ∼ 70 km depth, which we interpret as the maximum depth for the lithosphere–asthenosphere boundary. The 3-D model shows a structure that is mostly horizontally layered, with Vs isocontours at 3.5–4.5 km s−1 highlighting oscillations of the crustal structure with wavelengths of ∼25–30 km. These oscillations may be due to changes in the rate of mantle upwelling and magma supply rate.

Pinzón, Juan I. (ORCID:0000000338078068)↗

3D structure of anisotropic flow in small collision systems at energies available at the BNL Relativistic Heavy Ion Collider

Here, we present (3 + 1)-dimensional [(3 + 1)D] dynamical simulations of asymmetric nuclear collisions at the BNL Relativistic Heavy Ion Collider (RHIC). Employing a dynamical initial state model coupled to (3 + 1)D viscous relativistic hydrodynamics, we explore the rapidity dependence of anisotropic flow in the RHIC small system scan at 200 GeV center-of-mass energy. We calibrate parameters to describe central 3 He + Au collisions and make extrapolations to d + Au and p + Au collisions. Our calculations demonstrate that approximately 50% of the v 3 (p T ) difference between the measurements by the STAR and PHENIX Collaborations can be explained by the use of reference flow vectors from different rapidity regions. This emphasizes the importance of longitudinal flow decorrelation for anisotropic flow measurements in asymmetric nuclear collisions, and the need for (3 + 1)D simulations. We also present results for the beam energy dependence of particle spectra and anisotropic flow in d + Au collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Decentralized Voltage Control of Large-Scale Distribution System with PVs Based on MADRL

This paper proposes a model-free decentralized control framework for the voltage regulation of large-scale distribution systems through the coordinated control of PV inverters. This is achieved by developing a novel interaction mechanism between the surrogate model and the centralized training and decentralized execution multiagent deep reinforcement learning framework. Specifically, the sparse Gaussian processes regression method is first utilized to develop the surrogate model of the original distribution system for reward calculation during the training stage, where each agent represents a sub-region in the centralized fashion for coordination strategy learning. After that, the learned control rules are used to inform controllers within each sub-region for real-time decisions with only local measurements. Comparative tests among various methods on the EPRI Ckt5 test system demonstrate the effectiveness of the proposed method.

distribution system↗

Leaf turgor loss point shapes local and regional distributions of evergreen but not deciduous tropical trees

Summary The effects of climate change on tropical forests will depend on how diverse tropical tree species respond to drought. Current distributions of evergreen and deciduous tree species across local and regional moisture gradients reflect their ability to tolerate drought stress, and might be explained by functional traits. We measured leaf water potential at turgor loss (i.e. ‘wilting point’; π tlp ), wood density (WD) and leaf mass per area (LMA) on 50 of the most abundant tree species in central Panama. We then tested their ability to explain distributions of evergreen and deciduous species within a 50 ha plot on Barro Colorado Island and across a 70 km rainfall gradient spanning the Isthmus of Panama. Among evergreen trees, species with lower π tlp were associated with drier habitats, with π tlp explaining 28% and 32% of habitat association on local and regional scales, respectively, greatly exceeding the predictive power of WD and LMA. In contrast, π tlp did not predict habitat associations among deciduous species. Across spatial scales, π tlp is a useful indicator of habitat preference for tropical tree species that retain their leaves during periods of water stress, and holds the potential to predict vegetation responses to climate change.

59 BASIC BIOLOGICAL SCIENCES↗

pnnl/grid_prediction

Two datadriven predictive approaches, namely, {\em Koopman Operator Theoretic (KOT)-based} model and {\em Graph Neural Network (GNN)-based model}, to enable effective power system state predictions. The KOT-based approaches (Robust DMD, deepDMD) capture the power system evolution as a linear dynamical system on an abstract space. The GNNs model the spatio-temporal correlations using graph convolutional network and are called Spatio-Temporal Graph Convolutional Network (STGCN). These predictive models are trained, tested and compared rigorously based on their predictions of frequencies in the IEEE 68 bus system when subjected to a disturbance. GridSTAGE framework developed at Pacific Northwest National Laboratory is leveraged to generate multiple datasets (in the form of PMU measurements) for training and testing by strategically creating load changes across the spatial locations of the network

Nandanoor, Sai Pushpak↗

Utilizing a Storm-Generating Hotspot to Study Convective Cloud Transitions: The CACTI Experiment

The Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign was designed to improve understanding of orographic cloud life cycles in relation to surrounding atmospheric thermodynamic, flow, and aerosol conditions. The deployment to the Sierras de Córdoba range in north-central Argentina was chosen because of very frequent cumulus congestus, deep convection initiation, and mesoscale convective organization uniquely observable from a fixed site. The C-band Scanning Atmospheric Radiation Measurement (ARM) Precipitation Radar was deployed for the first time with over 50 ARM Mobile Facility atmospheric state, surface, aerosol, radiation, cloud, and precipitation instruments between October 2018 and April 2019. An intensive observing period (IOP) coincident with the RELAMPAGO field campaign was held between 1 November and 15 December during which 22 flights were performed by the ARM Gulfstream-1 aircraft. A multitude of atmospheric processes and cloud conditions were observed over the 7-month campaign, including numerous orographic cumulus and stratocumulus events; new particle formation and growth producing high aerosol concentrations; drizzle formation in fog and shallow liquid clouds; very low aerosol conditions following wet deposition in heavy rainfall; initiation of ice in congestus clouds across a range of temperatures; extreme deep convection reaching 21-km altitudes; and organization of intense, hail-containing supercells and mesoscale convective systems. These comprehensive datasets include many of the first ever collected in this region and provide new opportunities to study orographic cloud evolution and interactions with meteorological conditions, aerosols, surface conditions, and radiation in mountainous terrain.

54 ENVIRONMENTAL SCIENCES↗

Improving the Performance of NEML2 with Modern Graph Compilation Backends

NEML2 vectorizes constitutive-model evaluation for large-scale multiphysics simulation, using PyTorch as its tensor backend so that a batch of material-point updates runs on CPU or GPU through a single implementation. In the two prior reports in this series it was a C++-native library, deployed through TorchScript tracing and just-in-time (JIT) compilation; it has since been rewritten from the ground up into a Python-native library deployed through Ahead-of-Time Inductor (AOTInductor), a modern PyTorch graph-compilation backend. The rewrite is driven by a persistent tension, not a language preference: NEML2 composes constitutive models at runtime from a registry of small, independently-authored pieces, and that flexibility is difficult to reconcile with the compile-time knowledge an efficient GPU kernel needs. This report documents the rewrite and the investment that accompanied it: the AOTInductor export pipeline that turns a Python-authored model into a portable, Python-free compiled artifact loadable from pure C++; the eager and compiled runtimes and the new implicit solver layer built on them; a head-to-head benchmark of legacy JIT against AOTInductor; the physics-model catalog and its worked examples; the developer tooling; and the corresponding overhaul of MOOSE’s NEML2 integration that lets MOOSE consume it. A central objective is to examine whether modern PyTorch graph-compilation backends are effective for MOOSE GPU integration. The benchmark answers directly: AOTInductor outperforms legacy JIT on every GPU scenario measured, by 1.0–4.5×. Modern graph-compilation backends are effective for MOOSE GPU integration, and AOTInductor specifically – not compilation in the abstract – is why.

Hu, Gary (Tianchen) [Argonne National Laboratory (↗

Wind and Temperature Consensus at Horn Point, HU-Beltsville, Piney Run (Maryland) in support of CoURAGE

The Maryland Department of the Environment (MDE) operates a ground-based atmospheric profiling network consisting of collocated radar wind profilers (RWP) and radio acoustic sounding systems (RASS) as part of its Ambient Air Monitoring Program. This network provides continuous observations of wind and temperature structure in the lower troposphere to support air quality forecasting, regulatory analysis, and atmospheric research. The network currently includes three fixed sites across Maryland: Horn Point (HP, lower eastern shore) [38.587525°,-76.141006°], Howard University-Beltsville (HUB, central Maryland) [39.055277°, -76.878632°], and Piney Run (PR, western Maryland) [39.705950°, -79.012000°] The network is designed to capture regional variability in atmospheric transport and boundary-layer processes. These systems measure vertical profiles of horizontal wind speed and direction using Doppler radar techniques, with observations typically spanning from ~100 m above ground level up to approximately 2.5–4 km. Measurements are derived from the Doppler shift of backscattered electromagnetic signals, enabling retrieval of wind vectors at multiple altitudes with high temporal resolution (e.g., 30-minute averages reported every 6 minutes). Each radar wind profiler is paired with a Radio Acoustic Sounding System (RASS) to provide profiles of virtual temperature in the lower atmosphere (~100–200 m AGL) by measuring the propagation speed of acoustic waves. Together, the RWP/RASS system yields a coupled data set of thermodynamic and kinematic atmospheric structure, including additional parameters such as vertical velocity, radial velocity, signal-to-noise ratio, and spectral width for advanced analysis. There are two types of files for each station: wind data (files with a "w" prefix) and virtual temperature RASS data (files with a "t" prefix). The wind data files are in the format wYYDDD.cns, where YY is the 2-digit year and DDD is the day of the year. The RASS virtual temperature data files are in the format tYYDDD.cns. Each record has the following header structure: Line 1 : Station Name RASS files Line 2 : RASS rev DeTect_2.0, WINDS files Line 2 : WINDS rev ATI 5.1 Line 3 : N latitude, W longitude, and site elevation (m) Line 4 : Date and begin time of consensus: yy mm dd hh mn ss plus # minutes to add to get UTC Line 5 : Consensus averaging time (minutes); number of beams; number of range gates Line 6 : Number of records required to make consensus (num) total number of records (tot) and the consensus window size (m/s) in the format: num:tot (window) RASS files Line 7 : no. of coded cells, no. of spec, pulse width (ns), and inter-pulse period (µs), WINDS files Line 7 : No. of coded cells, no. of spectra, pulse width (ns), and inter-pulse period (µs), each with a pair of values: first value is for oblique beams, second for vertical RASS files Line 8 : Full scale Doppler value (m/s) Delay to first gate (ns) Number of gates Spacing of gates (ns), WINDS files Line 8 : Full scale Doppler velocity (m/s), oblique and vertical Vertical correction applied to oblique beams? (0 = no, 1 = yes) Delay to first gate (ns), oblique and vertical Number of gates, oblique and vertical Spacing of gates (ns), oblique and vertical Line 9 : Azimuth and elevation (9s indicate vertical beam not used) RASS files Line 10, values : HT = Height above ground (km), T = Uncorrected virtual temperature consensus (deg C), Tc = Corrected virtual temperature consensus (deg C), W = Vertical wind consensus (9s indicate vertical beam not used, w-component, positive upward, m/s), CNT = Number of records that made consensus (for the 3 values in same order), SNR = Average signal to noise ratio (dB) of records in consensus (same order) WINDS files Line 10, values : HT = Height above ground (km), SPD = Wind speed (m/s), DIR = Wind direction (deg E of N from N), RAD = Radial velocities for each beam (m/s) in order given in azimuth and elevation line (positive toward radar; 9s indicate vertical beam not used, CNT = Number of records that made consensus, SNR = Average signal to noise ratio (dB) of records in consensus

{"wind speed and direction",temperature}↗

Outflows in star-forming galaxies: Stacking analyses of resolved winds and the relation to their hosts’ properties

ABSTRACT Outflows form an integral component in regulating the gas cycling in and out of galaxies, although their impact on the galaxy hosts is still poorly understood. Here we present an analysis of 405 high mass (log M*/M⊙ ≥ 10), star-forming galaxies (excluding AGN) with low inclinations at z ∼ 0, using stacking techniques of the Na D λλ5889, 5895 Å neutral gas tracer in IFU observations from the MaNGA DR15 survey. We detect outflows in the central regions of 78/405 galaxies and determine their extent and power through the construction of stacked annuli. We find outflows are most powerful in central regions and extend out to ∼1Re, with declining mass outflow rates and loading factors as a function of radius. The stacking of spaxels over key galaxy quantities reveals outflow detections in regions of high ΣSFR (≳0.01 M⊙ yr−1 kpc−2) and $\Sigma _{M_{*}}$ (≳107 M⊙ kpc−2) along the resolved main sequence. Clear correlations with ΣSFR suggest it is the main regulator of outflows, with a critical threshold of ∼0.01 M⊙ yr−1 kpc−2 needed to escape the weight of the disc and launch them. Furthermore, measurements of the Hδ and Dn4000 indices reveal virtually identical star formation histories between galaxies with outflows and those without. Finally, through stacking of H i 21 cm observations for a subset of our sample, we find outflow galaxies show reduced H i gas fractions at central velocities compared to their non-detection control counterparts, suggestive of some removal of H i gas, likely in the central regions of the galaxies, but not enough to completely quench the host.

Roberts-Borsani, G. W.↗

Channelrhodopsin C1C2: Photocycle kinetics and interactions near the central gate

Channelrhodopsins (ChR) are light-sensitive cation channels used in optogenetics, a technique that applies light to control cells (e.g., neurons) that have been modified genetically to express those channels. Although mutations are known to affect pore kinetics, little is known about how mutations induce changes at the molecular scale. To address this issue, we first measured channel opening and closing rates of a ChR chimera (C1C2) and selected variants (N297D, N297V, and V125L). Then, we used atomistic simulations to correlate those rates with changes in pore structure, hydration, and chemical interactions among key gating residues of C1C2 in both closed and open states. Overall, the experimental results show that C1C2 and its mutants do not behave like ChR2 or its analogous variants, except V125L, making C1C2 a unique channel. Furthermore, our atomistic simulations confirmed that opening of the channel and initial hydration of the gating regions between helices I, II, III, and VII of the channel occurs with 1) the presence of 13-cis retinal; 2) deprotonation of a glutamic acid gating residue, E129; and 3) subsequent weakening of the central gate hydrogen bond between the same glutamic acid E129 and asparagine N297 in the central region of the pore. Also, an aspartate (D292) is the unambiguous primary proton acceptor for the retinal Schiff base in the hydrated channel.

59 BASIC BIOLOGICAL SCIENCES↗

A precise measurement of the jet energy scale derived from single-particle measurements and in situ techniques in proton–proton collisions at $\sqrt{s}=$ 13 TeV with the ATLAS detector

The jet energy calibration and its uncertainties are derived from measurements of the calorimeter response to single particles in both data and Monte Carlo simulation using proton–proton collisions at $\sqrt{s} = 13$ TeV collected with the ATLAS detector during Run 2 at the Large Hadron Collider. The jet calibration uncertainty for anti-$k_T$ jets with a jet radius parameter of R$_\textrm{jet} = 0.4$ and in the central jet rapidity region is about 2.5% for transverse momenta ($p_{\text {T}}$) of 20 $\text {GeV}$ , about 0.5% for $p_{\text {T}} = 300$ GeV and 0.7% for $p_{\text {T}} = 4$ TeV . Excellent agreement is found with earlier determinations obtained from -balance based in situ methods ($Z/\gamma$ +jets). The combination of these two independent methods results in the most precise jet energy measurement achieved so far with the ATLAS detector with a relative uncertainty of 0.3% at $p_\textrm{T} = 300$ GeV and 0.6% at 4 TeV. The jet energy calibration is also derived with the single-particle calorimeter response measurements separately for quark- and gluon-induced jets and furthermore for jets with R jet varying from 0.2 to 1.0 retaining the correlations between these measurements. Differences between inclusive jets and jets from boosted top-quark decays, with and without grooming the soft jet constituents, are also studied.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Combining Eddy Covariance Towers, Field Measurements, and the MEMS 2 Ecosystem Model Improves Confidence in the Climate Impacts of Bioenergy With Carbon Capture and Storage

ABSTRACT Carbon dioxide removal technologies such as bioenergy with carbon capture and storage (BECCS) are required if the effects of climate change are to be reversed over the next century. However, BECCS demands extensive land use change that may create positive or negative radiative forcing impacts upstream of the BECCS facility through changes to in situ greenhouse gas fluxes and land surface albedo. When quantifying these upstream climate impacts, even at a single site, different methods can give different estimates. Here we show how three common methods for estimating the net ecosystem carbon balance of bioenergy crops established on former grassland or former cropland can differ in their central estimates and uncertainty. We place these net ecosystem carbon balance forcings in the context of associated radiative forcings from changes to soil N 2 O and CH 4 fluxes, land surface albedo, embedded fossil fuel use, and geologically stored carbon. Results from long term eddy covariance measurements, a soil and plant carbon inventory, and the MEMS 2 process‐based ecosystem model all agree that establishing perennials such as switchgrass or mixed prairie on former cropland resulted in net negative radiative forcing (i.e., global cooling) of −26.5 to −39.6 fW m −2 over 100 years. Establishing these perennials on former grassland sites had similar climate mitigation impacts of −19.3 to −42.5 fW m −2 . However, the largest climate mitigation came from establishing corn for BECCS on former cropland or grassland, with radiative forcings from −38.4 to −50.5 fW m −2 , due to its higher plant productivity and therefore more geologically stored carbon. Our results highlight the strengths and limitations of each method for quantifying the field scale climate impacts of BECCS and show that utilizing multiple methods can increase confidence in the final radiative forcing estimates.

Falvo, Grant [Department of Plant, Soil and Microb↗

Measurements of the Hubble Constant with a Two-rung Distance Ladder: Two Out of Three Ain’t Bad

Abstract The three-rung distance ladder, which calibrates Type Ia supernovae (SNe Ia) through stellar distances linked to geometric measurements, provides the highest precision direct measurement of the Hubble constant. In light of the Hubble tension, it is important to test the individual components of the distance ladder. For this purpose, we report a measurement of the Hubble constant from 35 extragalactic Cepheid hosts measured by the SH0ES team, using their distances and redshifts at cz ≤ 3300 km s −1 , instead of any more distant SNe Ia, to measure the Hubble flow. The Cepheid distances are calibrated geometrically in the Milky Way, NGC 4258, and the Large Magellanic Cloud. Peculiar velocities are a significant source of systematic uncertainty at z ∼ 0.01, and we present a formalism for both mitigating and quantifying their effects, making use of external reconstructions of the density and velocity fields in the nearby universe. We identify a significant source of uncertainty originating from different assumptions about the selection criteria of this sample, whether distance or redshift limited, as it was assembled over three decades. Modeling these assumptions yields central values ranging from H 0 = 71.7 to 76.4 km s −1 Mpc −1 . Combining the four best-fitting selection models yields H 0 = 72.9 − 2.2 + 2.4 as a fiducial result, at 2.4 σ tension with Planck. While SNe Ia are essential for a precise measurement of H 0 , unknown systematics in these supernovae are unlikely to be the source of the Hubble tension.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurement of coherent exclusive J/ψ → μ+μ− production in ultraperipheral Pb+Pb collisions at sNN=5.36 TeV with the ATLAS detector

The ATLAS experiment has performed a measurement of coherent exclusive J/ψ → μ+μ− production in ultraperipheral Pb+Pb collisions at sNN=5.36$$ \sqrt{s_{\textrm{NN}}}=5.36 $$ TeV. The data was recorded at the Large Hadron Collider (LHC) during 2023, and corresponds to an integrated luminosity of 79 μb−1. Exclusive J/ψ candidates were selected with a dedicated track-sensitive trigger based on the ATLAS transition radiation tracker. The analysis involves reconstruction of the dimuon invariant mass based on muon tracks from the inner detector, as the muon transverse momentum range of interest precludes the use of the standard muon reconstruction and identification algorithms. Differential cross sections are measured as a function of J/ψ rapidity and are compared with theoretical predictions. After extrapolation to sNN=5.02$$ \sqrt{s_{\textrm{NN}}}=5.02 $$ TeV, they are also compared with previous measurements performed by other experiments using data from LHC Run 2. While the results agree reasonably well with theoretical predictions, they are in tension with previous Run-2 results for the central rapidity region.

Aad, G↗

Generalized quantum master equations can improve the accuracy of semiclassical predictions of multitime correlation functions

Multitime quantum correlation functions are central objects in physical science, offering a direct link between the experimental observables and the dynamics of an underlying model. While experiments such as 2D spectroscopy and quantum control can now measure such quantities, the accurate simulation of such responses remains computationally expensive and sometimes impossible, depending on the system’s complexity. A natural tool to employ is the generalized quantum master equation (GQME), which can offer computational savings by extending reference dynamics at a comparatively trivial cost. However, dynamical methods that can tackle chemical systems with atomistic resolution, such as those in the semiclassical hierarchy, often suffer from poor accuracy, limiting the credence one might lend to their results. By combining work on the accuracy-boosting formulation of semiclassical memory kernels with recent work on the multitime GQME, here we show for the first time that one can exploit a multitime semiclassical GQME to dramatically improve both the accuracy of coarse mean-field Ehrenfest dynamics and obtain orders of magnitude efficiency gains.

Chemistry↗

Classical 1/3 scaling of convection holds up to Ra = 10 15

Significance The heat transport law in turbulent convection remains central to current research in the field. Our present knowledge of the heat transport law for R a > 1 0 12 is inconclusive, where the Rayleigh number R a is a measure of the strength of convection. Massively parallel simulations of the three-dimensional convection have progressed to R a = 1 0 15 in slender cells. We resolve velocity gradients inside thin boundary layers and show that the turbulent heat transport continues to follow the classical 1 / 3 scaling law with no transition to the so-called “ultimate” state that is variously argued to have 1 / 2 scaling. Our work suggests that the boundary layers remain marginally stable and continue to act as the bottleneck for global heat transport.

Science & Technology - Other Topics↗

An efficient numerical model for predicting residual stress and strain in parts manufactured by laser powder bed fusion

Abstract Computational modeling of additively manufactured structures plays an increasingly important role in product design and optimization. For laser powder bed fusion processes, the accurate modeling of stress and distortion requires large amount of computational cost due to very localized heat input and evolving complex geometries. The current study takes advantage of a graphics processing unit accelerated explicit finite element analysis code and approximated heat conduction analysis to predict the macroscopic thermo-mechanical behavior in laser selective melting. Adjacent layers and tracks were lumped to reduce the number of time steps and elements in the finite element model. The effects of track and layer grouping on prediction accuracy and solution efficiency are investigated to provide a guidance for a cost-effective simulation. Thin-wall builds from Inconel alloy 625 (IN625) powders were simulated by applying the developed modeling approach to get the detailed residual stress and distortion at a computational speed 50 times higher than conventional approach. Under repeated heating and cooling cycles, a high tensile stress was produced near surfaces of a build due to a larger shrinkage on surface than that in central area. It is also shown that horizontal stresses concentrate near the root and top layers of the IN625 build. The predicted residual elastic strain distribution was validated by the experimental measurement using x-ray synchrotron diffraction.

36 MATERIALS SCIENCE↗

Galaxy properties as revealed by MaNGA – III. Kinematic profiles and stellar population gradients in S0s

ABSTRACT This is the third paper of a series where we study the stellar population gradients (SP; ages, metallicities, α-element abundance ratios, and stellar initial mass functions) of early-type galaxies (ETGs) at $z$ ≤ 0.08 from the Mapping Nearby Galaxies at APO Data Release 15 (MaNGA-DR15) survey. In this work, we focus on the S0 population and quantify how the SP varies across the population as well as with galactocentric distance. We do this by measuring Lick indices and comparing them to SP synthesis models. This requires spectra with high signal-to-noise ratio which we achieve by stacking in bins of luminosity (Lr) and central velocity dispersion (σ0). We find that: (1) there is a bimodality in the S0 population: S0s more massive than $3\times 10^{10}\, \mathrm{M}_\odot$ show stronger velocity dispersion and age gradients (age and σr decrease outwards) but little or no metallicity gradient, while the less massive ones present relatively flat age and velocity dispersion profiles, but a significant metallicity gradient (i.e. [M/H] decreases outwards). Above $2\times 10^{11}\, \mathrm{M}_\odot$, the number of S0s drops sharply. These two mass scales are also where global scaling relations of ETGs change slope. (2) S0s have steeper velocity dispersion profiles than fast-rotating elliptical galaxies (E-FRs) of the same luminosity and velocity dispersion. The kinematic profiles and SP gradients of E-FRs are both more similar to those of slow-rotating ellipticals (E-SRs) than to S0s, suggesting that E-FRs are not simply S0s viewed face-on. (3) At fixed σ0, more luminous S0s and E-FRs are younger, more metal rich and less α-enhanced. Evidently for these galaxies, the usual statement that ‘massive galaxies are older’ is not true if σ0 is held fixed.

Domínguez Sánchez, H.↗