Search NASA⌕ Search

SEARCH · Search NASA

Results for “phase space methods”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Phase-space methods for neutrino oscillations: Extension to multibeams

The phase-space approach (PSA), which was originally introduced in Lacroix [] to describe neutrino flavor oscillations for interacting neutrinos emitted from stellar objects is extended to describe arbitrary numbers of neutrino beams. The PSA is based on mapping the quantum fluctuations into a statistical treatment by sampling initial conditions followed by independent mean-field evolution. A new method is proposed to perform this sampling that allows treating an arbitrary number of neutrinos in each neutrino beams. We validate the technique successfully and confirm its predictive power on several examples where a reference exact calculation is possible. We show that it can describe many-body effects, such as entanglement and dissipation induced by the interaction between neutrinos. Due to the complexity of the problem, exact solutions can only be calculated for rather limited cases, with a limited number of beams and/or neutrinos in each beam. The PSA approach considerably reduces the numerical cost and provides an efficient technique to accurately simulate arbitrary numbers of beams. Examples of PSA results are given here, including up to 200 beams with time-independent or time-dependent Hamiltonians. We anticipate that this approach will be useful to bridge exact microscopic techniques with more traditional transport theories used in neutrino oscillations. It will also provide important reference calculations for future quantum computer applications where other techniques are not applicable to classical computers. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Demonstration of a Novel Phase-Space Painting Method in a Coupled Lattice to Mitigate Space Charge in High-Intensity Hadron Beams

Multiturn charge-exchange injection is the primary method of creating high-intensity hadron beams in circular accelerators, and phase space painting during injection enables tailoring of the accumulated phase space distribution. A technique we call eigenpainting allows injection of particles into a single mode of a coupled ring, providing full four-dimensional control of the phase space distribution. Under ideal conditions, uniform eigenpainting generates a linear-force equilibrium distribution in the transverse plane, with zero volume in four-dimensional transverse phase space, even including space charge. Here, we have implemented eigenpainting for the first time in the spallation neutron source accumulator ring. Injecting 8.8 μ⁢C of an 800 MeV beam, we obtain a final ratio of intrinsic transverse emittances of ≈2.4. We analyze the effect of space charge on the final distribution through comparison of the reconstructed phase space to particle-in-cell simulations.

Evans, Nicholas J. [Oak Ridge National Laboratory ↗
DOE OSTI2022

Data-Driven Operator Theoretic Methods for Phase Space Learning and Analysis

This paper uses data-driven operator theoretic approaches to explore the global phase space of a dynamical system. In this work, we defined conditions for discovering new invariant subspaces in the state space of a dynamical system starting from an invariant subspace based on the spectral properties of the Koopman operator. When the system evolution is known locally in several invariant subspaces in the state space of a dynamical system, a phase space stitching result is derived that yields the global Koopman operator. Additionally, in the case of equivariant systems, a phase space stitching result is developed to identify the global Koopman operator using the symmetry properties between the invariant subspaces of the dynamical system and time-series data from any one of the invariant subspaces. Finally, these results are extended to topologically conjugate dynamical systems; in particular, the relation between the Koopman tuple of topologically conjugate systems is established. The proposed results are demonstrated on several second-order nonlinear dynamical systems including a bistable toggle switch. Our method elucidates a strategy for designing discovery experiments: experiment execution can be done in many steps, and models from different invariant subspaces can be combined to approximate the global Koopman operator.

42 ENGINEERING↗
DOE OSTI2022

Role of non-Gaussian quantum fluctuations in neutrino entanglement

The flavor evolution of neutrinos in environments with large neutrino number densities is an open problem at the nexus of astrophysics and neutrino flavor physics. Among the many unanswered questions pertaining to this problem, it remains to be determined whether neutrino-neutrino coherent scattering can give rise to nontrivial quantum entanglement among neutrinos, and whether this can affect the flavor evolution in a meaningful way. To gain further insight into this question, here we study a simple system of two interacting neutrino beams and obtain the exact phase space explored by this system using the Husimi quasiprobability distribution. We observe that the entanglement induced by the coupling leads to strong delocalization in phase-space with largely non-Gaussian quantum fluctuations. The link between the neutrino entanglement and quantum fluctuations is illustrated using the one- and two-neutrino entropy. In addition, we propose an approximate phase-space method to describe the interacting neutrinos problem, where the exact evolution is replaced by a set of independent mean-field evolutions with a statistical sampling of the initial conditions. The phase-space approach provides a simple and accurate method to describe the gross features of the neutrino entanglement problem. Applications are shown using time-independent and time-dependent Hamiltonians in the nonadiabatic regime.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Role of non-gaussian quantum fluctuations in neutrino entanglement

The flavor evolution of neutrinos in environments with large neutrino number densities is an open problem at the nexus of astrophysics and neutrino flavor physics. Among the many unanswered questions pertaining to this problem, it remains to be determined whether neutrino-neutrino coherent scattering can give rise to nontrivial quantum entanglement among neutrinos, and whether this can affect the flavor evolution in a meaningful way. To gain further insight into this question, here we study a simple system of two interacting neutrino beams, and obtain the exact phase-space explored by this system using the Husimi quasi-probability distribution. We observe that the entanglement induced by the coupling leads to strong delocalization in phase-space with largely non-Gaussian quantum fluctuations. The link between the neutrino entanglement and quantum fluctuations is illustrated using the one- and two-neutrino entropy. In addition, we propose an approximate phase-space method to describe the interacting neutrinos problem, where the exact evolution is replaced by a set of independent mean-field evolutions with a statistical sampling of the initial conditions. The phase-space approach provides a simple and accurate method to describe the gross features of the neutrino entanglement problem. Applications are shown using time-independent and time-dependent Hamiltonians in the non-adiabatic regime.

Patwardhan, Amol↗
DOE OSTI2023

An interoperable implementation of collective‐variable based enhanced sampling methods in extended phase space within the OpenMM package

Collective variable (CV)-based enhanced sampling techniques are widely used today for accelerating barrier-crossing events in molecular simulations. A class of these methods, which includes temperature accelerated molecular dynamics (TAMD)/driven-adiabatic free energy dynamics (d-AFED), unified free energy dynamics (UFED), and temperature accelerated sliced sampling (TASS), uses an extended variable formalism to achieve quick exploration of conformational space. These techniques are powerful, as they enhance the sampling of a large number of CVs simultaneously compared to other techniques. Extended variables are kept at a much higher temperature than the physical temperature by ensuring adiabatic separation between the extended and physical subsystems and employing rigorous thermostatting. Here, in this work, we present a computational platform to perform extended phase space enhanced sampling simulations using the open-source molecular dynamics engine OpenMM. The implementation allows users to have interoperability of sampling techniques, as well as employ state-of-the-art thermostats and multiple time-stepping. This work also presents protocols for determining the critical parameters and procedures for reconstructing high-dimensional free energy surfaces. As a demonstration, we present simulation results on the high dimensional conformational landscapes of the alanine tripeptide in vacuo, tetra-N-methylglycine (tetra-sarcosine) peptoid in implicit solvent, and the Trp-cage mini protein in explicit water.

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

Energetic Particles - microturbulence interaction thrust (EP SIWG white paper)

This white paper outlines a strategic approach to Energetic Particle (EP) research, aimed at applying the results to the future burning devices such as ITER. Our goal is to deepen the understanding of the interplay between plasma background microturbulence and diverse species of energetic ions in burning plasma devices. This strategy is poised to advance our knowledge in this crucial area of fusion science, paving the way for significant breakthroughs in plasma physics and fusion energy research. We consider the confinement of energetic ions such as auxiliary heating beams, minority ions from ICRH, and fusion product alpha particles. Recent studies have underscored the potential of EP/microturbulence interaction to enhance fusion plasma performance by stabilizing microturbulence. However, it also poses risks, such as exacerbating Alfvénic eigenmode (AE) instabilities, which could compromise the sustainability of plasma discharge by causing EP losses. The interplay between EPs and background microturbulence is evident in the phenomenon of effective pitch angle scattering, a crucial aspect of the quasilinear (QL) theory. This theory is integral to developing numerically-efficient yet comprehensive and self-consistent approaches, recently employed to investigate the relaxation of energetic particle populations in the holistic modeling of fusion-grade plasmas. The primary objectives of this white paper thrust encompass a dual focus. First is that we need to delve into the mechanisms through which microturbulence engenders effective pitch angle scattering. This investigation will entail the examination of the presence of known or self-consistently predicted spectra of modes accountable for microturbulence. In these explorations, electrostatic microturbulence serves as a logical initial stride towards achieving the trust's objectives. This endeavor is anticipated to yield formulations expressing the parametric dependencies of the effective pitch angle scattering frequency on variables such as thermal plasma electron and/or ion temperatures, as well as their respective thermal conductivities. In our second objective, we aim to delve into the intricate formation of zonal flow (ZF) structures amidst the complex interplay of microturbulence and Alfvenic eigenmodes (AEs). This endeavor poses greater challenges as we seek to unravel the macroscopic manifestations influenced by microturbulence, termed as zonal structures (ZS), stemming from microturbulence-induced ZF. Drawing from this understanding, we anticipate employing a QL approach to yield comprehensive insights into the distribution function of energetic particles (EP) within phase space. This method entails resolving the dominant multidimensional phase space diffusion processes while effectively averaging over the rapid ballistic responses. Nevertheless, substantial strides remain imperative to realize a comprehensive whole-device modeling framework. This entails meticulous verification and validation exercises against experimental observations, as well as rigorous benchmarking against theoretical frameworks and numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ∼1/Np with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly non-equilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ~1/$\sqrt{N_p}$ with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly nonequilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Physics-constrained superresolution diffusion for six-dimensional phase space diagnostics

Adaptive physics-constrained superresolution diffusion is developed for noninvasive virtual diagnostics of the six-dimensional (6D) phase space density of charged particle beams. An adaptive variational autoencoder embeds initial beam condition images and scalar measurements to a low-dimensional latent space from which a 32 6 pixel 6D tensor representation of the beam's 6D phase space density is generated. Projecting from a 6D tensor generates physically consistent two-dimensional projections. Physics-guided superresolution diffusion transforms low-resolution images of the 6D density to high resolution 256 × 256 pixel images. Unsupervised adaptive latent space tuning enables tracking of time-varying beams without knowledge of time-varying initial conditions. The method is demonstrated with experimental data and multiparticle simulations at the HiRES UED. The general approach is applicable to a wide range of complex dynamic systems evolving in high-dimensional phase space. The method is shown to be robust to distribution shift without retraining. Published by the American Physical Society 2025

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Predicting beam transmission using 2-dimensional phase space projections of hadron accelerators

We present a method to compress the 2D transverse phase space projections from a hadron accelerator and use that information to predict the beam transmission. This method assumes that obtaining at least three projections of the 4D transverse phase space is possible and that an accurate simulation model is available for the beamline. Using a simulated model, we show that—a computer can train a convolutional autoencoder to reduce phase-space information which can later be used to predict the beam transmission. Finally, we argue that although using projections from a realistic nonlinear distribution produces less accurate results, the method still generalizes well.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Accelerating particle-in-cell kinetic plasma simulations via reduced-order modeling of space-charge dynamics using dynamic mode decomposition

We present a data-driven reduced-order modeling of the space-charge dynamics for electromagnetic particle-in-cell (EMPIC) plasma simulations based on dynamic mode decomposition (DMD). The dynamics of the charged particles in kinetic plasma simulations such as EMPIC is manifested through the plasma current density defined along the edges of the spatial mesh. We showcase the efficacy of DMD in modeling the time evolution of current density through a low-dimensional feature space. Not only do such DMD based predictive reduced-order models help accelerate EMPIC simulations, they also have the potential to facilitate investigative analysis and control applications. Here, we demonstrate the proposed DMD-EMPIC scheme for reduced-order modeling of current density and speedup in EMPIC simulations involving electron beam under the influence of magnetic field, virtual cathode oscillations, and backward wave oscillator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Realistic simulations of spin squeezing and cooperative coupling effects in large ensembles of interacting two-level systems

We describe an efficient numerical method for simulating the dynamics of interacting spin ensembles in the presence of dephasing and decay. The method builds on the discrete truncated Wigner approximation for isolated systems, which combines the mean-field dynamics of a spin ensemble with a Monte Carlo sampling of discrete initial spin values to account for quantum correlations. Here we show how this approach can be generalized for dissipative spin systems by replacing the deterministic mean-field evolution by a stochastic process, which describes the decay of coherences and populations while preserving the length of each spin. We demonstrate the application of this technique for simulating nonclassical spin-squeezing effects or the dynamics and steady states of cavity QED models with 10 5 interacting two-level systems. This opens up the possibility to perform accurate real-scale simulations of a diverse range of experiments in quantum optics or with solid-state spin ensembles under realistic laboratory conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Four-dimensional emittance measurement at the Spallation Neutron Source

A coasting hadron beam with an elliptical transverse profile, uniform charge density, and small transverse four-dimensional (4D) emittance could improve accelerator performance in several contexts. A phase space painting method to generate such a distribution is being tested in the Spallation Neutron Source (SNS) accumulator ring. A critical component of these efforts is to measure the 4D emittance of the beam during accumulation. The 4D emittance can be reconstructed from measured beam profiles in two ways: in the multi-optics method, the optics between a reconstruction and measurement location are varied; in the fixed-optics method, multiple measurement locations are used without modifying the optics. The fixed-optics method is faster but can lead to large uncertainty in the reconstructed 4D emittance. In this paper, we implement a variant of the multi-optics method using the four available wire-scanners near the SNS target. We also modify the wire-scanner region to reduce the uncertainty of the fixed-optics method. We then demonstrate the usefulness of the fixed-optics method by reconstructing the 4D emittance evolution during accumulation in the SNS ring.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Poisson_CCD: A dedicated simulator for modeling CCDs

A dedicated simulator, Poisson_CCD, has been constructed which models astronomical CCDs by solving Poisson’s equation numerically and simulating charge transport within the CCD. The potentials and free carrier densities within the CCD are self-consistently solved for, giving realistic results for the charge distribution within the CCD storage wells. The simulator has been used to model the CCDs that are being used to construct the LSST digital camera. The simulator output has been validated by comparing its predictions with several different types of CCD measurements, including astrometric shifts, brighter-fatter induced pixel–pixel covariances, saturation effects, and diffusion spreading. The code is open source and freely available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

On the merit of hot ion mode for tearing mode stabilization

The stabilization of tearing modes with rf driven current benefits from the cooperative feedback loop between rf power deposition and electron temperature within the island. This effect, termed rf current condensation, can greatly enhance and localize the current driven within magnetic islands. It has previously been shown that the condensation effect opens the possibility of passive stabilization with broad rf profiles, as would be typical of LHCD for steady state operation. Here, we show that this self-healing effect can be dramatically amplified by operation in a hot ion mode, due to the additional electron heat source provided by the hotter ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Metaplectic geometrical optics for ray-based modeling of caustics: Theory and algorithms

The optimization of radio frequency-wave (RF) systems for fusion experiments is often performed using ray-tracing codes, which rely on the geometrical-optics (GO) approximation. However, GO fails at caustics such as cutoffs and focal points, erroneously predicting the wave intensity to be infinite. This is a critical shortcoming of GO, since the caustic wave intensity is often the quantity of interest, e.g., RF heating. Full-wave modeling can be used instead, but the computational cost limits the speed at which such optimizations can be performed. Here, we have developed a less expensive alternative called metaplectic geometrical optics (MGO). Instead of evolving waves in the usual x (coordinate) or k (spectral) representation, MGO uses a mixed X$\equiv$Ax+Bk representation. By continuously adjusting the matrix coefficients A and B along the rays, one can ensure that GO remains valid in the X coordinates without caustic singularities. The caustic-free result is then mapped back onto the original x space using metaplectic transforms. Here, we overview the MGO theory and review algorithms that will aid the development of an MGO-based ray-tracing code. We show how using orthosymplectic transformations leads to considerable simplifications compared to previously published MGO formulas. We also prove explicitly that MGO exactly reproduces standard GO when evaluated far from caustics (an important property that until now has only been inferred from numerical simulations), and we relate MGO to other semiclassical caustic-removal schemes published in the literature. Finally this discussion is then augmented by an explicit comparison of the computed spectrum for a wave bounded between two cutoffs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Particle-in-cell modeling of plasma jet merging in the large-Hall-parameter regime

The merging process of magnetized plasma jets with parameters relevant to the plasma-jet-driven magneto-inertial fusion (PJMIF) design and the plasma liner experiment (PLX) is modeled by fully kinetic particle-in-cell (PIC) simulations in one and two spatial dimensions. Here, the modified two-stream instability is identified to be the main mechanism responsible for stopping the plasma jets and preventing species interpenetration. The electron and ion Hall parameters of the merged plasma are greater than unity, and the plasma β is close to unity, which is the desired characteristic of planned experiments at PLX. Our 2D PIC simulations validate the results of the radiation magneto-hydrodynamics code FLASH, which will be the primary tool for modeling various stages of future PJMIF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗