Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stability assessment”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Physics basis for the reference flat-top plasma scenario in the ST–E1 fusion power plant

As part of the U.S. Department of Energy’s Milestone-Based Fusion Energy Development Program, Tokamak Energy has completed the pre-concept design of the ST–E1 fusion power plant. ST–E1 is envisaged to operate in two phases: a pilot plant phase, targeting sustained net power production of 300 - 500 MWe for a duration >1 hr, followed by a commercial power plant phase targeting steady-state operations and a normalised overnight capital cost of ⩽12 000 $\$$/kWe. The design process adopted was highly iterative, integrating all major plant systems and progressing in a phased fidelity approach. At the pre-conceptual stage, the emphasis has been on exploring the design space, identifying the main system-level trade-offs, and making the key decisions that define the overall plant concept, rather than optimising a single operating point. This paper, part of a focused collection detailing the ST–E1 pre-concept design, addresses the development of a series of reference flat-top plasma operating points for the pilot plant phase. A modelling workflow was established to develop and assess candidate plasma design points and explore key dependencies. The workflow includes integrated core plasma modelling, magnetohydrodynamic (MHD) stability assessment, equilibrium generation, scrape-off-layer and exhaust modelling, heating & current drive design and optimisation, and turbulent transport modelling. Using this framework, the impact of several key parameters on the flat-top operating space was investigated, including the density limit, core radiation fraction and divertor power loading, level of external heating and curent drive power and assumed pedestal characteristics. The MHD stability, controllability and micro-stability characteristics of these plasmas were also analysed. These investigations informed the definition of a set of fully non-inductive, flat-top reference operating points that satisfy the high-level ST–E1 mission, including a low and high density case, a case that is stable to resistive wall modes and a case with reduced divertor power loading.

ST–E1↗
DOE OSTI2023

Deep Sigma Point Processes-Assisted Chance-Constrained Power System Transient Stability Preventive Control

Here this paper proposes a deep sigma point processes (DSPP)-assisted chance-constrained power system transient stability preventive control method to deal with uncertain renewable energy and loads-induced stability risk. The traditional transient stability-constrained preventive control is reformulated as a chance-constrained optimization problem. To deal with the computational bottleneck of the time-domain simulation-based probabilistic transient stability assessment, the DSPP is developed. DSPP is a parametric Bayesian approach that allows us to predict system transient stability with high computational efficiency while accurately quantifying the confidence intervals of the predictions that can be used to inform system instability risk. To this end, with a given preset confidence probability, we embed DSPP into the primal dual interior point method to help solve the chance-constrained preventive control problem, where the corresponding Jacobian and Hessian matrices are derived. Comparison results with other existing methods show that the proposed method can significantly speed up preventive control while maintaining high accuracy and convergence

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Operationalizing Synchrophasors for Enhanced Grid Reliability and Asset Utilization

The report discusses the development of production-grade synchrophasor software applications. The objective is to improve system reliability and asset utilization for offline as well as real-time considerations, as per the Department of Energy (DOE) guidelines specified by Funding Opportunity Announcement (FOA) DE-FOA-0001492. On the one hand, we have developed grid management applications on how to transform synchrophasor-based monitoring into operator guidance, and commercial grade tools addressing industry needs such as NERC Reliability Standards (MOD 26, MOD 27 and MOD 33). On the other hand, with the wide industrial adoption in mind, the developed applications have been fully discussed with utility partners and installed in the utility partners site including 6 field tests. The Model Validation and Calibration software module has been integrated into GE’s PhasorAnalyticsTM. The first commercial contract has been signed even before the close date of this project. The Angle-based Grid Management (AGM), Fast Voltage Security Assessment (FastVSA), Enhanced Islanding Monitoring(EIM) and Enhanced Disturbance Management (EDM) have been integrated into GE’s Grid Stability Assessment line of applications in e-terra Suite of EMS product.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Limitations of Thermal Stability Analysis via In-Situ TEM/Heating Experiments

This work highlights some limitations of thermal stability analysis via in-situ transmission electron microscopy (TEM)-annealing experiments on ultrafine and nanocrystalline materials. We provide two examples, one on nanocrystalline pure copper and one on nanocrystalline HT-9 steel, where in-situ TEM-annealing experiments are compared to bulk material annealing experiments. The in-situ TEM and bulk annealing experiments demonstrated different results on pure copper but similar output in the HT-9 steel. The work entails discussion of the results based on literature theoretical concepts, and expound on the inevitability of comparing in-situ TEM annealing experimental results to bulk annealing when used for material thermal stability assessment.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Physiochemical Machine Learning Models Predict Operational Lifetimes of CH3NH3PbI3 Perovskite Solar Cells

Halide perovskites are promising photovoltaic (PV) materials with the potential to lower the cost of electricity and greatly expand the penetration of PV if they can demonstrate long-term stability under illumination in the presence of moisture and oxygen. The solar cell service lifetime as quantified by the T80 (the time required for the power conversion efficiency to drop to 80% of its starting value) is a useful metric to assess stability. The T80 for utility, commercial, or residential PV systems needs to be several decades in order to yield low-cost electricity, and thus it is not practical to directly measure the T80. It would be useful if T80 could be predicted from the initial dynamics of a solar cell’s performance, but until now no models have been developed to forecast T80. In this work, we report the development of machine learning models to predict T80 of ITO/NiOx/CH3NH3PbI3/C60/BCP/Ag solar cells operating at maximum power point under 1-sun equivalent photon flux in air at varying temperatures and relative humidities. Efficiency losses are driven by short-circuit current and fill factor, indicating that chemical decomposition of the perovskite is a major contributor to degradation. Spatial patterns evident from in situ dark field optical microscopy suggest that the electric field gradient at device edges plays a significant role in perovskite decomposition, along with photochemical reactions with O2 and H2O. Models are trained using a menu of features from three distinct categories: (i) features based on measurements of the initial rates of change of device parameters, (ii) features based on the ambient conditions during operation (temperature, & partial pressure of H2O), and (iii) features based on underlying physics and chemistry. We show that a theory-based physiochemical feature derived from a model of the chemical reaction kinetics of the rate of degradation of the CH3NH3PbI3 is particularly valuable for prediction. This physiochemical feature was selected as the first or second most dominant feature in the best performing models. With a dataset consisting of 45 accelerated degradation experiments with T80 that range over a factor of 30, the model predicts T80 with an accuracy of about 40% (|predicted T80 - observed T80| / observed T80) on samples not used in training. This hybrid ML approach should be effective when applied to other compositions, device architectures, and advanced packaging schemes.

14 SOLAR ENERGY↗
DOE OSTI2022

Physiochemical Machine Learning Models Predict Operational Lifetimes of CH3NH3PbI3 Perovskite Solar Cells

Halide perovskites are promising photovoltaic (PV) materials with the potential to lower the cost of electricity and greatly expand the penetration of PV if they can demonstrate long-term stability under illumination in the presence of moisture and oxygen. The solar cell service lifetime as quantified by the T80 (the time required for the power conversion efficiency to drop to 80% of its starting value) is a useful metric to assess stability. The T80 for utility, commercial, or residential PV systems needs to be several decades in order to yield low-cost electricity, and thus it is not practical to directly measure the T80. It would be useful if T80 could be predicted from the initial dynamics of a solar cell’s performance, but until now no models have been developed to forecast T80. In this work, we report the development of machine learning models to predict T80 of ITO/NiOx/CH3NH3PbI3/C60/BCP/Ag solar cells operating at maximum power point under 1-sun equivalent photon flux in air at varying temperatures and relative humidities. Efficiency losses are driven by short-circuit current and fill factor, indicating that chemical decomposition of the perovskite is a major contributor to degradation. Spatial patterns evident from in situ dark field optical microscopy suggest that the electric field gradient at device edges plays a significant role in perovskite decomposition, along with photochemical reactions with O2 and H2O. Models are trained using a menu of features from three distinct categories: (i) features based on measurements of the initial rates of change of device parameters, (ii) features based on the ambient conditions during operation (temperature, & partial pressure of H2O), and (iii) features based on underlying physics and chemistry. We show that a theory-based physiochemical feature derived from a model of the chemical reaction kinetics of the rate of degradation of the CH3NH3PbI3 is particularly valuable for prediction. This physiochemical feature was selected as the first or second most dominant feature in the best performing models. With a dataset consisting of 45 accelerated degradation experiments with T80 that range over a factor of 30, the model predicts T80 with an accuracy of about 40% (|predicted T80 - observed T80| / observed T80) on samples not used in training. This hybrid ML approach should be effective when applied to other compositions, device architectures, and advanced packaging schemes.

14 SOLAR ENERGY↗
DOE OSTI2026

Small-Signal Stability of Grid-Forming Converters Under Fault Conditions

Threshold virtual impedance (TVI)-based current limiting for grid-forming converters (GFMs) has gained great interest due to its ability to maintain voltage source behaviour during faults. However, sequence component extraction (SCE) and negative-sequence control (NSC) are often overlooked in small-signal stability assessments during faults. This paper develops small-signal sequence impedance models for GFMs under four well-known SCE methods based on TVI current limiting control during symmetrical fault conditions. Using the developed impedance models, the impacts of SCE and NSC, and the voltage and current control loop bandwidths, on system stability during faults are investigated. Additionally, since negative-sequence TVI (TVI-) is typically added along with its positive-sequence counterpart, which is often inductive, inductive and resistive TVI- are examined. The findings suggest that a higher voltage or current control loop bandwidth has a negative impact on system stability, while SCE and NSC largely reduce the stable range for voltage and current control loop bandwidth during faults, and that the severity of such impacts is determined by the particular SCE method. Furthermore, it is observed that inductive TVI- significantly degrades system stability, while resistive TVI- can enhance stability when suitable SCE methods are appropriately selected and designed. Matlab/Simulink electromagnetic transient simulations validate these analytical results.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Small-Signal Stability Impacts of Load and Network Dynamics on Grid-Forming Inverters

This paper presents several stability analyses for grid-forming inverters and synchronous generators considering the dynamics of transmission lines and different load models. Load models are usually of secondary importance compared to generation source models, but as the results show, they play a crucial role in stability studies with the introduction of inverter-based resources. Given inverter control time scales, the implications of considering or neglecting electromagnetic transients of the network are very relevant in the stability assessments. In this paper, we perform eigenvalue analyses for inverter-based resources and synchronous machines connected to a load and explore the effects of multiples models under different network representations. We explore maximum loadability of inverter-based resources and synchronous machines, while analyzing the effects of load and network dynamic models on small-signal stability. The results show that the network representation plays a fundamental role in the stability of the system of different load models. The resulting stability regions are significantly different depending on the source and load model considered.

grid forming inverters↗
DOE OSTI2020

First-principles prediction of incipient order in arbitrary high-entropy alloys: exemplified in Ti 0.25 CrFeNiAl

High-temperature disordered multi-component alloys, including high-entropy alloys, experience either segregation or partially-ordered phases to reach low-temperature phases. For Ti 0.25 CrFeNiAl x (0 ≤ x ≤ 1), experiments suggest a partially-ordered B2 phase, whereas CALculation of PHAse Diagrams (CALPHAD) predicts a region of L21+B2 coexistence. We employ first-principles KKR-CPA to assess stability of phases with arbitrary order and a KKR-CPA linear-response theory to predict atomic short-range order (SRO) in the disorder phase that reveals the competing long-range ordered (LRO) phases in a given Bravais lattice. The favorable SRO provides a specific concentration-waves (site occupation probabilities and partially-ordered unit cells) and estimated energy gains that can then be assess directly by KKR-CPA formation enthalpies. Our results are in good agreement with experiments and CALPHAD in Al-poor regions (x ≤ 0.75) and with CALPHAD in Al-rich region (1≥ x >0.75). Our first-principles KKR-CPA and SRO-based concentration-wave analysis is shown to be a powerful and fast method to assess competing LRO phases in complex solid-solution alloys, and our results suggests more careful experiments in Al-rich region are needed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Critical role of oxide anchoring groups in organic electron transport layers for perovskite solar cell stability

Charge transport layers (CTLs) play a critical role in the performance and long-term stability of perovskite solar cells (PSCs) by facilitating efficient charge extraction and providing interfacial stabilization. For organic CTLs, different functional groups have been designed to modulate the interactions with the substrate. Anchoring groups are used to establish strong adhesion to transparent conductive oxides (TCOs) but are not routinely incorporated in organic CTL design. In this study, we investigate the influence of anchoring groups in organic electron transport layers (ETLs) on PSC performance and stability by synthesizing and evaluating two naphthalene diimide (NDI)-based ETLs: NDI-(PhPA) 2 , functionalized with polar phenyl phosphonic acid (PhPA) groups, and NDI-(PhBr) 2 , modified with nonpolar bromophenyl (PhBr) groups. X-ray photoelectron spectroscopy confirms that PhPA functionalization leads to NDIs anchored onto fluorine-doped tin oxide (FTO) following chemical bath deposition, while NDI-(PhBr) 2 shows no detectable surface presence after deposition. We show that NDI-(PhPA) 2 deposits as a thin film on FTO, exhibiting strong adhesion, high solvent resistance, and enhanced light and thermal stability, achieving a maximum device efficiency of 14.3%. Stability assessments conducted under continuous illumination at 25 °C for 200 hours and at 65 °C for 100 hours demonstrate performance on par with TiO 2 -based ETLs. In contrast, NDI-(PhBr) 2 , which lacks strong anchoring interactions, fails to adhere to the FTO surface, resulting in a rapid decrease in efficiency over a few hours, comparable to devices without an ETL. These findings highlight the importance of anchoring groups in organic CTLs for achieving long-term PSC stability under operational conditions.

Kim, Sanggyun [Georgia Institute of Technology, At↗
DOE OSTI2021

A Multi-modal Approach to Understanding Degradation of Organic Photovoltaic Materials

State-of-the-art organic photovoltaic (OPV) materials are composed of complex, chemically diverse polymeric and molecular structures that form highly intricate solid-state interactions, collectively yielding exceptional tunability in performance and aesthetics. These properties are especially attractive for semitransparent power-generating windows or shades in living environments, greenhouses, or other architectural integrations. However, before such a future is realized, a broader and deeper understanding of property stability must be acquired. Stability during operating and environmental conditions is critical, namely, material color steadfastness, optoelectronic performance retention, morphological rigidity, and chemical robustness. To date, no single investigation encompasses all four distinct, yet interconnected, metrics. Here, we present a multimodal strategy that captures a dynamic and interconnected evolution of each property during the course of an accelerated photobleaching experiment. We demonstrate this approach across relevant length scales (from molecular to visual macroscale) using X-ray photoelectron spectroscopy, grazing-incidence X-ray scattering, microwave conductivity, and time-dependent photobleaching spectroscopies for two high-performance semitransparent OPV blends—PDPP4T:PC 60 BM and PDPP4T:IEICO-4F, with comparisons to the stabilities of the individual components. We present direct evidence that specific molecular acceptor (fullerene vs nonfullerene) designs and the resulting donor–acceptor interactions lead to distinctly different mechanistic routes that ultimately arrive at what is termed “OPV degradation.” We directly observe a chemical oxidation of the cyano endcaps of the IEICO-4F that coincides with a morphological change and large loss in photoconductivity while the fullerene acceptor-containing blend demonstrates a significantly greater fraction of oxygen uptake but retains 55% of the photoconductivity. Finally, this experimental roadmap provides meaningful guidance for future high-throughput, multimodal studies, benchmarking the sensitivity of the different analytical techniques for assessing stability in printable active layers, independent of complete device architectures.

14 SOLAR ENERGY↗
DOE OSTI2022

A Laplace-Domain Circuit Model for Fault and Stability Analysis Considering Unbalanced Topology

For systems subject to unbalanced faults, analytical model building for stability assessment is a challenging task. This letter presents a straightforward modeling approach. A generalized dynamic circuit representation is achieved by use of the Laplacian transform variable s . Here, we translate the voltage and current relationship at the fault location into the relationship of three subsystems. The final circuit model is an interconnected sequence network with impedances in the Laplace domain. This circuit can be directly converted from a steady-state sequence network. This modeling procedure is illustrated by an example case of an induction motor served by a grid through a series compensated line. Electromagnetic transient simulation results demonstrate that sub-synchronous oscillations can be mitigated when a single-line to ground fault is applied at the motor terminal. Stability analysis results based on the dynamic circuit corroborate the simulation results. What's more, the derived circuit effortlessly reveals why unbalance can enhance stability.

42 ENGINEERING↗
DOE OSTI2025

Large-signal Stability Analysis of Grid-forming Inverters with Equivalent-circuit Models

Here, this paper proposes an energy function-based direct method for large-signal stability assessment of grid-forming (GFM) inverters leveraging an equivalent-circuit representation of all involved control- and physical-layer dynamics. Three different primary controls, a standard inner-current outer-voltage cascaded-control architecture, output LCL filter, and reference-current saturation limiting are featured in the modeling and analysis framework. A composite energy function for the GFM inverter is obtained by summing up individual energy contributions gleaned from the circuit representation. The approach can readily be generalized to different primary controls, output-filter arrangements, and current limiters since it is based on a circuit-theoretic foundation. Numerical simulations validate the efficacy of the approach in estimating the critical clearing time following a large-signal disturbance.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Leveraging National Laboratory Assets to Address Stability Challenges due to Declining Grid Inertia Using Geographically Distributed Electrical–Thermal Co-Emulation

Due to increased penetration of low-inertia resources into the electric grid, challenges are increasing for maintaining wide-area system stability. Grid stability assessment requires a faithful representation of the multiple-physics interaction at the system level, and timescales of interaction varying in orders of magnitude, from microseconds to seconds to several minutes. Along with the simulation-based techniques, hardware-in-the-loop (HIL), controller HIL, and power HIL techniques have been developed to better understand the emergent behavior of the system with emerging technologies. US National Laboratories have played a vital role in research and development to understand the behavior of individual technologies and devices integrated to the electric grid. Each national laboratory forwards a technological and strategic initiative tied core and enabling capabilities. Due to strategic, efficiency, and economic reasons, not all the labs have assets to conduct research on all technologies concomitantly, so it becomes crucial to integrate the labs across geographies to understand the interplay of different technologies together at the system level. This approach avoids duplication of the assets at different lab facilities and helps understand the integrated system behavior of various technologies representative of actual grid conditions by connecting multiple national labs. This article talks about techniques of connecting three national laboratories to enable co-emulation of electrical–mechanical–thermal characteristics of devices and systems. Such an approach can be used to understand the dynamic and transient interaction of multi-physics in a system level, at-scale emulation using real-time simulation tools and techniques.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Borehole breakout modeling in arkose and granite rocks

The existence of a deep borehole in the Earth’s crust disturbs the local stresses and creates a stress concentration that may result in breakout and damage to the borehole. Maintaining wellbore integrity mitigates environmental impacts such as groundwater contamination, gas leakage to the atmosphere, and fluid spills and seepage at the surface. In this paper, the stability of deep boreholes (5 km) is examined by laboratory experiments and numerical models in the context of nuclear waste disposal in Israel. Two rock types in southern Israel are considered: the crystalline basement (granite) and the Zenifim Formation (arkose). A series of room-temperature triaxial rock deformation experiments were conducted at different confining pressures. This mechanical characterization was then used to parameterize the elastic properties and damage behavior of the rocks. This facilitated modeling the stability of the deep boreholes by two different formulations of damage rheology: a dynamic-oriented formulation used to model deformation immediately after the creation of the open hole and a quasi-static formulation used to model longer stress corrosion regime. The calibrated modeling results indicate greater stability with Zenifim arkose than the crystalline granite for deep borehole conditions despite the granite having a greater triaxial compressive strength. Dissipation associated with dilation and porous compaction in the arkose during deformation plays a significant stabilizing role in the borehole compared to crystalline rocks. These results suggest that common strength-based borehole stability assessment may lead to inaccurate predictions. Three-dimensional modeling of bottom-hole stress conditions and the effects of transient borehole geometry show conventional two-dimensional analysis may not be conservative when predicting borehole damage.

58 GEOSCIENCES↗
DOE OSTI2023

High-Precision Calculation of Nanoparticle (Nanocrystal) Potentials of Mean Force and Internal Energies

Thermodynamic stability assessment of nanocrystal systems requires precise free energy calculations. This study highlights the importance of meticulous control over various factors, including the thermostat, time step, potential cutoff, initial configuration, sampling method, and overall simulation duration. Free energy computations in dry (solvent-free) systems are on the order of several hundred $k_BT$ but can be obtained with consistent accuracy. However, calculation of internal energies becomes challenging, as they are typically much larger in magnitude than free energies and exhibit significant noise and reduced reliability. To address this limitation, we propose a new internal energy estimate that drastically reduces the noise. Here, we also present formulas that enable the optimization of the parameters of the harmonic bias potential for optimal convergence. Finally, we discuss the implications of these findings for the computation of free energies in nanocrystal clusters and superlattices.

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

Paramagnetic Mn 8 Fe 4 -co-Polystyrene Nanobeads as a Potential T 1 –T 2 Multimodal Magnetic Resonance Imaging Contrast Agent with In Vivo Studies

In developing a cluster-nanocarrier design, as a magnetic resonance imaging contrast agent, we have investigated the enhanced relaxivity of a manganese and iron-oxo cluster grafted within a porous polystyrene nanobead with increased relaxivity due to a higher surface area. The synthesis of the cluster-nanocarrier for the cluster Mn 8 Fe 4 O 12 (O 2 CC 6 H 4 CH=CH 2 ) 16 (H 2 O) 4 , cross-linked with polystyrene (the nano-carrier), under miniemulsion conditions is described. By including a branched hydrophobe, iso-octane, the resulting nanobeads are porous and similar to 70 nm in diameter. The increased surface area of the nanobeads compared to nonporous nanobeads leads to an enhancement in relaxivity; r 1 increases from 3.8 to 5.2 ± 0.1 mM -1 s -1 , and r 2 increases from 11.9 to 50.1 ± 4.8 mM -1 s -1 , at 9.4 teslas, strengthening the potential for T 1 and T 2 imaging. Several metrics were used to assess stability, and the porosity produced no reduction in metal stability. Synchrotron X-ray fluorescence microscopy was used to demonstrate that the nanobeads remain intact in vivo . In depth, physicochemical characteristics were determined, including extensive pharmacokinetics, in vivo imaging, and systemic biodistribution analysis.

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

Candidate ferroelectrics via ab initio high-throughput screening of polar materials

Ferroelectrics are a class of polar and switchable functional materials with diverse applications, from microelectronics to energy conversion. Computational searches for new ferroelectric materials have been constrained by accurate prediction of the polarization and switchability with electric field, properties that, in principle, require a comparison with a nonpolar phase whose atomic-scale unit cell is continuously deformable from the polar ground state. For most polar materials, such a higher-symmetry nonpolar phase does not exist or is unknown. Here, we introduce a general high-throughput workflow that screens polar materials as potential ferroelectrics. We demonstrate our workflow on 1978 polar structures in the Materials Project database, for which we automatically generate a nonpolar reference structure using pseudosymmetries, and then compute the polarization difference and energy barrier between polar and nonpolar phases, comparing the predicted values to known ferroelectrics. Focusing on a subset of 182 potential ferroelectrics, we implement a systematic ranking strategy that prioritizes candidates with large polarization and small polar-nonpolar energy differences. To assess stability and synthesizability, we combine information including the computed formation energy above the convex hull, the Inorganic Crystal Structure Database id number, a previously reported machine learning-based synthesizability score, and ab initio phonon band structures. To distinguish between previously reported ferroelectrics, materials known for alternative applications, and lesser-known materials, we combine this ranking with a survey of the existing literature on these candidates through Google Scholar and Scopus databases, revealing ~130 promising materials uninvestigated as ferroelectric. Our workflow and large-scale high-throughput screening lays the groundwork for the discovery of novel ferroelectrics, revealing numerous candidates materials for future experimental and theoretical endeavors.

36 MATERIALS SCIENCE↗