Search NASA⌕ Search

SEARCH · Search NASA

Results for “strong scaling”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Preliminary Techno-Economic Assessment of Gas Switching Reforming (GSR) of Natural Gas for Pure Hydrogen Production and Power Generation with Integrated CO2 Capture

The increasing demand for hydrogen and the CO2 intensity of natural gas (NG) reforming motivate the development of low-carbon-emission hydrogen production technologies. Gas Switching Reforming (GSR) is an advanced auto-thermal reforming technology that produces hydrogen or syngas from natural gas. It integrates inherent CO2 capture by utilizing a specialized oxygen carrier in a single fluidized bed reactor, eliminating the need for complex, energy-intensive post-combustion separation. GSR technology builds upon Chemical Looping Reforming (CLR), an experimentally proven technology with strong potential for scaling up. In this study, select oxygen carriers (OC) (NiO/Al2O3, Fe2O3-CeO2/Al2O3, and magnetite) were tested in methane steam reforming in a fixed bed reactor to determine their relative reactivities under relevant conditions for GSR (800 °C, 7 bar total pressure). Process models were then developed to perform techno-economic analysis (TEA) of GSR for hydrogen production (GSR-H₂) and a combined cycle (GSR-CC) in which high-purity H₂ is fired in a gas turbine to produce electricity. Operating at 10 bar and 1100 °C and with the additional recovery steps implemented increased H₂ production by ~30% and improved efficiency relative to prior studies. For GSR-H₂, the levelized cost of hydrogen (LCOH) is 1.61–1.64 $/kg-H₂, competitive with a reference SMR case, though operating and maintenance costs are higher due to increased electricity demand. GSR-CC has a significantly higher levelized cost of electricity (LCOE) than its reference NGCC (natural gas combined cycle) plant, suggesting it is less competitive; however, increasing production scale could make it more attractive. Life-cycle results for GSR-H2 indicate NG consumption drives ~75% of total global warming impacts (~2.3 kg CO2 eq/kg H2). An environmental, health, and safety screening suggests iron-based carriers are comparatively safer, whereas NiO may pose greater risks. Overall, GSR-H2 is a scalable, competitive option for hydrogen production using nickel and non-nickel OC.

03 NATURAL GAS↗
DOE OSTI2025

Total Recycling of Copper Cable Scrap and Production of Carbon Using Fast Microwave Technology

The recycling of cable scrap, particularly from discarded electrical wiring, is gaining significant attention due to the rising demand for copper and the need for sustainable management of electronic waste. Traditionally, mechanical and thermal processings have been used to recover copper and plastic from cables. However, these approaches are often energy-intensive, time-consuming, and costly in terms of equipment and labor. In this study, we present a simple and effective method for recovering materials from cable scrap using a domestic microwave oven. Cable pieces (2–2.5 cm long) were exposed to 700 W of microwave irradiation under rotation for 30 s, enabling the rapid and efficient separation of high-quality copper metal from the core wire, and activated carbon from the carbonized plastic sheath. Microwaves facilitate this process through Ohmic heating, which induces electrical resistance in the metal, generating heat that mechanically loosens the metal and carbonized plastic components. The process demonstrates high efficiency, achieving an 80% reduction in energy consumption compared to conventional processings. This fast and energy-efficient method shows strong potential for scaling up to industrial recycling, offering a cost-effective and environmentally friendly way to recover high-quality materials for further use or repurposing.

Bourlinos, Athanasios B. (ORCID:0000000256165993)↗
DOE OSTI2026

Idealized simulations of wind farm interactions with intermittent turbulence in stable boundary layer conditions

Stable atmospheric boundary layer conditions typically correspond to weak turbulence levels, but intermittent periods of elevated turbulence can occur during otherwise quiescent conditions. The interaction between intermittent turbulence and wind turbines is not well understood because of sparse observations, as well as the difficulty in realistically resolving small-scale turbulence during strongly stable conditions with numerical simulations. In this study, an explicit filtering and reconstruction approach for large-eddy simulation (LES) is used to simulate weakly and strongly stable conditions, with surface cooling rates of −0.2 and −2.0 K h −1 , respectively. This approach can sustain resolved background turbulence at relatively coarse grid spacing and stronger stratification compared to conventional closures, permitting more realistic intermittent stable boundary layer (SBL) turbulence. The idealized LES capability of the Weather Research and Forecasting model is employed with turbine rotors parameterized using generalized actuator disks to examine (1) how the presence of turbine wakes affects SBL evolution and (2) the effect of intermittent turbulence on power production and wake recovery. Wakes increase mixing and deepen the SBL, with a stronger effect under strongly stable conditions, primarily because the SBL is shallower and closer to the top of the wind turbine rotor layer. Intermittent turbulence does not have a significant impact on mean power generation and wake recovery because the relevant intermittent turbulent structures in this study only affect the bottom half of the rotor disk. Power production is, however, more variable during periods of elevated turbulence, demonstrating the impact of SBL intermittency. This study uses an idealized configuration, focusing on LES model performance and physical understanding, with the goal of informing future simulations of the conditions observed during the American Wake Experiment.

Energy - Wind↗
DOE OSTI2020

Hydraulic fracturing experiments at 1500 m depth in a deep mine: Highlights from the kISMET project

In support of the U.S. DOE SubTER Crosscut initiative, we established a field test facility in a deep mine and designed and carried out in situ hydraulic fracturing experiments relevant to enhanced geothermal systems (EGS) in crystalline rock to characterize the stress field, understand the effects of rock fabric on fracturing, and gain experience in monitoring using geophysical methods. The project also included pre- and post-fracturing simulation and analysis, and laboratory measurements and experiments. The kISMET (permeability (k) and Induced Seismicity Management for Energy Technologies) site was established in the West Access Drift of the Sanford Underground Research Facility (SURF) 4757 ft (1450 m) below ground (on the 4850 ft level (4850L)) in phyllite of the Precambrian Poorman Formation. We drilled and continuously cored five near-vertical boreholes in a line on 3 m (10 ft) spacing, deviating the two outermost boreholes slightly to create a five-spot pattern around the test borehole centered in the test volume 40 m below the drift invert (floor) at a total depth of ~1490 m (4890 ft). Laboratory measurements of core from the center test borehole showed P-wave velocity heterogeneity along each core indicating strong, fine-scale (~1 cm or smaller) changes in the mechanical properties of the rock. Field measurements of the stress field by hydraulic fracturing showed that the minimum horizontal stress at the kISMET site averages 21.7 MPa (3146 psi) trending approximately N-S (356 degrees azimuth) and plunging slightly NNW at 12°. The vertical and horizontal maximum stresses are similar in magnitude at 42-44 MPa (6090-6380 psi) for the depths of testing, which averaged approximately 1530 m (5030 ft). Hydraulic fractures were remarkably uniform suggesting core-scale and larger rock fabric did not play a role in controlling fracture orientation. Analytical solutions suggest that the fracture radius of the large fracture (stimulation test) was more than 6 m (20 ft), depending on the unknown amount of leak-off.

Oldenburg, C↗
DOE OSTI2024

A Scalable Real-Time Data Assimilation Framework for Predicting Turbulent Atmosphere Dynamics

AI-based foundation models like FourCastNet, GraphCast are revolutionizing weather and climate predictions but are not yet ready for operational use. Their limitation lies in the absence of a data assimilation system to incorporate real-time Earth system observations, crucial for accurately forecasting events like tropical cyclones. To overcome these obstacles, we introduce a generic real-time data assimilation framework and demonstrate its end-to-end performance on the Frontier supercomputer. This framework comprises two primary modules: an ensemble score filter (EnSF), which significantly outperforms the state-of-the-art data assimilation method, and a vision transformer-based surrogate capable of real-time adaptation through the integration of observational data. We demonstrate both the strong and weak scaling of our framework up to 1024 GPUs on the Exascale supercomputer, Frontier. Our results not only illustrate the framework's exceptional scalability on high-performance computing systems, but also demonstrate the importance of supercomputers in real-time data assimilation for weather and climate predictions.

Lu, Dan↗
DOE OSTI2021

Relativistically transparent magnetic filaments: scaling laws, initial results and prospects for strong-field QED studies

Abstract Relativistic transparency enables volumetric laser interaction with overdense plasmas and direct laser acceleration of electrons to relativistic velocities. The dense electron current generates a magnetic filament with field strength of the order of the laser amplitude (>10 5 T). The magnetic filament traps the electrons radially, enabling efficient acceleration and conversion of laser energy into MeV photons by electron oscillations in the filament. The use of microstructured targets stabilizes the hosing instabilities associated with relativistically transparent interactions, resulting in robust and repeatable production of this phenomenon. Analytical scaling laws are derived to describe the radiated photon spectrum and energy from the magnetic filament phenomenon in terms of the laser intensity, focal radius, pulse duration, and the plasma density. These scaling laws are compared to 3D particle-in-cell (PIC) simulations, demonstrating agreement over two regimes of focal radius. Preliminary experiments to study this phenomenon at moderate intensity ( a 0 ∼ 30) were performed on the Texas Petawatt Laser. Experimental signatures of the magnetic filament phenomenon are observed in the electron and photon spectra recorded in a subset of these experiments that is consistent with the experimental design, analytical scaling and 3D PIC simulations. Implications for future experimental campaigns are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Polymer Macrocycles: A novel topology to control dynamics of rubbery materials (Final Report)

One of the current challenges in polymer physics and engineering is to understand the behavior of circular macromolecules. Though such molecules have been made by dilute solution ring closure methods in the past, such procedures have been limited to relatively small, unentangled rings as well as to small sample quantities. Initial studies were also hampered by issues surrounding purity of the rings, though this problem has been somewhat overcome by the use of liquid chromatography at critical conditions (LCCC) methods that separate rings and linear chains. Furthermore, in spite of the continued progress a full understanding of ring dynamics remains elusive. The present work has addressed multiple aspects of the behavior of ring molecules by using a novel synthetic route: reversible radical recombination redox polymerization (R3P) to produce poly(3,6-dioxa-1,8-octanedithiol) (PolyDODT) molecules. The synthesis conditions were chosen to produce linear (LDODT), cyclic (RDODT) and linear-cyclic mixtures (LRDODT) in order to investigate dilute solution and rheological responses of the materials. Importantly, the R3P synthesis makes larger quantities and higher molecular weights than do the dilute solution ring closure methods and also has the potential to be scaled to industrial quantities, thus becoming commercializable. We emphasize two important outcomes of the work. The first is that the rheological investigation shows that very large molecular weight the PolyDODT rings permitted the investigation of ring dynamics from glassy towards terminal flow for molecular sizes corresponding to approximately 300 entanglements of the linear counterpart. At the same time, these very large sized rings could be diluted in a solvent in order to reduce the molecular entanglement density to less than unity. The results demonstrate that ring viscosity follows a Rouse-like linear dependence on molecular weight or entanglement number to approximately 15 entanglements. This range of "Rouse-like" response is greater than what was achieved in all prior works where synthesis conditions have limited the molecular sizes achievable. After the entanglement number $Z_w$ increases beyond approximately 15, we find that the viscosity-molecular weight scaling undergoes a rapid change from the Rouse-like $η\sim{Z}^1_w$ to a very strong power-law scaling that is approximately $η\sim{Z}^{5.8}_w$. This aspect of the work strongly suggests that prior investigations on highly pure rings have invariably been performed on unentangled systems simply because ring molecules have very high entanglement thresholds. The results are further supported by plateau modulus measurements that show that in the "Rouse-like" viscosity regime, for the RDODT samples there is a corresponding lack of a rubbery entanglement plateau and once entanglement coupling is seen, the rubbery plateau also appears in the dynamic modulus data. Furthermore, upon appropriate scaling, we show that the unentangled polyDODT rings (i.e., the diluted high molecular weight rings) show virtually identical dynamics (upon appropriate scaling to account for the chain dilutions) to the dynamics of similarly unentangled LCCC fractionated polystyrene rings.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Search for the critical point of strongly-interacting matter in $^{40}$Ar + $^{45}$Sc collisions at 150A Ge V /c using scaled factorial moments of protons

The critical point of dense, strongly interacting matter is searched for at the CERN SPS in $^{40}$Ar + $^{45}$Sc collisions at 150A Ge V /c. The dependence of second-order scaled factorial moments of proton multiplicity distribution on the number of subdivisions of transverse momentum space is measured. The intermittency analysis is performed using both transverse momentum and cumulative transverse momentum. For the first time, statistically independent data sets are used for each subdivision number. The obtained results do not indicate any statistically significant intermittency pattern. An upper limit on the fraction of correlated proton pairs and the power of the correlation function is obtained based on a comparison with the Power-law Model developed for this purpose.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Local Spin Density Approximation Strongly Improved by a Better-Informed Local Scaling of Its Self-Interaction Correction

The Perdew−Zunger self-interaction correction (PZSIC) makes density functional approximations (DFAs) exact for all one-electron densities. However, it overcorrects in manyelectron regions, introducing errors for the uniform-density limit, where uncorrected DFAs are exact. The locally scaled PZSIC (LSIC), based on the iso-orbital indicator zσ [which distinguishes single-orbital and slowly varying density regions and is used with the local spin density approximation (LSDA)], restores the uniform-density limit and significantly improves results for many properties, including chemical reaction barrier heights, atomization energies, and ionization potentials. Yet, LSIC performs poorly for weakly bonded systems, leaving many unbound, due to limitations of its iso-orbital indicator. To correct this, in this work we propose a new local scaling, LSIC-α, based on the iso-orbital indicator ασ (which additionally identifies regions of overlapping density tails). A two-parameter scaling function of ασ is fitted to a subset of the nonbonded appropriate norms for the SCAN and r2SCAN meta- GGAs, and tested on many properties of main-group atoms, molecules, and molecular complexes. LSIC-α greatly improves the interaction energies of weakly bonded systems in the S22 data set while retaining LSIC’s accuracy for other properties. This work shows that the errors of LSDA (and presumably of higher-level DFAs) can be largely but not entirely repaired by a proper “do no harm” self-interaction correction.

Approximation↗
DOE OSTI2023

Cosmological structure formation and soliton phase transition in fuzzy dark matter with axion self-interactions

Herein we investigate cosmological structure formation in fuzzy dark matter (FDM) with the attractive self-interaction (SI) with numerical simulations. Such a SI would arise if the FDM boson were an ultra-light axion, which has a strong CP symmetry-breaking scale (decay constant). Although weak, the attractive SI may be strong enough to counteract the quantum ‘pressure’ and alter structure formation. We find in our simulations that the SI can enhance small-scale structure formation, and soliton cores above a critical mass undergo a phase transition, transforming from dilute to dense solitons.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

A filter-dependent granular temperature model from large-scale CFD-DEM data

The computational study of strongly-coupled, gas–solid flows at scales relevant to most environmental and engineering applications requires the use of ‘coarse-grained’ methodologies such as the two-fluid model, particle-in-cell approach or the multiphase Reynolds Averaged Navier–Stokes equations. While these strategies enable computations at desirable length- and time-scales, they rely heavily on models to capture important flow physics that occur at scales smaller than the mesh. To date, the models that do exist are based on a limited set of flow conditions, such as very dilute particle phase. To this end, we leverage a large-scale repository of CFD-DEM data to develop filter-size dependent models for the mean variance in particle volume fraction, a quantity commonly used to assess the degree of clustering, and the granular temperature, a key quantity for accurately predicting gas–solid flows. In conclusion, because of its filter-size dependence, the granular temperature model can be directly translated to coarse-grained approaches and tied directly to grid size.

AMReX↗
DOE OSTI2020

Kelvin Waves during GOAmazon and Their Relationship to Deep Convection

The 2014–15 Observations and Modeling of the Green Ocean Amazon (GOAmazon) field campaign over the central Amazon near Manaus, Brazil, occurred in coordination with the larger Cloud Processes of the Main Precipitation Systems in Brazil: A Contribution to Cloud-Resolving Modeling and to the Global Precipitation Measurement (CHUVA) project across Brazil. These programs provide observations of convection over the central Amazon on diurnal to annual time scales. Here, we address the question of how Kelvin waves, observed in satellite observations of deep cloud cover over the GOAmazon region during the 2014–15 time period, modulate the growth, type, and organization of convection over the central Amazon. The answer to this question has implications for improved predictability of organized systems over the region and representation of convection and its growth on local to synoptic scales in global models. Our results demonstrate that Kelvin waves are strong modulators of synoptic-scale low- to midlevel free-tropospheric moisture, integrated moisture convergence, and surface heat fluxes. These regional modifications of the environment impact the local diurnal cycle of convection, favoring the development of mesoscale convective systems. As a result, localized rainfall is also strongly modulated, with the majority of rainfall in the GOAmazon region occurring during the passage of these systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Intense infrared lasers for strong-field science

The advent of chirped-pulse amplification in the 1980s and femtosecond Ti:sapphire lasers in the 1990s enabled transformative advances in intense laser–matter interaction physics. Whereas most of experiments have been conducted in the limited near-infrared range of 0.8–1 μm, theories predict that many physical phenomena such as high harmonic generation in gases favor long laser wavelengths in terms of extending the high-energy cutoff. Significant progress has been made in developing few-cycle, carrier-envelope phase-stabilized, high-peak-power lasers in the 1.6–2 μm range that has laid the foundation for attosecond X ray sources in the water window. Even longer wavelength lasers are becoming available that are suitable to study light filamentation, high harmonic generation, and laser–plasma interaction in the relativistic regime. Long-wavelength lasers are suitable for sub-bandgap strong-field excitation of a wide range of solid materials, including semiconductors. In the strong-field limit, bulk crystals also produce high-order harmonics. In this review, we first introduce several important wavelength scaling laws in strong-field physics, then describe recent breakthroughs in short- (1.4–3 μm), mid- (3–8 μm), and long-wave (8–15 μm) infrared laser technology, and finally provide examples of strong-field applications of these novel lasers. Furthermore, some of the broadband ultrafast infrared lasers will have profound effects on medicine, environmental protection, and national defense, because their wavelengths cover the water absorption band, the molecular fingerprint region, as well as the atmospheric infrared transparent window.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Pore-Scale Dynamics of Liquid CO2–Water Displacement in 2D Axisymmetric Porous Micromodels Under Strong Drainage and Weak Imbibition Conditions: High-Speed μPIV Measurements

Resolving pore-scale transient flow dynamics is crucial to understanding the physics underlying multiphase flow in porous media and informing large-scale predictive models. Surface properties of the porous matrix play an important role in controlling such physics, yet interfacial mechanisms remain poorly understood, in part due to a lack of direct observations. This study reports on an experimental investigation of the pore-scale flow dynamics of liquid CO 2 and water in two-dimensional (2D) circular porous micromodels with different surface characteristics employing high-speed microscopic particle image velocimetry (μPIV). The design of the micromodel minimized side boundary effects due to the limited size of the domain. The high-speed μPIV technique resolved the spatial and temporal dynamics of multiphase flow of CO 2 and water under reservoir-relevant conditions, for both drainage and imbibition scenarios. When CO 2 displaced water in a hydrophilic micromodel (i.e., drainage), unstable capillary fingering occurred and the pore flow was dominated by successive pore-scale burst events (i.e., Haines jumps). When the same experiment was repeated in a nearly neutral wetting micromodel (i.e., weak imbibition), flow instability and fluctuations were virtually eliminated, leading to a more compact displacement pattern. Energy balance analysis indicates that the conversion efficiency between surface energy and external work is less than 30%, and that kinetic energy is a disproportionately smaller contributor to the energy budget. This is true even during a Haines jump event, which induces velocities typically two orders of magnitude higher than the bulk velocity. These novel measurements further enabled direct observations of the meniscus displacement, revealing a significant alteration of the pore filling mechanisms during drainage and imbibition. While the former typically featured burst events, which often occur only at one of the several throats connecting a pore, the latter is typically dominated by a cooperative filling mechanism involving simultaneous invasion of a pore from multiple throats. This cooperative filling mechanism leads to merging of two interfaces and releases surface energy, causing instantaneous high-speed events that are similar, yet fundamentally different from, burst events. Finally, pore-scale velocity fields were statistically analyzed to provide a quantitative measure of the role of capillary effects in these pore flows.

Li, Yaofa↗
DOE OSTI2021

Magnetic Field Amplification by a Nonlinear Electron Streaming Instability

Magnetic field amplification by relativistic streaming plasma instabilities is central to a wide variety of high-energy astrophysical environments as well as to laboratory scenarios associated with intense lasers and electron beams. We report on a new secondary nonlinear instability which arises for relativistic dilute electron beams after the saturation of the linear Weibel instability. This instability grows due to the transverse magnetic pressure associated with the beam current filaments, which cannot be quickly neutralized due to the inertia of background ions. We show that it can amplify the magnetic field strength and spatial scale by orders of magnitude, leading to large-scale plasma cavities with strong magnetic field and to very efficient conversion of the beam kinetic energy into magnetic energy. The instability growth rate, saturation level, and scale length are derived analytically and shown to be in good agreement with fully-kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Spectroscopic Characterization of Highly Active Fe–N–C Oxygen Reduction Catalysts and Discovery of Strong Interaction with Nafion Ionomer

Scaling up clean-energy applications necessitates the development of platinum group metal (PGM)-free fuel cell electrocatalysts with high activity, stability, and low cost. Here, X-ray absorption (XAS) at the Fe K-edge and Fe K β X-ray emission (XES) spectroscopies were used to study the electronic structure of Fe centers in highly active Fe–N–C oxygen reduction catalysts with significant commercial potential. X-ray absorption near-edge structure (XANES) analysis has shown that the majority (>95%) of Fe centers are in the Fe 3+ oxidation state, while extended X-ray absorption fine structure (EXAFS) detected a mixture of single site Fe–N 4 centers (>95%) and centers with short (~2.5 Å) Fe–Fe interactions of Fe metal and/or Fe-carbide nanoparticles (<5%) featuring the Fe o oxidation state. Surprisingly, addition of Nafion, the most widely used ionomer, resulted in pronounced changes in the XAS spectra, consistent with a strong catalyst–ionomer interaction where long Fe–Fe interactions at ~3.1 Å were shown to be a feature of Fe 3+ ions bound with the Nafion. We conclude that exposure to Nafion during the device formulation has a different effect from the aggressive acid leaching typically used in the preparation of Fe–N–C catalysts. Furthermore, it was hypothesized that the polymer interacts with single sites’ Fe 3+ centers, as well as with graphene layers protecting the Fe o nanoparticles, and extracts some Fe ions into the Nafion matrix.

25 ENERGY STORAGE↗
DOE OSTI2022

Large-scale reactive thermoset printing: Complex interactions between temperature evolution, viscosity, and cure shrinkage

Thermoset composites are strong candidates for large-scale material extrusion additive manufacturing (AM). However, uncured thermoset composites become unstable as print height increases. Here, chemical initiation of vinyl ester immediately before deposition was used to suppress collapse and enable large-scale thermoset printing. Structural stability was assessed by printing thin walls at various layer times and monitoring thermal gradients with an infrared camera. Stable printing was observed at a layer time of 4.50 min, approximately half the gel time of the material (8 min). Self-weight collapse occurred at short layer time (0.68 min), whereas warpage occurred at long layer time (6.50 min). Here, a new behavior was discovered at intermediate layer time (2.25 min) where the heat generated by the reaction causes new, un-gelled layers to flow. Forced convection eliminates this behavior, stabilizing the 2.25-min layer time print. These findings motivated development of a less exothermic material system, which was used to print a large-scale mold and demonstrate the value of this understanding. By presenting these challenges of large-scale reactive thermoset AM for the first time, this work motivates focused studies of the complex interplay between rheological, thermal, and chemical behaviors to improve the feasibility of large-scale thermoset AM. Tensile properties of the printed material were also measured. Longitudinal and transverse elastic moduli are 3.79 and 2.95 GPa, respectively, and corresponding tensile strengths are 36.11 and 18.83 MPa. The glass transition temperature is 93.13 °C.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The ellipticity parameterization for an NFW profile: An overlooked angular structure in strong lens modeling

Galaxy-scale gravitational lenses are often modeled with two-component mass profiles where one component represents the stellar mass and the second is a Navarro Frenk White (NFW) profile representing the dark matter. Outside of the spherical case, the NFW profile is costly to implement, and so it is approximated via two different methods; ellipticity can be introduced via the lensing potential (NFWp) or via the mass by approximating the NFW profile as a sum of analytical profiles (NFWm). While the NFWp method has been the default for lensing applications, it gives a different prescription of the azimuthal structure, which we show introduces ubiquitous gradients in ellipticity and boxiness in the mass distribution rather than having a constant elliptical shape. Because an unmodeled azimuthal structure has been shown to be able to bias lens model results, we explored the degree to which this azimuthal structure that was introduced can affect the model accuracy. We constructed input profiles using composite models using both the NFWp and NFWm methods and fit these mocks with a power-law elliptical mass distribution (PEMD) model with external shear. As a measure of the accuracy of the recovered lensing potential, we calculated the value of the Hubble parameter H 0 one would determine from the lensing fit. We found that the fits to the NFWp input return H 0 values that are systematically biased by about 3% lower than the NFWm counterparts. We explored whether such an effect is attributable to the mass sheet transformation (MST) by using an MST-independent quantity, ξ 2 . We show that, as expected, the NFWm mocks are degenerate with PEMD through an MST. For the NFWp, an additional bias was found beyond the MST due to the azimuthal structure exterior to the Einstein radius. We recommend modelers use an NFWm prescription in the future, such that the azimuthal structure can be introduced explicitly rather than implicitly.

79 ASTRONOMY AND ASTROPHYSICS↗