SEARCH · Search NASA

Results for “parallel data structures”

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 307 records · Page 17

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Enhancing Terminal Deoxynucleotidyl Transferase Activity on Substrates with 3′ Terminal Structures for Enzymatic De Novo DNA Synthesis

Enzymatic oligonucleotide synthesis methods based on the template-independent polymerase terminal deoxynucleotidyl transferase (TdT) promise to enable the de novo synthesis of long oligonucleotides under mild, aqueous conditions. Intermediates with a 3′ terminal structure (hairpins) will inevitably arise during synthesis, but TdT has poor activity on these structured substrates, limiting its usefulness for oligonucleotide synthesis. Here, we described two parallel efforts to improve the activity of TdT on hairpins: (1) optimization of the concentrations of the divalent cation cofactors and (2) engineering TdT for enhanced thermostability, enabling reactions at elevated temperatures. By combining both of these improvements, we obtained a ~10-fold increase in the elongation rate of a guanine-cytosine hairpin.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2018

The finale of a trilogy: comparing terpolymers and ternary blends with structurally similar backbones for use in organic bulk heterojunction solar cells

Building on our previous works that compared the efficacy of terpolymers vs. ternary blends in improving the performance of bulk heterojunction organic solar cells, the final piece of this series of studies focuses on comparing terpolymer and ternary blends constructed with two polymers with structurally similar backbones (monoCNTAZ and FTAZ) yet markedly different open circuit voltage (V oc ) values. Furthermore, terpolymers and ternary blends of five different ratios were studied and the results demonstrate that while the overall performance of both the systems is similar, the ternary blends exhibit higher short circuit current (J sc ) values, while the terpolymers exhibit higher V oc values. Investigation of the charge transfer state using low-energy external quantum efficiency (EQE) indicates that the ternary blends are governed by a parallel-like mechanism, while the terpolymer does not follow this mechanism. Furthermore, the key morphological difference between the systems, as elucidated by resonance soft X-ray scattering (RSoXS), is the slightly smaller size (~60 nm) of domains in the ternary blends compared to that of the terpolymer (~80 nm), which may affect exciton harvesting in the terpolymer system and lead to lower J sc values. In addition, a lower driving force for the formation of charge transfer (CT) state is also likely to contribute to the lower J sc values in the terpolymer system. All together, the data show that structurally similar (perhaps even miscible) polymers still exhibit key differences in performance when paired in terpolymers vs. ternary blends and allow us to further illuminate the underlying mechanisms of such complex systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Development of Steady-State and Dynamic Mass and Energy Constrained Neural Networks for Distributed Chemical Systems Using Noisy Transient Data

The paper presents the development of algorithms for mass and energy constrained neural network models that can exactly conserve the overall mass and energy of distributed chemical process systems, even though the noisy transient data used for optimal model training violate the same. In contrast to approximately satisfying mass and energy balance constraints of a system by soft penalization of objective function, algorithms have been developed for solving equality-constrained nonlinear optimization problems, thus providing the guarantee of exactly satisfying the system mass and energy conservation laws. For developing dynamic mass-energy constrained network models for distributed systems, hybrid series and parallel dynamic-static neural networks have been leveraged. The developed algorithms for solving both the training and forward problems are validated using both steady-state and dynamic data in the presence of various noise characteristics. The developed data-driven algorithms are flexible to exactly satisfy mass and energy balance constraints for dynamic chemical processes if the system holdup information is available. The proposed network structures and algorithms are applied to the development of data-driven lumped and distributed models of an adiabatic superheater/reheater system, a nonisothermal continuous stirred tank reactor, as well as an electrically heated plug-flow reactor system where one form of energy gets transformed to another. It has been observed that the mass-energy constrained neural networks yield a root mean squared error of <1% with respect to the system truth for the case studies evaluated in this work.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Parthenon—a performance portable block-structured adaptive mesh refinement framework

On the path to exascale the landscape of computer device architectures and corresponding programming models has become much more diverse. While various low-level performance portable programming models are available, support at the application level lacks behind. To address this issue, we present the performance portable block-structured adaptive mesh refinement (AMR) framework Parthenon, derived from the well-tested and widely used Athena++ astrophysical magnetohydrodynamics code, but generalized to serve as the foundation for a variety of downstream multi-physics codes. Parthenon adopts the Kokkos programming model, and provides various levels of abstractions from multidimensional variables, to packages defining and separating components, to launching of parallel compute kernels. Parthenon allocates all data in device memory to reduce data movement, supports the logical packing of variables and mesh blocks to reduce kernel launch overhead, and employs one-sided, asynchronous MPI calls to reduce communication overhead in multi-node simulations. Using a hydrodynamics miniapp, we demonstrate weak and strong scaling on various architectures including AMD and NVIDIA GPUs, Intel and AMD x86 CPUs, IBM Power9 CPUs, as well as Fujitsu A64FX CPUs. At the largest scale on Frontier (the first TOP500 exascale machine), the miniapp reaches a total of 1.7 × 10 13 zone-cycles/s on 9216 nodes (73,728 logical GPUs) at [Formula: see text] weak scaling parallel efficiency (starting from a single node). In combination with being an open, collaborative project, this makes Parthenon an ideal framework to target exascale simulations in which the downstream developers can focus on their specific application rather than on the complexity of handling massively-parallel, device-accelerated AMR.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Conceptual Designs for Irradiation Creep Testing of SiC in HFIR

Understanding irradiation creep of nuclear fuel cladding is important to properly size the initial fuel-cladding gap and understand when pellet-cladding contact is expected to occur due to a combination of fuel swelling and cladding creep-down. Irradiation creep also plays a role in relaxing stresses that develop in-pile. Silicon carbide fiber–reinforced silicon carbide matrix (SiC/SiC) composites are the leading long-term accident-tolerant fuel cladding concept for light-water reactors (LWRs). Although some limited data are available regarding irradiation creep of the individual constituents (fibers, matrix), data regarding irradiation creep of SiC/SiC composites are currently insufficient. Additional data regarding irradiation creep compliance and the rupture lifetime (combination of creep and slow crack growth) are needed to understand material limitations. This work describes the design and development of two irradiation vehicles that are being pursued for testing SiC/SiC concepts in the High Flux Isotope Reactor (HFIR). The first is a passive experiment, referred to as the PRECISE experiment, that leverages the constant coolant pressure of HFIR to compress a metallic bellows and provide a well-characterized load to drive creep in a SiC/SiC dog bone specimen. The total creep strain would be quantified post-irradiation by measuring dimensional changes of the specimen length as well as local dimensional changes within the gauge region. Non-stressed specimens would also be irradiated under the same conditions to provide an indication of dimensional changes due to radiation-induced swelling in the absence of creep. A second, more complex experiment, referred to as the INSITE experiment, is being designed in parallel that would use pneumatics to pressurize a metal bellows and linear variable differential transformers (LVDTs) to measure the specimen displacement in situ during irradiation. Such an experiment would provide significantly more data regarding the evolution of the creep compliance as a function of dose and applied stress within a single experiment but would require significantly more development time and cost to execute. The primary concern with the INSITE experiment is the accuracy, reliability, and expected lifetime of the LVDTs during irradiation at elevated temperatures. Efforts are being made to adjust the experiment design and operating procedure to limit LVDT temperatures and mitigate or otherwise compensate for uncertainties due to factors such as temperature fluctuations, creep in the surrounding structural materials, and drift of the LVDTs. This work describes the experiment designs, thermal and structural analysis that were performed to ensure that the desired temperature and stress conditions can be achieved, some initial sensitivity analyses to predict the evolution of the radiation-induced specimen displacements, and potential sources of uncertainty in the measurements. Out-of-pile testing is being performed in parallel to confirm that the test trains achieve the expected stress states in the specimens and do not result in prohibitive stress concentrators (e.g., in the grip regions) that might risk pre-mature failure. The PRECISE experiments are proceeding toward fabrication and assembly with HFIR insertion planned during fiscal year 2026. The INSITE experiment is progressing toward out-of-pile demonstrations, which will provide more conclusive evidence regarding the feasibility of executing these tests in HFIR or whether alternative displacement monitoring techniques may need to be considered.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2017

An alternate mode of oligomerization for E. coli SecA

SecA is the ATPase of preprotein translocase. SecA is a dimer in solution and changes in its oligomeric state may function in preprotein translocation. The SecA-N68 construct, in which the C-terminal helical domains of SecA are deleted, was used to investigate the mechanism of SecA oligomerization. SecA-N68 is in equilibrium between monomers, dimers, and tetramers. Subunit interactions in the SecA-N68 tetramer are mediated entirely by unstructured regions at its N- and C-termini: when the termini are deleted to yield SecA-N68ΔNC, the construct is completely monomeric. This monomeric construct yielded crystals diffracting to 2.6 Å that were used to solve the structure of SecA-N68, including the “preprotein crosslinking domain” (PPXD) that was missing from previous E. coli SecA structures. The SecA-N68 structure was combined with small angle X-ray scattering (SAXS) data to construct a model of the SecA-N68 tetramer that is consistent with the essential roles of the extreme N- and C-termini in oligomerization. This mode of oligomerization, which depends on binding of the extreme N-terminus to the DEAD motor domains, NBD1 and NBD2, was used to model a novel parallel and flexible SecA solution dimer that agrees well with SAXS data.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

The magnetic structure of long-wavelength magnetohydrodynamic modes in current-carrying stellarator plasmas

Magnetic field fluctuations are observed in current-carrying stellarator plasmas when the rotational transform is close to a rational value at the edge of the plasma. At low plasma pressure, these fluctuations are associated with perturbed currents parallel to the equilibrium magnetic field lines. A model for these magnetohydrodynamic modes in a low-β, three-dimensional stellarator equilibria has been developed. A set of helical current filaments are constrained to mimic the structure of magnetic field lines on rational surfaces derived from three-dimensional (3D) equilibrium reconstructions. Transformation to straight field line coordinates then allows fitting of the poloidal magnetic sensor data to a single harmonic function, which fixes the modeled toroidal mode structure via the field line flow geometry. The developed procedure accurately captures phase and amplitude variation for m/n = 3/2, 3/1, and 4/1 modes in the 3D equilibria of the compact toroidal hybrid experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Streamlined spatial and environmental expression signatures characterize the minimalist duckweed Wolffia australiana

Single-cell genomics permits a new resolution in the examination of molecular and cellular dynamics, allowing global, parallel assessments of cell types and cellular behaviors through development and in response to environmental circumstances, such as interaction with water and the light–dark cycle of the Earth. Here, we leverage the smallest, and possibly most structurally reduced, plant, the semiaquaticWolffia australiana, to understand dynamics of cell expression in these contexts at the whole-plant level. We examined single-cell-resolution RNA-sequencing data and foundWolffiacells divide into four principal clusters representing the above- and below-water-situated parenchyma and epidermis. Although these tissues share transcriptomic similarity with model plants, they display distinct adaptations thatWolffiahas made for the aquatic environment. Within this broad classification, discrete subspecializations are evident, with select cells showing unique transcriptomic signatures associated with developmental maturation and specialized physiologies. Assessing this simplified biological system temporally at two key time-of-day (TOD) transitions, we identify additional TOD-responsive genes previously overlooked in whole-plant transcriptomic approaches and demonstrate that the core circadian clock machinery and its downstream responses can vary in cell-specific manners, even in this simplified system. Distinctions between cell types and their responses to submergence and/or TOD are driven by expression changes of unexpectedly few genes, characterizingWolffiaas a highly streamlined organism with the majority of genes dedicated to fundamental cellular processes.Wolffiaprovides a unique opportunity to apply reductionist biology to elucidate signaling functions at the organismal level, for which this work provides a powerful resource.

Biochemistry & Molecular Biology↗

DOE OSTI2021

YAHFC: A Code Framework to Model Nuclear Reactions and Estimate Correlated Uncertainties

Reaction modeling is a key ingredient in designing experiments and interpreting their results, and is an essential component in the process of evaluating nuclear data and assembling nuclear data libraries used in nuclear technology applications. Typically, experimental data are available only for a handful of reaction channels and theory models are used to fill in the gaps. In addition, theory is often called upon as the arbitrator between discrepant data. Most importantly, theory and modeling are required for an accurate determination of uncertainties in the evaluated data and the correlations between the multiple channels. A fast, accurate, and flexible modeling capability has been developed at LLNL with the code system YAHFC (Yet Another Hauser-Feshbach Code). YAHFC is a Monte Carlo, Hauser-Feshbach code framework, making full use of dynamic memory allocation, derived types, and parallel computing. YAHFC can generate events to simulate experiments and is guiding experiments designed to measure inelastic neutron scattering from actinide targets. YAHFC is also being used to analyze decays from surrogate experiments, thereby enabling the inference of reaction cross sections inaccessible by direct measurement. Finally, by modeling nuclear reactions with constraints from experimental data, YAHFC can deliver complete nuclear data libraries, with evaluated uncertainties, using the modernized Generalized Nuclear Data Structure (GNDS).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Efficacy of the Cell Perturbation Method in Large-Eddy Simulations of Boundary Layer Flow over Complex Terrain

A challenge to simulating turbulent flow in multiscale atmospheric applications is the efficient generation of resolved turbulence motions over an area of interest. One approach is to apply small perturbations to flow variables near the inflow planes of turbulence-resolving simulation domains nested within larger mesoscale domains. While this approach has been examined in numerous idealized and simple terrain cases, its efficacy in complex terrain environments has not yet been fully explored. Here, we examine the benefits of the stochastic cell perturbation method (CPM) over real complex terrain using data from the 2017 Perdigão field campaign, conducted in an approximately 2-km wide valley situated between two nearly parallel ridges. Following a typical configuration for multiscale simulation using nested domains within the Weather Research and Forecasting (WRF) model to downscale from the mesoscale to a large-eddy simulation (LES), we apply the CPM on a domain with horizontal grid spacing of 150 m. At this resolution, spurious coherent structures are often observed under unstable atmospheric conditions with moderate mean wind speeds. Results from such an intermediate resolution grid are often nested down for finer, more detailed LES, where these spurious structures adversely affect the development of turbulence on the subsequent finer grid nest. We therefore examine the impacts of the CPM on the representation of turbulence within the nested LES domain under moderate mean flow conditions in three different stability regimes: weakly convective, strongly convective, and weakly stable. In addition, two different resolutions of the underlying terrain are used to explore the role of the complex topography itself in generating turbulent structures. We demonstrate that the CPM improves the representation of turbulence within the LES domain, relative to the use of high-resolution complex terrain alone. During the convective conditions, the CPM improves the rate at which smaller-scales of turbulence form, while also accelerating the attenuation of the spurious numerically generated roll structures near the inflow boundary. During stable conditions, the coarse mesh spacing of the intermediate LES domain used herein was insufficient to maintain resolved turbulence using CPM as the flow develops downstream, highlighting the need for yet higher resolution under even weakly stable conditions, and the importance of accurate representation of flow on intermediate LES grids.

17 WIND ENERGY↗

DOE OSTI2025

Randomized Preconditioned Solvers for Strong Constraint 4D-Var Data Assimilation

The Strong Constraint 4D Variational (SC-4DVAR) data assimilation method is widely used in climate and weather applications. SC-4DVAR involves solving a minimization problem to compute the maximum a posteriori estimate, which we tackle using the Gauss-Newton method. The computation of the descent direction is expensive since it involves the solution of a large-scale and potentially ill-conditioned linear system, solved using the preconditioned conjugate gradient (PCG) method. Here, to address this cost, we efficiently construct scalable preconditioners using three different randomization techniques, which all rely on a certain low-rank structure involving the Gauss-Newton Hessian. The proposed techniques come with theoretical guarantees on the condition number, and at the same time, are amenable to parallelization. We also develop an adaptive approach to estimate the sketch size and choose between the reuse or recomputation of the preconditioner. We demonstrate the performance and effectiveness of our methodology on two representative model problems—the Burgers and barotropic vorticity equation—showing a drastic reduction in both the number of PCG iterations and the number of Gauss-Newton Hessian products after including the preconditioner construction cost.

Gauss-Newton↗

DOE OSTI2023

Anticipating gelation and vitrification with medium amplitude parallel superposition (MAPS) rheology and artificial neural networks

Abstract Anticipating qualitative changes in the rheological response of complex fluids (e.g., a gelation or vitrification transition) is an important capability for processing operations that utilize such materials in real-world environments. One class of complex fluids that exhibits distinct rheological states are soft glassy materials such as colloidal gels and clay dispersions, which can be well characterized by the soft glassy rheology (SGR) model. We first solve the model equations for the time-dependent, weakly nonlinear response of the SGR model. With this analytical solution, we show that the weak nonlinearities measured via medium amplitude parallel superposition (MAPS) rheology can be used to anticipate the rheological aging transitions in the linear response of soft glassy materials. This is a rheological version of a technique called structural health monitoring used widely in civil and aerospace engineering. We design and train artificial neural networks (ANNs) that are capable of quickly inferring the parameters of the SGR model from the results of sequential MAPS experiments. The combination of these data-rich experiments and machine learning tools to provide a surrogate for computationally expensive viscoelastic constitutive equations allows for rapid experimental characterization of the rheological state of soft glassy materials. We apply this technique to an aging dispersion of Laponite ® clay particles approaching the gel point and demonstrate that a trained ANN can provide real-time detection of transitions in the nonlinear response well in advance of incipient changes in the linear viscoelastic response of the system.

Lennon, Kyle R. (ORCID:0000000212515461)↗

DOE OSTI2019

Picosecond Dynamical Response to a Pressure-Induced Break of the Tertiary Structure Hydrogen Bonds in a Membrane Chromoprotein

In this work, we used elastic incoherent neutron scattering (EINS) to find out if structural changes accompanying local hydrogen bond rupture are also reflected in global dynamical response of the protein complex. Chromatophore membranes from LH2-only strains of the photosynthetic bacterium Rhodobacter sphaeroides, with spheroidenone or neurosporene as the major carotenoids, were subjected to high hydrostatic pressure at ambient temperature. Optical spectroscopy conducted at high pressure confirmed rupture of tertiary structure hydrogen bonds. In parallel, we used EINS to follow average motions of the hydrogen atoms in LH2, which reflect the flexibility of this complex. A decrease of the average atomic mean square displacements of hydrogen atoms was observed up to a pressure of 5 kbar in both carotenoid samples due to general stiffening of protein structures, while at higher pressures a slight increase of the displacements was detected in the neurosporene mutant LH2 sample only. These data reflect a correlation between the local pressure-induced breakage of H-bonds, observed in optical spectra, with the altered protein dynamics monitored by EINS. The slightly higher compressibility of the neurosporene mutant sample shows that even subtle alterations of carotenoids are manifested on a larger scale and emphasize a close connection between the local structure and global dynamics of this membrane protein complex.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Impacts of Sub‐Grid Topographic Representations on Surface Energy Balance and Boundary Conditions in the E3SM Land Model: A Case Study in Sierra Nevada

Abstract Sub‐grid topographic heterogeneity has large impacts on surface energy balance and land‐atmosphere interactions. However, the impacts of representing sub‐grid topographic effects in land surface models (LSMs) on surface energy balance and boundary conditions remain unclear. This study analyzed and evaluated the impacts of sub‐grid topographic representations on surface energy balance, turbulent heat flux, and scalar (co‐)variances in the Energy Exascale Earth System Model (E3SM) land model (ELM). Three sub‐grid topographic representations in ELM were compared: (a) the default sub‐grid structure (D), (b) the recently developed sub‐grid topographic structure (T), and (c) high spatial resolution (1KM). Additionally, two different solar radiation schemes in ELM were compared: (a) the default plane‐parallel radiative transfer scheme (PP) and (b) the parameterization scheme (TOP) that accounts for sub‐grid topographic effects on solar radiation. A series of offline simulations with the three grid discretization structures (D, T, and 1KM) and two schemes of solar radiation (TOP and PP) were carried out using the Sierra Nevada, California. 1KM simulations with TOP well capture the spatial heterogeneity of surface fluxes compared to Moderate Resolution Imaging Spectroradiometer remote sensing data. There are significant differences between TOP and PP in the 1‐km simulated surface energy balance, but the differences in mean values and standard deviations become small when aggregated to the grid scale (i.e., 0.5°). The T configuration better mimics the 1KM simulations with TOP than the D configuration and better captures the sub‐grid topographic effects on surface energy balance and boundary conditions. These results underline the importance of representing sub‐grid topographic heterogeneities in LSMs and motivate future research to understand the sub‐grid topographic effects on land‐atmosphere interactions over mountainous areas.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2024

Low-cost sensor-based damage localization for large-area monitoring of FRP composites

Abstract In recent years, there has been growing interest in self-sensing structural materials across research and industry sectors. Detecting and locating structural damage typically requires numerous sensors wired to a data acquisition (DAQ) circuit, rendering implementation impractical in real structures. This paper proposes an innovative, cost-effective sensor network for damage detection and localization in fiber-reinforced polymer composites. The innovation encompasses three key elements: (1) utilizing carbon fiber tows within the composite as piezoresistive sensors, eliminating the need for additional foreign sensor devices; (2) introducing a novel sensor layout wherein sensor tow branches with varied resistance values are connected in parallel, reducing the number of connections to the DAQ circuit and cutting manufacturing costs significantly; (3) developing a practical sensor terminal fabrication technique to minimize manufacturing expenses. The proposed design methodology for the branch resistance values is first validated using a demonstration panel. Subsequently, the overall strategy is assessed by conducting impact tests on carbon and glass fiber-reinforced composite specimens. Results validate the sensor’s ability to accurately detect and locate structural damage.

Demo, Luke B. (ORCID:0000000208806522)↗

DOE OSTI2023

Phyllosphere exudates select for distinct microbiome members in sorghum epicuticular wax and aerial root mucilage

Phyllosphere exudates create specialized microhabitats that shape microbial community diversity. We explored the microbiome associated with two sorghum phyllosphere exudates, the epicuticular wax and aerial root mucilage. We assessed the microbiome associated with the wax from sorghum plants over two growth stages, and the root mucilage additionally from nitrogen-fertilized and non-fertilized plants. In parallel, we isolated and characterized hundreds of bacteria from wax and mucilage, and integrated data from cultivation-independent and cultivation-dependent approaches to gain insights into exudate diversity and bacterial phenotypes. We found that Sphingomonadaceae and Rhizobiaceae families were the major taxa in the wax regardless of water availability and plant developmental stage to plants. The cultivation-independent mucilage-associated bacterial microbiome contained Erwiniaceae, Flavobacteriaceae, Rhizobiaceae, Pseudomonadaceae, Sphingomonadaceae, and its structure was strongly influenced by sorghum development but only modestly influenced by fertilization. In contrast, the fungal community structure of mucilage was strongly affected by the year of sampling but not by fertilization or plant developmental stage, suggesting a decoupling of fungal-bacterial dynamics in the mucilage. Our bacterial isolate collection from wax and mucilage had several isolates that matched 100% to detected amplicon sequence variants, and were enriched on media that selected for phenotypes including phosphate solubilization, putative diazotrophy, resistance to desiccation, capability to grow on methanol as a carbon source, and ability to grow in the presence of linalool and β-caryophyllene (terpenes in sorghum wax). This work expands our understanding of the microbiome of phyllosphere exudates and supports our long-term goal to translate microbiome research to support sorghum cultivation.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2026

Atomic‐Scale Surface Imaging of Bulk Epitaxial CsPbBr 3 Perovskite Single Crystals on Mica Using Light Assisted Scanning Tunneling Microscopy at Low‐Temperature (80 K)

Epitaxial single-crystalline CsPbBr 3 perovskite films on mica, prepared ex situ, are explored using a low-temperature scanning tunneling microscope (STM) by probing the unoccupied electronic states of their surface in ultra-high vacuum (UHV) at 80 K. Light-assisted STM measurements under a broadband illumination with visible light were employed to enhance and stabilize surface conductivity. STM imaging across the surface of macroscopic bulk CsPbBr 3 films reveals large flat terraces characterized by a specific type of surface reconstruction, consisting of parallel rows of U-shaped atomic nanostructures. These structures are spaced by 12 Å and exhibit an internal feature size of 5.1 Å. Density functional theory (DFT) calculations reproduce the experimental observations and reveal a competition between different orthorhombic CsPbBr 3 (110) surface reconstructions: a Cs-rich structure, identified as the most energetically stable, and three alternative Pb–Br-rich reconstructions, which are slightly higher in energy yet remain consistent with the STM data. Additional analyses that explicitly account for the mica substrate exclude the cubic CsPbBr 3 phase and other orthorhombic surface orientations, while showing that variations in the mica surface termination do not alter the preferred CsPbBr 3 (110) reconstruction. In conclusion, this combined approach thereby confirms our assignment and resolves previous STM interpretations of CsPbBr 3 .

36 MATERIALS SCIENCE↗

DOE OSTI2017

Building Halo Merger Trees from the Q Continuum Simulation

Cosmological N-body simulations rank among the most computationally intensive efforts today. A key challenge is the analysis of structure, substructure, and the merger history for many billions of compact particle clusters, called halos. Effectively representing the merging history of halos is essential for many galaxy formation models used to generate synthetic sky catalogs, an important application of modern cosmological simulations. Generating realistic mock catalogs requires computing the halo formation history from simulations with large volumes and billions of halos over many time steps, taking hundreds of terabytes of analysis data. We present fast parallel algorithms for producing halo merger trees and tracking halo substructure from a single-level, density-based clustering algorithm. Merger trees are created from analyzing the halo-particle membership function in adjacent snapshots, and substructure is identified by tracking the "cores" of merging halos - sets of particles near the halo center. Core tracking is performed after creating merger trees and uses the relationships found during tree construction to associate substructures with hosts. The algorithms are implemented with MPI and evaluated on a Cray XK7 supercomputer using up to 16,384 processes on data from HACC, a modern cosmological simulation framework. We present results for creating merger trees from 101 analysis snapshots taken from the Q Continuum, a large volume, high mass resolution, cosmological simulation evolving half a trillion particles.

79 ASTRONOMY AND ASTROPHYSICS↗