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.
Inverse-sandwich samarium and ytterbium biphenyl complexes were synthesized and characterized by X-ray crystallography. Combined experimental and computational studies indicated that they have distinct electronic structures and bonding interactions.
The effects of migration direction, barriers, and defect-interactions on the strain-modified oxygen-ion conductivity in rare-earth substituted ceria are deconvoluted in this combined experimental and computational study.
Here we report a comprehensive computational study of unsupervised machine learning for extraction of chemically relevant information in X-ray absorption near edge structure (XANES) and in valence-to-core X-ray emission spectra (VtC-XES) for classification of a broad ensemble of sulphorganic molecules. By progressively decreasing the constraining assumptions of the unsupervised machine learning algorithm, moving from principal component analysis (PCA) to a variational autoencoder (VAE) to t-distributed stochastic neighbour embedding (t-SNE), we find improved sensitivity to steadily more refined chemical information. Surprisingly, when embedding the ensemble of spectra in merely two dimensions, t-SNE distinguishes not just oxidation state and general sulphur bonding environment but also the aromaticity of the bonding radical group with 87% accuracy as well as identifying even finer details in electronic structure within aromatic or aliphatic sub-classes. We find that the chemical information in XANES and VtC-XES is very similar in character and content, although they unexpectedly have different sensitivity within a given molecular class. We also discuss likely benefits from further effort with unsupervised machine learning and from the interplay between supervised and unsupervised machine learning for X-ray spectroscopies. Our overall results, i.e., the ability to reliably classify without user bias and to discover unexpected chemical signatures for XANES and VtC-XES, likely generalize to other systems as well as to other one-dimensional chemical spectroscopies.
Carbon fixing materials are a new class of self-healing, self-reinforcing materials we have introduced that utilize ambient CO 2 to chemically add to an ever extending carbon backbone. This class of materials can utilize biological or non-biological photocatalysts and support a wide range of potential backbone chemistries. However, there is no analysis to date that describes their fundamental limits in terms of chemical kinetics and mass transfer. In this computational study, we employ a reaction engineering, and materials science analysis to answer basic questions about the maximum growth rate, photocatalytic requirements and limits of applicable materials. Our proposed mathematical framework envelops three main functions required for carbon fixing materials: (1) adsorption of CO 2 from air, (2) photocatalytic reduction of CO 2 into selected monomers, and (3) polymerization of CO 2 -derived monomers. First, by performing a Damköhler number analysis, we derive criteria for the cross over from kinetic control to mass transfer limited growth, setting upper limits on performance of a potential photocatalyst. Next, we analyze photocatalytic reduction of CO 2 to single carbon products, using known catalytic pathways and kinetic data. We identify formaldehyde as a C1 intermediate having unique potential for incorporation into the material backbone of carbon fixing materials. As an example, we find that a diamond-shaped reaction network graph for CO 2 reduction, passing through CO and HCOOH intermediates, accurately describes kinetic data for cobalt-promoted TiO 2 nanoparticles at room temperature and 1 atm CO 2 . Finally, as an applied case study, we consider and analyze the example of a carbon fixing poly(oxymethylene) system with embedded catalyst promoting the photocatalytic reduction of CO 2 to formaldehyde. The latter then trimerizes to a trioxane monomer which subsequently polymerizes to polyoxymethylene. This reaction engineering analysis introduces benchmarks for carbon fixing materials with respect to achievable rates of photocatalysis, adsorption, and polymerization. These results should prove valuable for the design, evaluation, and benchmarking of this emergent and new class of environmentally sustainable materials.
Remdesivir (RDV) prodrug can be metabolized into a triphosphate form nucleotide analogue (RDV-TP) to bind and insert into the active site of viral RNA dependent RNA polymerase (RdRp) to further interfere with viral genome replication. Here in this work, we computationally studied how RDV-TP binds and inserts to the SARS-CoV-2 RdRp active site, in comparison with natural nucleotide substrate adenosine triphosphate (ATP). To do that, we first constructed atomic structural models of an initial binding complex (active site open) and a substrate insertion complex (active site closed), based on high-resolution cryo-EM structures determined recently for SARS-CoV-2 RdRp or non-structural protein (nsp) 12, in complex with accessory protein factors nsp7 and nsp8. By conducting all-atom molecular dynamics simulation with umbrella sampling strategies on the nucleotide insertion between the open and closed state RdRp complexes, our studies show that RDV-TP can initially bind in a comparatively stabilized state to the viral RdRp active site, as it primarily forms base stacking with the template uracil nucleotide (nt +1), which under freely fluctuations supports a low free energy barrier of the RDV-TP insertion (~1.5 kcal mol -1 ). In comparison, the corresponding natural substrate ATP binds initially to the RdRp active site in Watson–Crick base pairing with the template nt, and inserts into the active site with a medium low free energy barrier (~2.6 kcal mol -1 ), when the fluctuations of the template nt are well quenched. The simulations also show that the initial base stacking of RDV-TP with the template can be specifically stabilized by motif C-S759, S682 (near motif B) with the base, and motif G-K500 with the template backbone. Although the RDV-TP insertion can be hindered by motif F-R555/R553 interaction with the triphosphate, the ATP insertion seems to be facilitated by such interactions. The inserted RDV-TP and ATP can be further distinguished by specific sugar interaction with motif B-T687 and motif A-D623, respectively.
Here a benchmark computational study of K-edge core-ionization energies of third-row elements using relativistic delta-coupled-cluster (ΔCC) methods and a revised core-valence separation (CVS) scheme is reported. High-level relativistic (HLR) corrections beyond the spin-free exact two-component theory in its one-electron variant (SFX2C-1e), including the contributions from two-electron picture-change effects, spin–orbit coupling, the Breit term, and quantum electrodynamics effects, have been taken into account and demonstrated to play an important role. Relativistic ΔCC calculations are shown to provide accurate results for core-ionization energies of third-row elements. The SFX2C-1e-CVS-ΔCC results augmented with HLR corrections show a maximum deviation of less than 0.5 eV with respect to experimental values.
Transition metal carbides are increasingly used as catalysts for the transformation of CO 2 into useful chemicals. Recently, the effect of nanostructuring of such carbides has started to gain relevance in tailoring their catalytic capabilities. Catalytic materials based on molybdenum carbide nanoparticles (MoC y ) have shown a remarkable ability to bind CO 2 at room temperature and to hydrogenate it into oxygenates or light alkanes. However, the involved chemistry is largely unknown. In the present work, a systematic computational study is presented aiming to elucidate the chemistry behind the bonding of CO 2 with a representative set of MoC y nanoparticles of increasing size, including stoichiometric and non-stoichiometric cases. Furthermore, the obtained results provide clear trends to tune the catalytic activity of these systems and to move towards more efficient CO 2 transformation processes.
The confinement of water molecules is vital in fields from biology to nanotechnology. The conditions allowing confinement in small finite polycyclic aromatic hydrocarbons (PAHs) are unclear, yet are crucial for understanding confinement in larger systems. Here, we report a computational study of water cluster confinement within PAHs dimers. Our results serve as a model for larger carbon allotropes and for understanding molecular interactions in confined systems. We identified size and structural motifs allowing confinement and demonstrated the motifs in various PAHs systems. Here we show that optimal OH∙∙∙π interactions between water clusters and the PAH dimer permit optimal confinement to occur. However, the lack of such interactions leads to the formation of CH∙∙∙O interactions, resulting in less ideal confinement. Confinement of layered clusters is also possible, provided that the optimal OH∙∙∙π interactions are conserved.
The controlled functionalization of a single fluorine in a CF 3 group is difficult and rare. Photochemical C–F bond functionalization of the sp 3 -C–H bond in trifluorotoluene, PhCF 3 , is achieved using catalysts made from earth-abundant lanthanides, (Cp Me4 ) 2 Ln(2-O-3,5- t Bu 2 -C 6 H 2 )(1-C{N(CH) 2 N( i Pr)}) (Ln = La, Ce, Nd and Sm, Cp Me4 = C 5 Me 4 H). The Ce complex is the most effective at mediating hydrodefluorination and defluoroalkylative coupling of PhCF 3 with alkenes; addition of magnesium dialkyls enables catalytic C–F bond cleavage and C–C bond formation by all the complexes. Mechanistic experiments confirm the essential role of the Lewis acidic metal and support an inner-sphere mechanism of C–F activation. Computational studies agree that coordination of the C–F substrate is essential for C–F bond cleavage. The unexpected catalytic activity for all members is made possible by the light-absorbing ability of the redox non-innocent ligands. The results described herein underscore the importance of metal–ligand cooperativity, specifically the synergy between the metal and ligand in both light absorption and redox reactivity, in organometallic photocatalysis.
A spectroscopic and computational study reveals that bulky electron-withdrawing imido substituents stabilize a rare S = 1 iron( iv )–bisimido complex with a pseudo doubly degenerate ground electronic structure, accounting for its facile HAT reactivity.
Crystallographic, spectroscopic, and computational studies on diiron( ii ) hydride complexes reveal rapid hydride motions that cause major changes in the electronic structures of the iron sites.
The stability of metal–organic frameworks (MOFs) in water affects their ability to function as chemical catalysts, their capacity as adsorbents for separations in water vapor presence, and their usefulness as recyclable water harvesters. Here, we have examined water stability of four node-modified variants of the mesoporous MOF, NU-1000, namely formate-, Acac-, TFacac-, and Facac-NU-1000, comparing these with node-accessible NU-1000. These NU-1000 variants present ligands grafted to NU-1000's hexa-Zr(IV)-oxy nodes by displacing terminal aqua and hydroxo ligands. Facac-NU-1000, containing the most hydrophobic ligands, showed the greatest water stability, being able to undergo at least 20 water adsorption/desorption cycles without loss of water uptake capacity. Computational studies revealed dual salutary functions of installed Facac ligands: (1) enhancement of framework mechanical stability due to electrostatic interactions; and (2) transformation and shielding of the otherwise highly hydrophilic nodes from H-bonding interactions with free water, presumably leading to weaker channel-stressing capillary forces during water evacuation – consistent with trends in free energies of dehydration across the NU-1000 variants. Water harvesting and hydrolysis of chemical warfare agent simulants were examined to gauge the functional consequences of modification and mechanical stabilization of NU-1000 by Facac ligands. The studies revealed a harvesting capacity of ~1.1 L of water vapor per gram of Facac-NU-1000 per sorption cycle. They also revealed retention of catalytic MOF activity following 20 water uptake and release cycles. This study provides insights into the basis for node-ligand-engendered stabilization of wide-channel MOFs against collapse during water removal.
We propose the general selection guidelines of amorphous cathode coatings for lithium-ion batteries based on an extensive high-throughput computational study and detailed ionic diffusion analysis.
Determining the factors that influence and can help predict energetic material sensitivity has long been a challenge in the explosives community. Decades of literature reports identify a multitude of factors both chemical and physical that influence explosive sensitivity; however no unifying theory has been observed. Recent work by our team has demonstrated that the kinetics of “trigger linkages” (i.e., the weakest bonds in the energetic material) showed strong correlations with experimental drop hammer impact sensitivity. These correlations suggest that the simple kinetics of the first bonds to break are good indicators for the reactivity observed in simple handling sensitivity tests. Herein we report the synthesis of derivatives of the explosive pentaerythritol tetranitrate (PETN) in which one, two or three of the nitrate ester functional groups are substituted with an inert group. Experimental and computational studies show that explosive sensitivity correlates well with Q (heat of explosion), due to the change in the number of trigger linkages removed from the starting material. In addition, this correlation appears more significant than other observed chemical or physical effects imparted on the material by different inert functional groups, such as heat of formation, heat of explosion, heat capacity, oxygen balance, and the crystal structure of the material.
This computational study shows that the interactions between nanoparticles and vesicles of diverse shapes are influenced not just by particle size and the extent of wrapping, but also by their initial positions and the sequence of interactions.
The biphenyl molecule (C 12 H 10 ) acts as a fundamental molecular backbone in the stereoselective synthesis of organic materials due to its inherent twist angle causing atropisomerism in substituted derivatives and in molecular mass growth processes in circumstellar environments and combustion systems. Here, we reveal an unconventional low-temperature phenylethynyl addition–cyclization–aromatization mechanism for the gas-phase preparation of biphenyl (C 12 H 10 ) along with ortho-, meta-, and para-substituted methylbiphenyl (C 13 H 12 ) derivatives through crossed molecular beams and computational studies providing compelling evidence on their formation via bimolecular gas-phase reactions of phenylethynyl radicals (C 6 H 5 CC, X 2 A 1 ) with 1,3-butadiene-d 6 (C 4 D 6 ), isoprene (CH 2 C(CH 3 )CHCH 2 ), and 1,3-pentadiene (CH 2 CHCHCHCH 3 ). The dynamics involve de-facto barrierless phenylethynyl radical additions via submerged barriers followed by facile cyclization and hydrogen shift prior to hydrogen atom emission and aromatization to racemic mixtures (ortho, meta) of biphenyls in overall exoergic reactions. Furthermore, these findings not only challenge our current perception of biphenyls as high temperature markers in combustion systems and astrophysical environments, but also identify biphenyls as fundamental building blocks of complex polycyclic aromatic hydrocarbons (PAHs) such as coronene (C 24 H 12 ) eventually leading to carbonaceous nanoparticles (soot, grains) in combustion systems and in deep space thus affording critical insight into the low-temperature hydrocarbon chemistry in our universe.
The functionalisation of organic linkers in metal–organic frameworks (MOFs) to improve gas uptake is well-documented. Although the positive role of free carboxylic acid sites in MOFs for binding gas molecules has been proposed in computational studies, relatively little experimental evidence has been reported in support of this. Primarily this is because of the inherent synthetic difficulty to prepare MOF materials bearing free, accessible –COOH moieties which would normally bind to metal ions within the framework structure. Here, we describe the direct binding of CO 2 and C 2 H 2 molecules to the free –COOH sites within the pores of MFM-303(Al). MFM-303(Al) exhibits highly selective adsorption of CO 2 and C 2 H 2 with a high selectivity for C 2 H 2 over C 2 H 4 . In situ synchrotron X-ray diffraction and inelastic neutron scattering, coupled with modelling, highlight the cooperative interactions of adsorbed CO 2 and C 2 H 2 molecules with free –COOH and –OH sites within MFM-303(Al), thus rationalising the observed high selectivity for gas separation.
We analyzed the hydroxyl librational signatures of five structurally related aluminum (oxy)hydroxides, using inelastic neutron scattering (INS) and plane-wave lattice dynamics simulations. A clear trend across these aluminum-containing phases illustrates the relationship between hydrogen bonding, local atomic structure, and the spectral location and profile of the librational bands. The INS spectra have been compared to previous optical spectroscopy and computational studies, highlighting the complementary nature of the INS technique. Taking into account other structurally or chemically related material analogs, we have identified a correlation between a blueshift (to higher energy) of the upper librational band edge and the geometry of the hydrogen bond interactions, mirroring (with opposite correlation) the well-known redshift in the intramolecular O–H stretching energy with increasing hydrogen bond strength. For hydroxyl groups that do not participate in hydrogen bonding effectively, the bending librations occur at lower energies and hybridize with metal–oxygen lattice modes. Standard density functional theory approximations, including dispersion corrections, struggle to correctly predict vibrational frequencies of motions dominated by H but perform well for metal–oxygen modes, allowing us to make detailed mode assignments in several cases, including a demonstration of how layer-to-layer disorder in boehmite hydrogen bond orientations is reflected in the sharp but minor low energy peaks (at ∼70–80 meV) of the INS spectrum.