Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Complex chemical systems”

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

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

At least 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Methods and apparatus for gnotobiotic plant growth

The disclosure relates to a plant growth apparatus and related system to grow and maintain plants under controlled biotic conditions, for example to grow axenic (microbe-free) plants, gnotobiotic (defined microbiota) plants, and holoxenic (complex, or undefined microbiota) plants. This system allows aseptic bottom irrigation with water, soluble nutrients, chemicals, and/or microbiota. Plants can be inverted for dipping and/or vacuum infiltration. The system also allows for passive (gravity) drainage, thereby allowing for gas exchange and preventing root anoxia. A variety of plant growth substrates can be used within the plant growth apparatus as a plant growth medium. The plant growth apparatus, containing the growth substrate medium, can be completely flushed via the drainage port to remove potential toxic byproducts of the sterilization processes. The entire system is suitably constructed using autoclavable material.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

UO2 microstructural evolutions induced by Ni, Mo, and W dopants for intentional forensics

The concept of tagging nuclear fuel with a chemical barcode to enable forensics analysis across the nuclear fuel cycle is an area of active investigation, particularly to ensure fabrication viability without disrupting current fuel performance. This study explored the feasibility of using Ni, Mo, and W isotopic double-spikes as dopants in UO2 fuel from the perspective of fuel fabrication. Doped UO2 pellets were produced using conventional fuel fabrication processes, including powder mixing, sieving, pressing, and sintering in a reductive atmosphere. Two composition levels, 100 and 1000 ppm, were evaluated for each dopant element with isotopic double-spike configurations. For the Ni system, additional dopant concentrations of 250 and 500 ppm were produced with nonperturbed isotopic ratios. The results demonstrated that successful incorporation of Ni, Mo, and W double-spikes into UO2 pellets occurred with minimal shift in final density or dopant loss during pellet fabrication. Isotopic analysis confirmed the presence of the double-spike signature even when diluted with natural isotopic material in ratio of 1:5 in the fabrication process. Microstructural examinations revealed different impacts on grain size compared with undoped UO2. This study showed that Ni incorporation up to ∼500 ppm promoted moderate grain growth, whereas the Mo and W systems caused grain size reduction at all concentrations. Changes in the UO2 lattice parameter as a function of composition were detected exclusively for Ni up to 500 ppm, indicating that the Ni solid solution was the main factor for the observed grain growth. Insoluble (Mo and W) or supersaturated (Ni > 500 ppm) conditions produced grain size reduction. The Ni-doped pellets in the solution range resulted in a final microstructure within fuel specifications, demonstrating its potential benefits of employing complex dopant systems for potential nuclear forensic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Modeling aqueous association constants and mineral solubilities at subcritical and supercritical temperatures

The need for sustainable power generation has increased interest in the use of hydrothermal fluids for industrial applications. New high-enthalpy geothermal systems and biowaste-to-fuel processes are two relevant examples that employ supercritical fluids which require an in-depth understanding of complex chemical reactions occurring near the supercritical temperature of water (374 °C). As these processes operate in thermodynamic regimes that are not currently covered by a standard molar Gibbs energy of formation model, only empirical fits for single reaction systems are available which limit the use of multi-component phase equilibria calculations that are standard practice for less extreme environments. Here, in this paper, we advance a standard molar Gibbs energy of formation model able to operate in these otherwise inaccessible thermodynamic states to include species needed for key mineral solubility systems and ion association reactions. This work extends a model based on molecular statistical thermodynamics (MST) into four new systems (Na 3 PO 4 -H 2 O, LiOH-H 2 O, KOH-H 2 O, and BaSO 4 -H 2 O) by extending the model to cover 10 new species. For each of these systems, model predictions were consistently within the experimental uncertainties for the new systems covered. A breakdown of MST contributions to the model revealed that electrostatic and hard sphere contributions were key to reproducing density dependencies of standard molar Gibbs energy of formation values around the critical point of water.

42 ENGINEERING↗
DOE OSTI2020

Amphiphile Organization in Organic Solutions: An Alternative Explanation for Small-Angle X-ray Scattering Features in Malonamide/Alkane Mixtures

In this work, the role of different intermolecular interactions in the aggregation of amphiphiles in an organic solvent is studied for systems of relevance to liquid-liquid extraction (LLE), a chemical process used to selectively recover metals from complex mixtures. Of specific interest is the role, or lack thereof, of hydrogen bonding, which is often assumed to be a main driver of the organic phase structural organization that has been linked to separation efficacy. Toward that end, a series of malonamide extractants in n-dodecane have been studied in the absence of any extracted aqueous solutes, including water. The series of extractants includes N,N'-dimethyl-N,N'-dibutyltetradecylmalonamide (DMDBTDMA), two of its homologs, and N,N'-dimethyl-N,N'-dioctylhexylethoxymalonamide (DMDOHEMA). This simplified model LLE system enables systematic investigation of the role of dipole-dipole and alkyl tail steric interactions in amphiphile aggregation. Small-angle X-ray scattering (SAXS) profiles computed from molecular dynamics trajectories are in good agreement with experimental SAXS data. Molecular dynamics simulations show that malonamide aggregation results from dipole-driven self-association and lacks characteristic aggregate sizes. Mid-q correlation peaks in the SAXS profiles emerge at high concentration for each malonamide. In those densely packed solutions, the correlation peaks are observed to result from alkyl tail-induced spacing between electron-rich polar head groups, with peak positions determined by the different alkyl tail lengths present in the malonamide molecule. This explanation of the SAXS correlation peaks contrasts with the prevailing literature, which attributes mesoscale features observed in small-angle scattering to the formation of microemulsions. Instead, this work finds that these features are present in the absence of water or any reverse micellar organization of the malonamides. As such, molecular-scale malonamide self-association and packing, rather than microemulsion-based colloidal-scale descriptions, is a more appropriate framework for these LLE systems.

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

Electrochemical Characterization of Photo-Driven Hole-Scavenging by Cadmium Sulfide Quantum Dot-Nitrogenase Biohybrid Complexes

Light-driven biohybrid systems that couple semiconductor nanocrystals with enzymes offer a promising strategy for solar-to-chemical energy-conserving reduction reactions, yet are often limited by inefficient hole scavenging. Hole scavenging is critical for maintaining charge separation in the light-absorbing electron donor molecule and, thus, for sustaining catalytic turnover, making quantification essential to improving system efficiency. Here, we introduce a photoelectrochemical approach for real-time monitoring of hole scavenging in cadmium sulfide (CdS) quantum dot (QD)-nitrogenase MoFe protein biohybrids that catalyze the reduction of H+ and N2. Using hydroquinone (HQ) as the sacrificial electron donor (SED), oxidation of HQ to benzoquinone (BQ) by photogenerated valence-band holes is coupled to the electrochemical reduction of BQ at an electrode, enabling quantification of the hole-scavenging reaction via chronoamperometry. This approach provides a continuous, real-time readout of charge-transfer dynamics under illumination. Systematic variation of light intensity and donor concentration reveals that hole scavenging scales with photon flux and exhibits a non-linear dependence on SED concentration, with an optimum at intermediate HQ concentrations. These results establish chronoamperometry as a powerful tool for probing photoinduced charge transfer, advancing the quantitative understanding of hole-scavenging in biohybrid systems, and offering a generalizable framework for optimizing solar-driven biocatalysis.

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

Criticality Alarm Systems at the Idaho Chemical Processing Plant

The Idaho Chemical Processing Plant (ICPP) has a number of criticality alarm systems (CAS) in the various facilities. The complexity of the plant creates some unique CAS problems. A detailed review of the CAS has been made to determine the level of compliance with ANSI 8.3-1986. A major CAS upgrade project is currently in the design stage which will expand coverage and bring all CAS monitoring unshielded areas into full compliance with ANSI 8.3-1986. Because DOE Orders do not contain an exclusion that exempts the need for CAS coverage in shielded areas as provided in ANSI 8.3-1986, a future CAS upgrade project is in the planning stage which will also upgrade the CAS coverage in shielded portions of the plant.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Kinetics of amorphous defect phases measured through ultrafast nanocalorimetry

Recognition of the role of extended defects on local phase transitions has led to the conceptualization of the defect phase, localized thermodynamically stable interfacial states that have since been applied in a myriad of material systems to realize significant enhancements in material properties. Here, we explore the kinetics of grain boundary confined amorphous defect phases, utilizing the high temperature and scanning rates afforded by ultrafast differential scanning calorimetry to apply targeted annealing/quenching treatments at high rates capable of capturing the kinetic behavior. Four Al-based nanocrystalline alloys, including two binary systems, Al–Ni and Al–Y, and two ternary systems, Al–Mg–Y and Al–Ni–Y, are selected to probe the materials design space (enthalpy of mixing, enthalpy of segregation, chemical complexity) for amorphous defect phase formation and stability, with correlative transmission electron microscopy applied to link phase evolution and grain stability to nanocalorimetry signatures. A series of targeted isothermal annealing heat treatments is utilized to construct a Time–Temperature-Transformation curve for the Al–Ni system, from which a critical cooling rate of 2400 °C/s was determined for the grain boundary confined disordered-to-ordered transition. Finally, a thermal profile consisting of 1000 repeated annealing sequences was created to quantify the recovery of the amorphous defect phase following sequential annealing treatments, with results indicating remarkable microstructural stability after annealing at temperatures above 90% of the melting temperature. This work contributes to a deeper understanding of grain boundary localized thermodynamics and kinetics, with potential implications for the design and optimization of advanced materials with enhanced stability and performance.

Amorphous defect phases↗
DOE OSTI2024

Interpretable, extensible linear and symbolic regression models for charge density prediction using a hierarchy of many-body correlation descriptors

Here, density functional theory (DFT) is routinely used to make electronic structure predictions for high-throughput screening of materials and molecules for technologically relevant areas, like the identification of better catalysts, electronic materials, and drug discovery. However, the DFT formalism is limited by (a) its poor (quadratic-to-quartic) scaling, and (b) the need to perform repeated eigenvalue computations of the electronic Hamiltonian as part of its self-consistent field (SCF) iteration procedure to obtain the converged ground state electron density, ρ (r). Approaches that directly predict ρ (r) of a structure with high accuracy can accelerate conventional SCF calculations and can also be used in linearly scaling methods such as orbital-free DFT. To this end, we present a procedure to predict the ground state electron density of molecular and periodic three-dimensional systems directly from the atomic structure with a particular emphasis on physical interpretability. In our framework, ρ (r) is modeled using many-body correlation descriptors that accurately capture the effects of local atomic arrangements in the neighborhood of a grid point. Our use of a linear regression scheme to fit to charge density data enables transparent analysis of the relative contributions of various types of local atomic correlations. By systematically including increasingly complex correlations, our model is shown to accurately predict ρ (r) for a variety of chemically and electronically diverse systems — amorphous Ge, Al(001) slab, crystalline Ga 2 O 3 , molecular benzene, and polyethylene. We then demonstrate a symbolic regression-based protocol to construct easily computable, interpretable features from lower-order correlations that significantly improves our electron density predictions with effectively no increase in the computational cost.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Raman Spectroscopic In Situ Monitoring of Highly Turbid Media

The ability to run chemical processing more efficiently and cost effectively is a need that spans critical materials recovery and legacy nuclear waste cleanup. Sensors integrated to provide online monitoring are essential to addressing this need by providing near-real time feedback on process conditions, which can improve efficiency, aid in decision making, and reduce the need for grab sample measurements. Optical spectroscopy is well-suited for providing online chemical composition information and has been widely applied in varied chemical systems. However, applications in turbid matrices continue to represent substantial challenges to sensor performance, where absorption or scattering of excitation light can cause significant signal interference. Here, in this study, close-focus Raman probes are investigated for use in turbid media as a way to overcome the signal loss from the scattering of the Raman excitation source. This, paired with advanced data science techniques, allowed for the development of chemometric models for the accurate quantification of several analytes of interest (NO 3 – , NO 2 – , and PO 4 3– ) in highly turbid solutions with solids loadings of up to 20 wt %. This work focuses on offline sample measurement and characterization as an initial step toward the development of online monitoring capabilities. Chemical systems of interest were focused on nuclear waste at the Hanford Site, which represents highly complex matrices that could realize significant processing benefits through the integration of online monitoring.

Felmy, Heather M. [Pacific Northwest National Labo↗
DOE OSTI2021

Microscopic and Macroscopic Characterization of Grain Boundary Energy and Strength in Silicon Carbide via Machine-Learning Techniques

Predicting the properties of grain boundaries poses a challenge because of the complex relationships between structural and chemical attributes both at the atomic and continuum scales. Grain boundary systems are typically characterized by parameters used to classify local atomic arrangements in order to extract features such as grain boundary energy or grain boundary strength. The present work utilizes a combination of high-throughput atomistic simulations, macroscopic and microscopic descriptors, and machine-learning techniques to characterize the energy and strength of silicon carbide grain boundaries. Additionally, a diverse data set of symmetric tilt and twist grain boundaries are described using macroscopic metrics such as misorientation, the alignment of critical low-index planes, and the Schmid factor, but also in terms of microscopic metrics, by quantifying the local atomic structure and chemistry at the interface. These descriptors are used to create random-forest regression models, allowing for their relative importance to the grain boundary energy and decohesion stress to be better understood. Results show that while the energetics of the grain boundary were best described using the microscopic descriptors, the ability of the macroscopic descriptors to reasonably predict grain boundaries with low energy suggests a link between the crystallographic orientation and the resultant atomic structure that forms at the grain boundary within this regime. For grain boundary strength, neither microscopic nor macroscopic descriptors were able to fully capture the response individually. However, when both descriptor sets were utilized, the decohesion stress of the grain boundary could be accurately predicted. These results highlight the importance of considering both macroscopic and microscopic factors when constructing constitutive models for grain boundary systems, which has significant implications for both understanding the fundamental mechanisms at work and the ability to bridge length scales.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Remembering the Work of Phillip L. Geissler: A Coda to His Scientific Trajectory

Phillip L. Geissler made important contributions to the statistical mechanics of biological polymers, heterogeneous materials, and chemical dynamics in aqueous environments. He devised analytical and computational methods that revealed the underlying organization of complex systems at the frontiers of biology, chemistry, and materials science. In this retrospective we celebrate his work at these frontiers.

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

Architector for high-throughput cross-periodic table 3D complex building

Abstract Rare-earth and actinide complexes are critical for a wealth of clean-energy applications. Three-dimensional (3D) structural generation and prediction for these organometallic systems remains a challenge, limiting opportunities for computational chemical discovery. Here, we introduce Architector , a high-throughput in-silico synthesis code for s-, p-, d-, and f-block mononuclear organometallic complexes capable of capturing nearly the full diversity of the known experimental chemical space. Beyond known chemical space, Architector performs in-silico design of new complexes including any chemically accessible metal-ligand combinations. Architector leverages metal-center symmetry, interatomic force fields, and tight binding methods to build many possible 3D conformers from minimal 2D inputs including metal oxidation and spin state. Over a set of more than 6,000 x-ray diffraction (XRD)-determined complexes spanning the periodic table, we demonstrate quantitative agreement between Architector-predicted and experimentally observed structures. Further, we demonstrate out-of-the box conformer generation and energetic rankings of non-minimum energy conformers produced from Architector , which are critical for exploring potential energy surfaces and training force fields. Overall, Architector represents a transformative step towards cross-periodic table computational design of metal complex chemistry.

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

Spatially co-registered wide-field nonlinear optical imaging of living and complex biosystems in a total internal reflection geometry

Nonlinear optical microscopy that leverages an objective based total internal reflection (TIR) excitation scheme is an attractive means for rapid, wide-field imaging with enhanced surface sensitivity. Through select combinations of distinct modalities, one can, in principle, access complementary chemical and structural information for various chemical species near interfaces. Here, we report a successful implementation of such a wide-field nonlinear optical microscope system, which combines coherent anti-Stokes Raman scattering (CARS), two-photon fluorescence (TPF), second harmonic generation (SHG), and sum frequency generation (SFG) modalities on the same platform. The intense optical fields needed to drive these high order nonlinear optical processes are achieved through the use of femtosecond pulsed light in combination with the intrinsic field confinement induced by TIR over a large field of view. The performance of our multimodal microscope was first assessed through the experimental determination of its chemical fidelity, intensity and polarization dependences, and spatial resolution using a set of well-defined model systems. Subsequently, its unique capabilities were validated through imaging complex biological systems, including Hydrangea quercifolia pollen grains and Pantoea sp. YR343 bacterial cells. Specifically, the spatial distribution of different molecular groups in the former was visualized via vibrational contrast mechanisms of CARS, whereas co-registered TPF imaging allowed the identification of spatially localized intrinsic fluorophores. We further demonstrate the feasibility of our microscope for wide-field CARS imaging on live cells through independent characterization of cell viability using spatially co-registered TPF imaging. Overall, this approach to TIR enabled wide-field imaging is expected to provide new insights into bacterial strains and their interactions with other species in the rhizosphere in a time-resolved and chemically selective manner.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Gamma Radiolysis of Biphasic Acetohydroxamic Acid (AHA) Solvent Systems

Acetohydroxamic acid (AHA) has been proposed as a substitute for uranium(IV) and hydrazine as a plutonium complexant and neptunium reductant for simplified, single-cycle used nuclear fuel reprocessing flowsheets. However, the chemical behavior of AHA in an intense multi-component radiation field is poorly understood, especially under representative biphasic reprocessing solvent system conditions. In response to this critical knowledge gap, this study has investigated the gamma radiolytic integrity of AHA in aqueous nitric solutions in contact with an organic phase, comprising current and future reprocessing ligands (TBP, DEHBA, and DEHiBA) dissolved in n-dodecane diluent. Our data show negligible effect of the organic phase on the radiolytic behavior of AHA compared to complementary single-phase experiments.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Standoff Detection of Chemical Plumes Using Swept-ECQCL Remote Detection Technology: Results from EMRTC Complex Terrain Dual Tracer Experiment

As part of the Complex Terrain Dual Tracer Experiments performed in support of the LYNM program, a swept-wavelength external cavity quantum cascade laser (swept-ECQCL) remote detection system was deployed to measure chemical plume propagation with high temporal resolution. The swept-ECQCL remote detection technology was initially developed at PNNL, and is now under further development at Opticslah, LLC under SBIR funding from DOE (DOE SBIR Phase I Project DE-SC0019855, Remote Detection Technologies, Program Manager Chris Ramos). The swept-ECQCL technology provides a single-frequency infrared laser source which is rapidly scanned over a large wavelength range to measure the infrared absorption of gases, liquids, or solids. Analysis of the measured infrared absorption spectrum is used to detect, identify, and quantify the chemical species which are present. While most previous laser-based techniques used small wavelength tuning ranges to measure only a few isolated spectral lines of a few species, the swept-ECQCL provides the large tuning range needed to measure multiple species simultaneously, and to detect broad absorption features from large molecules or condensed phase materials. In addition, the swept-ECQCL systems have high enough spectral resolution and sensitivity to measure narrow lines of small molecules and determine isotope ratios. The remote detection system for gas plumes is based on the swept-ECQCL technology, and consists of a swept-ECQCL source, infrared detector, control electronics, and beam transmitting and alignment optics. The remote detection system directs the infrared output of the swept-ECQCL to a target retro reflector typically located at ~10-10,000 m distance. The swept-ECQCL wavelength is varied continuously over a large wavelength range (>1000 nm) at rates up to 1000 Hz. The swept ECQCL light reflected from the remote target is collected and focused onto an infrared photodetector. By measuring the detector signal continuously as the swept-ECQCL wavelength is varied, a time-series of absorption spectra are measured for the gases along the beam line of-sight. Based on analysis and fitting of the measured absorption spectra, the time-dependent chemical concentrations of species are determined. Multiple species and mixtures can be detected simultaneously, for all species with absorption features above the sensor noise floor. The sensor output is thus a continuous, high-speed record of identified chemicals along with their concentrations and isotope ratios of interest. The high-speed operation enables low-noise operation in conditions of atmospheric turbulence and allows tracking of fluctuating concentrations in transient chemical plumes as needed for reliable remote detection applications. All components operate with minimal cooling requirements and no cryogens are needed. In addition, the output infrared light is invisible with an intensity below the maximum permissible exposure (MPE) threshold of 100 mW/cm2 for these infrared wavelengths and thus eye safe. The measurement concept for plume detection using the swept-ECQCL remote detection system is as follows. The ECQCL beam is directed toward a remote retro-reflector, which reflects the beam back to the ECQCL system where it is focused onto an infrared photodetector. The line-of-sight from the ECQCL system to the retro-reflector defines the measurement path. The detected infrared signal is recorded continuously and detects the absorption of molecular species along the measurement path. Before a plume release, the system measures absorption from atmospheric constituents, which in this LWIR band consist primarily of lines from H 2 O and CO 2 , with additional weaker absorption from N 2 O, NH 3 , and O 3 possibly contributing. The system then measures any changes in absorption due to the tracer compounds when they enter the measurement path. The shape of absorption features is used to distinguish and identify the different tracer compounds, and the strength of the absorption is used to determine the path-integrated concentration – also called the column density.

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

Stoichiometrically-informed symbolic regression for extracting chemical reaction mechanisms from data

A data-driven computational method is introduced to extract chemical reaction mechanisms from time series chemical concentration data. It is realized through the use of dynamic symbolic regression in which a sparse analytical form for a dynamical system is discoverable from the underlying data. We specifically develop the stoichiometrically-informed symbolic regression (SISR) method to address a standing challenge in complex chemical reaction networks: given a time-series dataset of concentrations of several components, what is the mechanism and the associated rate constants? SISR finds the optimal mechanism, kinetic equations and rate constants by combining differential optimization with a genetic optimization approach that searches a symbolic space of possible reaction mechanisms. Use of SISR in several paradigmatic examples spanning linear and nonlinear reaction schemes results in excellent agreement between true and predicted mechanisms, including when the method is applied to noisy data. The advantages of a stoichiometrically-informed approach such as SISR to address reaction discovery is illustrated through comparison with the use of generic state-of-the-art data-driven approaches.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Critical analysis of velocimetry methods for particulate flows from synthetic data

Particle tracking methods that extract high-fidelity particle velocity data from high speed video of particle laden flows is a common experimental technique applied to chemical processes. These measurements are used to better understand the motion of particles and fluids in complex systems and create data against which computational models are validated. However, the methods, codes, and experimental setups all have limitations. It is imperative that practitioners verify the methods and their implementation as well as understand the limitations of experimental setups. This work focuses on quantifying the visible depth of field in a high particle concentration fluidized bed. Following a precedent set by the particle imaging velocimetry community, a particle velocity field is manufactured using a computational fluid dynamics and discrete element method simulation. Photo realistic high-speed videos are rendered based on the simulated data using the three-dimensional creation software Blender. Particle velocities are extracted from the synthetic high-speed videos using three variants of Particle Tracking Velocimetry and Optical Flow Velocimetry methodologies. Here, the tracked results are then compared to the known solution, quantifying the error associated with the assumed visible depth. The results indicate that at depth of one particle diameter, all three particle tracking codes give accurate measurements, largely within 5%. However, the error increases when the full bed video measurements are compared to the known solution at one particle diameter, i.e., mimicking a validation study. Finally, for some statistics the constant depth assumption only increases the error slightly, for others significantly.

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

Coupled cluster benchmarks of large noncovalent complexes: The L7 dataset as well as DNA–ellipticine and buckycatcher–fullerene

In this work, benchmark binding energies for dispersion-bound complexes in the L7 dataset, the DNA–ellipticine intercalation complex, and the buckycatcher–C60 complex with 120 heavy atoms using a focal-point method based on the canonical form of second-order Møller−Plesset theory (MP2) and the domain based local pair natural orbital scheme for the coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] extrapolated to the complete basis set (CBS) limit are reported. This work allows for increased confidence given the agreement with respect to values recently obtained using the local natural orbital CCSD(T) for L7 and the canonical CCSD(T)/CBS result for the coronene dimer (C2C2PD). Therefore, these results can be considered pushing the CCSD(T)/CBS binding benchmark to the hundred-atom scale. The disagreements between the two state-of-the-art methods, CCSD(T) and fixed-node diffusion Monte Carlo, are substantial with at least 2.0 (∼10%), 1.9 (∼5%), and 10.3 kcal/mol (∼25%) differences for C2C2PD in L7, DNA–ellipticine, and buckycatcher–C60, respectively. Such sizable discrepancy above “chemical accuracy” for large noncovalent complexes indicates how challenging it is to obtain benchmark binding interactions for systems beyond small molecules, although the three up-to-date density functionals, PBE0+D4, ωB97M-V, and B97M-V, agree better with CCSD(T) for these large systems. In addition to reporting these values, different basis sets and various CBS extrapolation parameters for Hartree–Fock and MP2 correlation energies were tested for the first time in large noncovalent complexes with the goal of providing some indications toward optimal cost effective routes to approach the CBS limit without substantial loss in quality.

Ballesteros, Francisco (ORCID:0000000325732307)↗