Cascade of Spin-State Transitions in the Intermetallic Marcasite FeP 2
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Methods for calibration of non-contact velocity measurements and systems for implementing the same are described. Generally, the method comprises inducing a shock wave into a sample at a stress intensity that varies across the sample's elastic limit, which corresponds to the elastic-plastic state transition of the sample. That transition state may be at the sample's Hugoniot elastic limit. The velocity of the sample is measured using a non-contact velocity measurement instrument such as a velocimeter. The measurement may be compared to a predicted velocity or a velocity measurement made by another system to determine calibration parameters.
Transit buses operate primarily in dense urban areas, where nearby populations face increased exposure to fine particulates, nitrogen oxides, and other harmful pollutants. Electrifying transit buses presents a clear opportunity to reduce greenhouse gas emissions and improve urban air quality. However, widespread adoption may pose significant energy and infrastructure challenges, which can be mitigated through proactive planning and investment. This report presents a robust modeling framework and an initial estimation of the hourly electricity demand at transit bus depots across the United States. The resulting depot-level dataset, available at data.nrel.gov/submissions/282, provides valuable insights for infrastructure planning and electricity demand forecasting, supporting the scalable electrification of transit bus fleets nationwide.
Monitoring nuclear reactors is an important safety and security task with growing requirements. We explore the possibility of using seismic and acoustic data for inferring the power level of an operating reactor. Continuous data recorded at a single seismoacoustic station that is located about 50 m away from a research reactor was visualized and analyzed. The data show a clear correlation between seismoacoustic features and reactor main operational states. We designed a workflow that includes two machine learning (ML) models to classify the reactor operational states (OFF, transition, and ON) and estimate reactor power levels (10%, 30%, 50%, 70%, and 90%). We applied and compared five ML algorithms for the reactor OFF-transition-ON and four approaches for the power level classification. We also compared the performance of ML models trained with seismic-only, acoustic-only, and both types of data. Five-fold cross validations were implemented to assure a thorough evaluation of the model performances. Additionally, the results show the extreme boosting gradient algorithm worked best for the first model, whereas random forests performed best for the second model. Combining seismic and acoustic data leads to better performance than using a single type of data. Seismic data contributed more than acoustic data for both models. We reached an accuracy of 0.98 for reactor OFF and ON. The accuracies for the transition state and power levels are less optimal with a minimum accuracy of 0.66. However, our results suggest seismic and acoustic data contain useful information about the transition state as well as power levels. Seismic and acoustic data could be integrated with other observations to improve monitoring performance.
A thermally sensitive ionic redox transistor comprises a channel, a reservoir layer, and an electrolyte layer disposed between the channel and the reservoir layer. A conductance of the channel is varied by changing concentration of ions in the channel layer. The electrolyte layer is configured to undergo a state change at a state transition temperature. Below the state transition temperature, ions in the electrolyte layer are substantially immobile. Above the state transition temperature, ions can move freely between the reservoir layer and the channel across the electrolyte layer in response to a voltage being applied between the channel and the reservoir layer. When the device is cooled below the state transition temperature or temperature range, the ions are trapped in one or more of the layers because the electrolyte layer loses its ionic conductivity. A state of the redox transistor can be read by measuring the conductance of the channel.
The net rate of a composite reaction is the difference between the forward and reverse reaction rates, which are kinetically distinct despite sharing elementary reaction steps and therefore have different rate-controlling transition states and species. Thus, degrees of rate control defined to identify rate-controlling transition states and species for the net rate confound contributions from the forward and reverse reactions. Herein, the forward and reverse degrees of rate control are defined to independently quantify the extent that species and transition states control the forward and reverse rates in reversible reactions. These degrees of rate control are defined as the relative change in the forward and reverse reaction rates per k B T decrease in the standard-state molecular free energies of transition states and species, and they are related to experimentally measurable quantities such as steady-state fractional coverages, reaction orders, and energies and entropies of activation of the forward and reverse reaction rates. The forward/reverse degrees of rate control represent stoichiometric coefficients for species and transition states in an equilibrium between the apparent initial states and the apparent transition states of the apparent rate-controlling steps of forward and reverse overall reactions. At equilibrium, the apparent transition states for the forward reaction and reverse reaction converge, and thus the forward and reverse rate-controlling steps combine to form a single apparent rate-controlling step at equilibrium. Here, this apparent rate-controlling step is comprised of an apparent initial state, transition state, and final state, where the apparent final state of the forward reaction is the apparent initial state of the reverse reaction. The apparent rate-controlling step behaves identically to an elementary step reaction at equilibrium with a pseudo-mass-action rate function given by the transition-state-theory (TST) form rate function (Foley and Bhan, 2020) with a stoichiometric number equal to the affinity-averaged stoichiometric number, $\bar{σ}$ The ratio of the forward and reverse TST-form rate functions is identical to the overall thermodynamic equilibrium relation, consistent with the principles of microscopic reversibility and detailed balance at equilibrium.
Systems, methods, and computer media for mitigating cybersecurity vulnerabilities of systems are provided herein. A current cybersecurity maturity of a system can be determined based on maturity criteria. The maturity criteria can be ranked based on importance. Solution candidates for increasing the cybersecurity maturity of the system can be determined based on the ranking. The solution candidates specify cybersecurity levels for the maturity criteria. A present state value reflecting the current cybersecurity maturity of the system can be calculated. For the solution candidates, an implementation state value and a transition state value can be determined. The implementation state value represents implementation of the maturity levels of the solution candidate, and the transition state value represents a transition from the present state value to the implementation state value. Based on the transition state values, a solution candidate can be selected for the system, and the system can be modified accordingly.
“Degree of rate control” (DRC) analysis provides a quantitative approach for analysing the kinetics of multi-step reaction mechanisms that has been widely applied to both heterogeneous and homogeneous catalysis research, as well as electrocatalysis. The DRC of any given transition state or intermediate is defined as a partial derivative such that it approximately equals the fractional increase in net rate to the product of interest per differential decrease in its standard-state free energy for that species (÷RT), holding constant the standard-state free energies of all other transition states and intermediates. Even very complex mechanisms usually have only a few species with non-zero DRCs and are thus the species whose interactions with the catalyst most strongly affect the net rate. These key DRC values thus offer a simple and intuitive route to optimize catalyst materials, especially with the assistance of computational methods like density functional theory (DFT). These high-DRC species are also the species whose energetics must be most accurately measured or calculated to achieve an accurate kinetic model for any reaction mechanism. In simple cases with a single “rate-determining step”, the DRC for its transition state (TS) is + 1. Catalyst-bound intermediates, on the other hand, often have negative DRCs equal to a small integer times their fractional occupancy of catalyst sites. The apparent activation energy equals a weighted average of the standard-state enthalpies (relative to reactants) of all the species (intermediates, transition states and products) in the reaction mechanism, each weighted by its DRC (+RT). It has been shown that the apparent transfer coefficient in electrocatalysis, an inverted form of the Tafel slope, is a weighted average of the number of electrons transferred to generate each intermediate or product species in the mechanism, weighted again by the DRC. Quantitative analysis of kinetic isotope effects (KIEs, or the ratio of net rates for different reactant isotopes) in complex mechanisms has shown that the logarithm of the KIE equals the weighted average over all species in the mechanism of the difference between the two isotopes in their standard-state free energies (÷RT), again weighted by the DRC. The reaction orders with respect to fluid-phase concentrations of reactants, products and intermediates have also been proven to be directly related to DRCs. Thus, there are numerous experimental observables which equate to short linear combinations of DRCs, so that combinations of experimental measurements might provide access to DRC values. Since its invention for transition states in 1994 and its generalization to include intermediates in 2009, the DRC has thus far mainly been calculated for microkinetic models based either entirely on DFT or on DFT with the key energies (i.e., those for high-DRC species) being fine-tuned to match experiments. The relationships summarized in this work provide new opportunities for using experiments earlier in the development and optimization of microkinetic models that require input from computational catalysis.
The mechanistic roles of free hydronium ions, surface hydrides, and interfacial protons during guaiacol hydrodeoxygenation (HDO) on ruthenium nanoparticles are established. As guaiacol adsorbs on Ru, it loses its strong aromaticity and undergoes a rapid H-shift from its hydroxyl to meta carbons (in relation to its hydroxyl group), leading enol and keto surface isomers to exist in chemical equilibrium. HDO occurs via a hydridic H-adatom (H*) attack to the enol, followed by a kinetically relevant C-O bond rupture step, during which water shuttles the hydroxyl proton, enabling its intramolecular attack to the methoxy, evolving a high charged [Ru(s)-(C6H5O-)…(H+)…OCH3]† transition state. The competing HYD begins with a rapid H* attack to the keto, before a second, kinetically relevant H* attack, without proton involvement. Water, despite shifting the thermodynamics towards the more polar surface keto, promotes HDO to a much greater extent than HYD, because of its dual catalytic roles—it mobilizes hydroxyl proton (Brønsted acid) to cleave the strong C-O bond, synchronizing with the Ru metal surface (base) function that stabilizes the resulting [Ru(s)-(C6H5O-)…(H+)…OCH3]† transition state, and the water layers solvate this charged transition state, further reducing its free energy. Free hydronium ions do catalyze a separate homogeneous enol-keto isomerization, but this reaction is kinetically unrelated to HDO catalysis. This mechanistic picture explains the strong effects of polar protic solvent in hydrodeoxygenation, highlighting (i) the requirements of surface hydrides and interfacial protons acting in tandem to complete a HDO turnover and (ii) the cooperative role of protic solvent and metal surface in breaking the aromaticity and stabilizing charged reactive precursors and transition states.
The mechanism for hydrocarbon ring growth in sooting environments is still the subject of considerable debate. The reaction of phenyl radical (C 6 H 5 ) with propargyl radical (H 2 CCCH) provides an important prototype for radical–radical ring-growth pathways. We studied this reaction experimentally over the temperature range of 300–1000 K and pressure range of 4–10 Torr using time-resolved multiplexed photoionization mass spectrometry. We detect both the C 9 H 8 and C 9 H 7 + H product channels and report experimental isomer-resolved product branching fractions for the C 9 H 8 product. We compare these experiments to theoretical kinetics predictions from a recently published study augmented by new calculations. Here, these ab initio transition state theory-based master equation calculations employ high-quality potential energy surfaces, conventional transition state theory for the tight transition states, and direct CASPT2-based variable reaction coordinate transition state theory (VRC-TST) for the barrierless channels. At 300 K only the direct adducts from radical–radical addition are observed, with good agreement between experimental and theoretical branching fractions, supporting the VRC-TST calculations of the barrierless entrance channel. As the temperature is increased to 1000 K we observe two additional isomers, including indene, a two-ring polycyclic aromatic hydrocarbon, and a small amount of bimolecular products C 9 H 7 + H. Our calculated branching fractions for the phenyl + propargyl reaction predict significantly less indene than observed experimentally. We present further calculations and experimental evidence that the most likely cause of this discrepancy is the contribution of H atom reactions, both H + indenyl (C 9 H 7 ) recombination to indene and H-assisted isomerization that converts less stable C 9 H 8 isomers into indene. Especially at low pressures typical of laboratory investigations, H-atom-assisted isomerization needs to be considered. Regardless, the experimental observation of indene demonstrates that the title reaction leads, either directly or indirectly, to the formation of the second ring in polycyclic aromatic hydrocarbons.
Pilgrim is a program written in Python and designed to use direct dynamics in the calculation of thermal rate constants of chemical reactions by the variational transition state theory (VTST), based on electronic structure calculations for the potential energy surface. Pilgrim can also simulate reaction mechanisms using kinetic Monte Carlo (KMC). For reaction processes with many elementary steps, the rate constant of each of these steps can be calculated by means of conventional transition state theory (TST) or by using VTST. In the current version, Pilgrim can evaluate thermal rates using the canonical version of reaction-path VTST, which requires the calculation of the minimum energy path (MEP) associated with each elementary step or transition structure. Multi-dimensional quantum effects can be incorporated through the small-curvature tunneling (SCT) approximation. These methodologies are available both for reactions involving a single structure of the reactants and the transition state and also for reactions involving flexible molecules with multiple conformations of the reactant and/or of the transition state. For systems with many conformers, the program can evaluate each of the elementary reaction rate constants by multipath canonical VTST or multi-structural VTST. Moreover, the reactant can be unimolecular or bimolecular. Torsional anharmonicity can be incorporated through either the MSTor or the Q2DTor programs. Dual-level calculations are also available in Pilgrim: automatic high-level single-point energies can be used to correct the energy of reactants, transition states, products, and MEP points using the interpolated single-point energies (ISPE) algorithm. When the rate constants of all the chemical processes of interest are known, by means of their calculation using Pilgrim or alternatively through analytical fits to the rate constants as functions of temperature, it is possible to simulate a multistep mechanism under specified laboratory conditions using KMC. Finally, this algorithm allows performing a kinetic simulation to monitor the evolution of each chemical species with time and obtain the product yields.
Brønsted acids mediate chemical transformations via proton transfer to bound species and interactions between the conjugate anion and bound cationic intermediates and transition states that are also stabilized by van der Waals forces within voids of molecular dimensions in inorganic hosts. This Feature Article describes the relevant descriptors of reactivity in terms of the properties of acids and molecules that determine their ability to donate and accept protons and to reorganize their respective charges to optimize their interactions at bound states. The deprotonation energy (DPE) of the acids and the protonation energy (Eprot) of the gaseous analogs of bound intermediates and transition states reflect their respective properties as species present at non-interacting distances. These properties accurately describe the reactivity of acids of a given type, such as polyoxometalates (POM) with a given type of addenda atom but different central atoms and heterosilicates, for different families of reactions. They do not fully capture, however, differences among acid types (e.g., Mo and W POM, heterosilicates, mineral acids) for diverse types of chemical transformations (e.g., elimination, isomerization, dimerization, condensation). The incompleteness of such descriptors reflects their inability to describe how protonated molecular species and conjugate anions restructure their respective charges when present as a binding pair at interacting distances. Such interaction energies represent electrostatic forces that depend on charge distributions in the cations and anions and the ability to reorganize the distributions to maximize the interactions. In the case of deprotonation, the electrostatic and charge reorganization components of DPE for various acids solely reflect the ability of the conjugate anion to accept and distribute the negative charge, a characteristic unique of each type of solid acid and specifically of the composition of its extended conjugate anion framework. The energy required to accept and rearrange the positive charge in bound intermediates and transition states reflects, in turn, their respective ability to recover the ionic and covalent components of DPE, the energy required to detach proton from conjugate anions. The DPE components and the recovery fractions together lead to a modified DPE, which captures only the part of DPE that remains unrecovered by the ion-pair interactions at bound intermediates and transition states, as the unifying descriptor for broad families of acids and reactions. The electrostatic and charge reorganization energies involved in these general descriptors are placed in historical context by assessing their connections to the heuristics of hard–soft acid–base displacements. Further development of these concepts requires benchmarking and extension of electrostatic and reorganization components of energies for a more diverse set of reaction types and acid families and advancement of methods for more efficient calculations of electrostatic interactions. Reactivity descriptors must also account for dispersive interactions between host cavities and guest molecules, requiring a framework analogous to the one described here for ion-pair interactions; these dispersive interactions depend on the fit between their shapes and sizes as well as their ‘‘structural stiffness’’ that determines the ability to modify the shapes of molecules and voids to minimize free energy. Entropy considerations and estimates of their dependence on properties of catalysts and molecules are also required for accurately determining Gibbs free energies that ultimately determine reaction rates.
The rate of hydride transfer from three Ir hydride complexes of the type Cp*Ir( R bpy)H + (Cp* = C 5 Me 5 ; R bpy = 4,4′-R-2,2′-bipyridine, R = OMe, H, CO 2 Me) to six N-methylacridinium ( R Acr + ) acceptors with electronically different substituents in the 2- or 2,7-positions were measured. Using the thermodynamic hydricity of the donors and the hydride affinity of the acceptors the thermodynamic driving forces for hydride transfer were determined. Brønsted plots, which correlate kinetic and thermodynamic hydricity, demonstrate distinct linear free energy relationships for each complex, with different Brønsted α values. Thus, at the same driving force hydride transfer from Cp*Ir( OMe bpy)H + is faster than for Cp*Ir(bpy)H + or Cp*Ir( CO2Me bpy)H + . Experimental and computational analyses are consistent with a concerted hydride transfer mechanism for all Ir complexes. As the thermodynamic driving force increases an earlier transition state is observed and all transition states also include π-stacking interactions between the donor and acceptor, which likely contribute to the different α values. The experimental data fits well to the Marcus model, enabling the determination of reorganization energies (λ) that range from 58 to 69 kcal mol -1 . These are lower than λ values for hydride transfer reactions involving organic donors and acceptors. This work provides a rare example of the correlation of kinetic and thermodynamic hydricity using only experimental data and shows that hydride transfer reactions involving metal hydrides can follow Marcus theory. Furthermore, the findings offer insight into controlling metal-catalyzed hydride transfer reactions, which is valuable for designing improved systems for a range of transformations.
Experimentally, the thermal gas-phase deazetization of 2,3-diazabicyclo[2.2.1]hept-2-ene (1) results in the loss of N2 and the formation of bicyclo products 3 (exo) and 4 (endo) in a nonstatistical ratio, with preference for the exo product. Here, we report unrestricted M06-2X quasiclassical trajectories initialized from the concerted N2 ejection transition state that were able to replicate the experimental preference to form 3. We found that the 3:4 ratio results from the relative amounts of very fast (ballistic) exotype trajectories versus trajectories that lead to the 1,3-diradical intermediate 2. These quasiclassical trajectories provided a set of transition-state vibrational, velocity, momenta, and geometric features for the machine learning analysis. Additionally, a selection of popular supervised classification algorithms (e.g., random forest) provided poor prediction of trajectory outcomes based on only transition-state vibrational quanta and energy features. However, these machine learning models provided more accurate predictions using atomic velocities and atomic positions, attaining ~70% accuracy using initial conditions and between 85 and 95% accuracy at later reaction time steps. This increased accuracy allowed the feature importance analysis to reveal that, at the later-time analysis, the methylene bridge out-of-plane bending is correlated with trajectory outcomes for the formation of either the exo product or toward the diradical intermediate. Possible reasons for the struggle of machine learning algorithms to classify trajectories based on transition-state features is the heavily overlapping feature values, the finite but very large possible vibrational mode combinations, and the possibility of chaos as trajectories propagate. We examined this chaos by comparing a set of nearly identical trajectories that differed by only a very small scaling of the kinetic energies resulting from the transition-state reaction coordinate.
Knowledge of the degradation kinetics for polymer materials is important for understanding thermal stability. In this study, evolved gas analysis–mass spectrometry and pyrolysis gas-chromatography-mass spectrometry were evaluated for the potential to deliver additional insight into thermal degradation kinetics of diglycidal ether of bisphenol F (DGEBF) epoxy thermoset under inert and oxidative atmospheres. Degradation products of selected precursor ions were evaluated for their uniqueness to the specific precursor using extracted ion thermographs. Unique mass peaks, solely attributed to a single reaction pathway of a specific product, were determined from extracted ion thermographs and used to determine both activation energy (E a ) and pre-exponential factors for the specific primary reaction pathways. These primary reaction pathways for DGEBF epoxy degradation were then evaluated in the context of transition state theory (TST) and related transition state enthalpies (ΔH ‡ ) and entropies (ΔS ‡ ) of activation to further elucidate the degradation process. It was determined under pyrolysis conditions, as suggested by the E a , the formation of bisphenol F monomer was the rate-limiting step toward the formation of xanthene and phenol. In contrast, under thermo-oxidative conditions, reactions involving oxygen containing species were identified as the rate-limiting step for all observed products based on the large negative ΔS ‡ calculated from TST. This work demonstrates a powerful combination of technique and theory that can provide new insight into the degradation of polymer materials.
Predicting and controlling nanostructure formation during nucleation can pave the way to synthesizing novel energy materials via crystallization. However, such control over nucleation and crystallization remains challenging due to an inadequate understanding of critical factors that govern evolving atomistic structures and dynamics. Here, in this paper, we utilize coordination number as a reaction coordinate and rate theory to investigate how sodium sulfate, commonly known as a phase-change energy material, nucleates in a supersaturated aqueous solution. In conjunction with ab initio and force field-based molecular dynamics simulation, the rate theoretical analysis reveals that sodium sulfate from an initially dissolved metastable state transits to a heterogeneous mixture of prenucleated clusters and finally to a large cylindrical zigzag morphology. Measurements of Raman spectra and their ab initio modeling confirm that this nucleated morphology contains a few waters for every sulfate. Rate processes such as solvent exchange and desolvation exhibit high sensitivity to the evolving prenucleation/nucleation structures, providing a means to distinguish between critical nucleation precursors. Desolvation and forming the first-shell interionic coordination structure via monomer-by-monomer addition around sulfates are found to explain the formation of large nuclei. Thus, a detailed understanding of the step-by-step structure formation across scales has been achieved. This can be leveraged to predict nucleation-related structures and dynamics and potentially control the synthesis of novel phase-change materials for energy applications.
Alloying is well-known to improve the dehydrogenation selectivity of pure metals, but there remains considerable debate about the structural and electronic features of alloy surfaces that give rise to this behavior. To provide molecular-level insights into these effects, a series of Pd intermetallic alloy catalysts with Zn, Ga, In, Fe and Mn promoter elements was synthesized, and the structures were determined using in situ X-ray absorption spectroscopy (XAS) and synchrotron X-ray diffraction (XRD). The alloys all showed propane dehydrogenation turnover rates 5–8 times higher than monometallic Pd and selectivity to propylene of over 90%. Moreover, among the synthesized alloys, Pd 3 M alloy structures were less olefin selective than PdM alloys which were, in turn, almost 100% selective to propylene. This selectivity improvement was interpreted by changes in the DFT-calculated binding energies and activation energies for C–C and C–H bond activation, which are ultimately influenced by perturbation of the most stable adsorption site and changes to the d-band density of states. Furthermore, transition state analysis showed that the C–C bond breaking reactions require 4-fold ensemble sites, which are suggested to be required for non-selective, alkane hydrogenolysis reactions. These sites, which are not present on alloys with PdM structures, could be formed in the Pd 3 M alloy through substitution of one M atom with Pd, and this effect is suggested to be partially responsible for their slightly lower selectivity.