SEARCH · Search NASA

Results for “ANALOG SIMULATION”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Commissioning a time-gated camera for fast neutron beamline spatial-energy characterization at LANSCE-WNR spallation source

An energy-resolved fast neutron beam imaging diagnostic has been successfully commissioned at the Weapons Neutron Research (WNR) spallation source within the Los Alamos Neutron Science Center (LANSCE) facility. This diagnostic replaces the existing analog phosphor image plates, which integrate across all neutron energies, as well as other particles, with a near-real-time energy-sensitive imaging capability. The system uses a fast plastic scintillator coupled with an intensified CCD camera. Specifically, the Teledyne Pi-MAX4 camera is coupled with either a 4 mm thick Eljen (EJ) 204 or 228 plastic scintillator. These scintillators are most sensitive to the fast neutrons (0.8-800 MeV) directly from the spallation source rather than low energy background radiation. Experimentally, these plastic scintillators were shown to have sufficiently fast decay to differentiate the bright gamma flash from the spallation neutrons. The spatial resolution is dominated by neutron beam divergence, with minimal additional contributions from scatter and light divergence. The system successfully resolved changes in neutron beam characteristics caused by intentional proton steering variations. Additionally, simulations of scintillator light yield as a function of thickness conducted using PHITS (with Scinful-QMD package) found that increasing scintillator thickness from 4 mm to 6 or 8 mm could potentially increase brightness ~ 3x. This may be explored if there is a need to reduce image acquisition time from several minutes to under one minute.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2025

Cost of emulating a small quantum annealing problem in the circuit model

Demonstrations of quantum advantage for certain sampling problems have generated considerable excitement for quantum computing and have further spurred the development of circuit-model quantum computers, which represent quantum programs as a sequence of quantum gates acting on a finite number of qubits. Amongst this excitement, analog quantum computation has become less prominent, with the expectation that circuit-model quantum computers will eventually be sufficient for emulating analog quantum computation and thus rendering analog quantum computation obsolete. In this work we explore the basic requirements for emulating a specific analog quantum computation in the circuit model: the preparation of a biased superposition of degenerate ground states of an Ising Hamiltonian using an adiabatic evolution. We show that the overhead of emulation is substantial even for this simple problem. This supports using analog quantum computation for solving time-dependent Hamiltonian dynamics in the short term and midterm, assuming analog errors can be made low enough and coherence times long enough to solve problems of practical interest.

Quantum algorithms & computation↗

DOE OSTI2026

Single-Molecule Electron Transport in Peptoids

Peptoids are structural analogs of peptides in which side chains are appended to the backbone nitrogen rather than the α-carbon. The sequence-defined modularity of peptoids enables precise control over structure−function relationships, enabling applications in energy storage and biomedical materials. Despite recent progress, the role of sequence and conformation on electron transport in peptoid molecules is not fully understood. Here, we synthesize a library of peptoid oligomers and characterize their molecular electronic properties using the scanning tunneling microscope-break junction (STM-BJ) technique. Our results show well-defined electron transport behavior for peptoid sequences containing aromatic side groups lacking hydrogen bonds (H-bonds) and without chemical substitutions at the N−C α position. This behavior fundamentally differs from electron transport in peptides, where H-bond interactions give rise to higher conductance states. All-atom molecular dynamics (MD) simulations are used to understand the conformational heterogeneity of peptoids, and molecular conformations obtained from MD simulations are used in quantum mechanical calculations based on the nonequilibrium Green’s function−density functional theory (NEGF-DFT) formalism. In all cases, computational results are in reasonable qualitative agreement with experiments. Our work demonstrates that the conductance behavior of peptoids depends on monomer identity, including side-chain aromaticity and substitution at the N−C α position. Overall, this work provides new insights into the structure−function relationships governing electron transport in peptoid-based materials and establishes design rules for peptoid-based molecular junctions.

Charge transport↗

DOE OSTI2024

Tailoring coherent charge transport in graphene by deterministic defect generation

Harnessing the wave-nature of charge carriers in solid state devices, electron optics investigates and exploits coherent phenomena, in analogy with optics and photonics. Typically, this requires complex electronic devices leveraging macroscopically coherent charge transport in two-dimensional electron gases and superconductors. Here, collective coherent effects are induced in a simple counterintuitive architecture by defect engineering. Deterministically introduced lattice defects in graphene enable the phase coherent charge transport by playing the role of potential barriers, instead of scattering centres as conventionally considered. Thus, graphene preserves its quasi-ballistic quantum transport and can support phase-matched charge carrier waves. Based on this approach, multiple electronic Fabry-Pèrot cavities are formed by creating periodically alternating defective and pristine nano-stripes through low energy electron-beam irradiation. Indeed, defective stripes behave as partially reflecting mirrors and resonantly confine the charge carrier waves within the pristine areas, giving rise to Fabry-Pèrot resonant modes. These modes experimentally manifest as sheet resistance oscillations, as also confirmed by Landauer-Büttiker simulations. Moreover, these coherent phenomena survive up to 30 K for both polarities of charge carriers, contrarily to traditional monopolar electrostatically created Fabry-Pèrot interferometers. Our study positions defective graphene as an innovative platform for coherent electronic devices, with potential applications in nano and quantum technologies.

Melchioni, Nicola↗

DOE OSTI2025

On the Stability Analysis of Astrophysical Cooling Functions

To model the temperature evolution of optically thin astrophysical environments at MHD scales, radiative and collisional cooling rates are typically either pretabulated or fit into a functional form and then input into MHD codes as a radiative loss function. Thermal balance requires estimates of the analogous heating rates, which are harder to calculate, and due to uncertainties in the underlying dissipative heating processes these rates are often simply parameterized. The resulting net cooling function defines an equilibrium curve that varies with density and temperature. Such cooling functions can make the gas prone to thermal instability (TI), which will cause departures from equilibrium. There has been no systematic study of thermally unstable parameter space for nonequilibrium states. Motivated by our recent finding that there is a related linear instability, catastrophic cooling instability, that can dominate over TI, here we carry out such a study. We show that Balbus instability criteria for TI can be used to define a critical cooling rate, Λ c , that permits a nonequilibrium analysis of cooling functions through the mapping of TI zones. We furthermore illustrate how thermal conduction modifies the shape of TI zones. Upon applying a Λ c -based stability analysis to coronal loop simulations, we find that loops undergoing periodic episodes of coronal rain formation are linearly unstable to catastrophic cooling instability, while TI is stabilized by thermal conduction.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

SAGAbg. III. Environmental Stellar Mass Functions, Self-quenching, and the Stellar-to-halo Mass Relation in the Dwarf Galaxy Regime

Recent efforts have extended our view of the number and properties of satellite galaxies beyond the Local Group firmly down to M ⋆ ∼ 10 6 M ⊙ . A similarly complete view of the field dwarf population has lagged behind. Using the background galaxy sample from the Satellites Around Galactic Analogs (SAGA) survey at z < 0.05, we take inventory of the dwarf population down to M ⋆ ∼ 5 × 10 6 M ⊙ using three metrics: the stellar mass function (SMF) as a function of environment, the stellar-to-halo mass relation (SHMR) of dwarf galaxies inferred via abundance matching, and the quenched fraction of highly isolated dwarfs. We find that the low-mass SMF shape shows minimal environmental dependence, with the field dwarf SMF described by a low-mass power-law index of α 1 = −1.44 ± 0.09 down to M ⋆ ∼ 5 × 10 6 M ⊙ , and that the quenched fraction of isolated dwarfs drops monotonically to f q ∼ 10 −3 at M ⋆ ∼ 10 8.5 M ⊙ . Though slightly steeper than estimates from H I kinematic measures, our inferred SHMR agrees with literature measurements of satellite systems, consistent with minimal environmental dependence of the SHMR in the probed mass range. Finally, although most contemporary cosmological simulations against which we compare accurately predict the SAGAbg-SMF SHMR, we find that big-box cosmological simulations largely overpredict isolated galaxy quenched fractions via a turnaround in f q (M ⋆ ) at 10 8 ≲ M ⋆ /M ⊙ ≲ 10 9 , underscoring the complexities in disentangling the drivers of galaxy formation and the need for systematic multidimensional observations of the dwarf population across environments.

Kado-Fong, Erin [Yale Univ., New Haven, CT (United↗

DOE OSTI2026

Graphical User Interface for a Biasing Board for the PSEC6

The PSEC6 is an application-specific integrated circuit (ASIC) designed for a readout system for a large area picosecond photodetector (LAPPD). The PSEC6 is currently in fabrication and pending testing. The testing system for the PSEC5, the previous iteration of the ASIC, required expensive and non-portable equipment, because the ASIC needs twelve adjustable reference voltages. The new testing system consists of an low-cost, open-source, cross-platform graphical user interface (GUI), a digital system, and a biasing board. The digital system is the interface between the GUI and biasing board, and can be implemented on a microcontroller or field-programmable gate array (FPGA). The biasing board contains twelve digital-to-analog converters (DACs) that are configurable via the GUI, which gives users the ability to write voltage values to all or specific DACs. The GUI was developed in C on Linux using the widget library GTK4 and cross-compiled for Windows compatibility. I2C and SPI protocols were implemented on an Adafruit Feather ESP32-S3 microcontroller to write commands to the DACs and PSEC6. A hardware implementation of the I2C protocol is in development on an FPGA. Since LAPPDs will be used by the Accelerator Neutrino Neutron Interaction Experiment (ANNIE) at Fermilab, the PSEC6 testing system in this internship project can potentially benefit future neutrino research. The project is relevant to the Department of Energy’s microelectronics mission, because the PSEC6 is an ASIC that will handle fast time signals arriving from the detector for readout. It also provided experience with building a cross-platform user interface, practicing digital design and implementation in hardware description language (HDL), and using simulations to inform new design iterations.

Guerrero, Sasha Camila [North Central Coll.]↗

DOE OSTI2026

Graphical User Interface for a Biasing Board for the PSEC6

Guerrero, Sasha Camila [North Central Coll.]↗

DOE OSTI2025

Engineering 2D Square Lattice Hubbard Models in 90° Twisted GeX/SnX (X =S, Se) Moiré Superlattices

Because of the large-period superlattices emerging in moiré two-dimensional (2D) materials, electronic states in such systems exhibit low energy flat bands that can be used to simulate strongly correlated physics in a highly tunable setup. While many investigations have thus far focused on moiré flat bands and emergent correlated electron physics in triangular, honeycomb, and quasi-one-dimensional lattices, tunable moiré realizations of square lattices subject to strong correlations remain elusive. Here, in this work, we propose a feasible scheme to construct moiré square lattice systems by twisting two or more layers of 2D materials in a rectangular lattice by 90°. We demonstrate the concept with twisted GeX/SnX (X =S, Se) moiré superlattices and calculate their electronic structures from first principles. We show that the lowest conduction flat band in these systems can be described by a square lattice Hubbard model with parameters which can be controlled by varying the choice of host materials, number of layers, and external electric fields. In particular, twisted double bilayer GeSe realizes a square lattice Hubbard model with strong frustration due to the next-nearest-neighbor hopping that could host unconventional superconductivity, in close analogy to the Hubbard model for copper-oxygen planes of cuprate high-temperature superconductors. The presented scheme uses 90° twisted 2D materials with rectangular unit cells as a promising platform for realizing the physical phenomena of square lattice Hubbard models, establishing a new route for studying its rich phase diagram of magnetism, charge order, and unconventional superconductivity in a highly tunable setting.

2-dimensional systems↗

DOE OSTI2024

SAGAbg. I. A Near-unity Mass-loading Factor in Low-mass Galaxies via Their Low-redshift Evolution in Stellar Mass, Oxygen Abundance, and Star Formation Rate

Measuring the relation between star formation and galactic winds is observationally difficult. In this work we make an indirect measurement of the mass-loading factor (the ratio between the mass outflow rate and star formation rate) in low-mass galaxies using a differential approach to modeling the low-redshift evolution of the star-forming main sequence and mass–metallicity relation. We use Satellites Around Galactic Analogs (SAGA) background galaxies, i.e., spectra observed by the SAGA Survey that are not associated with the main SAGA host galaxies, to construct a sample of 11,925 spectroscopically confirmed low-mass galaxies from 0.01 ≲ z ≤ 0.21 and measure auroral line metallicities for 120 galaxies. The crux of the method is to use the lowest-redshift galaxies as the boundary condition of our model, and to infer a mass-loading factor for the sample by comparing the expected evolution of the low-redshift reference sample in stellar mass, gas-phase metallicity, and star formation rate against the observed properties of the sample at higher redshift. We infer a mass-loading factor of η m = $0.92^{+1.76}_{-0.74}$, which is in line with direct measurements of the mass-loading factor from the literature despite the drastically different sets of assumptions needed for each approach. While our estimate of the mass-loading factor is in good agreement with recent galaxy simulations that focus on resolving the dynamics of the interstellar medium, it is smaller by over an order of magnitude than the mass-loading factor produced by many contemporary cosmological simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Playing Nonlocal Games across a Topological Phase Transition on a Quantum Computer

Many-body quantum games provide a natural perspective on phases of matter in quantum hardware, crisply relating the quantum correlations inherent in phases of matter to the securing of quantum advantage at a device-oriented task. In this Letter, we introduce a family of multiplayer quantum games for which topologically ordered phases of matter are a resource yielding quantum advantage. Unlike previous examples, quantum advantage persists away from the exactly solvable point and is robust to arbitrary local perturbations, irrespective of system size. We demonstrate this robustness experimentally on Quantinuum’s H1-1 quantum computer by playing the game with a continuous family of randomly deformed toric code states that can be created with constant-depth circuits leveraging midcircuit measurements and unitary feedback. We are thus able to tune through a topological phase transition—witnessed by the loss of robust quantum advantage—on currently available quantum hardware. This behavior is contrasted with an analogous family of deformed Greenberger-Horne-Zeilinger states, for which arbitrarily weak local perturbations destroy quantum advantage in the thermodynamic limit. Lastly, we discuss a topological interpretation of the game, which leads to a natural generalization involving an arbitrary number of players.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Velocity reconstruction in the era of DESI and Rubin/LSST. I. Exploring spectroscopic, photometric, and hybrid samples

Peculiar velocities of galaxies and halos can be reconstructed from their spatial distribution alone. This technique is analogous to the baryon acoustic oscillations reconstruction, using the continuity equation to connect density and velocity fields. The resulting reconstructed velocities can be used to measure imprints of galaxy velocities on the cosmic microwave background like the kinematic Sunyaev-Zel’dovich effect or the moving lens effect. As the precision of these measurements increases, characterizing the performance of the velocity reconstruction becomes crucial to allow unbiased and statistically optimal inference. In this paper, we quantify the relevant performance metrics: the variance of the reconstructed velocities and their correlation coefficient with the true velocities. We show that the relevant velocities to reconstruct for kSZ and moving lens are actually the halo—rather than galaxy—velocities. We quantify the impact of redshift-space distortions, photometric redshift errors, satellite galaxy fraction, incorrect cosmological parameter assumptions and smoothing scale on the reconstruction performance. Here, we also investigate hybrid reconstruction methods, where velocities inferred from spectroscopic samples are evaluated at the positions of denser photometric samples. We find that using exclusively the photometric sample is better than performing a hybrid analysis. The 2 Gpc/ℎ length simulations from abacussummit with realistic galaxy samples for DESI and Rubin LSST allow us to perform this analysis in a controlled setting. In the companion paper [B. Hadzhiyska, S. Ferraro, B. Ried Guachalla, and E. Schaan, companion paper, Phys. Rev. D 109, 103534 (2024).], we further include the effects of evolution along the light cone and give realistic performance estimates for DESI luminous red galaxies, emission line galaxies, and Rubin LSST-like samples.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

SAGAbg. II. The Low-mass Star-forming Sequence Evolves Significantly between 0.05 < z < 0.21

The redshift-dependent relation between galaxy stellar mass and star formation rate (SFR), known as the star-forming sequence (SFS), is a key observational yardstick for galaxy assembly. We use the SAGAbg-A sample of background galaxies from the Satellites Around Galactic Analogs (SAGA) Survey to model the low-redshift evolution of the low-mass SFS. The sample is comprised of 23,258 galaxies with Hα-based SFRs spanning 6 < log 10 (M * /[M ⊙ ]) < 10 and z < 0.21 (t < 2.5 Gyr). Although it is common to bin or stack galaxies at z ≲ 0.2 for galaxy population studies, the difference in lookback time between z = 0 and z = 0.21 is comparable to the time between z = 1 and z = 2. We develop a model to account for both the physical evolution of low-mass SFS and the selection function of the SAGA Survey, allowing us to disentangle redshift evolution from redshift-dependent selection effects across the SAGAbg-A redshift range. Our findings indicate significant evolution in the SFS over the last ∼2.5 Gyr, with a rising normalization: $\langle$SRF(M * = 10 8.5 M ⊙ )$\rangle$ (z) = 1.24$^{+0.25}_{–0.23}$z – 1.47$^{+0.03}_{–0.03}$. We also identify the redshift limit at which a static SFS is ruled out at the 95% confidence level, which is z = 0.05 based on the precision of the SAGAbg-A sample. Comparison with cosmological hydrodynamic simulations reveals that some contemporary simulations underpredict the recent evolution of the low-mass SFS. This demonstrates that the recent evolution of the low-mass SFS can provide new constraints on the assembly of the low-mass Universe and highlights the need for improved models in this regime.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

On the minimum number of radiation field parameters to specify gas cooling and heating functions

Fast and accurate approximations of gas cooling and heating functions are needed for hydrodynamic galaxy simulations. We use machine learning to analyze atomic gas cooling and heating functions in the presence of a generalized incident local radiation field computed by Cloudy. We characterize the radiation field through binned radiation field intensities instead of the photoionization rates used in our previous work. We find a set of 6 energy bins whose intensities exhibit relatively low correlation. We use these bins as features to train machine learning models to predict Cloudy cooling and heating functions at fixed metallicity. We compare the relative SHapley Additive exPlanation (SHAP) value importance of the features. From the SHAP analysis, we identify a feature subset of 3 energy bins (0.5-1, 1-4, and 13-16Ry) with the largest importance and train additional models on this subset. We compare the mean squared errors and distribution of errors on both the entire training data table and a randomly selected 20% test set withheld from model training. The machine learning models trained with 3 and 6 bins, as well as 3 and 4 photoionization rates, have comparable accuracy everywhere, with errors ≳10 times smaller than for the interpolation table of Gnedin and Hollon (2012). We conclude that 3 energy bins (or 3 analogous photoionization rates: molecular hydrogen photodissociation, neutral hydrogen HI, and fully ionized carbon CVI) are sufficient to characterize the dependence of the gas cooling and heating functions on our assumed incident radiation field model.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

On the minimum number of radiation field parameters to specify gas cooling and heating functions

Fast and accurate approximations of gas cooling and heating functions are needed for hydrodynamic galaxy simulations. We use machine learning to analyze atomic gas cooling and heating functions computed by Cloudy in the presence of a generalized incident local radiation field. We characterize the radiation field through binned radiation field intensities instead of the photoionization rates used in our previous work. We find a set of 6 energy bins whose intensities exhibit relatively low correlation. We use these bins as features to train machine learning models to predict Cloudy cooling and heating functions at fixed metallicity. We compare the relative SHapley Additive exPlanation (SHAP) value importance of the features. From the SHAP analysis, we identify a feature subset of 3 energy bins ($0.5-1, 1-4$, and $13-16 \, \mathrm{Ry}$) with the largest importance and train additional models on this subset. We compare the mean squared errors and distribution of errors on both the entire training data table and a randomly selected 20% test set withheld from model training. The machine learning models trained with 3 and 6 bins, as well as 3 and 4 photoionization rates, have comparable accuracy everywhere, with errors $\gtrsim 10$ times smaller than for the interpolation table of Gnedin and Hollon (2012). We conclude that 3 energy bins (or 3 analogous photoionization rates: molecular hydrogen photodissociation, neutral hydrogen HI, and fully ionized carbon CVI) are sufficient to characterize the dependence of the gas cooling and heating functions on our assumed incident radiation field model.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Development of a Data Overflow Protection System for Super-Kamiokande to Maximize Data from Nearby Supernovae

Neutrinos from very nearby supernovae, such as Betelgeuse, are expected to generate more than ten million events over 10 s in Super-Kamokande (SK). At such large event rates, the buffers of the SK analog-to-digital conversion board (QBEE) will overflow, causing random loss of data that are critical for understanding the dynamics of the supernova explosion mechanism. In order to solve this problem, two new data-acquisition (DAQ) modules were developed to aid in the observation of very nearby supernovae. The first of these, the SN module, is designed to save only the number of hit photomultiplier tubes during a supernova burst and the second, the Veto module, prescales the high-rate neutrino events to prevent the QBEE from overflowing based on information from the SN module. In the event of a very nearby supernova, these modules allow SK to reconstruct the time evolution of the neutrino event rate from beginning to end using both QBEE and SN module data. This paper presents the development and testing of these modules together with an analysis of supernova-like data generated with a flashing laser diode. We demonstrate that the Veto module successfully prevents DAQ overflows for Betelgeuse-like supernovae as well as the long-term stability of the new modules. During normal running the Veto module is found to issue DAQ vetos a few times per month resulting in a total dead-time less than 1 ms, and does not influence ordinary operations. Additionally, using simulation data we find that supernovae closer than 800 pc will trigger the Veto module, resulting in a prescaling of the observed neutrino data.

F20 Instrumentation and technique↗

DOE OSTI2025

A parallel-kinetic-perpendicular-moment model for magnetised plasmas

We describe a new model for the study of weakly collisional, magnetised plasmas derived from exploiting the separation of the dynamics parallel and perpendicular to the magnetic field. This unique system of equations retains the particle dynamics parallel to the magnetic field while approximating the perpendicular dynamics through a spectral expansion in the perpendicular degrees of freedom, analogous to moment-based fluid approaches. In so doing, a hybrid approach is obtained that is computationally efficient enough to allow for larger-scale modelling of plasma systems while eliminating a source of difficulty in deriving fluid equations applicable to magnetised plasmas. We connect this system of equations to historical asymptotic models and discuss advantages and disadvantages of this approach, including the extension of this parallel-kinetic-perpendicular moment beyond the typical region of validity of these more traditional asymptotic models. This paper forms the first of a multi-part series on this new model, covering the theory and derivation, alongside demonstration benchmarks of this approach that include shocks and magnetic reconnection.

astrophysical plasmas↗

DOE OSTI2026

Magnus force induced magnetic diode effect in skyrmion systems

We show that skyrmions can exhibit a “magnetic diode effect,” where there is a nonreciprocal response in the transport when the magnetic field is reversed. This effect can be achieved for skyrmions moving in channels with a sawtooth potential on one side and a reversed sawtooth potential on the other side. We consider the cases of both spin-transfer torque (STT) and spin–orbit torque (SOT). When the magnetic field is held fixed, the velocity response of the skyrmion is the same for current applied in either direction for both torques, so there is no current diode effect. When the magnetic field is reversed, under STT driving the velocity of the skyrmion reverses and its absolute value changes. Under SOT driving, the velocity remains in the same direction but drops to a much lower value, resulting in negative differential conductivity. For a fixed current, we find a nonreciprocal skyrmion velocity as a function of the applied field’s sign, in analogy to the velocity–current curves observed in the usual diode effect. The nonreciprocity is generated by the Magnus force, which causes skyrmions to interact preferentially with one side of the channel. Since the channel sides have opposite asymmetry, a positive magnetic field can cause the skyrmion to interact with the “hard” asymmetry side of the channel, while a negative magnetic field causes the skyrmion to interact with the “easy” asymmetry side. This geometry could be used to create new kinds of magnetic-field-induced diode effects that can be harnessed in new types of skyrmion-based devices.

36 MATERIALS SCIENCE↗