Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecules potential 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 145 records · Page 8

Diffusion Monte Carlo Study of the Structure and Spectroscopy of H 3 O –

A potential energy surface for H 3 O – has been developed based on the NN+(MOB-ML) approach we developed for studies of complexes of OH – with two or three water molecules. Unlike those systems, H 3 O – has two low-energy isomers, H – ·H 2 O and OH – ·H 2 , which differ in energy by less than 2.5 kcal mol –1 , and which are separated by a barrier of roughly 4.5 kcal mol –1 . We find that by training the NN+(MOB-ML) model using structures based on diffusion Monte Carlo (DMC) simulations initiated in the two potential minima, we are able to obtain a potential surface that describes both isomers. Using these potentials, the structure and spectra of H 3 O – and its deuterated analogues are investigated using DMC. These calculations show that the ground state wave function for H 3 O – is mainly localized in the H – ·H 2 O minimum in the potential, with a small amount of the probability amplitude (<5%) in the region of the OH – ·H 2 minimum. The delocalization of the wave function into the secondary minimum is lowered by deuteration of the water molecule, while replacing H – with D – increases the isomerization due to the shortening of the average distance between the hydride ion and the hydrogen atom in water to which it is bound. Introducing one quantum of excitation in the H – ···H 2 O stretching vibration increases the amount of isomerization, while the isomerization decreases with additional excitation of this mode. Excitation of the free OH stretch in the water molecule also increases the amount of isomerization, while excitation of the out-of-plane bending vibration suppresses the isomerization. Furthermore, the effects of partial deuteration on the frequencies for these vibrations are also explored.

Chemical structure↗

w19_OMEC Scientific Highlight [Slides]

Scientific Achievement: Revealed the effect of graphite substrate on H binding / hydrogen evolution reaction (HER) activity of organic molecular electrocatalysts (OMECs); and demonstrated the effect of molecular modification on calculated HER activity descriptor. Significance and Impact: Understanding of HER processes is important for discovering low-cost OMECs for hydrogen production and energy application; design of OMEC structures remains challenging due to the unexplored redox properties of organic molecules; and density functional theory (DFT) has been proposed to provide needed insights into HER activity for improved materials design. Research Details: DFT calculations on the binding energy of HER intermediates are used for predicting reduction potentials and and calculated H adsorption energy as activity descriptor; and LANL Institutional Computing (IC) resources are essential for successful project execution given the size of structures involved in modeling molecule-support interaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Double-Atom Catalysts Featuring Inverse Sandwich Structure for CO 2 Reduction Reaction: A Synergetic First-Principles and Machine Learning Investigation

Electrocatalytic CO 2 reduction reactions (CO 2 RR) based on scalable and highly efficient catalysis provide an attractive strategy for reducing CO 2 emissions. Here in this work, we combined first-principles density functional theory (DFT) and machine learning (ML) to comprehensively explore the potential of double-atom catalysts (DACs) featuring an inverse sandwich structure anchored on defective graphene (gra) to catalyze CO 2 RR to generate C 1 products. We started with five homonuclear M 2 ⊥gra (M = Co, Ni, Rh, Ir, and Pt), followed by 127 heteronuclear MM'⊥gra (M = Co, Ni, Rh, Ir, and Pt, M' = Sc–Au). Stable DACs were screened by evaluating their binding energy, formation energy, and dissolution potential of metal atoms, as well as conducting first-principles molecular dynamics simulations with and without solvent water molecules. Based on DFT calculations, Rh 2 ⊥gra DAC was found to outperform the other four homonuclear DACs and the Rh-based single- and double-atom catalysts of noninverse sandwich structures. Out of the 127 heteronuclear DACs, 14 were found to be stable and have good catalytic performance. An ML approach was adopted to correlate key factors with the activity and stability of the DACs, including the sum of radii of metal and ligand atoms (d M–M' , d M–C , and d M'–C ), the sum and difference of electronegativity of two metal atoms (P M + P M' , P M – P M '), the sum and difference of first ionization energy of two metal atoms (I M + I M' , I M – I M '), the sum and difference of electron affinity of two metal atoms (A M + A M' , A M – A M '), and the number of d-electrons of the two metal atoms (Nd). The obtained ML models were further used to predict 154 potential electrocatalysts out of 784 possible DACs featuring the same inverse sandwich configuration. Overall, this work not only identified promising CO 2 RR DACs featuring the reported inverse sandwich structure but also provided insights into key atomic characteristics associated with high CO 2 RR activity.

30 DIRECT ENERGY CONVERSION↗

Ab initio quantum scattering calculations and a new potential energy surface for the HCl( X 1Σ+)–O2(X3Σg−) system: Collision-induced line shape parameters for O2-perturbed R(0) 0–0 line in H35Cl

The remote sensing of abundance and properties of HCl—the main atmospheric reservoir of Cl atoms that directly participate in ozone depletion—is important for monitoring the partitioning of chlorine between “ozone-depleting” and “reservoir” species. Such remote studies require knowledge of the shapes of molecular resonances of HCl, which are perturbed by collisions with the molecules of the surrounding air. In this work, we report the first fully quantum calculations of collisional perturbations of the shape of a pure rotational line in H35Cl perturbed by an air-relevant molecule [as the first model system we choose the R(0) line in HCl perturbed by O2]. The calculations are performed on our new highly accurate HCl(X1Σ+)–O2(X3Σg−) potential energy surface. In addition to pressure broadening and shift, we also determine their speed dependencies and the complex Dicke parameter. This gives important input to the community discussion on the physical meaning of the complex Dicke parameter and its relevance for atmospheric spectra (previously, the complex Dicke parameter for such systems was mainly determined from phenomenological fits to experimental spectra and the physical meaning of its value in that context is questionable). We also calculate the temperature dependence of the line shape parameters and obtain agreement with the available experimental data. We estimate the total combined uncertainties of our calculations at 2% relative root-mean-square error in the simulated line shape at 296 K. This result constitutes an important step toward computational population of spectroscopic databases with accurate ab initio line shape parameters for molecular systems of terrestrial atmospheric importance.

Chemistry↗

Signatures of Non-universal Quantum Dynamics of Ultracold Chemical Reactions of Polar Alkali Dimer Molecules with Alkali Metal Atoms: Li( 2 S) + NaLi( a 3 Σ + ) → Na( 2 S) + Li 2 ( a 3 Σ u + )

Ultracold chemical reactions of weakly bound triplet-state alkali metal dimer molecules have recently attracted much experimental interest. Here we perform rigorous quantum scattering calculations with a new ab initio potential energy surface to explore the chemical reaction of spin-polarized NaLi(a 3 Σ + ) and Li( 2 S) to form Li 2 (a 3 Σ u + ) and Na( 2 S). The reaction is exothermic and proceeds readily at ultralow temperatures. Significantly, we observe strong sensitivity of the total reaction rate to small variations of the three-body part of the Li 2 Na interaction at short range, which we attribute to a relatively small number of open Li 2 (a 3 Σ u + ) product channels populated in the reaction. This provides the first signature of highly non-universal dynamics seen in rigorous quantum reactive scattering calculations of an ultracold exothermic insertion reaction involving a polar alkali dimer molecule, opening up the possibility of probing microscopic interactions in atom+molecule collision complexes via ultracold reactive scattering experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational studies reveal Fluorine based quinolines to be potent inhibitors for proteins involved in SARS-CoV-2 assembly

World is witnessing one of the worst pandemics of this century caused by SARS-CoV-2 virus which has affected millions of individuals. Despite rapid efforts to develop vaccines and drugs for COVID-19, the disease is still not under control. Chloroquine (CQ) and Hydroxychloroquine (HCQ) are two very promising inhibitors which have shown positive effect in combating the disease in preliminary experimental studies, but their use was reduced due to severe side-effects. Here, we performed a theoretical investigation of the same by studying the binding of the molecules with SARS-COV-2 Spike protein, the complex formed by Spike and ACE2 human receptor and a human serine protease TMPRSS2 which aids in cleavage of the Spike protein to initiate the viral activation in the body. Both the molecules had shown very good docking energies in the range of -6kcal/mol. Subsequently, we did a high throughput screening for other potential quinoline candidates which could be used as inhibitors. From the large pool of ligand candidates, we shortlisted the top three ligands (binding energy -8kcal/mol). We tested the stability of the docked complexes by running Molecular Dynamics (MD) simulations where we observed the stability of the quinoline analogues with the Spike-ACE2 and TMPRSS2 nevertheless the quinolines were not stable with the Spike protein alone. Thus, although the inhibitors bond very well with the protein molecules their intrinsic binding affinity depends on the protein dynamics. Moreover, the quinolines were stable when bound to electronegative pockets of Spike-ACE2 or TMPRSS2 but not with Viral Spike protein. We also observed that a Fluoride based compound: 3-[3-(Trifluoromethyl)phenyl]quinoline helps the inhibitor to bind with both Spike-ACE2 and TMPRSS2 with equal probability. The molecular details presented in this study would be very useful for developing quinoline based drugs for COVID-19 treatment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Supramolecular Self-Assembled Multi-Electron-Acceptor Organic Molecule as High-Performance Cathode Material for Li-Ion Batteries

Organic electrode materials possess many advantages such as low toxicity, sustainability, and chemical/structural tunability toward high energy density. However, to compete with inorganic-based compounds, crucial aspects such as redox potential, capacity, cycling stability, and electronic conductivity need to be improved. Herein, a comprehensive strategy on the molecular design of small organic electron-acceptor-molecule—hexaazatrianthranylene (HATA) embedded quinone (HATAQ) is reported. By introducing conjugated quinone moieties into the electron-deficient hexaazatriphenylene-derivative core, HATAQ with highly extended π-conjugation can yield extra-high capacity for lithium storage, delivering a capacity of 426 mAh g -1 at 200 mA g -1 (0.4C). At an extremely high rate of 10 A g -1 (19C), a reversible capacity of 209 mAh g -1 corresponding to nearly 85% retention is obtained after 1000 cycles. A unique network of unconventional lock-and-key hydrogen bonds in the solid-state facilitates favorable supramolecular 2D layered arrangement, enhancing cycling stability. To the best of the authors’ knowledge, the capacity and rate capability of HATAQ are found to be the best ever reported for organic small-molecule-based cathodes. Finally, these results together with density functional theory studies provide proof-of-concept that the design strategy is promising for the development of organic electrodes with exceptionally high energy density, rate capability, and cycling stability.

25 ENERGY STORAGE↗

Advancing technologies for lignin-based jet fuel production in aqueous phase

Integrating lignin into a cellulosic ethanol plant for the co-production of lignin-based jet fuel (LJF) in aqueous phase offers a significant opportunity to boost operational efficiency, economic viability, carbon conversion, and the overall sustainability of biofuel and chemical production. LJF is lignin-structure-based jet fuel blendstocks primarily composed of alkyl-substituted mono-, bi-, and tri-cyclohexanes. It exhibits high energy density, potential for low emissions, and favourable blend characteristics that comply with drop-in specifications. An overview of lignin feedstock, catalytic processes, LJF chemical compositions, fuel properties tests, and techno-economic analysis (TEA) and life cycle assessment (LCA) indicate that (1) the reactivity of lignin plays a crucial role in its structure transformation to LJF molecules; (2) catalytic processing of lignin to LJF can occur through a simultaneous depolymerization and hydrodeoxygenation process, bypassing the intermediate step of producing and upgrading lignin-derived oil; (3) LJF's uniqueness molecules making it more promising for high energy content and low emission jet fuel properties for next generation sustainable aviation fuel (SAF); and (4) TEA and LCA demonstrate that LJF is not only potentially cost-effective but also offers favourable carbon footprint compared to other SAFs. In conclusion, this review highlights the most recent advancements in LJF technology, along with the challenges and opportunities that lie ahead in fulfilling its potential.

09 BIOMASS FUELS↗

Molecular-orbital-based machine learning for open-shell and multi-reference systems with kernel addition Gaussian process regression

We introduce a novel machine learning strategy, kernel addition Gaussian process regression (KA-GPR), in molecular-orbital-based machine learning (MOB-ML) to learn the total correlation energies of general electronic structure theories for closed- and open-shell systems by introducing a machine learning strategy. The learning efficiency of MOB-ML(KA-GPR) is the same as the original MOB-ML method for the smallest criegee molecule, which is a closed-shell molecule with multi-reference characters. In addition, the prediction accuracies of different small free radicals could reach the chemical accuracy of 1 kcal/mol by training on one example structure. Accurate potential energy surfaces for the H10 chain (closed-shell) and water OH bond dissociation (open-shell) could also be generated by MOB-ML(KA-GPR). To explore the breadth of chemical systems that KA-GPR can describe, we further apply MOB-ML to accurately predict the large benchmark datasets for closed- (QM9, QM7b-T, and GDB-13-T) and open-shell (QMSpin) molecules.

Chemistry↗

A Stable, Low Permeable TEMPO Catholyte for Aqueous Total Organic Redox Flow Batteries

We report aqueous organic redox flow batteries (AORFBs) are highly attractive for large-scale energy storage because of their nonflammability, low cost, and sustainability. (2,2,6,6-Tetramethylpiperidin-1-yl)oxyl (TEMPO) derivatives, a class of redox active molecules bearing air-stable free nitroxyl radicals and high redox potential (>0.8 V vs NHE), has been identified as promising catholytes for AORFBs. However, reported TEMPO based molecules are either permeable through ion exchange membranes or not chemically stable enough for long-term energy storage. Herein, a new TEMPO derivative functionalized with a dual-ammonium dicationic group, N 1 , N 1 , N 1 , N 3 , N 3 , 2, 2, 6, 6-nonamethyl-N 3 -(piperidinyloxy)propane-1,3-bis(ammonium) dichloride (N 2 -TEMPO) as a stable, low permeable catholyte for AORFBs is reported. Ultraviolet–visible (UV–vis) and proton nuclear magnetic resonance ( 1 H-NMR) spectroscopic studies reveal its exceptional stability and ultra-low permeability (1.49 × 10 -12 cm 2 s -1 ). Coupled with 1,1'-bis[3-(trimethylammonio)propyl]-4,4'-bipyridinium tetrachloride ((NPr) 2 V) as an anolyte, a 1.35 V N2-TEMPO/(NPr)2V AORFB with 0.5 m electrolytes (9.05 Wh L -1 ) delivers a high power density of 114 mW cm -2 and 100% capacity retention for 400 cycles at 60 mA cm-2. At 1.0 m electrolyte concentrations, the N 2 -TEMPO/(NPr) 2 V AORFB achieves an energy density of 18.1 Wh L -1 and capacity retention of 90% for 400 cycles at 60 mA cm -2 .

TEMPO↗

Adsorption of ethylenediamine on Cu surfaces: attributes of a successful capping molecule using first-principles calculations

The shape-controlled synthesis of Cu nanocrystals can benefit a wide range of applications, though challenges exist in achieving high and selective yields to a particular shape. Capping agents play a pivotal role in controlling shape, but their exact role remains ambiguous. In this study, the adsorption of ethylenediamine (EDA) on Cu(100) and Cu(111) was investigated with quantum density functional theory (DFT) to reveal the complex roles of EDA in promoting penta-twinned Cu nanowire growth. We find EDA has stronger binding on Cu(100) than on Cu(111), which agrees the general expectation that penta-twinned Cu nanowires express facets with stronger capping-molecule binding. Despite this stronger binding, ab initio thermodynamics reveals the surface energy of EDA-covered Cu(111) is lower than that EDA-covered Cu(100) at all solution-phase EDA chemical potentials, so there is no thermodynamic driving force for penta-twinned nanowires. We also investigated the capability of EDA to protect Cu surfaces from oxidation in water by quantifying energy barriers for a water molecule to diffuse through EDA layers on Cu(100) and Cu(111). The energy barrier on Cu(100) is significantly lower, which supports observations of faster oxidation of Cu(100) in electrochemical experiments. Thus, we elucidate another possible function of a capping agent – to enable selective oxidation of crystal facets. This finding adds to the general understanding of successful attributes of capping agents for shape-selective nanocrystal growth.

36 MATERIALS SCIENCE↗

A Donor–Acceptor [2]Catenane for Visible Light Photocatalysis

Colored charge-transfer complexes can be formed by the association between electron-rich donor and electron-deficient acceptor molecules, bringing about the narrowing of HOMO–LUMO energy gaps so that they become capable of harnessing visible light. In an effort to facilitate the use of these widespread, but nonetheless weak, interactions for visible light photocatalysis, it is important to render the interactions strong and robust. Furthermore, we employ a well-known donor–acceptor [2]catenane—formed by the mechanical interlocking of cyclobis(paraquat-p-phenylene) and 1,5-dinaphtho[38]crown-10—in which the charge-transfer interactions between two 4,4'-bipyridinium and two 1,5-dioxynaphthalene units are enhanced by mechanical bonding, leading to increased absorption of visible light, even at low concentrations in solution. As a result, since this [2]catenane can generate persistent bipyridinium radical cations under continuous visible-light irradiation without the need for additional photosensitizers, it can display good catalytic activity in both photo-reductions and -oxidations, as demonstrated by hydrogen production—in the presence of platinum nanoparticles—and aerobic oxidation of organic sulfides, such as l-methionine, respectively. This research, which highlights the usefulness of nanoconfinement present in mechanically interlocked molecules for the reinforcement of weak interactions, can not only expand the potential of charge-transfer interactions in solar energy conversion and synthetic photocatalysis but also open up new possibilities for the development of active artificial molecular shuttles, switches, and machines.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Plasmon Dynamics Driven by Aggregation of Tris(2,2′-bipyridine)ruthenium(II)-Functionalized Gold Nanoparticles Probed by XANES and Transient Absorption Spectroscopy

Energy conversion dynamics is critical for advancing next-generation photovoltaics, optoelectronics, and light-harvesting technologies. Noble metal plasmonic nanoparticles play a pivotal role as nanoscale electromagnetic confinement structures, driving photon-induced chemical reactions. In this study, we explore the effects of [Ru(bpy) 3 ] 2+ functionalization and aggregation on citrate-capped gold nanoparticles (AuNPs) of 40 and 100 nm diameters, focusing on molecule-plasmon interactions and their influence on electronic and energy dissipation properties. X-ray absorption near-edge spectroscopy (XANES) revealed that [Ru(bpy) 3 ] 2+ functionalization induces controlled aggregation without altering the oxidation state of gold. A more pronounced white-line intensity is observed in 40 nm AuNPs, consistent with greater s–p–d hybridization and a higher density of surface states, likely influenced by both nanoparticle size and aggregation. Transient absorption (TA) spectroscopy highlights faster electron–phonon relaxation dynamics in aggregated 40 nm nanoparticles, which is attributed to increased electron delocalization and more efficient coupling to the phonon bath. In contrast, 100 nm nanoparticles exhibit minimal changes due to a lower degree of aggregation. Interestingly, we observe that enhanced electron–phonon coupling in aggregated nanoparticles coincides with a slowing of electron–electron scattering. These observations suggest a competitive interplay between the two relaxation pathways, where enhanced energy transfer to the lattice in aggregated systems can suppress electronic thermalization. In conclusion, these findings underscore the critical role of nanoparticle size, aggregation, and molecule–surface interactions in modulating plasmonic dynamics and excited-state lifetimes and further provide valuable insights into designing tailored plasmonic systems with transformative potential for sensing, catalysis, and energy conversion.

36 MATERIALS SCIENCE↗

Transgenic Sugarcane–Oilcane: An Alternative Feedstock for the Production of Drop-in Fuel and Value-Added Bioproducts

The utilization of plants and other agricultural produce can partly offset petroleum dependency for energy requirements and can potentially provide sustainable solutions to global environmental problems. To this end, synthetic biology has shown great potential in developing transgenic bioenergy grasses such as sugarcane, sorghum, miscanthus, and energy cane that hyperaccumulate energy-rich lipid molecules in their vegetative tissues, such as leaves, stems, and roots. These perennial high-biomass transgenic C4 grasses are not targeted to grow on prime agricultural land and can be dedicatedly used to produce biofuels and other value-added bioproducts. Recently, sugarcane has been metabolically engineered to sequestrate carbon from juice towards biosynthesis of lipid molecules in the vegetative tissues. Transgenic sugarcane is referred to as “oilcane.” Transgenic sugarcane–oilcane has improved energy density due to an elevated lipid content in the vegetative tissues. Transgenic bioenergy crops cater to both cellulosic sugars and vegetative lipids, hence resulting in higher biofuel yield (biodiesel/renewable diesel and bioethanol) per unit area of cultivable land. In this book chapter, the genetic engineering of sugarcane and its bioprocessing are discussed to illustrate its development and use as an alternative feedstock for the production of biofuels (bio-jet fuel, biodiesel, and bioethanol) and value-added bioproducts.

Maitra, Shraddha↗

Range-separated hybrid functionals for mixed dimensional heterojunctions: Application to phthalocyanines/MoS 2

We analyze the electronic structure and level alignment of transition-metal phthalocyanine (MPc) molecules adsorbed on two-dimensional MoS 2 employing density functional theory (DFT) calculations. We develop a procedure for multi-objective optimal tuning of parameters of range-separated hybrid functionals in these mixed-dimensional systems. Using this procedure, which leads to the asymptotically-correct exchange-correlation potential between molecule and two-dimensional material, we obtain electronic structures consistent with experimental photoemission results for both energy level alignment and electronic bandgaps, representing a significant advance compared to standard DFT methods. We elucidate the MoS 2 valence resonance with the transition-metal phthalocyanine non-frontier 3d orbitals and its dependence on the transition metal atomic number. Based on our calculations, we derive parameter-free, model self-energy corrections that quantitatively accounts for the effects of the heterogeneous dielectric environment on the electronic structure of these mixed-dimensional heterojunctions.

2D materials↗

Multielectron Organic Redoxmers for Energy-Dense Redox Flow Batteries

Redox flow battery is a highly promising stationary energy storage method but the limited energy density and high chemical cost are among the main barriers for commercialization. Multi-electron organic redoxmers represent a family of structurally tailorable candidates that can achieve multiplied energy density with decreased materials consumption, potentially resulting in a viable solution to address these challenges. In this work, the recent development of organic molecules with reversible multi-redox activities in both aqueous and nonaqueous electrolytes is reviewed. The major focus is on the fundamental correlation between the chemical structures and the functional properties of reported multi-electron organic molecules. Valuable insights are offered on rational structural design strategies for improving the relevant physico- and electro-chemical properties. Finally, the current challenges are discussed to suggest future research needs along the avenue of using the multi-electron approach to achieve energy-dense, stable, cost-effective redox flow batteries.

25 ENERGY STORAGE↗

Donor End-Capped Alkyl Chain Length Dependent Non-Radiative Energy Loss in All-Small-Molecule Organic Solar Cells

A critical bottleneck for further efficiency breakthroughs in organic solar cells (OSCs) is to minimize the non-radiative energy loss (eΔV nr ) while maximizing the charge generation. With the development of highly emissive low-bandgap non-fullerene acceptors, the design of high-performance donors becomes critical to enable the blend with the electroluminescence quantum efficiency to approach or surpass the pristine acceptor. Herein, by shortening the end-capped alkyl chains of the small-molecular donors from hexyl (MPhS-C6) to ethyl (MPhS-C2), the material obtained aggregation that was insensitive to thermal annealing (TA) along with condensed packing simultaneously. The former leads to small phase separation and suppressed upshifts of the highest occupied molecular orbital energy level during TA, and the latter facilitates its efficient charge-transport at aggregation-less packing. Hence, the ΔV nr decreases from 0.242 to 0.182 V, from MPhS-C6 to MPhS-C2 based OSCs. An excellent PCE of 17.11% is obtained by 1,8-diiodoctane addition due to almost unchanged high J sc (26.6 mA cm –2 ) and V oc (0.888 V) with improved fill factor, which is the record efficiency with the smallest energy loss (0.497 eV) and ΔV nr (0.192 V) in all-small-molecule OSCs. Finally, these results emphasize the potential material design direction of obtaining concurrent TA-insensitive aggregation and condensed packing to maximize the device performances with a super low ΔV nr .

36 MATERIALS SCIENCE↗

Effect of Collective Dynamics and Anharmonicity on Entropy in Heterogenous Catalysis: Building the Case for Advanced Molecular Simulations

In this paper, we present a perspective on the computational determination of entropy and its effects and consequences on heterogeneous catalysis. Special attention is paid to the role of anharmonicity and the deviations from the standard harmonic oscillator approximations which can fail to provide a reliable assessment of entropy. To address these challenges, advanced methodologies are needed that can explicitly account for these thermodynamic drivers through the appropriate statistical sampling of reactive free energy surfaces. We discuss where anharmonicity should be expected, where it has been observed from a theoretical perspective, and the methods currently employed to address it. We concentrate on three types of systems where we have observed major, nonnegligible anharmonic effects: (1) supported nanoparticles, where the migration of metal atoms, complexes, and entire clusters exhibit anharmonic behavior in their dynamic motion; (2) porous solids where confinement effects distort potential energy surfaces and hinder molecular motions resulting in large entropic terms; and (3) solid/liquid interfaces, where interactions between solvent molecules and adsorbed species can result in large solvent organization free energy and unique reactivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗