SEARCH · Engineering Papers

Results for “fractionalization”

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,315 records · Page 73

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Modeling Frequency–Independent Q Viscoacoustic Wave Propagation in Heterogeneous Media

Quantifying the attenuation of seismic waves propagating in the Earth interior is critical to study the subsurface structure. Previous studies have proposed fractional anelastic wave equations to model the frequency–independent Q seismic wave propagation. Such wave equations involve fractional derivatives that pose computational challenges for the numerical schemes in terms of accuracy and efficiency when dealing with heterogeneous Earth media. To tackle these challenges, here we derive a new viscoacoustic wave equation, where the power terms of the fractional Laplacian operators are spatially independent, thus accurate and efficient methods (e.g., the Fourier pseudospectral method) can be adopted. Our derivation enables the resultant equation to capture both amplitude and phase signatures of the anelastic wave propagation by matching the complex wave numbers for all the frequencies of interest. We verify the derivation by comparing the dispersion curves of both the attenuation factor and the phase velocity produced by the new wave equation with their theoretical values as well as the Pierre Shale in situ measurements. Following that, we use a synthetic attenuating gas chimney model to demonstrate the attenuation effects on seismic waveforms and then construct the Q–compensated reverse time migration to undo these effects for seismic image enhancement. Finally, we find that our forward modeling results can characterize the spatiotemporal attenuation effects revealed in the Frio–II CO2 injection time–lapse seismic monitoring data. Here, we expect this proposed equation to be useful to quantify the attenuation in seismic data to push the resolution limits of seismic imaging and inversion.

25 ENERGY STORAGE↗

DOE OSTI2020

Wildfire‐Smoke Aerosols Lead to Increased Light Use Efficiency Among Agricultural and Restored Wetland Land Uses in California's Central Valley

Abstract There are few observational studies measuring the ecosystem‐scale productivity effects of changes in incident diffuse photosynthetically active radiation (PAR diffuse ), especially related to wildfire smoke. Climate change‐induced increases to the duration and intensity of fire conditions have made smoke a common occurrence across western North America, with largely unquantified ecosystem feedbacks. Under equivalent amounts of radiation, increased atmospheric particulate matter could lead to a boost in productivity as scattering redistributes photons throughout multilayer canopies. In this work, we leverage a meso‐network of eddy covariance measurement sites across a unique array of managed and restored C 3 and C 4 canopy types to understand how recent wildfire smoke affected ecosystem productivity during the summer of 2018, an especially smoky year in the agriculturally productive Central Valley. We find that diffuse PAR diffuse increased by more than a third compared to the previous growing season, while total PAR was only slightly diminished. These conditions caused nearly a doubling of light use efficiency over the range of diffuse fraction observed, with the highest sensitivity to diffuse fraction exhibited by corn and alfalfa crops. We utilized an empirical model to assess the trade‐off between enhanced diffuse fraction and reduced total PAR. Under mean radiation conditions, daily integrated gross ecosystem productivity increased by 1.2–4.2% compared to the previous growing season. Finally, we explore the potential negative effect of heightened ozone, a copollutant often associated with wildfire. In addition to the effects of wildfire smoke, the results of this natural experiment can help validate future predictions of aerosol‐productivity feedbacks.

Hemes, Kyle S.↗

DOE OSTI2022

Load‐Balancing Intense Physics Calculations to Embed Regionalized High‐Resolution Cloud Resolving Models in the E3SM and CESM Climate Models

Abstract We design a new strategy to load‐balance high‐intensity sub‐grid atmospheric physics calculations restricted to a small fraction of a global climate simulation's domain. We show why the current parallel load balancing infrastructure of Community Earth System Model (CESM) and Energy Exascale Earth Model (E3SM) cannot efficiently handle this scenario at large core counts. As an example, we study an unusual configuration of the E3SM Multiscale Modeling Framework (MMF) that embeds a binary mixture of two separate cloud‐resolving model grid structures that is attractive for low cloud feedback studies. Less than a third of the planet uses high‐resolution (MMF‐HR; sub‐km horizontal grid spacing) relative to standard low‐resolution (MMF‐LR) cloud superparameterization elsewhere. To enable MMF runs with Multi‐Domain cloud resolving models (CRMs), our load balancing theory predicts the most efficient computational scale as a function of the high‐intensity work's relative overhead and its fractional coverage. The scheme successfully maximizes model throughput and minimizes model cost relative to precursor infrastructure, effectively by devoting the vast majority of the processor pool to operate on the few high‐intensity (and rate‐limiting) high‐resolution (HR) grid columns. Two examples prove the concept, showing that minor artifacts can be introduced near the HR/low‐resolution CRM grid transition boundary on idealized aquaplanets, but are minimal in operationally relevant real‐geography settings. As intended, within the high (low) resolution area, our Multi‐Domain CRM simulations exhibit cloud fraction and shortwave reflection convergent to standard baseline tests that use globally homogenous MMF‐LR and MMF‐HR. We suggest this approach can open up a range of creative multi‐resolution climate experiments without requiring unduly large allocations of computational resources.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Impacts of Synoptic-Scale Dynamics on Clouds and Radiation in High Southern Latitudes

High-latitudinal mixed-phase clouds significantly affect Earth's radiative balance. Observations of cloud and radiative properties from two field campaigns in the Southern Ocean and Antarctica were compared with two global climate model simulations. A cyclone compositing method was used to quantify “dynamics-cloud-radiation” relationships relative to the extratropical cyclone centers. Observations show larger asymmetry in cloud and radiative properties between western and eastern sectors at McMurdo compared with Macquarie Island. Most observed quantities at McMurdo are higher in the western (i.e., post-frontal) than the eastern (frontal) sector, including cloud fraction, liquid water path (LWP), net surface shortwave and longwave radiation (SW and LW), except for ice water path (IWP) being higher in the eastern sector. Here, the two models were found to overestimate cloud fraction and LWP at Macquarie Island but underestimate them at McMurdo Station. IWP is consistently underestimated at both locations, both sectors, and in all seasons. Biases of cloud fraction, LWP, and IWP are negatively correlated with SW biases and positively correlated with LW biases. The persistent negative IWP biases may have become one of the leading causes of radiative biases over the high southern latitudes, after correcting the underestimation of supercooled liquid water in the older model versions. By examining multi-scale factors from cloud microphysics to synoptic dynamics, this work will help increase the fidelity of climate simulations in this remote region.

58 GEOSCIENCES↗

DOE OSTI2025

Glacial‐Interglacial and Millennial‐Scale Changes in Nitrous Oxide Emissions Pathways and Source Regions

During the transition from the Last Glacial Maximum (LGM) to the Holocene, the atmospheric N 2 O mole fraction increased by 80 nmol mol −1 . Using ice core measurements of N 2 O isotopomer ratios, we show that this increase was driven by increases in both nitrification and denitrification, with the relative partitioning between both production pathways depending on the assumed isotopic end-member source signatures. Similarly, we also attribute a 35 nmol mol −1 N 2 O mole fraction increase during the Heinrich Stadial 4/Dansgaard Oeschger 8 (HS4/DO8) millennial-scale event to increases in both N 2 O production pathways. In contrast, the 25 nmol mol −1 N 2 O mole fraction decrease during the Younger Dryas was driven almost exclusively by a decrease in nitrification. The deglacial and HS4/DO8 increases in N 2 O production occurred in both marine and terrestrial environments, with the terrestrial source responding faster to warming by about two centuries. Constraints on changes in nitrification and denitrification emissions are robust and consistent with previous studies showing the sensitivity of N 2 O emissions to abrupt Northern Hemisphere warming. This study demonstrates for the first time the importance of both denitrification and nitrification pathways in driving source changes. Absolute emissions are more uncertain due to uncertainty about source isotopomer signatures. For instance, the contribution of denitrification to emissions at the LGM shifts from (65 ± 10) % to (91 ± 6) % when factoring in isotope enrichment due to partial reduction of N 2 O to N 2 during denitrification. Reducing uncertainty in source signatures will increase the power of ice core N 2 O isotope records in deducing environmental change.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Light Transfers Through a Koch Shape Cloud

Abstract Modeling radiative transfer in a 3D cloudy atmosphere is critical to climate projections. A recently developed fast 3D radiation parameterization scheme gains some success in quantifying horizontal radiative transfer through cloud sides using cloud area fraction. Based on 3D Monte Carlo simulations of radiative transfer through an idealized single‐layer cloud with Koch‐shaped fractal geometry edges, here we show that radiative energy transport through cloud sides correlates more significantly with cloud area fraction than with cloud perimeter length. The results exemplify the importance of accounting for the horizontal radiative energy exchanges between cloud‐free and cloudy regions with cloud area fraction. Results from additional sensitivity simulations show that increased cloud vertical extent often enhances cloud‐side sunlight leak more significantly than cloud‐side sunlight interception. At low sun elevations, cloud‐side sunlight interception is enhanced more than cloud‐side sunlight leak does with the increase of cloud mass.

58 GEOSCIENCES↗

DOE OSTI2026

Large Eddy Simulations of the Atmospheric Boundary Layer Over Satellite‐Sensed Sea Ice Maps

Surface heterogeneity in the marginal ice zone (MIZ) causes multiscale secondary atmospheric circulations that are challenging to model or observe. The absence or inadequate representation of these circulations in ocean‐atmosphere exchange schemes in climate models is partially responsible for the underestimation of Arctic sea ice loss. Observationally, such circulations obfuscate the interpretation of polar atmospheric chemistry measurements, among others. To address this open challenge, large‐eddy simulations are conducted over real‐world satellite‐sensed sea ice maps with an ice‐sea temperature contrast, as well as over idealizations of these maps that alter the ice pattern but conserve its fraction, showing that the ice fraction of a surface is not sufficient to predict the surface heat flux. In a second suite of simulations, three other heterogeneity metrics (representing the surface fragmentation, patch edge tortuosity, and patch size variability) are introduced to complement the ice fraction in describing the surface. Simulations varying these parameters suggest that they also significantly impact surface‐air interactions. A roughness contrast is then added to the surface temperature contrast, indicating that the contribution of roughness changes to the resulting atmospheric circulations is less pronounced than that of thermal heterogeneity. Based on these results, we illustrate, using a multi‐linear regression on these map features, that MIZ surface parameterizations in Earth Systems Models can be improved if they account for these various characteristics of the sea‐ice patterns.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Physics or Knob‐Tuning? Tropical Anvil Peak Is Captured by GCMs

Tropical anvil clouds peak near 200 hPa and significantly impact Earth's climate, yet its physical realism in coarse‐resolution General Circulation Models (GCMs) remains debated. We examine anvil cloud formation by performing simulations with a GCM with a hierarchy of cloud fraction schemes ranging from a complex prognostic Tiedtke scheme to a simple binary scheme. All schemes consistently reproduce the anvil peak. The robust anvil peak arises because extremely cold temperatures at the upper troposphere facilitate frequent saturation events, producing clouds that disproportionately influence mean cloud fraction. Sensitivity experiments with enhanced evaporation of cloud condensate unexpectedly show increased anvil coverage, highlighting how slight evaporative moistening reinforces local saturation in cold upper‐tropospheric conditions. These results demonstrate that the tropical anvil cloud peak emerges from fundamental thermodynamic constraints, rather than specific cloud fraction parameterization choices.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Giant phonon anomalies in the proximate Kitaev quantum spin liquid α-RuCl 3

The Kitaev quantum spin liquid epitomizes an entangled topological state, for which two flavors of fractionalized low-energy excitations are predicted: the itinerant Majorana fermion and the Z 2 gauge flux. Here, it was proposed recently that fingerprints of fractional excitations are encoded in the phonon spectra of Kitaev quantum spin liquids through a novel fractional-excitation-phonon coupling. Here, we detect anomalous phonon effects in α-RuCl 3 using inelastic X-ray scattering with meV resolution. At high temperature, we discover interlaced optical phonons intercepting a transverse acoustic phonon between 3 and 7 meV. Upon decreasing temperature, the optical phonons display a large intensity enhancement near the Kitaev energy, J K ~8 meV, that coincides with a giant acoustic phonon softening near the Z 2 gauge flux energy scale. These phonon anomalies signify the coupling of phonon and Kitaev magnetic excitations in α-RuCl 3 and demonstrates a proof-of-principle method to detect anomalous excitations in topological quantum materials.

36 MATERIALS SCIENCE↗

DOE OSTI2025

PAVC: The foundation for a Pan-Arctic Vegetation Cover database

Field-measured Arctic vegetation cover data is essential for creating accurate, high-quality vegetation structure and composition maps. Extrapolating field data into high-resolution cover maps provides detailed, function-specific information for use in Earth System Models, vegetation classifications, and monitoring vegetation change over time and space. However, field campaigns that collect plant cover vary substantially in scope, method, and purpose, which makes them difficult to unify across data stores, and they are often not designed to meet remote sensing needs. In this work, we synthesized and harmonized field-based fractional cover data from various data stores to create a high-quality, consistent repository schema for remote sensing-based vegetation cover mapping applications. We developed a reproducible workflow for synthesizing visual estimate and point-intercept fractional cover data. The resultant Pan-Arctic Vegetation Cover (PAVC) database contains synthesized fractional cover at both the species and plant functional type levels. The latter includes absolute foliar cover for deciduous shrubs and trees, evergreen shrubs and trees, forbs, graminoids, lichen, bryophytes, and “other” vegetation, as well as absolute cover for litter and top cover for water and bare ground.

Steckler, Morgan R. [Oak Ridge National Laboratory↗

DOE OSTI2021

In situ synchrotron diffraction and modeling of non-equilibrium solidification of a MnFeCoNiCu alloy

The solidification mechanism and segregation behavior of laser-melted Mn 35 Fe 5 Co 20 Ni 20 Cu 20 was firstly investigated via in situ synchrotron x-ray diffraction at millisecond temporal resolution. The transient composition evolution of the random solid solution during sequential solidification of dendritic and interdendritic regions complicates the analysis of synchrotron diffraction data via any single conventional tool, such as Rietveld refinement. Therefore, a novel approach combining a hard-sphere approximation model, thermodynamic simulation, thermal expansion measurement and microstructural characterization was developed to assist in a fundamental understanding of the evolution of local composition, lattice parameter, and dendrite volume fraction corresponding to the diffraction data. This methodology yields self-consistent results across different methods. Via this approach, four distinct stages were identified, including: (I) FCC dendrite solidification, (II) solidification of FCC interdendritic region, (III) solid-state interdiffusion and (IV) final cooling with marginal diffusion. It was found out that in Stage I, Cu and Mn were rejected into liquid as Mn 35 Fe 5 Co 20 Ni 20 Cu 20 solidified dendritically. During Stage II, the lattice parameter disparity between dendrite and interdendritic region escalated as Cu and Mn continued segregating into the interdendritic region. After complete solidification, during Stage III, the lattice parameter disparity gradually decreases, demonstrating a degree of composition homogenization. The volume fraction of dendrites slightly grew from 58.3 to 65.5%, based on the evolving composition profile across a dendrite/interdendritic interface in diffusion calculations. Postmortem metallography further confirmed that dendrites have a volume fraction of 64.7 ± 5.3% in the final microstructure.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Crystallization kinetics of rapid spray plasma processed multiple cation perovskites in open air

Rapid Spray Plasma Processing (RSPP) is a high throughput, scalable, and open-air route toward manufacturing perovskite solar modules. The plasma exposure is dosed such that the perovskite precursor solution is cured in milliseconds using a combination of rapid heating, reactive species, and UV photons. We implemented in situ wide angle X-ray scattering (WAXS) to characterize the crystal growth in multiple cation, mixed halide perovskite thin films—Cs x (MA z FA 1-z ) 1-x Pb(Br 0.17 I 0.83 ) 3 for 0.05 ≤ x ≤ 0.25—fabricated in open air. The use of synchrotron radiation and high detector speeds resolved the formation and dissolution of a transient intermediate crystalline phase during the tens of milliseconds when the perovskite crystallized. Increasing the mole fraction of cesium and methylammonium resulted in a decrease in the intermediate phase and an increase in the average grain diameter. Increasing the cesium fraction suppressed phase segregation, as observed in photoluminescence. Calculations showed the refractive index decreased and the extinction coefficient increased with cesium fraction. Based on the rapid crystallization kinetics, RSPP perovskite films exhibit residual stress values >5× lower than spin coated films. Understanding the kinetics of perovskite formation enables the design of perovskite films with improved thermomechanical and operational stability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Multicolor polymeric carbon dots: synthesis, separation and polyamide-supported molecular fluorescence

Multicolor carbon dots (CDs) have been developed recently and demonstrate great potential in bio-imaging, sensing, and LEDs. However, the fluorescence mechanism of their tunable colors is still under debate, and efficient separation methods are still challenging. Herein, we synthesized multicolor polymeric CDs through solvothermal treatment of citric acid and urea in formamide. Automated reversed-phase column separation was used to achieve fractions with distinct colors, including blue, cyan, green, yellow, orange and red. This work explores the physicochemical properties and fluorescence origins of the red, green, and blue fractions in depth with combined experimental and computational methods. Three dominant fluorescence mechanism hypotheses were evaluated by comparing time-dependent density functional theory and molecular dynamics calculation results to measured characteristics. We find that blue fluorescence likely comes from embedded small molecules trapped in carbonaceous cages, while pyrene analogs are the most likely origin for emission at other wavelengths, especially in the red. Also important, upon interaction with live cells, different CD color fractions are trafficked to different sub-cellular locations. Super-resolution imaging shows that the blue CDs were found in a variety of organelles, such as mitochondria and lysosomes, while the red CDs were primarily localized in lysosomes. These findings significantly advance our understanding of the photoluminescence mechanism of multicolor CDs and help to guide future design and applications of these promising nanomaterials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Controlled site coverage of strong metal–support interaction (SMSI) on Pd NP catalysts

Here, strong metal–support interaction catalysts have been shown to improve desired product selectivity at the cost of fractional rates due to active site coverage. The goal of this study was to determine if the active site coverage of metallic nanoparticles could be controlled to lower levels than have been previously reported in SMSI catalysts with the aim of improving the rate while maintaining high selectivity. 2Pd–XTi/SiO 2 (2 wt% Pd, X wt% Ti) strong metal–support interaction (SMSI) catalysts with Ti loadings between 0–1.0 wt% were synthesized to control Pd nanoparticle coverage. Calcination at 450 °C and reduction at 550 °C were sufficient for forming ~2 nm sized Pd particles in all catalysts. Increasing the Ti loading from 0.1 to 1.0 wt% increased the surface coverage from 40 to 85% at a fixed reduction temperature of 550 °C. The IR spectra of the SMSI catalysts were similar with a high fraction of linear bonded CO which was much higher than that of Pd nanoparticles of similar size. The SMSI overlayer could be removed by oxidation at 350 °C and re-reduction at 200 °C. EXAFS of the oxidized catalysts indicates that nearly full oxidation of the metallic nanoparticle was required to remove the SMSI overlayer. Oxidation temperatures from 30 to 300 °C partially oxidized the Pd nanoparticles and subsequent re-reduction at 200 °C partially decreases the SMSI coverage. The fractional surface coverage was determined by measuring the rate of propylene hydrogenation with and without the SMSI overlayer. Increasing the reduction temperature from 200 to 550 °C increased the SMSI coverage from 0 to 85% depending on the Ti loading and temperature. In conclusion, after reduction at 550 °C and oxidation at 350 °C, the range of coverages varied between ~10% with 0.1 wt% Ti after re-reduction at 300 °C and ~85% with 1 wt% Ti after reduction at 550 °C.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Classroom aerosol dispersion modeling: experimental assessment of a low-cost flow simulation tool

The purpose of this study was to assess the utility of a low-cost flow simulation tool for an indoor air modeling application by comparing its outputs with the results of a physical experiment, as well as those from a more advanced computational fluid dynamics (CFD) software package. In this study, five aerosol dispersion tests were performed in two different classrooms by releasing a CO 2 tracer gas from six student locations. Resultant steady-state concentrations were monitored at 13 locations around the periphery of the room. Subsequently, the experiments were modeled using both a low-cost tool (SolidWorks Flow Simulation) and a more sophisticated tool (STAR-CCM+). Models were evaluated based on their ability to predict the experimentally measured concentrations at the 13 monitoring locations by calculating four performance parameters commonly used in the evaluation of dispersion models: fractional mean bias (FB), normalized mean-square error (NMSE), fraction of predicted value within a factor of two (FAC2), and normalized absolute difference (NAD). The more sophisticated model performed better in 15 of the 20 possible cases (five tests at four parameters each), with parameters meeting acceptance criteria in 19 of 20 cases. However, the lower-cost tool was only slightly worse, with parameters meeting acceptance criteria in 18 of 20 cases, and it performed better than the other tool in 3 of 20 cases. Because it provides useful results at a fraction of the monetary and training cost and is already widely accessible to many institutions, such a tool may be worthwhile for many indoor aerosol dispersion applications, especially for students or researchers just beginning CFD modeling.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Multi-scale computational screening and mechanistic insights of cyclic amines as solvents for improved lignocellulosic biomass processing

A computational screening workflow for the efficient deconstruction of cellulose, lignin and hemicellulose fractions of lignocellulosic biomass using cyclic amines as solvents. Lignocellulosic biomass is a promising feedstock for production of affordable fuels and chemicals from renewable resources. Effective solubilization and subsequent deconstruction of its cellulose, hemicellulose, and lignin fractions is essential for the viability of future biorefineries. This study used quantum chemistry-based equilibrium thermodynamics methods to evaluate the potential of 650 cyclic amines to solubilize cellulose, hemicellulose, and lignin. The activity coefficients of solvent - biopolymer interactions were predicted using the COSMO-RS (COnductor-like Screening MOdel for Real Solvents) method and used to identify cyclic amines that can efficiently dissolve and extract selective fractions of biopolymers during biomass pretreatment. Among the 650 cyclic amines, 1-piperazineethanmaine was predicted to be an effective solvent for extracting all three polymers and was experimentally shown to achieve the highest lignin removal (97.1%). Non-covalent interaction, reduced density gradient and quantum chemical calculations were performed to elucidate the dissolution mechanism of lignin, cellulose and hemicellulose and gain further molecular level insights into the interactions between the cyclic amines and biomass polymers that promote efficient solubilization and extraction. These analyses indicated that 1-piperazineethanmaine and 1-methylimidazole make noncovalent van der Waals, electrostatic interactions and hydrogen bonding with lignin, leading to enhanced lignin removal, while the strong intramolecular hydrogen bonding interactions in cellulose and hemicellulose result in weaker solvent-biopolymer interactions. Overall, the computational approach provided an efficient method for identifying cyclic amines tailored for optimal biomass pretreatment and resulted in the identification of a potential new class of solvents for effective biomass pretreatment.

Kumar, Nikhil↗

DOE OSTI2020

Dust and gas absorption in the high mass X-ray binary IGR J16318-4848

With an absorption column density on the order of 10 24 cm -2 , IGR J16318-4848 is one of the most extreme cases of a highly obscured high mass X-ray binary. In addition to the overall continuum absorption, the source spectrum exhibits a strong iron and nickel fluorescence line complex at 6.4 keV. Previous empirical modeling of these features and comparison with radiative transfer simulations raised questions about the structure and covering fraction of the absorber and the profile of the fluorescence lines. We aim at a self-consistent description of the continuum absorption, the absorption edges, and the fluorescence lines to constrain the properties of the absorbing material, such as ionization structure and geometry. We further investigate the effects of dust absorption on the observed spectra and the possibility of fluorescence emission from dust grains. We used XMM-Newton and NuSTAR spectra to first empirically constrain the incident continuum and fluorescence lines. Next we used XSTAR to construct a customized photoionization model where we vary the ionization parameter, column density, and covering fraction. In the third step we modeled the absorption and fluorescence in a dusty olivine absorber and employed both a simple analytical model for the fluorescence line emission and a Monte Carlo simulation of radiative transfer that generates line fluxes, which are very close to the observational data. Our empirical spectral modeling is in agreement with previous works. Our second model, the single gas absorber does not describe the observational data. In particular, irrespective of the ionization state or column density of the absorber, a much higher covering fraction than previously estimated is needed to produce the strong fluorescence lines and the large continuum absorption. A dusty, spherical absorber (modeled as consisting of olivine dust, although the nature of dust cannot be constrained) is able to produce the observed continuum absorption and edges. A dense, dusty absorber in the direct vicinity of the source consisting of dust offers a consistent description of both the strong continuum absorption and the strong emission features in the X-ray spectrum of IGR J16318-4848. In particular, for low optical depth of individual grains, which is the case for typical volume densities and grain size distribution models, the dust will contribute significantly to the fluorescence emission.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

J-PLUS: Spectral evolution of white dwarfs by PDF analysis

We estimated the spectral evolution of white dwarfs with effective temperature using the Javalambre Photometric Local Universe Survey (J-PLUS) second data release (DR2), which provides 12 photometric optical passbands over 2176 deg 2 . We analyzed 5926 white dwarfs with r ≤ 19.5 mag in common between a white dwarf catalog defined from Gaia EDR3 and J-PLUS DR2. We performed a Bayesian analysis by comparing the observed J-PLUS photometry with theoretical models of hydrogen- and helium-dominated atmospheres. We estimated the probability distribution functions for effective temperature (T eff ), surface gravity, parallax, and composition; and the probability of having a H-dominated atmosphere (p H ) for each source. We applied a prior in parallax, using Gaia EDR3 measurements as a reference, and derived a self-consistent prior for the atmospheric composition as a function of T eff . We described the fraction of white dwarfs with a He-dominated atmosphere (f He ) with a linear function of the effective temperature at 5000 < T eff < 30 000 K. We find f He = 0.24 ± 0.01 at T eff = 10 000 K, a change rate along the cooling sequence of 0.14 ± 0.02 per 10 kK, and a minimum He-dominated fraction of 0.08 ± 0.02 at the high-temperature end. We tested the obtained p H by comparison with spectroscopic classifications, finding that it is reliable. We estimated the mass distribution for the 351 sources with distance d < 100 pc, mass M > 0.45 M ⊙ , and T eff > 6000 K. The result for H-dominated white dwarfs agrees with previous studies, with a dominant M = 0.59 M ⊙ peak and the presence of an excess at M ~ 0.8 M ⊙ . This high-mass excess is absent in the He-dominated distribution, which presents a single peak. The J-PLUS optical data provide a reliable statistical classification of white dwarfs into H- and He-dominated atmospheres. We find a 21 ± 3% increase in the fraction of He-dominated white dwarfs from T eff = 20 000 K to T eff = 5000 K.

79 ASTRONOMY AND ASTROPHYSICS↗