Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distributed Energy Program”

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 937 records · Page 52

Cosmic-Ray Extremely Distributed Observatory

The Cosmic-Ray Extremely Distributed Observatory (CREDO) is a newly formed, global collaboration dedicated to observing and studying cosmic rays (CR) and cosmic-ray ensembles (CRE): groups of at least two CR with a common primary interaction vertex or the same parent particle. The CREDO program embraces testing known CR and CRE scenarios, and preparing to observe unexpected physics, it is also suitable for multi-messenger and multi-mission applications. Perfectly matched to CREDO capabilities, CRE could be formed both within classical models (e.g., as products of photon–photon interactions), and exotic scenarios (e.g., as results of decay of Super-Heavy Dark Matter particles). Their fronts might be significantly extended in space and time, and they might include cosmic rays of energies spanning the whole cosmic-ray energy spectrum, with a footprint composed of at least two extensive air showers with correlated arrival directions and arrival times. As the CRE are predominantly expected to be spread over large areas and, due to the expected wide energy range of the contributing particles, such a CRE detection might only be feasible when using all available cosmic-ray infrastructure collectively, i.e., as a globally extended network of detectors. Thus, with this review article, the CREDO Collaboration invites the astroparticle physics community to actively join or to contribute to the research dedicated to CRE and, in particular, to pool together cosmic-ray data to support specific CRE detection strategies.

79 ASTRONOMY AND ASTROPHYSICS↗

Data Centers and Digital Assurance Workshop 3 – Mitigations for Digital Assurance Risks

The third session of the TADA (Technical Assistance for Digital Assurance) Data Centers Cohort, held on November 18, 2025, focused on developing mitigation strategies for digital assurance risks identified in previous workshops. Hosted by Idaho National Laboratory (INL) and ScottMadden, the session emphasized the application of Cyber-Informed Engineering (CIE) to data center infrastructure, particularly at the utility–data center interface. Participants revisited and ranked key digital assurance risks, including architecture and interface weaknesses, governance gaps, and AI-enabled threats. The workshop introduced the 12 principles of CIE, advocating for consequence-focused design, engineered controls, and secure information architecture to proactively reduce cyber-physical vulnerabilities. These principles were applied to critical data center systems such as power distribution, UPS, cooling, SCADA/BMS, and grid-forming batteries. The session also addressed governance challenges at the interconnection boundary, highlighting the need for clear roles in telemetry sharing, firmware management, and trip settings. Special attention was given to emerging risks from behind-the-meter (BTM) generation, including reverse-power flow and the integration of small modular reactors (SMRs), which shift data centers from large loads to complex generation nodes. Participants explored how interconnection agreements can serve as enforceable instruments for digital assurance, and reviewed gaps in current standards such as NERC CIP, IEC 62443, and IEEE 1547. The workshop concluded with pathways to standardization, including model agreement language, state-level programs, and expanded NERC guidance. INL also presented tools and frameworks for secure procurement and supplier risk management, reinforcing the need for integrated engineering and policy solutions to secure the evolving data center–grid ecosystem. Session 3 of 3.

24 - POWER TRANSMISSION AND DISTRIBUTION↗

Energy Northwest — Advanced Grid Interactive Load Efficiency (AGILE) : A Techno-economic Assessment

The Advanced Grid Interactive Load Efficiency project team, led by Energy Northwest, was awarded funding as part of the fourth round of the Clean Energy Fund program to investigate, co-create, and complete a preliminary design for grid-interactive efficient buildings (GEBs) for a number of schools served by Grays Harbor PUD. The project team has evaluated different combinations of systems at selected school buildings and the potential to respond to grid and customer needs, and develop a preliminary system technical design, including specifying desired GEB technology options, determining communication and electrical interface requirements, assessing control software sensor requirements, and interfacing with building management systems. PNNL led the techno-economic assessment efforts, leveraging its advanced modeling and analytical methods and tools developed for building-to-grid integration. This report documents the economic assessment for different GEB design options, encompassing an in-depth examination of use cases and value propositions, assumptions and inputs, models, methods, case studies, as well as key findings.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Partner Baseline Assessment Guide: Clean Energy for Low Income Communities

The U.S. Department of Energy Better Buildings Clean Energy for Low Income Communities Accelerator aims to lower energy bills in low to moderate income communities through expanded installation of energy efficiency and distributed renewables. Upon joining the Accelerator, partners committed to developing Action Plans after a year. As shown in the diagram below, baselining is a useful approach to informing your planning efforts, and should involve input from a broad range of stakeholders. The following guide was developed by DOE to assist partners with baselining, and identifying needs and gaps to address in action plans.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Nuclear Physics Network Requirements Review (Final Report)

The Energy Sciences Network (ESnet) is the high-performance network user facility for the US Department of Energy (DOE) Office of Science (SC) and delivers highly reliable data transport capabilities optimized for the requirements of data-intensive science. In essence, ESnet is the circulatory system that enables the DOE science mission by connecting all its laboratories and facilities in the US and abroad. ESnet is funded and stewarded by the Advanced Scientific Computing Research (ASCR) program and managed and operated by the Scientific Networking Division at Lawrence Berkeley National Laboratory (LBNL). ESnet is widely regarded as a global leader in the research and education networking community. ESnet interconnects DOE national laboratories, user facilities, and major experiments so that scientists can use remote instruments and computing resources as well as share data with collaborators, transfer large datasets, and access distributed data repositories. ESnet is specifically built to provide Between July 2023 and October 2023, ESnet and the Nuclear Physics program (NP) of the DOE SC organized an ESnet requirements review of NP-supported activities. Preparation for these events included identification of key stakeholders: program and facility management, research groups, and technology providers. Each stakeholder group was asked to prepare formal case study documents about its relationship to the NP program to build a complete understanding of the current, near-term, and long-term status, expectations, and processes that will support the science going forward.

97 MATHEMATICS AND COMPUTING↗

Requirements Description of GAMSOR and GAMSRC

This report presents the modeling and simulation capabilities of Argonne National Laboratory’s GAMSOR and GAMSRC code that is used in present reactor analysis activities. These capabilities will then be used to establish the set of verification tasks necessary to verify GAMSOR for use on commercial projects. A similar approach was taken for the PERSENT, REBUS and DIF3D software packages. The GAMSOR program is a modified version of DIF3D used to obtain the coupled neutron gamma flux distributions for a given reactor system. GAMSOR was first developed in the 1980s and its work has progressed in parallel with DIF3D throughout its various development phases. While GAMSOR is compatible with DIF3D-FD and DIF3D-Nodal, DIF3D-VARIANT is most commonly used today and will be the focus of the present verification effort. As covered in [1], GAMSOR was built to provide more accurate heating distributions by explicitly solving for both the gamma and neutron heating distributions. Without this consideration, gamma heat is normally assumed to be deposited locally at the site of gamma production. This is acceptable for modeling some phenomena because the total energy production is preserved, but it can lead to inaccurate approximations of fluence-dependent behaviors such as fuel performance, swelling, cladding embrittlement, radiation induced creep, and power peaking. The outputs from GAMSOR are contained in multiple intermediate and final output files: FIXSRC, VARSRC, GTFLUX, GHFLUX, PWDINT, NPDINT, and GPDINT. DIF3D normally produces a PWDINT file containing power density results, with the assumption that gamma rays are absorbed where they are emitted (i.e., energy deposition is local with no modeled transport). With GAMSOR, neutron power density (NPDINT) and gamma power density (GPDINT) are both produced, with their summed power densities replacing the results in PWDINT for total heating (i.e., recoverable) power density.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Accelerating Ocean-Based Renewable Energy Educational Opportunities to Achieve a Clean Energy Future

The United Nations has named 2021–2030 the Decade of Ocean Science for Sustainable Development with goals to 'strengthen the international cooperation needed to develop the scientific research and innovative technologies that can connect ocean science with the needs of society' (IOC 2019 The science we need for the ocean we want: the United Nations decade of ocean science for sustainable development (2021–2030) (Paris) p 24). Important actions that have been identified in support of sustainable development goals include capacity-building, training, and education. This includes educational opportunities for ocean-based renewable energy development in support of a healthy planet and ocean. Offshore wind is experiencing rapid development globally and, while the U.S. offshore wind market is still nascent, it is on the brink of exponential growth based on large cost reductions driven largely by European development and technology advances. Growth is also expected in wave and tidal energy with significant opportunities identified for various distributed markets through the Powering the Blue Economy™ initiative, with longer-term implications for expansion at the utility scale (LiVecchi et al 2019 Powering the blue economy; exploring opportunities for marine renewable energy in maritime markets p 207). In order to expedite progress and maximize benefits to the national, state, and local economies, these development actions will require a broad, diverse, and appropriately trained workforce. The ocean-based renewable energy workforce needs engineers and scientists to develop cost-effective technologies, as well as trade and maritime workers to eventually deploy the technologies at scale. In the United States, educational institutions, state governments, and private developers are taking action to understand job skills and capability requirements and to develop educational and training programs to meet offshore workforce needs; most are focused on offshore wind power, with gaining interest in marine energy. This article explores the workforce requirements of the growing ocean-based renewable energy industry and the current state of education and training programs to meet those requirements in order to identify gaps and make recommendations for further workforce development activities and initiatives. An international view needs to be adopted that incorporates the education and skill needs of early-stage marine energy technologies and evolving offshore wind technologies together with more market-ready offshore renewable energy markets. By accelerating educational development opportunities in ocean-based renewable energy, these growing blue economy markets can deliver significant economic and social benefits.

50 EE - Wind and Water Power Program - Water (EE-4↗

Enhancement of PyARC for Westinghouse Electric Company’s Lead Fast Reactor Design and Modeling (Final TCF Report)

Westinghouse Electric Company is a nuclear reactor vendor headquartered in the U.S. that is developing advanced reactor technology for the U.S. and global markets. Westinghouse has been relying on the neutronics Argonne Reactor Codes (ARC) executed through the NEAMS Workbench and its PyARC module that are developed under the DOE-NE Nuclear Energy Advanced Modeling and Simulation (NEAMS) and Advanced Reactor Technology (ART) – Fast Reactor programs. Through this user experience, Westinghouse identified several enhancements that would benefit the ARC codes’ usability by the US industry and therefore its commercialization potential. The enhancements were proposed to deliver both improvements in workflow and analysis capabilities to better support effective fast reactor core design and analysis to the nuclear industry. The PyARC workflow was extended in this project by integrating non-neutronic ARC codes DASSH and NUBOW-3D. The Ducted Assembly Steady-State Heat equation (DASSH) code is developed at ANL to perform steady-state thermal hydraulic sub-channel analysis in liquid metal fast reactor assemblies to determine optimized coolant flow and temperature distributions, which in this project was updated and validated for lead fast reactor (LFR) applications. The interface between REBUS and NUBOW-3D were improved in this project to assess the impact of the core restraint design and thermal induced expansion effects on the reactivity of the core, and to model the deformations of the fuel assemblies induced by temperature and irradiation. Finally, the ARC models that were extensively verified and validated through various SFR-based modeling benchmarks are extended in this project through code-to-code comparison on relevant LFR-specific neutronics benchmarks against Monte-Carlo neutronic solutions. Overall, this work enables verification of the capability of the ARC codes for a wide range of Generation-IV reactor designs. The outcome of this project is the release of a comprehensive modeling toolkit of validated, robust and efficient codes, as well as their user interface, that enables industry to perform a wide range of fast reactor analyses for design and licensing of their concepts.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Life Cycle Analysis of Refinery Fuel Products Air Emissions Using Refinery Emissions Inventory Data

With the transportation sector being a major source of air pollutant emissions, this study seeks to evaluate criteria pollutant emissions associated with the use of petroleum fuels in vehicles using the GREET® (Greenhouse gases, Regulated Emissions, and Energy use in Transportation) well-to-wheel model. The life cycle of petroleum fuels includes the major stages of crude recovery, crude transportation, crude refining, fuel distribution and fuel combustion (vehicle operation). This study focuses on quantifying refinery emissions, in terms of both refinery throughput and various refinery fuel products, by matching a refinery emissions inventory dataset with refinery operation information, at various levels (national, regional, facility, and process unit). The refinery emissions in terms refinery throughput (per refinery crude input, in g/mmbtu crude) are developed at a regional level for various Petroleum Administration for Defense Districts (PADD), and also aggregated at the national level. To allocate refinery emissions to refinery products, a methodology is developed using a refinery operation optimization Linear Programming (LP) model for individual refineries. The results from individual refinery facilities are then aggregated to obtain statistical national average (capacity weighted) emissions for each product. Because refineries vary by complexity of operation, configuration, crude quality, and product slate, the derived emissions demonstrate large variation for each product and per unit of crude input. Relative to previous estimates used in the GREET model, the present research shows reductions in refinery (and life cycle) emissions for most criteria pollutants. The reduction is especially significant for sulfur dioxide (SOx) emissions.

Allocation↗

Macrophysical properties of continental cumulus clouds from active and passive remote sensing

Cloud amount is an essential and extensively used macrophysical parameter of cumulus clouds. It is commonly defined as a cloud fraction (CF) from zenith-pointing ground-based active and passive remote sensing. However, conventional retrievals of CF from the remote sensing data with very narrow field-of-view (FOV) may not be representative of the surrounding area. Here we assess its representativeness using an integrated dataset collected at the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) program's Southern Great Plains (SGP) site in Oklahoma, USA. For our assessment with focus on selected days with single-layer cumulus clouds (2005-2016), we include the narrow-FOV ARM Active Remotely Sensed Clouds Locations (ARSCL) and large-FOV Total Sky Imager (TSI) cloud products, the 915-MHz Radar Wind Profiler (RWP) measurements of wind speed and direction, and also high-resolution satellite images from Landsat and the Moderate Resolution Imaging Spectroradiometer (MODIS). We demonstrate that a root-mean-square difference (RMSD) between the 15-min averaged ARSCL cloud fraction (CF) and the 15-min averaged TSI fractional sky cover (FSC) is large (up to 0.3). We also discuss how the horizontal distribution of clouds can modify the obtained large RMSD using a new uniformity metric. The latter utilizes the spatial distribution of the FSC over the 100° FOV TSI images obtained with high temporal resolution (30 sec sampling). We demonstrate that cases with more uniform spatial distribution of FSC show better agreement between the narrow-FOV CF and large-FOV FSC, reducing the RMSD by up to a factor of 2.

ARM Active Remotely Sensed Cloud Locations (ARSCL)↗

Exploring baryon resonances with transition generalized parton distributions: status and perspectives

QCD gives rise to a rich spectrum of excited baryon states. Understanding their internal structure is important for many areas of nuclear physics, such as nuclear forces, dense matter, and neutrino-nucleus interactions. Generalized parton distributions (GPDs) are an established tool for characterizing the QCD structure of the ground-state nucleon. They are used to create 3D tomographic images of the quark/gluon structure and quantify the mechanical properties such as the distribution of mass, angular momentum, and forces in the system. Transition GPDs extend these concepts to N → N* transitions and can be used to characterize the 3D structure and mechanical properties of baryon resonances. They can be probed in high-momentum-transfer exclusive electroproduction processes with resonance transitions e + N → e' + M + N*, such as deeply-virtual Compton scattering (M = γ) or meson production (M = π, K, etc.), and in related photon/hadron-induced processes. This White Paper describes a research program aiming to explore baryon resonance structure with transition GPDs. This includes the properties and interpretation of the transition GPDs, theoretical methods for structures and processes, first experimental results from JLab 12 GeV, future measurements with existing and planned facilities (JLab detector and energy upgrades, COMPASS/AMBER, EIC, EicC, J-PARC, LHC ultraperipheral collisions), and the theoretical and experimental developments needed to realize this program.

Experimental Nuclear Physics↗

Dynamic Probabilistic Risk Assessment Based Response Surface Approach for FLEX and Accident Tolerant Fuels for Medium Break LOCA Spectrum

After the Fukushima Daiichi Accident, the safety features such as accident tolerant fuel (ATF) and diverse and flexible coping strategies (FLEX) for existing nuclear fleets are being investigated by the US Department of Energy under the Light Water Reactor Sustainability Program. This research is being conducted to quantify the risk-benefit of these safety features. Dynamic probabilistic risk assessment (DPRA)-based response-surface approach has been presented to quantify the FLEX and ATF benefits by estimating the risk associated with each option. ATFs with multilayered silicon carbide (SiC), iron-chromium-aluminum, and chromium-coated zirconium cladding were considered in this study. While these ATF candidates perform better than the current zirconium cladding (Zr), they may introduce additional failure modes in some operating conditions. The fuel failure analysis modules (FAMs) were developed to investigate ATF performance. The dynamic risk assessments were performed using RAVEN, a DPRA tool, coupled with RELAP5 and FAMs. A cumulative distribution function-based index provided a mean of comparing the benefits of safety enhancements. For medium break loss of coolant accidents, FLEX operational timing window for each fuel type was estimated. Among these ATF candidates, SiC-type ATF was the most beneficial candidate for an increased safety margin than Zr-based fuel and was found to complement FLEX strategies in terms of risk and coping time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

HERSCHEL OBSERVATIONS OF DUSTY DEBRIS DISKS

We present results from two Herschel observing programs using the Photodetector Array Camera and Spectrometer. During three separate campaigns, we obtained Herschel data for 24 stars at 70, 100, and 160 μ m. We chose stars that were already known or suspected to have circumstellar dust based on excess infrared (IR) emission previously measured with the InfraRed Astronomical Satellite ( IRAS ) or Spitzer and used Herschel to examine long-wavelength properties of the dust. Fifteen stars were found to be uncontaminated by background sources and possess IR emission most likely due to a circumstellar debris disk. We analyzed the properties of these debris disks to better understand the physical mechanisms responsible for dust production and removal. Seven targets were spatially resolved in the Herschel images. Based on fits to their spectral energy distributions, nine disks appear to have two temperature components. Of these nine, in three cases, the warmer dust component is likely the result of a transient process rather than a steady-state collisional cascade. The dust belts at four stars are likely stirred by an unseen planet and merit further investigation.

79 ASTRONOMY AND ASTROPHYSICS↗

Global transpiration data from sap flow measurements: the SAPFLUXNET database

Plant transpiration links physiological responses of vegetation to water supply and demand with hydrological, energy, and carbon budgets at the land–atmosphere interface. However, despite being the main land evaporative flux at the global scale, transpiration and its response to environmental drivers are currently not well constrained by observations. Here we introduce the first global compilation of whole-plant transpiration data from sap flow measurements (SAPFLUXNET, https://sapfluxnet.creaf.cat/, last access: 8 June 2021). We harmonized and quality-controlled individual datasets supplied by contributors worldwide in a semi-automatic data workflow implemented in the R programming language. Datasets include sub-daily time series of sap flow and hydrometeorological drivers for one or more growing seasons, as well as metadata on the stand characteristics, plant attributes, and technical details of the measurements. SAPFLUXNET contains 202 globally distributed datasets with sap flow time series for 2714 plants, mostly trees, of 174 species. SAPFLUXNET has a broad bioclimatic coverage, with woodland/shrubland and temperate forest biomes especially well represented (80 % of the datasets). The measurements cover a wide variety of stand structural characteristics and plant sizes. The datasets encompass the period between 1995 and 2018, with 50 % of the datasets being at least 3 years long. Accompanying radiation and vapour pressure deficit data are available for most of the datasets, while on-site soil water content is available for 56 % of the datasets. Many datasets contain data for species that make up 90 % or more of the total stand basal area, allowing the estimation of stand transpiration in diverse ecological settings. SAPFLUXNET adds to existing plant trait datasets, ecosystem flux networks, and remote sensing products to help increase our understanding of plant water use, plant responses to drought, and ecohydrological processes. SAPFLUXNET version 0.1.5 is freely available from the Zenodo repository (https://doi.org/10.5281/zenodo.3971689; Poyatos et al., 2020a). The “sapfluxnetr” R package – designed to access, visualize, and process SAPFLUXNET data – is available from CRAN.

58 GEOSCIENCES↗

Optimization of Energy Storage System Economics and Controls by Incorporating Battery Degradation Costs in REopt

The use of stationary electrochemical energy storage systems utilizing lithium-ion batteries has increased rapidly as the production scale and price for lithium-ion batteries has decreased. These energy storage systems are crucial for maintaining grid resiliency, especially for grids operating with high penetration of renewable energy generation assets or for with a variety of distributed energy generation and storage systems. One challenging factor for the development of battery energy storage systems is estimating the proper sizing, in terms of both power and energy, that minimizes total costs over the lifetime of the systems; this calculation is difficult in simple cases, where a battery is costed independently, but is extremely challenging when building loads and electrical generation by photovoltaic resources are also considered. REopt is a techoeconomic optimization tool developed by NREL to address these challenges. Previously, battery degradation has been priced by simply assuming a 10-year replacement schedule for battery systems. However, this does not account for varying degradation trends observed across real-world batteries, or allow for batteries to be operated in a degradation-aware manner that optimizes battery dispatch based on operating costs. This work incorporates a battery life model into REopt. This battery life model is simple, so that it may be solvable within the constrains of a mixed-integer linear optimization problem, but is fit to accelerated aging data recorded in the lab. To achieve the best possible accuracy for lifetime estimates given these constraints, parameters for the battery life model in REopt are estimated by fitting 20-year simulations of battery life after identifying state-space battery degradation model from accelerated aging data. Comparisons of battery life predicted in REopt and from the state-space battery degradation model to ensure validity of lifetime estimates made by REopt. Battery life and cost is optimized by controlling three decision to minimize system life cost: battery sizing, daily state-of-charge, and daily energy-throughput. The cost of battery degradation as a function of these control variables is then estimated assuming two possible maintenance strategies: replacement, where the entire battery system is replaced if cell reach an end-of-life capacity threshold; and augmentation, which establishes a fund to pay for continual purchase of new batteries to maintain the initial energy capacity of the system. These two strategies offer conservative (for replacement) and optimistic (for augmentation) bounds for total system cost. The degradation cost incurred by these strategies is then used to control battery dispatch decisions, operating the battery in a degradation-aware manner that maximizes battery lifetime while also providing energy when favorable. Because the mixed-integer linear program has perfect foresight of future energy needs, batteries with degradation costs are always operated using 'just-in-time' charging, which is unrealistic, as no energy is left in the storage system to perform other energy services or to serve as emergency back-up power. To combat this, an inequality constraint on the average annual state-of-charge is imposed, and the sensitivity of system cost to average stored energy, e.g., the cost of system resiliency, can be quantified. Analysis of results has several conclusions, for instance, oversizing of battery storage systems is not a cost burden when battery storage is an optimal solution, as any additional battery capacity can simply be utilized to avoid costs of purchasing energy from a utility.

battery↗

Key insights from US Department of Energy Better Plants workforce development bootcamps (2022–2025)

This study examines the effectiveness of the US Department of Energy’s Better Plants Program Bootcamps, which are designed to enhance participants’ technical skills in improving energy efficiency and optimizing operations in manufacturing facilities. Through the analysis of survey data collected from 529 participants across 9 bootcamps, the research investigates the motivations, benefits, and demographic trends of attendees. The findings reveal that skill acquisition and improvement are primary drivers for participation, with key benefits including hands-on training on diagnostic equipment and software tools, networking opportunities, and access to technical resources. The analysis shows strong participation from sectors characterized by high energy consumption and employment, such as chemical and transportation equipment manufacturing. Over 50% of participants have job titles that include “EHS” or “Energy” showing their key roles in leading energy efficiency and energy management efforts in manufacturing. Furthermore, the analysis highlights the distribution of participants across managerial, engineering, and technical roles, revealing a higher representation of managers and engineers. This observation suggests a need for targeted outreach to engage technicians, equipment operators, maintenance staff, and floor workers to ensure comprehensive workforce development. The post-bootcamp survey showed that the participants highly valued the opportunities for peer learning and idea exchange, and the benefits they gained from them. This research contributes to the advancement of manufacturing education by demonstrating the efficacy of specialized training in addressing critical industry challenges and fostering a more competent and empowered workforce.

Energy efficiency↗

A Massively Parallel Implementation of the CCSD(T) Method Using the Resolution-of-the-Identity Approximation and a Hybrid Distributed/Shared Memory Parallelization Model

In this work, a parallel algorithm is described for the coupled-cluster singles and doubles method augmented with a perturbative correction for triple excitations [CCSD(T)] using the resolution-of-the-identity (RI) approximation for two-electron repulsion integrals (ERIs). The algorithm bypasses the storage of four-center ERIs by adopting an integral-direct strategy. The CCSD amplitude equations are given in a compact quasi-linear form by factorizing them in terms of amplitude-dressed three-center intermediates. A hybrid MPI/OpenMP parallelization scheme is employed, which uses the OpenMP-based shared memory model for intranode parallelization and the MPI-based distributed memory model for internode parallelization. Parallel efficiency has been optimized for all terms in the CCSD amplitude equations. Two different algorithms have been implemented for the rate-limiting terms in the CCSD amplitude equations that entail and -scaling computational costs, where N O and N V denote the number of correlated occupied and virtual orbitals, respectively. One of the algorithms assembles the four-center ERIs requiring N V 4 and N O 2 N V 2 -scaling memory costs in a distributed manner on a number of MPI ranks, while the other algorithm completely bypasses the assembling of quartic memory-scaling ERIs and thus largely reduces the memory demand. It is demonstrated that the former memory-expensive algorithm is faster on a few hundred cores, while the latter memory-economic algorithm shows a better strong scaling in the limit of a few thousand cores. The program is shown to exhibit a near-linear scaling, in particular for the compute-intensive triples correction step, on up to 8000 cores. The performance of the program is demonstrated via calculations involving molecules with 24–51 atoms and up to 1624 atomic basis functions. As the first application, the complete basis set (CBS) limit for the interaction energy of the π-stacked uracil dimer from the S66 data set has been investigated. This work reports the first calculation of the interaction energy at the CCSD(T)/aug-cc-pVQZ level without local orbital approximation. The CBS limit for the CCSD correlation contribution to the interaction energy was found to be -8.01 kcal/mol, which agrees very well with the value -7.99 kcal/mol reported by Schmitz, Hättig, and Tew [ Phys. Chem. Chem. Phys. 2014 , 16 , 22167-22178]. The CBS limit for the total interaction energy was estimated to be -9.64 kcal/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗