Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “close”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Hydrodynamics and Load Shedding Behavior of a Variable-Geometry Oscillating Surge Wave Energy Converter (OSWEC)

In order to improve their long-term viability, wave energy converters (WECs) need to be able to shed loads when a threshold wave condition is exceeded. As shown by Tom et al. (2016) [1], provision of adjustable flaps within the body of an oscillating surge wave energy converter (OSWEC) allows wave energy to pass through the device. A control system may then be able to open and close the flaps when waves approaching the device exceed preset thresholds. The variable-geometry OSWEC (VG-OSWEC) concept studied in this paper is a bottom-hinged, rectangular wave paddle with five flaps of elliptical cross-section embedded into the face of the paddle. System ID tests were conducted on this VG-OSWEC device at a 1:14 scale in a wave basin. Free decay tests showed that the damping was distinctly nonlinear when the flaps were fully open, and the natural frequency increased by 40% when compared with the flaps in a fully closed configuration. Tests with regular wave conditions were used to develop the response amplitude operator for the rotational motion about the hinge. These response amplitude operator results when compared with numerical simulations run using WEC-Sim/WAMIT and ANSYS AQWA, show strong agreement with the flap open and closed conditions. The regular-wave condition measurements also show that the wave excitation moment about the hinge was reduced by up to 60% when the flaps were fully open. The experiments serve to demonstrate the potential of the variable geometry design to shed loads and survive harsh ocean environments.

load shedding↗
DOE OSTI2022

Identifying the Gate-Opening Mechanism in the Flexible Metal–Organic Framework UTSA-300

Atomic-level understanding of the gate-opening phenomenon in flexible porous materials is an important step toward learning how to control, design, and engineer them for applications such as the separation of gases from complex mixtures. Here, we report such mechanistic insight through an in-depth study of the pressure-induced gate-opening phenomenon in our earlier reported metal–organic framework (MOF) Zn(dps) 2 (SiF 6 ) (dps = 4,4'-dipyridylsulfide), also called UTSA-300, using isotherm and calorimetry measurements, in situ infrared spectroscopy, and ab initio simulations. UTSA-300 is shown to selectively adsorb acetylene (C 2 H 2 ) over ethylene (C 2 H 4 ) and ethane (C 2 H 6 ) and undergoes an abrupt gate-opening phenomenon, making this framework a highly selective gas separator of this complex mixture. The selective adsorption is confirmed by pressure-dependent in situ infrared spectroscopy, which, for the first time, shows the presence of multiple C 2 H 2 species with varying strengths of bonding. A rare energetic feature at the gate-opening condition of the flexible MOF is observed in our differential heat energies, directly measured by calorimetry, showcasing the importance of this tool in adsorption property exploration of flexible frameworks and offering an energetic benchmark for further energy-based fundamental studies. Based on the agreement of this feature with ab initio-based adsorption energies of C 2 H 2 in the closed-pore structure UTSA-300a (“a” refers to the activated form), this feature is assigned to the weakening of the H-bond C–H···F formed between C 2 H 2 and fluorine of the MOF. Furthermore, our analysis identifies the weakening of this H-bond, the expansion of the closed-pore MOF upon successive C 2 H 2 coadsorption until its volume is close to that of the open-pore MOF, and the spontaneous gate opening to energetically favor C 2 H 2 adsorption in the open-pore structure as crucial steps in the gate-opening mechanism in this system.

36 MATERIALS SCIENCE↗
DOE OSTI2022

All-Atom Simulations Uncover Structural and Dynamical Properties of STING Proteins in the Membrane System

Recent studies have shown that the stimulator of interferon gene (STING) protein plays a central role in the immune system by facilitating the production of type I interferons in cells. The STING signaling pathway is also a prominent activator of cancer-killing T cells that initiate a powerful adaptive immune response. Since biomolecular signaling pathways are complicated and not easily identified through traditional experiments, molecular dynamics (MD) has often been used to study structural and dynamical responses of biological pathways. In this work, we carried out MD simulations for full-length chicken and human STING (chSTING and hSTING) proteins. Specifically, we investigated ligand-bound closed (holo) and ligand-unbound open (apo) forms of STING in the membrane system by comparing their conformational and dynamical differences. Our research provides clues for understanding the mechanism of the STING signaling pathway by uncovering detailed insights for the examined systems: the residues from each chain in the binding pocket are strongly correlated to one another in the open STING structure compared with those in the closed STING structure. Ligand-bound closed STING displays ~174° rotation of the ligand-binding domain (LBD) relative to the open STING structure. The dynamical analysis of residue Cys148 located in the linker region of hSTING does not support the earlier hypothesis that Cys148 can form disulfide bonds between adjacent STING dimers. We also demonstrate that using the full-length proteins is critical, since the MD simulations of the LBD portion alone cannot properly describe the global conformational properties of STING.

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

Directing Polymorphism of Colloid Crystals Using Conformational Entropy of Polymer Chains

Controlling the polytypes of close-packed structures of spherical colloids is still a challenging problem despite their wide occurrence. Here, in this work, we show that systematic engineering of the polytype structures of close-packed colloids is possible by using the conformational entropy of polymer chains confined in the interstitial space of colloid crystals. Our interstitial space analysis shows that the hexagonal close-packed (HCP) structures offer larger local interstitial space domains, and the structure director chains in favor of HCP counteract the entropic advantages of the face-centered cubic (FCC) lattices. Using model block copolymer colloids and the known lattice entropy of FCC, a proportionality parameter, β CP = 1.91 × 10 –3 ± 3.67 × 10 –4 , for quantifying the conformational entropy contribution toward HCP structures is extracted. This work demonstrates that the interstitial space of colloid crystals serves as a new structure engineering tool for the self-assembly of colloids.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A Polycrystalline Pd Surface Studied by Two-Dimensional Surface Optical Reflectance during CO Oxidation: Bridging the Materials Gap

Industrial catalysts are complex materials systems operating in harsh environments. The active parts of the catalysts are nanoparticles that expose different facets with different surface orientations at which the catalytic reactions occur. However, these facets are close to impossible to study in detail under industrially relevant operating conditions. Instead, simpler model systems, such as single crystals with a well-defined surface orientation, have been successfully used to study gas–surface interactions such as adsorption and desorption, surface oxidation, and oxidation/reduction reactions. To more closely mimic the many facets exhibited by nanoparticles and thereby close the so-called materials gap, there has also been a recent move toward using polycrystalline surfaces and curved crystals. However, these studies are limited either by the pressure or spatial resolution at realistic pressures or by the number of surfaces studied simultaneously. In this work, we demonstrate the use of reflectance microscopy to study a vast number of catalytically active surfaces simultaneously under realistic and identical reaction conditions. As a proof of concept, we have conducted an operando experiment to study CO oxidation over a Pd polycrystal, where the polycrystalline surface acts as a collection of many single-crystal surfaces. Finally, we visualized the resulting data by plotting the reflectivity as a function of surface orientation. We think the techniques and visualization methods introduced in this work will be key toward bridging the materials gap in catalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Formic Acid Electrooxidation on Pt or Pd Monolayer on Transition-Metal Single Crystals: A First-Principles Structure Sensitivity Analysis

Here, we present a density functional theory analysis of trends for the electrooxidation of formic acid (FAO) on a single Pt or Pd monolayer supported on the close-packed (111) facet of transition metals (Pt*/M or Pd*/M): Au, Ag, Cu, Pt, Pd, Ir, and Rh, the close-packed (0001) facet of Os, Ru, and Re, and the open (100) facet of Au, Ag, Pt, Pd, Ir, and Rh. We show that the deposition of Pt or Pd pseudomorphic monolayers on these single crystals modifies the electronic structure of the Pt or Pd monolayer. Specifically, we found a direct correlation between the d-band center of the Pt and Pd monolayer and the free energy of adsorbed CO*, the latter being a reactivity descriptor for FAO. Together with the free energy of adsorbed OH* as a second reactivity descriptor, we depict the thermochemistry of the reaction network as phase diagrams showing calculated free energies across regions of rate-determining steps. We found that FAO is structure-sensitive on most surfaces studied. Pt*/Au(111) is predicted to be the most active among all Pt*/M(111/0001) surfaces studied, despite binding CO*, the strongest among the close-packed facets. This is the case because of its superiority in activating water to OH*, thus removing CO* at lower potentials than other surfaces. On similar grounds, Pt*/Pd(100), Pd*/Re(0001), Pd*/Au(111), Pd*/Ag(111), and Pd*/Pt(111) are predicted to show higher FAO activity than the corresponding monometallic Pt and Pd surfaces.

(100) facets↗
DOE OSTI2023

Solvent Dynamics Are Critical to Understanding Carbon Dioxide Dissolution and Hydration in Water

Simulations of carbon dioxide (CO 2 ) in water may aid in understanding the impact of its accumulation in aquatic environments and help advance technologies for carbon capture and utilization (via, e.g., mineralization). Quantum mechanical (QM) simulations based on static molecular models with polarizable continuum solvation poorly reproduce the energetics of CO 2 hydration to form carbonic acid in water, independent of the level of QM theory employed. Only with density-functional-theory-based molecular dynamics and rare-event sampling, followed by energy corrections based on embedded correlated wavefunction theory (in conjunction with density functional embedding theory), can a close agreement between theory and experiment be achieved. Such multilevel simulations can serve as benchmarks for simpler, less costly models, giving insight into potential errors of the latter. The strong influence of sampling/averaging over dynamical solvent configurations on the energetics stems from the difference in polarity of both the transition state and product (both polar) versus the reactant (nonpolar). When a solute undergoes a change in polarity during reaction, affecting its interaction with the solvent, careful assessment of the energetic contribution of the solvent response to this change is critical. Here, we show that static models (without structural sampling) that incorporate three explicit water molecules can yield far superior results than models with more explicit water molecules because fewer water molecules yield less configurational artifacts. Static models intelligently incorporating both explicit (molecules directly participating in the reaction) and implicit solvation, along with a proper QM theory, e.g., CCSD(T) for closed-shell systems, can close the accuracy gap between static and dynamic models.

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

Mechanochemically Enabled Metastable Niobium Tungsten Oxides

Metastable compounds have greatly expanded the synthesizable compositions of solid-state materials and have attracted enormous amounts of attention in recent years. Especially, mechanochemically enabled metastable materials synthesis has been very successful in realizing cation-disordered materials with highly simple crystal structures, such as rock salts. Application of the same strategy for other structural types, especially for non-close-packed structures, is peculiarly underexplored. Niobium tungsten oxides (NbWOs), a class of materials that have been under the spotlight because of their diverse structural varieties and promising electrochemical and thermoelectric properties, are ideally suited to fill such a knowledge gap. In this work, we develop a new series of metastable NbWOs and realize one with a fully cation-disordered structure. Furthermore, we find that metastable NbWOs transform to a cation-disordered cubic structure when applied as a Li-ion battery anode, highlighting an intriguing non-close-packed–close-packed conversion process, as evidenced in various physicochemical characterizations, in terms of diffraction, electronic, and vibrational structures. Lastly, by comparing the cation-disordered NbWO with other trending cation-disordered oxides, we raise a few key structural features for cation disorder and suggest a few possible research opportunities for this field.

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

Autonomous Nanoparticle Synthesis Guided by In Situ Multiscale Structural Characterization

Autonomous synthesis platforms promise rapid exploration of vast parameter spaces; yet, integrating in situ structural characterization in closed-loop synthesis optimization remains challenging. We demonstrate a realization of such a closed-loop platform coupled with a droplet-flow microreactor, in situ X-ray scattering methods (SAXS/WAXS), and Gaussian process optimization to synthesize citrate-reduced Au nanoparticles with targeted characteristics. The system efficiently explored ∼19,000 synthesis recipes through 365 experiments, achieving precise control over size (4–60 nm) and polydispersity (σ < 0.11) across large citrate/gold ratios, exceeding traditional synthesis boundaries (1–10). Beyond confirming classical Turkevich–Frens trends, partial-dependence analysis revealed strong nonlinear coupling among precursor, citrate, and pH effects. Combining quantitative SAXS/WAXS analysis with electron microscopy characterization, we uncovered that crystallite size (d c ) and particle size (d) follow d c = 0.18d + β, where synthesis chemistry controls the intercept β while maintaining a universal slope. This parallel-band structure enables independent tuning of crystallite domain size at fixed particle diameter through a combination of chloride, gold precursor, citrate, and pH contributions (cross-validated Spearman ρ = 0.7 ± 0.1). High-resolution electron microscopy shows multiple lattice-fringe orientations within single particles, directly confirming polycrystalline domains and the ability to tune d c at the fixed d. The platform’s validation includes indistinguishable static versus flowing measurements, stable droplet transport at 100 °C, and <5% run-to-run variation, establishing a robust framework for mapping and controlling multiscale nanoparticle structure across expansive chemical spaces. In conclusion, the developed closed-loop platform can be applied to a borad range of nanosyntheis processes.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Adsorbate-Induced Adatom Formation on Lithium, Iron, Cobalt, Ruthenium, and Rhenium Surfaces

Recent experimental and theoretical studies have demonstrated the reaction-driven metal–metal bond breaking in metal catalytic surfaces even under relatively mild conditions. Here, we construct a density functional theory (DFT) database for the adsorbate-induced adatom formation energy on the close-packed facets of three hexagonal close-packed metals (Co, Ru, and Re) and two body-centered cubic metals (Li and Fe), where the source of the ejected metal atom is either a step edge or a close-packed surface. For Co and Ru, we also considered their metastable face-centered cubic structures. We studied 18 different adsorbates relevant to catalytic processes and predicted noticeably easier adatom formation on Li and Fe compared to the other three metals. The NH 3 - and CO-induced adatom formation on Fe(110) is possible at room temperature, a result relevant to NH 3 synthesis and Fischer-Tropsch synthesis, respectively. There also exist other systems with favorable adsorbate effects for adatom formation relevant to catalytic processes at elevated temperatures (500–700 K). Our results offer insight into the reaction-driven formation of metal clusters, which could play the role of active sites in reactions catalyzed by Li, Fe, Co, Ru, and Re catalysts.

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

Evolution of Droplet Size Distributions During the Transition of an Ultraclean Stratocumulus Cloud System to Open Cell Structure: An LES Investigation Using Lagrangian Microphysics

Here, a state-of-the-art Lagrangian microphysics scheme is used in a large-eddy simulation to investigate the stratocumulus transition from closed to open cell structure. Processes controlling precipitation development, which is a key to the transition, are analyzed by leveraging unique benefits of Lagrangian microphysics, particularly the ability to track computational drops in the flow. Sufficient time is needed for coalescence growth of cloud drops to drizzle within the updraft-downdraft cycle of large eddies. This favors broad drop size distributions (DSDs) and drizzle growth in downdrafts, where drops are typically much older than in updrafts. During the closed cell stage, mean cloud drop radius is too small, and the DSDs are too narrow, so that the timescale for coalescence is much longer than the large eddy turnover time and drizzle growth is limited. The closed-to-open cell transition occurs when these timescales become comparable and the precipitation flux increases sharply.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Model-based evaluation of school- and non-school-related measures to control the COVID-19 pandemic

The role of school-based contacts in the epidemiology of SARS-CoV-2 is incompletely understood. We use an age-structured transmission model fitted to age-specific seroprevalence and hospital admission data to assess the effects of school-based measures at different time points during the COVID-19 pandemic in the Netherlands. Our analyses suggest that the impact of measures reducing school-based contacts depends on the remaining opportunities to reduce non-school-based contacts. If opportunities to reduce the effective reproduction number (R e ) with non-school-based measures are exhausted or undesired and R e is still close to 1, the additional benefit of school-based measures may be considerable, particularly among older school children. As two examples, we demonstrate that keeping schools closed after the summer holidays in 2020, in the absence of other measures, would not have prevented the second pandemic wave in autumn 2020 but closing schools in November 2020 could have reduced R e below 1, with unchanged non-school-based contacts.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Neutralizing antibodies induced in immunized macaques recognize the CD4-binding site on an occluded-open HIV-1 envelope trimer

Broadly-neutralizing antibodies (bNAbs) against HIV-1 Env can protect from infection. We characterize Ab1303 and Ab1573, heterologously-neutralizing CD4-binding site (CD4bs) antibodies, isolated from sequentially-immunized macaques. Ab1303/Ab1573 binding is observed only when Env trimers are not constrained in the closed, prefusion conformation. Fab-Env cryo-EM structures show that both antibodies recognize the CD4bs on Env trimer with an ‘occluded-open’ conformation between closed, as targeted by bNAbs, and fully-open, as recognized by CD4. The occluded-open Env trimer conformation includes outwardly-rotated gp120 subunits, but unlike CD4-bound Envs, does not exhibit V1V2 displacement, 4-stranded gp120 bridging sheet, or co-receptor binding site exposure. Inter-protomer distances within trimers measured by double electron-electron resonance spectroscopy suggest an equilibrium between occluded-open and closed Env conformations, consistent with Ab1303/Ab1573 binding stabilizing an existing conformation. Studies of Ab1303/Ab1573 demonstrate that CD4bs neutralizing antibodies that bind open Env trimers can be raised by immunization, thereby informing immunogen design and antibody therapeutic efforts.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

An autoinhibited state of 53BP1 revealed by small molecule antagonists and protein engineering

The recruitment of 53BP1 to chromatin, mediated by its recognition of histone H4 dimethylated at lysine 20 (H4K20me2), is important for DNA double-strand break repair. Using a series of small molecule antagonists, we demonstrate a conformational equilibrium between an open and a pre-existing lowly populated closed state of 53BP1 in which the H4K20me2 binding surface is buried at the interface between two interacting 53BP1 molecules. In cells, these antagonists inhibit the chromatin recruitment of wild type 53BP1, but do not affect 53BP1 variants unable to access the closed conformation despite preservation of the H4K20me2 binding site. Thus, this inhibition operates by shifting the conformational equilibrium toward the closed state. Our work therefore identifies an auto-associated form of 53BP1—autoinhibited for chromatin binding—that can be stabilized by small molecule ligands encapsulated between two 53BP1 protomers. Such ligands are valuable research tools to study the function of 53BP1 and have the potential to facilitate the development of new drugs for cancer therapy.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Prediction of solubility parameters of lignin and ionic liquids using multi-resolution simulation approaches

The solubility parameter (SP) of a molecular species is a vital feature that indicates polarity and quantifies the ‘like-seeks-like’ principle, which is used in chemistry to screen solvents for dissolution. Recent studies have demonstrated that ionic liquids (ILs) and deep eutectic solvents (DESs) efficiently solubilize lignocellulosic biomass and promote enzymatic saccharification into sugars used for the production of biofuels and value-added chemicals. Understanding the solubility of plant biopolymers, particularly lignin, in ILs and DESs is critical for selecting candidate ILs and DESs for biomass pretreatment; however, experimentally measuring SPs is challenging. Thus, the present study investigates lignin dissolution mechanisms in IL/DES and prediction of the solubility parameters (Hildebrand and Hansen) of lignin, ILs, and DESs using multi-resolution simulation approaches. Here, solubility parameters of the studied compounds were predicted using molecular dynamics (MD) simulations, and the SP of lignin was determined to be 23–27 MPa 1/2 , which was close to the polymeric lignin solubility parameter (24.3–25.5 MPa 1/2 ). The SPs of ILs namely [Ch][Lys], [Ch][Oct], and [Emim][Lys] were predicted to be ~26 MPa 1/2 , which is close to lignin's SPs and resulted in increased biomass delignification. The MD simulated SPs were validated by both the COSMO-RS model and experimental investigations, with the results showing a close agreement between the predicted and experimentally obtained SPs. In addition, the enthalpy of vaporization (ΔH vap ) of ILs/DESs was predicted based on the potential energy of the system, and the ΔH vap of ILs/DESs was around 40–65 kcal mol –1 , which is 5–8 times higher than that of traditional organic solvents.

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

Radial velocity variability and the evolution of hot subdwarf stars

Hot subdwarf stars represent a late and peculiar stage in the evolution of low-mass stars, since they are likely formed by close binary interactions. In this work, we perform a radial velocity (RV) variability study of a sample of 646 hot subdwarfs with multi-epoch radial velocities based on spectra from Sloan Digital Sky Survey (SDSS) and Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAMOST). The atmospheric parameters and RVs were taken from the literature. For stars with archival spectra but without literature values, we determined the parameters by fitting model atmospheres. In addition, we redetermined the atmospheric parameters and RVs for all the He-enriched sdO/Bs. This broad sample allowed us to study RV-variability as a function of the location in theT eff - log g- and T eff - log n (He)/n (H) diagrams in a statistically significant way. We used the fraction of RV-variable stars and the distribution of the maximum RV variations ΔRV max as diagnostics. Both indicators turned out to be quite inhomogeneous across the studied parameter ranges. A striking feature is the completely dissimilar behaviour of He-poor and He-rich hot subdwarfs. While the former have a high fraction of close binaries, almost no significant RV variations could be detected for the latter. This has led us to the conclusion that there is likely no evolutionary connection between these subtypes. On the other hand, intermediate He-rich- and extreme He-rich sdOB/Os are more likely to be related. Furthermore, we conclude that the vast majority of this population is formed via one or several binary merger channels. Hot subdwarfs with temperatures cooler than ~24 000 K tend to show fewer and smaller RV-variations. These objects might constitute a new subpopulation of binaries with longer periods and late-type or compact companions. The RV-variability properties of the extreme horizontal branch (EHB) and corresponding post-EHB populations of the He-poor hot subdwarfs match and confirm the predicted evolutionary connection between them. Stars found below the canonical EHB at somewhat higher surface gravities show large RV variations and a high RV variability fraction. These properties are consistent with most of them being low-mass EHB stars or progenitors of low-mass helium white dwarfs in close binaries.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Fossil group origins

Aims: We analyse the large-scale structure out to 100 Mpc around a sample of 16 confirmed fossil systems using spectroscopic information from the Sloan Digital Sky Survey Data Release 16. Methods: We computed the distance between our fossil groups (FGs) and the centres of filaments and nodes from the literature. We also studied the density of bright galaxies, since this parameter is thought to be a good mass tracers, as well as the projected over-densities of galaxies. Finally, we applied a friends-of-friends (FoF) algorithm to detect virialised structures around our FGs and obtain an estimate of the mass available in their surroundings. Results: We find that FGs are mainly located close to filaments, with a mean distance of 3.7 ± 1.1 R 200 and a minimum distance of 0.05 R 200 . On the other hand, none of our FGs were found close to intersections, with a mean and minimum distance of 19.3 ± 3.6 and 6.1 R 200 , respectively. There is a correlation that indicates FGs at higher redshifts are found in denser regions, when we use bright galaxies as tracers of the mass. At the same time, FGs with the largest magnitude gaps (Δm 12 > 2.5) are found in less dense environments and tend to host (on average) smaller central galaxies. Conclusions: Our results suggest that FGs formed in a peculiar position within the cosmic web, close to filaments and far from nodes, whereby their interaction with the cosmic web itself may be limited. We deduce that FGs with brightest central galaxies (BCGs) that are relatively faint, high values of Δm 12 , and low redshifts could, in fact, be systems that are at the very last stage of their evolution. Moreover, we confirm theoretical predictions that systems with the largest magnitude gap are not massive.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Explicitly correlated coupled cluster method for accurate treatment of open-shell molecules with hundreds of atoms

We present a near-linear scaling formulation of the explicitly correlated coupled-cluster singles and doubles with the perturbative triples method [CCSD(T)F12¯] for high-spin states of open-shell species. The approach is based on the conventional open-shell CCSD formalism [M. Saitow et al., J. Chem. Phys. 146, 164105 (2017)] utilizing the domain local pair-natural orbitals (DLPNO) framework. The use of spin-independent set of pair-natural orbitals ensures exact agreement with the closed-shell formalism reported previously, with only marginally impact on the cost (e.g., the open-shell formalism is only 1.5 times slower than the closed-shell counterpart for the C160H322 n-alkane, with the measured size complexity of ≈1.2). Evaluation of coupled-cluster energies near the complete-basis-set (CBS) limit for open-shell systems with more than 550 atoms and 5000 basis functions is feasible on a single multi-core computer in less than 3 days. The aug-cc-pVTZ DLPNO-CCSD(T)F12¯ contribution to the heat of formation for the 50 largest molecules among the 348 core combustion species benchmark set [J. Klippenstein et al., J. Phys. Chem. A 121, 6580–6602 (2017)] had root-mean-square deviation (RMSD) from the extrapolated CBS CCSD(T) reference values of 0.3 kcal/mol. For a more challenging set of 50 reactions involving small closed- and open-shell molecules [G. Knizia et al., J. Chem. Phys. 130, 054104 (2009)], the aug-cc-pVQ(+d)Z DLPNO-CCSD(T)F12¯ yielded a RMSD of ∼0.4 kcal/mol with respect to the CBS CCSD(T) estimate.

Kumar, Ashutosh (ORCID:0000000175896030)↗