Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Operator approximation”

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

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

At least 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Empirical Study of the Effect of Thermal Loading on the Heating Efficiency of Variable-Speed Air Source Heat Pumps

Heating buildings with air source heat pumps (ASHPs) has the potential to save energy compared to utilizing conventional heat sources. Accurate understanding of the efficiency of ASHPs is important to maximize the energy savings. While it is well understood that, in general, ASHP efficiency decreases with decreasing outdoor temperature, it is not well understood how the ASHP efficiency changes with different levels of thermal loading, even though it is an important consideration for sizing and controlling ASHPs. The goal of this study was to create an empirical model of the ASHP efficiency as a function of two independent variables–outside temperature and level of thermal loading. Four ductless mini-split ASHPs were evaluated in a cold chamber where the temperature (representing the outdoor temperature) was varied over a wide range. For each temperature, the ASHP performance data were collected at several levels of thermal loading. The data for all four ASHPs were combined and approximated with an analytical function that can be used as a general model for the ASHP steady-state efficiency as a function of the outside temperature and level of thermal loading. To the knowledge of the authors, no such empirical model that is solely based on third-party test data has been published before. While limitations exist, the model can be used to help guide future selection and operation of ASHPs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Idaho National Laboratory’s Mobile Hot Cell Transportation: Engineering Solutions for Global Disused Sealed Radioactive Sources.

Title: Idaho National Laboratory’s Mobile Hot Cell Transportation: Engineering Solutions for Global Disused Sealed Radioactive Sources. Abstract: The Mobile Hot Cell (MHC), currently under development by Idaho National Laboratory (INL) for the Off-Site Source Recovery Project (OSRP), is designed to help international partners meet the unique challenges of end-of-life radioactive material management. The MHC will provide a critical resource for countries that require assistance securing and disposing of Disused Sealed Radioactive Sources (DSRS) and orphaned sources in challenging environments, allowing these sources to be secured against misuse and nefarious activities. The MHC is a rapidly deployable system for conditioning and preparing end-of-life radioactive sources for transportation or storage. It is designed to handle sources of up to 1,000 Ci Co-60 equivalent while maintaining full radiological and biological containment. It will be deployable within 48 hours of an alert, making it ideal for emergency situations. The MHC features an operational suite for control, support racks for electronics, pneumatics, and welding systems, and a modular robust steel structure providing radiological shielding and internal robotic support. This design allows configurations for multiple device types to be conditioned and the ability to safely manage routine issues such as leaking or damaged sources. The MHC has been designed with the transportation challenges of rapid deployment to difficult environments in mind. The system weighs approximately 150,000 pounds, with individual systems breaking down into pieces not exceeding 20,000 pounds. Components are to be transportable on standard 20ft ISO containers, with shielding shells on 20ft flat racks. It is estimated that a total of eight containers and flat racks will be required. The use of 20ft containers, as opposed to 40ft containers, minimizes the impact on less developed road infrastructures, enabling the MHC to be positioned in constrained environments such as hospital parking lots. The system’s modularity also allows for deployment using smaller equipment, such as a 10-ton boom truck or forklift, which is crucial given the potential logistical challenges in different countries. This transportation strategy, evaluated in collaboration with Utah State University, ensures the MHC can be deployed via ground, rail, sea, or air, addressing the primary concern of international transport logistics.

99 - GENERAL AND MISCELLANEOUS↗
DOE OSTI2020

Investigation of fouling mechanisms for diesel engine exhaust heat recovery

The effect of diesel engine exhaust fouling on heat and mass exchangers for waste heat recovery is investigated. Experiments are conducted for a range of operating conditions relevant to the use of the waste heat from a diesel engine for powering a 2.71 kW capacity absorption heat pump. Thermal resistances of fouling layers as well as their effects on pressure drop are determined, and the underlying mechanisms for deposition and removal are investigated. Fouling layer effects are most severe for the lowest coolant temperatures and 100% generator load. At these conditions, the fouling thermal resistance is approximately 70% of the total resistance, and the ratio of the steady-state pressure drop to the baseline pressure drop is 3.25. Flow induced shear is found to be the primary growth inhibiting mechanism at high exhaust velocities. Findings from this study can guide heat exchanger design for a variety of waste heat recovery systems.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2024

Bayesian reduced-order deep learning surrogate model for dynamic systems described by partial differential equations

We propose a reduced-order deep-learning surrogate model for dynamic systems described by time-dependent partial differential equations. This method employs space–time Karhunen–Loève expansions (KLEs) of the state variables and space-dependent KLEs of space-varying parameters to identify the reduced (latent) dimensions. Subsequently, a deep neural network (DNN) is used to map the parameter latent space to the state variable latent space. An approximate Bayesian method is developed for uncertainty quantification (UQ) in the proposed KL-DNN surrogate model. The KL-DNN method is tested for the linear advection–diffusion and nonlinear diffusion equations, and the Bayesian approach for UQ is compared with the deep ensembling (DE) approach, commonly used for quantifying uncertainty in DNN models. It was found that the approximate Bayesian method provides a more informative distribution of the PDE solutions in terms of the coverage of the reference PDE solutions (the percentage of nodes where the reference solution is within the confidence interval predicted by the UQ methods) and log predictive probability. The DE method is found to underestimate uncertainty and introduce bias. For the nonlinear diffusion equation, we compare the KL-DNN method with the Fourier Neural Operator (FNO) method and find that KL-DNN is 10% more accurate and needs less training time than the FNO method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Spontaneous Polarization Suppression of Exciton–Exciton Annihilation in Rhombohedral-Stacked Bilayer Molybdenum Disulfide

Rapid exciton–exciton annihilation (EEA) in two-dimensional semiconductors limits access to high-density excitonic regimes essential for efficient optoelectronic operation under strong excitation. Here, we show that EEA is strongly suppressed in rhombohedral (3R)-stacked MoS 2 bilayers relative to nonpolar 2H bilayers, and that this suppression is quantitatively consistent with repulsive dipole–dipole interactions between layer-polarized excitons induced by spontaneous polarization. Using ultrafast pump–probe spectroscopy, we measure an EEA rate of γ EEA = (5.03 ± 0.99) × 10 –3 cm 2 s –1 in 3R bilayers, which is approximately 18.2-fold smaller than that in monolayers and 2.9-fold smaller than that in nonpolar 2H bilayers. Despite the higher exciton diffusivity recently reported for 3R relative to 2H bilayers, the reduced EEA rate in 3R indicates a rate-limited regime governed by the close-encounter annihilation probability rather than diffusion. A rate-limited annihilation model incorporating a dipole–dipole repulsive potential captures the observed ratio γ EEA,3R /γ EEA,2H ≈ 0.35 for an exciton–exciton encounter distance of ∼1.3 nm, consistent with the bilayer exciton Bohr radius. These results indicate that spontaneous polarization in 3R-stacked bilayers can suppress nonlinear excitonic losses and provide a route toward high-density excitonics.

dipolar excitons↗
DOE OSTI2022

Scaling quantum approximate optimization on near-term hardware

The quantum approximate optimization algorithm (QAOA) is an approach for near-term quantum computers to potentially demonstrate computational advantage in solving combinatorial optimization problems. However, the viability of the QAOA depends on how its performance and resource requirements scale with problem size and complexity for realistic hardware implementations. Here, we quantify scaling of the expected resource requirements by synthesizing optimized circuits for hardware architectures with varying levels of connectivity. Assuming noisy gate operations, we estimate the number of measurements needed to sample the output of the idealized QAOA circuit with high probability. We show the number of measurements, and hence total time to solution, grows exponentially in problem size and problem graph degree as well as depth of the QAOA ansatz, gate infidelities, and inverse hardware graph degree. These problems may be alleviated by increasing hardware connectivity or by recently proposed modifications to the QAOA that achieve higher performance with fewer circuit layers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Modeling of particle transport, neutrals and radiation in magnetically-confined plasmas with Aurora

In this work, we present Aurora, an open-source package for particle transport, neutrals and radiation modeling in magnetic confinement fusion plasmas. Aurora's modern multi-language interface enables simulations of 1.5D impurity transport within high-performance computing frameworks, particularly for the inference of particle transport coefficients. A user-friendly Python library allows simple interaction with atomic rates from the Atomic Data and Atomic Structure database as well as other sources. This enables a range of radiation predictions, both for power balance and spectroscopic analysis. We discuss here the superstaging approximation for complex ions, as a way to group charge states and reduce computational cost, demonstrating its wide applicability within the Aurora forward model and beyond. Aurora also facilitates neutral particle analysis, both from experimental spectroscopic data and other simulation codes. Leveraging Aurora's capabilities to interface SOLPS-ITER results, we demonstrate that charge exchange is unlikely to affect the total radiated power from the ITER core during high performance operation. Finally, we describe the ImpRad module in the one modeling framework for integrated task framework, developed to enable experimental analysis and transport inferences on multiple devices using Aurora.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Color Transformations of Photometric Measurements of Galaxies in Optical and Near-infrared Wide-field Imaging Surveys*

Abstract Over the past 2 decades, wide-field photometric surveys in optical and infrared domains reached a nearly all-sky coverage thanks to numerous observational facilities operating in both hemispheres. However, subtle differences among exact realizations of Johnson and SDSS photometric systems require one to convert photometric measurements into the same system prior to analysis of composite data sets originating from multiple surveys. It turns out that the published photometric transformations lead to substantial biases when applied to integrated photometry of galaxies from the corresponding catalogs. Here we present photometric transformations based on piece-wise linear approximations of integrated photometry of galaxies in the optical surveys SDSS, DECaLS, BASS, MzLS, DES, DELVE, KiDS, VST ATLAS, and the near-infrared surveys UKIDSS, UHS, VHS, and VIKING. We validate our transformations by constructing k -corrected color–magnitude diagrams of non-active galaxies and measuring the position and tightness of the “red sequence.” We also provide transformations for aperture magnitudes and show how they are affected by the image quality difference among the surveys. We present the implementation of the derived transformations in python and idl and also a web-based color transformation calculator for galaxies. By comparing DECaLS and DES, we identified systematic issues in DECaLS photometry for extended galaxies, which we attribute to the photometric software package used by DECaLS. As an application of our method, we compiled two multi-wavelength photometric catalogs for over 200,000 low- and intermediate-redshift galaxies originating from CfA FAST and Hectospec spectral archives.

Astronomy & Astrophysics↗
DOE OSTI2021

Hyper Suprime-Cam Legacy Archive

Abstract We present the launch of the Hyper Suprime-Cam Legacy Archive (HSCLA), a public archive of processed, science-ready data from the Hyper Suprime-Cam (HSC). HSC is an optical wide-field imager installed at the prime focus of the Subaru Telescope and has been in operation since 2014. While ∼1/3 of the total observing time of HSC has been used for the Subaru Strategic Program (SSP), the remainder of the time is used for Principal Investigator (PI)-based programs. We have processed the data from these PI-based programs and make the processed, high-quality data available to the community through HSCLA. The current version of HSCLA includes data taken in the first year of science operation, 2014. We provide both individual and coadd images as well as photometric catalogs. The photometric catalog from the coadd is loaded to the database, which offers a fast access to the large catalog. There are other online tools such as an image browser and an image cutout tool and they will be useful for science analyses. The coadd images reach 24–27th magnitudes at 5σ for point sources and cover approximately 580 square degrees in at least one filter with 150 million objects in total. We perform extensive quality assurance tests and verify that the photometric and astrometric quality of the data is good enough for most scientific explorations. However, the data are not without problems and users are referred to the list of known issues before exploiting the data for science. All the data and documentations can be found at the data release site, 〈https://hscla.mtk.nao.ac.jp/〉.

Tanaka, Masayuki↗
DOE OSTI2022

PlasMod0 v.1.0

SAND2022-12609 O PlasMod0 is a 0-dimensional plasma modeling which approximates a plasma as a well-stirred reactor. The model uses two temperatures, one for electrons and one for ions and neutrals. This model was originally developed for laser produced plasma in single-particle aerosol mass spectrometers. As a result, there are some specialized capabilities that include spherical expansion into vacuum at a user-specified rate and photon reactions. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Lietz, Amanda↗
DOE OSTI2020

Imagery Applications for Advanced Event Analytics (WBS 24.3.1.3.3-IDC FY202 Final Project Report)

Accurate event locations and replicability of location analyses are essential for assessing the nature of an event, its context, ambient site conditions, and proximity to relevant facilities and infrastructure. Additionally, accurate event locations provide valuable information that reduce uncertainties, improve confidence in event analyses, and inform in-field verification activities. However, event location/relocation and replicability are difficult due to a number of factors, including spatially-sparse network coverage in some areas of the globe and variability in seismic data processing. This team proposed that the incorporation of high-fidelity imagery as a data backbone to the analytical assessment of a suspected underground explosion and/or an advanced seismic event bulletin produced by the International Data Centre (IDC) of the Preparatory Commission for the Comprehensive Nuclear-Test-Ban Treaty Organization (CTBTO PrepCom) could reduce uncertainties and improve confidence in analyses. Specifically, temporally-separated images can reduce uncertainty by identifying areas where change has occurred (e.g., building construction or demolition, road or facilities improvements). The primary goal of this project was to develop an automated geospatial processing script for imagery change detection to better reflect needs of the technical community (including the IDC) and to make the use of such a tool accessible in a variety of settings across platforms. Technical experts at Los Alamos National Laboratory successfully built GAIA: the Geospatial Automated Imagery Analysis tool, to fill this need. GAIA combines five tool components to produce orthorectified time separated imagery and imagery change detection maps. Our toolkit (1) reduces error by providing a standardized workflow for image analyses and (2) significantly reduces processing time from between 7 and 24+ hours to approximately 5 minutes. Technical experts at Sandia National Laboratories supported GAIA via beta-testing and by introducing a web-based system approach for increased applicability. To test the function, performance, broad application, and ease-of-use of GAIA, we applied it to four separate test cases. The results of this preliminary investigation show promise in reducing uncertainty in seismic event locations: if satellite imagery can show regions where operations that produce seismic activity likely occurred, then pursuing imagery to locate epicenters of seismic nuclear events could reduce the time needed to find the true epicenter location.

58 GEOSCIENCES↗
DOE OSTI2020

Imagery Applications for Advanced Event Analytics (WBS 24.3.1.3.3-IDC FY2020 Final Project Report)

Accurate event locations and replicability of location analyses are essential for assessing the nature of an event, its context, ambient site conditions, and proximity to relevant facilities and infrastructure. Additionally, accurate event locations provide valuable information that reduce uncertainties, improve confidence in event analyses, and inform in-field verification activities. However, event location/relocation and replicability are difficult due to a number of factors, including spatially-sparse network coverage in some areas of the globe and variability in seismic data processing. This team proposed that the incorporation of high-fidelity imagery as a data backbone to the analytical assessment of a suspected underground explosion and/or an advanced seismic event bulletin produced by the International Data Centre (IDC) of the Preparatory Commission for the Comprehensive Nuclear-Test-Ban Treaty Organization (CTBTO PrepCom) could reduce uncertainties and improve confidence in analyses. Specifically, temporally-separated images can reduce uncertainty by identifying areas where change has occurred (e.g., building construction or demolition, road or facilities improvements). The primary goal of this project was to develop an automated geospatial processing script for imagery change detection to better reflect needs of the technical community (including the IDC) and to make the use of such a tool accessible in a variety of settings across platforms. Technical experts at Los Alamos National Laboratory successfully built GAIA: the Geospatial Automated Imagery Analysis tool, to fill this need. GAIA combines five tool components to produce orthorectified time-separated imagery and imagery change detection maps. Our toolkit (1) reduces error by providing a standardized workflow for image analyses and (2) significantly reduces processing time from between 7 and 24+ hours to approximately 5 minutes. Technical experts at Sandia National Laboratories supported GAIA via beta-testing and by introducing a web-based system approach for increased applicability. To test the function, performance, broad application, and ease-of-use of GAIA, we applied it to four separate test cases. The results of this preliminary investigation show promise in reducing uncertainty in seismic event locations: if satellite imagery can show regions where operations that produce seismic activity likely occurred, then pursuing imagery to locate epicenters of seismic nuclear events could reduce the time needed to find the true epicenter location.

58 GEOSCIENCES↗
DOE OSTI2021

Oak Ridge National Laboratory Roane and Anderson Counties, Tennessee Historic Architectural Resource Survey June 2018

From February through September 2017, Cultural Resource Analysts, Inc. personnel completed a historic architectural resource survey of the Oak Ridge National Laboratory in Roane and Anderson Counties, Tennessee. The survey was completed at the request of UT-Battelle, LLC, the private not-for-profit company that manages and operates Oak Ridge National Laboratory for the U.S. Department of Energy. The purpose of the survey was to assist the Department of Energy in complying with the National Historic Preservation Act of 1966, as amended, by providing updated recommendations regarding the National Register of Historic Places eligibility of historic architectural resources at Oak Ridge National Laboratory. The survey included the entirety of Oak Ridge National Laboratory’s main campus that occupies approximately 1,100 acres in Bethel and Melton Valleys and on Chestnut Ridge. Based on Facilities Information Management System data and guidance provided by UT-Battelle, LLC, preservation and modernization staff, the survey team developed a list of resources for inclusion in the survey. The list included all substantial buildings constructed in or before 1980 in order to provide sufficient coverage of buildings that will turn 50 years of age within the next decade. Only those buildings with recorded area over 100 sq ft were included in the survey. Additionally, the UT-Battelle team identified 15 buildings for more intensive survey because such buildings are central to Oak Ridge National Laboratory’s future mission needs and are expected to undergo renovations in the foreseeable future. These building were subject to both exterior and interior documentation with particular focus on identifying character-defining features, and thus were subject to more detailed analysis than the other resources included in the survey. The findings of this historic resource survey build on the conclusions of the 1994 survey by DuVall & Associates, Inc., as well as the survey updates completed by Thomason and Associates in 2004 and 2015 in support of the Oak Ridge National Laboratory Historic Preservation Plan. Within 12 distinct areas of Oak Ridge National Laboratory, as well as outlying areas, the survey recorded 195 buildings and structures constructed in or before 1980. In consideration of the role and significance of each area and each individual resource within the context of Oak Ridge National Laboratory, and in light of their current integrity and condition, each recorded resource was evaluated to determine if it is eligible for listing in the National Register of Historic Places.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Field Laboratory for Emerging Stacked Unconventional Plays in Central Appalachia

The goal of the Field Laboratory for Emerging Stacked Unconventional Plays (ESUP) in Central Appalachia project was to investigate and characterize the resource potential for multi-play production of emerging unconventional reservoirs in Central Appalachia. Project activities included drilling, logging and coring of a vertical characterization well drilled to basement to approximately 15,000 feet in depth. The data from drilling and well logs as well recovered core samples were used by the research team to support characterization of the geology and potential pay zones within the Nora Gas Field of southwestern Virginia and the greater Central Appalachian Basin. This project was led by Dr. Nino Ripepi of the Virginia Center for Coal and Energy Research (VCCER) in close collaboration with EnerVest Operating LLC (EnerVest). VCCER is housed in the Department of Mining and Minerals Engineering at Virginia Tech and was created by an Act of the Virginia General Assembly on March 30, 1977, as an interdisciplinary study, research, information, and resource facility for the Commonwealth of Virginia. EnerVest is a top tier, low-cost oil and natural gas company with a long history of traditional operating relationships with institutional investors as well as public and private companies. This report summarizes two major research activities: 1) Log and Core Analysis and 2) Numerical Modeling. Core analysis comprises the characterization efforts of the ESUP Field Laboratory and includes a discussion of the well logs run in the deep characterization wells, an initial analysis of those, the number, size and location of cores recovered, core testing and analysis done to date as well as plans for future analysis. Numerical modeling is comprised of seven major modeling efforts, including (1) Multiphysics shale transport modeling work, (2) the developed in-house compositional simulator studies, (3) reservoir modeling studies using GEM software, (4) fracture modeling work using EFRAC3D software and (5) numerical modeling efforts with ABAQUS software, (6) Numerical modeling efforts using FLAC3D and PFC3D, and (7) Automatic Machine Learning (AutoML) Studies.

02 PETROLEUM↗
DOE OSTI2025

RTN-095: The Vera C. Rubin Observatory Data Preview 1

We present Rubin Data Preview 1 (DP1), the first release of data from the NSF-DOE Vera C. Rubin Observatory, consisting of raw and calibrated single-epoch images, coadds, difference images, detection catalogs, and other derived data products. DP1 is based on 1792 science-grade optical/near-infrared exposures acquired over 48 distinct nights by the Rubin Commissioning Camera, LSSTComCam, on the Simonyi Survey Telescope at the Summit Facility on Cerro Pachón, Chile during the first on-sky commissioning campaign in late 2024. DP1 covers a total of ~15 sq. deg. over seven roughly equally-sized non-contiguous fields, each independently observed in six broad photometric bands, ugrizy, spanning a range of stellar densities and latitudes and overlapping with external reference datasets. The median image quality across all bands, measured by the FWHM of the point-spread function, is approximately 1.13 arcseconds, with the sharpest images reaching about 0.65 arcseconds. DP1 contains approximately 2.3 million distinct astrophysical objects, of which 1.6 million are extended in at least one band, and 431 solar system objects, of which 93 are new discoveries. DP1 is approximately 3.5 TB in size and available to Rubin data rights holders via the Rubin Science Platform, a cloud-based environment for the analysis of petascale astronomical data. While small compared to future LSST releases, its high quality and diversity of data support a broad range of early science investigations across all four LSST themes, providing a valuable opportunity to engage with Rubin data ahead of the start of full operations in late 2025.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Dynamic System Scaling Application to Accelerated Nuclear Fuel Testing

The development of nuclear fuel and materials requires a continuous effort to investigate the acute and prolonged effects of irradiation, thermal-material stress, chemical change, or other conceivable damage mechanics acquired during normal operation or accident scenarios throughout its lifetime. As in-core fuel property measurement techniques advance to support in the real-time, non-invasive, and enhanced accuracy realm, it is the future of fuel development to pursue to higher degrees of control, predictability of integral test behavior via separate effects test (SET), and shorter test time intervals. The fuel development life cycle from initial concept to commercial licensing is approximated to be 20 years and current literature suggests by optimizing fuel performance codes with SETs, the process could possibly be compressed to 5 to 10 years. Recently, a research group in the Idaho National Laboratory (INL) is testing reduced scale fuel rods and increased power density to accelerate evolution of fuel phenomena in metallic fuels. In support of nuclear fuel rod development, compressing fuel test process, and accelerating fuel phenomena, it was the purpose of this study to investigate nuclear fuel performance phenomena via literature review and effectively scale the initial conditions, boundary conditions, and geometric properties to describe the time-dependent response including to fuel burnup, thermal and mechanical stress, transmutations and inter-diffusion, and other relevant observed phenomena. The study was based on BISON simulations of historic EBR-II metallic fuel experiments and Dynamical System Scaling (DSS) method are utilized to assess effects of scaling activity including fuel phenomena acceleration and calculations of time-dependent distortions. The research successfully scaled metallic fuel phenomena, accelerated fuel testing, and assessed the distortions for each scaled case derived.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Monte Carlo Explicitly Correlated Second-Order Many-Body Green’s Function Calculations of Semiconductor Band Gaps

A systematically converging series of ab initio, post-density-functional, size-consistent, electron-correlated approximations is desired for predictive computing of felectronic band structures of insulating, semiconducting, and metallic solids. A series that meets all of these desiderata (except the applicability to metals) is ab initio many-body Green's function theory based on Gaussian-type-orbital (GTO) basis sets. Here, its leading-order approximation, the second-order Green's function (GF2) method in the diagonal and frequency-independent approximations with the aug-cc-pVDZ basis set, is applied to the fundamental band gaps of three semiconductors (diamond, silicon, and silicon carbide in the zincblende structure) using cluster models. Corrections are made to the basis-set-incompleteness errors by the explicit-correlation (F12) ansatz (GF2-F12) for the valence band edges. The crystals are modeled as surface-passivated clusters of increasing sizes, whose wave functions are expanded by up to 2709 GTO basis functions. Immense computational costs of these calculations are overcome by the highly scalable stochastic algorithm of the Monte Carlo GF2-F12 method, whose operation cost per state increases only as a cubic power of system size, which has a tiny memory footprint and easily achieves near-perfect parallel efficiency on thousands of CPUs or on hundreds of GPUs. The correlated, F12-corrected highest-occupied and lowest-unoccupied molecular-orbital energy (HOMO-LUMO) gap is 5.78 ± 0.07 eV for C 87 H 76 as compared with the experimental value of the fundamental (indirect) band gap of bulk diamond at 5.48 eV. The correlated, F12-corrected HOMO-LUMO gaps for Si 75 H 76 and Si 32 C 43 H 76 are 2.56 ± 0.15 eV and 3.50 ± 0.12 eV, respectively, which are expected to decrease further with increasing cluster sizes. As a result, the experimental fundamental (indirect) band gaps of bulk silicon and silicon carbide are 1.17 eV and 2.42 eV, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗