SEARCH · Engineering Papers

Results for “3D systems”

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

SourceResearch connectionFrom yearThrough year

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

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2024

A comparative study on deep learning models for condition monitoring of advanced reactor piping systems

Advanced nuclear reactors offer innovative applications due to their portability, reliability, resiliency, and high capacity factors. To operate them on a wider scale, reducing maintenance life-cycle costs while ensuring their integrity is essential. Autonomous operations in advanced nuclear reactors using augmented Digital Twin (DT) technology can serve as a cost-effective solution by increasing awareness about the system’s health. A key component of nuclear DT frameworks is the condition monitoring of safety systems, such as piping-equipment systems, which involves acquiring and monitoring the plant’s sensor data. Here, this research proposes a condition monitoring methodology utilizing deep learning algorithms, such as multilayer perceptions (MLP) and convolutional neural networks (CNNs), to detect degradation and its severity in nuclear piping-equipment systems. Sensor signals are processed to obtain the power spectral density and the Short-Time Fourier transform, and feature extraction methodologies are proposed to develop degradation-sensitive data repositories. The performance of MLP, one-dimensional (1D) CNN, and 2D CNN within the proposed condition monitoring framework is compared using a finite element model of a 3D piping system subjected to seismic loads as the application case study. Various approaches, such as dropout, k-Fold validation, regularization, and early stopping of training the network, are investigated to avoid overfitting the models to the input sensor data. The predictive capability and computational capacity of the deep learning algorithms are also compared to detect degradation in the Z-pipe system of the Experimental Breeder Reactor II (EBRII). The Z-pipe system is subjected to harmonic excitations that represent normal operating loads, such as pump-induced vibrations. The findings of the study indicate that the proposed artificial intelligence (AI)-driven condition monitoring framework demonstrates superior prediction accuracies with a 2D CNN, whereas the MLP exhibits higher computational efficiency.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS2017

Imaging Systems for Size Measurements of Debrisat Fragments

The overall objective of the DebriSat project is to provide data to update existing standard spacecraft breakup models. One of the key sets of parameters used in these models is the physical dimensions of the fragments (i.e., length, average-cross sectional area, and volume). For the DebriSat project, only fragments with at least one dimension greater than 2 mm are collected and processed. Additionally, a significant portion of the fragments recovered from the impact test are needle-like and/or flat plate-like fragments where their heights are almost negligible in comparison to their other dimensions. As a result, two fragment size categories were defined: 2D objects and 3D objects. While measurement systems are commercially available, factors such as measurement rates, system adaptability, size characterization limitations and equipment costs presented significant challenges to the project and a decision was made to develop our own size characterization systems. The size characterization systems consist of two automated image systems, one referred to as the 3D imaging system and the other as the 2D imaging system. Which imaging system to use depends on the classification of the fragment being measured. Both imaging systems utilize point-and-shoot cameras for object image acquisition and create representative point clouds of the fragments. The 3D imaging system utilizes a space-carving algorithm to generate a 3D point cloud, while the 2D imaging system utilizes an edge detection algorithm to generate a 2D point cloud. From the point clouds, the three largest orthogonal dimensions are determined using a convex hull algorithm. For 3D objects, in addition to the three largest orthogonal dimensions, the volume is computed via an alpha-shape algorithm applied to the point clouds. The average cross-sectional area is also computed for 3D objects. Both imaging systems have automated size measurements (image acquisition and image processing) driven by the need to quickly and accurately measure tens of thousands of debris fragments. Moreover, the automated size measurement reduces potential fragment damage/mishandling and ability for accuracy and repeatability. As the fragment characterization progressed, it became evident that the imaging systems had to be revised. For example, an additional view was added to the 2D imaging system to capture the height of the 2D object. This paper presents the DebriSat project's imaging systems and calculation techniques in detail; from design and development to maturation. The experiences and challenges are also shared.

Shiotani, B.↗

DOE OSTI2021

Origin of the quasi-quantized Hall effect in ZrTe 5

The quantum Hall effect (QHE) is traditionally considered to be a purely two-dimensional (2D) phenomenon. Recently, however, a three-dimensional (3D) version of the QHE was reported in the Dirac semimetal ZrTe 5 . It was proposed to arise from a magnetic-field-driven Fermi surface instability, transforming the original 3D electron system into a stack of 2D sheets. Here, we report thermodynamic, spectroscopic, thermoelectric and charge transport measurements on such ZrTe 5 samples. The measured properties: magnetization, ultrasound propagation, scanning tunneling spectroscopy, and Raman spectroscopy, show no signatures of a Fermi surface instability, consistent with in-field single crystal X-ray diffraction. Instead, a direct comparison of the experimental data with linear response calculations based on an effective 3D Dirac Hamiltonian suggests that the quasi-quantization of the observed Hall response emerges from the interplay of the intrinsic properties of the ZrTe 5 electronic structure and its Dirac-type semi-metallic character.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

A GPU-based Approach for Path Planning Optimization via Travel Length Reduction

Typically, before constructing an object with an additive manufacturing system, the 3D object must be sent through a process called slicing. Slicing converts a 3D object commonly in the form of an STL file into a set of layers by horizontally intersecting a plane with the object at various heights. At each height, called a layer, multiple 2D polygons can be generated. Each polygon represents a boundary for solid geometry and is called an island. Each island is then comprised of multiple path types in an attempt to optimally fill the polygon. To move between each island and each islands’ paths, travels are inserted. Travels are simply motion by the system to move from one area of construction to another. Travels do not contribute to the construction of the object, and so, are considered wasted motion. In large-scale additive manufacturing, objects can be quite large and the distance between islands can be large as well. As a result, these travels can waste a significant amount of time. Ideally, travels would be as short as possible, however, computing global minimal travel paths is computationally expensive. To combat this problem, researchers at Oak Ridge National Lab developed a GPU-based approach to travel insertion based on a unique factoradic representation. This representation was then utilized by the GPU to solve the Traveling Salesman Problem (TSP). This algorithm was able to compute global minimal travel paths quickly resulting in faster object construction. A general investigation was also carried out to determine when a GPU vs CPU implementation would be beneficial.

Borish, Michael↗

DOE OSTI2021

Flow boiling transient critical heat flux tests with stainless steel and FeCrAl: Transient correlation implementation, model calibration, and sensitivity analysis

In this study, steady-state and transient internal-flow critical heat flux (CHF) experiments were carried out under two flow conditions, at atmospheric pressure, on 316 stainless steel (316-SS) and iron-chromium-aluminum (FeCrAl) tubes. Slow power transients at a low mass flow (300 kg/m 2 - s) with very low subcooling (X e =-0.0054) generated premature CHFs, which were prevented by faster power transients and high mass flow (1000 kg/m 2 - s). The measured transient CHFs increased linearly with increasing power transients, compared with the steady-state CHF. Yet, the wall superheat at the CHFs decreased with faster power transients. Transient CHF correlations highlighting heterogeneous spontaneous nucleation were calibrated to the measured CHFs and compared with other existing correlations. Transient CHF multipliers were acquired from pool and flow boiling empirical CHF correlations that were generated. The multipliers were applied to the RELAP5-3D nuclear system code to analyze the discrepancies between the measured data and the predicted CHF and post-CHF behavior, which improved peak cladding temperature predictions by 24.6%. A variance-based global sensitivity study perturbing the experimental uncertainties and cladding material thermal properties showed the diminishing influence of flow boiling heat transfer with increasing power transients, highlighting the significance of the volumetric heat capacity for the cladding integrity during transients. Transient CHF correlations were applied to the most limiting design basis accident: a hot-zero-power reactivity-initiated accident on a generic pressurized water reactor RELAP5-3D model.

36 MATERIALS SCIENCE↗

DOE OSTI2025

SiC Receiver/Reactor by Additive Manufacturing for Concentrated Solar Thermocatalysis with Thermal Energy Storage (Final Technical Report - Public)

The direct use of solar thermal energy provides opportunities for low-cost heating sources for a variety of applications. Ultra-high temperatures around 1000°C are high value and highly useful for energy-demanding industries. Many materials cannot withstand these conditions. In the area of Sustainable Chemicals, further limitations on material stability exist. Combining state-of-the-art materials with new designs provides a promising pathway for harvesting solar thermal energy and performing high temperature chemical processes. However, conventional manufacturing limits the potential for design flexibility. In this project, Additive Manufacturing was combined with advanced materials and new chemical reactor designs. In addition, 24/7 energy is necessary for chemical processing, and designs for ultra-high temperature thermal storage were devised. Specifically, preliminary design of a novel solar thermal receiver was developed in this project and designed to work with thermocatalytic reactors for producing sustainable chemicals and fuels. An ultra-high temperature particle storage system and heat exchangers were proposed to transport ultra-hot air as thermal fluid for the system. On a broader scale, this system could be used to tap solar thermal energy for a centralized facility with capability of transferring that heat to various segments at a full range of temperatures to 1000°C. The project pushed the temperature boundaries past those in current use, and Additive Manufacturing was envisaged for fabricating the receiver to meet requirements of extreme environments. An extensive analysis of silicon carbide additive manufacturing was performed to compare the thermal and mechanical properties of complex geometries compared to conventional material and those manufactured via other methods. The Additive Manufacturing via Binder-Jet printing was optimized and characterized to provide high quality and reproducible components capable of withstanding the proposed extreme environments. The designs for the concentrating solar thermal cavity with ultra-hot air thermal fluid showed high performance in simulations, attributable to the complex optimized geometries of the 3D printed systems. The bright future of Additive Manufacturing with advanced materials developments should provide more options and even higher quality as the technology further develops. Current costs for Additive Manufacturing of advanced ceramics is relatively low, however post-processing of the materials for extreme environments is currently high. There is little industrial-scale infrastructure for these, but it is growing as niche applications become more mainstream. The results of the project can be translated into similar extreme environments for concentrating solar thermal energy as well as its integration with ultra-hot air thermal fluids. A number of industries that require ultra-high temperatures need to electrify or otherwise decarbonize for climate goals, and this project showed that theoretically there is a pathway to do so with direct concentrated solar thermal power.

10 SYNTHETIC FUELS↗

NASA NTRS2017

Habitable Moist Atmospheres on Terrestrial Planets near the Inner Edge of the Habitable Zone Around M Dwarfs

Terrestrial planets in the habitable zones (HZs) of low-mass stars and cool dwarfs have received significant scrutiny recently. Transit spectroscopy of such planets with the James Webb Space Telescope (JWST) represents our best shot at obtaining the spectrum of a habitable planet within the next decade. As these planets are likely tidally locked, improved 3D numerical simulations of such planetary atmospheres are needed to guide target selection. Here we use a 3D climate system model, updated with new water-vapor absorption coefficients derived from the HITRAN 2012 database, to study ocean-covered planets at the inner edge of the HZ around late M to mid- K stars (2600 K less than or equal to T(sub eff) less than or equal to 4500 K). Our results indicate that these updated water-vapor coefficients result in significant warming compared to previous studies, so the inner HZ around M dwarfs is not as close as suggested by earlier work. Assuming synchronously rotating Earth-sized and Earth-mass planets with background 1 bar N2 atmospheres, we find that planets at the inner HZ of stars with T(sub eff) greater than 3000 K undergo the classical ''moist greenhouse'' (H2O mixing ratio greater than 10(exp -3) in the stratosphere) at significantly lower surface temperature (approx. 280 K) in our 3D model compared with 1D climate models (approx. 340 K). This implies that some planets around low-mass stars can simultaneously undergo water loss and remain habitable. However, for stars with T(sub eff) less than or equal to 3000 K, planets at the inner HZ may directly transition to a runaway state, while bypassing the moist greenhouse water loss entirely. We analyze transmission spectra of planets in a moist greenhouse regime and find that there are several prominent H2O features, including a broad feature between 5 and 8 micron, within JWST MIRI instrument range. Thus, relying only on standard Earth-analog spectra with 24 hr rotation period around M dwarfs for habitability studies will miss the strong H2O features that one would expect to see on synchronously rotating planets around M dwarf stars, with JWST.

Kopparapu, Ravi kumar↗

DOE OSTI2022

Quasiparticle characteristics of the weakly ferromagnetic Hund metal MnSi

Hund metals are multiorbital systems with 3d or 4d electrons exhibiting both an itinerant character and local moments, and they feature Kondo-like screenings of local orbital and spin moments, with suppressed coherence temperatures driven by Hund's coupling J H . They often exhibit magnetic order at low temperature, but how the interaction between the Kondo-like screening and long-range magnetic order is manifested in the quasiparticle spectrum remains an open question. Here, we present the spectroscopic signature of such an interaction in a Hund metal candidate MnSi exhibiting weak ferromagnetism. Our photoemission measurements reveal renormalized quasiparticle bands near the Fermi level with strong momentum dependence: The ferromagnetism manifests through possibly exchange-split bands (Q1) below T C , while the spin/orbital screenings lead to the gradual development of quasiparticles (Q2) upon cooling. Our results demonstrate how the characteristic spin/orbital coherence in a Hund metal could coexist and compete with the magnetic order to form a weak itinerant ferromagnet, via quasiparticle bands that are well separated in momentum space and exhibit distinct temperature dependence. Here, our results imply that the competition between the spin/orbital screening and the magnetic order in a Hund metal bears interesting similarities to the Kondo lattice systems

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1992

Virtual environment display for a 3D audio room simulation

The development of a virtual environment simulation system integrating a 3D acoustic audio model with an immersive 3D visual scene is discussed. The system complements the acoustic model and is specified to: allow the listener to freely move about the space, a room of manipulable size, shape, and audio character, while interactively relocating the sound sources; reinforce the listener's feeling of telepresence in the acoustical environment with visual and proprioceptive sensations; enhance the audio with the graphic and interactive components, rather than overwhelm or reduce it; and serve as a research testbed and technology transfer demonstration. The hardware/software design of two demonstration systems, one installed and one portable, are discussed through the development of four iterative configurations.

Chapin, William L.↗

NASA NTRS2019

Astrobiology Survey of a Lava Cave at Lava Beds National Monument by a Rover Carrying a Remote Sensing Instrument Payload

We report here on a survey of a lava tube cave by a rover that is instrumented for astrobiology missions. The NASA Ames testbed rover, CaveR, was deployed in Valentine Cave in Lava Beds National Monument (N. CA, USA) during August of 2018. The rover carried an instrument package consisting of Near Infrared and Visible Spectrometer System (NIRVSS) a point spectrometer operating in 1590-3400 nm range, sensitive to H2O and -OH bearing minerals, pyroxenes, and carbonates (Roush, et al 2018); the bore sighted Drill Operations Camera (DOC), a monochrome imager illuminated by LEDs at 410, 540, 640, 740, 905 and 940 nm; a Realsense ™ depth sensor system for 3D model generation; and a high resolution DSLR stereo camera. The payload was mounted on a tiltable instrument platform attached to the left side of the rover. The rover was driven manually in the cave by field operators, following instructions from a remote science operations team, and simulating a mission concept with science-guided autonomy. A simulated mission took place for 3 days with a team of 3 scientists selecting targets and interpreting data from the payload. To begin the mission, the rover drove along one wall of the cave imaging continuously with the Realsense in 20 m cave segments, three total. At the start of each day, the images from a 20m segment and a panorama stitched from them were provided to the science team to examine. The science team used these data to prioritize specific points along the cave wall for the collection of NIRVSS, DOC, and DSLR data. The objective of the data collection was to identify and study putative biological and mineralogical features in the cave. The data were delivered in xGDS, a customized mapping, planning, and data base management software developed at NASA Ames (Lee, et al 2013). Once the targets for further observations were selected, a plan for collecting the observations (positions in the cave and pointing for each requested observation) was constructed using xGDS and delivered to a rover team to execute the science data collection plan. Acquired data were delivered back to the science team for analysis. Preliminary results from the experiment illustrate the utility of the system (rover plus payload) to study the cave geology and mineralogy and its potential for identifying biomineral features.

Wong, Uland↗

DOE OSTI2021

System and method for stimulated emission depletion projection stereolithography

The present disclosure relates to a system for forming a three dimensional (3D) part. The system may incorporate a beam delivery subsystem for generating optical signals, and a mask subsystem that receives the optical signals and generates optical images therefrom. A first one of the optical images activates a polymerization species of a photo-sensitive resin in accordance with illuminated areas thereof, to thus cause polymerization of select portions of the photo-sensitive resin to help form a layer of the 3D part. A second one of the optical images causes stimulated emission depletion of subportions of the polymerization species, simultaneously, over various areas of the layer, to enhance resolution of at least one subportion of the select portions of the photo-sensitive resin.

Hensleigh, Ryan↗

NASA NTRSDate not supplied

Pterodactyl: Effects of 3D Thermal Analysis on Thermal Protection System Design for a Flap Control System

NASA’s Pterodactyl project has investigated the use of a novel multi-flap control system to facilitate precision targeting of Deployable Entry Vehicles (DEVs) during atmospheric re-entry [1]. DEVs can be folded to fit within the limiting cross-sectional area of current launch systems. Once deployed, the vehicle expands and settles into a pre-determined blunt body shape. While DEVs provide a more efficient solution to increased payload sizes, the absence of a back shell does not allow for easy integration of reaction control systems, which have historically been used for guidance and control (G&C) of rigid aeroshells during entry. One DEV solution, called the Adaptable, Deployable Entry Placement Technology (ADEPT), employs mechanically deployed gores. This provided the Pterodactyl project with the opportunity to incorporate a rib-mounted, 8-flap control system for G&C. Each flap can deflect independently in and out of the hypersonic flow, requiring TPS for the flaps. Previous work within the Pterodactyl project utilized 1D thermal analysis to design the TPS. However, the Pterodactyl project was concerned that the base level tools for 1D thermal analysis were not accounting for the 3D effects of the large heating gradients along and across the flaps, the in-depth in-plane conductivity of the TPS and internal structure, and the effects of the small edge radii at the neck of the flaps. This paper discusses the methodology and results of a 3D thermal analysis of the flap control system. Additionally, the results from the 1D and 3D thermal analyses are compared. It is found that elements of the original design that resulted from the 1D analysis may be overly conservative, and implementation of a 3D thermal analysis indicated a reduction in TPS thickness is feasible. This result helps to reduce the mechanical integration complexity of the flap at the rib tip and realizes potential mass savings.

Sarah N D'Souza↗

DOE OSTI2022

Physical Models from Physical Templates Using Biocompatible Liquid Crystal Elastomers as Morphologically Programmable Inks For 3D Printing

Advanced manufacturing has received considerable attention as a tool for the fabrication of cell scaffolds however, finding ideal biocompatible and biodegradable materials that fit the correct parameters for 3D printing and guide cells to align remain a challenge. Herein, a photocrosslinkable smectic-A (Sm-A) liquid crystal elastomer (LCE) designed for 3D printing is presented, that promotes cell proliferation but most importantly induces cell anisotropy. The LCE-based bio-ink allows the 3D duplication of a highly complex brain structure generated from an animal model. Vascular tissue models are generated from fluorescently stained mouse tissue spatially imaged using confocal microscopy and subsequently processed to create a digital 3D model suitable for printing. The 3D structure is reproduced using a Digital Light Processing (DLP) stereolithography (SLA) desktop 3D printer. Synchrotron Small-Angle X-ray Diffraction (SAXD) data reveal a strong alignment of the LCE layering within the struts of the printed 3D scaffold. The resultant anisotropy of the LCE struts is then shown to direct cell growth. This study offers a simple approach to produce model tissues built within hours that promote cellular alignment.

3D cell culture↗

DOE OSTI2025

Results of an interlaboratory study on the working curve in vat photopolymerization II: Towards a standardized method

The working curve measurement in photopolymer additive manufacturing is a ubiquitous measure of the cure depth of a printing resin as a function of radiant exposure of light. The fit parameters from this measurement (the depth of light penetration D p and the critical exposure E c ) are used to evaluate and report a resin’s printability, optimize processing parameters, and inform print and resin quality control. Despite its widespread use, the working curve lacks a standard measurement method. Here, following up on our paper “Results of an Interlaboratory Study on the Working Curve in Vat Photopolymerization” from last year, an interlaboratory study on the working curve was performed using calibrated, reproducible, bandpass filtered light sources. With these light sources, the variability between labs in measured working curves was dramatically reduced from the initial interlaboratory study. Aggregate data from this experiment produced reliable D p and E c measurements at 385 nm of 39.2 ± 3.7 µm and 12.3 ± 3.0 mJ cm −2 , respectively. At 405 nm the values of D p and E c are 69.3 ± 3.8 µm and 17.9 ± 2.3 mJ cm −2 , respectively. The results are agnostic to the thickness measurement tool utilized by participants, ensuring broad applicability across laboratories. We also tested the generalizability of the proposed method of using a filtered light source by filtering a commercial 405 nm light source and obtaining a working curve in agreement with the aggregate data from the interlaboratory study. This interlaboratory study provides a basis for a documentary standard for the working curve, so that the entire photopolymer additive manufacturing industry can share reproducible and interoperable working curve data.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Mapping and Characterizing the Visual Impacts of the Existing US Wind Turbine Fleet

Visual impacts of wind turbines have been a persistent concern for wind energy development in the United States (US) for decades and remain a major source of project delays and cancellations. Assessments of visual impacts are frequently performed at a local scale for individual projects, but a comprehensive understanding of broader geographic patterns in visual impacts across the US is lacking. This paper presents a visual impact assessment of the existing land-based wind turbine fleet of the contiguous United States (CONUS). The assessment combines geographic information systems and 3D simulation methods to account for key factors driving the visual magnitude of impacts from the installed turbines. The results indicate that, despite the deployment of approximately 70,000 turbines and over 144 gigawatts of land-based wind in the CONUS, the visual impacts are very small when measured as a proportion of land area, population, and sensitive visual resources. Nonetheless, visual impacts are not evenly distributed: people experience a concentrated share in a small number of natural settings, primarily including plains, prairies, and deserts. Finally, we find that although increased density of wind development consistently leads to visual impacts across a greater proportion of land, it does not always lead to impacts to a greater share of the population. These findings suggest that visual impacts from wind energy are generally well-mitigated across the CONUS to date but also highlight the need for a deeper understanding of landscape sensitivity and individual perceptions of wind turbines in the most heavily impacted natural settings.

3D simulation↗

DOE OSTI2023

Insights Into Lightning K‐Leader Initiation and Development From Three Dimensional Broadband Interferometric Observations

Abstract We report detailed observations of K‐leaders and the activity between them with the three‐dimensional Broadband Interferometric Mapping and Polarization system (BIMAP‐3D) at Los Alamos National Laboratory. It is found that K‐leaders have a general propagation trend of initial acceleration and then gradual deceleration, and the corresponding very high frequency (VHF) radiation power is exponentially correlated with the leader speed. Based on the 3D development and simultaneous electric field change measurement, some simple K‐leaders can be modeled with time‐evolving point charges at the propagating leader tip and at the stationary origin. We found that the charge magnitude increases during the initial acceleration stage and stays relatively constant for the rest of the development. K‐leaders are observed to interact with other branches; the branches affect the leader's propagation speed, and may affect the charge transfer. After the occurrence of a K‐leader, VHF emissions are quenched for several milliseconds. VHF sources then reappear in an impulsive and scattered manner as “twinkling,” and these sources are found not uniquely on the so‐called needles, but also on the main channel. These twinkling sources start near the apparent positive leader tip, and extend back toward the direction of the flash origin at about 10 5 m/s, while the apparent positive tip continues to extend forward at about 10 4 m/s. The twinkling extending toward the direction of flash origin appears to initiate the following K‐leader, although it may be interrupted by a K‐leader along a different branch, or simply die out without more K‐leader activity.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Lattice distortions and metal-insulator transition in hexagonal FeS

Using total x-ray scattering and density functional theory calculations based on experimental data, we study the metal-insulator transition (MIT) in the strongly correlated system FeS with hexagonal symmetry. We find that it arises from the combined effect of static lattice distortions and antiferromagnetic ordering, while electron-electron correlations appear to play a minor role. The presence of local lattice distortions above the MIT explains well the evolution of transport and magnetic properties with changing temperature or pressure. Intrinsically, lattice distortions in hexagonal FeS act as lattice degrees of freedom that bridge competing electronic phases, thus facilitating transitions between them. The mechanism of MIT operating in h-FeS appears to be common to strongly correlated binary systems involving 3d transition metals, warranting further investigations on the lattice distortions-property relationship in this broad class of materials using the approach adopted here.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Exchange-correlation functional challenges in modeling quaternary chalcogenides

The development of next-generation quaternary chalcogenides, such as Cu2ZnSnS4 (CZTS) and Cu2ZnGeS4 (CZGS), for solar energy and thermoelectric applications hinges upon both careful experimentation and accurate quantum mechanical modeling. To address the latter, many have turned to density functional theory (DFT), which offers several choices for the approximate treatment of electron exchange and correlation (XC). Popular XC functionals include the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA) and the recently developed strongly constrained and appropriately normed (SCAN) meta-GGA. Extensions of DFT functionals, such as adding a Hubbard U correction and introducing a fraction of the Fock exchange (hybrid functionals), have been used widely to model systems containing 3d metal ions. However, no studies yet have compared comprehensively PBE(+U) and SCAN(+U) in the quality of their predictions of the bulk and defect thermodynamics of quaternary chalcogenides, which play a critical role in device fabrication and performance. Hence, here we calculate the (i) 0 K formation energies of bulk Ge compounds and (ii) neutral defect formation energies including charge-balanced (e.g., CuZn + ZnCu) and charge-imbalanced (e.g., CuSn) combinations of antisites and vacancies in CZTS and CZGS using the PBE, PBE +U, SCAN, SCAN +U, and the hybrid Heyd-Scuseria-Ernzerhof XC frameworks. We find that the formation energies of charge-imbalanced defects are more sensitive to the choice of the XC functional than those of charge-balanced defects, which can be explained by the differences in the extent of penalization of defect-generated delocalized electrons/holes by PBE, PBE +U, SCAN, and SCAN +U. Additionally, our results show that SCAN systematically underbinds Ge-containing compounds, thus highlighting the need for even further improvement of XC functionals. Based on our findings, we recommend the use of SCAN for modeling quaternary chalcogenides because its errors are systematic, and it has the firmest theoretical underpinning. Our work provides guidance for future modeling of quaternary chalcogenides.

14 SOLAR ENERGY↗