Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Random Effects”

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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Synthesis and property of EPEG ‐based polycarboxylate ether superplasticizers via RAFT polymerization

Abstract Polycarboxylate ether superplasticizers (PCEs) with comb‐shaped molecular structures have attracted considerable attention in the area of admixtures for construction engineering. However, there are some contradictions in studying the relationship between the molecular structure and property of PCE due to its random structure based on traditional free radical polymerization, which is confused for designing and developing novel PCEs. In this article, a series of reversible addition fragmentation chain transfer (RAFT) polymerization‐based polycarboxylate ether superplasticizers (RPCs) were prepared by the copolymerization of acrylic acid (AA) and ethylene‐glycol monovinyl polyethylene glycol (EPEG). The chemical structure and composition of RPCs were characterized and their adsorption behavior and dispersion properties were conducted. The results demonstrate that the molecular weights of RPCs do not vary significantly with prolonging reaction time, and the actual monomer ratio in RPCs displays a significant change to the feed one due to the monomer reactivity and conversion at higher reaction temperature, which is different from conventional free radical polymerization. The dispersion and adsorption properties of RPCs to cement were evaluated by slurry fluidity tests and total organic carbon (TOC) experiments. The effect of polymerization temperature and time on the dispersion ability has been systematically studied. The results indicate that the dispersion ability and adsorption capacity of RPCs are closely related to their molecular structures. The dispersion ability of RPCs were enhanced with increasing carboxyl group density in a certain range of monomer molar ratios (≤6:1). This study will provide a new platform for designing and developing high‐performance and rational cost‐effective superplasticizers with well‐defined structure and controlled molecular weight.

Chen, Xiaodong↗
DOE OSTI2024

Microbial inoculum effects on the rumen epithelial transcriptome and rumen epimural metatranscriptome in calves

Manipulation of the rumen microbial ecosystem in early life may affect ruminal fermentation and enhance the productive performance of dairy cows. The objective of this experiment was to evaluate the effects of dosing three different types of microbial inoculum on the rumen epithelium tissue (RE) transcriptome and the rumen epimural metatranscriptome (REM) in dairy calves. For this objective, 15 Holstein bull calves were enrolled in the study at birth and assigned to three different intraruminal inoculum treatments dosed orally once weekly from three to six weeks of age. The inoculum treatments were prepared from rumen contents collected from rumen fistulated lactating cows and were either autoclaved (control; ARF), processed by differential centrifugation to create the bacterial-enriched inoculum (BE), or through gravimetric separation to create the protozoal-enriched inoculum (PE). Calves were fed 2.5 L/d pasteurized waste milk 3x/d from 0 to 7 weeks of age and texturized starter until euthanasia at 9 weeks of age, when the RE tissues were collected for transcriptome and microbial metatranscriptome analyses, from four randomly selected calves from each treatment. The different types of inoculum altered the RE transcriptome and REM. Compared to ARF, 9 genes were upregulated in the RE of BE and 92 in PE, whereas between BE and PE there were 13 genes upregulated in BE and 114 in PE. Gene ontology analysis identified enriched GO terms in biological process category between PE and ARF, with no enrichment between BE and ARF. The RE functional signature showed different KEGG pathways related to BE and ARF, and no specific KEGG pathway for PE. We observed a lower alpha diversity index for RE microbiome in ARF (observed genera and Chao1 (p < 0.05)). Five microbial genera showed a significant correlation with the changes in host gene expression: Roseburia (25 genes), Entamoeba (two genes); Anaerosinus, Lachnospira, and Succiniclasticum were each related to one gene. sPLS-DA analysis showed that RE microbial communities differ among the treatments, although the taxonomic and functional microbial profiles show different distributions. Co-expression Differential Network Analysis indicated that both BE and PE had an impact on the abundance of KEGG modules related to acyl-CoA synthesis, type VI secretion, and methanogenesis, while PE had a significant impact on KEGGs related to ectoine biosynthesis and D-xylose transport. Our study indicated that artificial dosing with different microbial inocula in early life alters not only the RE transcriptome, but also affects the REM and its functions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A Robust Method to Secure Multi-Inverter Grid Tied PV and Battery Energy Storage Systems Against Cyber Intrusions

This paper details a robust method to secure a multi-inverter grid tied system that interfaces photovoltaic (PV) and battery energy storage against potential cyber-attacks. The method can be applied to any third-party inverter systems without a need to modify their internal controls. A small random private excitation signal termed "watermark" is injected into the DC input voltage terminals (via a series transformer) connected to the PV/battery inverter system. An external robust cyber intrusion detector (CID) hardware consisting of a digital signal processor (DSP) generates the "watermark" and also receives the sensor signals that control the setpoints of the PV/battery grid tied system. The CID algorithm is shown to detect all possible cyber intrusions (such as false data injection(FDI)) on external sensor signals such as P and Q measured by a smart meter that control the overall system operation. The proposed CID computes online system ID and two variance tests in real time on each sensor signal and is able pinpoint intrusion location in a multi-inverter system. Results on a hardware in the loop (HIL) of a two-inverter grid connected system demonstrate effectiveness of the proposed CID system for FDI and unobservable FDI. Test results on a laboratory prototype will be discussed in the conference presentation.

Ibrahim, Hasan↗
DOE OSTI2024

Effects of spatial variability in vegetation phenology, climate, landcover, biodiversity, topography, and soil property on soil respiration across a coastal ecosystem

Coastal terrestrial-aquatic interfaces (TAIs) are crucial contributors to global biogeochemical cycles and carbon exchange. A systematic evaluation of the interaction between coastal catchment properties and carbon dioxide (CO2) emission by soil respiration is significant for assessing carbon dynamics and predicting the future trajectory of atmospheric CO2 concentrations in coastal TAIs. The soil CO2 efflux in these transition zones is however poorly understood due to the high spatiotemporal dynamics of TAIs, as various sub-ecosystems in this region are compressed and expanded by complex influences of tides, changes in river levels, climate, and land use. We focus on the Chesapeake Bay region to (i) investigate the spatial heterogeneity of the coastal ecosystem and identify spatial zones with similar environmental characteristics based on the spatial data layers, including vegetation index (kNDVI), climate, landcover, diversity, topography, soil property, and relative tidal elevation; (ii) understand the primary driving factors affecting soil respiration within sub-ecosystems of the coastal ecosystem. Specifically, we employed hierarchical clustering analysis to identify spatial regions with distinct environmental characteristics, followed by the determination of main driving factors using Random Forest regression and SHapley Additive exPlanations. Maximum and minimum temperature are the main drivers common to all sub-ecosystems, while each region also has additional unique major drivers that differentiate them from one another. Precipitation exerts an influence on vegetated lands, while soil pH value holds importance specifically in forested lands. In croplands characterized by high clay content and low sand content, the significant role is attributed to bulk density. Wetlands demonstrate the importance of both elevation and sand content, with clay content being more relevant in non-inundated wetlands than in inundated wetlands. The topographic wetness index significantly contributes to the mixed vegetation areas, including shrub, grass, pasture, and forest. Additionally, our research reveals that dense vegetation land covers and urban/developed areas exhibit distinct soil property drivers. Overall, there is no one-size-fits-all approach to modeling carbon fluxes in coastal TAIs, and our study highlights the importance of further research and monitoring practices to improve our understanding of carbon dynamics and promote the sustainable management of coastal TAIs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Parallelized domain decomposition for multi-dimensional Lagrangian random walk mass-transfer particle tracking schemes

Lagrangian particle tracking schemes allow a wide range of flow and transport processes to be simulated accurately, but a major challenge is numerically implementing the inter-particle interactions in an efficient manner. This article develops a multi-dimensional, parallelized domain decomposition (DDC) strategy for mass-transfer particle tracking (MTPT) methods in which particles exchange mass dynamically. We show that this can be efficiently parallelized by employing large numbers of CPU cores to accelerate run times. In order to validate the approach and our theoretical predictions we focus our efforts on a well-known benchmark problem with pure diffusion, where analytical solutions in any number of dimensions are well established. In this work, we investigate different procedures for “tiling” the domain in two and three dimensions (2-D and 3-D), as this type of formal DDC construction is currently limited to 1-D. An optimal tiling is prescribed based on physical problem parameters and the number of available CPU cores, as each tiling provides distinct results in both accuracy and run time. We further extend the most efficient technique to 3-D for comparison, leading to an analytical discussion of the effect of dimensionality on strategies for implementing DDC schemes. Increasing computational resources (cores) within the DDC method produces a trade-off between inter-node communication and on-node work. For an optimally subdivided diffusion problem, the 2-D parallelized algorithm achieves nearly perfect linear speedup in comparison with the serial run-up to around 2700 cores, reducing a 5 h simulation to 8 s, while the 3-D algorithm maintains appreciable speedup up to 1700 cores.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

The Effect of the Prior and the Experimental Design on the Inference of the Precision Matrix in Gaussian Chain Graph Models

Here, we investigate whether (and how) experimental design could aid in the estimation of the precision matrix in a Gaussian chain graph model, especially the interplay between the design, the effect of the experiment and prior knowledge about the effect. Estimation of the precision matrix is a fundamental task to infer biological graphical structures like microbial networks. We compare the marginal posterior precision of the precision matrix under four priors: flat, conjugate Normal-Wishart, Normal-MGIG and a general independent. Under the flat and conjugate priors, the Laplace-approximated posterior precision is not a function of the design matrix rendering useless any efforts to find an optimal experimental design to infer the precision matrix. In contrast, the Normal-MGIG and general independent priors do allow for the search of optimal experimental designs, yet there is a sharp upper bound on the information that can be extracted from a given experiment. We confirm our theoretical findings via a simulation study comparing (i) the KL divergence between prior and posterior and (ii) the Stein’s loss difference of MAPs between random and no experiment. Our findings provide practical advice for domain scientists conducting experiments to better infer the precision matrix as a representation of a biological network.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Searching for quasar candidates with periodic variations from the Zwicky Transient Facility: results and implications

ABSTRACT We conduct a systematic search for quasars with periodic variations from the archival photometric data of the Zwicky Transient Facility by cross-matching with the quasar catalogues of the Sloan Digital Sky Survey and Véron-Cetty and Véron. We first select out 184 primitive periodic candidates using the generalized Lomb–Scargle periodogram and autocorrelation function and then estimate their statistical significance of periodicity based on two red-noise models, i.e. damped random walk (DRW) and single power-law (SPL) models. As such, we finally identify 106 (DRW) and 86 (SPL) candidates with the most significant periodic variations out of 143 700 quasars. We further compare DRW and SPL models using Bayes factors, which indicate a relative preference of the SPL model for our primitive sample. We thus adopt the candidates identified with SPL as the final sample and summarize its basic properties. We extend the light curves of the selected candidates by supplying other archival survey data to verify their periodicity. However, only three candidates (with 6–8 cycles of periods) meet the selection criteria. This result clearly implies that, instead of being strictly periodic, the variability must be quasi-periodic or caused by stochastic red-noise. This exerts a challenge to the existing search approaches and calls for developing new effective methods.

Chen, Yong-Jie (ORCID:0000000342807673)↗
DOE OSTI2023

Development of an interatomic potential for the W–Ta system

A physics-inspired, data-driven interatomic potential framework for multi-element system is presented. This potential is based on the generalization of the embedded atom method potential and systematically incorporates two-, three- and many-body effects. Different atomic environments are described by atom-centered Gaussian basis sets which provides sufficient flexibility to capture diverse atomic environments and allows for easy optimization of the free parameters. An interatomic potential model for the tungsten–tantalum (W–Ta) system is developed using this framework by training on data from ab initio density functional theory (DFT) calculations. The potential is thoroughly tested at various compositions by comparing bulk and defect properties from experimental and DFT data. It is shown that the potential predicts the elastic constants, defect properties, such as vacancy and interstitial formation energies, core structures of dislocations, and melting points for both pure and tantalum and tungsten with an accuracy comparable to other single element machine learning based potentials. In conclusion, the potential model is used to investigate the formation energies of ordered alloys and vacancy formation energies for chemically random W–Ta alloys.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

First estimation of high-resolution solar photovoltaic resource maps over China with Fengyun-4A satellite and machine learning

Fengyun-4A (FY-4A), which is the latest-generation Chinese geostationary meteorological satellite, measures solar reflection and thermal emission with high temporal, spatial, and spectral resolutions. It is expected to be highly beneficial for solar resource assessment and forecasting in China. This study is the first to estimate, using FY-4A and a random forest model, the global horizontal irradiance (GHI) at a 4-km–15-min spatio-temporal resolution over China, as a means to arrive at a solar photovoltaic (PV) resource map. In terms of GHI estimates, the root mean square error and mean bias error between hourly measured and retrieved values are 147.02 (35.2%), –5.64 W/m 2 (–1.4%), respectively, whereas the values of daily estimates are 29.20 (18.0%), –2.97 W/m 2 (–1.3%). The retrieval accuracy is found much better for instances with solar zenith angles smaller than 60°. Relatively larger errors are found at locations in the Sichuan Basin and northeastern China, which can be attributed to bright surfaces and/or strong cloud transients. With the retrieved irradiance, PV resource is derived through a physical model chain. The annual mean PV resource map suggests that, over most of the west areas, the annual mean effective irradiance exceeds 1700 kWh/m 2 , with the highest value found in Tibet (around 2000 kWh/m 2 per annum). Eastern China has an annual effective irradiance of only 1300–1500 kWh/m 2 . In conclusion, the region with poorest solar resource is the Sichuan Basin (less than 1100 kWh/m 2 per annum).

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Effects of Size and Shape on the Tolerances for Misalignment and Probabilities for Successful Oriented Attachment of Nanoparticles

Oriented attachment (OA) of nanoparticles is an important pathway of crystal growth, but tools for quantitatively modeling OA are lacking. Here we present several simple models that relate the probability of achieving OA to basic geometric parameters such as particle size, shape, and lattice periodicity. A Moiré-domain model is applied to understand twist-misorientations between parallel surfaces, and it predicts that the range of twist angles yielding perfect OA is inversely related to the width of the contact area. This is confirmed and further developed using a surface functional model, which predicts how crystallographic registration forces drive the emergence of complex orientational energy landscapes. The energy landscapes are predicted to possess local minima that can trap particles in imperfect alignments, and these local minima become deeper and more numerous as the contact area increases, making OA more challenging for large particles. Further, a second set of models is presented to understand the sequence of events by which two crystallographic faces become co-planer after collision. We use a ‘central force approximation’ to quantitatively predict the odds of attaining coalignment between various faces when particles collide with random misalignments, and we show that in the absence of biasing forces, the probability of attaining alignment on a given face is roughly proportional to its solid angle as viewed from the center of the particle. The model predicts that OA is most favorable between well-faceted particles and becomes exceedingly unlikely for large spherical particles that express many microfacets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Extended Interacting Boson Model with Configuration Mixing Description For 98 Mo

An extended version of the interacting boson model (IBM), encompassing two-particle–two-hole (2p-2h) excitations originating from below a valence shell and one-particle–one-hole (1p-1h) excitations extending to the upper next-nearest major shell, is developed to investigate cross-shell excitation effects in the transitional nucleus 98 Mo. This extended IBM framework is utilized to reproduce experimentally observed properties such as low-lying positive-parity level energies, 𝐵⁡(𝐸⁢2) values, and the electric quadrupole moment of the 2$^{+}_{1}$ state, 𝜌 2 ⁡(𝐸⁢0, 0$^{+}_{2}$ → 0$^{+}_{g}$), as well as the strength functions of the isoscalar giant monopole resonance (ISGMR) and the isoscalar giant quadrupole resonance (ISGQR) in 98 Mo. It is shown that the extended consistent-𝑄 IBM, incorporating configuration mixing of both 2p-2h excitations from below the valence shell and 1p-1h excitations to the upper next-nearest major shell, describes the low-lying structural properties and the ISGMR and ISGQR strength distributions of 98 Mo quite well, with the ISGMR and the ISGQR strength distributions comparable to the results of the self-consistent quasiparticle random-phase approximation with Skyrme interactions.

90 ≤ A ≤ 149↗
DOE OSTI2022

Disentangling direct from indirect relationships in association networks

Networks are vital tools for understanding and modeling interactions in complex systems in science and engineering, and direct and indirect interactions are pervasive in all types of networks. However, quantitatively disentangling direct and indirect relationships in networks remains a formidable task. Here, we present a framework, called iDIRECT (Inference of Direct and Indirect Relationships with Effective Copula-based Transitivity), for quantitatively inferring direct dependencies in association networks. Using copula-based transitivity, iDIRECT eliminates/ameliorates several challenging mathematical problems, including ill-conditioning, self-looping, and interaction strength overflow. With simulation data as benchmark examples, iDIRECT showed high prediction accuracies. Application of iDIRECT to reconstruct gene regulatory networks in Escherichia coli also revealed considerably higher prediction power than the best-performing approaches in the DREAM5 (Dialogue on Reverse Engineering Assessment and Methods project, #5) Network Inference Challenge. In addition, applying iDIRECT to highly diverse grassland soil microbial communities in response to climate warming showed that the iDIRECT-processed networks were significantly different from the original networks, with considerably fewer nodes, links, and connectivity, but higher relative modularity. Further analysis revealed that the iDIRECT-processed network was more complex under warming than the control and more robust to both random and target species removal ( P < 0.001). As a general approach, iDIRECT has great advantages for network inference, and it should be widely applicable to infer direct relationships in association networks across diverse disciplines in science and engineering.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Resource Occurrence and Productivity in Existing and Proposed Wind Energy Lease Areas on the Northeast US Shelf

States in the Northeast United States have the ambitious goal of producing more than 22 GW of offshore wind energy in the coming decades. The infrastructure associated with offshore wind energy development is expected to modify marine habitats and potentially alter the ecosystem services. Species distribution models were constructed for a group of fish and macroinvertebrate taxa resident in the Northeast US Continental Shelf marine ecosystem. These models were analyzed to provide baseline context for impact assessment of lease areas in the Middle Atlantic Bight designated for renewable wind energy installations. Using random forest machine learning, models based on occurrence and biomass were constructed for 93 species providing seasonal depictions of their habitat distributions. We developed a scoring index to characterize lease area habitat use for each species. Subsequently, groups of species were identified that reflect varying levels of lease area habitat use ranging across high, moderate, low, and no reliance on the lease area habitats. Among the species with high to moderate reliance were black sea bass ( Centropristis striata ), summer flounder ( Paralichthys dentatus ), and Atlantic menhaden ( Brevoortia tyrannus ), which are important fisheries species in the region. Potential for impact was characterized by the number of species with habitat dependencies associated with lease areas and these varied with a number of continuous gradients. Habitats that support high biomass were distributed more to the northeast, while high occupancy habitats appeared to be further from the coast. There was no obvious effect of the size of the lease area on the importance of associated habitats. Model results indicated that physical drivers and lower trophic level indicators might strongly control the habitat distribution of ecologically and commercially important species in the wind lease areas. Therefore, physical and biological oceanography on the continental shelf proximate to wind energy infrastructure development should be monitored for changes in water column structure and the productivity of phytoplankton and zooplankton and the effects of these changes on the trophic system.

17 WIND ENERGY↗
DOE OSTI2021

Thermodynamics of order and randomness in dopant distributions inferred from atomically resolved imaging

Abstract Exploration of structure-property relationships as a function of dopant concentration is commonly based on mean field theories for solid solutions. However, such theories that work well for semiconductors tend to fail in materials with strong correlations, either in electronic behavior or chemical segregation. In these cases, the details of atomic arrangements are generally not explored and analyzed. The knowledge of the generative physics and chemistry of the material can obviate this problem, since defect configuration libraries as stochastic representation of atomic level structures can be generated, or parameters of mesoscopic thermodynamic models can be derived. To obtain such information for improved predictions, we use data from atomically resolved microscopic images that visualize complex structural correlations within the system and translate them into statistical mechanical models of structure formation. Given the significant uncertainties about the microscopic aspects of the material’s processing history along with the limited number of available images, we combine model optimization techniques with the principles of statistical hypothesis testing. We demonstrate the approach on data from a series of atomically-resolved scanning transmission electron microscopy images of Mo x Re 1- x S 2 at varying ratios of Mo/Re stoichiometries, for which we propose an effective interaction model that is then used to generate atomic configurations and make testable predictions at a range of concentrations and formation temperatures.

25 ENERGY STORAGE↗
DOE OSTI2022

Geochemistry and Multiomics Data Differentiate Streams in Pennsylvania Based on Unconventional Oil and Gas Activity

Unconventional oil and gas (UOG) extraction is increasing exponentially around the world, as new technological advances have provided cost-effective methods to extract hard-to-reach hydrocarbons. While UOG has increased the energy output of some countries, past research indicates potential impacts in nearby stream ecosystems as measured by geochemical and microbial markers. Here, we utilized a robust data set that combines 16S rRNA gene amplicon sequencing (DNA), metatranscriptomics (RNA), geochemistry, and trace element analyses to establish the impact of UOG activity in 21 sites in northern Pennsylvania. These data were also used to design predictive machine learning models to determine the UOG impact on streams. We identified multiple biomarkers of UOG activity and contributors of antimicrobial resistance within the order Burkholderiales. Furthermore, we identified expressed antimicrobial resistance genes, land coverage, geochemistry, and specific microbes as strong predictors of UOG status. Of the predictive models constructed (n = 30), 15 had accuracies higher than expected by chance and area under the curve values above 0.70. The supervised random forest models with the highest accuracy were constructed with 16S rRNA gene profiles, metatranscriptomics active microbial composition, metatranscriptomics active antimicrobial resistance genes, land coverage, and geochemistry (n = 23). The models identified the most important features within those data sets for classifying UOG status. These findings identified specific shifts in gene presence and expression, as well as geochemical measures, that can be used to build robust models to identify impacts of UOG development.

16S rRNA↗
DOE OSTI2022

Hot Forging Options for Thick Castings

Uranium alloyed with 10 wt% molybdenum (U-10Mo) is a monolithic nuclear fuel relevant to the National Nuclear Security Administration’s nonproliferation efforts. Research has been underway to optimize processing techniques for the U-10Mo fuel. This study investigated the use of hot compression or “hot forging” on a thick (~1") cast and homogenized U-10Mo plate before standard hot and cold rolling procedures. After plates were cast and homogenized, six samples were cut and forged at 700°C at a strain rate of either 0.10 s –1 or 0.01 s –1 at six levels of reduction. After forging, all samples underwent hot and cold rolling followed by annealing to achieve a final foil thickness of about 0.0085". Samples were taken at each stage of the casting and thermomechanical processing to assess the microstructural evolution. Chemical composition, microstructure, and uranium carbide morphology are presented and assessed in this study. Upon hot forging, dislocations accumulate along the grain boundaries, which serve as nucleation sites for randomly oriented, strain-free grains during subsequent annealing steps. Hot forging and subsequent annealing produced very heterogeneous grain sizes. However, no molybdenum segregation was observed after forging. No obvious trend was observed between forging conditions (strain rate and reduction percentage) and the microstructure after final thermomechanical processing. Upon hot rolling to 0.04" and annealing (700°C for 45 min), the average grain diameters from OM was 17 ± 2 µm across the six different forged samples. The subsequent cold rolling to 0.0085" and then annealing (700°C for 45 min) resulted in an average of 13 ± 2 µm between the six samples. Thus, the starting, as-forged microstructure did not appear to significantly influence the final microstructure of the cold-rolled foil. These results will help with understanding and expanding hot working capabilities for thicker U-10Mo castings. They also provide useful information on the effects of hot forging and its potential use to minimize defects that can arise during subsequent hot and cold rolling procedures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Optimization of Thermal Conductance at Interfaces Using Machine Learning Algorithms

We report optimization of thermal transport across the interface of two different materials is critical to micro-/nanoscale electronic, photonic, and phononic devices. Although several examples of compositional intermixing at the interfaces having a positive effect on interfacial thermal conductance (ITC) have been reported, an optimum arrangement has not yet been determined because of the large number of potential atomic configurations and the significant computational cost of evaluation. On the other hand, computation-driven materials design efforts are rising in popularity and importance. Yet, the scalability and transferability of machine learning models remain as challenges in creating a complete pipeline for the simulation and analysis of large molecular systems. In this work we present a scalable Bayesian optimization framework, which leverages dynamic spawning of jobs through the Message Passing Interface (MPI) to run multiple parallel molecular dynamics simulations within a parent MPI job to optimize heat transfer at the silicon and aluminum (Si/Al) interface. We found a maximum of 50% increase in the ITC when introducing a two-layer intermixed region that consists of a higher percentage of Si. Because of the random nature of the intermixing, the magnitude of increase in the ITC varies. We observed that both homogeneity/heterogeneity of the intermixing and the intrinsic stochastic nature of molecular dynamics simulations account for the variance in ITC.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Entropy Stabilization Effects and Ion Migration in 3D “Hollow” Halide Perovskites

A recently discovered new family of 3D halide perovskites with the general formula (A) 1–x (en) x (Pb) 1–0.7x (X) 3–0.4x (A = MA, FA; X = Br, I; MA = methylammonium, FA = formamidinium, en = ethylenediammonium) is referred to as “hollow” perovskites owing to extensive Pb and X vacancies created on incorporation of en cations in the 3D network. The “hollow” motif allows fine tuning of optical, electronic, and transport properties and bestowing good environmental stability proportional to en loading. To shed light on the origin of the apparent stability of these materials, we performed detailed thermochemical studies, using room temperature solution calorimetry combined with density functional theory simulations on three different families of “hollow” perovskites namely en/FAPbI 3 , en/MAPbI 3 , and en/FAPbBr 3 . We found that the bromide perovskites are more energetically stable compared to iodide perovskites in the FA-based hollow compounds, as shown by the measured enthalpies of formation and the calculated formation energies. The least stable FAPbI 3 gains stability on incorporation of the en cation, whereas FAPbBr 3 becomes less stable with en loading. This behavior is attributed to the difference in the 3D cage size in the bromide and iodide perovskites. Configurational entropy, which arises from randomly distributed cation and anion vacancies, plays a significant role in stabilizing these “hollow” perovskite structures despite small differences in their formation enthalpies. With the increased vacancy defect population, we have also examined halide ion migration in the FA-based “hollow” perovskites and found that the migration energy barriers become smaller with the increasing en content.

36 MATERIALS SCIENCE↗