Search NASA⌕ Search

SEARCH · Search NASA

Results for “fractional occupation number”

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.

Fractional occupation numbers and self-interaction correction-scaling methods with the Fermi-Löwdin orbital self-interaction correction approach

In this work, we present a new assessment of the Fermi-Löwdin orbital self-interaction correction (FLO-SIC) approach with an emphasis on its performance for predicting energies as a function of fractional occupation numbers (FONs) for various multielectron systems. Our approach is implemented in the massively parallelized NWChem quantum chemistry software package and has been benchmarked on the prediction of total energies, atomization energies, and ionization potentials of small molecules and relatively large aromatic systems. Within our study, we also derive an alternate expression for the FLO-SIC energy gradient expressed in terms of gradients of the Fermi-orbital eigenvalues and revisit how the FLO-SIC methodology can be seen as a constrained unitary transformation of the canonical Kohn–Sham orbitals. Finally, we conclude with calculations of energies as a function of FONs using various SIC-scaling methods to test the limits of the FLO-SIC formalism on a variety of multielectron systems. We find that these relatively simple scaling methods do improve the prediction of total energies of atomic systems as well as enhance the accuracy of energies as a function of FONs for other multielectron chemical species.

fractional occupation number↗

Mixed Precision Fermi-Operator Expansion on Tensor Cores from a Machine Learning Perspective

Here we present a second-order recursive Fermi-operator expansion scheme using mixed precision floating point operations to perform electronic structure calculations using tensor core units. A performance of over 100 teraFLOPs is achieved for half-precision floating point operations on Nvidia’s A100 tensor core units. The second-order recursive Fermi-operator scheme is formulated in terms of a generalized, differentiable deep neural network structure, which solves the quantum mechanical electronic structure problem. We demonstrate how this network can be accelerated by optimizing the weight and bias values to substantially reduce the number of layers required for convergence. We also show how this machine learning approach can be used to optimize the coefficients of the recursive Fermi-operator expansion to accurately represent the fractional occupation numbers of the electronic states at finite temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Density-Matrix Based Extended Lagrangian Born–Oppenheimer Molecular Dynamics

Extended Lagrangian Born–Oppenheimer molecular dynamics [ Phys. Rev. Lett. 2008, 100, 123004] is presented for Hartree–Fock theory, where the extended electronic degrees of freedom are represented by a density matrix, including fractional occupation numbers at elevated electronic temperatures. In contrast to regular direct Born–Oppenheimer molecular dynamics simulations, no iterative self-consistent field optimization is required prior to the force evaluations. To sample regions of the potential energy landscape where the gap is small or vanishing, which leads to particular convergence problems in regular direct Born–Oppenheimer molecular dynamics simulations, an adaptive integration scheme for the extended electronic degrees of freedom is presented. The integration scheme is based on a tunable, low-rank approximation of a fourth-order kernel, K, that determines the metric tensor, T ≡ K T K, used in the extended harmonic oscillator of the Lagrangian that generates the dynamics of the electronic degrees of freedom. Here, the formulation and algorithms provide a general guide to implement extended Lagrangian Born–Oppenheimer molecular dynamics for quantum chemistry, density functional theory, and semiempirical methods using a density matrix formalism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Symmetry Breaking with the SCAN Density Functional Describes Strong Correlation in the Singlet Carbon Dimer

The SCAN (strongly constrained and appropriately normed) meta-generalized gradient approximation (meta-GGA), which satisfies all 17 exact constraints that a meta-GGA can satisfy, accurately describes equilibrium bonds that are normally correlated. With symmetry breaking, it also accurately describes some sd equilibrium bonds that are strongly correlated. While sp equilibrium bonds are nearly always normally correlated, the C 2 singlet ground state is known from correlated wave function theory to be a rare case of strong correlation in an sp equilibrium bond. Earlier work that calculated atomization energies of the molecular sequence B 2 , C 2 , O 2 , and F 2 in the local spin density approximation (LSDA), the Perdew–Burke–Ernzerhof (PBE) GGA, and the SCAN meta-GGA, without symmetry breaking in the molecule, found that only SCAN was accurate enough to reveal an anomalous under-binding for C 2 . Here, this work shows that spin symmetry breaking in singlet C 2 , which involves the appearance of net up- and down-spin densities on opposite sides (not ends) of the bond, corrects that underbinding, with a small SCAN atomization-energy error more like that of the other three molecules, suggesting that symmetry breaking with an advanced density functional might reliably describe strong correlation. This article also discusses some general aspects of symmetry breaking and the insights into strong correlation that symmetry breaking can bring. The normally correlated low-lying triplet excited state has the right vertical excitation energy in SCAN but not in LSDA or PBE, where the triplet is a false ground state. Fractional occupation numbers are found only for the symmetry-unbroken singlet and only in LSDA and PBE GGA.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Graph-based quantum response theory and shadow Born–Oppenheimer molecular dynamics

Graph-based linear scaling electronic structure theory for quantum-mechanical molecular dynamics simulations [A. M. N. Niklasson et al., J. Chem. Phys. 144, 234101 (2016)] is adapted to the most recent shadow potential formulations of extended Lagrangian Born–Oppenheimer molecular dynamics, including fractional molecular-orbital occupation numbers [A. M. N. Niklasson, J. Chem. Phys. 152, 104103 (2020) and A. M. N. Niklasson, Eur. Phys. J. B 94, 164 (2021)], which enables stable simulations of sensitive complex chemical systems with unsteady charge solutions. The proposed formulation includes a preconditioned Krylov subspace approximation for the integration of the extended electronic degrees of freedom, which requires quantum response calculations for electronic states with fractional occupation numbers. For the response calculations, we introduce a graph-based canonical quantum perturbation theory that can be performed with the same natural parallelism and linear scaling complexity as the graph-based electronic structure calculations for the unperturbed ground state. Further, the proposed techniques are particularly well-suited for semi-empirical electronic structure theory, and the methods are demonstrated using self-consistent charge density-functional tight-binding theory both for the acceleration of self-consistent field calculations and for quantum-mechanical molecular dynamics simulations. Graph-based techniques combined with the semi-empirical theory enable stable simulations of large, complex chemical systems, including tens-of-thousands of atoms.

74 ATOMIC AND MOLECULAR PHYSICS↗

DESI mock challenge: Halo and galaxy catalogues with the bias assignment method

We present a novel approach to the construction of mock galaxy catalogues for large-scale structure analysis based on the distribution of dark matter halos obtained with effective bias models at the field level. We aim to produce mock galaxy catalogues capable of generating accurate covariance matrices for a number of cosmological probes that are expected to be measured in current and forthcoming galaxy redshift surveys (e.g. two- and three-point statistics). The construction of the catalogues shown in this paper is part of a mock-comparison project within the Dark Energy Spectroscopic Instrument (DESI) collaboration. We use the bias assignment method ( BAM ) to model the statistics of halo distribution through a learning algorithm using a few detailed N-body simulations, and approximated gravity solvers based on Lagrangian perturbation theory. We introduce cosmic-web-dependent corrections to modelling redshift-space distortions at the N-body level – both in the halo and galaxy distributions –, as well as a multi-scale approach for accurate assignment of halo properties. Using specific models of halo occupation distributions to populate halos, we generate galaxy mocks with the expected number density and central-satellite fraction of emission-line galaxies, which are a key target of the DESI experiment. BAM generates mock catalogues with per cent accuracy in a number of summary statistics, such as the abundance, the two- and three-point statistics of halo distributions, both in real and redshift space. In particular, the mock galaxy catalogues display ~3%-10% accuracy in the multipoles of the power spectrum up to scales of k ~ 0.4 h -1 Mpc. We show that covariance matrices of two- and three-point statistics obtained with BAM display a similar structure to the reference simulation. BAM offers an efficient way to produce mock halo catalogues with accurate two- and three-point statistics and is able to generate a variety of multi-tracer catalogues with precise covariance matrices of several cosmological probes. We discuss future developments of the algorithm towards mock production in DESI and other galaxy-redshift surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Imaging two-dimensional generalized Wigner crystals

The Wigner crystal has fascinated condensed matter physicists for nearly 90 years. Signatures of two-dimensional (2D) Wigner crystals were first observed in 2D electron gases under high magnetic field, and recently reported in transition metal dichalcogenide moiré superlattices. Direct observation of the 2D Wigner crystal lattice in real space, however, has remained an outstanding challenge. Conventional scanning tunnelling microscopy (STM) has sufficient spatial resolution but induces perturbations that can potentially alter this fragile state. Here we demonstrate real-space imaging of 2D Wigner crystals in WSe 2 /WS 2 moiré heterostructures using a specially designed non-invasive STM spectroscopy technique. This employs a graphene sensing layer held close to the WSe 2 /WS 2 moiré superlattice. Local STM tunnel current into the graphene layer is modulated by the underlying Wigner crystal electron lattice in the WSe 2 /WS 2 heterostructure. Different Wigner crystal lattice configurations at fractional electron fillings of n = 1/3, 1/2 and 2/3, where n is the electron number per site, are directly visualized. The n = 1/3 and n = 2/3 Wigner crystals exhibit triangular and honeycomb lattices, respectively, to minimize nearest-neighbour occupations. The n = 1/2 state spontaneously breaks the original C3 symmetry and forms a stripe phase. Here our study lays a solid foundation for understanding Wigner crystal states in WSe 2 /WS 2 moiré heterostructures and provides an approach that is generally applicable for imaging novel correlated electron lattices in other systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

How well do one-electron self-interaction-correction methods perform for systems with fractional electrons?

Recently developed locally scaled self-interaction correction (LSIC) is a one-electron SIC method that, when used with a ratio of kinetic energy densities (z σ ) as iso-orbital indicator, performs remarkably well for both thermochemical properties as well as for barrier heights overcoming the paradoxical behavior of the well-known Perdew–Zunger self-interaction correction (PZSIC) method. In this work, we examine how well the LSIC method performs for the delocalization error. Our results show that both LSIC and PZSIC methods correctly describe the dissociation of $H$$^{+}_{2}$ and $H$$^{+}_{2}$ but LSIC is overall more accurate than the PZSIC method. Likewise, in the case of the vertical ionization energy of an ensemble of isolated He atoms, the LSIC and PZSIC methods do not exhibit delocalization errors. For the fractional charges, both LSIC and PZSIC significantly reduce the deviation from linearity in the energy vs number of electrons curve, with PZSIC performing superior for C, Ne, and Ar atoms while for Kr they perform similarly. The LSIC performs well at the endpoints (integer occupations) while substantially reducing the deviation. The dissociation of LiF shows both LSIC and PZSIC dissociate into neutral Li and F but only LSIC exhibits charge transfer from Li + to F – at the expected distance from the experimental data and accurate ab initio data. Overall, both the PZSIC and LSIC methods reduce the delocalization errors substantially.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantum spin liquids bootstrapped from Ising criticality in Rydberg arrays

Arrays of Rydberg atoms constitute a highly tunable, strongly interacting venue for the pursuit of exotic states of matter. We develop a strategy for accessing a family of fractionalized phases known as quantum spin liquids in two-dimensional Rydberg arrays. We specifically use effective field theory methods to study arrays assembled from Rydberg chains tuned to an Ising phase transition that famously hosts emergent fermions propagating within each chain. This highly entangled starting point allows us to naturally access spin liquids familiar from Kitaev's honeycomb model—albeit from an entirely different framework. In particular, we argue that finite-range repulsive Rydberg interactions, which frustrate nearby symmetry-breaking orders, can enable coherent propagation of emergent fermions between the chains in which they were born. Delocalization of emergent fermions across the full two-dimensional Rydberg array yields a gapless Z 2 spin liquid with a single massless Dirac cone. Here, the Rydberg occupation numbers exhibit universal power-law correlations that provide a straightforward experimental diagnostic of this phase. We further show that explicitly breaking symmetries perturbs the gapless spin liquid into gapped, topologically ordered descendants: Breaking lattice symmetries generates toric-code topological order, whereas introducing Floquet-mediated chirality generates non-Abelian Ising topological order. In the toric-code phase, we analytically construct microscopic incarnations of non-Abelian defects, which can be created and transported by dynamically controlling the atom positions in the array. Our work suggests that appropriately tuned Rydberg arrays provide a cold-atoms counterpart of solid-state “Kitaev materials” and, more generally, it spotlights a different angle for pursuing experimental platforms for Abelian and non-Abelian fractionalization.

36 MATERIALS SCIENCE↗

Analysis and prediction of reaction kinetics using the degree of rate control

“Degree of rate control” (DRC) analysis provides a quantitative approach for analysing the kinetics of multi-step reaction mechanisms that has been widely applied to both heterogeneous and homogeneous catalysis research, as well as electrocatalysis. The DRC of any given transition state or intermediate is defined as a partial derivative such that it approximately equals the fractional increase in net rate to the product of interest per differential decrease in its standard-state free energy for that species (÷RT), holding constant the standard-state free energies of all other transition states and intermediates. Even very complex mechanisms usually have only a few species with non-zero DRCs and are thus the species whose interactions with the catalyst most strongly affect the net rate. These key DRC values thus offer a simple and intuitive route to optimize catalyst materials, especially with the assistance of computational methods like density functional theory (DFT). These high-DRC species are also the species whose energetics must be most accurately measured or calculated to achieve an accurate kinetic model for any reaction mechanism. In simple cases with a single “rate-determining step”, the DRC for its transition state (TS) is + 1. Catalyst-bound intermediates, on the other hand, often have negative DRCs equal to a small integer times their fractional occupancy of catalyst sites. The apparent activation energy equals a weighted average of the standard-state enthalpies (relative to reactants) of all the species (intermediates, transition states and products) in the reaction mechanism, each weighted by its DRC (+RT). It has been shown that the apparent transfer coefficient in electrocatalysis, an inverted form of the Tafel slope, is a weighted average of the number of electrons transferred to generate each intermediate or product species in the mechanism, weighted again by the DRC. Quantitative analysis of kinetic isotope effects (KIEs, or the ratio of net rates for different reactant isotopes) in complex mechanisms has shown that the logarithm of the KIE equals the weighted average over all species in the mechanism of the difference between the two isotopes in their standard-state free energies (÷RT), again weighted by the DRC. The reaction orders with respect to fluid-phase concentrations of reactants, products and intermediates have also been proven to be directly related to DRCs. Thus, there are numerous experimental observables which equate to short linear combinations of DRCs, so that combinations of experimental measurements might provide access to DRC values. Since its invention for transition states in 1994 and its generalization to include intermediates in 2009, the DRC has thus far mainly been calculated for microkinetic models based either entirely on DFT or on DFT with the key energies (i.e., those for high-DRC species) being fine-tuned to match experiments. The relationships summarized in this work provide new opportunities for using experiments earlier in the development and optimization of microkinetic models that require input from computational catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-redshift millennium and astrid galaxies in effective field theory at the field level

Effective field theory (EFT) modeling is expected to be a useful tool in the era of future higher-redshift galaxy surveys such as DESI-II and Spec-S5 due to its robust description of various large-scale structure tracers. However, large values of EFT bias parameters of higher-redshift galaxies could jeopardize the convergence of the perturbative expansion. Here, in this paper we measure the bias parameters and other EFT coefficients from samples of two types of star-forming galaxies in the state-of-the-art MilleniumTNG and astrid hydrodynamical simulations. Our measurements are based on the field-level EFT forward model that allows for precision EFT parameter measurements by virtue of cosmic variance cancellation. Specifically, we consider approximately representative samples of Lyman-break galaxies (LBGs) and Lyman-𝛼 emitters (LAEs) that are consistent with the observed (angular) clustering and number density of these galaxies at 𝑧 = 3. Reproducing the linear biases and number densities observed from existing LAE and LBG data, we find quadratic bias parameters that are roughly consistent with those predicted from the halo model coupled with a simple halo occupation distribution model. We also find nonperturbative velocity contributions (fingers of God) of a similar size for LBGs to the familiar case of luminous red galaxies. However, these contributions are quite small for LAEs despite their large satellite fraction values of up to ∼ 30%. Our results indicate that the effective momentum reach 𝑘 max at 𝑧 = 3 for LAEs (LBGs) will be in the range 0.3−0.6⁢ℎ Mpc −1 (0.2−0.8⁢ℎ Mpc −1 ), suggesting that EFT will perform well for high-redshift galaxy clustering. This work provides the first step toward obtaining realistic simulation-based priors on EFT parameters for LAEs and LBGs.

Sullivan, James M. [Massachusetts Inst. of Technol↗

The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: exploring the halo occupation distribution model for emission line galaxies

We study the modelling of the halo occupation distribution (HOD) for the eBOSS DR16 emission line galaxies (ELGs). Motivated by previous theoretical and observational studies, we consider different physical effects that can change how ELGs populate haloes. We explore the shape of the average HOD, the fraction of satellite galaxies, their probability distribution function (PDF), and their density and velocity profiles. Our baseline HOD shape was fitted to a semi-analytical model of galaxy formation and evolution, with a decaying occupation of central ELGs at high halo masses. We consider Poisson and sub/super-Poissonian PDFs for satellite assignment. We model both Navarro–Frenk–White and particle profiles for satellite positions, also allowing for decreased concentrations. We model velocities with the virial theorem and particle velocity distributions. Additionally, we introduce a velocity bias and a net infall velocity. We study how these choices impact the clustering statistics while keeping the number density and bias fixed to that from eBOSS ELGs. The projected correlation function, w p , captures most of the effects from the PDF and satellites profile. The quadrupole, ξ 2 , captures most of the effects coming from the velocity profile. We find that the impact of the mean HOD shape is subdominant relative to the rest of choices. Here, we fit the clustering of the eBOSS DR16 ELG data under different combinations of the above assumptions. The catalogues presented here have been analysed in companion papers, showing that eBOSS RSD+BAO measurements are insensitive to the details of galaxy physics considered here. These catalogues are made publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗

Technical Background and Validation Report on the Residential Water Inhalation Risk Calculator Presented in the Risk Assessment Information System

Indoor air quality (IAQ) is critical for human health. Poor IAQ is linked to respiratory issues, cardiovascular diseases, and cancer. Indoor pollutants are emitted by typical household items such as cleaning products, personal care items, building materials, and tap water - an understudied volatile organic compound (VOC) source. This document presents the Residential Water Inhalation Risk Calculator (RWIRC), which estimates daily VOC exposure concentrations from various household water uses, such as showering and dishwashing, to assess exposure risks for the most vulnerable occupant. Integrated into the Risk Assessment Information System (RAIS) and sponsored by the US Department of Energy (DOE), the calculator divides a house into three compartments: shower, bathroom, and other spaces, accounting for daily water usage patterns and calculating VOC concentrations. Exposure data generated using the calculator can assist health assessors in estimating excess lifetime cancer risk (ELCR) and hazard index (HI) from VOC inhalation. Unlike traditional exposure models that utilize Andelman’s constant, the RWIRC continuously assesses variability in VOC concentrations and environmental conditions using differential equations to track VOC concentrations and air exchange between compartments. The calculator also provides unique volatilization fractions for each chemical and appliance, enhancing accuracy of the exposure concentration estimation. The RWIRC is accessible online and allows users to customize parameters (i.e., number of bathrooms, water temperature, and exhaust fan conditions) and input VOC characteristics (i.e., tap water and ambient air concentrations). This document provides a step-by-step guide on implementing the calculator. It also provides comparisons with the ATSDR-SHOWER calculator, using eight VOCs with varying physicochemical properties to reveal differences in algorithms and output concentrations. Simulations also assess how bathroom door positions and exhaust fan usage affect VOC exposure. The calculator results can enhance EPA risk screening levels for inhalation exposure to VOCs from tap water, offering a sophisticated tool for assessing inhalation risks and improving public health protection.

54 ENVIRONMENTAL SCIENCES↗

Dirty bosons on the Cayley tree: Bose-Einstein condensation versus ergodicity breaking

Building on large-scale quantum Monte Carlo simulations, we investigate the zero-temperature phase diagram of hard-core bosons in a random potential on site-centered Cayley trees with branching number K=2. In order to follow how the Bose-Einstein condensate (BEC) is affected by the disorder, we focus on both the zero-momentum density, probing the quantum coherence, and the one-body density matrix (1BDM) whose largest eigenvalue monitors the off-diagonal long-range order. We further study its associated eigenstate which brings useful information about the real-space properties of this leading eigenmode. Upon increasing randomness, we find that the system undergoes a quantum phase transition at finite disorder strength between a long-range ordered BEC state, fully ergodic at large scale, and a new disordered Bose glass regime showing conventional localization for the coherence fraction while the 1BDM displays a nontrivial algebraic vanishing BEC density together with a nonergodic occupation in real space. These peculiar properties can be analytically captured by a simple phenomenological description on the Cayley tree which provides a physical picture of the Bose glass regime.

36 MATERIALS SCIENCE↗

Readout electronics for low occupancy High-Pressure Gas TPCs

High-Pressure Gas Time Projection Chambers (HPgTPCs) have benefits such as low energy thresholds, magnetisability, and 4π acceptance, making them ideal for neutrino experiments such as DUNE. We present the design of an FPGA-based solution optimised for Gaseous Argon Near Detector (ND-GAr), which is part of the Phase-II more capable near detector for DUNE. These electronics reduce the cost significantly compared to using collider readout electronics, which are typically designed for much higher occupancy and therefore, for example, need much larger numbers of FPGAs and power per channel. We demonstrate the performance of our electronics with the Teststand for an Overpressurised Argon Detector (TOAD) at Fermilab in the US at a range of pressures and gas mixtures up to 4.5 barA, reading out ∼10 000 channels from a Multi-Wire Proportional Chamber (MWPC). The operation took place between April and July of 2024. We measure the noise characteristics of the system to be sufficiently low, and we identify sources of noise that can be further mitigated in the next iteration. We also note that the cooling scheme used in the test requires improvement before full-scale deployment. Despite these necessary improvements, we show that the system can fulfil the needs of a HPgTPC for a fraction of the price of collider readout electronics.

Data acquisition concepts↗

On the Road to Achieving Safety Excellence at the Idaho Cleanup Project: The Initial Challenges, Solutions, and Positive Results of Striving to be the Best and Safest in the DOE Complex - 20518

When Fluor Idaho, LLC (Fluor Idaho) assumed the Idaho Cleanup Project Core contract (ICP Core) in mid-2016, it inherited the successful multi-faceted, employee-owned safety programs of two prior contractors supporting cleanup operations at the Idaho National Laboratory Site. Fluor Idaho incorporated the best of both safety programs along with Fluor Corporate's 106-year legacy of safety initiatives and successes. Fluor Idaho integrated the separate programs and cultures under a One Fluor approach - one project, one mission, One Fluor. Job changes, management changes, and procedures, processes, and system changes presented multiple challenges. Those challenges combined with the severe Winter of 2016/2017 resulted in numerous slips, trips, and falls first aids and reportable injuries. The road to safety excellence started off bumpy and stayed that way for some time. Since then, Fluor Idaho has successfully implemented several key initiatives including employee engagement; management engagement and time in the field; communications improvements; safety culture sustainment; injury/illness prevention; organizational learning; and DOE Voluntary Protection Program (VPP) readiness and preparation. Management became more engaged by making more workplace visits and increasing their supervision of 'life critical work' or complex activities. Additional training became mandatory for first-line supervisors and front-line workers. Most recently, Fluor Idaho increased its injury and illness prevention initiatives extensively. In the last year, things have turned around in a dramatic way. This shift is directly attributable to the significant number of initiatives launched in late 2017 and 2018. These initiatives addressed safety performance in the five key areas of the DOE VPP and promoted sustainability of the company's safety culture. The five key areas are: management leadership, employee involvement, work-site analysis, hazard prevention and control, and safety and health training. Employees and managers alike recommitted to making safety the highest value. The result of this laser-like focus on safety and accident and injury prevention is that injuries and first aids are a fraction of what they were. During FY 2019, employees worked more than 4.2 million hours without a serious injury or lost work-day and 1 million hours without a recordable injury. Fluor Idaho is proud of its safety trend. The company invested significant resources in improving its safety record, and that investment is paying off. Other contractors in the DOE Complex can learn and benefit from this discussion of the challenges Fluor Idaho faced, as well as its innovative safety initiatives that are producing results that are being recognized by our U.S. Department of Energy (DOE) customer. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Characterization of the Emissions and Crystalline Silica Content of Airborne Dust Generated from Grinding Natural and Engineered Stones

Abstract In this study, we systematically characterized the airborne dust generated from grinding engineered and natural stone products using a laboratory testing system designed and operated to collect representative respirable dust samples. Four stone samples tested included two engineered stones consisting of crystalline silica in a polyester resin matrix (formulations differed with Stones A having up to 90wt% crystalline silica and Stone B up to 50wt% crystalline silica), an engineered stone consisting of recycled glass in a cement matrix (Stone C), and a granite. Aerosol samples were collected by respirable dust samplers, total dust samplers, and a Micro-Orifice Uniform Deposit Impactor. Aerosol samples were analyzed by gravimetric analysis and x-ray diffraction to determine dust generation rates, crystalline silica generation rates, and crystalline silica content. Additionally, bulk dust settled on the floor of the testing system was analyzed for crystalline silica content. Real-time particle size distributions were measured using an Aerodynamic Particle Sizer. All stone types generated similar trimodal lognormal number-weighted particle size distributions during grinding with the most prominent mode at an aerodynamic diameter of about 2.0-2.3 μm, suggesting dust formation from grinding different stones is similar. Bulk dust from Stone C contained no crystalline silica. Bulk dust from Stone A, Stone B, and granite contained 60, 23, and 30wt% crystalline silica, respectively. In Stones A and B, the cristobalite form of crystalline silica was more plentiful than the quartz form. Only the quartz form was detected in granite. The bulk dust, respirable dust, and total dust for each stone had comparable amounts of crystalline silica, suggesting that crystalline silica content in the bulk dust could be representative of that in respirable dust generated during grinding. Granite generated more dust per unit volume of material removed than the engineered stones, which all had similar normalized dust generation rates. Stone A had the highest normalized generation rates of crystalline silica, followed by granite, Stone B, and Stone C (no crystalline silica), which likely leads to the same trend of respirable crystalline silica (RCS) exposure when working with these different stones. Manufacturing and adoption of engineered stone products with formulations such as Stone B or Stone C could potentially lower or eliminate RCS exposure risks. Combining all the effects of dust generation rate, size-dependent silica content, and respirable fraction, the highest normalized generation rate of RCS consistently occurs at 3.2-5.6 µm for all the stones containing crystalline silica. Therefore, removing particles in this size range near the generation sources should be prioritized when developing engineering control measures.

Public, Environmental & Occupational Health↗

How galaxies populate halos in extreme density environments: An analysis of the halo occupation distribution in SDSS

Recent works have shown that the properties of galaxy populations in simulated dark matter halos vary with large-scale environments. These results suggest a variation in the halo occupation distribution (HOD) in extreme density environments since the dynamical and astrophysical conditions prevailing in these regions may significantly affect the formation and evolution of their halos and residing galaxies, influencing the mean number of satellite galaxies. To analyse these effects, we identify cosmic voids and future virialised structures (FVS) in the Sloan Digital Sky Server Data Release 12 (SDSS-DR12) and estimate the HOD within these super-structures using group catalogues as dark matter halo proxies. Our goal is to use observational galaxy data to characterise the HOD within cosmic voids and FVS, explore the different properties of these galaxies populations and compare them with the general results outside these super-structures. We use a publicly available observational galaxy catalogue with information on redshifts, positions, magnitudes and other astrophysical features to build a volume complete galaxy sample and identify cosmic voids and FVS. Using a publicly available galaxy group catalogue as a proxy to dark matter halos, we compute the HOD within both types of super-structures for different absolute magnitude thresholds. We also study the dependence on the results on the main void and FVS properties, density and volume. We also analysed the main characteristics of the stellar content of galaxies inside these extreme-density regions such as the mean stellar age and the galaxy light concentration index. In all cases, we compare the results with those derived from the Field sample, defined by objects outside both types of environments. Inside cosmic voids, we find a strong decrease in HOD concerning the Field results. In the most extreme cases, the mean number of satellites fall to ~ 50%. Inside FVS, the HOD shows a significant increase to the Field, with a ~ 40% excess in the mean number of satellites. These results are present for the different galaxy luminosity ranges explored. In both environments, the differences with respect to the Field increases for the extreme values of the density environments. However, we obtain no signs of variations related to intrinsic characteristics of the super-structures, indicating that the effects depend mainly on the density of the large-scale environment. In addition, we find that the cumulative distribution of the mean age of stars of the central galaxy also varies in the different regions, this suggests that the history of the formation of the dark matter halos may be different. Finally, we explore the HOD for the 25% youngest (oldest) galaxies, based on the mean age of their stars. We find that for the low-mass groups the youngest galaxies are only present inside voids and are generally central galaxies. On the other hand, for the high-mass groups the FVS environments show the same increase in the HOD concerning the Field as previously mentioned. We find that cosmic voids lack a significant fraction of galaxies with the oldest stellar population.

79 ASTRONOMY AND ASTROPHYSICS↗