SEARCH · Search NASA

Results for “model structural error”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 1,225 records · Page 68

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Extraction of mechanical properties of materials through deep learning from instrumented indentation

Instrumented indentation has been developed and widely utilized as one of the most versatile and practical means of extracting mechanical properties of materials. This method is particularly desirable for those applications where it is difficult to experimentally determine the mechanical properties using stress–strain data obtained from coupon specimens. Such applications include material processing and manufacturing of small and large engineering components and structures involving the following: three dimensional (3D) printing, thin-film and multilayered structures, and integrated manufacturing of materials for coupled mechanical and functional properties. Here, we utilize the latest developments in neural networks, including a multifidelity approach whereby deep-learning algorithms are trained to extract elastoplastic properties of metals an d alloys from instrumented indentation results using multiple datasets for desired levels of improved accuracy. We have established algorithms for solving inverse problems by recourse to single, dual, and multiple indentation and demonstrate that these algorithms significantly outperform traditional brute force computations and function-fitting methods. Moreover, we present several multifidelity approaches specifically for solving the inverse indentation problem which 1) significantly reduce the number of high-fidelity datasets required to achieve a given level of accuracy, 2) utilize known physical and scaling laws to improve training efficiency and accuracy, and 3) integrate simulation and experimental data for training disparate datasets to learn and minimize systematic errors. The predictive capabilities and advantages of these multifidelity methods have been assessed by direct comparisons with experimental results for indentation for different commercial alloys, including two wrought aluminum alloys and several 3D printed titanium alloys.

3D printed materials, stress-strain behavior, mult↗

DOE OSTI2018

The Role of Proton Transport in Gating Current in a Voltage Gated Ion Channel, as Shown by Quantum Calculations

Over two-thirds of a century ago, Hodgkin and Huxley proposed the existence of voltage gated ion channels (VGICs) to carry Na + and K + ions across the cell membrane to create the nerve impulse, in response to depolarization of the membrane. The channels have multiple physiological roles, and play a central role in a wide variety of diseases when they malfunction. The first channel structure was found by MacKinnon and coworkers in 1998. Subsequently, the structure of a number of VGICs was determined in the open (ion conducting) state. This type of channel consists of four voltage sensing domains (VSDs), each formed from four transmembrane (TM) segments, plus a pore domain through which ions move. Understanding the gating mechanism (how the channel opens and closes) requires structures. One TM segment (S4) has an arginine in every third position, with one such segment per domain. It is usually assumed that these arginines are all ionized, and in the resting state are held toward the intracellular side of the membrane by voltage across the membrane. They are assumed to move outward (extracellular direction) when released by depolarization of this voltage, producing a capacitive gating current and opening the channel. We suggest alternate interpretations of the evidence that led to these models. Measured gating current is the total charge displacement of all atoms in the VSD; we propose that the prime, but not sole, contributor is proton motion, not displacement of the charges on the arginines of S4. It is known that the VSD can conduct protons. Quantum calculations on the K v 1.2 potassium channel VSD show how; the key is the amphoteric nature of the arginine side chain, which allows it to transfer a proton. This appears to be the first time the arginine side chain has had its amphoteric character considered. We have calculated one such proton transfer in detail: this proton starts from a tyrosine that can ionize, transferring to the NE of the third arginine on S4; that arginine’s NH then transfers a proton to a glutamate. The backbone remains static. A mutation predicted to affect the proton transfer has been qualitatively confirmed experimentally, from the change in the gating current-voltage curve. The total charge displacement in going from a normal closed potential of -70 mV across the membrane to 0 mV (open), is calculated to be approximately consistent with measured values, although the error limits on the calculation require caution in interpretation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Generalized framework for likelihood-based field-level inference of growth rate from velocity and density fields

Measuring the growth rate of large-scale structures ( f ) as a function of redshift has the potential to break degeneracies between modified gravity and dark energy models, when combined with expansion-rate probes. Direct estimates of peculiar velocities of galaxies have attracted interest as a means of estimating fσ 8 . In particular, field-level methods can be used to fit the field nuisance parameter along with cosmological parameters simultaneously. This article aims to provide the community with a unified framework for the theoretical modeling of the likelihood-based field-level inference by performing fast field covariance calculations for velocity and density fields. Our purpose is to lay the foundations for a nonlinear extension of the likelihood-based method at the field level. We have developed a generalized framework, implemented in the dedicated software flip to perform a likelihood-based inference of fσ 8 . We derived a new field covariance model, which includes wide-angle corrections. We also included the models previously described in the literature inside our framework. We compared their performance against ours, and we validated our model by comparing it with the two-point statistics of a recent N-body simulation. The tests we performed have allowed us to validate our software and determine the appropriate wavenumber range to integrate our covariance model and its validity in terms of separation. Our framework allows for a wider wavenumber coverage to be used in our calculations than in previous works, which is particularly interesting for nonlinear model extensions. Finally, our generalized framework allows us to efficiently perform a survey geometry-dependent Fisher forecast of the fσ 8 parameter. We show that the Fisher forecast method we developed gives an error bar that is 30% closer to a full likelihood-based estimation than a standard volume Fisher forecast.

Ravoux, Corentin↗

DOE OSTI2021

The clustering of the SDSS-IV extended Baryon Oscillation Spectroscopic Survey quasar sample: testing observational systematics on the Baryon Acoustic Oscillation measurement

ABSTRACT Baryon Acoustic Oscillations (BAO) are considered to be a very robust standard ruler against various systematics. This premise has been tested against observational systematics, but not to the level required for the next generation of galaxy surveys such as the Dark Energy Spectroscopic Instrument (DESI) and Euclid. In this paper, we investigate the effect of observational systematics on the BAO measurement of the final sample of quasars from the extended Baryon Oscillation Spectroscopic Survey Data Release 16 in order to prepare and hone a similar analysis for upcoming surveys. We employ catalogues with various treatments of imaging systematic effects using linear and neural network-based non-linear approaches and consider how the BAO measurement changes. We also test how the variations to the BAO fitting model respond to the observational systematics. As expected, we confirm that the BAO measurements obtained from the DR16 quasar sample are robust against imaging systematics well within the statistical error, while reporting slightly modified constraints that shift the line-of-sight BAO signal by less than 1.1 per cent. We use realistic simulations with similar redshift and angular distributions as the DR16 sample to conduct statistical tests for validating the pipeline, quantifying the significance of differences, and estimating the expected bias on the BAO scale in future high-precision data sets. Although we find a marginal impact for the eBOSS QSO data, the work presented here is of vital importance for constraining the nature of dark energy with the BAO feature in the new era of big data cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Development of High-Performance Cast Crankshafts. Final Technical Report

The objective of this project was to develop technologies that would enable the production of cast crankshafts that can replace high performance forged steel crankshafts. To achieve this, the Ultimate Tensile Strength (UTS) of the new material needs to be 850 MPa with a desired minimum Yield Strength (YS; 0.2% offset) of 615 MPa and at least 10% elongation. Perhaps more challenging, the cast material needs to be able to achieve sufficient local fatigue properties to satisfy the durability requirements in today’s high performance gasoline and diesel engine applications. The project team focused on the development of cast steel alloys for application in crankshafts to take advantage of the higher stiffness over other potential material choices. The material and process developed should be able to produce high-performance crankshafts at no more than 110% of the cost of current production cast units, perhaps the most difficult objective to achieve. To minimize costs, the primary alloy design strategy was to design compositions that can achieve the required properties with minimal alloying and post-casting heat treatments. An Integrated Computational Materials Engineering (ICME) based approach was utilized, rather than relying only on traditional trial-and-error methods, which has been proven to accelerate alloy development time. Prototype melt chemistries designed using ICME were cast as test specimens and characterized iteratively to develop an alloy design within a stage-gate process. Standard characterization and material testing was done to validate the alloy performance against design targets and provide feedback to material design and manufacturing process models. Finally, the project called for Caterpillar and General Motors (GM) to develop optimized crankshaft designs using the final material and manufacturing processing path developed. A multi-disciplinary effort was to integrate finite element analyses by engine designers and geometry-specific casting simulations with existing materials models to optimize crankshaft cost and performance. Prototype crankshafts of the final design were to be produced and validated using laboratory bench testing and on-engine durability testing. ICME process simulation tools were used to investigate a broad range of processing concepts. These concepts included casting orientation, various mold and core materials, and various filling and feeding strategies. Each crankshaft was first simulated without gating and risers, which is termed natural solidification. The natural solidification results were used as a baseline for strategy development of each concept. Casting process simulations and ICME tools were proven to be reasonable predictors of real world results. Potential alloys were developed that could meet the project material property goals with appropriate normalization and temper treatments. For the alloys considered, post-normalization temper treatments proved to be necessary to achieve the desired yield strengths and elongations and appropriate heat treatments were designed using ICME tools. The experimental data of all the alloys were analyzed in combination with ICME tools to establish chemistry-process-structure relations. Several GM small gas engine (SGE) crankshafts were successfully cast in sand molds using two different sprue, runner, gate, riser, chill designs. These crankshafts were cast in two different steel alloys developed during the project, but casting finishing (e.g. riser removal) remains a cost challenge. A long list of future work was left unfinished when this project was unexpectedly terminated.

36 MATERIALS SCIENCE↗

DOE OSTI2014

SU-E-CAMPUS-I-02: Estimation of the Dosimetric Error Caused by the Voxelization of Hybrid Computational Phantoms Using Triangle Mesh-Based Monte Carlo Transport

Purpose: Computational voxel phantom provides realistic anatomy but the voxel structure may result in dosimetric error compared to real anatomy composed of perfect surface. We analyzed the dosimetric error caused from the voxel structure in hybrid computational phantoms by comparing the voxel-based doses at different resolutions with triangle mesh-based doses. Methods: We incorporated the existing adult male UF/NCI hybrid phantom in mesh format into a Monte Carlo transport code, penMesh that supports triangle meshes. We calculated energy deposition to selected organs of interest for parallel photon beams with three mono energies (0.1, 1, and 10 MeV) in antero-posterior geometry. We also calculated organ energy deposition using three voxel phantoms with different voxel resolutions (1, 5, and 10 mm) using MCNPX2.7. Results: Comparison of organ energy deposition between the two methods showed that agreement overall improved for higher voxel resolution, but for many organs the differences were small. Difference in the energy deposition for 1 MeV, for example, decreased from 11.5% to 1.7% in muscle but only from 0.6% to 0.3% in liver as voxel resolution increased from 10 mm to 1 mm. The differences were smaller at higher energies. The number of photon histories processed per second in voxels were 6.4×10{sup 4}, 3.3×10{sup 4}, and 1.3×10{sup 4}, for 10, 5, and 1 mm resolutions at 10 MeV, respectively, while meshes ran at 4.0×10{sup 4} histories/sec. Conclusion: The combination of hybrid mesh phantom and penMesh was proved to be accurate and of similar speed compared to the voxel phantom and MCNPX. The lowest voxel resolution caused a maximum dosimetric error of 12.6% at 0.1 MeV and 6.8% at 10 MeV but the error was insignificant in some organs. We will apply the tool to calculate dose to very thin layer tissues (e.g., radiosensitive layer in gastro intestines) which cannot be modeled by voxel phantoms.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

FLASSH 1.0: Thermal scattering law evaluation and cross section generation for reactor physics applications

The Full Law Analysis Scattering System Hub (FLASSH) is a modern, advanced code which evaluates the thermal scattering law (TSL) along with accompanying cross sections. FLASSH features generalized methods which accommodate any material structure. Historical approximations including the incoherent and cubic approximations have been removed. Instead, the latest release of FLASSH features advanced physics options including distinct corrections (1-phonon contributions) and non-cubic formulations. The non-cubic elastic and inelastic contributions are necessary to accurately evaluate 1-phonon contributions. Both non-cubic and 1-phonon calculations require high-density sampling of the various scattering directions. Optimization and parallelization of these routines were therefore necessary to produce results in a reasonable timeframe. With these notable improvements to the generalized TSL, FLASSH 1.0 meets benchmark requirements, demonstrating noticeable agreement with experiment for both TSLs and the resulting integrated cross sections. Additional features including a graphical user interface (GUI), plotting diagnostics, and formatted output options including ACE files allow users to complete a TSL evaluation with minimal input and maximum flexibility. The user GUI creates input files for FLASSH, reducing user error and also providing built-in error checks. Autofill options and suggested input values help make TSL evaluation accessible to novice users. The FLASSH code is compiled to run on both Windows and Linux platforms with automatic parallelization. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2016

Validation of the Modelling of Corrosion Driven Expansion in Grouted Waste-forms - 16125

One of the current preferred solutions for the disposal of Intermediate Level radioactive waste in the UK is grout encapsulation inside steel boxes for long term storage. These boxes are then to be entombed in a Geological Disposal Facility once one is available. However, both prior to and during geological disposal, waste contents can react with the water present in the encapsulant grout and begin to corrode. This corrosion product occupies a larger volume (by several hundred percent) than the original material and the corrosion occurs over the course of decades and centuries, long after the encapsulant grout has set and hardened. The expanding waste then creates stresses in the encapsulant causing both it and the surrounding box to deform and, in extreme circumstances, and after long times, for the structural integrity of the box to fail. This behaviour needs to be quantified before the box can be approved for use as any excessive deformation could lead to a loss of containment or, since the boxes are stacked on top of each other, to a significant accident in the stores. Due to the long timescales involved, physical testing of the phenomenon in full scale is impossible. Instead, a combination of chemical modelling and finite element analysis is used to predict the extent of the corrosion in the waste and its effect of the integrity of the box. There are a number of difficulties in modelling the corrosion-driven expansion and subsequent deformation of the boxes. Although the corrosion occurs over a very long time scale, the event is highly non-linear and dynamic. All of the material models, including the steel and grout are non-linear. The deformation is in the order of 10s of centimetres, well beyond the static elastic range. Finally, the components must be allowed to move relative to each other which involve a lot of part to part contacts. This means that simple static modelling, which would solve quite quickly, cannot be used. Instead a dynamic analysis code, LS-DYNA, has been used to expand the product over time to simulate the corrosion-driven expansion. This software is generally used for modelling car crashes or scenarios where waste containers have been dropped. The analyses therefore include a significant inertial component. This paper will detail the NNL's approach to modelling this long-term, quasi-static phenomenon in the dynamic, impact oriented software code LS-DYNA. It will cover the problems inherent in this approach, the errors that it creates as well as the tests and solutions that NNL worked on in order to overcome them and be confident that the solution could be validated. This work will ultimately be used to determine the maximum inventory which can be stored in each box based on the waste content and the potential extent of expansion as well as any modification that could be made to the box design or the processing streams that may improve the storage solution. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2015

Predictive Engineering Tools for Injection-Molded Long-Carbon-Fiber Thermoplastic Composites (FY 2015 Third Quarterly Report)

During the third quarter of FY 2015, the following technical progress has been made toward project milestones: 1) Magna oversaw the tool build and prepared the molding plan for the complex part of Phase II. 2) PlastiComp hosted a visit by Magna and Toyota on April 23rd to finalize the molding scope and schedule. The plan for molding trials including selection of molding parameters for both LFT and D-LFT for the U-shape complex part was established. 3) Toyota shipped the U-shape complex part tool to Magna on May 28th, 2015. 4) Plasticomp provided 30wt% LCF/PP and 30wt% LCF/PA66 compounded pellets to Magna for molding the complex part. 5) Magna performed preliminary molding trials on June 2nd, 2015 to validate wall thickness, fill profile, tool temperature and shot size requirements for the complex part. 6) Magna performed the first complex part run on June 16th and 17th, 2015 at Magna’s Composite Centre of Excellence in Concord, ON, Canada. Dale Roland of Plasticomp, and Umesh Gandhi of Toyota also attended the molding. 7) Magna discussed and finalized the plan with PNNL and the team for cutting samples from molded parts at selected locations for fiber orientation and length measurements. 8) Magna provided the computer-aided design (CAD) files of the complex parts with and without ribs to PNNL and Autodesk to build the corresponding ASMI models for injection molding simulations. Magna also provided the actual parameters used. 9) Plasticomp’s provided knowledge and experience of molding LCF materials essential to the successful molding of the parts including optimization of fill speed, tool temperatures, and plasticizing conditions for the 30wt% LCF/PP and 30wt% LCF/PA66 materials in both rib and non-rib versions. 10) Magna molded additional parts for evaluation of mechanical property testing including torsional stiffness on June 29th and 30th, 2015 at Magna’s Composite Center of Excellence. 11) Toyota began preparation for the torsion test of the specimens. Preparation of a computer-aided engineering (CAE) model to predict the performance is in progress. 12) Autodesk fixed an error in the implementation of the proper orthogonal decomposition (POD) calculation of fiber length that had caused the ASMI solution to crash and provided an updated build of ASMI containing the fix. 13)Autodesk reviewed and provided feedback for the complex part molding and measurement locations. 14) Autodesk provided support to set up the workflow for ASMI-ABAQUS® analysis, and provided a fix and workaround for a bug in the ASMI-ABAQUS® output command. 15) Autodesk helped build ASMI analysis models for the complex parts with and without ribs. 16) Autodesk worked on improving the orientation prediction accuracy in the shearing layer for 3D meshes based on comparison to measured data of the plaque moldings. 17) PNNL installed a new ASMI version received from Autodesk and performed comparative analyses to assess mid-plane versus 3D fiber length predictions using the full fiber length model and the reduced-order model (ROM) using POD. 18) PNNL presented the project scope, accomplishments, significant results and future plans to DOE and the USCAR Materials Tech Team on June 3rd, 2015. 19) PNNL discussed the cutting of samples from molded parts and finalized a plan with Magna and the team suggesting the sample size, locations and number of samples per location. 20) PNNL and Autodesk built ASMI models for the complex parts with and without ribs, and preliminary analyses of the part with ribs were conducted using the actual molding parameters received from Magna. 21) PNNL worked on a procedure to extract fiber orientation and length results from a 3D ASMI analysis to a 3D ABAQUS model. This procedure is essential to import ASMI fiber orientation and length to a 3D ABAQUS model of the part allowing future part structural analysis for weight reduction study.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Timelike Compton Scattering from a longitudinally polarised target with CLAS12 at Jefferson Lab

Explorations into the internal dynamics of hadrons are constantly evolving, and the requirement for experimental results to verify theoretical models of hadron structure is paramount. A key area in this field is the study of Generalised Parton Distributions (GPDs), which are functions used to model the momenta of quarks and gluons within hadrons, and the methods to access GPDs experimentally. One such scattering process that allows access to these is Timelike Compton Scattering. TCS complements existing Deeply Virtual Compton Scattering experiments and allows investigation into the universality of GPDs through access to the real and imaginary parts of the parton helicity independent GPD Hq via beam spin asymmetries (BSA), and it provides novel access to the real and imaginary parts of the parton helicity dependent GPD ˜Hq through target polarisation asymmetries (TSA). This thesis work presents a comparative study with the first published BSA for TCS at the Thomas Jefferson National Accelerator Facility (JLab), alongside a first time extraction of a Target Spin Asymmetry with the Summer 2022 data taking run. JLab hosts the Continuous Electron Beam Accelerator Facility (CEBAF) which provides a 12 GeV electron beam to four experimental halls. Hall-B contains the CEBAF Large Acceptance Spectrometer, which took data across three run periods on a longitudinally polarised NH3 and ND3 fixed target from 2022-2023, to extract measurements of electron-proton scattering, from which a TCS signal could be extracted. The thesis discusses work done to understand and eliminate contributions from the non-/low-polarised nuclear background, testing pre-established cuts to eliminate pion background from a dilepton (e+e-) final state and modifying them as needed for the new experimental run, and attempts to hone in on a clean TCS signal from which to extract the two asymmetry observables. A comparison with existing BSA results was performed; however, the statistical errors are too large to draw a significant conclusion as to whether there is agreement across each bin. More data is needed for a multidimensionally binned extraction. A proof of principle was achieved in the TSA measurements, with two out of four kinematic bins showing preliminary agreement in shape with theoretical values. Again, the errors are significant due to the contributions from the nuclear background. To support these conclusions, a further study was done, which takes into account an estimate of the asymmetries with the full available dataset (this thesis is based only on data taken in the summer set; at the time of writing processing was still being conducted for the final two datasets), as well as an estimate including additional future experiment days that were awarded in July 2024. Additional work was done on a secondary project exploring the feasibility of measuring TCS at the upcoming Electron Ion Collider, supporting the design proposal for the detector for the first interaction region and giving a positive outlook for the future of these types of measurements beyond JLab.

Gates, Kayleigh [Univ. of Glasgow, Scotland (Unite↗

DOE OSTI2021

On the computational solution of vector-density based continuum dislocation dynamics models: A comparison of two plastic distortion and stress update algorithms

Continuum dislocation dynamics models of mesoscale plasticity consist of dislocation transport-reaction equations coupled with crystal mechanics equations. The coupling between these two sets of equations is such that dislocation transport gives rise to the evolution of plastic distortion (strain), while the evolution of the latter fixes the stress from which the dislocation velocity field is found via a mobility law. Earlier solutions of these equations employed a staggered solution scheme for the two sets of equations in which the plastic distortion was updated via time integration of its rate, as found from Orowan's law. In this work, we show that such a direct time integration scheme can suffer from accumulation of numerical errors. We introduce an alternative scheme based on field dislocation mechanics that ensures consistency between the plastic distortion and the dislocation content in the crystal. The new scheme is based on calculating the compatible and incompatible parts of the plastic distortion separately, and the incompatible part is calculated from the current dislocation density field. Stress field and dislocation transport calculations were implemented within a finite element based discretization of the governing equations, with the crystal mechanics part solved by a conventional Galerkin method and the dislocation transport equations by the least squares method. A simple test was first performed to show the accuracy of the two schemes for updating the plastic distortion, which shows that the solution method based on field dislocation mechanics is more accurate. This method then was used to simulate an austenitic steel crystal under uniaxial loading and multiple slip conditions. By considering dislocation interactions caused by junctions, a hardening rate similar to discrete dislocation dynamics simulation results was obtained. Finally, the simulations show that dislocations exhibit some self-organized structures as the strain is increased.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Deep nonparametric estimation of intrinsic data structures by chart autoencoders: Generalization error and robustness

Autoencoders have demonstrated remarkable success in learning low-dimensional latent features of high-dimensional data across various applications. Assuming that data are sampled near a low-dimensional manifold, we employ chart autoencoders, which encode data into low-dimensional latent features on a collection of charts, preserving the topology and geometry of the data manifold. Our paper establishes statistical guarantees on the generalization error of chart autoencoders, and we demonstrate their denoising capabilities by considering n noisy training samples, along with their noise-free counterparts, on a d-dimensional manifold. By training autoencoders, we show that chart autoencoders can effectively denoise the input data with normal noise. We prove that, under proper network architectures, chart autoencoders achieve a squared generalization error in the order of n–$\frac{2}{d+2}$log 4 n, which depends on the intrinsic dimension of the manifold and only weakly depends on the ambient dimension and noise level. We further extend our theory on data with noise containing both normal and tangential components, where chart autoencoders still exhibit a denoising effect for the normal component. As a special case, our theory also applies to classical autoencoders, as long as the data manifold has a global parametrization. Furthermore, our results provide a solid theoretical foundation for the effectiveness of autoencoders, which is further validated through several numerical experiments.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

RedMaPPer: Evolution and Mass Dependence of the Conditional Luminosity Functions of Red Galaxies in Galaxy Clusters

In this work we characterize the luminosity distribution, halo-mass dependence, and redshift evolution of red galaxies in galaxy clusters using the SDSS Data Release 8 redMaPPer cluster sample. We propose a simple prescription for the relationship between the luminosity of both red central and red satellite galaxies and the mass of their host halos, and show that this model is well fit by the data. Using a larger galaxy cluster sample than previously employed in the literature, we find that the luminosities of red central galaxies scale as $\langle \mathrm{log}L\rangle \propto {A}_{L}\mathrm{log}({M}_{200b})$, with A L = 0.39 ± 0.04, and that the scatter of the red central galaxy luminosity at fixed M 200b (${\sigma }_{\mathrm{log}L| M}$) is ${0.23}_{-0.04}^{+0.05}$ dex, with the error bar including systematics due to mis-centering of the cluster finder, photometry, and photometric redshift estimation. Our data prefers a positive correlation between the luminosity of red central galaxies and the observed richness of clusters at a fixed halo mass, with an effective correlation coefficient ${d}_{\mathrm{eff}}={0.36}_{-0.16}^{+0.17}$. The characteristic luminosity of red satellites becomes dimmer from z = 0.3 to z = 0.1 by ~20% after accounting for passive evolution. We estimate the fraction of galaxy clusters where the brightest red galaxy is not the central to be P BNC ~ 20%. We discuss implications of these findings in the context of galaxy evolution and the galaxy–halo connection.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

BOSS DR12 full-shape cosmology: Λ CDM constraints from the large-scale galaxy power spectrum and bispectrum monopole

Here, we present a full Λ CDM analysis of the BOSS DR12 dataset, including information from the power spectrum multipoles, the real-space power spectrum, the reconstructed power spectrum and the bispectrum monopole. This is the first analysis to feature a complete treatment of the galaxy bispectrum, including a consistent theoretical model and without large-scale cuts. Unlike previous works, the statistics are measured using window-free estimators: this greatly reduces computational costs by removing the need to window-convolve the theory model. Our pipeline is tested using a suite of high-resolution mocks and shown to be robust and precise, with systematic errors far below the statistical thresholds. Inclusion of the bispectrum yields consistent parameter constraints and shrinks the σ 8 posterior by 13% to reach < 5 % precision; less conservative analysis choices would reduce the error bars further. Our constraints are broadly consistent with Planck: in particular, we find H 0 = 69.6$^{+1.1}_{–1.3}$ km s –1 Mpc –1 , σ 8 = 0.692$^{+0.035}_{–0.041}$ and n s = 0.870$^{+0.067}_{–0.064}$, including a BBN prior on the baryon density. When n s is set by Planck, we find H 0 = 68.31$^{+0.83}_{–0.86}$ km s –1 Mpc –1 and σ 8 = 0.722$^{+0.032}_{–0.036}$. Our S 8 posterior, 0.751 ± 0.039, is consistent with weak lensing studies, but lower than Planck. Constraints on the higher-order bias parameters are significantly strengthened from the inclusion of the bispectrum, and we find no evidence for deviation from the dark matter halo bias relations. These results represent the most complete full-shape analysis of BOSS DR12 to-date, and the corresponding spectra will enable a variety of beyond-Λ CDM analyses, probing phenomena such as the neutrino mass and primordial non-Gaussianity.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Baryon acoustic oscillation theory and modelling systematics for the DESI 2024 results

This paper provides a comprehensive overview of how fitting of baryon acoustic oscillations (BAO) is carried out within the upcoming Dark Energy Spectroscopic Instrument’s (DESI) 2024 results using its DR1 data set, and the associated systematic error budget from theory and modelling of the BAO. We derive new results showing how non-linearities in the clustering of galaxies can cause potential biases in measurements of the isotropic (⁠α iso ⁠) and anisotropic (⁠α ap ) BAO distance scales, and how these can be effectively removed with an appropriate choice of reconstruction algorithm. We then demonstrate how theory leads to a clear choice for how to model the BAO and develop, implement, and validate a new model for the remaining smooth-broad-band (i.e. without BAO) component of the galaxy clustering. Finally, we explore the impact of all remaining modelling choices on the BAO constraints from DESI using a suite of high-precision simulations, arriving at a set of best practices for DESI BAO fits, and an associated theory and modelling systematic error. Overall, our results demonstrate the remarkable robustness of the BAO to all our modelling choices and motivate a combined theory and modelling systematic error contribution to the post-reconstruction DESI BAO measurements of no more than 0.1 per cent (0.2 per cent) for its isotropic (anisotropic) distance measurements. We expect the theory and best practices laid out to here to be applicable to other BAO experiments in the era of DESI and beyond.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Description of reaction and vibrational energetics of CO 2 –NH 3 interaction using quantum computing algorithms

CO 2 capture is critical to solving global warming. Amine-based solvents are extensively used to chemically absorb CO 2 . Thus, it is crucial to study the chemical absorption of CO 2 by amine-based solvents to better understand and optimize CO 2 capture processes. Here, we use quantum computing algorithms to quantify molecular vibrational energies and reaction pathways between CO 2 and a simplified amine-based solvent model—NH 3 . Molecular vibrational properties are important to understanding kinetics of reactions. However, the molecule size correlates with the strength of anharmonicity effect on vibrational properties, which can be challenging to address using classical computing. Quantum computing can help enhance molecular vibrational calculations by including anharmonicity. We implement a variational quantum eigensolver (VQE) algorithm in a quantum simulator to calculate ground state vibrational energies of reactants and products of the CO 2 and NH 3 reaction. The VQE calculations yield ground vibrational energies of CO 2 and NH 3 with similar accuracy to classical computing. In the presence of hardware noise, Compact Heuristic for Chemistry (CHC) ansatz with shallower circuit depth performs better than Unitary Vibrational Coupled Cluster. The “Zero Noise Extrapolation” error-mitigation approach in combination with CHC ansatz improves the vibrational calculation accuracy. Excited vibrational states are accessed with quantum equation of motion method for CO 2 and NH 3 . Using quantum Hartree–Fock (HF) embedding algorithm to calculate electronic energies, the corresponding reaction profile compares favorably with Coupled Cluster Singles and Doubles while being more accurate than HF. Our research showcases quantum computing applications in the study of CO 2 capture reactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2013

Spatial frequency spectrum of the x-ray scatter distribution in CBCT projections

Purpose: X-ray scatter is a source of significant image quality loss in cone-beam computed tomography (CBCT). The use of Monte Carlo (MC) simulations separating primary and scattered photons has allowed the structure and nature of the scatter distribution in CBCT to become better elucidated. This work seeks to quantify the structure and determine a suitable basis function for the scatter distribution by examining its spectral components using Fourier analysis.Methods: The scatter distribution projection data were simulated using a CBCT MC model based on the EGSnrc code. CBCT projection data, with separated primary and scatter signal, were generated for a 30.6 cm diameter water cylinder [single angle projection with varying axis-to-detector distance (ADD) and bowtie filters] and two anthropomorphic phantoms (head and pelvis, 360 projections sampled every 1°, with and without a compensator). The Fourier transform of the resulting scatter distributions was computed and analyzed both qualitatively and quantitatively. A novel metric called the scatter frequency width (SFW) is introduced to determine the scatter distribution's frequency content. The frequency content results are used to determine a set basis functions, consisting of low-frequency sine and cosine functions, to fit and denoise the scatter distribution generated from MC simulations using a reduced number of photons and projections. The signal recovery is implemented using Fourier filtering (low-pass Butterworth filter) and interpolation. Estimates of the scatter distribution are used to correct and reconstruct simulated projections.Results: The spatial and angular frequencies are contained within a maximum frequency of 0.1 cm{sup −1} and 7/(2π) rad{sup −1} for the imaging scenarios examined, with these values varying depending on the object and imaging setup (e.g., ADD and compensator). These data indicate spatial and angular sampling every 5 cm and π/7 rad (∼25°) can be used to properly capture the scatter distribution, with reduced sampling possible depending on the imaging scenario. Using a low-pass Butterworth filter, tuned with the SFW values, to denoise the scatter projection data generated from MC simulations using 10{sup 6} photons resulted in an error reduction of greater than 85% for the estimating scatter in single and multiple projections. Analysis showed that the use of a compensator helped reduce the error in estimating the scatter distribution from limited photon simulations by more than 37% when compared to the case without a compensator for the head and pelvis phantoms. Reconstructions of simulated head phantom projections corrected by the filtered and interpolated scatter estimates showed improvements in overall image quality.Conclusions: The spatial frequency content of the scatter distribution in CBCT is found to be contained within the low frequency domain. The frequency content is modulated both by object and imaging parameters (ADD and compensator). The low-frequency nature of the scatter distribution allows for a limited set of sine and cosine basis functions to be used to accurately represent the scatter signal in the presence of noise and reduced data sampling decreasing MC based scatter estimation time. Compensator induced modulation of the scatter distribution reduces the frequency content and improves the fitting results.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2014

Neutrino constraints: what large-scale structure and CMB data are telling us?

We discuss the reliability of neutrino mass constraints, either active or sterile, from the combination of different low redshift Universe probes with measurements of CMB anisotropies. In our analyses we consider WMAP 9-year or Planck Cosmic Microwave Background (CMB) data in combination with Baryonic Acoustic Oscillations (BAO) measurements from BOSS DR11, galaxy shear measurements from CFHTLenS, SDSS Ly α forest constraints and galaxy cluster mass function from Chandra observations. At odds with recent similar studies, to avoid model dependence of the constraints we perform a full likelihood analysis for all the datasets employed. As for the cluster data analysis we rely on to the most recent calibration of massive neutrino effects in the halo mass function and we explore the impact of the uncertainty in the mass bias and re-calibration of the halo mass function due to baryonic feedback processes on cosmological parameters. We find that none of the low redshift probes alone provide evidence for massive neutrino in combination with CMB measurements, while a larger than 2σ detection of non zero neutrino mass, either active or sterile, is achieved combining cluster or shear data with CMB and BAO measurements. Yet, the significance of the detection exceeds 3σ if we combine all four datasets. For a three active neutrino scenario, from the joint analysis of CMB, BAO, shear and cluster data including the uncertainty in the mass bias we obtain ∑ m{sub ν} =0.29{sup +0.18}{sub -0.21} eV and ∑ m{sub ν} =0.22{sup +0.17}{sub -0.18} eV 95%CL) using WMAP9 or Planck as CMB dataset, respectively. The preference for massive neutrino is even larger in the sterile neutrino scenario, for which we get m{sub s}{sup eff}=0.44{sup +0.28}{sub -0.26} eV and Δ N{sub eff}=0.78{sup +0.60}{sub -0.59} 95%CL) from the joint analysis of Planck, BAO, shear and cluster datasets. For this data combination the vanilla ΛCDM model is rejected at more than 3σ and a sterile neutrino mass as motivated by accelerator anomaly is within the 2σ errors. Conversely, the Ly α data favour vanishing neutrino masses and from the data combination Planck+BAO+Ly α we get the tight upper limits ∑ m{sub ν} <0.14 eV and m{sub s}{sup eff}<0.22 eV—Δ N{sub eff}<1.11 95%CL) for the active and sterile neutrino model, respectively. Finally, results from the full data combination reflect the tension between the σ{sub 8} constraints obtained from cluster and shear data and that inferred from Ly α forest measurements; in the active neutrino scenario for both CMB datasets employed, the full data combination yields only an upper limits on ∑ m{sub ν}, while assuming an extra sterile neutrino we still get preference for non-vanishing mass, m{sub s}{sup eff}=0.26{sup +0.22}{sub -0.24} eV, and dark contribution to the radiation content, Δ N{sub eff}=0.82±0.55.

79 ASTRONOMY AND ASTROPHYSICS↗