Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Level Set Method”

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 OSTI2021

Automatic 3D image registration for nano-resolution chemical mapping using synchrotron spectro-tomography

Nano-resolution synchrotron X-ray spectro-tomography has been demonstrated as a powerful tool for probing the three-dimensional (3D) structural and chemical heterogeneity of a sample. By reconstructing a number of tomographic data sets recorded at different X-ray energy levels, the energy-dependent intensity variation in every given voxel fingerprints the corresponding local chemistry. The resolution and accuracy of this method, however, could be jeopardized by non-ideal experimental conditions, e.g. instability in the hardware system and/or in the sample itself. Herein is presented one such case, in which unanticipated sample deformation severely degrades the data quality. Here, to address this issue, an automatic 3D image registration method is implemented to evaluate and correct this effect. The method allows the redox heterogeneity in partially delithiated Li x Ta 0.3 Mn 0.4 O 2 battery cathode particles to be revealed with significantly improved fidelity.

3D image registration↗
DOE OSTI2024

GDB-9-Ex_EOM-CCSD: Dataset containing Equation of Motion Coupled Cluster (EOM-CCSD) calculations for organic molecules of the GDB-9-Ex dataset.

This dataset contains data-intensive quantum chemical electronic structure calculations for 80,593 organic molecules of the GDB-9-Ex dataset. Calculations were performed using the Equation of Motion Coupled Cluster (EOM-CCSD) first principles method using the ORCA software. It provides UV-vis spectra calculations of molecules with a high level of accuracy. The optical spectra behavior was collected based on the optimized molecular geometries in the DFTB method with 3ob parameters. All calculations utilized the def2-TZVP basis sets with the auxiliary def2/J and def2-TZVP/C basis sets. The similarity-transformed EOM-CCSD method that used domain-based local pair natural orbitals (DLPNO) approximation which constitutes the STEOM-DLPNO-CCSD method was used. This method is based on the STEOM approach and was found to make accurate predictions of transition energies for organic molecules. For the excitation energy calculations, the lowest 50 excitation states were calculated.

AI dataset↗
DOE OSTI2024

GDB-9-Ex_TD-DFT-PBE0: Dataset containing Time Dependent Density Functional Theory (TDDFT) calculations for organic molecules of the GDB-9-Ex dataset.

This dataset contains data-intensive quantum chemical electronic structure calculations for 96,766 organic molecules of the GDB-9-Ex dataset. Calculations were performed using the Time Dependent Density Functional Theory (TDDFT) first principles method using the ORCA software. It provides UV-vis spectra calculations of molecules with a high level of accuracy. The optical spectra behavior was collected based on the optimized molecular geometries in the DFTB method with 3ob parameters. All calculations utilized the def2-TZVP basis sets with the auxiliary def2/J and def2-TZVP/C basis sets. The time-dependent density-functional theory (TDDFT) approach with the PBE0 exchange-correlation functional and ORCAs default integration grid was employed. For the excitation energy calculations, the lowest 50 excitation states were calculated.

AI dataset↗
DOE OSTI2022

TDCOSMO VIII. A key test of systematics in the hierarchical method of time-delay cosmography

The largest source of systematic errors in the time-delay cosmography method likely arises from the lens model mass distribution, where an inaccurate choice of model could in principle bias the value of $H$ 0 . A Bayesian hierarchical framework has been proposed which combines lens systems with kinematic data, constraining the mass profile shape at a population level. The framework has been previously validated using a small sample of lensing galaxies drawn from hydro-simulations. The goal of this work is to expand the validation to a more general set of lenses consistent with observed systems, as well as confirm the capacity of the method to combine two lens populations: one which has time delay information and one which lacks time delays and has systematically different image radii. For this purpose, we generated samples of analytic lens mass distributions made of baryons+dark matter and fit the subsequent mock images with standard power-law models. Corresponding kinematics data were also emulated. The hierarchical framework applied to an ensemble of time-delay lenses allowed us to correct the $H$ 0 bias associated with model choice to find $H$ 0 within 1.5σ of the fiducial value. We then combined this set with a sample of corresponding lens systems which have no time delays and have a source at lower $z$, resulting in a systematically smaller image radius relative to their effective radius. The hierarchical framework has successfully accounted for this effect, recovering a value of $H$ 0 which is both more precise (σ ~ 2%) and more accurate (0.7% median offset) than the time-delay set alone. This result confirms that non-time-delay lenses can nonetheless contribute valuable constraining power to the determination of $H$ 0 via their kinematic constraints, assuming they come from the same global population as the time-delay set.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Performance Results for Sensor Assignment Problem as Solved on a Multi-Node Cluster

An earlier report described a procedure for optimal sensor set selection and its implementation on a computational cluster. This new and innovative capability was developed to facilitate a reduction in operations staffing levels to improve plant economics. By automating surveillance and maintenance tasks through early detection of degrading sensors and equipment, staff can be more efficiently deployed. The method uses automated reasoning and domain knowledge in the form of the conservation equations to infer from plant measurements the state of equipment health. Inclusion of domain knowledge addresses the problem that exists with pure data-driven methods that there are no rigorous guidelines for determining what constitutes an adequate sensor set. Formalizing the procedure for sensor set selection as we have done results in a more reliable and explainable diagnosis of plant equipment health. Importantly, from the standpoint of the plant owner, personnel are provided with an early and explicit diagnosis of an equipment problem. That in principle automates the process and eliminates having to send personnel into the plant to find the cause as typically occurs when a data-driven method detects an anomaly. In this report we describe first results obtained using a computational cluster to solve the sensor set selection problem as framed above. The case described addresses the problem of equipment health monitoring in the high-pressure (HP) feedwater system of a pressurized light water reactor as seen through the eyes of our collaborating utility partner. Maintenance of this system can amount to millions of dollars per year if equipment health issues go undiagnosed and lead to loss of function. On examining the potential that is inherent in the installed sensor set for diagnosing equipment health degradation, it was found that greater fault resolution capability can be achieved using a sensor set that is 20 percent fewer in number. The take-away is that compared to the installed sensor set there exists a more strategic assignment of sensors that will furnish better health monitoring capability and with fewer sensors. Where the problem defies solution by manual inspection, as is the case here, one can be found by an algorithm. The solution was obtained in four hours using 30 computational cores. The HP feedwater problem as posed above illustrates the added value of approaching the sensor selection problem as one amenable to algorithmic solution. This problem is of interest to advanced reactor designers and to utilities that are setting up remote monitoring and diagnostic centers.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Multi-fidelity learning for interatomic potentials: low-level forces and high-level energies are all you need

The promise of machine learning interatomic potentials (MLIPs) has led to an abundance of public quantum mechanical (QM) training datasets. The quality of an MLIP is directly limited by the accuracy of the energies and atomic forces in the training dataset. Unfortunately, most of these datasets are computed with relatively low-accuracy QM methods, e.g. density functional theory with a moderate basis set. Due to the increased computational cost of more accurate QM methods, e.g. coupled-cluster theory with a complete basis set (CBS) extrapolation, most high-accuracy datasets are much smaller and often do not contain atomic forces. The lack of high-accuracy atomic forces is quite troubling, as training with force data greatly improves the stability and quality of the MLIP compared to training to energy alone. Because most datasets are computed with a unique level of theory, traditional single-fidelity (SF) learning is not capable of leveraging the vast amounts of published QM data. In this study, we apply multi-fidelity learning (MFL) to train an MLIP to multiple QM datasets of different levels of accuracy, i.e. levels of fidelity. Specifically, we perform three test cases to demonstrate that MFL with both low-level forces and high-level energies yields an extremely accurate MLIP—far more accurate than a SF MLIP trained solely to high-level energies and almost as accurate as a SF MLIP trained directly to high-level energies and forces. Therefore, MFL greatly alleviates the need for generating large and expensive datasets containing high-accuracy atomic forces and allows for more effective training to existing high-accuracy energy-only datasets. Indeed, low-accuracy atomic forces and high-accuracy energies are all that are needed to achieve a high-accuracy MLIP with MFL.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Bond Dissociation Energies of the Actinide Halides AnX, An = Ac–Lr and X = F–I, Utilizing Relativistic Composite Coupled Cluster Approaches

Bond dissociation energies (BDEs) have been calculated for the set of actinide halides AnX with An=Ac, Pa, and Np-Lr and X=F-I. Two composite thermochemistry methods based on the Feller-Peterson-Dixon (FPD) approach have been utilized, one involving spinor-based relativistic CCSD(T) calculations where spin-orbit (SO) was included at the orbital level and another using scalar relativistic CCSD(T) with a posteriori SO contributions based on 2-component multireference configuration interaction calculations. The method that was chosen for a given actinide halide was based on which representation yielded the best single determinant reference determinant for the coupled cluster calculation. The spinor-based method was chosen for all cases except for AmX, CmX, and BkX. Both composite approaches included contributions accounting for basis set truncation, outer-core correlation, the Gaunt interaction, and QED. The scalar FPD results, as well as the spinor-based calculations for AcF, also included higher order electron correlation up through CCSDT(Q). In addition to BDEs, CCSD(T) equilibrium bond lengths, harmonic frequencies, and vibrational anharmonicity constants are reported for all species. Last, the FPD BDEs for the fluorides were used to confirm the trend across the actinide series previously predicted by Gibson using bonding models based atomic promotion energies that provide a single 6d electron for bonding. In particular the local minimum in the BDEs at AmF is confirmed in the present calculations. Furthermore, the BDEs for LrX are predicted to be slightly larger than those of AcX, making them the largest in the actinide halide series.

Actinides↗
DOE OSTI2024

Random coordinate descent: A simple alternative for optimizing parameterized quantum circuits

Variational quantum algorithms rely on the optimization of parameterized quantum circuits in noisy settings. The commonly used back-propagation procedure in classical machine learning is not directly applicable in this setting due to the collapse of quantum states after measurements. Thus, gradient estimations constitute a significant overhead in a gradient-based optimization of such quantum circuits. This paper introduces a random coordinate descent algorithm as a practical and easy-to-implement alternative to the full gradient descent algorithm. This algorithm only requires one partial derivative at each iteration. Motivated by the behavior of measurement noise in the practical optimization of parameterized quantum circuits, this paper presents an optimization problem setting that is amenable to analysis. Under this setting, the random coordinate descent algorithm exhibits the same level of stochastic stability as the full gradient approach, making it as resilient to noise. The complexity of the random coordinate descent method is generally no worse than that of the gradient descent and can be much better for various quantum optimization problems with anisotropic Lipschitz constants. Theoretical analysis and extensive numerical experiments validate our findings. Published by the American Physical Society 2024

Ding, Zhiyan (ORCID:000000018863403X)↗
DOE OSTI2022

An Assessment of Hydropower Potential at National Conduits

In support of the mission of the US Department of Energy Water Power Technologies Office to promote domestic next-generation hydropower growth, this study provides a first-of-its-kind assessment of the hydropower development potential on existing water conduits nationwide. This study evaluates the potential for conduit hydropower development in the municipal, agricultural, and industrial sectors across the United States. The research team developed systematic methods to evaluate conduit hydropower capacity potential (MW) and energy generation potential (GWh/year) at four categories of conduits across the United States. They include the following: (1) Water supply pipelines for municipal and industrial uses, (2) Wastewater discharge conduits from municipal and industrial systems, (3) Agricultural water conduits including irrigation canals and ditches in the 17 western states that rely heavily on irrigation, and (4) Thermoelectric power plant cooling water discharge conduits. For each type of conduit, the team developed and implemented a method to estimate hydraulic head and annual water flows—and the hydropower potential—based on analyses of satellite imagery, topography, and existing data sets on water systems and power plants. These methods support a consistent, replicable estimate of conduit hydropower potential across 3 sectors and all 50 states, aggregated at the county and state levels. The assessment was conducted at the reconnaissance level, considering resources that could be available for development at the state and national levels using present-day assumptions about conduit hydropower technology.

13 HYDRO ENERGY↗
DOE OSTI2020

Assessing United States County-Level Exposure for Research on Tropical Cyclones and Human Health

Tropical cyclone epidemiology can be advanced through exposure assessment methods that are comprehensive and consistent across space and time, as these facilitate multiyear, multistorm studies. Further, an understanding of patterns in and between exposure metrics that are based on specific hazards of the storm can help in designing tropical cyclone epidemiological research. a) Provide an open-source data set for tropical cyclone exposure assessment for epidemiological research; and b) investigate patterns and agreement between county-level assessments of tropical cyclone exposure based on different storm hazards. We created an open-source data set with data at the county level on exposure to four tropical cyclone hazards: peak sustained wind, rainfall, flooding, and tornadoes. The data cover all eastern U.S. counties for all land-falling or near-land Atlantic basin storms, covering 1996–2011 for all metrics and up to 1988–2018 for specific metrics. We validated measurements against other data sources and investigated patterns and agreement among binary exposure classifications based on these metrics, as well as compared them to use of distance from the storm’s track, which has been used as a proxy for exposure in some epidemiological studies. Our open-source data set was typically consistent with data from other sources, and we present and discuss areas of disagreement and other caveats. Over the study period and area, tropical cyclones typically brought different hazards to different counties. Therefore, when comparing exposure assessment between different hazard-specific metrics, agreement was usually low, as it also was when comparing exposure assessment based on a distance-based proxy measurement and any of the hazard-specific metrics. Our results provide a multihazard data set that can be leveraged for epidemiological research on tropical cyclones, as well as insights that can inform the design and analysis for tropical cyclone epidemiological research.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

VARI3D & PERSENT: Perturbation and Sensitivity Analysis

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code. This manuscript serves as a single manual for two separate codes: VARI3D and PERSENT. The VARI3D code (VARIational 3D) is based upon the classic finite difference diffusion theory solver available in DIF3D. The PERSENT code (PERturbation and SENitivity for Transport) is based upon the variational nodal method employed in DIF3D termed VARIANT. The VARIANT solver was added to DIF3D in 1995 and has seen continued development and use for the last 18 years. Because VARI3D primarily uses deprecated coding practices, rather than incorporating the perturbation and sensitivity treatments for transport within VARI3D, a new coding development was built using modern Fortran coding. The primary purpose of this manual is to describe the theory behind PERSENT (and by convenience, that of VARI3D) and discuss the input and output of PERSENT along with giving potential users an idea of how to use it. While this manuscript does describe the input and output of VARI3D, the PERSENT code is intended to be the replacement capability of VARI3D as PERSENT can generate nearly identical (if not superior) diffusion theory results. In this manuscript, the relevant aspects of generalized perturbation theory and exact perturbation theory that apply to both VARI3D and PERSENT are covered. The input and output of VARI3D is displayed by excerpting several of the example problems. Similarly, the input and output of PERSENT is displayed along with tips on how best to use the code. Note that the input and output of the inhomogeneous solver wrapped around DIF3D (DIF3D_IFS) is also discussed as it is needed to carry out some of the sensitivities in PERSENT such as reaction rate ratios. This manuscript describes several perturbation and sensitivity problems, and the results computed using PERSENT. From these sections, potential users should find that PERSENT provides not only the typical tables of numbers desired in perturbation and sensitivity analysis work, but also can visually plot the result for a more thorough understanding of the space and energy distribution (Section 5). Overall, PERSENT is observed to produce accurate reactivity worths and sensitivities for the displayed set of test problems and clearly demonstrates the need to have a transport-based sensitivity capability as evident from the thousands of percent errors observed in the 21-group hexagonal fast reactor problem (covered in Section 7). The uncertainty calculation capability is described in Section 3 and demonstrated in Section 7.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Complete collision data set for electrons scattering on molecular hydrogen and its isotopologues: I. Fully vibrationally-resolved electronic excitation of H 2 $X^1Σ^+_g)$

Here, we present a comprehensive set of vibrationally-resolved cross sections for electron-impact electronic excitation of molecular hydrogen suitable for implementation in collisional-radiative models. The adiabatic-nuclei molecular convergent close-coupling method is used to calculate cross sections for excitation of all bound vibrational levels and dissociative excitation of the B 1 Σ u + , C 1 Π u , E F 1 Σ g + , B ′ 1 Σ u + , G K 1 Σ g + , I 1 Π g , J 1 Δ g , D 1 Π u , H 1 Σ g + , b 3 Σ u + , c 3 Π u , a 3 Σ g + , e 3 Σ u + , d 3 Π u , h 3 Σ g + , g 3 Σ g + , i 3 Π g , and j 3 Δ g electronic states from all –14 bound vibrational levels of the ground electronic ($\text{X}$ $^1&#x3A3^+_g$) state. The data set consists of cross sections from threshold to 500 eV for over 5000 transitions, representing all possible electronic and vibrational transitions between the state and the –3 singlet and triplet states (where refers to the united-atoms-limit principle quantum number). The cross sections are presented in graphical form and provided as both numerical values and analytic fit functions in supplementary data files.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Custom Accessors: Enabling Scalable Data Ingestion, (Re-)Organization, and Analysis on Distributed Systems

The emerging class of high velocity and high volume data analytic workflows comprise interwoven data ingestion, organization, and processing stages, with ingestion and organization steps often contributing comparable or even higher computational costs than actual processing steps. Since complex workflows consist of a variety of phases that view and use data differently, being able to construct efficient, scalable, distributed data structures (arrays, vectors, sets, maps, and multi-maps) is essential and requires custom methods to extend and shrink containers, analyze and position data, and, maintain globallyconsistent meta-data. In this paper, we propose a novel datastructure access paradigm based on the concept of Accessors. At a high level, accessors are customizable callable objects that can modify the behavior of insert, read, update, and delete operations for distributed containers while preserving atomicity guarantees. Accessors provide a very clean and natural way to implement a variety of programming patterns, e.g., conditional insertion/deletion and cascading computations, which would be otherwise hard (or even impossible) to express in parallel and distributed settings without using locks. We demonstrate the practicality and usefulness of our approach with two representative use cases and study the performance of these applications on a distributed High-Performance Computing system. Our analysis highlights that our proposed abstraction allows for an effective overlapping and concurrent execution of different workflow steps (e.g., data ingestion and analysis), which in a conventional analytics pipeline would execute sequentially, contributing cumulatively to the overall latency.

Castellana, Vito G. [BATTELLE (PACIFIC NW LAB)] (O↗
DOE OSTI2025

Engineering Polyketide Stereocenters with Ketoreductase Domain Exchanges

Polyketide synthases (PKSs) are versatile biosynthetic megasynthases capable of producing a diverse range of natural products with many applications, including in pharmaceuticals. The stereochemical precision of PKSs makes them a powerful tool for engineering tailored, unnatural polyketides; however, modifying the stereocenters of a PKS product while maintaining production levels remains a significant challenge. In this study, we systematically tested and evaluated strategies for ketoreductase (KR) domain exchanges, the domain responsible for setting stereocenters of polyketide products. After first optimizing the method for KR exchanges, we then performed 44 KR domain exchanges on three different PKSs to obtain high production of all four stereoisomers in vivo. By testing both one- and two-module PKS systems, we investigated how downstream modules process intermediates with altered stereochemistry and found that the configuration of the α-substituents was critical for gatekeeping by the ketosynthase (KS). To overcome this constraint, we investigated two different strategies for altering the KS domain, including introducing targeted mutations in the downstream KS, and exploring boundaries in exchanging the entire functional unit from the donor PKS. Both strategies successfully modified the KS stereocontrol with distinct trade-offs; the functional unit exchange resulted in higher titer improvements, though it was more likely to break the entire PKS. This study demonstrates a comprehensive approach to successfully engineering all four stereochemical configurations in multiple PKS systems, advancing our understanding of and ability to rationally modify polyketide stereochemistry through multiple engineering strategies.

Keiser, Leah S. [Joint BioEnergy Institute (JBEI),↗
DOE OSTI2023

Fringing Analysis and Simulation for the Vera C. Rubin Observatory’s Legacy Survey of Space and Time

Abstract The presence of fringing in astronomical Charge-Coupled Device (CCD) images will have an impact on photometric quality and measurements. Yet its impact on the Vera C. Rubin Observatory’s Legacy Survey of Space and Time (LSST) has not been fully studied. We present a detailed study on fringing for CCDs already implemented on the Rubin Observatory LSST Camera’s focal plane. After making physical measurements and knowing the compositions, we developed a model for the e2v CCDs. We present a method to fit for the internal height variation of the epoxy layer within the sensors based on fringing measurements in a laboratory setting. This method is generic enough that it can be easily modified to work for other CCDs. Using the derived fringing model, we successfully reproduce comparable fringing amplitudes that match the observed levels in images taken by existing telescopes with different optical designs. This model is then used to forecast the expected level of fringing in a single LSST y -band sky background exposure with Rubin telescope optics in the presence of a realistic time-varying sky spectrum. The predicted fringing amplitude in LSST images ranges from 0.04% to 0.2% depending on the location of a CCD on the focal plane. We find that the predicted variation in surface brightness caused by fringing in LSST y -band sky background images is about 0.6 μ Jy arcsec − 2 , which is 40 times larger than the current measurement error. We conclude that it is necessary to include fringing correction in the Rubin’s LSST image processing pipeline.

Astronomy & Astrophysics↗
DOE OSTI2024

The AGORA High-resolution Galaxy Simulations Comparison Project. VI. Similarities and Differences in the Circumgalactic Medium

Abstract We analyze the circumgalactic medium (CGM) for eight commonly-used cosmological codes in the AGORA collaboration. The codes are calibrated to use identical initial conditions, cosmology, heating and cooling, and star formation thresholds, but each evolves with its own unique code architecture and stellar feedback implementation. Here, we analyze the results of these simulations in terms of the structure, composition, and phase dynamics of the CGM. We show properties such as metal distribution, ionization levels, and kinematics are effective tracers of the effects of the different code feedback and implementation methods, and as such they can be highly divergent between simulations. This is merely a fiducial set of models, against which we will in the future compare multiple feedback recipes for each code. Nevertheless, we find that the large parameter space these simulations establish can help disentangle the different variables that affect observable quantities in the CGM, e.g., showing that abundances for ions with higher ionization energy are more strongly determined by the simulation’s metallicity, while abundances for ions with lower ionization energy are more strongly determined by the gas density and temperature.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Lagrangian–Eulerian multidensity topology optimization with the material point method

Abstract In this paper, a hybrid Lagrangian–Eulerian topology optimization (LETO) method is proposed to solve the elastic force equilibrium with the Material Point Method (MPM). LETO transfers density information from freely movable Lagrangian carrier particles to a fixed set of Eulerian quadrature points. This transfer is based on a smooth radial kernel involved in the compliance objective to avoid the artificial checkerboard pattern. The quadrature points act as MPM particles embedded in a lower‐resolution grid and enable a subcell multidensity resolution of intricate structures with a reduced computational cost. A quadrature‐level connectivity graph‐based method is adopted to avoid the artificial checkerboard issues commonly existing in multiresolution topology optimization methods. Numerical experiments are provided to demonstrate the efficacy of the proposed approach.

Li, Yue↗
DOE OSTI2026

Concordant Mode Approach (CMA): Vibrational Analysis of New and Upgraded Intermolecular Benchmarks for Noncovalent Bonding

The Concordant Mode Approach (CMA) is a novel method that offers tremendous potential for increasing the system size and the level of theory attainable in quantum chemical computations of molecular vibrational frequencies. To investigate the extension of CMA to intermolecular vibrations, computations with coupled cluster singles and doubles with perturbative triples theory [CCSD(T)] using two augmented correlation-consistent polarized-valence triple-ζ basis sets (aug-cc-pVTZ or h-aug-ccpVTZ) were performed on 17 prototypical loosely bound complexes of hydrogen-bonded, dispersion, and mixed character. These Level A results provide new and upgraded benchmarks for noncovalent bonding and a severe test for CMA vibrational analyses. The Level A target frequencies were recovered remarkably well using second-order Møller−Plesset perturbation theory (MP2) with h-aug-cc-pVTZ for generating the underlying (Level B) normal modes of the CMA scheme. Employing this Level B within the lowest-rung CMA-0A method reproduces the 435 benchmark frequencies with a mean absolute error (MAE) of 0.23 cm −1 and a corresponding standard deviation (σ) of 0.84 cm −1 ; strikingly, the corresponding subset of 106 interfragment frequencies exhibits MAE = 0.34 cm −1 and σ = 0.90 cm −1 . Subsequent application of the higher-rung CMA-2A scheme eliminates all outliers and reduces the overall MAE to a minuscule 0.08 cm −1 with the inclusion of only 3.0% of the off-diagonal couplings not accounted for by CMA-0A. Accordingly, the highly efficient CMA methodology proves to be robust even for vibrations on flat potential energy surfaces.

Aromatic compounds↗