Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dissociation”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Cis-Divacant Octahedral Fe(II) in a Dimensionally Reduced Family of 2-(Pyridin-2-yl)pyrrolide Complexes

Four-coordinate transition metal complexes can adopt a diverse array of coordination geometries, with square planar and tetrahedral coordination being the most prevalent. Previously, we reported the synthesis of a trinuclear Fe(II) complex, Fe 3 TPM 2 , supported by a three-fold symmetric 2-pyridylpyrrolide ligand (i.e., tris(5-(pyridin-2-yl)-1H-pyrrol-2-yl)methane), that featured a rare cis-divacant octahedral (CDO) geometry at each Fe(II) center. Here, a series of truncated 2-pyridylpyrrolide ligands is described that support mono- and binuclear Fe(II) complexes that also exhibit CDO geometries. Metallation of tetradentate ligand bis(5-(pyridin-2-yl)-1H-pyrrol-2-yl)methane (H 2 BPM) in THF results in a binuclear complex Fe 2 (BPM) 2 (THF) 2 in which both Fe(II) ions are octahedrally coordinated. The coordinated THF solvent ligands are labile: THF dissociation leads to Fe 2 (BPM) 2 , which features five-coordinate Fe(II) ions. The Fe–Fe distance in these binuclear complexes can be elongated by ligand methylation. Metalation of bis(5-(6-methylpyridin-2-yl)-1H-pyrrol-2-yl)methane (H 2 BPM Me ) in THF leads to the formation of four-coordinate, CDO Fe(II) centers in Fe(BPM Me ) 2 . Further ligand truncation affords bidentate ligands 2-(1H-pyrrol-2-yl)pyridine (PyrPyrrH) and 2-methyl-6-(1H-pyrrol-2-yl)pyridine (Pyr Me PyrrH). Metalation of these ligands in THF affords six-coordinate complexes Fe(PyrPyrr) 2 (THF) 2 and Fe(Pyr Me Pyrr) 2 (THF) 2 . Dissociation of labile solvent ligands provides access to fourcoordinate Fe(II) complexes. Ligand disproportionation at Fe(PyrPyrr) 2 results in the formation of Fe(PyrPyrr) 3 and Fe(0). Ligand methylation suppresses this disproportionation and enables isolation of Fe(Pyr Me Pyrr) 2 , which is rigorously CDO. Complete ligand truncation, by separating the 2-pyridylpyrrolide ligands into the constituent monodentate pyridyl and pyrrolide donors, affords Fe(Pyr) 2 (Pyrr) 2 in which the Fe(II) is tetrahedrally coordinated. Computational analysis indicates that the potential energy surface that dictates the coordination geometry in this family of four-coordinate complexes is fairly flat in the vicinity of CDO coordination. Furthermore, these synthetic studies provide the structural basis to explore the implications of CDO geometry on Fe-catalyzed reactions.

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

Cobalt(II) Phthalocyanine Substituents Tune the Electrocatalytic CO 2 Conversion to Methanol

Cobalt phthalocyanine (Co(II)Pc) and its derivatives are promising molecular electrocatalysts for the electrochemical reduction of CO 2 to CO and methanol (CH 3 OH). Despite increasing interest, a detailed mechanistic understanding of how ligand substituents influence catalytic activity, selectivity, and efficiency remains limited. In this study, we employ density functional theory (DFT) to systematically investigate the influence of electron-donating groups (EDGs) and electron-withdrawing groups (EWGs) on the electronic structure and redox properties of the Co(II)Pc electrocatalyst, and to elucidate CO 2 RR mechanistic pathways. Our results reveal that EWGs cause a positive shift in the reduction potentials, favor CO 2 binding over protonation of the Co metal center and promote downstream methanol formation at mild potentials. EDGs show opposite trends including favorable protonation steps, promoting a negative shift in the reduction potential, and facilitating the hydrogen evolution reaction (HER), which competes with the desired CO 2 RR pathway. Notably, CO dissociation is thermodynamically and kinetically unfavorable across all systems, positioning the redox potential versus CO dissociation energy as a key factor for methanol selectivity. Furthermore, these insights provide a predictive framework for rational catalyst design and underscore the critical role of electronic tuning in advancing molecular electrocatalysts for sustainable CO 2 conversion.

Alcohols↗
DOE OSTI2025

Synthesis, Properties, and Electrochemical Proton Reduction of a Homoleptic Tetrathiolato Ni-Site Model of [NiFe]-Hydrogenase

[NiFe]-hydrogenase enzymes process H 2 at a nonplanar tetracysteinato-Ni site, the sole participator in proton binding/redox chemistry during turnover. With the objective of assessing whether a simple tetrahedral/ tetrathiolato-Ni 2+ core could promote H 2 evolution reaction (HER), we synthesized (Et 4 N) 2 [Ni(S-p-CF 3 −Ph) 4 ] (1) employing para-trifluoromethylbenzenethiolate ( − S-p-CF 3 −Ph) as a Ni-site analog of [NiFe]-hydrogenase. Spectroscopic measurements and X-ray crystallography confirm the distorted tetrahedral geometry of 1. Dissolution of 1 results in partial thiolate dissociation and formation of S,S-bridged complexes such as (Et 4 N) 2 [Ni 2 (S-p-CF 3 −Ph) 6 ] (3) among other ill-defined species. Dissociation is further accelerated in the presence of Brønsted acids, complicating the assessment of 1 for proton reduction. However, this dissociation/proton instability is suppressed in the presence of additional thiolate ligand to ensure tetrahedral/tetrathiolato 1 persists in solution. Electrochemical HER activity was evaluated by monitoring the current response of an MeCN solution of 1/ excess thiolate after sequential titration with a weak Brønsted acid (acetic acid). The results suggest that 1 is a modest electrocatalyst for the HER with a turnover frequency of 14.5 ± 3.6 s −1 and an overpotential of 0.72 ± 0.02 V. Control experiments and supplementary DFT computations indicate that 1, or a species derived from 1, is responsible for the HER and suggest an ECCE-type mechanism.

Evolution reactions↗
DOE OSTI2024

QC-GN 2 oMS 2 : a Graph Neural Net for High Resolution Mass Spectra Prediction

Predicting the mass spectrum of a molecular ion is often accomplished via three generalized approaches: rules-based methods for bond breaking, deep learning, or quantum chemical (QC) modeling. Rules-based approaches are often limited by the conditions for different chemical subspaces and perform poorly under chemical regimes with few defined rules. QC modeling is theoretically robust but requires significant amounts of computational time to produce a spectrum for a given target. Among deep learning techniques, graph neural networks (GNNs) have performed better than previous work with fingerprint-based neural networks in mass spectra prediction. To explore this technique further, we investigate the effects of including quantum chemically derived information as edge features in the GNN to increase predictive accuracy. The models we investigated include categorical bond order, bond force constants derived from extended tight-binding (xTB) quantum chemistry, and acyclic bond dissociation energies. Throughout this work, we evaluated these models against a control GNN with no edge features in the input graphs. Bond dissociation enthalpies yielded the best improvement with a cosine similarity score of 0.462 relative to the baseline model (0.437). In this work we also apply dynamic graph attention which improves performance on benchmark problems and supports the inclusion of edge features. Between implementations, we investigate the nature of the molecular embedding for spectra prediction and discuss the recognition of fragment topographies in distinct chemistries for further development in tandem mass spectrometry prediction.

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

Including Physics-Informed Atomization Constraints in Neural Networks for Reactive Chemistry

Machine learning interatomic potentials (MLIPs) have emerged as powerful tools for investigating atomistic systems with high accuracy and a relatively low computational cost. However, a common and unaddressed challenge with many current neural network (NN) MLIP models is their limited ability to accurately predict the relative energies of systems containing isolated or nearly isolated atoms, which appear in various reactive processes. To address this limitation, we present a mathematical technique for modifying any existing atom-centered NN architecture to account for the energies of isolated atoms. The result produces a consistent prediction of the atomization energy (AE) of a system using minimal constraints on the model. Using this technique, we build a model architecture that we call hierarchically interacting particle neural network (HIP-NN)-AE, an AE-constrained version of the HIP-NN, as well as ANI-AE, the AE-constrained version of the accurate NN engine for molecular energies (ANI). Our results demonstrate AE consistency of AE-constrained models, which drastically improves the AE predictions for the models. We compare the AE-constrained approach to unconstrained models as well as models from the literature in other scenarios, such as bond dissociation energies, bond dissociation pathways, and extensibility tests. These results show that the constraints improve the model performance in some of these tasks and do not negatively affect the performance on any tasks. The AE constraint approach thus offers a robust solution to the challenges posed by isolated atoms in energy prediction tasks.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Stability Analysis of Substituted Cobaltocenium [Bis(cyclopentadienyl)cobalt(III)] Employing Chemistry-Informed Neural Networks

Cobaltocenium derivatives are promising components of the anion exchange membranes due to their excellent thermal and alkaline stability under the operating conditions of a fuel cell. Here we present an efficient modeling approach of assessing the chemical stability of substituted cobaltocenium CoCp 2 + based on the computed electronic structure enhanced by machine learning techniques. Within the aqueous environment, the positive charge of the metal cation is balanced by the hydroxide anion through formation of the CoCp 2 + OH¯ complexes, whose dissociation is studied within the implicit solvent employing density functional theory. The data set of about 118 species based on 42 substituent groups characterized by a range of electron- donating (ED) and electron-withdrawing (EW) properties is constructed and analyzed. Given 12 carefully chosen chemistry-informed descriptors of the complexes and relevant fragments, the stability of the complexes is found to strongly correlate with the energies of the highest occupied and lowest unoccupied molecular orbitals, modulated by a switching function of the Hirshfeld charge. The latter is used as a measure of the electron withdrawing-donating character of the substituents. Based on this observation from the conventional regression analysis, two fully connected, feed-forward neural network (FNN) models with different unit structures, called the chemistry-informed (CINN) and the quadratic (QNN) neural networks, are developed. Both models predict the bond dissociation energies of the cobaltocenium complexes with mean relative errors less than 5.40% and average absolute errors less than 0.94 kcal/mol. The results show the potential of QNN to efficiently capture more complex relationships. Here, the concept of incorporating the domain (chemical) knowledge/insight into the neural network structure paves the way to applications of machine learning techniques with small data sets, ultimately leading to better predictive models compared to the conventional regression analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Improving Excited-State Potential Energy Surfaces via Optimal Orbital Shapes

In this work, we demonstrate that, rather than resorting to high-cost dynamic correlation methods, qualitative failures in excited-state potential energy surface predictions can often be remedied at no additional cost by ensuring that optimal molecular orbitals are used for each individual excited state. This approach also avoids the weighting choices required by state-averaging and dynamic weighting and obviates their need for expensive wave function response calculations when relaxing excited-state geometries. Although multistate approaches are of course preferred near conical intersections, other features of excited-state potential energy surfaces can benefit significantly from our single-state approach. In three different systems, including a double bond dissociation, a biologically relevant amino hydrogen dissociation, and an amino-to-ring intramolecular charge transfer, we show that state-specific orbitals offer qualitative improvements over the state-averaged status quo.

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

Formation of Phenanthrene via Recombination of Indenyl and Cyclopentadienyl Radicals: A Theoretical Study

This work presents quantum chemical G3(MP2,CC)//B2PLYPD3/6-311G(d,p) calculations of the potential energy surface for the indenyl (C 9 H 7 ) + cyclopentadienyl (C 5 H 5 ) reaction followed by unimolecular decomposition of the C 14 H 11 radicals formed as the primary products, as well as the Rice–Ramsperger–Kassel–Marcus master equation (RRKM-ME) calculations to predict temperature- and pressure-dependent reaction rate constants and product branching ratios. The reaction begins with the barrierless recombination of indenyl and cyclopentadienyl forming a C 14 H 12 molecule with a new C–C bond connecting two five-membered rings, which subsequently dissociates to C 14 H 11 radicals by H losses. The primary products of the C 9 H 7 + C 5 H 5 → C 14 H 11 + H reaction can directly decompose by another H loss to benzofulvalene, and this pathway is most favorable in terms of the entropy factor and hence is preferable at higher temperatures. Otherwise, the initial C 14 H 11 isomers can undergo significant structural rearrangements before eliminating an H atom and producing phenanthrene, anthracene, or benzoazulenes, among which the formation of phenanthrene via the “spiran” pathway is clearly preferred. The calculated barriers along the computed favorable dissociation pathways are relatively low, in the ~30–40 kcal/mol range, making the C 14 H 11 radicals unstable at temperatures above 1000–1250 K at 1 atm. The results of RRKM-ME calculations show that, under typical combustion conditions, the decomposition of the C 14 H 11 radicals predominantly leads to benzofulvalene. However, the latter can be rapidly converted to phenanthrene via H-assisted isomerization with the rate constant for the benzofulvalene + H → phenanthrene + H reaction being close to 10 –11 cm 3 molecule –1 s –1 at 1000–1500 K and 1 atm. Here, the results provide further support for the hypothesis that recombination of two π radicals containing five-membered rings can lead to a growth of PAH with the formation of two fused six-membered rings, but the reaction mechanism may not be direct and is likely to involve two consecutive H atom losses leading to a fulvalene-like product, with subsequent H-assisted isomerization of the latter to a benzenoid PAH.

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

Characterization of Uranyl Coordinated by Equatorial Oxygen: Oxo in UO 3 versus Oxyl in UO 3 +

Uranium trioxide, UO 3 , has a T-shaped structure with bent uranyl, UO 2 2+ , coordinated by an equatorial oxo, O 2– . The structure of cation UO 3 + is similar but with an equatorial oxyl, O •– . Neutral and cationic uranium trioxide coordinated by nitrates were characterized by collision induced dissociation (CID), infrared multiple-photon dissociation (IRMPD) spectroscopy, and density functional theory. CID of uranyl nitrate, [UO 2 (NO 3 ) 3 ] – (complex A1), eliminates NO 2 to produce nitrate-coordinated UO 3 + , [UO 2 (O • )(NO 3 ) 2 ]– (B1), which ejects NO 3 to yield UO 3 in [UO 2 (O)(NO 3 )] – (C1). Finally, C1 associates with H 2 O to afford uranyl hydroxide in [UO 2 (OH) 2 (NO 3 )] – (D1). IRMPD of B1, C1, and D1 confirms uranyl equatorially coordinated by nitrate(s) along with the following ligands: (B1) radical oxyl O •– ; (C1) oxo O 2– ; and (D1) two hydroxyls, OH – . As the nitrates are bidentate, the equatorial coordination is six in A1, five in B1, four in D1, and three in C1. Ligand congestion in low-coordinate C1 suggests orbital-directed bonding. Hydrolysis of the equatorial oxo in C1 epitomizes the inverse trans influence in UO 3 , which is uranyl with inert axial oxos and a reactive equatorial oxo. The uranyl ν 3 IR frequencies indicate the following donor ordering: O 2– [best donor] >> O •– > OH – > NO 3 – .

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

Pressure Effects on the Relaxation of an Excited Ethane Molecule in High-Pressure Bath Gases

Here, we use molecular dynamics to calculate the rotational and vibrational energy relaxation of C 2 H 6 in Ar, Kr, and Xe bath gases over a pressure range of 10 to 400 atm and at temperatures of 300 K and 800 K. The C 2 H 6 is instantaneously excited by 80 kcal/mol randomly distributed into both vibrational and rotational modes. The computed relaxation rates show little sensitivity to the identity of the noble gas in the bath. Vibrational relaxation rates show a non-linear pressure dependence at 300 K. At 800 K the reduced range of bath gas densities covered by the range of pressures do not yet show any non-linearity in the pressure dependence. Rotational relaxation is characterized with two relaxation rates. The slower rate is comparable to the vibrational relaxation rate. The faster rate has a linear pressure dependence at 300 K but an irregular, nonlinear pressure dependence at 800 K. To understand this, a model was developed based on approximating the periodic box used in the molecular dynamics simulations by an equal-volume collection of cubes where each cube is sized to allow only single occupancy by the noble gas or the molecule. Combinatorial statistics then leads to a pressure and temperature dependent analytic distribution of the bath gas species the molecule encounters in a collision. This distribution, the dissociation energy of molecule/bath gas complexes and bath gas clusters, and the computed energy release per collision combine to show that only at 300 K is the energy release sufficient to dissociate likely complexes and clusters. This suggests that persistent and pressure-dependent clusters and complexes at 800 K may be responsible for the non-linear pressure dependence of rotational relaxation.

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

Thermodynamics of the Metal Carbonates and Bicarbonates of Mn, Co, Ni, Cu, and Zn Relevant to Mineral Energetics

The gas phase heats of formation of ground-state MCO 3 , M(HCO 3 ) 2 , and M(HCO 3 )(OH), where M = Mn, Co, Ni, Cu, and Zn, have been predicted using the correlated molecular orbital theory at the CCSD(T) level extrapolated to the complete basis set limit using the Feller–Peterson–Dixon (FPD) approach. Cohesive energies of the carbonates were predicted based on the calculated gas phase and experimental solid heats of formation. Coulombic dissociation energies (CDEs) between metal cations and anions show a near-linear correlation with Shannon metal cation atomic radii, yet no correlation is found with the hardness of these cations. The total reaction dissociation energies (TRDEs) of transition metals are higher than their CDEs for the di-bicarbonates, in contrast to those for Mg and Ca based on our prior work. In addition to differences in the energies needed to prepare the transition metal dications, electron donation from the ligands to the 3d orbitals of open-shell transition metal dications from lone pairs of adjacent O atoms also plays a role. No electron donation from the ligands to the fully occupied 3d orbitals of Zn and Cd was found. Decomposition energies for generating MO, CO 2 , and/or H 2 O were calculated. Gas phase metal exchange energies only partially correlate with the electrochemical series for M(s) → M 2+( aq). The FPD heats of formation were used to benchmark a range of density functional theory exchange–correlation functionals, including those commonly used in solid-state mineral calculations. None of the functionals provided chemical accuracy agreement (±1 kcal/mol) with the FPD results. So, the best agreement with the FPD results is predicted for the $\tau$-HCTH functional with an average unsigned error of 8.3 kcal/mol.

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

Computational Predictions of the Hydrolysis of 2,4,6-Trinitrotoluene (TNT) and 2,4-Dinitroanisole (DNAN)

Hydrolysis is a common transformation reaction that can affect the environmental fate of many organic compounds. In this study, three proposed mechanisms of alkaline hydrolysis of 2,4,6-trinitrotoluene (TNT) and 2,4-dinitroaniline (DNAN) were investigated with plane-wave density functional theory (DFT) combined with ab initio and classical molecular dynamics (AIMD/MM) free energy simulations, Gaussian basis set DFT calculations, and correlated molecular orbital theory calculations. Most of the computations in this study were carried out using the Arrows web based tools. For each mechanism: Meisenheimer complex formation, nucleophilic aromatic substitution, and proton abstraction reaction energies and activation barriers were calculated for the reaction at each relevant site. For TNT, it was found that the most kinetically favorable first hydrolysis steps involve Meisenheimer complex formation by attachment of OH$^{-}$ at the C1 and C3 arene carbons and proton abstraction from the methyl group. The nucleophilic aromatic substitution reactions at the C2 and C4 arene carbons were found to be thermodynamically favorable. However, the calculated activation barriers were slightly lower than previous studies, but still found to be $\Delta$G$^{\ddagger} \approx 18$ kcal/mol using PBEo AIMD/MM free energy simulations, suggesting the reactions are not kinetically significant. For DNAN, the barriers of nucleophilic aromatic substitution were even greater ($\Delta$G$^{\ddagger} > 29$ kcal/mol PBEo AIMD/MM). Further, the most favorable hydrolysis reaction for DNAN was found to be a two-step process in which the hydroxyl first attacks the C1 carbon to form a Meisenheimer complex at the C1 arene carbon C1-(OCH$_3$)OH$^{-}$, subsequently, the methoxy anion (-OCH$_3$) at the C1 arene carbon dissociates and the proton shuttles from the C1-OH to the dissociated methoxy group, resulting in methanol and an aryloxy anion.

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

Association Kinetics for Perfluorinated n -Alkyl Radicals

Radical-radical reaction channels are important in the pyrolysis and oxidation chemistry of perfluoroalkyl substances (PFAS). In particular, unimolecular dissociation reactions within unbranched n-perfluoroalkyl chains, and their corresponding reverse barrierless association reactions, are expected to be significant contributors to the gas-phase thermal decomposition of families of species such as perfluorinated carboxylic acids and perfluorinated sulfonic acids. Unfortunately, experimental data for these reactions are scarce and uncertain. Furthermore, obtaining reliable theoretical predictions for such reactions is a laborious and computationally intensive task. Here, in this work, the chemical kinetics of the various association/decomposition reactions producing/decomposing the C 2 -C 4 series of unbranched n-perfluoroalkanes (C 2 F 6 , C 3 F 8 , and C 4 F 10 ) are examined using state-of-the-art ab initio transition-state-theory-based master-equation calculations. The variable-reaction-coordinate transition-state theory (VRC-TST) formalism is employed in computing the microcanonical and canonical rates for the association reactions. Reaction thermochemistry is obtained via composite quantum chemistry calculations and the laddering of error-canceling reaction schemes via a connectivity-based hierarchy approach employing ANL1/ANL0-style reference energies. Lennard-Jones collision model parameters for the considered systems were estimated by a direct dynamics approach, and collisional energy transfer parameters were obtained from analogies to systems of similar size and heavy-atom connectivity. A one-dimensional master equation approach was used to convert the microcanonical rate coefficients from the VRC-TST analysis into temperature- and pressure-dependent rate constants for the association reactions and the reverse dissociation reactions. The data are reported in standardized formats for usage in comprehensive chemical kinetic models for PFAS thermal destruction.

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

Quantifying the Decomposition Kinetics of Linear C 2 –C 4 Perfluoroalkanes

The pyrolysis of C 2 –C 4 linear perfluoroalkanes was investigated using laser schlieren densitometry (LS) in a diaphragmless shock tube. Mixtures of 1, 2, and 4% perfluoroethane (PFE), perfluoropropane (PFP), and perfluorobutane (PFB) in excess krypton were shock-heated to 1400–2500 K at pressures of 60, 120, and 240 Torr. The initiation step for all perfluoroalkanes involved carbon–carbon bond fission. The rate constants for unimolecular decomposition were computed using microcanonical rate theory, and these rate constants were combined into a detailed chemical kinetic mechanism. For PFE and PFP, there was one dissociation channel, while PFB dissociated via two competing bond-fission reactions. The perfluoroalkyl radicals subsequently recombined or underwent further unimolecular decomposition. In conclusion, the kinetic model accurately simulated the density gradient profiles of all perfluoroalkanes.

Guzman, Eduardo H. [Brown Univ., Providence, RI (U↗
DOE OSTI2020

Diffusivity and Structure of Room Temperature Ionic Liquid in Various Organic Solvents

Room-temperature ionic liquids (RTILs) hold promise for applications in electric double layer capacitors (EDLCs), owing to a much wider potential window, lower vapor pressure, and better thermal and chemical stabilities compared to conventional aqueous and organic electrolytes. However, because the low diffusivity of ions in neat RTILs negates the EDLCs’ advantage of high power density, the ionic liquids are often used in mixture with organic solvents. In this study, we measured the diffusivity of cations and anions in RTIL, 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl) ([BMIM + ][TFSI – ]), mixed with 10 organic solvents, by using the pulsed-field gradient NMR method. The ion diffusivity was found to follow that of neat solvents and in most studied solvents showed an excellent agreement with the predicted values reported in the recent molecular dynamics (MD) study [Thompson, M. W.; J. Phys. Chem. B 2019, 123, 1340-1347]. In two solvents consisting of long-chain molecules, however, the MD simulations predictions slightly underestimated the ionic diffusivities. The degree of ion dissociation was also estimated for each solvent by comparing the ionic conductivity with the molar conductivity derived from the diffusion measurements. The degree of ion dissociation and the hydrodynamic radius of ions suggest that the ions are coordinated by ~1 solvent molecule. Finally, the scarcity of solvent–ion interactions explains the fact that the diffusivity of ions in the mixture significantly depends on the viscosity of the solvent.

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

Transition Path Sampling Study of the Feruloyl Esterase Mechanism

Serine hydrolases cleave peptide and ester bonds and are ubiquitous in nature, with applications in biotechnology, in materials, and as drug targets. The serine hydrolase two-step mechanism employs a serine–histidine–aspartate/glutamate catalytic triad, where the histidine residue acts as a base to activate poor nucleophiles (a serine residue or a water molecule) and as an acid to allow the dissociation of poor leaving groups. This mechanism has been the subject of debate regarding how histidine shuttles the proton from the nucleophile to the leaving group. To elucidate the reaction mechanism of serine hydrolases, we employ quantum mechanics/molecular mechanics-based transition path sampling to obtain the reaction coordinate using the Aspergillus niger feruloyl esterase A (AnFaeA) as a model enzyme. The optimal reaction coordinates include terms involving nucleophilic attack on the carbonyl carbon and proton transfer to, and dissociation of, the leaving group. During the reaction, the histidine residue undergoes a reorientation on the time scale of hundreds of femtoseconds that supports the “moving histidine” mechanism, thus calling into question the “ring flip” mechanism. We find a concerted mechanism, where the transition state coincides with the tetrahedral intermediate with the histidine residue pointed between the nucleophile and the leaving group. Moreover, motions of the catalytic aspartate toward the histidine occur concertedly with proton abstraction by the catalytic histidine and help stabilize the transition state, thus partially explaining how serine hydrolases enable poor nucleophiles to attack the substrate carbonyl carbon. Rate calculations indicate that the second step (deacylation) is rate-determining, with a calculated rate constant of 66 s –1 . Altogether, these results reveal the pivotal role of active-site dynamics in the catalytic mechanism of AnFaeA, which is likely similar in other serine hydrolases.

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

A Model for Atomic Precision p-Type Doping with Diborane on Si(100)-2×1

Diborane (B 2 H 6 ) is a promising molecular precursor for atomic precision p-type doping of silicon that has recently been experimentally demonstrated. In this study, we use density functional theory (DFT) calculations to determine the reaction pathway for diborane dissociating into a species that will incorporate as electrically active substitutional boron after adsorbing onto the Si(100)-2×1 surface. Our calculations indicate that diborane must overcome an energy barrier to adsorb, explaining the experimentally observed low sticking coefficient (<1 × 10 –4 at room temperature) and suggesting that heating can be used to increase the adsorption rate. Upon sticking, diborane has an ≈50% chance of splitting into two BH 3 fragments versus merely losing hydrogen to form a dimer such as B 2 H 4 . As boron dimers are likely electrically inactive, whether this latter reaction occurs is shown to be predictive of the incorporation rate. The dissociation process proceeds with significant energy barriers, necessitating the use of high temperatures for incorporation. Using the barriers calculated from DFT, we parameterize a Kinetic Monte Carlo model that predicts the incorporation statistics of boron as a function of the initial depassivation geometry, dose, and anneal temperature. Our results suggest that the dimer nature of diborane inherently limits its doping density as an acceptor precursor and furthermore that heating the boron dimers to split before exposure to silicon can lead to poor selectivity on hydrogen and halogen resists. This suggests that, while diborane works as an atomic precision acceptor precursor, other non-dimerized acceptor precursors may lead to higher incorporation rates at lower temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗