Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “difficult”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Coupled Cluster Theory for Nonadiabatic Dynamics: Nuclear Gradients and Nonadiabatic Couplings in Similarity Constrained Coupled Cluster Theory

Coupled cluster theory is one of the most accurate electronic structure methods for predicting ground and excited state chemistry. However, the presence of numerical artifacts at electronic degeneracies, such as complex energies, has made it difficult to apply the method in nonadiabatic dynamics simulations. While it has already been shown that such numerical artifacts can be fully removed by using similarity constrained coupled cluster (SCC) theory [J. Phys. Chem. Lett. 2017, 8(19), 4801–4807], simulating dynamics requires efficient implementations of gradients and nonadiabatic couplings. Here, we present an implementation of nuclear gradients and nonadiabatic derivative couplings at the similarity constrained coupled cluster singles and doubles (SCCSD) level of theory, thereby making possible nonadiabatic dynamics simulations using a coupled cluster theory that provides a correct description of conical intersections between excited states. We present a few numerical examples that show good agreement with literature values and discuss some limitations of the method.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Revisiting Artifacts of Kohn–Sham Density Functionals for Biosimulation

We revisit the problem of unphysical charge density delocalization/fractionalization induced by the self-interaction error of common approximate Kohn–Sham (KS) density functional theory functionals on simulation of small to medium-sized proteins in a vacuum. Aside from producing unphysical electron densities and total energies, the vanishing of the HOMO–LUMO gap associated with the unphysical charge delocalization leads to an unphysical low-energy spectrum and catastrophic failure of most popular solvers for the KS self-consistent field (SCF) problem. We apply a robust quasi-Newton SCF solver to obtain solutions for some of these difficult cases. The anatomy of the charge delocalization is revealed by the natural deformation orbitals obtained from the density matrix difference between the Hartree–Fock and KS solutions; the charge delocalization not only can occur between charged fragments (such as in zwitterionic polypeptides) but also involves neutral fragments. The vanishing-gap phenomenon and troublesome SCF convergence are both attributed to the unphysical KS Fock operator eigenspectra of molecular fragments (e.g., amino acids or their side chains). Analysis of amino acid pairs suggests that the unphysical charge delocalization can be partially ameliorated by the use of some range-separated hybrid functionals but not by semilocal or standard hybrid functionals. Last, we demonstrate that solutions without the unphysical charge delocalization can be located even for semilocal KS functionals highly prone to such defects, but such solutions have non-Aufbau character and are unstable with respect to mixing of the non-overlapping “frontier” orbitals. Caution should be exercised when unexpectedly small (or vanishing) HOMO–LUMO gaps and atypical SCF convergence patterns (e.g., oscillatory) are observed in KS DFT simulations in any context (bio or otherwise).

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

Nonempirical Prediction of the Length-Dependent Ionization Potential in Molecular Chains

The ionization potential of molecular chains is well-known to be a tunable nanoscale property that exhibits clear quantum confinement effects. State-of-the-art methods can accurately predict the ionization potential in the small molecule limit and in the solid-state limit, but for intermediate, nanosized systems prediction of the evolution of the electronic structure between the two limits is more difficult. Recently, optimal tuning of range-separated hybrid functionals has emerged as a highly accurate method for predicting ionization potentials. This was first achieved for molecules using the ionization potential theorem (IPT) and more recently extended to solid-state systems, based on an ansatz that generalizes the IPT to the removal of charge from a localized Wannier function. Here, we study one-dimensional molecular chains of increasing size, from the monomer limit to the infinite polymer limit using this approach. By comparing our results with other localization-based methods and where available with experiment, we demonstrate that Wannier-localization-based optimal tuning is highly accurate in predicting ionization potentials for any chain length, including the nanoscale regime.

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

Data-Driven Kinetic Reaction Networks for Separation Chemistry

Understanding complex, multistep chemical reactions at the molecular level is a major challenge whose solution would greatly benefit the design and optimization of numerous chemical processes. The separation of rare-earth (4f) and actinide (5f) elements is an example where improving our chemical understanding is important for designing and optimizing new chemistries, even with a limited number of observations. Here, in this work, we leverage data-driven artificial intelligence and machine-learning approaches to develop kinetic reaction networks that describe the liquid–liquid extraction mechanism of uranium using N,N-di-2-ethylhexyl-isobutyramide (DEHiBA). Specifically, we compare and contrast the properties of two classes of models: (1) purely data-driven models that are regularized using chemistry-agnostic, L1 regression and (2) chemistry-informed models that are regularized using relative reaction energies provided by quantum mechanical calculations. We observe that purely data-driven models are unbiased, simple, and accurate in their predictions of experimental measurements when provided with sufficient data but are difficult to fully constrain and interpret. In contrast, chemistry-informed models exhibit significantly improved chemical interpretability and consistency, providing a detailed description of the separation process while achieving high accuracy through ensemble averaging. Overall, the dominant species predicted to be extracted into the organic phase is UO 2 (NO 3 ) 2 (DEHiBA) 2 , agreeing with experimental slope analysis, thermodynamic modeling, EXAFS, and crystal structures. This work demonstrates that leveraging the fundamental structure of the problem can lead to efficient learning schemes that provide both accurate predictions and chemical insights at a low computational cost.

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

Scaling Field-Theoretic Simulation for Multicomponent Mixtures with Neural Operators

Multicomponent polymer mixtures are ubiquitous in biological self-organization but are notoriously difficult to study computationally. Plagued by both slow single molecule relaxation times and slow equilibration within dense mixtures, molecular dynamics simulations are typically infeasible at the spatial scales required to study the stability of mesophase structure. Polymer field theories offer an attractive alternative, but analytical calculations are only tractable for mean-field theories and nearby perturbations, constraints that become especially problematic for fluctuation-induced effects such as coacervation. Here, we show that a recently developed technique for obtaining numerical solutions to partial differential equations based on operator learning, neural operators, lends itself to a highly scalable training strategy by parallelizing per-species operator maps. We illustrate the efficacy of our approach on six-component mixtures with randomly selected compositions and that it significantly outperforms the state-of-the-art pseudospectral integrators for field-theoretic simulations, especially as polymer lengths become long.

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

A Multiple Filter Based Neural Network Approach to the Extrapolation of Adsorption Energies on Metal Surfaces for Catalysis Applications

Computational catalyst discovery involves the development of microkinetic reactor models based on estimated parameters determined from density functional theory (DFT). For complex surface chemistries, the number of reaction intermediates can be very large, and the cost of calculating the adsorption energies by DFT for all surface intermediates even for one active site model can become prohibitive. In this paper, we have identified appropriate descriptors and machine learning models that can be used to predict a significant part of these adsorption energies given data on the rest of them. Furthermore, our investigations also included the case when the species data used to train the predictive model are of different size relative to the species the model tries to predict—this is an extrapolation in the data space which is typically difficult with regular machine learning models. Due to the relative size of the available data sets, we have attempted to extrapolate from the larger species to the smaller ones in the current work. Here, we have developed a neural network based predictive model that combines an established additive atomic contribution based model with the concepts of a convolutional neural network that, when extrapolating, achieves a statistically significant improvement over the previous models.

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

Robust, Accurate, and Efficient: Quantum Embedding Using the Huzinaga Level-Shift Projection Operator for Complex Systems

Using wave function (WF) in density functional theory (DFT) embedding methods provides a framework for performing localized, high-accuracy WF calculations on a system, while not incurring the full computational cost of the WF calculation on the full system. To effectively partition a system into localized WF and DFT subsystems, we utilize the Huzinaga level-shift projection operator within an absolutely localized basis. In this work, we study the ability of the absolutely localized Huzinaga level-shift projection operator method to study complex WF and DFT partitions, including partitions between multiple covalent bonds, a double bond, and transition-metal–ligand bonds. We find that our methodology can accurately describe all of these complex partitions. Additionally, we study the robustness of this method with respect to the WF method, specifically where the embedded systems were described using a multiconfigurational WF method. We found that the method is systematically improvable with respect to both the number of atoms in the WF region and the size of the basis set used, with energy errors less than 1 kcal/mol. Additionally, we calculated the adsorption energy of H2 to a model of an iron metal–organic framework (Fe-MOF-74) to within 1 kcal/mol compared to CASPT2 calculations performed on the full model while incurring only a small fraction of the full computational cost. In conclusion, this work demonstrates that the absolutely localized Huzinaga level-shift projection operator method is applicable to very complex systems with difficult electronic structures.

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

Synthesis of Perdeuterated Alkyl Amines/Amides with Pt/C as Catalyst under Mild Conditions

A convenient method for the synthesis of perdeuterated alkyl amides/amines is disclosed. Perdeuterated acetyl amides can be achieved by a hydrogen–deuterium (H/D) exchange protocol with Pt/C as a catalyst and D 2 O as a deuterium source under mild conditions. Further, after removal or reduction of the acetyl group, this protocol can provide perdeuterated primary, secondary, and tertiary amines, which are difficult to achieve via other methods.

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

Structure–Property Effects in the Generation of Transient Aqueous Benzoic Acid Anhydrides by Carbodiimide Fuels

The design of dissipative systems, which operate out-of-equilibrium by consuming chemical fuels, is difficult. As yet, there are a few examples of privileged fuel chemistries that can be broadly applied in abiotic systems in the same way that ATP hydrolysis is exploited throughout biochemistry. The key issue is that designing nonequilibrium systems is inherently about balancing the relative rates of coupled reactions. The use of carbodiimides as fuels to generate transient aqueous carboxylic anhydrides has recently been used in examples of new nonequilibrium materials and supramolecular assemblies. Here, we explore the kinetics of formation and decomposition of a series of benzoic anhydrides generated from the corresponding acids and EDC under typical conditions (EDC = N-(3-dimethylaminopropyl)-N'-ethylcarbodiimide hydrochloride). The reactions can be described by a simple mechanism that merges known behavior for the two processes independently. Structure–property effects in these systems are dominated by differences in the anhydride decomposition rate. The kinetic parameters allow trends in concentration-dependent properties to be simulated, such as reaction lifetimes, peak anhydride concentrations, and yields. For key properties, there are diminishing returns with the addition of increasing amounts of fuel. These results should provide useful guidelines for the design of functional systems making use of this chemistry.

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

High-Pressure Equation of State of 1,3,5-triamino-2,4,6-trinitrobenzene: Insights into the Monoclinic Phase Transition, Hydrogen Bonding, and Anharmonicity

The high-pressure equation of state (EOS) of energetic materials (EMs) is important for continuum and mesoscale models of detonation performance and initiation safety. Additionally, obtaining a high-fidelity EOS of the insensitive EM 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) has proven to be difficult because of challenges in experimental characterization at high pressures (HPs). In this work, powder X-ray diffraction patterns were fitted using the recently discovered monoclinic I2/a phase above 4 GPa, which shows that TATB is less compressible than when indexed with the triclinic $P\bar{1}$ phase. First-principles calculations were performed with Perdew–Burke–Ernzerhof (PBE) and PBE0 functionals including thermal effects using the $P\bar{1}$ phase. PBE0 improves the description of hydrogen bonding and thus predicts accurate planar a and b lattice parameters under ambient conditions. However, discrepancies in the predicted lattice parameters above 4–10 GPa compared with experimental measurements indexed with $P\bar{1}$ are further evidence of a structural modification at high pressure. Layer sliding defects are formed during molecular dynamics simulations, which induces an anharmonic effect on the thermal expansion of the c lattice parameter. In short, the results provide several insights into determining high-fidelity EOS parameters for TATB and other molecular crystals.

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

Inferring the Energetics of CO 2 –Aniline Adduct Formation from Vibrational Spectroscopy

Control of atmospheric CO 2 is an important contemporary scientific and engineering challenge. Towards this goal, the reaction of CO 2 with amines to form carbamate bonds is an established method for CO 2 capture. However, controllable reversal of this reaction remains difficult and requires tuning the energetics of the carbamate bond. Through IR spectroscopy, we show that a characteristic frequency observed upon carbamate formation varies as a function of the substituent’s Hammett parameter for a family of para- substituted anilines. We present computational evidence that the vibrational frequency of the adducted CO 2 serves as a predictor of the energy of formation of the carbamate. Electron donating groups typically enhance the driving force of carbamate formation by transferring more charge to the adducted CO 2 and thus increasing the occupancy of the anti-bonding orbital in the carbon-oxygen bonds. Increased occupancy of the anti-bonding orbital within adducted CO 2 indicates a weaker bond, leading to a red shift in the characteristic carbamate frequency. As a result, our work serves the large field of CO 2 capture research where spectroscopic observables, such as IR frequencies, are more easily obtainable and can stand in as a descriptor of driving forces.

Amines↗
DOE OSTI2023

Re-Examination of the N 2 O + O Reaction

The reaction of N 2 O with O is a key step in consumption of nitrous oxide in thermal processes. It has two product channels, NO + NO (R2) and N 2 + O 2 (R3). The rate constant for R2 has been measured both in the forward and the reverse direction at elevated temperature and is well established. However, the rate constant for the N 2 + O 2 channel (R3) has been difficult to quantify and has significant error limits. Here, the direct reaction on the triplet surface has a barrier of around 40 kcal mol –1 , and it is too slow for the N 2 + O 2 channel to have any practical significance. Recently, Pham and Lin (2022) suggested an alternative low activation energy reaction path that involves intersystem crossing and reaction on the singlet surface. In the present work, we re-examined a wide range of experiments relevant for the N 2 O + O reaction through kinetic modeling, paying attention to the impact of artifacts such as impurities and surface reactions. Experimental results from shock tubes and batch reactors on the final NO yield in N 2 O decomposition, covering temperatures of 973–2200 K and pressures of 0.013–11.5 atm, support k 3 ~ 0, consistent with the high activation energy for reaction on the triplet surface and a low probability of ISC.

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

Atomic and Electronic Structure in MgO–SiO 2

Understanding disordered structure is difficult due to insufficient information in experimental data. Here, we overcome this issue by using a combination of diffraction and simulation to investigate oxygen packing and network topology in glassy (g-) and liquid (l-) MgO–SiO 2 based on a comparison with the crystalline topology. We find that packing of oxygen atoms in Mg 2 SiO 4 is larger than that in MgSiO 3 , and that of the glasses is larger than that of the liquids. Moreover, topological analysis suggests that topological similarity between crystalline (c)- and g-(l-) Mg 2 SiO 4 is the signature of low glass-forming ability (GFA), and high GFA g-(l-) MgSiO 3 shows a unique glass topology, which is different from c-MgSiO 3 . We also find that the lowest unoccupied molecular orbital (LUMO) is a free electron-like state at a void site of magnesium atom arising from decreased oxygen coordination, which is far away from crystalline oxides in which LUMO is occupied by oxygen’s 3s orbital state in g- and l-MgO–SiO 2 , suggesting that electronic structure does not play an important role to determine GFA. We finally concluded the GFA of MgO–SiO 2 binary is dominated by the atomic structure in terms of network topology.

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

Comparison of Computational Strategies for the Calculation of the Electronic Coupling in Intermolecular Energy and Electron Transport Processes

Electronic couplings in intermolecular electron and energy transfer processes calculated by six different existing computational techniques are compared to nonorthogonal configuration interaction for fragments (NOCI-F) results. The paper addresses the calculation of the electronic coupling in diketopyrrolopyrol, tetracene, 5,5'-difluoroindigo, and benzene–Cl for hole and electron transport, as well as the local exciton and singlet fission coupling. NOCI-F provides a rigorous computational scheme to calculate these couplings, but its computational cost is rather elevated. The here-considered ab initio Frenkel–Davydov (AIFD), Dimer projection (DIPRO), transition dipole moment coupling, Michl–Smith, effective Hamiltonian, and Mulliken–Hush approaches are computationally less demanding, and the comparison with the NOCI-F results shows that the NOCI-F results in the couplings for hole and electron transport are rather accurately predicted by the more approximate schemes but that the NOCI-F exciton transfer and singlet fission couplings are more difficult to reproduce.

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

Realistic Ion Dynamics through Charge Renormalization in Nonaqueous Electrolytes

While many practically important electrolytes contain lithium ions, interactions of these ions are particularly difficult to probe experimentally because of their small X-ray and neutron scattering cross sections and large neutron absorption cross sections. Molecular dynamics (MD) is a powerful tool for understanding the properties of nonaqueous electrolyte solutions from the atomic level, but the accuracy of this computational method crucially depends on the physics built into the classical force field. In this study, we demonstrate that several force fields for lithium bistriflimide (LiTFSI) in acetonitrile yield a solution structure that is consistent with the neutron scattering experiments, yet these models produce dramatically different ion dynamics in solution. Such glaring discrepancies indicate that inadequate representation of long-range interactions leads to excessive ionic association and ion-pair clustering. We show that reasonable agreement with the experimental observations can be achieved by renormalization of the ion charges using a “titration” method suggested herewith. This simple modification produces realistic concentration dependencies for ionic diffusion and conductivity in <2 M solutions, without loss in quality for simulation of the structure.

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

Screening Ionic Liquids Based on Ionic Volume and Electrostatic Potential Analyses

Ionic liquids are known to have tunable solvation properties, based on the pairing of different anions and cations, but the compositional landscape is vast and challenging to navigate efficiently. Some computational screening protocols are available, but they can be either time-consuming or difficult to implement. Herein, we perform a detailed investigation of the fundamental role of electrostatic interactions in these systems. We effectively develop a bridge between the previous volume-based approach (VBA) with a quantum structure–property relationship (QSPR) approach to create fast, simple screening guidelines. We propose a new parameter that is applicable to both monovalent and multivalent ions, the ionic polarity index (IPI), which is defined as the ratio of the average electrostatic surface potential ($\bar{V}$) of the ion to the net charge of the ion (q). The IPI correlation has been tested on a diverse dataset of 121 ions and reliable predictions can be obtained within a homologous series of ionic liquid compounds.

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

First-Principles Simulations of CuCl in High-Temperature Water Vapor

Experimental data suggest that the solubility of copper in high-temperature water vapor is controlled by the formation of hydrated clusters of the form CuCl(H 2 O) n , where the average number of water molecules in the cluster generally increases with increasing density [Migdisov, A. A.; et al. Geochim. Cosmochim. Acta2014, 129, 33–53]. However, the precise nature of these clusters is difficult to probe experimentally. Moreover, there are some discrepancies between experimental estimates of average cluster size and prior simulation work [Mei, Y. Geofluids2018, 2018, 4279124]. Here, we have performed first-principles Monte Carlo (MC) and molecular dynamics (MD) simulations to explore these clusters in finer detail. We find that molecular dynamics is not the most appropriate technique for studying aggregation in vapor phases, even at relatively high temperatures. Specifically, our MD simulations exhibit substantial problems in adequately sampling the equilibrium cluster size distribution. In contrast, MC simulations with specialized cluster moves are able to accurately sample the phase space of hydrogen-bonding vapors. At all densities, we find a stable, slightly distorted linear H 2 O–Cu–Cl structure, which is in agreement with the earlier simulations, surrounded by a variable number of water molecules. The surrounding water molecules do not form a well-defined second solvation shell but rather a loose network of hydrogen-bonded water with molecular CuCl on the outside edge of the water cluster. We also find a broad distribution of hydration numbers, especially at higher densities. In contrast to previous simulation work but in agreement with experimental data, we find that the average hydration number substantially increases with increasing density. Moreover, the value of the hydration number depends on the choice of cluster definition.

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

Ion Pairing and Molecular Orientation at Liquid/Liquid Interfaces: Self-Assembly and Function

We report that molecular orientation plays a pivotal role in defining the functionality and chemistry of interfaces, yet accurate measurements probing this important feature are few, due, in part, to technical and analytical limitations in extracting information from molecular monolayers. For example, buried liquid/liquid interfaces, where a complex and poorly understood balance of inter- and intramolecular interactions impart structural constraints that facilitate the formation of supramolecular assemblies capable of new functions, are difficult to probe experimentally. Here, we use vibrational sum-frequency generation spectroscopy, numerical polarization analysis, and atomistic molecular dynamics simulations to probe molecular orientations at buried oil/aqueous interfaces decorated with amphiphilic oligomers. We show that the orientation of self-assembled oligomers changes upon the addition of salts in the aqueous phase. The evolution of these structures can be described by competitive ion effects in the aqueous phase altering the orientations of the tails extending into the oil phase. These specific anionic effects occur via interfacial ion pairing and associated changes in interfacial solvation and hydrogen-bonding networks. These findings provide more quantitative insight into orientational changes encountered during self-assembly and pave the way for the design of functional interfaces for chemical separations, neuromorphic computing applications, and related biomimetic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗