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.
Density functional theory and microkinetic modeling show that an Au@NiO core–shell nanostructure enhances the oxidative esterification of aldehydes by electron transfer from NiO to Au, stabilizing intermediates and lowering activation energies.
Rigidified heterocycles encompassing two dibenzo[b,f]oxazepine scaffolds and a 1,4-dihydropyrrolo[3,2-b]pyrrole core show intense blue emission and low turn-on voltage (≈3.5 V) in OLED devices.
This work investigates the high-pressure behavior of α-furil. At ∼1.6 GPa α-furil transforms into the previously theorizedtrans-planar conformation before undergoing a Diels–Alder reaction to form chemically homogeneous “double-core” nanothreads.
All-solid-state lithium batteries (ASSLBs) are emerging as promising next-generation batteries for electric vehicles owing to their high energy densities and safety features. However, challenges such as inadequate material percolation and low cathode utilization often hinder their potential. This paper presents a core–shell approach to optimize the cathode active material (CAM) utilization. The resultant CAM composite showed high ionic conductivity, a highly dense microstructure with <10% porosity, and minimal stack pressure changes during electrochemical cycling. The maximum CAM utilization was achieved while effectively mitigating electrochemo-mechanical side reactions by applying a uniformly coated Li 6 PS 5 Cl solid electrolyte layer (≈500 nm) and a LiNbO 3 buffer layer (≈10 nm) onto LiNi 0.8 Mn 0.1 Co 0.1 O 2 particles (LPSCl@LNO@NMC). The engineered LPSCl@LNO@NMC composites, which incorporated a 5 wt% LPSCl coating on LNO@NMC powders, exhibited a dense microstructure that enhanced the mechanical stability at the cathode. Sulfide-based solid electrolyte (SSE)/SSE contact provided better ionic pathways within the composite and increased CAM utilization. Thus, an enhanced reversible capacity (197 mA h g -1 ) and exceptional high-rate cycling performance (86.3% capacity retention after 1000 cycles at 2C) were observed. These findings pave the way for the advancement and commercialization of high-performance ASSLBs.
Presented molecules combine two features, the bent-core rigidity with the dimeric character, and reveal variety of specific phases. Longer homologues exhibit tilted smectic phases, short homologues manifest a propensity for helical structures.
Low-exergy waste heat, which constitutes the majority of industrial-scale thermal losses, remains largely unrecoverable with conventional technologies. Thermoelectrics offer a solid-state solution for converting this hard-to-access energy into electricity, making them attractive for decentralized power generation and sensor applications. High-entropy materials (HEMs) have gained traction as a strategy for better-performing thermoelectrics, but the mechanisms driving their benefits require further exploration. This article highlights key insights for heat and electronic transport in HEMs. For heat transport, we argue that reduced, and often ultralow, lattice thermal conductivity in HEMs—with respect to ordered counterparts—can be taken for granted, emerging naturally as a fifth core effect of high-entropy systems. While band convergence is often considered beneficial for electronic transport, its impact depends strongly on the electronic structure. We summarize the scenarios where it can be detrimental to thermoelectric performance. These insights motivate strategies that align seamlessly with advancements in artificial intelligence and data-driven approaches, helping accelerate the discovery of next-generation thermoelectric materials.
Context: Recent, nonlinear simulations of wave generation and propagation in full-star models have been carried out in the anelastic approximation using spectral methods. Although it makes long time steps possible, this approach excludes the physics of sound waves completely and requires rather high artificial viscosity and thermal diffusivity for numerical stability. A direct comparison with observations is thus limited. Aims: In this study, we explore the capabilities of our compressible multidimensional Seven-League Hydro (SLH) code to simulate stellar oscillations. Methods: We compare some fundamental properties of internal gravity and pressure waves in 2D SLH simulations to linear wave theory using two test cases: (1) an interval gravity wave packet in the Boussinesq limit and (2) a realistic 3 M ⊙ stellar model with a convective core and a radiative envelope. Oscillation properties of the stellar model are also discussed in the context of observations. Results: Our tests show that specialized low-Mach techniques are necessary when simulating oscillations in stellar interiors. Basic properties of internal gravity and pressure waves in our simulations are in good agreement with linear wave theory. As compared to anelastic simulations of the same stellar model, we can follow internal gravity waves of much lower frequencies. The temporal frequency spectra of velocity and temperature are flat and compatible with the observed spectra of massive stars. Conclusion: The low-Mach compressible approach to hydrodynamical simulations of stellar oscillations is promising. Our simulations are less dissipative and require less luminosity boosting than comparable spectral simulations. The fully-compressible approach allows for the coupling of gravity and pressure waves in the outer convective envelopes of evolved stars to be studied in the future.
There is strong observational evidence that the convective cores of intermediate-mass and massive main sequence stars are substantially larger than those predicted by standard stellar-evolution models. However, it is unclear what physical processes cause this phenomenon or how to predict the extent and stratification of stellar convective boundary layers. Convective penetration is a thermal-timescale process that is likely to be particularly relevant during the slow evolution on the main sequence. We use our low-Mach-number S EVEN -L EAGUE H YDRO code to study this process in 2.5D and 3D geometries. Starting with a chemically homogeneous model of a 15 M ⊙ zero-age main sequence star, we construct a series of simulations with the luminosity increased and opacity decreased by the same factor, ranging from 10 3 to 10 6 . After reaching thermal equilibrium, all of our models show a clear penetration layer; its thickness becomes statistically constant in time and it is shown to converge upon grid refinement. The penetration layer becomes nearly adiabatic with a steep transition to a radiative stratification in simulations at the lower end of our luminosity range. This structure corresponds to the adiabatic ‘step overshoot’ model often employed in stellar-evolution calculations. The simulations with the highest and lowest luminosity differ by less than a factor of two in the penetration distance. The high computational cost of 3D simulations makes our current 3D data set rather sparse. Depending on how we extrapolate the 3D data to the actual luminosity of the initial stellar model, we obtain penetration distances ranging from 0.09 to 0.44 pressure scale heights, which is broadly compatible with observations.
Radio-mode feedback associated with the active galactic nuclei (AGNs) at the cores of galaxy clusters injects a large amount of energy into the intracluster medium (ICM), offsetting radiative losses through X-ray emission. This mechanism prevents the ICM from rapidly cooling down and fueling extreme starburst activity as it accretes onto the central galaxies, and it is therefore a key ingredient in the evolution of galaxy clusters. However, the influence and mode of feedback at high redshifts (z ∼ 1) remains largely unknown. Low-frequency sub-arcsecond-resolution radio observations taken with the International LOFAR Telescope have demonstrated their ability to assist X-ray observations with constraining the energy output from the AGNs (or “cavity power”) in galaxy clusters, thereby enabling research at higher redshifts than before. In this pilot project, we tested this hybrid method on a high-redshift (0.6 < z < 1.3) sample of 13 galaxy clusters for the first time with the aim of verifying the performance of this method at these redshifts and providing the first estimates of the cavity power associated with the central AGN for a sample of distant clusters. We were able to detect clear radio lobes in three out of 13 galaxy clusters at redshifts of 0.7 < z < 0.9, and we used these detections in combination with ICM pressures surrounding the radio lobes obtained from standard profiles to calculate the corresponding cavity powers of the AGNs. Combining our results with the literature, the current data appear to suggest that the average cavity power peaked at a redshift ofz ∼ 0.4 and slowly decreases toward higher redshifts. However, we require more and tighter constraints on the cavity volume and a better understanding of our observational systematics to confirm any deviation of the cavity power trend from a constant level.
Context. While supernova remnants (SNRs) have been considered the most relevant Galactic cosmic ray (CR) accelerators for decades, core-collapse supernovae (CCSNe) could accelerate particles during the earliest stages of their evolution and hence contribute to the CR energy budget in the Galaxy. Some SNRs have indeed been associated with TeV γ -rays, yet proton acceleration efficiency during the early stages of an SN expansion remains mostly unconstrained. Aims. The multi-wavelength observation of SN 2023ixf, a Type II supernova (SN) in the nearby galaxy M 101 (at a distance of 6.85 Mpc), opens the possibility to constrain CR acceleration within a few days after the collapse of the red super-giant stellar progenitor. With this work, we intend to provide a phenomenological, quasi-model-independent constraint on the CR acceleration efficiency during this event at photon energies above 100 MeV. Methods. We performed a maximum-likelihood analysis of γ -ray data from the Fermi Large Area Telescope up to one month after the SN explosion. We searched for high-energy, non-thermal emission from its expanding shock, and estimated the underlying hadronic CR energy reservoir assuming a power-law proton distribution consistent with standard diffusive shock acceleration. Results. We do not find significant γ -ray emission from SN 2023ixf. Nonetheless, our non-detection provides the first limit on the energy transferred to the population of hadronic CRs during the very early expansion of a CCSN. Conclusions. Under reasonable assumptions, our limits would imply a maximum efficiency on the CR acceleration of as low as 1%, which is inconsistent with the common estimate of 10% in generic SNe. However, this result is highly dependent on the assumed geometry of the circumstellar medium, and could be relaxed back to 10% by challenging spherical symmetry. Consequently, a more sophisticated, inhomogeneous characterisation of the shock and the progenitor’s environment is required before establishing whether or not Type II SNe are indeed efficient CR accelerators at early times.
Next generation High-Energy Physics (HEP) experiments are presented with significant computational challenges, both in terms of data volume and processing power. Using compute accelerators, such as GPUs, is one of the promising ways to provide the necessary computational power to meet the challenge. The current programming models for compute accelerators often involve using architecture-specific programming languages promoted by the hardware vendors and hence limit the set of platforms that the code can run on. Developing software with platform restrictions is especially unfeasible for HEP communities as it takes significant effort to convert typical HEP algorithms into ones that are efficient for compute accelerators. Multiple performance portability solutions have recently emerged and provide an alternative path for using compute accelerators, which allow the code to be executed on hardware from different vendors. We apply several portability solutions, such as Kokkos, SYCL, C++17 std::execution::par, Alpaka, and OpenMP/OpenACC, on two mini-apps extracted from the mkFit project: p2z and p2r. These apps include basic kernels for a Kalman filter track fit, such as propagation and update of track parameters, for detectors at a fixed z or fixed r position, respectively. The two mini-apps explore different memory layout formats. We report on the development experience with different portability solutions, as well as their performance on GPUs and many-core CPUs, measured as the throughput of the kernels from different GPU and CPU vendors such as NVIDIA, AMD and Intel.
Valence-to-core x-ray emission spectroscopy (VtC XES) combines the sample flexibility and element specificity of hard x-rays with the chemical environment sensitivity of valence spectroscopy. Here, we extend this technique to study geometric and electronic structural changes induced by photoexcitation in the femtosecond time domain via laser-pump, x-ray probe experiments using an x-ray free electron laser. The results of time-resolved VtC XES on a series of ferrous complexes [Fe(CN) 2n (2, 2'-bipyridine) 3-n ] -2n+2 , n = 1, 2, 3, are presented. Comparisons of spectra obtained from ground state density functional theory calculations reveal signatures of excited state bond length and oxidation state changes. An oxidation state change associated with a metal-to-ligand charge transfer state with a lifetime of less than 100 fs is observed, as well as bond length changes associated with metal-centered excited states with lifetimes of 13 ps and 250 ps.
In this work, we test the efficacy of excited state mean field theory and its excited-state-specific perturbation theory on the prediction of K-edge positions and x-ray peak separations. We find that the mean field theory is surprisingly accurate, even though it contains no accounting of differential electron correlation effects. In the perturbation theory, we test multiple core-valence separation schemes and find that, with the mean field theory already so accurate, electron-counting biases in one popular separation scheme become a dominant error when predicting K-edges. Happily, these appear to be relatively easy to correct for, leading to a perturbation theory for K-edge positions that is lower scaling and more accurate than coupled cluster theory and competitive in accuracy with recent high-accuracy results from restricted open-shell Kohn-Sham theory. For peak separations, our preliminary data show excited state mean field theory to be exceptionally accurate, but more extensive testing will be needed to see how it and its perturbation theory compare to coupled cluster peak separations more broadly.
Two existing particle-in-cell gyrokinetic codes, GEM for the core region and XGC for the edge region, have been successfully coupled with a spatial coupling scheme at the interface in a toroidal geometry. Additionally, a mapping technique is developed for transferring data between GEM's structured and XGC's unstructured meshes. Two examples of coupled simulations are presented to demonstrate the coupling scheme. The optimization of GEM for graphics processing unit is also presented.
Two-dimensional transition metal carbides, carbonitrides, and nitrides called MXenes exhibit high metallic conductivity, ion intercalation capability and reversible redox activity, prompting their applications in energy storage and conversion, electromagnetic interference (EMI) shielding, and electronics, among many other fields. It has been shown that replacement of about 50% of carbon atoms in the most popular MXene family member, titanium carbide (Ti 3 C 2 T x ), by nitrogen atoms, forming titanium carbonitride (Ti 3 CNT x ), leads to drastically different properties, such as very high negative charge in solution and extreme EMI shielding effectiveness, exceeding all known materials, even metals at comparable thicknesses. Here, by using ultraviolet photoemission spectroscopy (UPS), the electronic structures of Ti 3 CNT x and Ti 3 C 2 T x are systematically investigated and compared as a function of charge carrier density. We observe that, in contrast to Ti 3 C 2 T x , the Ti 3p core-level of Ti 3 CNT x exhibits a counterintuitive shift to a lower binding energy of up to approximately 250 meV upon increasing the electron density, which is a spectroscopic signature of negative electronic compressibility (NEC). These experimentally measured chemical potential shifts are well-captured by the density functional theory (DFT) calculation. The DFT results also further suggest that the hybridization of titanium-nitrogen bonding in Ti 3 CNT x helps promoting the available states of Ti atoms for receiving more electron above the Fermi level and leads to the observed NEC. Furthermore, our findings explain the differences in electronic properties between the two very important and widely studied MXenes and also suggest a new strategy to apply the NEC effect of Ti 3 CNT x in energy and charge storage applications.
Magnetization of clusters is often simulated using atomistic spin dynamics for a fixed lattice. Coupled spinlattice dynamics simulations of the magnetization of nanoparticles have up to now neglected the change in the size of the atomic magnetic moments near surfaces. We show that the introduction of variable magnetic moments leads to a better description of experimental data for the magnetization of small Fe nanoparticles. To this end we divide atoms into a surface-near shell and a core with bulk properties. Here, it is demonstrated that both the magnitude of the shell magnetic moment and the exchange interactions need to be modified to obtain a fair representation of the experimental data. This allows for a reasonable description of the average magnetic moment versus cluster size, and also the cluster magnetization versus temperature.
With the increase in collected data volumes, either from experimental measurements or high fidelity simulations, there is an ever-growing need to develop computationally efficient tools to process, analyze, and interpret these datasets. Modal analysis techniques have gained great interest due to their ability to identify patterns in the data and extract valuable information about the system being considered. Dynamic mode decomposition (DMD) relies on elements of the Koopman approximation theory to compute a set of modes, each associated with a fixed oscillation frequency and a decay/growth rate. Extracting these details from large datasets can be computationally expensive due to the need to implement singular value decomposition of the input data matrix. Sketching algorithms have become popular in numerical linear algebra where statistical theoretic approaches are utilized to reduce the cost of major operations. A sketch of a matrix is another matrix, which is significantly smaller, but still sufficiently approximates the original system. We put forth an efficient DMD framework, SketchyDMD, based on a core sketching algorithm that captures information about the range and corange (their mutual relationship) of input data. The proposed sketching-based framework can accelerate various portions of the DMD routines, compared to classical methods that operate directly on the raw input data. We conduct numerical experiments using the spherical shallow water equations as a prototypical model in the context of geophysical flows. In conclusion, we show that the proposed SketchyDMD is superior to existing randomized DMD methods that are based on capturing only the range of the input data.
Three-dimensional effects play a crucial role during the hot-spot formation in inertial confinement fusion (ICF) implosions. A data analysis technique for 3D hot-spot reconstruction from experimental observables has been developed to characterize the effects of low modes on 3D hot-spot formations. In nuclear measurements, the effective flow direction, governed by the maximum eigenvalue in the velocity variance of apparent ion temperatures, has been found to agree with the measured hot-spot flows for implosions dominated by mode ℓ = 1. Asymmetries in areal-density (ρR) measurements were found to be characterized by a unique cosine variation along the hot-spot flow axis. In x-ray images, a 3D hot-spot x-ray emission tomography method was developed to reconstruct the 3D hot-spot plasma emissivity using a generalized spherical-harmonic Gaussian function. The gradient-descent algorithm was used to optimize the mapping between the projections from the 3D hot-spot emission model and the measured x-ray images along multiple views. Furthermore, this work establishes a platform to analyze 3D low-mode core asymmetries in ICF.