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.
Not Available
Monte Carlo (MC) numerical simulation is carried out for the intensity of multiply backscattered radiation as a function of the source–receiver distance for the Henyey–Greenstein (HG) and Rayleigh–Gans (RG) anisotropic scattering models. It is found that, in spite of the standard diffusive regime assumption on the universality of the description of scattering in terms of the extinction coefficient and the transport length, the multiple scattering intensity depends on the form of the phase function and on the degree of anisotropy even for systems with identical extinction coefficients and transport lengths. In this case, the results obtained by the RG phase function are more sensitive to the scattering ansisotropy than those obtained by the HG function.
In this paper, we study the topological Hall (TH) effect in a diffusive ferromagnetic metal thin film by solving a Boltzmann transport equation in the presence of spin-flip scattering. A generalized spin-diffusion equation is derived which contains an additional source term associated with the gradient of the emergent magnetic field that arises from skyrmions. Because of the source term, spin accumulation may build up in the vicinity of the skyrmions. This gives rise to a spin-polarized diffusion current that in general suppresses the bulk TH current. Only when the spin-diffusion length is much smaller than the skyrmion size does the TH resistivity approach the value derived by Bruno et al. [Phys. Rev. Lett. 93, 096806 (2004)]. Finally, we derive a general expression of the TH resistivity that applies to thin-film geometries with spin-flip scattering, and show that the corrections to the TH resistivity become large when the size of room temperature skyrmions is further reduced to tens of nanometers.
Electron beam lithography (EBL) of chemically amplified resists (CARs) faces fundamental challenges, including stochastic electron scattering and acid diffusion, that limit resolution and reproducibility. Using SU-8 as a model CAR, this study systematically investigated complementary strategies to address these challenges, combining multipass exposure, proximity effect correction (PEC) with midrange correction factors, base quencher incorporation, and post-exposure bake (PEB) suppression. Monte Carlo simulations and calibrated PEC modeling revealed that extending the point spread function to include a midrange scattering component significantly improved critical dimension (CD) control across varying pattern densities, correcting deviations that conventional two-term PEC failed to capture. Multipass exposure, particularly 4-pass writing with a 25% offset, redistributed the dose to average stochastic beam and scattering fluctuations, reducing line-width roughness by more than 50% and yielding more uniform nanoscale features. Photoacid confinement was investigated by adding urea as a base quencher, which successfully reduced acid diffusion but introduced substantial sensitivity penalties without improving ultimate resolution or Z-factor performance, underscoring the trade-offs of chemical versus physical confinement. Suppressing PEB most directly minimized acid diffusion, resulting in improved Z-factors and reproducible 30 nm half-pitch dense line/space patterns. Overall, these results demonstrated that PEC with midrange correction, multipass strategies, quencher additives, and PEB-free processing addresses different aspects of the EBL process window and that their integration provides a comprehensive framework for managing stochastic scattering, diffusion, and chemical amplification effects. This framework advances dense nanoscale patterning in CARs and establishes guiding principles for optimizing resist design and process strategies in high-resolution EBL and potentially other advanced lithographies, such as extreme ultraviolet (EUV) lithography.
Here, we investigated the impact of Li stoichiometry and host flexibility on Li + diffusion processes in LiAlGeO 4 at the microscopic level using quasielastic neutron scattering (QENS) and ab initio molecular dynamics (AIMD) simulations. Using sufficiently long AIMD trajectories, we could simulate the observed QENS signal and identify the localized dynamics of Li in crystalline LiAlGeO 4 . Such information is vital to identify the bottleneck of diffusion processes and design materials for battery application. Our AIMD simulations in LiAlGeO 4 reveal that the Li + conductivity can be significantly improved by manipulating the Li stoichiometry and/or host flexibility via amorphization. We determined that excess Li stoichiometry enhances the Coulomb repulsion of neighboring Li sites and softens the host structure to enable faster Li + diffusion along the hexagonal c-axis. In the amorphous structure, random orientations of AlO 4 and GeO 4 polyhedral units create a wide distribution of intersite distances and significantly soften the host structure, greatly enhancing the Li + diffusion. The simulations are used to understand the nature of diffusion, especially the role of the host structure, the possible hopping pathways, and the diffusion behavior in various structures of LiAlGeO 4 .
Highlights: • Diffuse X-ray scattering collected from 4 compositions of Na{sub ½}Bi{sub ½}TiO{sub 3}-x%BaTiO{sub 3}. • Oxygen octahedral tilts, atomic displacements and modulations existed. • All local disordered structures evolved across the morphotropic phase boundary. - Abstract: Disordered structures contribute to the relaxor ferroelectric properties of the lead-free system Na{sub ½}Bi{sub ½}TiO{sub 3}-(x)%BaTiO{sub 3}. Such structures are difficult to characterize with conventional structural analysis techniques. Here, diffuse X-ray scattering was used to analyze the disordered structures from four single-crystals of Na{sub ½}Bi{sub ½}TiO{sub 3}-(x)%BaTiO{sub 3}, with compositions of (x)% = 0%, 4%, 9% and 17%. Such a range of compositions allowed a distinction between the structures that occur across the morphotropic phase boundary, where enhanced electro-mechanical properties have been reported. Throughout the compositional range, the resulting diffuse scattering showed the modification of correlated atomic displacements, transformation of planar defects in the oxygen octahedral tilt-systems, and the presence of modulations of approximately 10 unit cells. The results reported here demonstrate an evolution of the short-range structures across the morphotropic phase boundary that may underpin the electro-mechanical properties of Na{sub ½}Bi{sub ½}TiO{sub 3}-(x)%BaTiO{sub 3}. This information will ultimately assist in tailoring disordered structures for next-generation piezoelectric materials.
We have performed quasielastic neutron scattering (QENS) experiments up to 1243 K and ab initio molecular dynamics (AIMD) simulations to investigate the Na diffusion in various phases of NaAlSiO4 (NASO), namely, low-carnegieite (L-NASO; trigonal), high-carnegieite (H-NASO; cubic) and nepheline (N-NASO; hexagonal) phases. The QENS measurements reveal Na ions localized diffusion behavior in L-NASO and N-NASO, but long-range diffusion behavior in H-NASO. Furthermore, the AIMD simulation supplemented the QENS measurements and showed that excess Na ions in H-NASO enhance the host network flexibility and activate the AlO 4 /SiO 4 tetrahedra rotational modes. These framework modes enable the long-range diffusion of Na across a pathway of interstitial sites. The simulations also show Na diffusion in Na-deficient N-NASO through vacant Na sites along the hexagonal c-axis.
The determination of effective diffusion coefficients of gases or solutes in the water–saturated pore space of mudrocks is time consuming and technically challenging. Yet, reliable values of effective diffusion coefficients are important to predict migration of hydrocarbon gases in unconventional reservoirs, dissipation of (explosive) gases through clay barriers in radioactive waste repositories, mineral alteration of seals to geological CO 2 storage reservoirs and contaminant migration through aquitards. In this study, small angle and very small angle neutron scattering techniques have been utilized to determine a range of transport properties in mudrocks, including porosity, pore size distributions and surface and volume fractal dimensions of pores and grains, from which diffusive transport parameters can be estimated. Using a fractal model derived from Archie's Law, we calculate effective diffusion coefficients from these parameters and compare them to laboratory–derived effective diffusion coefficients for CO 2 , H 2 , CH 4 and HTO on either the same or related mudrock samples. The samples include Opalinus Shale from the underground laboratory in Mont Terri, Switzerland; Boom Clay from a core drilled in Mol, Belgium and a marine claystone cored in Utah, USA. The predicted values were compared to laboratory diffusion measurements. The measured and modeled diffusion coefficients show good agreement, differing generally by less than factor 5. Furthermore, neutron or X–ray scattering analysis is therefore proposed as a novel method for fast, accurate estimation of effective diffusion coefficients in mudrocks, together with simultaneous measurement of multiple transport parameters including porosity, pore size distributions and surface areas, important for (reactive) transport modelling.
We explore the possibility for orientation recovery in single-molecule coherent diffractive imaging with diffusion map. This algorithm approximates the Laplace-Beltrami operator, which we diagonalize with a metric that corresponds to the mapping of Euler angles onto scattering images. While suitable for images of objects with specific properties we show why this approach fails for realistic molecules. Here, we introduce a modification of the form factor in the scattering images which facilitates the orientation recovery and should be suitable for all recovery algorithms based on the distance of individual images. (C) 2016 Optical Society of America
We derive analytical expressions for the spectral moments of the dynamical response functions of the Hubbard model using the high-temperature series expansion. We consider generic dimension d as well as the infinite- d limit, arbitrary electron density n, and both finite and infinite repulsion U. We use moment-reconstruction methods to obtain the one-electron spectral function, the self-energy, and the optical conductivity. They are all smooth functions at high temperature and, at large U, they are featureless with characteristic widths of the order of the lattice hopping parameter t. In the infinite-d limit, we compare the series expansion results with accurate numerical renormalization group and interaction expansion quantum Monte Carlo results. We find excellent agreement down to surprisingly low temperatures, throughout most of the bad-metal regime, which applies for T≳ (1 - n)D, the Brinkman-Rice scale. The resistivity increases linearly in T at high temperature without saturation. This results from the 1/T behavior of the compressibility or kinetic energy, which play the role of the effective carrier number. In contrast, the scattering time (or diffusion constant) saturates at high T. We find that σ( n,T) ≈ (1 - n)σ(n = 0,T) to a very good approximation for all n, with σ(n = 0,T) ∝ t/T at high temperatures. The saturation at small n occurs due to a compensation between the density dependence of the effective number of carriers and that of the scattering time. The T dependence of the resistivity displays a kneelike feature which signals a crossover to the intermediate-temperature regime where the diffusion constant (or scattering time) starts increasing with decreasing T. At high temperatures, the thermopower obeys the Heikes formula, while the Wiedemann-Franz law is violated with the Lorenz number vanishing as 1/T 2 . The relevance of our calculations to experiments probing high-temperature transport in materials with strong electronic correlations or ultracold atomic gases in an optical lattice is briefly discussed.
A novel, drag-diffusion Langevin method of hybrid, kinetic-fluid Coulomb scattering in plasmas is presented. Unlike previous methods, the frictional drag is always applied in the simulation frame of reference. The velocity-space diffusion is performed in the stationary-fluid frame of reference when anisotropic, and in the laboratory frame when isotropic. While the general method is mass-ratio independent, we focus on interactions of kinetic-ions and fluid-electrons to show first-order modifications to the electron velocity distribution function that are an important correction for the accurate calculation of electric resistivity. Inclusion of sub-cycling and a limit to the maximum collision frequency is shown to allow for arbitrarily large timesteps without numerical failure. Here the Langevin method is compared with a grid-based binary method and found to require a much less restrictive timestep in cases of ion–electron slowing and temperature equilibration; this finding differs from previous work and is dependent on the mass ratio.
The microscopic information of solute and solute-solvent motions can be measured using ultrafast diffuse x-ray scattering, and compared to molecular dynamics simulations, as demonstrated in studies of solvated molecular systems done in X-ray free-electron lasers. However, for the typical photon energies used in these experiments, the photon momentum transfer range was limited to lower Q values where the signals arising from different parts of the studied system overlap significantly. Recently, this limitation was greatly alleviated by extending the photon energy at the Linac Coherent Light-Source (LCLS), significantly increasing the accessible momentum transfer range. This improvement is transformative for ultrafast diffuse x-ray scattering measurements, enabling for the first time to disentangle the solute scattering difference signal, that persists to higher Q values, from the bulk solvent and solute-solvent cross-terms difference signals. In addition, it enhances the prediction ability of modeling and simulation methods such as the hybrid QM/MM approach. The extended Q range also opens the way to resolve in real-space details regarding coherent wavepackets motions beyond their center-of-mass positions. Here, we present the first results on high-energy (18keV) diffuse x-ray scattering in solution, demonstrating high fidelity time-resolved scattering and analysis of the photoexcited model photocatalysts PtPOP ([Pt2(POP)4]4-), and IrDimen ([Ir2(dimen)4]2+) in several solvents. These complexes provide ideal systems for demonstrating the ability of high-Q ultrafast scattering and QM/MM simulations to decompose ultrafast solvation dynamics into specific changes in the solute-solvent pair distribution function with a particular focus on how electronically excited states change the interaction between the solvent and photo-catalytically active metal sites. We introduce an approach to further utilize the high-energy capability and demonstrate a single-shot ultra-wide-angle X-ray Scattering modality using two perpendicular detectors, spanning a scattering angle range of more than 100 degrees, allowing to extend the accessible momentum transfer range up to Q~14 Å-1. We develop a model-free method to invert the x-ray scattering signals and enable the recovery of multiple pair-density motions that happen simultaneously, allowing us to directly measure in real-space nuclear wavepacket motions. . [1] Natan, Adi. "Real-Space Inversion and Super-Resolution of Ultrafast X-ray Scattering using Natural Scattering Kernels." arXiv preprint arXiv:2107.05576 (2021)
To characterize in detail the charge density wave (CDW) transition of 1𝑇−VSe 2 , its electronic structure and lattice dynamics are comprehensively studied by means of x-ray diffraction, muon spectroscopy, angle resolved photoemission (ARPES), diffuse and inelastic x-ray scattering, and state-of-the-art first-principles density functional theory calculations. Resonant elastic x-ray scattering does not show any resonant enhancement at either V or Se, indicating that the CDW peak at the 𝐾 edges describes a purely structural modulation of the electronic ordering. ARPES experiments identify (i) a pseudogap at 𝑇 > 𝑇 CDW , which leads to a depletion of the density of states in the ML-M'L' plane at 𝑇 < 𝑇 CDW , and (ii) anomalies in the electronic dispersion reflecting a sizable impact of phonons on it. A diffuse scattering precursor, characteristic of soft phonons, is observed at room temperature (RT) and leads to the full collapse of the low-energy phonon (𝜔 1 ) with propagation vector (0.25 0 −0.3) r.l.u. Here, we show that the frequency and linewidth of this mode are anisotropic in momentum space, reflecting the momentum dependence of the electron-phonon interaction (EPI), hence demonstrating that the origin of the CDW is, to a much larger extent, due to the momentum dependent EPI with a small contribution from nesting. The pressure dependence of the 𝜔 1 soft mode remains nearly constant up to 13 GPa at RT, with only a modest softening before the transition to the high-pressure monoclinic C2/m phase. The wide set of experimental data is well captured by our state-of-the art first-principles anharmonic calculations with the inclusion of van der Waals corrections in the exchange-correlation functional. The comprehensive description of the electronic and dynamical properties of VSe 2 reported here adds important pieces of information to the understanding of the electronic modulations in the family of transition-metal dichalcogenides.
Sodium iron phosphate, Na x FePO 4 , is currently one of the most promising cathode materials for sodium ion batteries. Here we employ a systematic approach using three neutron scattering techniques to investigate temperature dependent properties with a focus on the charge transfer processes. Two distinct Na + diffusive motions are observed with quasielastic neutron scattering, including room temperature motion. These processes correspond to rapid double jumps through extrinsic vacancy sites of ∼6 Å along with the more fundamental nearest neighbor jumps of 3.1 Å. Experimental observations of room temperature ionic diffusion are highly desirable, as the diffusion pathways are naturally well suited for efficient transport. Inelastic neutron scattering was used to probe the vibrational density of states that assist in the ionic and electronic (polaronic) transport processes, which are further supported with density functional theory. Neutron diffraction examines the disordering of the sodium sublattice upon transition to the solid-solution phase, which is believed to enhance the overall kinetics during sodium insertion/extraction. Finally, 57 Fe Mössbauer spectroscopy measurements probes polaron hopping and the local structural evolution that modifies the electron density surrounding the redox active iron sites. This study sheds light on the temperature induced atomic scale dynamics that occur in olivine Na x FePO 4 .
ABSTRACT Understanding the transport of energetic cosmic rays belongs to the most challenging topics in astrophysics. Diffusion due to scattering by electromagnetic fluctuations is a key process in cosmic ray transport. The transition from a ballistic to a diffusive-propagation regime is presented in direct numerical calculations of diffusion coefficients for homogeneous magnetic field lines subject to turbulent perturbations. Simulation results are compared with theoretical derivations of the parallel diffusion coefficient’s dependences on the energy and the fluctuation amplitudes in the limit of weak turbulence. The present study shows that the widely used extrapolation of the energy scaling for the parallel diffusion coefficient to high turbulence levels predicted by quasi-linear theory does not provide a universally accurate description in the resonant-scattering regime. It is highlighted here that the numerically calculated diffusion coefficients can be polluted for low energies due to missing resonant interaction possibilities of the particles with the turbulence. Five reduced-rigidity regimes are established, which are separated by analytical boundaries derived in this work. Consequently, a proper description of cosmic ray propagation can only be achieved by using a turbulence-level-dependent diffusion coefficient and can contribute to solving the Galactic cosmic ray gradient problem.
The collapse or merging of individual plumes of direct-injection gasoline injectors is of fundamental importance to engine performance because of its impact on fuel–air mixing. But, the mechanisms of spray collapse are not fully understood and are difficult to predict. The purpose of this work is to study the aerodynamics in the inter-spray region, which can potentially lead to plume collapse. High-speed (100 kHz) particle image velocimetry is applied along a plane between plumes to observe the full temporal evolution of plume interaction and potential collapse, resolved for individual injection events. Supporting information along a line of sight is obtained using simultaneous diffused back illumination and Mie-scatter techniques. Experiments are performed under simulated engine conditions using a symmetric eight-hole injector in a high-temperature, high-pressure vessel at the “Spray G” operating conditions of the engine combustion network. Indicators of plume interaction and collapse include changes in counter-flow recirculation of ambient gas toward the injector along the axis of the injector or in the inter-plume region between plumes. Furthermore, the effect of ambient temperature and gas density on the inter-plume aerodynamics and the subsequent plume collapse are assessed. Increasing ambient temperature or density, with enhanced vaporization and momentum exchange, accelerates the plume interaction. Plume direction progressively shifts toward the injector axis with time, demonstrating that the plume interaction and collapse are inherently transient.
Here, we investigate the gradual diffusion of energetic electrons from the inner edge of the outer radiation belt into the slot region. The Van Allen Probes observed slow inward diffusion and decay of ~200–600 keV electrons following the intense geomagnetic storm that occurred on 17 March 2013. During the 10 day nondisturbed period following the storm, the peak of electron fluxes gradually moved from L ~ 2.7 to L ~ 2.4, and the flux levels decreased by a factor of ~2–4 depending on the electron energy. We simulated the radial intrusion and decay of electrons using a three–dimensional diffusion code, which reproduced the energy–dependent transport of electrons from ~100 keV to 1 MeV in the slot region. At energies of 100–200 keV, the electrons experience fast transport across the slot region due to the dominance of radial diffusion; at energies of 200–600 keV, the electrons gradually diffuse and decay in the slot region due to the comparable rate of radial diffusion and pitch angle scattering by plasmaspheric hiss; at energies of E > 700 keV, the electrons stopped diffusing near the inner edge of outer radiation belt due to the dominant pitch angle scattering loss. In addition to plasmaspheric hiss, magnetosonic waves and VLF transmitters can cause the loss of high pitch angle electrons, relaxing the sharp “top–hat” shaped pitch angle distributions created by plasmaspheric hiss. Our simulation indicates the importance of balance between radial diffusion and loss through pitch angle scattering in forming the diffusive intrusion of energetic electrons across the slot region.