Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phenol”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Developing a new ethylene glycol/H 2 O pretreatment system to achieve efficient enzymatic hydrolysis of sugarcane bagasse cellulose and recover highly active lignin: Countercurrent extraction

Improving pretreatment efficiency is a critical premise in achieving efficient biomass conversion, and obtaining high-performance natural polymers is the guarantee of high-value conversion of biomass. Here, in this study, a new pilot-scale continuous countercurrent pretreatment reaction unit about ethylene glycol-alkali solution was designed for pretreating sugarcane bagasse in order to achieve efficient separation of the three major components of lignocellulose when expanding the scale of pretreatment, reduce lignin deposition on the fiber surface, and obtain highly active lignin and excellent enzymatic hydrolysis efficiency of cellulose. X-ray diffractometer (XRD), X-ray photoelectron spectrometer (XPS), brunauer-emmett-teller (BET) and scanning electron microscope (SEM) methods are used to analyze the structural properties of sugarcane bagasse before and after pretreatment, and high-performance liquid chromatography (HPLC) is used to analyze the monosaccharide components in the enzymatic solution. In addition, the structural properties of the recovered lignin are analyzed by gel permeation chromatography (GPC), 31 P NMR and 2D-HSQC-NMR methods. The results indicate that the system can gain a high cellulose recovery of 92.99% along with a lignin removal of 95.33%, and recovered lignin has low lignin carbohydrate complexes, low condensation, and rich in phenolic hydroxyl groups for 1.95 mmol/g. Meanwhile, the countercurrent pretreatment system can effectively reduce the deposition of lignin on the cellulose surface, which is evidently superior to the non-countercurrent pretreatment and facilitates the efficiency of enzymatic saccharification of substrate, achieving a high glucose yield of 99% as well as a total sugar yield of 91.11%. The method efficiently separates biomass in a green manner, and solid residues are easily hydrolyzed, showing potential for industrial-scale production.

09 BIOMASS FUELS↗
DOE OSTI2020

Antecedent conditions determine the biogeochemical response of coastal soils to seawater exposure

Coastal landscapes are increasingly exposed to seawater due to sea level rise and extreme weather events. The biogeochemical responses of these vulnerable ecosystems are poorly understood, limiting our ability to predict how their role in local and global biogeochemical cycles will shift under future conditions. Here we evaluate how antecedent conditions influence the biogeochemical response of soil to seawater inundation events based on a 42-day laboratory incubation experiment with soils collected from a natural salinity gradient across a coastal floodplain. We quantified influences of seawater inundation on intact soil cores through high-frequency carbon dioxide (CO 2 ) and methane (CH 4 ) gas fluxes measurements as well as ultrahigh resolution characterization of organic matter chemistry and metabolites. Mean CO 2 and CH 4 fluxes were higher after inundation compared to control cores for soils that had low in situ salinity. Soils with low in situ salinity also exhibited significant shifts in organic matter profiles after inundation, with surficial soils (0-7.5 cm) becoming more enriched in phenolic compounds, compared to deeper soils (7.5-15 cm). The number of biochemical transformations inferred from mass spectrometry increased significantly after inundation for soils with low in situ salinity, and a trend of decreasing thermodynamic favorability with increasing salinity was observed, suggesting an increase in microbial activity following inundation. Our results suggest that seawater inundation of low-salinity terrestrial environments can lead to increased microbial activity, preferential depletion of more thermodynamically favorable compounds, and increased soil carbon release. Finally, we conclude that the biogeochemical impacts of future seawater exposure will be modulated by antecedent conditions associated with landscape position within coastal watersheds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Alkylation of monomeric, dimeric, and polymeric lignin models through carbon-hydrogen activation using Ru-catalyzed Murai reaction

In this study, we have assessed directed carbon-hydrogen activation (CHA) for alkylation of monomeric, dimeric, and polymeric lignin models using Murai's catalyst [RuH 2 (CO)(PPh 3 ) 3 ]. Based on related work from our laboratory showing that isolated organosolv lignin bears benzylic directing groups ideal for CHA reactions, this approach could offer new methodology for the valorization of biorefinery lignin. Monomeric and dimeric models bearing a keto group at the benzylic position undergo Ru-catalyzed alkylation in good to excellent yield. Similarly, models bearing a benzylic OH group also undergo alkylation via a tandem oxidation/alkylation process enabled by the Ru catalyst. Polymeric models show low levels of functionalization as a result of the poor solubility of the starting polymer. With unsymmetrical models, functionalization occurs first at the least sterically hindered ortho-site, but a subsequent alkylation, leading to disubstituted products can occur at the more sterically hindered site, leading to hexasubstituted arenes. The reaction shows sensitivity to free phenolic OH groups, which appears to reduce the yield in some reactions, and is also a contributing factor to the low yields observed with polymeric lignin models. Finally, combining CHA methodology with lignin isolation technology able to introduce appropriate directing groups for catalytic functionalization will form the basis for improved conversion of lignin to high value chemical products.

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

Enhancing physicochemical, bioactive, and nutritional properties of sweet potatoes: Ultrasonic contact drying with slot jet nozzles compared to hot-air drying and freeze drying

Sweet potatoes are a rich source of nutrients and bioactive compounds, but their quality can be impacted by the drying process. This study investigates the impact of slot jet reattachment (SJR) nozzle and ultrasound (US) combined drying (SJR + US) on sweet potato quality, compared to freeze-drying (FD), SJR drying, and hot air drying (HAD). SJR + US drying at 50 °C closely resembled FD in enhancing quality attributes and outperformed HAD and SJR in key areas such as rehydration, shrinkage ratios, and nutritional composition. Notably, SJR + US at 50 °C produced the highest total starch (36.84 g/100 g), total dietary fiber (8.48 g/100 g), total phenolic content (158.19 mg GAE/100 g), total flavonoid content (119.08 mg QE/g), DPPH antioxidant activity (6.44 μmol TE/g), β-carotene (31.98 mg/100 g), and vitamin C (5.27 mg/100 g). It also exhibited higher glass transition temperatures (Tg: 14.49 °C), indicating better stability at room temperature. The hardness values for SJR + US samples were similar to FD, while HAD samples had the highest hardness. SJR + US at 50 °C resulted in the lowest total color changes (ΔE), indicating minimal impact on appearance. Additionally, FTIR analysis revealed that peaks in specific spectral regions indicated superior preservation of bioactive compounds in SJR + US samples compared to other methods, which was also confirmed by principal component analysis (PCA) and heatmap visualization. Overall, these findings suggest that SJR + US is an effective alternative to conventional drying techniques, significantly improving the quality of dried sweet potatoes.

Color↗
DOE OSTI2026

Organic colloid composition in variable-redox porewaters within a mountainous floodplain

Redox gradients, often driven by changes in sediment moisture levels in porous, heterogeneous groundwater systems, create dynamic conditions that may promote the production and transport of colloids within natural waters. While much research has focused on the inorganic composition of colloids, the organic composition remains less well understood. Organic matter (OM) in colloids may associate with minerals, complex metal ions, and serve as an electron donor for microbial respiration; therefore, its composition is of high interest. We examined the composition of porewater OM along a redox gradient in a riparian soil located along the Slate River in Crested Butte, Colorado, USA as a function of depth (90, 130, 200, and 350 cm below ground surface). All depths were oxic to suboxic, except 200 cm, where the products of iron and sulfate reduction were observed concomitant with an increase in dissolved and/or colloidal OM, pH, alkalinity, and conductivity. We investigated the composition of porewater using correlated scanning transmission X-ray microscopy and transmission electron microscopy. We observed a change in the OM chemistry from carboxylate-rich at the 200 cm depth (where it was also enmeshed with non-crystalline iron) to phenol- and substituted-aromatic-rich at other depths. Radiocarbon dating revealed carbon in porewater at 200 cm was younger than depths above and below. Soil porewater can flow down into the underlaying gravel bed during baseflow conditions, thus we speculate whether riparian porewater could transport iron- and carboxylate-rich organic colloids into surrounding surface waters through the gravel bed conduit.

Alpine floodplain↗
DOE OSTI2024

Tribological behavior of hybrid Aluminum-TiB 2 metal matrix composites for brake rotor applications

Here, this research investigates the feasibility of hybrid aluminum metal matrix composites (MMC) incorporating TiB 2 particles for brake rotor applications. The composites were produced by incorporating both in-situ, submicron-sized TiB 2 particles and regular micron-sized TiB 2 powders via stir and squeeze casting techniques into A206 aluminum alloy matrix. Systematic adjustments in the fractions of in-situ and ex-situ TiB 2 particles were conducted to evaluate their impact on wear behavior and mechanisms. Combination of both particle types allowed composites with up to 10 vol% of reinforcements. Composites with higher proportions of ex-situ particles demonstrated increased wear resistance compared to those solely composed of in-situ particles, control specimens without TiB 2 , and conventional cast iron counterparts. The lowest wear rate for the hybrid composites sliding against phenolic brake pads was 1.1 × 10 -5 mm 3 /Nm, signifying a 3-fold reduction relative to cast iron sliding against the same pads. Wear analysis elucidated distinctive mechanisms within the hybrid composites, characterized by mild fragmented abrasive wear, adhesive wear, and plastic deformation, accompanied by the formation of an intermixed tribo-oxide layer.

42 ENGINEERING↗
DOE OSTI2023

Expanded Coverage of Phytocompounds by Mass Spectrometry Imaging Using On-Tissue Chemical Derivatization by 4-APEBA

Probing the entirety of any species metabolome is an analytical grand challenge, especially at a cellular scale. Where spatial metabolomics, completed primarily by matrix-assisted laser desorption/ionization (MALDI), has limited molecular coverage for several reasons. To expand the scope of spatial metabolomics, we developed an on-tissue chemical derivatization (OTCD) workflow using 4-APEBA for confident identification of several dozen elusive phytocompounds, including several phytohormones, which have various roles within stress responses and cellular communication. Superiority of 4-APEBA is established in comparison to other derivatization agents with (1) broad specificity towards carbonyls, (2) low background, and (3) introduction of bromine isotopes, where the latter two facilitate confident bioinformatics. In conclusion, the outlined workflow trailblazes a path towards spatial hormonomics within plant samples, enhancing detection of carboxylates, aldehydes, ketones, and plausibly phenols.

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

Co-Hydrotreatment of Yellow Greases and the Water-Insoluble Fraction of Pyrolysis Oil. Part I: Experimental Design to Increase Kerosene Yield and Reduce Coke Formation

Here, this paper reports the co-hydrotreatment of Biomass Technology Group commercial pyrolysis oil water-insoluble (WIS) phase (also known as pyrolytic lignin) and yellow greases (waste cooking oil), aiming to produce sustainable aviation fuels (SAFs). We use a sulfided NiMo/Al 2 O 3 blend with 16 wt % WIS and a central composite experimental design to identify processing conditions increasing kerosene yield and reducing coke formation. The input variables were: (1) reaction temperature (320, 350, and 380 °C), (2) initial hydrogen pressure (5, 6, and 7 MPa), and (3) amount of catalyst (0.7, 1.0, and 1.3 g). The hydrotreated oily phases were distilled to obtain gasoline (<150 °C), kerosene (150–250 °C), diesel (250–350 °C), and residual oil (>350 °C). The reaction temperature is the main factor affecting the yield of gaseous, solid, and liquid products. Meanwhile, a higher initial hydrogen pressure and catalyst loading increased the yield of kerosene and other distillates and decreased the coke formation. A high temperature correlated with a lower content of oxygenates in kerosene cuts. Based on our experimental results, we propose to conduct hydrotreatment studies at 380 °C, initial H 2 pressure of 7 MPa, and 1.3 g of catalyst. Under the identified conditions, it was possible to improve the kerosene yield to more than 20 wt % and reduce the yield of coke to close to 2.0 wt %. The chemical composition and fuel properties of the gasoline, kerosene, and diesel cuts were thoroughly analyzed. The content of aromatics and phenols in the kerosene fraction produced at the conditions identified in this project exceeded the recommended values for SAFs. New strategies (such as blending and more intense hydrotreatment to remove oxygenated compounds) need to be implemented to reduce the content of these molecules in our final product.

09 BIOMASS FUELS↗
DOE OSTI2023

Elucidating Biomass-Derived Pyrolytic Lignin Structures from Demethylation Reactions through Density Functional Theory Calculations

Pyrolytic lignin is a fraction of pyrolysis oil that contains a wide range of phenolic compounds that can be used as intermediates to produce fuels and chemicals. However, the characteristics of the raw lignin structure make it difficult to establish a pyrolysis mechanism and determine pyrolytic lignin structures. Herein this study proposes dimer, trimer, and tetramer structures based on their relative thermodynamic stability for a hardwood lignin model in pyrolysis. Different configurations of oligomers were evaluated by varying the positions of the guaiacyl (G) and syringyl (S) units and the bonds βO4 and β5 in the hardwood model lignin through electronic structure calculations. The homolytic cleavage of βO4 bonds is assumed to occur and generate two free radical fragments. These can stabilize by taking hydrogen radicals that may be in solution during the intermediate liquid (pathway 1) formation before the thermal ejection. An alternative pathway (pathway 2) could occur when the radicals use intramolecular hydrogen, turning themselves into stable products. Subsequently, a demethylation reaction can take place, thus generating a methane molecule and new oligomeric lignin-derived molecules. The most probable resulting structures were studied. We used FTIR and NMR spectra of selected model compounds to evaluate our calculation approach. Thermophysical properties were calculated using group contribution methods. The results give insights into the lignin oligomer structures and how these molecules are formed. They also provide helpful information for the design of pyrolysis oil separation and upgrading equipment.

09 BIOMASS FUELS↗
DOE OSTI2024

Liquid–Liquid Equilibrium Prediction in Fast Pyrolysis Bio-Oil Systems: A Framework for Incorporating Bio-Oil Complexity

The study of mixtures of bio-oil, water and organic solvents in different proportions can serve as a cost-effective analysis of its content due to the formation of immiscible phases. This manuscript attempts to replicate experimentally determined partition coefficients (K OW ) of relevant species present in fast pyrolysis bio-oil (FPBO). A commercial flowsheeting simulator with surrogate bio-oil model representation is used. Concurrently, pyrolytic lignins in FPBO (‘pyrolignin’) do not have an agreed-upon structural representation, and the literature is ripe with wide variations of said representations. Thus, during the description of FPBO, this pyrolignin fraction was modeled using 20 possible structures (phenolic dimers to tetramers), with the goal of determining the structures for which the experimental data are best described. Two cases were considered: Case 1 normalized the reported experimental mass balance, while Case 2 included the unreported fraction in the mass balance to the total pyroligin. Please, add here a comment on the prediction of the Water oil equilibrium. The best KOW predictions for levoglucosan (LVG) were obtained when the system was modeled with no pyrolignin, presenting an MRE under 10% for both systems WO and BO. Among the possible structures, D2 (dimer), F1 (trimer) and I1, and I3 (tetramers) presented MRE ≤ 13% for both cases.

09 BIOMASS FUELS↗
DOE OSTI2021

Elucidating the Role of Single-Atom Pd for Electrocatalytic Hydrodechlorination

In this study, we loaded Pd catalysts onto a reduced graphene oxide (rGO) support in an atomically dispersed fashion [i.e., Pd single-atom catalysts (SACs) on rGO or Pd 1 /rGO] via a facile and scalable synthesis based on anchor-site and photoreduction techniques. The as-synthesized Pd 1 /rGO significantly outperformed the Pd nanoparticle (Pdnano) counterparts in the electrocatalytic hydrodechlorination of chlorinated phenols. Downsizing Pd nano to Pd 1 leads to a substantially higher Pd atomic efficiency (14 times that of Pdnano), remarkably reducing the cost for practical applications. The unique single-atom architecture of Pd 1 additionally affects the desorption energy of the intermediate, suppressing the catalyst poisoning by Cl – , which is a prevalent challenge with Pdnano. Characterization and experimental results demonstrate that the superior performance of Pd 1 /rGO originates from (1) enhanced interfacial electron transfer through Pd–O bonds due to the electronic metal–support interaction and (2) increased atomic H (H*) utilization efficiency by inhibiting H 2 evolution on Pd 1 . Furthermore, this work presents an important example of how the unique geometric and electronic structure of SACs can tune their catalytic performance toward beneficial use in environmental remediation applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Molecular Characterization of Water-Soluble Brown Carbon Chromophores in Snowpack from Northern Xinjiang, China

This study reports molecular-level characterization of brown carbon (BrC) attributed to water-soluble organic carbon in six snowpack samples collected from northern Xinjiang, China. Here, the molecular composition and light-absorbing properties of BrC chromophores were unraveled by application of high-performance liquid chromatography (HPLC) coupled to a photodiode array (PDA) detector and high-resolution mass spectrometry. The chromophores were classified into five major types, that is, (1) phenolic/lignin-derived compounds, (2) flavonoids, (3) nitroaromatics, (4) oxygenated aromatics, and (5) other chromophores. Identified chromophores account for ~23–64% of the total light absorption measured by the PDA detector in the wavelength range of 300–370 nm. In the representative samples from urban and remote areas, oxygenated aromatics and nitroaromatics dominate the absorption in the wavelengths below and above 320 nm, respectively. The highly polluted urban sample shows the most complex HPLC-PDA chromatogram, and more other chromophores contribute to the bulk absorption. Phenolic/ lignin-derived compounds are the most light-absorbing species in the soil-influenced sample. Chromophores in two remote samples exhibit ultraviolet–visible features distinct from other samples, which are attributed to flavonoids. Identification of individual chromophores and quantitative analysis of their optical properties are helpful for elucidating the roles of BrC in snow radiative balance and photochemistry.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Effect of Redox-Inactive Metal Ion–Nickel(III) Interactions on the Redox Properties and Proton-Coupled Electron Transfer Reactivity

We report Mononuclear nickel(II) and nickel(III) complexes of a bisamidate-bisalkoxide ligand, (NMe 4 ) 2 [Ni- II (HMPAB)] (1) and (NMe 4 )[Ni III (HMPAB)] (2), respectively, have been synthesized and characterized by various spectroscopic techniques including X-ray crystallography. The reaction of redox-inactive metal ions (M n+ = Ca 2+ , Mg 2+ , Zn 2+ , Y 3+ , and Sc 3+ ) with 2 resulted in 2-M n+ adducts, which was assessed by an array of spectroscopic techniques including X-ray absorption spectroscopy (XAS), electron paramagnetic resonance (EPR), and reactivity studies. The X-ray structure of Ca 2+ coordinated to Ni(III) complexes, 2Ca 2+ T, was determined and exhibited an average Ni-Ca distance of 3.1253 angstrom, close to the metal ions' covalent radius. XAS analysis of 2-Ca 2+ and 2-Y 3+ in solution further revealed an additional coordination to Ca and Y in the 2-M n+ adducts with shortened Ni-M distances of 2.15 and 2.11 Å, respectively, implying direct bonding interactions between Ni and Lewis acids (LAs). Such a short interatomic distance between Ni(III) and M is unprecedented and was not observed before. EPR analysis of 2 and 2-M n+ species, moreover, displayed rhombic signals with g av > 2.12 for all complexes, supporting the +III oxidation state of Ni. The Ni III /Ni II redox potential of 2 and 2-M n+ species was determined, and a plot of E 1/ 2 of 2-M n+ versus pK a of [M(H 2 O) n ] m+ exhibited a linear relationship, implying that the Ni III /Ni II potential of 2 can be tuned with different redox-inactive metal ions. Reactivity studies of 2 and 2-M n+ with different 4-X-2,6-ditert-butylphenol (4-XDTBP) and other phenol derivatives were performed, and based on kinetic studies, we propose the involvement of a proton-coupled electron transfer (PCET) pathway. Analysis of the reaction products after the reaction of 2 with 4-OMe-DTBP showed the formation of a Ni(II) complex (1a) where one of the alkoxide arms of the ligand is protonated. A pK a value of 24.2 was estimated for 1a. The reaction of 2-M n+ species was examined with 4-OMe-DTBP, and it was observed that the k 2 values of 2-M n+ species increase by increasing the Lewis acidity of redox-inactive metal ions. However, the obtained k 2 values for 2-M n+ species are much lower compared to the k 2 value for 2. Such a variation of PCET reactivity between 2 and 2-M n+ species may be attributed to the interactions between Ni(III) and LAs. Our findings show the significance of the secondary coordination sphere effect on the PCET reactivity of Ni(III) complexes and furnish important insights into the reaction mechanism involving high-valent nickel species, which are frequently invoked as key intermediates in Ni-mediated enzymatic reactions, solar-fuel catalysis, and biomimetic/synthetic transformation reactions.

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

Distinguishing Desirable and Undesirable Reactions in Multicomponent Systems for Redox Activation of the Uranyl Ion

Although it has been established that covalent functionalization of the U–O bonds in the uranyl dication (UO 2 2+ ) generally requires use of strong reductants and electrophiles, little work has examined how interactions between the individual reaction components could affect final outcomes in solution. Here, the patterns of such reactivity have been studied in a UO 2 2+ -containing model system supported by a workhorse pentadentate ligand, 2,2′-[(methylimino)bis(2,1-ethanediylnitrilomethylidyne)]bis-phenol. Oxo activation and functionalization have been tested with (i) electrochemical and chemical reduction, and (ii) coordinating and noncoordinating solvents. In acetonitrile, uranyl reduction was achieved cleanly, but treatment of the reduced species with tris(pentafluorophenyl)borane (BCF) resulted in a mixture of products arising from direct electron transfer to BCF. In dichloromethane (CH 2 Cl 2 ), electrochemical reduction of uranyl was achieved cleanly, but clean chemical reactivity was inaccessible. Despite these challenges, one trinuclear and oxo-deficient uranium-containing product was crystallized from CH 2 Cl 2 solution and characterized; thus, desirable electrophilic reactivity can proceed to some degree in CH 2 Cl 2 with BCF. Computational studies were used to investigate the properties of the trinuclear uranium product and the changes that could be inducible by further reduction. Here, taken together, the reactivity patterns identified here could inform design of improved systems for actinyl oxo functionalization.

cations↗
DOE OSTI2022

At-Line Sampling and Characterization of Pyrolytic Vapors from Biomass Feedstock Blends Using SPME-GC/MS-PCA: Influence of Char on Fast Pyrolysis

Solid-phase microextraction (SPME) coupled with GC-MS analysis was used for at-line sampling of pyrolytic vapors produced during fast pyrolysis of biomass. The pure and binary blends of switchgrass (SWG) and pine harvest residues (PT6) were used as feedstock. Sequential SPME sampling allowed for monitoring of changes in the pyrolysis vapors as char accumulated in the fluid bed. The concentration and composition of the vapors desorbed from the SPME fibers were investigated using GC-MS, and the data sets were then analyzed using principal component analysis (PCA) to compare the composition of the pyrolysis vapors over the course of the pyrolysis run. The chemical composition of both carbohydrate and lignin fragments varied as the char builds up in the reactor bed. Fragments derived from cellulose and xylan included anhydrosugars, furans and light oxygenated compounds. Lignin fragments included methoxyphenols, phenolic ketones and aldehydes, low molecular weight aromatics. The composition of the carbohydrate fragments changed more than the lignin fragments as the char build-up in the fluid bed. This combination of SPME-GC/MS-PCA were a novel, easy and effective method for measuring the composition and changes in the composition of pyrolysis vapors during the fast pyrolysis process. Here, this work also highlighted the effect of char build-up on the composition of the overall pyrolysis vapors.

09 BIOMASS FUELS↗
DOE OSTI2023

Nuclear–Electronic Orbital QM/MM Approach: Geometry Optimizations and Molecular Dynamics

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods allow simulations of chemical reactions in atomistic solvent and heterogeneous environments such as proteins. Herein, the nuclear–electronic orbital (NEO) QM/MM approach is introduced to enable the quantization of specified nuclei, typically protons, in the QM region using a method such as NEO-density functional theory (NEO-DFT). This approach includes proton delocalization, polarization, anharmonicity, and zero-point energy in geometry optimizations and dynamics. Expressions for the energies and analytical gradients associated with the NEO-QM/MM method, as well as the previously developed polarizable continuum model (NEO-PCM), are provided. Geometry optimizations of small organic molecules hydrogen bonded to water in either dielectric continuum solvent or explicit atomistic solvent illustrate that aqueous solvation can strengthen hydrogen-bonding interactions for the systems studied, as indicated by shorter intermolecular distances at the hydrogen-bond interface. We then performed a real-time direct dynamics simulation of a phenol molecule in explicit water using the NEO-QM/MM method. Furthermore, these developments and initial examples provide the foundation for future studies of nuclear–electronic quantum dynamics in complex chemical and biological environments.

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

Analytic Nuclear Gradients for Complete Active Space Linearized Pair-Density Functional Theory

Accurately modeling photochemical reactions is difficult due to the presence of conical intersections and locally avoided crossings, as well as the inherently multiconfigurational character of excited states. As such, one needs a multistate method that incorporates state interaction in order to accurately model the potential energy surface at all nuclear coordinates. The recently developed linearized pair-density functional theory (L-PDFT) is a multistate extension of multiconfiguration PDFT, and it has been shown to be a cost-effective post-MCSCF method (as compared to more traditional and expensive multireference many-body perturbation methods or multireference configuration interaction methods) that can accurately model potential energy surfaces in regions of strong nuclear–electronic coupling in addition to accurately predicting Franck–Condon vertical excitations. Here, in this paper, we report the derivation of analytic gradients for L-PDFT and their implementation in the PySCF-forge software, and we illustrate the utility of these gradients for predicting ground- and excited-state equilibrium geometries and adiabatic excitation energies for formaldehyde, s-trans-butadiene, phenol, and cytosine.

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

Hybrid Quantum Mechanical, Molecular Mechanical, and Machine Learning Potential for Computing Aqueous-Phase Adsorption Free Energies on Metal Surfaces

Performing reliable computer simulations of elementary processes occurring at metal–water interfaces is pivotal for novel catalyst design in sustainable energy applications. Computational catalyst design hinges on the ability to reliably and efficiently compute the potential energy surface (PES) of the system. Here, due to the large system sizes needed for studying processes at liquid water–metal interfaces, these systems can currently not be described using density functional theory (DFT). In this work, we used a hybrid quantum mechanical, molecular mechanical, and machine learning potential for studying the adsorption behavior of phenol, atomic hydrogen, 2-butanol, and 2-butanone on the (0001) facet of Ru under reducing conditions when Ru is not oxidized. Specifically, we describe the adsorbate and the surrounding metal atoms at the DFT level of theory. Here, we also considered the electrostatic field effect of the water molecules on adsorbate–metal interactions. Next, for the water–water and water–adsorbate interactions, we used established classical force fields. Finally, for the water–Ru surface interaction, for which no reliable force fields have been published, we used Behler–Parrinello high-dimensional neural network potentials (HDNNPs). Employing this setup, we used our explicit solvation for metal surface (eSMS) approach to compute the aqueous-phase effect on the low-coverage adsorption of selected molecules and atoms on the (0001) facet of Ru. In agreement with previous experimental and computational studies of oxygenated molecules over transition metal facets, we found that liquid water destabilizes the tested adsorbates on Ru(0001). Interestingly, our findings indicate that adsorbates on Ru are less affected by the presence of an aqueous phase than on other transition metals (e.g., Pt), highlighting the necessity of experimental investigations of Ru-based catalytic systems in liquid water.

Adsorption↗