Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Surface Energy”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Ab Initio Bulk Free Energy Surface of Proper Ferroelectrics

We report a systematic and accurate approach for deriving the bulk free energy surface (FES), a function of temperature, polarization, and strain, from the first-principles density functional theory (DFT) of proper ferroelectrics. The core of our approach is the metadynamics algorithm that extracts the polarization dependence of the FES from all-atom molecular dynamics simulations without an a priori ansatz. The rest of the FES is derived from the metadynamics trajectories that span the relevant phase space. We demonstrate our approach in the case of lead titanate. The errors across the phase transition, due to DFT numerics, all-atom molecular dynamics, and free energy evaluation by enhanced sampling, can be systematically controlled and are of the order of 1 meV/atom. The accuracy of the resulting ab initio FES is only limited by the adopted functional approximation of DFT.

Xie, Pinchen [Lawrence Berkeley National Laborator↗
DOE OSTI2026

Surface Energy Modulated Self‐Assembly of Percolative Gallium Nanodroplet Films for Stretchable Electronics

Advancements in soft robotics and wearable technologies have led to an increased need for flexible and stretchable electronic systems. Liquid metals like Gallium (Ga) and its alloys offer high potential in the field of stretchable electronics and photonics due to their inherent high conductivity and malleability. However, their high surface tension and native oxide formation hinder uniform wetting and continuous film formation to get conducting films. Here, we report a fabrication strategy to realize electrically conductive Ga thin films with effective thicknesses around 100 nm on polymeric and elastomeric substrates. By enhancing the surface energy of the elastomer, we create uniform Ga nanodroplets capable of sustaining stretchability up to 200% strain. Our devices exhibit electronic conduction primarily through tunneling and display strain tunable piezoresistive characteristics, which are explained using Kirchhoff's resistive network and Simmon's analytical tunneling current model. To showcase the device's capability in fields of stretchable electronics, we demonstrate several functional prototypes, including flexible devices integrated with Light Emitting Diodes (LEDs), curvature and pressure sensors, and soft robotic grippers capable of sensing object curvature. The demonstrated devices highlight the potential of this approach for scalable integration of stretchable electronics in soft robotic and wearable platforms.

42 ENGINEERING↗
DOE OSTI2025

Interplay of surface energy and rheology in biopolymer soil enhancement

Biopolymers such as xanthan gum (XG) and locust bean gum (LBG) hold great potential as eco-friendly alternative soil binders. In this work, we investigated the impact of XG/LBG mixtures on the unconfined compressive strength (UCS) of sand. The high strength of dry biopolymer/sand arises from the cohesion between solid polymer films and sand particles which supported by work of adhesion calculation and soil mechanics measurement. LBG exhibits much lower sand reinforcement efficacy because polymers unevenly distributed within sand matrix. The formation of a core-shell structure in LBG/sand is an interplay of surface free energy and viscoelastic properties of polymer solutions. This structure is altered when LBG mixed with XG at varying ratios as those physical properties changed due to the complexity of polymer chains association. By probing these factors, we aim to elucidate the role of surface energies and polymer physics in governing the strength of the sand/polymer network, thereby contributing to a more comprehensive understanding polymer-sand interface. The low strength of gels (G’ ∼10Pa) cannot solely account for the increased UCS of wet sand over 10 kPa. Instead, the high strength of biopolymer/sand is more likely derived from the granular particles with biopolymers as solid glue.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Quantum scattering of HC 5 N and para -H 2 on a new potential energy surface

In the interstellar medium (ISM), non-local thermodynamic equilibrium situations are common due to low density, and one needs to consider the effect of molecular collisions in order to interpret the observations. Among the species detected in the ISM, cyanopolyynes, with the general molecular formula HC 2n+1 N (n = 1, 2, …), are characterized by large dipole moments and small rotational constants and constitute an indispensable class of candidates for the sensitive tracers of local density and temperature. We present a study of the collisional (de-) excitation of HC 5 N by para -H 2 (p-H 2 ) in its ground rotational state, namely HC 5 N ( j 1 ) + H 2 ( j 2 = 0) → HC 5 N (j$_1^′$) + H2 (j$_2^′$ = 0), where j 1 (or j$_1^′$) and j 2 (or j$_2^′$) denote the initial (or final) rotational quantum numbers of HC 5 N and H 2 , respectively. We performed the quantum scattering calculations at low collision energy using a new four-dimensional ab initio potential energy surface. In the regime where p-H 2 remains in its rotational ground state, converged cross sections did not require including excited rotational states of p-H 2 in the rotational basis. State-to-state cross sections were computed by means of the quantum-mechanical close-coupling (CC) method and the coupled states (CS) approximation, and rate coefficients for the first 61 levels of HC 5 N were computed for the first time up to 20 K with the CC approach and up to 50 K with the CS method. CC and CS results were found to agree well at temperatures up to 20 K. Finally, these data should allow a more accurate derivation of the HC 5 N abundance in molecular clouds.

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

Pair wave function symmetry in UTe2 from zero-energy surface state visualization

Although nodal spin-triplet topological superconductivity appears probable in uranium ditelluride (UTe 2 ), its superconductive order parameter Δ k remains unestablished. In theory, a distinctive identifier would be the existence of a superconductive topological surface band, which could facilitate zero-energy Andreev tunneling to an s-wave superconductor and also distinguish a chiral from a nonchiral Δ k through enhanced s-wave proximity. In this study, we used s-wave superconductive scan tips and detected intense zero-energy Andreev conductance at the UTe 2 (0-11) termination surface. Imaging revealed subgap quasiparticle scattering interference signatures with a-axis orientation. The observed zero-energy Andreev peak splitting with enhanced s-wave proximity signifies that Δ k of UTe 2 is a nonchiral state: B 1u , B 2u , or B 3u . However, if the quasiparticle scattering along the a axis is internodal, then a nonchiral B 3u state is the most consistent for UTe 2 .

Gu, Qiangqiang [Cornell Univ., Ithaca, NY (United ↗
DOE OSTI2025

Observations of surface energy fluxes and meteorology in the seasonally snow-covered high-elevation East River watershed during SPLASH, 2021–2023

Abstract. From autumn 2021 through summer 2023, scientists from the National Oceanic and Atmospheric Administration (NOAA) and partners conducted the Study of Precipitation, the Lower Atmosphere, and Surface for Hydrometeorology (SPLASH) campaign in the East River watershed of Colorado. One objective of SPLASH was to observe the transfer of energy between the atmosphere and the surface, which was done at several locations. Two remote sites were chosen that did not have access to power utilities. These were along the valley floor near the East River in the vicinity of the unincorporated town of Gothic, Colorado. Energy balance measurements were made at these locations using autonomous, single-level flux towers referred to as atmospheric surface flux stations (ASFSs). The ASFSs were deployed on 28 September 2021 at the Kettle Ponds Annex site and on 12 October 2021 at the Avery Picnic site and operated until 19 July and 21 June 2023, respectively. Measurements included basic meteorology; upward and downward longwave and shortwave radiative fluxes and subsurface conductive flux, each at 1 min resolution; 3-D winds from a sonic anemometer and H2O/CO2 from an open-path gas analyzer, both at 20 Hz from which sensible, latent heat, and CO2 fluxes were derived; and profiles of soil properties in the upper 0.5 m (both sites) and temperature profiles through the snow (at Avery Picnic), each reported between 10 min and 6 h. The system uptime was 97 % (Kettle Ponds) and 90 % (Avery Picnic), and collectively 1184 d of data was obtained between the stations. The purpose of this article is to document the ASFS deployment at SPLASH, to document the data acquisition and post-processing of measurements, and to serve as a guide for interested users of the data sets, which are archived at Zenodo (https://doi.org/10.5281/zenodo.10313363, Cox et al., 2023b; https://doi.org/10.5281/zenodo.10327409, Cox et al., 2023c; https://doi.org/10.5281/zenodo.10313894, Cox et al., 2023d; https://doi.org/10.5281/zenodo.10307825, Cox et al., 2023e; https://doi.org/10.5281/zenodo.10310520, Cox et al., 2023f) with the Creative Commons Attribution 4.0 International license.

Cox, Christopher J. (ORCID:0000000322037173)↗
DOE OSTI2025

Cloud radiative effect dominates variabilities of surface energy budget in the dark Arctic

Climate models simulate a wide range of temperatures in the Arctic. Here we investigate one of the main drivers of changes in surface temperature: the net surface heat flux in the models. We show that in the winter months of the dark Arctic, there is a more than two-fold difference in the net surface heat fluxes among the models, and this difference is dominated by the downward infrared radiation from clouds. Owing to the small amount of water vapor in the winter Arctic, infrared radiation from clouds transmits more easily to the surface in the Arctic than at other latitudes, resulting in large cloud radiative effect at the surface. The dominant role of the cloud effect is also found in the transient variability of the net surface heat flux. Results demonstrate that accurate simulation of clouds is crucial for determining the net surface heat flux, which in turn affects surface temperature and sea ice properties in the Arctic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Basin-Size Mapping: Prediction of Metastable Polymorph Synthesizability Across TaC–TaN Alloys

The sizes of the basins of attraction on the potential energy surface are helpful indicators in determining the experimental synthesizability of metastable phases. In principle, these basins can be controlled with changes in thermodynamic conditions such as composition, pressure, and surface energy. Herein, we use random structure sampling to computationally study how alloying smoothly perturbs basin of attraction sizes. The TaC 1-x N x pseudobinary is an ideal test system given the structural and polymorphic contrast of its parent compounds and their technological relevance as epitaxial substrates for Al 1-x Ga x N. While we find limited thermodynamic stability across all computationally observed phases, random structure sampling shows a significant composition region where the rocksalt basin dominates. As such, we predict the potential for the nonequilibrium synthesis of metastable rocksalt TaC 1-x N x alloys as substrates for Al 1-x Ga x N. At higher nitrogen concentrations, other low-energy metastable polymorphs emerge that continue to retain the hexagonal close packing suitable for III-N growth. Confidence in these trends was established through uncertainty quantification of the basin sizes and energy distributions; such analysis utilized the Beta and Dirichlet distributions. In conclusion, we also find (a) polymorph basin sizes can be rationalized in terms of energetic preferences for different coordination environments; and (b) basin sizes universally shrink with increasing nitrogen content, making the system more prone to amorphous growth.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Review of High-Throughput Surface Treatments for Microlens Arrays

Microlenses are increasingly being integrated into modern manufactured devices. From printed security devices and screens to solar panels and microscopes, these optical materials offer high control over light focusing. Thus, understanding how to treat the surfaces of these fragile, transparent devices on an integrated manufacturing line is essential. Here, in this study, we review the surface treatments for the following application categories: cleaning, increasing surface energy, decreasing surface energy, and tunable surface modifications. This overview describes methods available for the large-scale manufacturing of microlens arrays and the potential impact of those treatments on common optical surfaces. Objectives and qualitative compatibility parameters are compared, and outlooks are provided for further study to aid in streamlining the method selection and process optimization for microlenses and similar optical components.

lens manufacturing↗
DOE OSTI2024

Numerical analysis of coalescence-induced bubble departure for enhanced boiling heat transfer

Boiling heat transfer plays a crucial role in a wide range of applications, such as power generation, refrigeration, electronics cooling, and pharmaceutics. Among the various factors that influence boiling heat transfer, the dynamics of vapor bubble nucleation, growth, and departure from the heated surface stand out as particularly important. An emerging phenomenon that can promote the departure of bubbles smaller than the Fritz diameter is coalescence-induced departure. If the dynamics of this process are fully understood, then surfaces can be engineered to promote faster bubble departure and substantially increase the performance of boiling heat transfer. Further, this work expands on published results by presenting a detailed numerical analysis of bubble coalescence and departure for a range of initial bubble diameters and size ratios between coalescing bubbles. Analysis of the results is focused on explaining how the release of surface energy and bubble surface dynamics lead to bubble departure, as well as fundamentally distinguishing capillary–inertial jumping and buoyant–inertial departure mechanisms across different bubble sizes and size ratios. The results show that both the initial sizes of the coalescing bubbles and the ratio between their sizes can determine whether the merged bubble will leave the surface through capillary–inertial jumping or buoyant departure. Below a certain bubble size, the release of surface energy by the merger is not sufficient to propel the merged bubble from the surface.

42 ENGINEERING↗
DOE OSTI2026

Protein Adhesion on Semi-Fluorinated Polystyrene Surfaces in Static and Dynamic Measurements

Reducing protein adhesion is a critical strategy in fouling-resistant material innovation, with broad applications spanning biomedical and healthcare devices, biosensors, industrial and environmental systems, and other important technological domains. Here, in this study, we elucidated protein adhesion behavior on polystyrene-based thin films by neutron reflectometry (NR) and quartz crystal microbalance with dissipation (QCM-D), using both lysozyme and bovine serum albumin (BSA) as model proteins. To this end, semifluorinated polystyrene thin films with gradient wettability and surface energy were fabricated through dry processing using plasma oxidation and gas-phase deposition. Although it is believed that a fully fluorinated alkyl chain offers extremely low surface energy, thus rejecting foulants, and has been used in many fouling-resistant surface designs, enhanced protein–surface interactions were observed consistently in NR and QCM-D results, due to the combined effects of surface morphology and chemistry. On the contrary, depositing shorter fluorinated silane onto a hydrophilic PS surface contributed to a more homogeneous nanoscale fluorine coating, resulting in less initial protein adsorption and improved surface recovery. Comparative analysis of proteins with different sizes on the nanopatterned semifluorinated surface revealed the influence of molecular characteristics on surface interactions. Lysozyme, being smaller and more compact, showed faster adsorption kinetics and higher surface coverage but largely reversible binding, whereas BSA, with its larger and more flexible structure, formed broader and more stable interfacial layers. This study fills the gap in understanding protein adhesion within the range of hydrophobicity (water contact angle ∼90°), as current strategies often associate with extreme hydrophilic and superhydrophobic surfaces due to hydration or low-surface-energy rejection mechanisms, respectively. It also provides in-depth insights into current combinatorial fouling-resistant surface design.

Yuan, Yue [Oak Ridge National Laboratory (ORNL), O↗
DOE OSTI2024

ELECTRONIC STRUCTURE METHODS AND PROTOCOLS WITH APPLICATION TO DYNAMICS, KINETICS AND THERMOCHEMISTRY

Hydrocarbon combustion involves the reaction dynamics of a tremendous number of species beginning with many-component fuel mixtures and proceeding via a complex system of intermediates to form primary and secondary products. Combustion conditions corresponding to new advanced engines and/or alternative fuels rely increasingly on autoignition and low-temperature-combustion chemistry. In these regimes various transient radical species such as HO2, ROO·, ·QOOH, HCO, NO2, HOCO, and Criegee intermediates play important roles in determining the detailed as well as more general dynamics. A clear understanding and accurate representation of these processes is needed for effective modeling. Given the difficulties associated with making reliable experimental measurements of these systems, computation can play an important role in developing these energy technologies. Accurate calculations have their own challenges since even within the simplest dynamical approximations such as transition state theory, the rates depend exponentially on critical barrier heights and these may be sensitive to the level of quantum chemistry. Moreover, it is well-known that in many cases it is necessary to go beyond statistical theories and consider the dynamics. Quantum tunneling, resonances, radiative transitions, and non-adiabatic effects governed by spin-orbit or derivative coupling can be determining factors in those dynamics. Building upon progress made during a period of prior support through the DOE Early Career Program, this project combines developments in the areas of potential energy surface (PES) fitting and multistate multireference quantum chemistry to allow spectroscopically and dynamically/kinetically accurate investigations of key molecular systems (such as those mentioned above), many of which are radicals with strong multireference character and have the possibility of multiple electronic states contributing to the observed dynamics. An ongoing area of investigation is to develop general strategies for robustly convergent electronic structure theory for global multichannel reactive surfaces including diabatization of energy and other relevant surfaces such as dipole transition. Combining advances in ab initio methods with automated interpolative PES fitting allows the construction of high-quality PESs (incorporating thousands of high-level data) to be done rapidly through parallel processing on high-performance computing (HPC) clusters. In addition, new methods and approaches to electronic structure theory will be developed and tested through applications. This project will explore limitations in traditional multireference calculations (e.g., MRCI) such as those imposed by internal contraction, lack of high-order correlation treatment and poor scaling. Methods such as DMRG-based extended active-space CASSCF and various Quantum Monte Carlo (QMC) methods will be applied (including VMC/DMC and FCIQMC). Insight into the relative significance of different orbital spaces and the robustness of application of these approaches on leadership class computing architectures will be gained. Synergy with other components of this research program such as automated PES fitting and multireference quantum chemistry will be used to address challenges encountered by the standard approaches to computational thermochemistry (those being single-reference quantum chemistry and perturbative treatments of the anharmonic vibrational energy, which break down for some cases of electronic structure or floppy strongly coupled vibrational modes).

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Evaluating the Structural Response of Amphiphilic Monolayers to Environmental Stimuli

Amphiphilic monolayers composed of end groups with distinct polar and nonpolar functional groups offer rapid and reversible interfacial adaptation in response to environmental stimuli such as a change in interfacial medium polarity. We have synthesized and characterized a suite of monolayers with functional groups of competing polarity designed to reconfigure their interfacial chemical composition in response to solvent polarity. In these films, the end group is designed to be able to reorient and expose the functional groups that minimize the interfacial free energy between the film and the environment. Using a combination of spectroscopic, computational, and wettability studies, we have investigated the responsive interfacial behavior of different end groups upon exposure to environments with varying polarities. Contact angle measurements across a series of polar and dispersive probe liquids reveal trends that reflect the underlying molecular flexibility and composition. Vibrational sum frequency generation (SFG) spectroscopy and atomistic molecular dynamics (MD) simulations confirm solvent-driven reorientation of the end groups, with restructuring observed at the interface. To quantify these effects, we have developed a surface energy calculation model that incorporates solvent-induced surface rearrangements into the estimations. Our findings reveal a strong dependence of surface energy and switching behavior on the length and flexibility of the functionalities in the end group, which affects the exposure of certain interfacial compositions under different solvents. These results offer new insights into the design of adaptive monolayers and provide a framework for evaluating solvent-responsive surfaces.

functional groups↗
DOE OSTI2026

Pressure-dependent kinetics of cyclopentene consumption

Recent experiments have uncovered deficiencies in the ability of current kinetic models for cyclopentane combustion to predict species profiles as a function of O2 concentration. One likely source of these discrepancies flows from the limited availability of relevant kinetic data for key functionalized intermediates, requiring the use of distant analogies to estimate their rates. This work fills this gap for cyclopentene, the most prominent C5 intermediate in low-temperature cyclopentane combustion. Ab initio transition-state-theory based master equation methods are used to calculate pressure-dependent rate constants across four potential energy surfaces (C5H8 + OH, C5H8 + HO2, C5H7, and C5H7 + O2). Potential energy surfaces are characterized at the CCSD(T)-F12/cc-pVDZ-F12//revDSD-PBEP86-D3BJ/def2-TZVPP level of theory. Pressure-dependent rate constants and branching fractions are determined via 1D master equation calculations. The results are used to analyze the competition for radicals between cyclopentene and cyclopentane in the intermediate temperature region, and to identify the dominant pathways flowing from these initiation channels. On the OH initiation surface, we find good agreement with recently published experimental rate constants at high temperature, and we elucidate the complex temperature dependence below 800 K, where 𝜋𝜋 addition becomes prominent. On the HO2 initiation surface, we find an interesting alternative route to bicyclic ether formation via well-skipping over the HO2 adduct well. On the cyclopentenyl radical surface, the major product is cyclopentadiene plus H. On the cyclopentenyl + O2 surface, we find that resonance stabilization largely prevents the allylic radical from undergoing oxidation, whereas the alkylic radical has a somewhat elevated propensity for oxidation pathways relative to cyclopentyl.

ab initio calculations↗
DOE OSTI2025

Potentials of mean force fail to describe chemical bond-breaking in solution

Many liquid phase studies assume that the potential energy surfaces of reacting molecules are the same as in the gas phase, neglecting complex solvent dynamics that can completely alter the nature of chemical reactivity. Even studies that include solvent effects typically only consider them in an average, equilibrium way as part of a potential of mean force (PMF). In this work, we use mixed quantum/classical simulations to compare how equilibrium and non-equilibrium solvent motions affect the photodissociation of a simple diatomic molecule, NaK + , in liquid tetrahydrofuran. A PMF analysis shows that as the excited-state molecule dissociates with the solvent at equilibrium, the bonding electron remains associated with K + at short bond distances but eventually localizes on Na + at the end of dissociation. When we examine non-equilibrium dynamical photodissociation trajectories, however, we find that they fall into three distinct categories: about a quarter of them have the bonding electron mainly associated with Na + , another quarter stay mainly associated with K + , and about half have the bonding electron shared roughly equally between the two ions. The results show that equilibrium PMFs cannot accurately describe the dynamics of bond-breaking chemical reactions in solution because there is insufficient time for the solvent to reach equilibrium on the time scale over which bond dissociation occurs. Furthermore, our analysis shows that the solvent coupling between the electronic energy surfaces is similar at and away from equilibrium, suggesting that other factors, such as solute velocity-driven solvent memory effects, play a more important role in explaining the failure of the equilibrium PMF to predict the non-equilibrium dynamics.

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

Ab-initio simulation of spin-vibronic spectra of methoxy radical

Despite the fact that experimental and theoretical work on the spectrum of methoxy has stretched from the microwave to the ultraviolet and proceeded for nearly 50 years, parts of the spectrum have remained a challenge to simulate theoretically and make reliable line-by-line assignments. The spectral complexity arises because the radical has a non-zero electron spin and significant vibronic coupling between the two elec- tronic components of the ground state due to the presence of a conical intersection. This work describes a completely ab initio effort to understand and assign the spin- vibronic levels of the X 2E state from 0 to above 3000 cm−1, a region that includes the fundamental transitions of the C-H symmetric and asymmetric stretches that have not previously been identified uniquely. A potential energy surface for methoxy was calculated at the EOM-CCSDT/ANO1 level of theory. Subsequently this potential energy surface was fit to a quartic power series expansion of all nine vibrational nor- mal coordinates (as determined at the minimum of the conical intersection) by the use of a machine-learning-based algorithm. After the addition of spin-orbit coupling, the spin-vibronic problem was solved using both the Krylov-Schur and Lanczos algorithms with the SOCJT3 software to converge eigenvalues up to 3500 cm−1 and their eigen- vectors. The latter were used, in conjunction with the calculated dipole moment and its derivatives (calculated using finite differences at EOM-CCSDT/ANO1 level), to determine spectral intensities for the spin-vibronic spectra. The calculated transition frequencies and intensities were used to simulate and assign the observed transitions of the spin-vibronic spectra of the radical. The credibility of the assignments and their significance is discussed in detail.

Sharma, Ketan [University of Florida, Gainesville,↗
DOE OSTI2025

Hyperfine-to-rotational energy transfer in ultracold atom–molecule collisions of Rb and KRb

Energy transfer between different mechanical degrees of freedom in atom–molecule collisions has been studied and largely understood. However, systems involving spins remain less explored. Here, in this study, we directly observed energy transfer from atomic hyperfine to molecular rotation in the 87 Rb ($| {F}_{a},{M}_{{F}_{a}}\rangle =| 2,2\rangle$) + 40 K 87 Rb (X 1 Σ + , rotational state N = 0) ⟶ Rb ($| 1,1\rangle$) + KRb ( N = 0, 1, 2) collision with state-to-state precision. We also performed quantum scattering calculations that rigorously included the coupling between spin and rotational degrees of freedom at short range under the assumption of rigid-rotor KRb monomers moving along a single potential energy surface. The calculated product rotational state distribution deviates from the observations even after extensive tuning of the atom–molecule potential energy surface. In addition, our ab initio calculations indicate that spin–rotation coupling is enhanced close to a conical intersection that is energetically accessible at short range. This, together with the deviation, suggests that vibrational degrees of freedom and conical intersections play an important part in the coupling. Our observations confirm that spin is coupled to mechanical rotation at short range and establish a benchmark for future theoretical studies.

chemical physics↗