Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “space computing”

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 793 records · Page 44

Control of monodomain polymer-stabilized cuboidal nanocrystals of chiral nematics by confinement.

Liquid crystals are important components of optical technologies. Cuboidal crystals consisting of chiral liquid crystals-the so-called blue phases (BPs), are of particular interest due to their crystalline structures and fast response times, but it is critical that control be gained over their phase behavior as well as the underlying dislocations and grain boundaries that arise in such systems. Blue phases exhibit cubic crystalline symmetries with lattice parameters in the 100 nm range and a network of disclination lines that can be polymerized to widen the range of temperatures over which they occur. Here, we introduce the concept of strain-controlled polymerization of BPs under confinement, which enables formation of strain-correlated stabilized morphologies that, under some circumstances, can adopt perfect single-crystal monodomain structures and undergo reversible crystal-to-crystal transformations, even if their disclination lines are polymerized. We have used super-resolution laser confocal microscopy to reveal the periodic structure and the lattice planes of the strain and polymerization stabilized BPs in 3D real space. Our experimental observations are supported and interpreted by relying on theory and computational simulations in terms of a free energy functional for a tensorial order parameter. Simulations are used to determine the orientation of the lattice planes unambiguously. The findings presented here offer opportunities for engineering optical devices based on single-crystal, polymer-stabilized BPs whose inherent liquid nature, fast dynamics, and long-range crystalline order can be fully exploited.

confinement↗

Attosecond response of molecules to impulsive ionization

When matter interacts with energetic radiation it can undergo sudden, or impulsive, ionization. This process can drive chemical change and occurs widely in space and planetary atmospheres, yet its comprehensive description challenges our current theoretical and computational capabilities as it requires advanced treatment of electron correlation and nonadiabatic dynamics beyond the Born–Oppenheimer approximation. Here, in this study, we measure the response of the para-aminophenol molecule to sudden ionization. Using attosecond X-ray absorption spectroscopy, we resolve the ultrafast dynamics of the ionized molecule with atomic precision. A subfemtosecond decay corresponds to states undergoing non-radiative decay, whereas few-femtosecond oscillatory signatures are associated with electronic wavepacket motion in stable cation states that later couple to nuclear motion. We compare our measurement with state-of-the-art computational modelling, qualitatively reproducing the observed response across multiple timescales. These results provide a benchmark for computational models of sudden ionization and ultrafast charge motion in matter.

Driver, Taran [SLAC National Accelerator Laborator↗

Simulation of adiabatic quantum computing for molecular ground states

Quantum computation promises to provide substantial speedups in many practical applications with a particularly exciting one being the simulation of quantum many-body systems. Adiabatic state preparation (ASP) is one way that quantum computers could recreate and simulate the ground state of a physical system. In this paper, we explore a novel approach for classically simulating the time dynamics of ASP with high accuracy and with only modest computational resources via an adaptive sampling configuration interaction scheme for truncating the Hilbert space to only the most important determinants. We verify that this truncation introduces negligible error and use this new approach to simulate ASP for sets of small molecular systems and Hubbard models. Furthermore, we examine two approaches to speeding up ASP when performed on quantum hardware: (i) using the complete active space configuration interaction (CASCI) wave function instead of the Hartree–Fock initial state and (ii) a nonlinear interpolation between the initial and target Hamiltonians. We find that starting with a CASCI wave function with a limited active space yields substantial speedups for many of the systems examined, while nonlinear interpolation does not. In additional, we observe interesting trends in the minimum gap location (based on the initial state) as well as how state preparation time can depend on certain molecular properties, such as the number of valence electrons. Importantly, we find that the required state preparation times do not show an immediate exponential wall that would preclude an efficient run of ASP on actual hardware.

Kremenetski, Vladimir↗

Kilopower Reactor Using Stirling TechnologY (KRUSTY) Component Critical Experiments

A series of critical experiments were conducted at the National Criticality Experiments Research Center (NCERC) in Nevada to evaluate the operational performance of a compact reactor that eventually will resemble the flight unit the National Aeronautics and Space Administration will use for deep space exploration. The results from the experiments are compared to preliminary results from computational models using MCNP and ENDF/B-7.1 neutron cross-section data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Simulating Diverse HEP Workflows on Heterogeneous Architectures

We evaluate key patterns and estimate throughput bounds of simulated transformation of conventional high energy physics (HEP) data processing workflows to heterogeneous equivalents. The simulation parameter space includes the number of offloaded tasks, CPU/accelerator ratios of intra-task computations, offload latencies, and run time efficiency of offloaded computations. The simulation is performed for a diverse set of state-of-the-art event reconstruction scenarios from ATLAS, LHCb and CMS - the frontier HEP experiments of the Large Hadron Collider project.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Design of detectors at the electron ion collider with artificial intelligence

Abstract Artificial Intelligence (AI) for design is a relatively new but active area of research across many disciplines. Surprisingly when it comes to designing detectors with AI this is an area at its infancy. The electron ion collider is the ultimate machine to study the strong force. The EIC is a large-scale experiment with an integrated detector that extends for about ±35 meters to include the central, far-forward, and far-backward regions. The design of the central detector is made by multiple sub-detectors, each in principle characterized by a multidimensional design space and multiple design criteria also called objectives. Simulations with Geant4 are typically compute intensive, and the optimization of the detector design may include non-differentiable terms as well as noisy objectives. In this context, AI can offer state of the art solutions to solve complex combinatorial problems in an efficient way. In particular, one of the proto-collaborations, ECCE, has explored during the detector proposal the possibility of using multi-objective optimization to design the tracking system of the EIC detector. This document provides an overview of these techniques and recent progress made during the EIC detector proposal. Future high energy nuclear physics experiments can leverage AI-based strategies to design more efficient detectors by optimizing their performance driven by physics criteria and minimizing costs for their realization.

Instruments & Instrumentation↗

Stability of the peroxide group in BaO 2 under high pressure

Alkaline earth metal peroxides are typical examples of ionic compounds containing polyanions. We herein report a stable BaO 2 phase at high pressure up to 130 GPa found via a first-principles computational structure searches and high-pressure experimental investigations. The identified monoclinic structure (space group C2/m) can be derived by sublattice distortions of Ba atoms and peroxide groups associated with the phonon mode softening of the lower-pressure Cmmm structure. Contrary to the previous expectation of polymerization of the peroxide group at elevated pressure, the new phase retains the peroxide group and interestingly exhibits an insulating behavior demonstrating an increase of the band gap under compression. Our synchrotron XRD measurements could not distinguish between Cmmm and C2/m BaO 2 definitively because the difference in XRD patterns is very subtle. However, our data do not show any sign of polymerization transition up to 120 GPa. Raman spectra of the O-O peroxide vibration show a small anomaly in frequency at 110 GPa, which is qualitatively similar to that predicted theoretically due to the Cmmm to C2/m phase transition thus supporting the predicted transformation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Topology invisible to eigenvalues in obstructed atomic insulators

Here we consider the extent to which symmetry eigenvalues reveal the topological character of bands. Specifically, we compare distinct atomic limit phases (band representations) that share the same irreducible representations (irreps) at all points in the Brillouin zone and, therefore, appear equivalent in a classification based on eigenvalues. We derive examples where such “irrep-equivalent” phases can be distinguished by a quantized Berry phase or generalization thereof. These examples constitute a generalization of the Su-Schrieffer-Heeger chain: neither phase is topological, in the sense that localized Wannier functions exist, yet there is a topological obstruction between them. We refer to two phases as “Berry obstructed atomic limits” if they have the same irreps, but differ by Berry phases. This is a distinct notion from eigenvalue obstructed atomic limits, which differ in their symmetry irreps at some point in the Brillouin zone. We compute exhaustive lists of elementary band representations that are irrep equivalent, in all space groups, with and without time-reversal symmetry and spin-orbit coupling, and use group theory to derive a set of necessary conditions for irrep equivalence. Finally, we conjecture, and in some cases prove, that irrep-equivalent elementary band representations that are not equivalent can be distinguished by a topological invariant.

36 MATERIALS SCIENCE↗

Long-distance nuclear matrix elements for neutrinoless double-beta decay from lattice QCD

Neutrinoless double-beta ( 0 ν β β ) decay is a heretofore unobserved process which, if observed, would imply that neutrinos are Majorana particles. Interpretations of the stringent experimental constraints on 0 ν β β -decay half-lives require calculations of nuclear matrix elements. This work presents the first lattice quantum chromodynamics (LQCD) calculation of the matrix element for 0 ν β β decay in a multinucleon system, specifically the n n → p p e e transition, mediated by a light left-handed Majorana neutrino propagating over nuclear-scale distances. This calculation is performed with quark masses corresponding to a pion mass of m π = 806 MeV at a single lattice spacing and volume. The statistically cleaner Σ − → Σ + e e transition is also computed in order to investigate various systematic uncertainties. The prospects for matching the results of LQCD calculations onto a nuclear effective field theory to determine a leading-order low-energy constant relevant for 0 ν β β decay with a light Majorana neutrino are investigated. This work, therefore, sets the stage for future calculations at physical values of the quark masses that, combined with effective field theory and nuclear many-body studies, will provide controlled theoretical inputs to experimental searches of 0 ν β β decay. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quantum Complexity of the Kronecker Coefficients

Whether or not the Kronecker coefficients of the symmetric group count some set of combinatorial objects is a longstanding open question. In this work we show that a given Kronecker coefficient is proportional to the rank of a projector that can be measured efficiently using a quantum computer. In other words a Kronecker coefficient counts the dimension of the vector space spanned by the accepting witnesses of a QMA verifier, where QMA is the quantum analogue of NP . This implies that approximating the Kronecker coefficients to within a given relative error is not harder than a certain natural class of problems that captures the complexity of estimating thermal properties of quantum many-body systems. A second consequence is that deciding positivity of Kronecker coefficients is contained in QMA , complementing a recent NP -hardness result of Ikenmeyer, Mulmuley, and Walter. We obtain similar results for the related problem of approximating row sums of the character table of the symmetric group. Finally, we discuss an efficient quantum algorithm that approximates normalized Kronecker coefficients to inverse-polynomial additive error. Published by the American Physical Society 2024

Physics↗

Advanced Models and Fast Simulation Algorithms for MARS Plants

With the successful commercialization of modular multilevel converters (MMC) and rapid increase in photovoltaic (PV) penetration and energy storage systems (ESS), one of the research areas being explored dynamically is integration of PV systems and ESS systems to MMCs. On this front, an integrated development of PV plants, ESS, and alternating- and direct-current (ac/dc) systems through Multi-Port Autonomous Reconfigurable Solar (MARS) plant is proposed and being researched upon. However, simulation of MARS systems on traditional existing software takes very long time (order of days) due to the presence of a large number of states and non-linear devices. Therefore, development of fast high fidelity models is of paramount importance to test the hardware design, control design and planning of power system expansions. In this context, two research objectives are addressed in this paper: firstly, an ultra fast single- or multi- CPU simulation algorithm to simulate the MARS system based on state-space models, hybrid discretization algorithm with a relaxation technique that reduces the imposed computational burden and a multi-rate method is presented. Secondly, a reduced order model of MARS is also presented. Both the developed techniques are validated with respect to reference PSCAD/EMTDC model.

Marthi, Phani Ratna Vanamali↗

PAGANI: a parallel adaptive GPU algorithm for numerical integration

We present a new adaptive parallel algorithm for the challenging problem of multi-dimensional numerical integration on massively parallel architectures. Adaptive algorithms have demonstrated the best performance, but efficient many-core utilization is difficult to achieve because the adaptive work-load can vary greatly across the integration space and is impossible to predict a priori. Existing parallel algorithms utilize sequential computations on independent processors, which results in bottlenecks due to the need for data redistribution and processor synchronization. Our algorithm employs a high-throughput approach in which all existing sub-regions are processed and sub-divided in parallel. Repeated sub-region classification and filtering improves upon a brute-force approach and allows the algorithm to make efficient use of computation and memory resources. A CUDA implementation shows orders of magnitude speedup over the fastest open-source CPU method and extends the achievable accuracy for difficult integrands. Our algorithm typically outperforms other existing deterministic parallel methods.

Sakiotis, Ioannis↗

Parallel Simulation of Beam Dynamics in Particle Accelerators [Slides]

Particle accelerators are among the most versatile and important tools of scientific discovery. The Nation's accelerators are responsible for a wealth of advances in materials science, chemistry, the biosciences, particle physics, and nuclear physics. They also have important applications to national security, the environment, energy, medicine, and on the quality of people's lives. LANL has a long history of making pioneering contributions to Accelerator Science including key contributions to the field of Computational Accelerator Physics. These include the development of early beam dynamics codes with space charge (such as PARMILA and PARMELA), the development of rf cavity codes and magnet codes (including Poisson and Superfish), and the development and distribution of codes to the accelerator community through the Los Alamos Accelerator Code Group. LANL researchers also helped pioneer the development of massively parallel space-charge codes. In project t22_accelsim we have moved beyond electrostatic models of collective effects (i.e., solving the Poisson equation in the bunch frame) to fully electromagnetic models based on the Lienard-Wiechert formalism. This approach enables the large-scale simulation of radiation production and collective effects in high brightness electron beams. This is highly relevant to LANL given its future goal of developing an X-ray Free Electron Laser (XFEL). It also directly impacts a LANL LDRD project to develop an undulator-based non-invasive beam profile monitor for beams created in laser-plasma accelerator systems.

43 PARTICLE ACCELERATORS↗

Nucleon valence quark distribution functions from Lattice QCD

We present results on the nucleon valence quark distribution extracted from Lattice QCD simulations, using a gauge ensemble of Nf=2+1 Wilson-Clover fermions with a pion mass of m?=350 MeV and lattice spacing of a=0.091 fm. We obtain reduced Ioffe Time Distributions (rITDs) by computing appropriate matrix elements on the lattice, and elaborate on the extraction of the desired quark distributions from the rITDs following the pseudo-PDF approach. A set of techniques are considered in order to ensure ground state dominance. Theoretical and experimental implications of our calculation are discussed.

Kallidonis, Christos↗

Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM SciDAC) (Technical Final Report)

Runaway electrons can severely damage the plasma facing components on ITER during a major disruption and pose a major risk for tokamak fusion. It has been recognized that an adequate disruption mitigation system (DMS) is essential for the safe operation of ITER. The United States is responsible for the design and implementation of the disruption mitigation system on ITER, and in July 2016 the Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM) was launched by DOE, in a joint Fusion Energy Sciences (FES) and Advanced Scientific Computing Research (ASCR) collaboration. SCREAM was a comprehensive theory and simulation SciDAC center that provided physics guidance in the avoidance and mitigation of runaway electrons, and in tandem with domestic and international experiments, helped establish the qualitative and quantitative bases for safe operational scenarios and viable mitigation techniques. The SCREAM center assembled a national team of experts in runaway electron physics, tokamak disruptions, magnetohydrodynamic (MHD) simulation, and advanced algorithms and computing. The team combined advanced simulation and analysis capability facilitated by direct participation of ASCR SciDAC institutes with theoretical models and code development by FES scientists to focus on the runaway risk for ITER and tokamaks in general. The research scope was focussed on integrated simulations of kinetic runaway electrons, including MHD and fluid models of impurity transport, within a research plan guided by theory. The specific research tasks were (1) establish the fundamental physics of runaway generation, saturation, and dynamical evolution in a tokamak; (2) examine the critical path toward runaway avoidance; and (3) investigate the viability and effectiveness of the leading candidate schemes for runaway mitigation. In all three areas, members of the team carried out scoping studies that established the readiness for rapid and critical advances, especially in the deployment and further development of large-to extreme-scale simulation tools. Our multi-pronged computational approach included (1) relativistic Fokker-Planck solvers with discretization in phase space, (2) self-consistent particle-in-cell techniques, (3) particle-based Monte-Carlo, and (4) MHD-particle hybrid simulations. Cross-check between these different methods provided an additional means for verification and further bolstered the fidelity of our physics prediction. Validation against experimental results brings confidence to the predictive capability for ITER and frequently leads to new ideas for understanding and mitigating the thermal quench driven runaway electron phenomenon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Qubit Regularization of Quantum Field Theories

To study quantum field theories on a quantum computer, we must begin with Hamiltonians defined on a finite-dimensional Hilbert space and then take appropriate limits. This approach can be seen as a new type of regularization for quantum field theories, which we refer to as qubit regularization. A related finite-dimensional regularization, known as the D-theory approach, was proposed long ago as a general framework for all quantum field theories. In this framework, the dimensionality of the local Hilbert space at each spatial point can increase as needed through an additional flavor index. To reproduce asymptotically free QFTs, most studies assume that qubit-regularized theories require extending the local Hilbert space to infinity. However, contrary to this common belief, recent discoveries in (1+1) dimensions have revealed two examples where asymptotic freedom appears to emerge within a strictly finite-dimensional local Hilbert space through a novel renormalization group (RG) flow. These findings motivate further investigation into whether asymptotically free gauge theories could also emerge within a strictly finite-dimensional local Hilbert space. To support these explorations, we propose an orthonormal basis called the monomer-dimer-tensor-network (MDTN) basis and use it to construct new types of qubit-regularized lattice gauge theories.

Chandrasekharan, Shailesh [Duke Univ., Durham, NC ↗

The Energetics of the Central Engine in the Powerful Quasar 3C 298

The compact steep-spectrum radio source 3C 298 (redshift of 1.44) has the largest 178 MHz luminosity in the Third Cambridge Revised Catalogue (3CR); its radio lobes are among the most luminous in the universe. The plasma state of the radio lobes is modeled with the aid of interferometric radio observations (in particular, the new Low Frequency Array observation and archival MERLIN data) and archival single-station data. It is estimated that the long-term time-averaged jet power required to fill these lobes with leptonic plasma is $\overline{Q}\,\approx 1.28\pm 0.51\times {10}^{47}\,\mathrm{erg}\,{{\rm{s}}}^{-1}$, rivaling the largest time-averaged jet powers from any quasar. Supporting this notion of extraordinary jet power is a 0.5–10 keV luminosity of ≈ 5.2 × 10 46 erg s –1 , comparable to luminous blazars, yet there is no other indication of strong relativistic beaming. We combine two new high signal-to-noise ratio optical spectroscopic observations from the Hobby-Eberly Telescope with archival Hubble Space Telescope, Two Micron All Sky Survey, and Galaxy Evolutionary Explorer data to compute a bolometric luminosity from the accretion flow of L bol ≈ 1.55 ± 0.15 × 10 47 erg s –1 . The ratio, $\overline{Q}/{L}_{\mathrm{bol}}\approx 1$, is the approximate upper limit for quasars. Characteristic of a large $\overline{Q}/{L}_{\mathrm{bol}}$, we find an extreme-ultraviolet (EUV) spectrum that is very steep (the "EUV deficit" of powerful radio quasars relative to radio-quiet quasars), and this weak ionizing continuum is likely a contributing factor to the relatively small equivalent widths of the broad emission lines in this quasar.

79 ASTRONOMY AND ASTROPHYSICS↗

Physical Yukawa couplings in heterotic string compactifications

One of the challenges of heterotic compactification on a CalabiYau threefold is to determine the physical (27) 3 Yukawa couplings of the resulting four-dimensional $\mathcal{N}$ = 1 theory. In general, the calculation necessitates knowledge of the Ricci-flat metric. However, in the standard embedding, which references the tangent bundle, we can compute normalized Yukawa couplings from the Weil-Petersson metric on the moduli space of complex structure deformations of the Calabi-Yau manifold. In various examples (the Fermat quintic, the intersection of two cubics in $\mathbb{P}$ 5 , and the TianYau manifold), we calculate the normalized Yukawa couplings for (2,1)-forms using the Weil-Petersson metric obtained from the Kodaira-Spencer map. In cases where $h^{1,1}$ = 1 , this is compared to a complementary calculation based on performing period integrals. A third expression for the normalized Yukawa couplings is obtained from a machine learned approximate Ricci-flat metric making use of explicit harmonic representatives. Finally, the excellent agreement between the different approaches opens the door to precision string phenomenology.

Butbaia, Giorgi↗