Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STRESS FUNCTION”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

CHGNet as a pretrained universal neural network potential for charge-informed atomistic modelling

Large-scale simulations with complex electron interactions remain one of the greatest challenges for atomistic modelling. Although classical force fields often fail to describe the coupling between electronic states and ionic rearrangements, the more accurate ab initio molecular dynamics suffers from computational complexity that prevents long-time and large-scale simulations, which are essential to study technologically relevant phenomena. Here we present the Crystal Hamiltonian Graph Neural Network (CHGNet), a graph neural network-based machine-learning interatomic potential (MLIP) that models the universal potential energy surface. CHGNet is pretrained on the energies, forces, stresses and magnetic moments from the Materials Project Trajectory Dataset, which consists of over 10 years of density functional theory calculations of more than 1.5 million inorganic structures. The explicit inclusion of magnetic moments enables CHGNet to learn and accurately represent the orbital occupancy of electrons, enhancing its capability to describe both atomic and electronic degrees of freedom. We demonstrate several applications of CHGNet in solid-state materials, including charge-informed molecular dynamics in Li x MnO 2 , the finite temperature phase diagram for Li x FePO 4 and Li diffusion in garnet conductors. We highlight the significance of charge information for capturing appropriate chemistry and provide insights into ionic systems with additional electronic degrees of freedom that cannot be observed by previous MLIPs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Non-dilute helium-related defect interactions in the near-surface region of plasma-exposed tungsten

We report a systematic energetic analysis of helium-related defect interactions that mediate helium (He) segregation on surfaces of plasma-exposed tungsten at different levels of He ion implantation. We focus on high He fluences that increase the He content in the plasma-exposed material well beyond the dilute limit of He concentration and employ atomic configurations generated by large-scale molecular dynamics simulations of He-implanted tungsten. We perform systematic molecular statics computations of cluster– defect interaction energetics in the highly defect-rich near-surface region of plasma-exposed tungsten for small mobile helium clusters as a function of the cluster distance from the surface. In this region, mobile clusters are also subjected to the stress fields generated by defects such as helium bubbles and other clusters, which govern cluster–defect interactions in addition to the cluster–surface interaction. Based on systematic investigation of individual cluster–defect interactions, we develop a mathematical framework to describe the interaction energy landscapes consisting of elastic interaction potential perturbations to finite-width square-well potentials, where the potential well accounts for cluster trapping by the defect at close range and subsequent coalescence and the perturbation potential is parameterized according to elastic inclusion theory. Superposition of all the relevant interaction potentials provides a comprehensive description of the interaction energy landscape that would be experienced by a small mobile cluster along its migration path toward the plasma-exposed surface at high He fluence. Such descriptions are particularly important for developing atomistically-informed, hierarchical multi-scale models of helium cluster dynamics in plasma-facing materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Cross slip of extended dislocations in face-centered cubic metals through phase-field modeling

Cross slip is a dislocation mechanism that significantly impacts the mechanical behavior of engineering alloys. Here, in this work, we advance a 3D phase-field dislocation dynamics (PFDD) mesoscale technique to simulate cross slip across a broad range of face-centered cubic (FCC) metals. The formulation incorporates elastic anisotropy, an FCC numerical grid, and a high-fidelity representation of the entire γ -surface from density functional theory for eight FCC metals and no adjustable parameters or rules. The relaxed core structures under zero stress for all metals are predicted to extend in plane. The analytical model for stacking fault width agrees well with the PFDD result under the assumption of elastic isotropy but overestimates it under elastic anisotropy, when the degree of anisotropy is large. The dynamic simulations are designed to elucidate the material parameters that influence the propensity for cross slip. Whether cross slip occurs under a non-Schmid stress or to bypass a hard obstacle, the critical stress to cross slip scales strongly with the anisotropic energy coefficient for a screw dislocation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Severe declines in hydraulic capacity and associated carbon starvation drive mortality in seawater exposed Sitka-spruce (Picea sitchensis) trees

Abstract Sea-level rise is causing widespread tree mortality of coastal forests, with large consequences on the Earth system as a result of these forests’ importance in carbon and nutrient export. The mechanisms of mortality under these conditions are, however, poorly tested. We used wood anatomy traits, wood δ 13 C, and tree radial growth to retrospectively assess the physiological process of seawater effects on whole tree xylem hydraulic capacity, gas exchange, and radial growth. During the latter stages of mortality (2018–2019), we directly measured metrics of water use and carbon metabolism across trees having crowns ranging from fully foliated to completely defoliated to investigate mortality processes at the sub-annual scale. Upon seawater exposure, soil salinity increased and allocation to hydraulic function declined, resulting in a dramatic reduction in water supply to the crown, increased crown-level water stress, and subsequent crown foliage loss. Simultaneously, leaf-level photosynthetic capacity declined steeply with increasing salinity. The combined loss of crown foliage area and photosynthetic rates per unit leaf area promoted carbon starvation, while no evidence of hydraulic failure was observed. These results elucidate mechanisms of coastal forest death under seawater exposure, enabling more accurate modeling in the future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The Ancient Salicoid Genome Duplication Event: A Platform for Reconstruction of De Novo Gene Evolution in Populus trichocarpa

Orphan genes are characteristic genomic features that have no detectable homology to genes in any other species and represent an important attribute of genome evolution as sources of novel genetic functions. Here, we identified 445 genes specific to Populus trichocarpa. Of these, we performed deeper reconstruction of 13 orphan genes to provide evidence of de novo gene evolution. Populus and its sister genera Salix are particularly well suited for the study of orphan gene evolution because of the Salicoid whole-genome duplication event which resulted in highly syntenic sister chromosomal segments across the Salicaceae. We leveraged this genomic feature to reconstruct de novo gene evolution from intergenera, interspecies, and intragenomic perspectives by comparing the syntenic regions within the P. trichocarpa reference, then P. deltoides, and finally Salix purpurea. Furthermore, we demonstrated that 86.5% of the putative orphan genes had evidence of transcription. Additionally, we also utilized the Populus genome-wide association mapping panel, a collection of 1,084 undomesticated P. trichocarpa genotypes to further determine putative regulatory networks of orphan genes using expression quantitative trait loci (eQTL) mapping. Functional enrichment of these eQTL subnetworks identified common biological themes associated with orphan genes such as response to stress and defense response. We also identify a putative cis-element for a de novo gene and leverage conserved synteny to describe evolution of a putative transcription factor binding site. Overall, 45% of orphan genes were captured in trans-eQTL networks.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Editorial: Plant-microbe omics

Omics-based studies have evolved over the past two decades to provide rich datasets from which deeper system-level understanding can be unraveled. Omics technologies (e.g., genomics, transcriptomics, proteomics, and metabolomics) can provide specific molecular insights into plant and microbe fitness, disease states, stress, cell signaling/cell-cell communication, and nutrient exchange, and these molecular observations can be correlated with system phenotype and function. The overall health of plant-microbe systems and their surrounding environments can be assessed, and a greater understanding of processes related to nutrient exchange/cycling, plant disease, and ecosystem homeostasis can now be achieved via the latest approaches to plant-microbe omics analysis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Compaction of crushed salt for safe containment – overview of the KOMPASS project

Abstract. In Germany, rock salt formations are possible host rock candidates for a repository for heat-emitting radioactive waste. The safety concept of a repository in salt bases on a multibarrier system consisting mainly of the geological barrier salt and geotechnical seals ensuring safe containment. Crushed salt will be used for backfilling of cavities and sealing measures in drifts and shafts due to its favourable properties and its easy availability (mined-off material). The creep of the rock salt leads to crushed salt compaction with time. Thereby, the crushed salts' porosity is reduced from the initial porosity of 30 %–40 % to a value comparable to the porosity of undisturbed rock salt (≤1 %). In such low porosity ranges, technical impermeability is assumed. The compaction behaviour of crushed salt is rather complex and involves several coupled THM processes (Kröhn et al., 2017; Hansen et al., 2014). It is influenced by internal properties like humidity and grain size distribution, as well as boundary conditions such as temperature, compaction rate or stress state. However, the current process understanding has some important gaps referring to the material behaviour, experimental database and numerical modelling. It needs to be extended and validated, especially in the low porosity range. The objective of the KOMPASS project was development of methods and strategies for the reduction of deficits in the prediction of crushed salt compaction leading to an improvement of the prognosis quality. Key results are as follows (KOMPASS Phase 1, 2020): selection of an easily available and permanently producible synthetic crushed salt mixture, acting as a reference material for generic investigations; development and proof of different techniques for producing pre-compacted samples for further investigations; establishment of a tool of microstructure investigation methods to demonstrate the comparability of grain structures of pre-compacted samples with in-situ compacted material for future investigations; execution of various laboratory experiments using pre-compacted samples, e.g. long-term creep tests which deliver reliable information about time- and stress-dependent compaction behaviour; development of a complex experimental investigation strategy to derive necessary model parameters considering individual functional dependencies. Its technical feasibility was successfully verified; benchmarking with various existing numerical models using datasets from three different triaxial long-term tests. The result was not entirely satisfactory; however, the number of influencing factors is small and further validation work has to be done. Overall, the KOMPASS project has made significant progress in the approaches to solving the outstanding question, building the basis for further investigations.

Friedenberg, Larissa↗
DOE OSTI2021

SAA1 is transcriptionally activated by STAT3 and accelerates renal interstitial fibrosis by inducing endoplasmic reticulum stress

Renal interstitial fibrosis (RIF) is the common irreversible pathway by which chronic kidney disease (CKD) progresses to the end stage. The transforming growth factor-β (TGF-β)/signal transducer and activator of transcription 3 (STAT3) signaling pathway is a common factor leading to inflammation-mediated RIF, but its downstream regulatory mechanism is still unclear. Bioinformatics analysis predicted that serum amyloid A protein 1 (SAA1) was one of the target genes for transcriptional activation of STAT3 signaling. As an acute phase reaction protein, SAA1 plays an important role in many inflammatory reactions, and research has suggested that SAA1 is significantly elevated in the serum of patients with CKD. In this research, multiple experiments were performed to investigate the role of SAA1 in the process of RIF. SAA1 was abnormally highly expressed in kidney tissue from individuals who underwent unilateral ureteral obstruction (UUO) and TGF-β-induced HK2 cells, and the abnormal expression was directly related to the transcriptional activation of STAT3. Additionally, SAA1 can directly target and bind valosin-containing protein (VCP)-interacting membrane selenoprotein (VIMP) to inhibit the function of the Derlin-1/VCP/VIMP complex, preventing the transportation and degradation of the misfolded protein, resulting in endoplasmic reticulum (ER) stress characterized by an increase in glucose-regulated protein 78 (GRP78) levels and ultimately promoting the occurrence and development of RIF.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Shear Stress Dependence of Force Networks in 3D Dense Suspensions

The geometric organization and force networks of 3D dense suspensions that exhibit both shear thinning and thickening have been examined as a function of varying strength of interparticle attractive interactions using lubrication flow discrete element simulations. Significant rearrangement of the geometric topology does not occur at either the local or global scale as these systems transition across the shear thinning and shear thickening regimes. In contrast massive rearrangements in the balance of attractive, lubrication, and contact forces are observed with interesting behavior of network growth and competition. In agreement with prior work, in shear thinning regions the attractive force is dominant, however as the shear thickening region is approached there is growth of lubrication forces. Lubrication forces oppose the attraction forces, but as viscosity continues to increase under increasing shear stress, the lubrication forces are dominated by contact forces that also resist attraction. Contact forces are the dominant interactions during shear thickening and are an order of magnitude higher than their values in the shear-thinning regime. At high attractive interaction strength, contact networks can form even under shear thinning conditions, however high shear stress is still required before contact networks become the driving mechanism of shear thickening. Analysis of the contact force network during shear thickening generally indicate a uniformly spreading network that rapidly forms across empty domains, however the growth patterns have structure that are significantly dependent upon the strength of interparticle interactions, indicating subtle variations in the mechanism of shear thickening.

Edens, Lance E.↗
DOE OSTI2022

Machine Learning of All Mycobacterium tuberculosis H37Rv RNA-seq Data Reveals a Structured Interplay between Metabolism, Stress Response, and Infection

Mycobacterium tuberculosis is one of the most consequential human bacterial pathogens, posing a serious challenge to 21st century medicine. A key feature of its pathogenicity is its ability to adapt its transcriptional response to environmental stresses through its transcriptional regulatory network (TRN). While many studies have sought to characterize specific portions of the M. tuberculosis TRN, and some studies have performed system-level analysis, few have been able to provide a network-based model of the TRN that also provides the relative shifts in transcriptional regulator activity triggered by changing environments. Here, we compiled a compendium of nearly 650 publicly available, high quality M. tuberculosis RNA-sequencing data sets and applied an unsupervised machine learning method to obtain a quantitative, top-down TRN. It consists of 80 independently modulated gene sets known as “iModulons,” 41 of which correspond to known regulons. These iModulons explain 61% of the variance in the organism’s transcriptional response. We show that iModulons (i) reveal the function of poorly characterized regulons, (ii) describe the transcriptional shifts that occur during environmental changes such as shifting carbon sources, oxidative stress, and infection events, and (iii) identify intrinsic clusters of regulons that link several important metabolic systems, including lipid, cholesterol, and sulfur metabolism. This transcriptome-wide analysis of the M. tuberculosis TRN informs future research on effective ways to study and manipulate its transcriptional regulation and presents a knowledge-enhanced database of all published high-quality RNA-seq data for this organism to date.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Insights into dual-functional modification for water stability enhancement of mesoporous zirconium metal–organic frameworks

The stability of metal–organic frameworks (MOFs) in water affects their ability to function as chemical catalysts, their capacity as adsorbents for separations in water vapor presence, and their usefulness as recyclable water harvesters. Here, we have examined water stability of four node-modified variants of the mesoporous MOF, NU-1000, namely formate-, Acac-, TFacac-, and Facac-NU-1000, comparing these with node-accessible NU-1000. These NU-1000 variants present ligands grafted to NU-1000's hexa-Zr(IV)-oxy nodes by displacing terminal aqua and hydroxo ligands. Facac-NU-1000, containing the most hydrophobic ligands, showed the greatest water stability, being able to undergo at least 20 water adsorption/desorption cycles without loss of water uptake capacity. Computational studies revealed dual salutary functions of installed Facac ligands: (1) enhancement of framework mechanical stability due to electrostatic interactions; and (2) transformation and shielding of the otherwise highly hydrophilic nodes from H-bonding interactions with free water, presumably leading to weaker channel-stressing capillary forces during water evacuation – consistent with trends in free energies of dehydration across the NU-1000 variants. Water harvesting and hydrolysis of chemical warfare agent simulants were examined to gauge the functional consequences of modification and mechanical stabilization of NU-1000 by Facac ligands. The studies revealed a harvesting capacity of ~1.1 L of water vapor per gram of Facac-NU-1000 per sorption cycle. They also revealed retention of catalytic MOF activity following 20 water uptake and release cycles. This study provides insights into the basis for node-ligand-engendered stabilization of wide-channel MOFs against collapse during water removal.

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

Diverse fates of ancient horizontal gene transfers in extremophilic red algae

Abstract Horizontal genetic transfer (HGT) is a common phenomenon in eukaryotic genomes. However, the mechanisms by which HGT‐derived genes persist and integrate into other pathways remain unclear. This topic is of significant interest because, over time, the stressors that initially favoured the fixation of HGT may diminish or disappear. Despite this, the foreign genes may continue to exist if they become part of a broader stress response or other pathways. The conventional model suggests that the acquisition of HGT equates to adaptation. However, this model may evolve into more complex interactions between gene products, a concept we refer to as the ‘Integrated HGT Model’ (IHM). To explore this concept further, we studied specialized HGT‐derived genes that encode heavy metal detoxification functions. The recruitment of these genes into other pathways could provide clear examples of IHM. In our study, we exposed two anciently diverged species of polyextremophilic red algae from the Galdieria genus to arsenic and mercury stress in laboratory cultures. We then analysed the transcriptome data using differential and coexpression analysis. Our findings revealed that mercury detoxification follows a ‘one gene‐one function’ model, resulting in an indivisible response. In contrast, the ars H gene in the arsenite response pathway demonstrated a complex pattern of duplication, divergence and potential neofunctionalization, consistent with the IHM. Our research sheds light on the fate and integration of ancient HGTs, providing a novel perspective on the ecology of extremophiles.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

A Novel Mode of Photoprotection Mediated by a Cysteine Residue in the Chlorophyll Protein IsiA

Oxygenic photosynthetic organisms have evolved a multitude of mechanisms for protection against high-light stress. IsiA, a chlorophyll a-binding cyanobacterial protein, serves as an accessory antenna complex for photosystem I. Intriguingly, IsiA can also function as an independent pigment protein complex in the thylakoid membrane and facilitate the dissipation of excess energy, providing photoprotection. The molecular basis of the IsiA-mediated excitation quenching mechanism remains poorly understood. In this study, we demonstrate that IsiA uses a novel cysteine-mediated process to quench excitation energy. The single cysteine in IsiA in the cyanobacterium Synechocystis sp. strain PCC 6803 was converted to a valine. Ultrafast fluorescence spectroscopic analysis showed that this single change abolishes the excitation energy quenching ability of IsiA, thus providing direct evidence of the crucial role of this cysteine residue in energy dissipation from excited chlorophylls. Under stress conditions, the mutant cells exhibited enhanced light sensitivity, indicating that the cysteine-mediated quenching process is critically important for photoprotection.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Tailoring a Global Iron Regulon to a Uropathogen

Pathogenicity islands and plasmids bear genes for pathogenesis of various Escherichia coli pathotypes. Although there is a basic understanding of the contribution of these virulence factors to disease, less is known about variation in regulatory networks in determining disease phenotypes. Here, we dissected a regulatory network directed by the conserved iron homeostasis regulator, ferric uptake regulator (Fur), in uropathogenic E. coli (UPEC) strain CFT073. Comparing anaerobic genome-scale Fur DNA binding with Fur-dependent transcript expression and protein levels of the uropathogen to that of commensal E. coli K-12 strain MG1655 showed that the Fur regulon of the core genome is conserved but also includes genes within the pathogenicity/genetic islands. Unexpectedly, regulons indicative of amino acid limitation and the general stress response were also indirectly activated in the uropathogen fur mutant, suggesting that induction of the Fur regulon increases amino acid demand. Using RpoS levels as a proxy, addition of amino acids mitigated the stress. In addition, iron chelation increased RpoS to the same levels as in the fur mutant. The increased amino acid demand of the fur mutant or iron chelated cells was exacerbated by aerobic conditions, which could be partly explained by the O 2 -dependent synthesis of the siderophore aerobactin, encoded by an operon within a pathogenicity island. Taken together, these data suggest that in the iron-poor environment of the urinary tract, amino acid availability could play a role in the proliferation of this uropathogen, particularly if there is sufficient O 2 to produce aerobactin. IMPORTANCE Host iron restriction is a common mechanism for limiting the growth of pathogens. We compared the regulatory network controlled by Fur in uropathogenic E. coli (UPEC) to that of nonpathogenic E. coli K-12 to uncover strategies that pathogenic bacteria use to overcome iron limitation. Although iron homeostasis functions were regulated by Fur in the uropathogen as expected, a surprising finding was the activation of the stringent and general stress responses in the uropathogen fur mutant, which was rescued by amino acid addition. This coordinated global response could be important in controlling growth and survival under nutrient-limiting conditions and during transitions from the nutrient-rich environment of the lower gastrointestinal (GI) tract to the more restrictive environment of the urinary tract. The coupling of the response of iron limitation to increased demand for amino acids could be a critical attribute that sets UPEC apart from other E. coli pathotypes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Multiaxial Plastic Deformation of Zircaloy-4 Nuclear Fuel Cladding Tubes

Here, this work is motivated by the desire to devise an internal pressure test that can mimic a displacement-controlled loading scenario and demonstrate how to apply the multiaxial stress and strain data from the test to develop an elastic/plastic constitutive model for a thin-walled tubular component. This is achieved by conducting simultaneous measurements of tangential and axial strain during the pressure test and integrating these strain measures into a feedback loop with the pressure controller. It is shown how data from such a test can be used to develop a large mechanical property data set relevant to biaxial loading conditions. The data obtained have high confidence evidenced by their low variability and alignment with other literature studies. Additionally, data from these internal pressure tests combined with full-tube axial tensile tests allow for the derivation of the Hill anisotropic yield function. The developed Hill yield function is validated by comparing the plastic strain ratios from the full-tube tension tests and by comparing the predicted yield stress in the tangential direction with measured values from ring tension tests in a previous study.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Vegetation responses to air dryness amplify future land surface warming

Temperature exerts a first-order control on vegetation photosynthesis and transpiration. Yet most studies investigating temperature impacts on plants rely on near-surface air temperature, rather than canopy temperature—the temperature plants actually experience. Because canopy temperature directly regulates ecosystem function, it provides a more accurate measure of vegetation–climate interactions. Combining Earth System Model (ESM) simulations and satellite observations in a dual emergent constraint, here we show that canopy temperature is projected to increase substantially more than air temperature (~0.11-degrees more or a 16% increase in their difference) over the 21st century. The ESM ensemble median fails to capture these stronger increases in the majority of vegetated regions. We find that the largest projected increases in the difference between canopy and air temperature are predicted to occur in regions where elevating moisture stress—particularly rising vapor pressure deficit—increasingly constrains vegetation growth and transpiration. This implies that future warming will impose stronger constraints on plant function than currently estimated. Relying on air temperature alone will therefore lead to systematic underestimation of temperature effects on photosynthesis, vegetation growth, and the land carbon sink. Accurate representation of canopy temperature in ESMs is thus essential to improve projections of ecosystem responses and feedbacks to climate change.

Green, Julia K↗
DOE OSTI2026

Cross-slip and easy-glide CRSS in titanium: Theoretical predictions and in-situ TEM measurements

This study investigates the mechanics of prismatic and first-order pyramidal $\langle$a$\rangle$ slip in titanium (Ti), elucidating the physics of easy-glide and cross-slip through a combination of theory and experiments. Screw-character prismatic (Pr) dislocations in Ti are of particular interest because their complex cores can be stable or unstable, leading to activation by either cross-slip or planar glide. To investigate these mechanisms, site-specific micro tensile samples were prepared using focused ion beam (FIB) milling and mounted on a push-to-pull (PTP) device for in-situ transmission electron microscopy (TEM) tensile testing. The in-situ experiments provide direct observations of the onset of dislocation motion and the precise determination of the critical resolved shear stress (CRSS) for the activated mechanisms, and its evolution with load cycling. A comprehensive theory has been developed to predict the CRSS values for easy glide, cross-slip, and multiplication of dislocations. Predicted critical stresses for pyramidal (π)-to-Pr and reverse cross-slip agree closely with the experimental measurements. The latter cross-slip stress is a factor of two higher than that of unobstructed planar slip. The model accounts for overlapping dislocation cores and employs a Wigner-Seitz based cell to evaluate misfit energies. By combining ab initio density functional theory (DFT) with anisotropic elasticity, the framework identifies minimum energy pathways for dislocation glide, which can be intermittent and zig-zag. A simplified expression utilizing (π) and Pr Schmid factor ratios is proposed for critical stress corresponding to (π)-to-Pr cross-slip transition. The results are strongly dependent on crystal orientation, underscoring non-Schmid behavior. Overall, this study explores key critical stress parameters essential for informing higher-scale simulations of plasticity in Ti.

Critical stress↗
DOE OSTI2020

Guided search for desired functional responses via Bayesian optimization of generative model: Hysteresis loop shape engineering in ferroelectrics

Advances in theoretical modeling across multiple disciplines have yielded generative models capable of high veracity in predicting macroscopic functional responses of materials emerging as a result of complex non-local interactions. Correspondingly, of interest is the inverse problem of finding the model parameter that will yield desired macroscopic responses, such as stress–strain curves, ferroelectric hysteresis loops, etc. Here, we suggest and implement Gaussian process based methods that allow to effectively sample the degenerate parameter space of a complex non-local model to output regions of parameter space which yield desired functionalities. In this work, we discuss the specific adaptation of the acquisition function and sampling function to make the process efficient and balance the efficient exploration of parameter space for multiple possible minima and exploitation to densely sample the regions of interest where target behaviors are optimized. Additionally, this approach is illustrated via the hysteresis loop engineering in ferroelectric materials but can be adapted to other functionalities and generative models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗