SEARCH · Search NASA
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.
Multi-Scale modeling of mercury (Hg) geochemical speciation and reactions has been performed by integrating atomistic quantum chemical calculations with continuum scale speciation models. Major progress has been made in the improvement of quantum chemical models to calculate critical thermodynamic data for Hg complexes in aquatic environments. Rapid and reliable quantum chemical approaches have been developed for calculating acid dissociation constants (pK a ) and stability constants (log K), with calculated mean unsigned errors of 0.5 and 1.5 log units, respectively for ligand molecules and Hg complexes. At the continuum scale, systematic analysis of uncertainty propagation in mercury (Hg) speciation modeling has been conducted and was used to identify environmental conditions under which thermodynamic constant uncertainties are significant and recommended to be accounted for. The integrated framework for multi-scale modeling of mercury geochemistry is open to the research community through the web-based multiscale modeling aqueous speciation resource, AQUA-MER. The improved quantum chemical approaches for thermodynamic constant calculations are accessible through AQUA-MER and can be used to provide the missing constants in the continuum scale speciation calculations. In addition to low molecular mass Hg complex speciation, modeling natural aquatic environments also involve the transport of high molecular weight dissolved organic matter (DOM) in reactive flows simultaneously with equilibrium and kinetic reactions. To this end, atomistic MD simulations were performed to capture the details of aggregation, mechanisms and distribution of functional groups in DOM at the molecular level. The elemental composition and calculated bulk properties of the DOM models are in close agreement with experimental measurements. A travel-time based reactive transport model in the hyporheic zone of stream corridors was established for the multicomponent Hg-DOM-S system and implemented through PFLOTRAN.
The Fourier transform microwave spectrum of 4-methylacetophenone recorded from 8 GHz to 18 GHz under jet-cooled conditions has revealed large tunneling splittings arising from a low barrier to internal rotation of the ring methyl group and small splittings from a high torsional barrier of the acetyl methyl group. The large splittings are especially challenging to model, while the small splittings are difficult to analyze due to the resolution limit of 120 kHz. The combination of two methyl groups undergoing internal rotations caused each rotational transition to split into five torsional species, which were resolved and fitted using a modified version of the XIAM code and the newly devel- oped ntop code to a root-mean-square deviation close to measurement accuracy, providing an estimate of the V 3 potential barriers of about 22 cm -1 and 584–588 cm -1 for the ring and the acetyl methyl groups, respectively. The assignment was aided by separately fitting the five torsional species using odd-power order operators. Only one conformer in which all heavy atoms are located on a symmetry plane could be identified in the spectrum, in agreement with results from conformation analysis using quantum chemical calculations.
Computer code was developed for calculating chemical quantities and qualities in equilibrium.
A pair of bimetallic compounds featuring Fe–Fe bonds, [Fe( i PrNPPh 2 ) 3 FeR] (R = PMe 3 , ≡N t Bu), have been investigated using High-Frequency Electron Paramagnetic Resonance (HFEPR) as well as field- and temperature-dependent 57 Fe nuclear γ resonance (Mössbauer) spectroscopy. To gain insight into the local site electronic structure, we have concurrently studied a compound containing a single Fe(II) in a geometry analogous to that of one of the dimer sites. Our spectroscopic studies have allowed for the assessment of the electronic structure via the determination of the zero-field splitting and 57 Fe hyperfine parameters for the entire series. We also report on our efforts to correlate structure with physical properties in metal–metal bonded systems using ligand field theory guided by quantum chemical calculations. Through the insight gained in this study, we discuss strategies for the design of single-molecule magnets based on polymetallic compounds linked via direct metal–metal bonds.
Conformational preferences of substituted, saturated cyclic ketones may deviate from those of their parent analogues. 1-Phenylpiperidin-4-one was found to exist in several forms differing by not only position (axial or equatorial) of the aryl-substituent but also by a configuration of the piperidine cycle itself: chair-Eq, chair-Ax and twist-Tw. Here, the ratios of the co-existing species contributions in gas phase predicted by various quantum chemical calculations were confirmed by synchronous gas-phase electron diffraction and mass spectrometric experiments carried out for the vapor at 337 K demonstrating a best fit with the B3LYP-D3/cc-pVTZ combination: Eq:Ax:Tw=55(13):22(9):23(10) vs. 40:35:25.
In the foreseeable future, an expedition may be undertaken to explore the planet Mars. Some of the power source options being considered for such a mission are photovoltaics, regenerative fuel cells and nuclear reactors. In addition to electrical power requirements, environmental conditions en route to Mars, in the planetary orbit and on the Martian surface must be simulated and studied in order to anticipate and solve potential problems. Space power systems components such as photovoltaic arrays, radiators, and solar concentrators may be vulnerable to degradation in the Martian environment. Natural characteristics of Mars which may pose a threat to surface power systems include high velocity winds, dust, ultraviolet radiation, large daily variation in temperature, reaction to components of the soil, atmosphere and atmospheric condensates as well as synergistic combinations. Most of the current knowledge of the characteristics of the Martian atmosphere and soil composition was obtained from the Viking 1 and 2 missions in 1976. A theoretical study is presented which was used to assess the effects of the Martian atmospheric conditions on the power systems components. A computer program written at NASA-Lewis for combustion research that uses a free energy minimization technique was used to calculate chemical equilibrium for assigned thermodynamic states of temperature and pressure. The power system component materials selected for this study include: silicon dioxide, silicon, carbon, copper, and titanium. Combinations of environments and materials considered include: (1) Mars atmosphere with power surface material, (2) Mars atmosphere and dust component with power surface material, and (3) Mars atmosphere and hydrogen peroxide or superoxide or superoxide with power system material. The chemical equilibrium calculations were performed at a composition ratio (oxidant to reactant) of 100. The temperature for the silicon dioxide material and silicon, which simulate photovoltaic cells, were 300 and 400 K; for carbon, copper and titanium, which simulate radiator surfaces, 300, 500, and 1000 K. All of the systems were evaluated at pressures of 700, 800, and 900 Pa, which stimulate the Martian atmosphere.
Electronically excited NdO is a possible product of the chemistry associated with the release of Nd into the ionosphere, and emission from these states may contribute to the observations following such experiments. To better characterize the energetics and spectroscopy of NdO, we report a combined experimental and theoretical study using slow photoelectron velocity-map imaging spectroscopy of cryogenically cooled NdO – anions (cryo-SEVI) supplemented by wave function-based quantum-chemical calculations. Using cryo-SEVI, we measure the electron affinity of NdO to be 1.0091(7) eV and resolve numerous transitions to low-lying electronic and vibrational states of NdO that are assigned with the aid of the electronic structure calculations. Additionally, temperature-dependent data suggest contributions from the (2)4.5 state of NdO – residing 2350 cm –1 above the ground anion state. Finally, photodetachment to higher-lying excited states of NdO is also reported, which may help to clarify observations from prior release experiments.
Quantum chemical calculations on atomistic systems have evolved into a standard approach to studying molecular matter. These calculations often involve a significant amount of manual input and expertise, although most of this effort could be automated, which would alleviate the need for expertise in software and hardware accessibility. Here, we present the AutoRXN workflow, an automated workflow for exploratory high-throughput electronic structure calculations of molecular systems, in which (i) density functional theory methods are exploited to deliver minimum and transition-state structures and corresponding energies and properties, (ii) coupled cluster calculations are then launched for optimized structures to provide more accurate energy and property estimates, and (iii) multi-reference diagnostics are evaluated to back check the coupled cluster results and subject them to automated multi-configurational calculations for potential multi-configurational cases. All calculations are carried out in a cloud environment and support massive computational campaigns. Key features of all components of the AutoRXN workflow are autonomy, stability, and minimum operator interference. We highlight the AutoRXN workflow with the example of an autonomous reaction mechanism exploration of the mode of action of a homogeneous catalyst for the asymmetric reduction of ketones.
Quantum chemical calculations of the C 6 H 5 O 2 potential energy surface (PES) were carried out to study the mechanism of the phenoxy + O( 3 P) and phenyl + O 2 reactions. CASPT2(15e,13o)/CBS//CASSCF(15e,13o)/DZP multi-reference calculations were utilized to map out the minimum energy path for the entrance channels of the phenoxy + O( 3 P) reaction. Stationary points on the C 6 H 5 O 2 PES were explored at the CCSD(T)-F12/cc-pVTZ-f12//B3LYP/6-311++G** level for the species with a single-reference character of the wave function and at the CASPT2(15e,13o)/CBS//B3LYP/6-311++G** level of theory for the species with a multi-reference character of the wave function. Conventional, variational, and variable reaction coordinate transition state theories were employed in Rice-Ramsperger-Kassel-Marcus Master Equation calculations to assess temperature- and pressure-dependent phenomenological rate constants and product branching ratios. Here, the main bimolecular product channels of the phenoxy + O( 3 P) reaction are concluded to be para/ortho-benzoquinone + H, 2,4-cyclopentadienone + HCO, and, at high temperatures, also phenyl + O 2 . The main bimolecular product channels of the phenyl + O 2 reaction include 2,4-cyclopentadienone + HCO at lower temperatures, and phenoxy + O( 3 P) at higher temperatures. Both for the phenoxy + O( 3 P) and phenyl + O 2 reactions, the collisional stabilization of peroxybenzene at low temperatures and high pressures competes with the bimolecular product channels.
Explore the source record for details and available documents.
Here, empirical, highly accurate non-relativistic electronic total atomization energies (eTAEs) are established by combining experimental or computationally converged treatments of the nuclear motion and relativistic contributions with the total atomization energies of HF, CO, N 2 , and H 2 O obtained from the Active Thermochemical Tables. These eTAEs, which have estimated (2σ) uncertainties of less than 10 cm -1 (0.12 kJ mol -1 ), form the basis for an analysis of high-level ab initio quantum chemical calculations that aim at reproducing these eTAEs for the title molecules. The results are then employed to analyze the performance of the high-accuracy extrapolated ab initio thermochemistry, or High-Accuracy Extrapolated Ab Initio Thermochemistry (HEAT), family of theoretical methods. The method known as HEAT-345(Q), in particular, is found to benefit from fortuitous error cancellation between its treatment of the zero-point energy, extrapolation errors in the Hartree-Fock and coupled cluster contributions, neglect of post-(T) core-correlation, and the basis-set error involved in higher-level correlation corrections. In addition to shedding light on a longstanding curiosity of the HEAT protocol—where the cheapest HEAT-345(Q) performs comparably to the theoretically more complete HEAT-456QP procedure—this study lays the foundation for extended HEAT variants that offer substantial improvements in accuracy relative to the established approaches.
This work presents a generalized approach for analytical method optimization that branches the gap between techniques historically employed and accurate modern optimization techniques suitable for various applications. The novelty of the described strategy is the utilization of multivariate, multiobjective optimization with Karush-Kuhn-Tucker conditions to bound the optimization space to solutions within the physical limitations of instrumentation. Briefly, the basic steps outlined in this paper are to (1) determine the objective(s) that should be maximized or minimized based on the goals of the analytical application, (2) conduct a screening experiment, (3) perform ANOVA to determine the parameters which have a statistically significant effect on the objective, (4) conduct an experiment (e.g., Box-Behnken design) to collect data for fitting the objective equation, and (5) determine the physical constraints of the parameters and solve the Lagrangian to determine the optimal method parameters. A broad approach to optimization target selection allows for robust method tuning to develop improved data sets amenable for chemometrics and machine learning algorithm development. Gas chromatography-mass spectrometry was selected as a use case due to its broad use across scientific fields and time-consuming method development involving numerous parameters. In conclusion, this strategy can reduce the cost of research, improve data quality, and enable the rapid development of new analytical technique.
The rate coefficient for the ion-molecule reaction NH3(+) + H2 yields NH4(+) + H has been calculated as a function of temperature with the use of the statistical phase space approach. The potential surface and reaction complex and transition state parameters used in the calculation have been taken from ab initio quantum chemical calculations. The calculated rate coefficient has been found to mimic the unusual temperature dependence measured in the laboratory, in which the rate coefficient decreases with decreasing temperature until 50-100 K and then increases at still lower temperatures. Quantitative agreement between experimental and theoretical rate coefficients is satisfactory given the uncertainties in the ab initio results and in the dynamics calculations. The rate coefficient for the unusual three-body process NH3(+) + H2 + He yields NH4(+) + H + He has also been calculated as a function of temperature and the result found to agree well with a previous laboratory determination.
High-pressure multiplexed photoionization mass spectrometry (MPIMS) with tunable vacuum ultraviolet (VUV) ionization radiation from the Lawrence Berkeley Labs Advanced Light Source is used to investigate the oxidation of diethyl ether (DEE). Kinetics and photoionization (PI) spectra are simultaneously measured for the species formed. Several stable products from DEE oxidation are identified and quantified using reference PI cross-sections. In addition, we directly detect and quantify three key chemical intermediates: peroxy (ROO˙), hydroperoxyalkyl peroxy (˙OOQOOH), and ketohydroperoxide (HOOP[double bond, length as m-dash]O, KHP). These intermediates undergo dissociative ionization (DI) into smaller fragments, making their identification by mass spectrometry challenging. Furthermore, with the aid of quantum chemical calculations, we identify the DI channels of these key chemical species and quantify their time-resolved concentrations from the overall carbon atom balance at T = 450 K and P = 7500 torr. This allows the determination of the absolute PI cross-sections of ROO˙, ˙OOQOOH, and KHP into each DI channel directly from experiment. The PI cross-sections in turn enable the quantification of ROO˙, ˙OOQOOH, and KHP from DEE oxidation over a range of experimental conditions that reveal the effects of pressure, O 2 concentration, and temperature on the competition among radical decomposition and second O 2 addition pathways.
It is known that the flow field of the redesigned solid rocket motor (RSRM) is very complicated due to the complex characteristics of turbulent multi-phase flow, chemical reaction, particle combustion, evaporation, breakup and agglomeration etc. It requires multi-phase calculations, chemical reaction simulation, and particle combustion, evaporation, and breakup models to obtain a better understanding of thermophysics for the RSRM design using numerical methods. Also, the slag buildup due to the molten particles is another factor affecting the performance of the RSRM. To achieve this goal, the volume of fluid (VOF) method is used to capture the free surface motion so as to simulate the accumulation of the molten particles (slag) of the RSRM. A finite rate chemistry model is used to simulate the chemical reaction effects. For multi-phase calculations, the Hermsen combustion model is used for the aluminum particle combustion analysis and the Taylor Analogy Breakup (TAB) model is used for the particle breakup analysis. An interphase mas-exchange model introduced by Spalding is used for the evaporation calculation. The particle trajectories are calculated using a one-step implicit method for several groups of particle sizes by which the drag forces and heat fluxes are then coupled with the gas phase equations. The preliminary results predicted a reasonable physical simulation of the particle effects using a simple two dimensional solid rocket motor configuration. It shows that the AL/AL2O3 particle sizes are reduced due to the combustion, evaporation, and breakup. The flow field is disturbed by the particles. Mach number distributions in the nozzle are deformed due to the effect of particle concentrations away from the center line.
State preparation for quantum algorithms is crucial for achieving high accuracy in quantum chemistry and competing with classical algorithms. The localized active space–unitary coupled cluster (LAS–UCC) algorithm iteratively loads a fragment-based multireference wave function onto a quantum computer. Here, in this study, we compare two state preparation methods, quantum phase estimation (QPE) and direct initialization (DI), for each fragment. We test the two state preparation methods on three systems, ranging from a model system, a set of interacting hydrogen molecules, to more realistic chemical problems, like the C–C double bond breaking in transbutadiene and the spin ladder in a bimetallic system. We analyze the impact of QPE parameters, such as the number of ancilla qubits and Trotter steps, on the prepared state. We find a trade-off between the methods, where DI requires fewer resources for smaller fragments, while QPE is more efficient for larger fragments. Our resource estimates highlight the benefits of system fragmentation in state preparation for subsequent quantum chemical calculations. These findings have broad applications for preparing multireference quantum chemical wave functions on quantum circuits that can be used for realistic chemical applications.