Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex solids”

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 1,045 records · Page 58

Autonomous platform for solution processing of electronic polymers

The manipulation of electronic polymers’ solid-state properties through processing is crucial in electronics and energy research. Yet, efficiently processing electronic polymer solutions into thin films with specific properties remains a formidable challenge. We introduce Polybot, an artificial intelligence (AI) driven automated material laboratory designed to autonomously explore processing pathways for achieving high-conductivity, low-defect electronic polymers films. Leveraging importance-guided Bayesian optimization, Polybot efficiently navigates a complex 7-dimensional processing space. In particular, the automated workflow and algorithms effectively explore the search space, mitigate biases, employ statistical methods to ensure data repeatability, and concurrently optimize multiple objectives with precision. The experimental campaign yields scale-up fabrication recipes, producing transparent conductive thin films with averaged conductivity exceeding 4500 S/cm. Feature importance analysis and morphological characterizations reveal key design factors. This work signifies a significant step towards transforming the manufacturing of electronic polymers, highlighting the potential of AI-driven automation in material science.

Wang, Chengshi [Argonne National Laboratory (ANL),↗

Deviatoric stress-induced metallization, layer reconstruction and collapse of van der Waals bonded zirconium disulfide

In contrast to two-dimensional (2D) monolayer materials, van der Waals layered transition metal dichalcogenides exhibit rich polymorphism, making them promising candidates for novel superconductor, topological insulators and electrochemical catalysts. Here, we highlight the role of hydrostatic pressure on the evolution of electronic and crystal structures of layered ZrS2. Under deviatoric stress, our electrical experiments demonstrate a semiconductor-to-metal transition above 30.2 GPa, while quasi-hydrostatic compression postponed the metallization to 38.9 GPa. Both X-ray diffraction and Raman results reveal structural phase transitions different from those under hydrostatic pressure. Under deviatoric stress, ZrS2 rearranges the original ZrS6 octahedra into ZrS8 cuboids at 5.5 GPa, in which the unique cuboids coordination of Zr atoms is thermodynamically metastable. The structure collapses to a partially disordered phase at 17.4 GPa. These complex phase transitions present the importance of deviatoric stress on the highly tunable electronic properties of ZrS2 with possible implications for optoelectronic devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Advancements in NbTiN based circuits for Superconducting Digital Logic

Superconducting (SC) electronics have emerged as a promising platform for high-speed, energy-efficient computing and quantum information processing. This work, centered on NbTiN, presents recent advances in material science and fabrication methods leading to significant improvements in performance, scalability and vertical integration. We specifically report on fabrication and characterization of key components, including Josephson junctions (JJs), flux trapping structures and SC interconnects. Together, these efforts represent critical steps towards realizing practical, complex, dense and large-scale SC integrated circuits.

Pokhrel, A. [Imec,Heverlee,Belgium]↗

Analysis of Planar-Cavity Receiver Reactor for Solar Thermochemical Dry-Reforming

Concentrating solar thermal (CST) systems can be leveraged to provide not only heat for power generation, but also for energy storage and thermochemical fuel production. Such solar thermochemical processes have been studied conceptually, from solar thermochemical hydrogen production (STCH) and thermochemical energy storage (TCES), to gasification, reforming, and fuel upgrading by various means. The solar receiver and reactor are critical components in the conversion of solar energy into chemical energy in the form of 'solar fuels'. For effective conversion of solar energy within a coupled solar receiver-reactor, extremely high temperatures are required, thereby demanding a high solar concentration ratio (CR) for efficient operation. This creates a design challenge for the receiver-reactor, as many thermochemical processes involve gas or gas-solid systems that are limited by low heat transfer coefficients. A unique receiver design is proposed that has the potential to incorporate various high-temperature thermochemical processes such as TCES-assisted power generation, methane reforming, or STCH processes. Modeling this receiver and its potential applications requires a full three-dimensional model that accurately captures the interconnected effects of receiver geometry, spatial solar irradiance, complex radiation, reaction kinetics, fluid dynamics, and heat transfer. In this paper we analyze a CST system integrated with this unique planar-cavity receiver-reactor design using the developed model. The model presented in this paper showed where improved thermal management was needed to achieve suitable receiver performance when a dry-methane reforming process is implemented.

14 SOLAR ENERGY↗

Quarterly Research Performance Progress Report (Q8)

As part of Task 1, we have started by testing our modeling capabilities by reproducing isothermal DFIT simulations presented in the literature. Once satisfied with the results we have started by targeting the modeling of the DFITs at conducted at well 58-32. We have a identified a specific test (cycle 4 in zone 2) as the most interesting to be model with GEOS hydraulic fracturing module. Thus, we have first produced results with an isothermal model and adjusted model parameters to get a satisfying match with field pressure data. The, we have added thermal effects and compared the modeling results with and without thermal effects to estimate how thermal effects may influence test interpretation. Models seem to suggest that, for small volumes of fluid, thermal effects are moderate. In Task 2, we have adapted GEOS phase-field formulation to be able to simulate near-wellbore hydraulic fracture nucleation and propagation. We have devised a novel formulation that, compared to other existing ones, incorporates rock strengths. We have submitted a journal publication about our work. We are currently employing this phase-field formulation to model the experiments taking place at U Pitt and help us understand the effect of various parameters. In Task 3, we have built a model of the region surround well 16A and have started modeling stage 3 stimulation because of its simpler planar geometry. After calibrating simulation parameters using known analytical solutions, we have simulated the stage 3 stimulation using our isothermal hydraulic fracturing module, varying the permeability field, the stress conditions including different physics to get a better understanding of the numerical challenges and of the effects of varying these parameters on the simulation results. In Task 4 laboratory experimentation, a set of specialized drilling and injection tools has been customized and constructed to accommodate an inclined well with an orientation of up to 30 degrees relative to material anisotropy or principal stress axes. These inclined samples have also undergone thermal stress and hydraulic fracturing at a temperature of 190 degrees Celsius. Furthermore, both vertical and deviated sampling testing setups enable an extended analysis of post-peak pressure behaviors, facilitating post-test pressure analyses such as the G-function, step rate, and fracture reopening measurements. Thus, the key components of in-situ stress estimation can be extracted and validated through our experiment, providing a solid foundation for validating existing in-situ stress estimation theories or proposing new ones. Simultaneously, we are integrating computer vision techniques with traditional experimental fracture observation methods such as multi-overcore/slicing and water-penetration fracture observation. This combination will prove beneficial in populating the hydraulic fracture patterns database, generated under challenging EGS conditions. This approach aims to deepen our understanding of the complexities in EGS reservoirs and pave the way for future data-driven investigations. Additionally, PITT has also equipped the ELE International compression machine, which is now prepared for conducting indirect tensile and fracture toughness tests. These tests will aid in characterizing how rock fabrics influence the resulting fracture patterns. Additionally, we have completed the required personnel training and gained access to Scanning Electron Microscopy (SEM) and Energy Dispersive Spectroscopy (EDS) for conducting more detailed characterization and analysis of rock fabrics, as well as the examination of thermal and hydraulically induced fracture patterns. Thus, the PITT team has effectively demonstrated the capabilities of our experimental apparatuses in exploring the thermal effects, well deviation angles, material anisotropy, and operational choices (such as circulation rate, injection fluid viscosity, and injection rate) and their impact on pressure responses and fracture trajectories under the Utah FORGE conditions.

15 GEOTHERMAL ENERGY↗

Mercury Accumulation in a Stream Ecosystem: Linking Labile Mercury in Sediment Porewaters to Bioaccumulative Mercury in Trophic Webs

Mercury (Hg) deposition and accumulation in the abiotic and biotic environments of a stream ecosystem were studied. This study aimed to link labile Hg in porewater to bioaccumulative Hg in biota. Sediment cores, porewaters, and biota were sampled from four sites along the Fourmile Branch (SC, USA) and measured for total Hg (THg) and methyl-Hg (MHg) concentrations. Water quality parameters were also measured at the sediment–water interface (SWI) to model the Hg speciation. In general, Hg concentrations in porewaters and bulk sediment were relatively high, and most of the sediment Hg was in the solid phase as non-labile species. Surface sediment presented higher Hg concentrations than the medium and bottom layers. Mercury methylation and MHg production in the sediment was primarily influenced by sulfate levels, since positive correlations were observed between sulfate and Hg in the porewaters. The majority of Hg species at the SWI were in non-labile form, and the dominant labile Hg species was complexed with dissolved organic carbon. MHg concentrations in the aquatic food web biomagnified with trophic levels (biofilm, invertebrates, and fish), increasing by 3.31 times per trophic level. Based on the derived data, a modified MHg magnification model was established to estimate the Hg bioaccumulation at any trophic level using Hg concentrations in the abiotic environment (i.e., porewater).

54 ENVIRONMENTAL SCIENCES↗

Simulation of Particulate Transport for Delivery of Solid Amendments into the Subsurface: FY24 Status Report

For particulate-based amendments to be viable for field-scale remediation at the Hanford Site (e.g., 200 DV-1 Operable Unit), particles need to be delivered a sufficient radial distance from an injection well and retained at concentrations high enough for effective treatment. An accurate description of the particle radius of influence (ROI) is critical for developing an overall remediation strategy. However, field-scale particle simulations are currently limited due to insufficient simulation capabilities and a lack of experimental data to validate and parameterize particle transport models. To help build toward field-scale deployment, this fiscal year (FY) we have (1) developed a pre screening tool to estimate particle transport, (2) implemented particle transport models within PFLOTRAN, and (3) conducted preliminary estimations of particle ROI. While field-scale numerical simulations will ultimately be necessary before remedy design and field implementation, we have developed a pre-screening tool that offers valuable estimations of expected particle injectability and ROI in a 1-D system. The advantage of the tool is that it does not require extensive laboratory experiments and instead makes predictions based solely on routine laboratory measurements. This tool can assist in down-selection and decision-making by identifying which particle amendment systems are worth pursuing in future laboratory experiments, such as 1-D column tests and beyond. With any system, scaling up from the lab to the field presents challenges. Currently, there is no field data available for model calibration or validation. However, the theoretical particle models being developed herein are the best tools available to guide progress toward field deployment. To help bridge this gap and verify model predictions, larger-scale lab experiments are being proposed. To advance simulation capabilities, six particle transport models are being integrated into the reactive transport simulator PFLOTRAN. These include colloid filtration theory (CFT) and five additional particle transport models (M1-M5). Each model, from M1 to M5, progressively incorporates additional particle transport and retention processes. Ultimately, the simplest model capable of accurately describing 1-D column data will be selected and parameterized. During FY24, the CFT and M1 model have been fully implemented within PFLOTRAN. Using an existing 1 D column experiment, the two currently implemented particle transport models (CFT and M1), and associated parameters, were fit to this experiment. While simpler model formulations are helpful for estimations, these formulations could not fully describe particle transport and retention behavior in the previous 1-D column experiment. Thus, additional complexities will need to be considered, which will be accounted for in the M2-M5 model formulations. Additionally, because a viscous, shear thinning fluid was required to keep particles in suspension, considerations for flow will also need to also be accounted for. Therefore, a new immiscible two-phase flow mode is currently being implemented in PFLOTRAN. With some modifications, this new flow module could also support simulation of non-Newtonian liquid amendments, foams, and emulsions. We also estimated the expected ROI of solid amendments using 1-D simulations. The average predicted ROI was approximately 15 ft for micron-sized zero valent iron (mZVI) suspended in xanthan gum (XG). Using the pre screening tool and ROI estimates, additional amendment-delivery laboratory characterization and experiments are proposed. The results from additional experiments can be used to validate and parametrize particulate transport model formulations, which will ultimately provide predictive capabilities for field amendment-delivery systems.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Reversible solid oxide cells: Early performance and microstructural evolution during electrolysis and switched mode operation

Incorporating more renewables into our energy mix will require a diversity of technological solutions. Reversible solid oxide cell systems (RSOCs) operating reversibly between fuel cell and electrolysis modes with high round-trip efficiencies offer one such solution. However, questions remain about the stability of the oxygen and fuel electrodes under reversible operation. Here, in this study, single cells with the configuration Ni-YSZ|YSZ|GDC10|NNO-NDC50|NNO are operated by reversible cycling under potentiostatic conditions. Test conditions are chosen to be representative of average stack conditions that are likely to be encountered during operation of a stack constructed of such cells. The results of reversible cycling over 500 h (21 cycles) are compared with operating the cells under potentiostatic conditions in the electrolysis-only mode (i.e., without cycling). Voltage-current density (DC), AC complex impedance spectroscopy, distribution of relaxation times (DRT) analysis, and scanning electron microscopy (SEM) indicate excellent reversibility and stability of the neodymium nickelate-based oxygen electrodes in both electrolysis and reversible cycling modes. Degradation behavior is attributed to the loss and coarsening of connected Ni particles in the fuel electrode, particularly in the active layer. However, reversible operation mitigates this degradation compared to operation in the electrolysis-only (SOEC-only) mode.

25 ENERGY STORAGE↗

Radiation Damage Mechanisms in Oxide-Based Solid Electrolytes

Ion implantation is a promising approach to mitigate dendrite formation in metal anode solid-state batteries by modifying the chemical and mechanical properties of oxide-based electrolyte surfaces such as Li 7 La 3 Zr 2 O 12 (LLZO). Yet, the detailed atomistic mechanisms and potentially adverse side effects of ion implantation remain poorly understood. In this study, we used molecular dynamics simulations to elucidate the evolution of radiation-induced damage in LLZO. The results reveal that radiation damage is primarily driven by recoils from heavy ions in the lattice, producing antisite defect clusters whose density decreases as recoil energy increases, a counterintuitive finding that highlights complex cascade fragmentation. Radiation-induced disruption to the crystal sublattice connectivity, especially in the La–O network critical for Li-ion transport, leads to degradation in ionic pathways at low energies, while higher energy cascades can form new conduction pathways via enhanced Zr–O network connectivity. These findings provide a mechanistic foundation for optimizing ion implantation strategies that balance disruption of the lattice for surface modification for dendrite resistance with preservation of essential ionic conductivity.

Defects↗

Alkene adsorption and cracking on acidic zeolites – A gradual process of understanding

Cracking of light alkenes on solid acids and, in particular on zeolites, is ubiquitous in refineries. Different zeolite framework topologies in the catalysts used provide spatial constraints for reactants and products, leading to discrete reaction routes and varying selectivities. The conversion of alkenes in cracking proceeds either via monomolecular or bimolecular reaction pathways. In both routes, the reacting alkene adsorbs first on Brønsted acid sites forming the initial state (e.g., π-complex, carbenium ion or alkoxide), followed by β-cleavage in the carbenium transition state, generating a smaller alkene and a carbenium ion (monomolecular route) or forming a dimer by alkene addition that cracks subsequently (bimolecular route). The current state of understanding alkene sorption and cracking is critically discussed in this paper and four questions are formulated that may be the central questions from our perspective in further developing new generations of catalysts and processes for alkene conversion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

1 H chemical shift anisotropy: a high sensitivity solid-state NMR dynamics probe for surface studies?

We report dynamics play significant roles in chemistry and biochemistry-molecular motions impact both large- and small-scale chemical reactions in addition to biochemical processes. In many systems, including heterogeneous catalysts, the characterization of dynamics remains a challenge. The most common approaches involve the solid-state NMR measurement of anisotropic interactions, in particular 2 H quadrupolar coupling and 1 H–X dipolar coupling, which generally require isotope enrichment. Due to the high sensitivity of 1 H NMR, 1 H chemical shift anisotropy (CSA) is a particularly enticing, and underexplored, dynamics probe. We carried out 1 H CSA and 1 H– 13 C dipolar coupling measurements in a series of model supported complexes to understand how 1 H CSA can be leveraged to gain dynamic information for heterogeneous catalysts. Mathematical descriptions are given for the dynamic averaging of the CSA tensor, and its dependence on orientation and asymmetry. The variability of the orientation of the tensor in the molecular frame, in addition to its magnitude and asymmetry, negatively impacts attempts to extract quantitative dynamic information. Nevertheless, 1 H CSA measurements can reveal useful qualitative insights into the motions of a particularly dilute site, such as from a surface species.

1H chemical shift anisotropy↗

Deciphering the local structure of Prussian blue analogue cathodes with Raman spectroscopy for sodium-ion batteries

Iron-based Prussian blue analogues (PBAs) have gained attention as low-cost, relatively higher-potential cathodes for sodium-ion batteries, due to their open 3D-framework structures. However, understanding the local structural changes is critical to unveiling the intercalation pathways and degradation mechanisms. We employ here operando Raman spectroscopy to probe the changes in the cyanide vibrational modes during cycling and degradation after cycling, which are not adequately resolved by X-ray diffraction due to the low structure factor and limited X-ray sensitivity of the cyanide groups. Vibrational spectroscopy has thus proven essential for deciphering these complex materials. Additionally, we implement pre-sodiation strategies to assess the impact of sodium inventory loss by pairing PBAs with a hard carbon anode in a full-cell configuration. Operando galvanostatic electrochemical impedance spectroscopy (EIS) and ex situ X-ray photoelectron spectroscopy (XPS) further elucidate the interface evolution and the role of water molecules in forming the cathode-electrolyte interphase (CEI) and solid-electrolyte interphase (SEI). The insights gained advance the understanding of PBAs and enhance their practical viability.

25 ENERGY STORAGE↗

High Temperature Erosion In Particle Based CSP Systems

Particle based concentrated solar power plants are increasingly being considered as an alternative to molten salts. Although these particle based systems provide the inherent safety and enhanced lifetime from reduced corrosion compared to molten salt systems, they bring in complexities associated with surface erosion from falling particle impact and sliding motion along the surfaces. Past research on particle related erosion has been limited to high velocity applications and doesn’t necessarily coincide with the operating conditions in CSP systems. In our work, we evaluate characterize the rate of erosion for particles and containment materials at 800ºC. Three different types of erosion resulting from a) impact of solid particles on receiver, particle storage, and heat exchanger walls; b) abrasion erosion from particle sliding motion along walls; and c) attrition erosion as the particles breakdown from particle-to-particle contact and particle-to-wall interactions are evaluated and compared across a varied spectrum of particle and containment material spectrum. From the results, It is noted that temperature plays a key role in the rate of erosion experienced by any material. With the increase in temperature, materials become more ductile and are more readily sheared away from particle impinging or sliding motion. Furthermore, the particles themselves have higher susceptibility to breakdown which also lowers particle thermal absorptance and thus efficiency as heat transfer fluid.

14 SOLAR ENERGY↗

High Temperature Erosion In Particle Based CSP Systems

Particle based concentrated solar power plants are increasingly being considered as an alternative to molten salts. Although these particle based systems provide the inherent safety and enhanced lifetime from reduced corrosion compared to molten salt systems, they bring in complexities associated with surface erosion from falling particle impact and sliding motion along the surfaces. Past research on particle related erosion has been limited to high velocity applications and doesn’t necessarily coincide with the operating conditions in CSP systems. In our work, we evaluate characterize the rate of erosion for particles and containment materials at 800ºC. Three different types of erosion resulting from a) impact of solid particles on receiver, particle storage, and heat exchanger walls; b) abrasion erosion from particle sliding motion along walls; and c) attrition erosion as the particles breakdown from particle-to-particle contact and particle-to-wall interactions are evaluated and compared across a varied spectrum of particle and containment material spectrum. From the results, It is noted that temperature plays a key role in the rate of erosion experienced by any material. With the increase in temperature, materials become more ductile and are more readily sheared away from particle impinging or sliding motion. Furthermore, the particles themselves have higher susceptibility to breakdown which also lowers particle thermal absorptance and thus efficiency as heat transfer fluid.

14 SOLAR ENERGY↗

Post Irradiation Examination Results of Irradiated Yttrium Hydride

Department of Energy’s (DOE’s) Microreactor program (MRP) aims to provide the fundamental data to enable the development of microreactors. As such, material property data of critical materials for microreactor technologies are researched. Because substoichiometric yttrium dihydride (YHx, where x<2) is considered as a potential solid neutron moderator, its material property data has been combined in the Advanced Moderator Material Handbook which includes thermodynamic and thermophysical properties of YHx with the exception of irradiated material’s properties due to limited PIE. To fill the knowledge gap for the irradiated YHx, specimens and irradiation capsules were prepared at Los Alamos National Laboratory (LANL). Specimens were irradiated in the Advanced Test Reactor (ATR) at Idaho National Laboratory’s (INL’s). post-irradiation examination (PIE) was performed at INL’s Materials and Fuels Complex (MFC). This report compiles the PIE results of irradiated YHx specimens through fiscal years 2022 and 2023 (FY22-23) . The PIE data will directly be incorporated into the newer version of the Advanced Moderator Material Handbook. The main takeaways include that (i) the geometrical stability and mechanical integrity of YHx was intact with couple exceptions after high-temperature irradiations (600-800°C), (ii) hydrogen content variation due to manufacturing or irradiation in YHx caused visible surface discoloration, that is also related to the microstructural changes, (iii) qualitative comparisons of PIE methods implied that H retention was significantly higher at 600°C as compared to 800°C, (iv) thermal properties included signatures correlated with the H loss or re-gain, (v) importance of manufacturing readiness and initial as-manufactured specimens history was emphasized, (vi) the needs of targeted irradiations focusing on temperature and time parameters and very targeted PIE were specified.

08 HYDROGEN↗

Report on the Incorporation of Post-Irradiation Examination Results of Yttrium Hydride for Advanced Moderator Handbook

Department of Energy’s (DOE’s) Microreactor program (MRP) aims to provide the fundamental data to enable the development of microreactors. As such, material property data of critical materials for microreactor technologies are researched. Because substoichiometric yttrium dihydride (YHx, where x<2) is considered as a potential solid neutron moderator, its material property data has been combined in the Advanced Moderator Material Handbook which includes thermodynamic and thermophysical properties of YHx with the exception of irradiated material’s properties due to limited PIE. To fill the knowledge gap for the irradiated YHx, specimens and irradiation capsules were prepared at Los Alamos National Laboratory (LANL). Specimens were irradiated in the Advanced Test Reactor (ATR) at Idaho National Laboratory’s (INL’s). post-irradiation examination (PIE) was performed at INL’s Materials and Fuels Complex (MFC). This report compiles the PIE results of irradiated YHx specimens through fiscal years 2022 and 2023 (FY22-23) . The PIE data will directly be incorporated into the newer version of the Advanced Moderator Material Handbook. The main takeaways include that (i) the geometrical stability and mechanical integrity of YHx was intact with couple exceptions after high-temperature irradiations (600-800°C), (ii) hydrogen content variation due to manufacturing or irradiation in YHx caused visible surface discoloration, that is also related to the microstructural changes, (iii) qualitative comparisons of PIE methods implied that H retention was significantly higher at 600°C as compared to 800°C, (iv) thermal properties included signatures correlated with the H loss or re-gain, (v) importance of manufacturing readiness and initial as-manufactured specimens history was emphasized, (vi) the needs of targeted irradiations focusing on temperature and time parameters and very targeted PIE were specified.

08 HYDROGEN↗

Rapid Viscoelastic Deformation Slows Marine Ice Sheet Instability at Pine Island Glacier

The ice sheets of the Amundsen Sea Embayment (ASE) are vulnerable to the marine ice sheet instability (MISI), which could cause irreversible collapse and raise sea levels by over a meter. The uncertain timing and scale of this collapse depend on the complex interaction between ice, ocean, and bedrock dynamics. The mantle beneath the ASE is likely less viscous (~10 18 Pa s) than the Earth's average mantle (~10 21 Pa s). In this paper we show that an effective equilibrium between Pine Island Glacier's retreat and the response of a weak viscoelastic mantle can reduce ice mass lost by almost 30% over 150 years. Other components of solid Earth response—purely elastic deformations and geoid perturbations—provide less stability than the viscoelastic response alone. Uncertainties in mantle rheology, topography, and basal melt affect how much stability we expect, if any. Our study indicates the importance of considering viscoelastic uplift during the rapid retreat associated with MISI.

54 ENVIRONMENTAL SCIENCES↗

Ab initio molecular dynamics with enhanced sampling in heterogeneous catalysis

Ab initio molecular dynamics simulations combined with enhanced sampling techniques are becoming widespread methods to investigate chemical phenomena in catalytic systems. These techniques automatically include finite temperature effects, anharmonicity, and collective dynamics in their robust description of enthalpic and entropic contributions, which can have significant impact on reaction free energy landscapes. This contrasts with standard ab initio static approaches that are based on assessing reaction free energies from various coarse-grained descriptions of the reaction potential energy surface. Enhanced sampling ab initio molecular dynamics opens the way to first principles simulations of systems of increasing complexity like solid/liquid catalytic interfaces. Here, we aim at guiding the reader through the basis of these techniques, summarizing their fundamental theoretical and practical aspects, and reviewing the relevant literature in the field. After a brief introduction to the problem, we will illustrate the advantage of using molecular simulations to include finite temperature effects, examine the most common ab initio techniques currently in use, describe their application to solid state heterogeneous catalysts, and finally critically review the most popular enhanced sampling techniques used in computational catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗