Molecular latent space simulators
Latent space simulators learn kinetic models for atomistic simulations and generate novel trajectories at six orders of magnitude lower cost.
SEARCH · Search NASA
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.
Latent space simulators learn kinetic models for atomistic simulations and generate novel trajectories at six orders of magnitude lower cost.
Recent discoveries of hyperbolic surface plasmon polaritons (SPPs) in 1T′ WTe 2 have garnered significant attention in 2D materials and nanophotonics. In this study, we employ finite-element simulations to investigate the real-space characteristics of hyperbolic SPPs in thin WTe 2 flakes. Our results show that the SPPs exhibit a pronounced sensitivity to excitation energy and sample thickness. By analyzing the plasmonic field patterns, we extract key plasmonic parameters including plasmon wavelengths, hyperbolic angles, and plasmonic figures of merit. In addition, we examine SPP modes in stacked WTe 2 flakes with varying twist angles, demonstrating that the plasmonic field patterns can be effectively tuned by adjusting the twist angle. Notably, as the twist angle increases, the SPPs undergo a topological transition from open hyperbolic modes to closed elliptic modes. This twist-engineering capability offers promising potential for the development of tunable plasmonic devices based on WTe 2 .
Not Available
Not Available
The CO2 capture efficiency in solvent-based carbon capture systems (CCSs) critically depends on the gas-solvent interfacial area (IA), making maximization of IA a foundational challenge in CCS design. While the IA associated with a particular CCS design can be estimated via a computational fluid dynamics (CFD) simulation, using CFD to derive the IAs associated with numerous CCS designs is prohibitively costly. However, previous works such as Deep Fluids (Kim et al., 2019) show that large simulation speed-ups and low error are achievable by replacing CFD simulators with neural-network (NN) surrogates that mimic the CFD simulation process. This raises the possibility for a fast, accurate replacement for a CFD simulator, and thus computationally feasible IA-based CCS-design optimization. As such, here, we explore whether an existing NN-surrogate approach (and variants we develop) can successfully be applied to our complex carbon-capture CFD simulations, with the ultimate goal of obtaining a fast and accurate simulator for our CCS-design application. Our experiments build on the Deep Fluids approach and find that resulting surrogates can produce large speed ups (4000x) while maintaining IA relative errors as low as 4% on unseen CCS configurations (interpolating between configurations seen during training). Thus, despite less faithfulness to the underlying physics of the problem, NN surrogates may be a promising tool for our CCS design optimization problem. Notably, though, the Deep Fluids approach has limitations for our application (such as model non-transferability to CCS-packing changes), which we discuss. We conclude with potential directions for future work that may, like innovations we introduced here (e.g., transformer-based dynamics prediction), improve performance on our complicated dataset of CCS CFD simulations.
Explore the source record for details and available documents.
A 2.5 MeV proton injector is being constructed for the IOTA ring at Fermilab to study the interaction of nonlinear integrable optics (NIO) with high space charge beams. Space charge in the transport line from the RFQ to the injection location has a significant current dependent effect on the phase space. Simulation studies to support efficient injection of intense bunches into IOTA are presented, included schemes to inject directly into NIO lattices.
A 2.5 MeV proton injector is being constructed for the IOTA ring at Fermilab to study the interaction of nonlinear integrable optics (NIO) with high space charge beams. Space charge in the transport line from the RFQ to the injection location has a significant current dependent effect on the phase space. Simulation studies to support efficient injection of intense bunches into IOTA are presented, included schemes to inject directly into NIO lattices.
Injection simulations of space charge dominated proton beams in IOTA ABSTRACT = A 2.5 MeV proton injector has been constructed for the IOTA ring at Fermilab to study the interaction of nonlinear integrable optics (NIO) with high space charge beams. Space charge in the transport line from the RFQ to the injection location has a significant current dependent effect on the phase space. Simulation studies to support efficient injection of intense bunches into IOTA are presented, included schemes to inject directly into NIO lattices.
The recent development of differentiable simulation codes for particle accelerators has enabled gradient-based workflows that promise finer control and more realistic modeling of accelerator facilities. However, when using reverse-mode automatic differentiation, the memory usage continuously increases during the simulation, and can potentially exceed the available hardware memory - especially when costly space charge computation is included. To study the memory requirements for differentiable simulations, we have implemented space charge in Cheetah, a PyTorch-based beam tracking code that supports reverse-mode differentiation. We find that the memory usage for reverse-mode differentiation grows linearly with the number of macroparticles and cells, and that it is proportional to the number of space charge kicks involved in the simulation. This general scaling can be used to evaluate whether a given differentiable simulation is feasible given hardware memory constraints.
The future high energy physics program at Fermilab requires that the proton complex operate with beam bunch intensities four times larger than is currently handled. At these intensities space charge nonlinear defocussing effects cause unacceptable particle losses especially in the low energy rapid-cycling-synchrotron (RCS) Booster. Focusing electron lens elements may offer a solution by providing partial space charge compensation but there is a need for detailed simulations as this technique has not been demonstrated. We report on high fidelity numerical 6D space charge simulations in a model accelerator lattice with a record high space charge tune shift approaching unity.
The role and mitigation of space charge effects are important aspects of the beam physics research to be performed at the Integrable Optics Test Accelerator (IOTA) at Fermilab. The impact of nonlinear integrability (partial and complete) on space charge driven incoherent and coherent resonances will be a new feature of this accelerator and has the potential to influence the design of future low energy proton synchrotrons. In this report we model transverse space charge effects using the PIC code pyORBIT. First we benchmark the single particle tracking results against MADX with checks on symplecticity, tune footprints, and dynamic aperture in a partially integrable lattice realized with a special octupole string insert. Our space charge calculations begin with an examination of the 4D symplecticity. Short term tracking is done first with the initial bare lattice and then with a matching of the rms values with space charge. Next, we explore slow initialization of charge as a technique to establish steady state and reduce emittance growth and beam loss following injection into a mismatched lattice. We establish values of space charge simulation parameters so as to ensure numerical convergence. Finally, we compare the simulated space charge induced tune shifts and footprints against theory.
Integrable optics is an innovation in particle accelerator design that enables strong nonlinear focusing without generating parametric resonances. We use a Synergia space-charge simulation to investigate the application of integrable optics to a high-intensity hadron ring that could replace the Fermilab Booster. We find that incorporating integrability into the design suppresses the beam halo generated by a mismatched KV beam. Our integrable rapid cycling synchrotron (iRCS) design includes other features of modern ring design such as low momentum compaction factor and harmonically canceling sextupoles. Experimental tests of high-intensity beams in integrable lattices will take place over the next several years at the Fermilab Integrable Optics Test Accelerator (IOTA) and the University of Maryland Electron Ring (UMER).
Symplectic simulation of space-charge effects is crucial for the design and operation of high-intensity particle accelerators. Traditional methods for simulating these effects are often computationally expensive, resulting in significant overhead. In this work, we introduce a generative model based on a U-Net architecture within a generative adversarial network framework to efficiently simulate space-charge effects. The model is trained to predict the transverse multiparticle space-charge Hamiltonian, which can be physically computed using a gridless spectral method. The one-step symplectic transverse transfer map for the particles is then obtained by differentiating the predicted Hamiltonian. Benchmarking results demonstrate that this generative model achieves an order of magnitude higher computational efficiency compared to the spectral method, providing a highly efficient alternative for simulating space-charge effects with a large number of particles. By maintaining symplecticity, the model effectively preserves the phase-space structure and mitigates nonphysical errors in long-term simulations. This model has been integrated into jutrack, a novel autodifferentiable accelerator modeling code developed in the julia programming language.
The results obtained by experimentally studying gamma rays emitted by samples prepared as analogs of planetary matter and irradiated with thermal neutrons are presented. The intensities of spectral lines of gamma rays emitted by such samples differing in chemical composition are compared. These results will be used in processing data on gamma-ray spectra of the Moon and Mercury from measurements performed onboard spacecrafts with the aim of studying the composition of the surface of these celestial bodies.
We present a systematic attempt to study magnetic null points and the associated magnetic energy conversion in kinetic particle-in-cell simulations of various plasma configurations. We address three-dimensional simulations performed with the semi-implicit kinetic electromagnetic code iPic3D in different setups: variations of a Harris current sheet, dipolar and quadrupolar magnetospheres interacting with the solar wind, and a relaxing turbulent configuration with multiple null points. Spiral nulls are more likely created in space plasmas: in all our simulations except lunar magnetic anomaly (LMA) and quadrupolar mini-magnetosphere the number of spiral nulls prevails over the number of radial nulls by a factor of 3–9. We show that often magnetic nulls do not indicate the regions of intensive energy dissipation. Energy dissipation events caused by topological bifurcations at radial nulls are rather rare and short-lived. The so-called X-lines formed by the radial nulls in the Harris current sheet and LMA simulations are rather stable and do not exhibit any energy dissipation. Energy dissipation is more powerful in the vicinity of spiral nulls enclosed by magnetic flux ropes with strong currents at their axes (their cross sections resemble 2D magnetic islands). These null lines reminiscent of Z-pinches efficiently dissipate magnetic energy due to secondary instabilities such as the two-stream or kinking instability, accompanied by changes in magnetic topology. Current enhancements accompanied by spiral nulls may signal magnetic energy conversion sites in the observational data.
Defocusing repulsive forces due to self space charge fields lead to degradation of high-intensity particle beams. Being of particular concern for low- and medium-energy proton beams, they result in emittance growth, beam halo formation, and beam loss. They set stringent limits on the intensity of frontier accelerators; therefore, the mitigation of space charge effects is a crucial challenge to improve proton beam intensity. The space charge effects in a positively charged proton beam can be effectively compensated using negatively charged electron columns. In this paper, we present the results of simulations using Synergia of the Electron Column lattice for IOTA. Beam loss due to space charge effects and aperture restrictions have been studied, as well as bunch formation and matching using an adiabatic ramp of the RF cavity. The results show the need for space charge compensation, and provide the basis for integration of the Synergia and Warp codes in order to form a complete simulation of space charge compensation using an Electron Column in IOTA.
PPMLR-MHD is a new magnetohydrodynamics (MHD) model used to simulate the interactions of the solar wind with the magnetosphere, which has been proved to be the key element of the space weather cause-and-effect chain process from the Sun to Earth. Compared to existing MHD methods, PPMLR-MHD achieves the advantage of high order spatial accuracy and low numerical dissipation. However, the accuracy comes at a cost. On one hand, this method requires more intensive computation. On the other hand, more boundary data is subject to be transferred during the process of simulation. In this work, we present a parallel hybrid solution of the PPMLR-MHD model implemented using the computing capabilities of both CPUs and GPUs. We demonstrate that our optimized implementation alleviates the data transfer overhead by using GPU Direct technology and can scale up to 151 processes and achieve significant performance gains by distributing the workload among the CPUs and GPUs on Titan at Oak Ridge National Laboratory. The performance results show that our implementation is fast enough to carry out highly accurate MHD simulations in real time.