Search NASA⌕ Search

SEARCH · Search NASA

Results for “High Order”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Two-dimensional phase cartography for high-harmonic spectroscopy

Extremely nonlinear spectroscopy based on high-order-harmonic generation has become a powerful investigation method for attosecond dynamics in gas and solid targets. In particular, the phase of harmonic emission was shown to carry profound insight into atomic and molecular structure and dynamics. However, current techniques offer phase measurements only along specific directions, thus providing partial characterization. Here we report on a new approach combining optical and quantum interferometers measuring along two dimensions the intensity and phase of harmonic emission from aligned molecules in the exact same experimental conditions. This two-dimensional cartography technique measures the phase with no arbitrary offset and no uncertainty on its sign. Measurements along different dimensions can be combined in two ways: either a single mapping or a redundant mapping allowing high-precision phase recovery using a Shack–Hartmann-like algorithm. We demonstrate both methods in a nitrogen test case, which allows disentangling structural and dynamical effects. Two-dimensional phase cartography paves the way to high-resolution high-harmonic spectroscopy for applications such as quantum orbital tomography and attosecond charge migration in molecules.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Slimmer Geminals For Accurate F12 Electronic Structure Models

The Slater-type F12 geminal length scales originally tuned for the second-order Mo̷ller-Plesset F12 method are too large for higher-order F12 methods formulated using the SP (diagonal fixed-coefficient spin-adapted) F12 ansatz. The new geminal parameters reported herein reduce the basis set incompleteness errors (BSIEs) of absolute coupled-cluster singles and doubles F12 correlation energies by a significant─and increase with the cardinal number of the basis─margin. The effect of geminal reoptimization is especially pronounced for the cc-pVXZ-F12 basis sets (specifically designed for use with F12 methods) relative to their conventional aug-cc-pVXZ counterparts. The BSIEs of relative energies are less affected, but substantial reductions can be obtained, especially for atomization energies and ionization potentials with the cc-pVXZ-F12 basis sets. The new geminal parameters are therefore recommended for all applications of high-order F12 methods, such as coupled-cluster F12 methods and transcorrelated F12 methods.

Powell, Samuel R. [Virginia Polytechnic Inst. and ↗
DOE OSTI2021

On the possibility of footprint compression with one lens in nonlinear accelerator lattice

Electromagnetic interaction of colliding beams along with other nonlinear fields often limits the beams' lifetimes and luminosities. Nonlinearities result in the spread of betatron frequencies (footprint) and, thus, may enhance dynamic diffusion of particles due to high order resonances. One of the possible ways to eliminate nonlinearities and overcome the corresponding difficulties is compensation of nonlinear forces, but, in practice, it is hardly possible to obtain exact linearity of the system. The compensation with a single nonlinear lens cannot cope with distributed nonlinearities, nonlinearities due to parasitic crossings, etc. Here, we present a method to compute parameters of nonlinear element (lens) that eliminates both the footprint and resonance strength without achieving full compensation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Single-Crystal Polycationic Polymers Obtained by Single-Crystal-to-Single-Crystal Photopolymerization

The efficient preparation of single-crystalline ionic polymers and fundamental understanding of their structure–property relationships at the molecular level remains a challenge in chemistry and materials science. Herein, we observe the single-crystal structure of a highly ordered polycationic polymer (polyelectrolyte) and its proton conductivity. The polyelectrolyte single crystals can be prepared on a gram-scale in quantitative yield, by taking advantage of an ultraviolet/sunlight-induced topochemical polymerization, from a tricationic monomer—a self-complementary building block possessing a preorganized conformation. A single-crystal-to-single-crystal photopolymerization was revealed unambiguously by in situ single-crystal X-ray diffraction analysis, which was also employed to follow the progression of molecular structure from the monomer, to a partially polymerized intermediate, and, finally, to the polymer itself. Collinear polymer chains are held together tightly by multiple Coulombic interactions involving counterions to form two-dimensional lamellar sheets (1 nm in height) with sub-nanometer pores (5 Å). The polymer is extremely stable under 254 nm light irradiation and high temperature (above 500 K). The extraordinary mechanical strength and environmental stability—in combination with its impressive proton conductivity (~3 × 10 –4 S cm –1 )—endow the polymer with potential applications as a robust proton-conducting material. By marrying supramolecular chemistry with macromolecular science, the outcome represents a major step toward the controlled synthesis of single-crystalline polyelectrolyte materials with perfect tacticity.

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

Architected mesoporous crystalline magnesium silicates with ordered pore structures

Novel approaches to harness earth abundant silicates as building blocks for carbon dioxide removal, capture, utilization, and storage are gaining increasing attention in the context of sustainable and low carbon energy and resource recovery. Advancing a calibrated understanding of these fluid-silicate interactions is essential for developing scalable processes. One of the challenges in developing predictive controls over these interactions is the compositional and morphological heterogeneity of naturally occurring, heterogeneous magnesium silicate minerals. To address this challenge, the synthesis of architected mesoporous crystalline magnesium silicate (Mg 2 SiO 4 ) is proposed. While synthesis routes for producing amorphous mesoporous magnesium silicates have been developed via sol-gel methods, approaches to synthesize crystalline magnesium silicates with well-controlled pore size distributions have not been explored. The conventional approaches of converting matter that is amorphous to crystalline states at elevated temperatures results in a heterogeneous pore size distribution. To develop controls on pore size distribution, amorphous mesoporous magnesium silicates are coated with carbon. This approach preserves the pore size distributions during the amorphous to crystalline transformations of Mg-silicates at elevated temperatures. The carbon coating is removed on heating. Magnesium silicate particles produced using this approach have highly ordered pores around 2.58 nm and a specific surface area of 124.25 m 2 /g. In this study, we report the chemical compositions, morphologies, phase transitions, and pore structures of the intermediate materials produced during the synthesis of crystalline mesoporous magnesium silicates. The transitions, and pore structures of the intermediate materials produced during the synthesis of crystalline mesoporous magnesium silicates. Furthermore, the synthesis routes discussed in this study can be applied translationally to produce metal silicates with ordered mesoporous structures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Direct, simple and efficient computation of all components of the virtual-casing magnetic field in axisymmetric geometries with Kapur–Rokhlin quadrature

In a recent publication (Toler et al., J. Plasma Phys., vol. 89, issue 2, 2023, p. 905890210), we demonstrated that for axisymmetric geometries, the Kapur–Rokhlin quadrature rule provided an efficient and high-order accurate method for computing the normal component, on the plasma surface, of the magnetic field due to the toroidal current flowing in the plasma, via the virtual-casing principle. The calculation was indirect, as it required the prior computation of the magnetic vector potential from the virtual-casing principle, followed by the computation of its tangential derivative by Fourier differentiation, to obtain the normal component of the magnetic field. Our approach did not provide the other components of the virtual-casing magnetic field. In this letter, we show that a more direct and more general approach is available for the computation of the virtual-casing magnetic field. The Kapur–Rokhlin quadrature rule accurately calculates the principal value integrals in the expression for all the components of the magnetic field on the plasma boundary, and the numerical error converges at a rate nearly as high as the indirect method we presented previously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Continuum electromagnetic gyrokinetic simulations of turbulence in the tokamak scrape-off layer and laboratory devices

Here, we present algorithms and results from Gkeyll, a full-f continuum, electromagnetic gyrokinetic code, designed to study turbulence in the edge region of fusion devices. The edge is computationally very challenging, requiring robust algorithms that can handle large-amplitude fluctuations and stable interactions with plasma sheaths. We present an energy-conserving high-order discontinuous Galerkin scheme that solves gyrokinetic equations in Hamiltonian form. Efficiency is improved by a careful choice of basis functions and automatically generated computation kernels. Previous verification tests were performed in the straight-field-line large plasma device [Shi et al., J. Plasma Phys. 83, 905830304 (2017)] and the Texas Helimak, a simple magnetized torus [Bernard et al., Phys. Plasmas 26, 042301 (2019)], including the effect of end-plate biasing on turbulence. Results for the scrape-off layer for NSTX parameters with a model helical magnetic geometry with bad curvature have been obtained [Shi et al., Phys. Plasmas 26, 012307 (2019)]. In this paper, we present algorithms for the two formulations of electromagnetic gyrokinetics: the Hamiltonian and the symplectic. We describe each formulation and show results of benchmark tests. Although our scheme works for the Hamiltonian formulation, the presence of spurious numerical modes for high-β and large $k^2_⊥ρ^2_s$ regimes shows that the symplectic formulation is more robust. We then review our recent algorithm for the symplectic formulation [Mandell et al., J. Plasma Phys. 86, 905860109 (2020)], along with example application of this new capability. Maintaining positivity of the distribution function can be challenging, and we describe a new and novel exponential recovery based algorithm to address this.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Large-Eddy Simulation of a wind turbine using a Filtered Actuator Line Model

When dealing with multirotor devices such as quadcopters or wind farms, the cost of blade-resolved large-eddy simulation (LES) becomes prohibitive. Combining LES with a family of lower-fidelity models, called actuator line models (ALMs), has grown in popularity in the past decade. ALM replaces full blade resolution with an array of actuator points or lines parameterized by aerodynamic lift/drag polar plots along the blades. Body forces computed based on these actuator points are then projected onto the LES flow mesh, mimicking the effect of rotating blades on the flow. However, the optimal projection radius and the associated LES grid size is often too restrictive for multirotor simulations. Recently, a new tip-correction-based filtered ALM (F-ALM) was proposed by Martinez-Tossas and Meneveau (2019), which allows coarser-than-optimal grids by avoiding the associated overprediction of thrust. In this work, F-ALM is implemented into a high-order, in-house LES code to simulate National Renewable Energy Laboratory Phase VI wind turbine. It is then followed by a comparison between the baseline ALM and the newly implemented F-ALM in terms of instantaneous and time-averaged flow fields and blade loads, revealing the advantages of F-ALM in preventing the overprediction of power on coarse grids. Finally, this encourages accurate and affordable simulations of multirotor devices in the future.

17 WIND ENERGY↗
DOE OSTI2020

A chemically stabilized sulfur cathode for lean electrolyte lithium sulfur batteries

Lithium sulfur batteries (LSBs) are promising next-generation rechargeable batteries due to the high gravimetric energy, low cost, abundance, nontoxicity, and high sustainability of sulfur. However, the dissolution of high-order polysulfide in electrolytes and low Coulombic efficiency of Li anode require excess electrolytes and Li metal, which significantly reduce the energy density of LSBs. Quasi-solid-state LSBs, where sulfur is encapsulated in the micropores of carbon matrix and sealed by solid electrolyte interphase, can operate under lean electrolyte conditions, but a low sulfur loading in carbon matrix (<40 wt %) and low sulfur unitization (<70%) still limit the energy density in a cell level. Here, we significantly increase the sulfur loading in carbon to 60 wt % and sulfur utilization to ∼87% by dispersing sulfur in an oxygen-rich dense carbon host at a molecular level through strong chemical interactions of C–S and O–S. In an all-fluorinated organic lean electrolyte, the C/S cathode experiences a solid-state lithiation/delithiation reaction after the formation of solid electrolyte interphase in the first deep lithiation, completely avoiding the shuttle reaction. The chemically stabilized C/S composite retains a high reversible capacity of 541 mAh⋅g −1 (based on the total weight of the C/S composite) for 200 cycles under lean electrolyte conditions, corresponding to a high energy density of 974 Wh⋅kg −1 . The superior electrochemical performance of the chemical bonding-stabilized C/S composite renders it a promising cathode material for high-energy and long-cycle-life LSBs.

25 ENERGY STORAGE↗
DOE OSTI2023

Real-space solution to the electronic structure problem for nearly a million electrons

We report a Kohn–Sham density functional theory calculation of a system with more than 200 000 atoms and 800 000 electrons using a real-space high-order finite-difference method to investigate the electronic structure of large spherical silicon nanoclusters. Our system of choice was a 20 nm large spherical nanocluster with 202 617 silicon atoms and 13 836 hydrogen atoms used to passivate the dangling surface bonds. To speed up the convergence of the eigenspace, we utilized Chebyshev-filtered subspace iteration, and for sparse matrix–vector multiplications, we used blockwise Hilbert space-filling curves, implemented in the PARSEC code. For this calculation, we also replaced our orthonormalization + Rayleigh–Ritz step with a generalized eigenvalue problem step. We utilized all of the 8192 nodes (458 752 processors) on the Frontera machine at the Texas Advanced Computing Center. We achieved two Chebyshev-filtered subspace iterations, yielding a good approximation of the electronic density of states. Our work pushes the limits on the capabilities of the current electronic structure solvers to nearly 106 electrons and demonstrates the potential of the real-space approach to efficiently parallelize large calculations on modern high-performance computing platforms.

Chemistry↗
DOE OSTI2020

Ultrashort XUV pulse absorption spectroscopy of partially oxidized cobalt nanoparticles

High-order harmonic generation (HHG) based transient extreme ultraviolet (XUV) absorption spectroscopy is an emerging technique to trace photoinduced charge carrier dynamics in condensed phase materials with femtosecond and even attosecond temporal resolution and elemental specificity. However, its application to nanoparticulate samples that are relevant, for example, for novel photocatalytic light harvesting concepts, has been limited. This is in part due to the challenge to produce residual-free samples on ultrathin, XUV-transparent substrates as well as a widespread understanding that sparsely distributed nanoparticles do not provide sufficient contrast for XUV absorption measurements. Here, we present static XUV absorption spectra of partially oxidized Co nanowire-structures with diameters of approximately 4.5 nm and lengths between 10 and 40 nm, recorded with an ultrashort pulse HHG light source. Nanoparticles are synthesized by the agglomeration of Co atoms inside superfluid helium droplets, followed by surface deposition and oxidation in ambient air. The method is uniquely suited for residual-free synthesis of transition metal nanowires and their deposition on ultrathin substrates. Analysis by high-resolution transmission electron microscopy reveals the formation of CoO nanowires with regions of unoxidized Co in their interior. The nanoparticle samples are investigated in an HHG-driven ultrafast XUV absorption setup. Despite the low surface coverage of only 23%, the recorded spectrum exhibits a distinct absorption feature at the Co M 2,3 (2p) edge near 60 eV with a peak height of about 40 mOD. Finally, the results support the feasibility of table-top ultrafast transient XUV absorption studies of photoinduced dynamics in transition metal oxide nanoparticles with sub-monolayer surface coverage.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Solving the electronic structure problem for over 100000 atoms in real space

Using a real-space high-order finite-difference approach, we investigate the electronic structure of large spherical silicon nanoclusters. Within Kohn-Sham density functional theory and using pseudopotentials, we report the self-consistent field convergence of a system with over 100000 atoms: a Si 107,641 ⁢H 9,084 nanocluster with a diameter of 16 nm. Our approach uses Chebyshev-filtered subspace iteration to speed up the convergence of the eigenspace, and blockwise Hilbert space-filling curves to speed up sparse matrix-vector multiplications, all of which are implemented in the parsec code. For the largest system, we utilized 2048 nodes (114 688 cores) on the Frontera machine in the Texas Advanced Computing Center. Our quantitative analysis of the electronic structure shows how it gradually approaches its bulk counterpart as a function of nanocluster size. The band gap is enlarged due to quantum confinement in nanoclusters, but decreases as the system size increases, as expected. In conclusion, our work serves as a proof of concept for the capacity of the real-space approach in efficiently parallelizing very large calculations using high-performance computer platforms, which can straightforwardly be replicated in other systems with more than 10 5 atoms.

0-dimensional systems↗
DOE OSTI2023

Resonance Effect in Brunel Harmonic Generation in Thin Film Organic Semiconductors

Abstract Nonlinear optical properties of organic semiconductors (OSCs) have been extensively investigated in the perturbative regime, while strong light induced high‐order processes in solid‐state OSCs are less studied. Here, below‐threshold harmonic generation is examined, both experimentally and theoretically, in two solid‐state thin film OSCs, that is, tetraphenylporphyrin and zinc tetraphenylporphyrin. Results show that the π–π* excitations of the porphyrin ring system generate the harmonic emission. The contribution of the Brunel harmonic to the 5th harmonic emission is uncovered, where the resonant 5‐photon transition ( S 0 → S 2 transition) is found to lead to an early onset of non‐perturbative behavior. A similar resonance effect is expected in Brunel harmonic generation in other organic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Direct numerical simulations of turbulent reacting flows with shock waves and stiff chemistry using many-core/GPU acceleration

Compressible reacting flows may display sharp spatial variation related to shocks, contact discontinuities or reactive zones embedded within relatively smooth regions. The presence of such phenomena emphasizes the relevance of shock-capturing schemes such as the weighted essentially non-oscillatory (WENO) scheme as an essential ingredient of the numerical solver. However, these schemes are complex and have more computational cost than the simple high-order compact or non-compact schemes. In this paper, we present the implementation of a seventh-order, minimally-dissipative mapped WENO (WENO7M) scheme in a newly developed direct numerical simulation (DNS) code called KAUST Adaptive Reactive Flows Solver (KARFS). In order to make efficient use of the computer resources and reduce the solution time, without compromising the resolution requirement, the WENO routines are accelerated via graphics processing unit (GPU) computation. The performance characteristics and scalability of the code are studied using different grid sizes and block decomposition. Furthermore, the performance portability of KARFS is demonstrated on a variety of architectures including NVIDIA Tesla P100 GPUs and NVIDIA Kepler K20X GPUs. In addition, the capability and potential of the newly implemented WENO7M scheme in KARFS to perform DNS of compressible flows is also demonstrated with model problems involving shocks, isotropic turbulence, detonations and flame propagation into a stratified mixture with complex chemical kinetics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Observing neutron alignment at high angular momenta in heavy deformed actinides

An outstanding puzzle in the structure of very heavy nuclei lies in the rotation response of deformed transuranic isotopes. While the rotation alignment of pairs of 𝑖 13/2 protons is ubiquitous and well understood, observing the analogous alignment of the 𝑗 15/2 neutrons has been elusive. In this work we present a new quantitative analysis approach: from their moments of inertia as a function of frequency, we compare the alignments of rotational bands with both “low” and high-𝑗 intrinsic structures in the 𝑁 = 151 isotones of Pu, Cm, Cf, and No, thereby distinguishing between the alignments of neutron versus proton pairs through “blocking” effects. In this way, the alignment of 𝑗 15/2 neutrons can be clearly identified for the first time. The protons and neutrons begin aligning at a very similar frequency. The contribution of 𝑗 15/2 neutrons to the total alignment is found to grow with increasing 𝑍 compared to 𝑖 13/2 protons. To explain these observations, it is necessary to adopt a self-consistent set of single-particle energies corresponding to the nuclear shape that minimizes the mass. The shape description needs to include high-order deformations up to at least 𝛽 6 .

A ≥ 220↗
DOE OSTI2022

Filming movies of attosecond charge migration in single molecules with high harmonic spectroscopy

Electron migration in molecules is the progenitor of chemical reactions and biological functions after light-matter interaction. Following this ultrafast dynamics, however, has been an enduring endeavor. Here we demonstrate that, by using machine learning algorithm to analyze high-order harmonics generated by two-color laser pulses, we are able to retrieve the complex amplitudes and phases of harmonics of single fixed-in-space molecules. These complex dipoles enable us to construct movies of laser-driven electron migration after tunnel ionization of N 2 and CO 2 molecules at time steps of 50 attoseconds. Moreover, the angular dependence of the migration dynamics is fully resolved. By examining the movies, we observe that electron holes do not just migrate along the laser polarization direction, but may swirl around the atom centers. Our result establishes a general scheme for studying ultrafast electron dynamics in molecules, paving a way for further advance in tracing and controlling photochemical reactions by femtosecond lasers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Quantum-Path-Resolved Attosecond High-Harmonic Spectroscopy

Strong-field ionization of molecules releases electrons which can be accelerated and driven back to recombine with their parent ion, emitting high-order harmonics. This ionization also initiates attosecond electronic and vibrational dynamics in the ion, evolving during the electron travel in the continuum. Revealing this subcycle dynamics from the emitted radiation usually requires advanced theoretical modeling. Here we show that this can be avoided by resolving the emission from two families of electronic quantum paths in the generation process. The corresponding electrons have the same kinetic energy, and thus the same structural sensitivity, but differ by the travel time between ionization and recombination-the pump-probe delay in this attosecond self-probing scheme. We measure the harmonic amplitude and phase in aligned CO 2 and N 2 molecules and observe a strong influence of laser-induced dynamics on two characteristic spectroscopic features: a shape resonance and multichannel interference. This quantum-path-resolved spectroscopy thus opens wide prospects for the investigation of ultrafast ionic dynamics, such as charge migration.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Novel Concepts for High Gradient Acceleration (Final Technical Report)

We have conducted an intensive, pioneering program to demonstrate novel concepts for achieving high gradient acceleration at frequencies from the conventional microwave bands up to the millimeter wave /THz bands. High gradient accelerators hold the promise of smaller and less costly accelerators for applications ranging from the largest scale accelerators used for discovery science down to the smallest accelerators used for industrial, homeland security and medical applications. The research consisted of two major research thrusts: 1.) Structure-based wakefield accelerator (SWFA) research in collaboration with the Argonne Wakefield Accelerator research group and 2.) Millimeter Wave / THz high gradient acceleration in collaboration with SLAC. The specific goals of the research program were: Design novel metallic metamaterial structures that increase the beam-wave coupling for the accelerator mode and reduce the effect of high order modes; test novel metamaterial structures to achieve higher output power, > 1 GW at X-Band (11.7 GHz), in test at the Argonne Wakefield Accelerator (AWA); determine the experimental breakdown threshold for nanosecond-scale pulses at X-Band in testing at the AWA; test a 110 GHz accelerator structure with a field emission electron gun, built at SLAC, using pulses from a 1 MW, 110 GHz gyrotron; design, build and test a 110 GHz quasi-optical, resonant-ring pulse compressor to compress microsecond pulses from the 1 MW gyrotron into > 20 MW, 5 ns output pulses for accelerator structure testing. The proposed research program built on our successes in our research program including: Generation of 510 MW, 2.1 ns (FWHM) pulses at 11.7 GHz from a metallic metamaterial structure in test at the Argonne Wakefield Accelerator using a train of eight 65 MeV electron bunches spaced at 1.3 GHz with a total charge of 280 nC. The metamaterial structure consisted of 100 copper unit cells each consisting of a “wagon-wheel” plate and a spacer plate with a total structure length of 0.2 m. The 510 MW pulse generated an on-axis wakefield of 130 MV/m that could be used to accelerate a trailing witness bunch. Demonstration of coupling of an unprecedented rf power level of 575kW into a 110 GHz accelerator structure using a quasi-optical setup. The standing structure consisted of a central copper cavity located between two matching cavities fed by a TM01 mode. The 6 ns input pulses were sliced from 3 microsecond pulses from the gyrotron using a laser-driven silicon switch. We obtained an unprecedented high gradient up to 230MV/m corresponding to a peak surface electric field of more than 520 MV/m.

43 PARTICLE ACCELERATORS↗