Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HISTOGRAM”

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

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

At least 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

The spatial structure of young stellar clusters. II. Total young stellar populations

We investigate the intrinsic stellar populations (estimated total numbers of OB and pre-main-sequence stars down to 0.1 M{sub ⊙}) that are present in 17 massive star-forming regions (MSFRs) surveyed by the MYStIX project. The study is based on the catalog of >31,000 MYStIX Probable Complex Members with both disk-bearing and disk-free populations, compensating for extinction, nebulosity, and crowding effects. Correction for observational sensitivities is made using the X-ray luminosity function and the near-infrared initial mass function—a correction that is often not made by infrared surveys of young stars. The resulting maps of the projected structure of the young stellar populations, in units of intrinsic stellar surface density, allow direct comparison between different regions. Several regions have multiple dense clumps, similar in size and density to the Orion Nebula Cluster. The highest projected density of ∼34,000 stars pc{sup −2} is found in the core of the RCW 38 cluster. Histograms of surface density show different ranges of values in different regions, supporting the conclusion of Bressert et al. that no universal surface-density threshold can distinguish between clustered and distributed star formation. However, a large component of the young stellar population of MSFRs resides in dense environments of 200–10,000 stars pc{sup −2} (including within the nearby Orion molecular clouds), and we find that there is no evidence for the B10 conclusion that such dense regions form an extreme “tail” of the distribution. Tables of intrinsic populations for these regions are used in our companion study of young cluster properties and evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

The statistics of BAT-to-XRT flux ratio in GRBs: Evidence for a characteristic value and its implications

We present the statistics of the ratio, R, between the prompt and afterglow “plateau” fluxes of gamma-ray bursts (GRBs). We define this as the ratio of the mean prompt energy flux in Swift BAT and the Swift XRT one, immediately following the steep transition between these two states and the beginning of the afterglow stage referred to as the “plateau”. Like the distribution of many other GRB observables, the histogram of R is log-normal with maximum at a value R{sub m}≃2000, FWHM of about two decades, and with the entire distribution spanning about five decades in the value of R. We note that the peak of the distribution is close to the proton-to-electron mass ratio (R{sub m}≃m{sub p}/m{sub e}=1836), as proposed to be the case in an earlier publication, on the basis of a specific model of the GRB dissipation process. It therefore appears that, in addition to the values of the energy of peak luminosity E{sub pk}∼m{sub e}c{sup 2}, GRBs present us with one more quantity with an apparent characteristic value. The fact that the values of both these quantities (E{sub pk} and R) are consistent with the same specific model invoked to account for the efficient conversion of their relativistic proton energies to electrons argues favorably for its underlying assumptions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

An atom probe perspective on phase separation and precipitation in duplex stainless steels

Here, three-dimensional chemical imaging of Fe–Cr alloys showing Fe-rich (α)/Cr-rich (α') phase separation is reported using atom probe tomography techniques. The extent of phase separation, i.e., amplitude and wavelength, has been quantitatively assessed using the Langer-Bar-on-Miller, proximity histogram, and autocorrelation function methods for two separate Fe–Cr alloys, designated 2101 and 2205. Although the 2101 alloy possesses a larger wavelength and amplitude after annealing at 427 °C for 100–10 000 h, it exhibits a lower hardness than the 2205 alloy. In addition to this phase separation, ultra-fine Ni–Mn–Si–Cu-rich G-phase precipitates form at the α/α' interfaces in both alloys. For the 2101 alloy, Cu clusters act to form a nucleus, around which a Ni–Mn–Si shell develops during the precipitation process. For the 2205 alloy, the Ni and Cu atoms enrich simultaneously and no core–shell chemical distribution was found. This segregation phenomenon may arise from the exact Ni/Cu ratio inside the ferrite. After annealing for 10 000 h, the number density of the G-phase within the 2205 alloy was found to be roughly one order of magnitude higher than in the 2101 alloy. The G-phase precipitates have an additional deleterious effect on the thermal embrittlement, as evaluated by the Ashby–Orowan equation, which explains the discrepancy between the hardness and the rate of phase separation with respect to annealing time (Gladman T 1999 Mater. Sci. Tech. Ser. 15 30–36).

36 MATERIALS SCIENCE↗
DOE OSTI2025

CaloChallenge 2022: a community challenge for fast calorimeter simulation

Here, we present the results of the ‘Fast Calorimeter Simulation Challenge 2022’—the CaloChallenge. We study state-of-the-art generative models on four calorimeter shower datasets of increasing dimensionality, ranging from a few hundred voxels to a few tens of thousand voxels. The 31 individual submissions span a wide range of current popular generative architectures, including variational autoencoders (VAEs), generative adversarial networks (GANs), normalizing flows, diffusion models, and models based on conditional flow matching. We compare all submissions in terms of quality of generated calorimeter showers, as well as shower generation time and model size. To assess the quality we use a broad range of different metrics including differences in one-dimensional histograms of observables, KPD/FPD scores, AUCs of binary classifiers, and the log-posterior of a multiclass classifier. The results of the CaloChallenge provide the most complete and comprehensive survey of cutting-edge approaches to calorimeter fast simulation to date. In addition, our work provides a uniquely detailed perspective on the important problem of how to evaluate generative models. As such, the results presented here should be applicable for other domains that use generative AI and require fast and faithful generation of samples in a large phase space.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

Estimating the density of states of frustrated spin systems

Abstract Estimating the density of states (DOS) of systems with rugged free energy landscapes is a notoriously difficult task of the utmost importance in many areas of physics ranging from spin glasses to biopolymers. DOS estimation has also recently become an indispensable tool for the benchmarking of quantum annealers when these function as samplers. Some of the standard approaches suffer from a spurious convergence of the estimates to metastable minima, and these cases are particularly hard to detect. Here, we introduce a sampling technique based on population annealing enhanced with a multi-histogram analysis and report on its performance for spin glasses. We demonstrate its ability to overcome the pitfalls of other entropic samplers, resulting in some cases in large scaling advantages that can lead to the uncovering of new physics. The new technique avoids some inherent difficulties in established approaches and can be applied to a wide range of systems without relevant tailoring requirements. Benchmarking of the studied techniques is facilitated by the introduction of several schemes that allow us to achieve exact counts of the degeneracies of the tested instances.

Barash, Lev (ORCID:000000022298785X)↗
DOE OSTI2022

Organizing Large Data Sets for Efficient Analyses on HPC Systems

Upcoming exascale applications could introduce significant data management challenges due to their large sizes, dynamic work distribution, and involvement of accelerators such as graphical processing units, GPUs. In this work, we explore the performance of reading and writing operations involving one such scientific application on two different supercomputers. Our tests showed that the Adaptable Input and Output System, ADIOS, was able to achieve speeds over 1TB/s, a significant fraction of the peak I/O performance on Summit. We also demonstrated the querying functionality in ADIOS could effectively support common selective data analysis operations, such as conditional histograms. In tests, this query mechanism was able to reduce the execution time by a factor of five. More importantly, ADIOS data management framework allows us to achieve these performance improvements with only a minimal amount of coding effort.

Gu, Junmin↗
DOE OSTI2014

Development of Bayesian analysis program for extraction of polarisation observables at CLAS

At the mass scale of a proton, the strong force is not well understood. Various quark models exist, but it is important to determine which quark model(s) are most accurate. Experimentally, finding resonances predicted by some models and not others would give valuable insight into this fundamental interaction. Several labs around the world use photoproduction experiments to find these missing resonances. The aim of this work is to develop a robust Bayesian data analysis program for extracting polarisation observables from pseudoscalar meson photoproduction experiments using CLAS at Jefferson Lab. This method, known as nested sampling, has been compared to traditional methods and has incorporated data parallelisation and GPU programming. It involves an event-by-event likelihood function, which has no associated loss of information from histogram binning, and results can be easily constrained to the physical region. One of the most important advantages of the nested sampling approach is that data from different experiments can be combined and analysed simultaneously. Results on both simulated and previously analysed experimental data for the K{sup +}Λ channel will be discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Open access to high-level data and analysis tools in the CMS experiment at the LHC

The CMS experiment, in recognition of its commitment to data preservation and open access as well as to education and outreach, has made its first public release of high-level data under the CC0 waiver: up to half of the proton-proton collision data (by volume) at 7 TeV from 2010 in CMS Analysis Object Data format. CMS has prepared, in collaboration with CERN and the other LHC experiments, an open-data web portal based on Invenio. The portal provides access to CMS public data as well as to analysis tools and documentation for the public. The tools include an event display and histogram application that run in the browser. In addition a virtual machine containing a CMS software environment along with XRootD access to the data is available. Within the virtual machine the public can analyse CMS data, example code is provided. As a result, we describe the accompanying tools and documentation and discuss the first experiences of data use.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2025

Tools for unbinned unfolding

Machine learning has enabled differential cross section measurements that are not discretized. Going beyond the traditional histogram-based paradigm, these unbinned unfolding methods are rapidly being integrated into experimental workflows. Here, in order to enable widespread adaptation and standardization, we develop methods, benchmarks, and software for unbinned unfolding. For methodology, we demonstrate the utility of boosted decision trees for unfolding with a relatively small number of high-level features. This complements state-of-the-art deep learning models capable of unfolding the full phase space. To benchmark unbinned unfolding methods, we develop an extension of existing dataset to include acceptance effects, a necessary challenge for real measurements. Additionally, we directly compare binned and unbinned methods using discretized inputs for the latter in order to control for the binning itself. Lastly, we have assembled two software packages for the OmniFold unbinned unfolding method that should serve as the starting point for any future analyses using this technique. One package is based on the widely-used RooUnfold framework and the other is a standalone package available through the Python Package Index (PyPI).

47 OTHER INSTRUMENTATION↗
DOE OSTI2014

IDENTIFYING POTENTIAL MARKERS OF THE SUN'S GIANT CONVECTIVE SCALE

Line-of-sight magnetograms from the Helioseismic and Magnetic Imager (HMI) of the Solar Dynamics Observatory (SDO) are analyzed using a diagnostic known as the magnetic range of influence (MRoI). The MRoI is a measure of the length over which a photospheric magnetogram is balanced and so its application gives the user a sense of the connective length scales in the outer solar atmosphere. The MRoI maps and histograms inferred from the SDO/HMI magnetograms primarily exhibit four scales: a scale of a few megameters that can be associated with granulation, a scale of a few tens of megameters that can be associated with super-granulation, a scale of many hundreds to thousands of megameters that can be associated with coronal holes and active regions, and a hitherto unnoticed scale that ranges from 100 to 250 Mm. We infer that this final scale is an imprint of the (rotationally driven) giant convective scale on photospheric magnetism. This scale appears in MRoI maps as well-defined, spatially distributed concentrations that we have dubbed ''g-nodes''. Furthermore, using coronal observations from the Atmospheric Imaging Assembly on SDO, we see that the vicinity of these g-nodes appears to be a preferred location for the formation of extreme-ultraviolet (and likely X-Ray) brightpoints. These observations and straightforward diagnostics offer the potential of a near real-time mapping of the Sun's largest convective scale, a scale that possibly reaches to the very bottom of the convective zone.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Molecular counting by photobleaching in protein complexes with many subunits: best practices and application to the cellulose synthesis complex

The constituents of large, multisubunit protein complexes dictate their functions in cells, but determining their precise molecular makeup in vivo is challenging. One example of such a complex is the cellulose synthesis complex (CSC), which in plants synthesizes cellulose, the most abundant biopolymer on Earth. In growing plant cells, CSCs exist in the plasma membrane as six-lobed rosettes that contain at least three different cellulose synthase (CESA) isoforms, but the number and stoichiometry of CESAs in each CSC are unknown. To begin to address this question, we performed quantitative photobleaching of GFP-tagged AtCESA3-containing particles in living Arabidopsis thaliana cells using variable-angle epifluorescence microscopy and developed a set of information-based step detection procedures to estimate the number of GFP molecules in each particle. The step detection algorithms account for changes in signal variance due to changing numbers of fluorophores, and the subsequent analysis avoids common problems associated with fitting multiple Gaussian functions to binned histogram data. The analysis indicates that at least 10 GFP-AtCESA3 molecules can exist in each particle. In conclusion, these procedures can be applied to photobleaching data for any protein complex with large numbers of fluorescently tagged subunits, providing a new analytical tool with which to probe complex composition and stoichiometry.

09 BIOMASS FUELS↗
DOE OSTI2020

Global adjoint tomography—model GLAD-M25

SUMMARY Building on global adjoint tomography model GLAD-M15, we present transversely isotropic global model GLAD-M25, which is the result of 10 quasi-Newton tomographic iterations with an earthquake database consisting of 1480 events in the magnitude range 5.5 ≤ Mw ≤ 7.2, an almost sixfold increase over the first-generation model. We calculated fully 3-D synthetic seismograms with a shortest period of 17 s based on a GPU-accelerated spectral-element wave propagation solver which accommodates effects due to 3-D anelastic crust and mantle structure, topography and bathymetry, the ocean load, ellipticity, rotation and self-gravitation. We used an adjoint-state method to calculate Fréchet derivatives in 3-D anelastic Earth models facilitated by a parsimonious storage algorithm. The simulations were performed on the Cray XK7 ‘Titan’ and the IBM Power 9 ‘Summit’ at the Oak Ridge Leadership Computing Facility. We quantitatively evaluated GLAD-M25 by assessing misfit reductions and traveltime anomaly histograms in 12 measurement categories. We performed similar assessments for a held-out data set consisting of 360 earthquakes, with results comparable to the actual inversion. We highlight the new model for a variety of plumes and subduction zones.

58 GEOSCIENCES↗
DOE OSTI2021

encore : an O ( N g2) estimator for galaxy N -point correlation functions

ABSTRACT We present a new algorithm for efficiently computing the N-point correlation functions (NPCFs) of a 3D density field for arbitrary N. This can be applied both to a discrete spectroscopic galaxy survey and a continuous field. By expanding the statistics in a separable basis of isotropic functions built from spherical harmonics, the NPCFs can be estimated by counting pairs of particles in space, leading to an algorithm with complexity $\mathcal {O}(N_\mathrm{g}^2)$ for Ng particles, or $\mathcal {O}(N_\mathrm{FFT}\log N_\mathrm{FFT})$ when using a Fast Fourier Transform with NFFT grid-points. In practice, the rate-limiting step for N > 3 will often be the summation of the histogrammed spherical harmonic coefficients, particularly if the number of radial and angular bins is large. In this case, the algorithm scales linearly with Ng. The approach is implemented in the encore code, which can compute the 3PCF, 4PCF, 5PCF, and 6PCF of a BOSS-like galaxy survey in ${\sim}100$ CPU-hours, including the corrections necessary for non-uniform survey geometries. We discuss the implementation in depth, along with its GPU acceleration, and provide practical demonstration on realistic galaxy catalogues. Our approach can be straightforwardly applied to current and future data sets to unlock the potential of constraining cosmology from the higher point functions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Scaling pair count to next galaxy surveys

ABSTRACT Counting pairs of galaxies or stars according to their distance is at the core of real-space correlation analyses performed in astrophysics and cosmology. Upcoming galaxy surveys (LSST, Euclid) will measure properties of billions of galaxies challenging our ability to perform such counting in a minute-scale time relevant for the usage of simulations. The problem is only limited by efficient access to the data, hence belongs to the big data category. We use the popular Apache Spark framework to address it and design an efficient high-throughput algorithm to deal with hundreds of millions to billions of input data. To optimize it, we revisit the question of non-hierarchical sphere pixelization based on cube symmetries and develop a new one dubbed the ‘Similar Radius Sphere Pixelization’ (SARSPix) with very close to square pixels. It provides the most adapted indexing over the sphere for all distance-related computations. Using LSST-like fast simulations, we compute autocorrelation functions on tomographic bins containing between a hundred million to one billion data points. In each case, we achieve the construction of a standard pair-distance histogram in about 2 min, using a simple algorithm that is shown to scale, over a moderate number of nodes (16–64). This illustrates the potential of this new techniques in the field of astronomy where data access is becoming the main bottleneck. They can be easily adapted to other use-cases as nearest-neighbours search, catalogue cross-match or cluster finding. The software is publicly available from https://github.com/astrolabsoftware/SparkCorr.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The role of magnetic fields in the stability and fragmentation of filamentary molecular clouds: two case studies at OMC-3 and OMC-4

We present the stability analysis of two regions, OMC-3 and OMC-4, in the massive and long molecular cloud complex of Orion A. We obtained 214 $\mu$m HAWC + /SOFIA polarization data, and we make use of archival data for the column density and C 18 O (1–0) emission line. We find clear depolarization in both observed regions and that the polarization fraction is anticorrelated with the column density and the polarization-angle dispersion function. We find that the filamentary cloud and dense clumps in OMC-3 are magnetically supercritical and strongly subvirial. This region should be in the gravitational collapse phase and is consistent with many young stellar objects (YSOs) forming in the region. Our histogram of relative orientation (HRO) analysis shows that the magnetic field is dynamically sub-dominant in the dense gas structures of OMC-3. We present the first polarization map of OMC-4. We find that the observed region is generally magnetically subcritical except for an elongated dense core, which could be a result of projection effect of a filamentary structure aligned close to the line of sight. The relative large velocity dispersion and the unusual positive shape parameters at high column densities in the HROs analysis suggest that our viewing angle may be close to axes of filamentary substructures in OMC-4. The dominating strong magnetic field in OMC-4 is unfavourable for star formation and is consistent with much fewer YSOs than in OMC-3.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Relative alignment between magnetic fields and molecular gas structure in molecular clouds

Here, we compare the structure of synthetic dust polarization with synthetic molecular line emission from radiative transfer calculations using a three-dimensional, turbulent collapsing-cloud magnetohydrodynamics simulation. The histogram of relative orientation (HRO) technique and the projected Rayleigh statistic (PRS) are considered. In our trans-Alfvénic (more strongly magnetized) simulation, there is a transition to perpendicular alignment at densities above ~4 × 10 3 cm –3 . This transition is recovered in most of our synthetic observations of optically thin molecular tracers; however, for 12 CO it does not occur and the PRS remains in parallel alignment across the whole observer space. We calculate the physical depth of the optical depth τ = 1 surface and find that for 12 CO it is largely located in front of the cloud mid-plane, suggesting that 12 CO is too optically thick and instead mainly probes low-volume density gas. In our super-Alfvénic simulation, the magnetic field becomes significantly more tangled, and all observed tracers tend towards no preference for perpendicular or parallel alignment. An observable difference in alignment between optically thin and optically thick tracers may indicate the presence of a dynamically important magnetic field, though there is some degeneracy with viewing angle. We convolve our data with a Gaussian beam and compare it with HRO results of the Vela C molecular cloud. We find good agreement between these results and our sub-Alfvénic simulations when viewed with the magnetic field in the plane of the sky (especially when sensitivity limitations are considered), though the observations are also consistent with an intermediately inclined magnetic field.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Multiparticle cumulant mapping for Coulomb explosion imaging: Calculations and algorithm

We present a versatile cumulant mapping algorithm for analyzing correlated particle emission, offering insights into complex electronic and nuclear dynamics. Recently, we have demonstrated the use of cumulant mapping to extract information-rich correlations between the momenta of multiple fragments produced in Coulomb explosion imaging experiments [C. Cheng et al., Phys. Rev. Lett. 130, 093001 (2023)]. We define cumulant mapping in terms of histograms, enabling fast computation of linear (additive) observables. However, applying the same algorithm to nonlinear (nonadditive) observables poses challenges, as the computation time of conventional estimators scales nonlinearly with data size. To overcome this, we develop estimators and an accompanying algorithm to enable computationally efficient estimation of the cumulant of interest. Comparisons of computation times and signal-to-noise ratios reveal the superior performance of our approach. This method is demonstrated on the (D+, D+, C+, O+) dissociation channel of CD 2 ⁢O 4+ produced in a strong-field ionization experiment. Additionally, Poisson statistics are used to simulate the two methods and provide insights into the efficiency of our algorithm. The proposed methodology unlocks efficient computation of cumulant mapping for a broader range of complex systems and observables, such as the laser pulse dependence of ionization dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Orientational ordering, buckling, and dynamic transitions for vortices interacting with a periodic quasi-one-dimensional substrate

In this study, we examine the statics and dynamics of vortices in the presence of a periodic quasi-one-dimensional substrate, focusing on the limit where the vortex lattice constant is smaller than the substrate lattice period. As a function of the substrate strength and filling factor, within the pinned state, we observe a series of order-disorder transitions associated with buckling phenomena in which the number of vortex rows that fit between neighboring substrate maxima increases. These transitions coincide with steps in the depinning threshold, jumps in the density of topological defects, and changes in the structure factor. At the buckling transition, the vortices are disordered, while between the buckling transitions the vortices form a variety of crystalline and partially ordered states. In the weak substrate limit, the buckling transitions are absent and the vortices form an ordered hexagonal lattice that undergoes changes in its orientation with respect to the substrate as a function of vortex density. At intermediate substrate strength, certain ordered states appear that are correlated with peaks in the depinning force. Under an applied drive, the system exhibits a rich variety of distinct dynamical phases, including plastic flow, a density-modulated moving crystal, and moving floating solid phases. We also find a dynamic smectic-to-smectic transition in which the smectic ordering changes from being aligned with the substrate to being aligned with the external drive. The different dynamical phases can be characterized using velocity histograms and the structure factor. We discuss how these results are related to recent experiments on vortex ordering in thin films with periodic thickness modulations. In conclusion, our results should also be relevant for other types of systems such as ions, colloids, or Wigner crystals interacting with periodic quasi-one-dimensional substrates.

36 MATERIALS SCIENCE↗