Search NASA⌕ Search

SEARCH · Search NASA

Results for “Time parallelization”

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 1,513 records · Page 84

Scalable training of graph convolutional neural networks for fast and accurate predictions of HOMO-LUMO gap in molecules

Abstract Graph Convolutional Neural Network (GCNN) is a popular class of deep learning (DL) models in material science to predict material properties from the graph representation of molecular structures. Training an accurate and comprehensive GCNN surrogate for molecular design requires large-scale graph datasets and is usually a time-consuming process. Recent advances in GPUs and distributed computing open a path to reduce the computational cost for GCNN training effectively. However, efficient utilization of high performance computing (HPC) resources for training requires simultaneously optimizing large-scale data management and scalable stochastic batched optimization techniques. In this work, we focus on building GCNN models on HPC systems to predict material properties of millions of molecules. We use HydraGNN, our in-house library for large-scale GCNN training, leveraging distributed data parallelism in PyTorch. We use ADIOS, a high-performance data management framework for efficient storage and reading of large molecular graph data. We perform parallel training on two open-source large-scale graph datasets to build a GCNN predictor for an important quantum property known as the HOMO-LUMO gap. We measure the scalability, accuracy, and convergence of our approach on two DOE supercomputers: the Summit supercomputer at the Oak Ridge Leadership Computing Facility (OLCF) and the Perlmutter system at the National Energy Research Scientific Computing Center (NERSC). We present our experimental results with HydraGNN showing (i) reduction of data loading time up to 4.2 times compared with a conventional method and (ii) linear scaling performance for training up to 1024 GPUs on both Summit and Perlmutter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evaluation of pre-training large language models on leadership-class supercomputers

Large language models (LLMs) have arisen rapidly to the center stage of artificial intelligence as the foundation models applicable to many downstream learning tasks. However, how to effectively build, train, and serve such models for many high-stake and first-principle-based scientific use cases are both of great interests and of great challenges. Moreover, pre-training LLMs with billions or even trillions of parameters can be prohibitively expensive not just for academic institutions, but also for well-funded industrial and government labs. Furthermore, the energy cost and the environmental impact of developing LLMs must be kept in mind. Here, in this work, we conduct a first-of-its-kind performance analysis to understand the time and energy cost of pre-training LLMs on the Department of Energy (DOE)’s leadership-class supercomputers. Employing state-of-the-art distributed training techniques, we evaluate the computational performance of various parallelization approaches at scale for a range of model sizes, and establish a projection model for the cost of full training. Our findings provide baseline results, best practices, and heuristics for pre-training such large models that should be valuable to HPC community at large. We also offer insights and optimization strategies for using the first exascale computing system, Frontier, to train models of the size of GPT-3 and beyond.

97 MATHEMATICS AND COMPUTING↗

Application of an Evolutionary Algorithm for Parameter Optimization in a Gully Erosion Model

Herein we demonstrate how to use model optimization to determine a set of best-fit parameters for a landform model simulating gully incision and headcut retreat. To achieve this result we employed the Covariance Matrix Adaptation Evolution Strategy (CMA-ES), an iterative process in which samples are created based on a distribution of parameter values that evolve over time to better fit an objective function. CMA-ES efficiently finds optimal parameters, even with high-dimensional objective functions that are non-convex, multimodal, and non-separable. We ran model instances in parallel on a high-performance cluster, and from hundreds of model runs we obtained the best parameter choices. This method is far superior to brute-force search algorithms, and has great potential for many applications in earth science modeling. We found that parameters representing boundary conditions tended to converge toward an optimal single value, whereas parameters controlling geomorphic processes are defined by a range of optimal values.

97 MATHEMATICS AND COMPUTING↗

Ultrafast laser-induced spin–lattice dynamics in the van der Waals antiferromagnet CoPS 3

CoPS 3 stands out in the family of the van der Waals antiferromagnets XPS 3 (X = Mn, Ni, Fe, and Co) due to the unquenched orbital momentum of the magnetic Co 2+ ions, which is known to facilitate the coupling of spins to both electromagnetic waves and lattice vibrations. Here, using a time-resolved magneto-optical pump–probe technique, we experimentally study the ultrafast laser-induced dynamics of mutually correlated spins and lattice. It is shown that a femtosecond laser pulse acts as an ultrafast heater and, thus, results in the melting of the antiferromagnetic order. At the same time, the resonant pumping of the 4 T 1g → 4 T 2g electronic transition in Co 2+ ions effectively changes their orbital momentum, giving rise to a mechanical force that moves the ions in the direction parallel to the orientation of their spins, thus generating a coherent B g phonon mode at the frequency of about 4.7 THz.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dissecting protein tyrosine phosphatase signaling by engineered chemogenetic control of its activity

Protein tyrosine phosphatases (PTPases) are critical mediators of dynamic cell signaling. A tool capable of identifying transient signaling events downstream of PTPases is essential to understand phosphatase function on a physiological time scale. We report a broadly applicable protein engineering method for allosteric regulation of PTPases. This method enables dissection of transient events and reconstruction of individual signaling pathways. Implementation of this approach for Shp2 phosphatase revealed parallel MAPK and ROCK II dependent pathways downstream of Shp2, mediating transient cell spreading and migration. Furthermore, we show that the N-SH2 domain of Shp2 regulates MAPK-independent, ROCK II-dependent cell migration. Engineered targeting of Shp2 activity to different protein complexes revealed that Shp2-FAK signaling induces cell spreading whereas Shp2-Gab1 or Shp2-Gab2 mediates cell migration. We identified specific transient morphodynamic processes induced by Shp2 and determined the role of individual signaling pathways downstream of Shp2 in regulating these events. Broad application of this approach is demonstrated by regulating PTP1B and PTP-PEST phosphatases.

59 BASIC BIOLOGICAL SCIENCES↗

Magnetic field dependence of low-energy magnons, anisotropic heat conduction, and spontaneous relaxation of magnetic domains in the cubic helimagnet ZnCr 2 Se 4

Here, anisotropic low-temperature properties of the cubic spinel helimagnet ZnCr 2 Se 4 in the single-domain spin-spiral state are investigated by a combination of neutron scattering, thermal conductivity, ultrasound velocity, and dilatometry measurements. In an applied magnetic field, neutron spectroscopy shows a complex and nonmonotonic evolution of the spin-wave spectrum across the quantum-critical point that separates the spin-spiral phase from the field-polarized ferromagnetic phase at high fields. A tiny spin gap of the pseudo-Goldstone magnon mode, observed at wave vectors that are structurally equivalent but orthogonal to the propagation vector of the spin helix, vanishes at this quantum critical point, restoring the cubic symmetry in the magnetic subsystem. The anisotropy imposed by the spin helix has only a minor influence on the lattice structure and sound velocity but has a much stronger effect on the heat conductivities measured parallel and perpendicular to the magnetic propagation vector. The thermal transport is anisotropic at T≲ 2K, highly sensitive to an external magnetic field, and likely results directly from magnonic heat conduction. We also report long-time thermal relaxation phenomena, revealed by capacitive dilatometry, which are due to magnetic domain motion related to the destruction of the single-domain magnetic state, initially stabilized in the sample by the application and removal of magnetic field. Our results can be generalized to a broad class of helimagnetic materials in which a discrete lattice symmetry is spontaneously broken by the magnetic order.

36 MATERIALS SCIENCE↗

Visualizing an Exascale Data Center Digital Twin: Considerations, Challenges and Opportunities

Digital twins are an excellent tool to model, visualize, and simulate complex systems, to understand and optimize their operation. In this work, we present the technical challenges of real-time visualization of a digital twin of the Frontier supercomputer.We show the initial prototype and current state of the twin and highlight technical design challenges of visualizing such a large High Performance Computing (HPC) system. The goal is to understand the use of augmented reality as a primary way to extract information and collaborate on digital twins of complex systems. This leverages the spatio-temporal aspect of a 3D representation of a digital twin, with the ability to view historical and real-time telemetry, triggering simulations of a system state and viewing the results, which can be augmented via dashboards for details. Finally, we discuss considerations and opportunities for augmented reality of digital twins of large-scale, parallel computers.

Maiterth, Matthias↗

SU-G-TeP1-15: Toward a Novel GPU Accelerated Deterministic Solution to the Linear Boltzmann Transport Equation

Purpose: To develop a Graphic Processor Unit (GPU) accelerated deterministic solution to the Linear Boltzmann Transport Equation (LBTE) for accurate dose calculations in radiotherapy (RT). A deterministic solution yields the potential for major speed improvements due to the sparse matrix-vector and vector-vector multiplications and would thus be of benefit to RT. Methods: In order to leverage the massively parallel architecture of GPUs, the first order LBTE was reformulated as a second order self-adjoint equation using the Least Squares Finite Element Method (LSFEM). This produces a symmetric positive-definite matrix which is efficiently solved using a parallelized conjugate gradient (CG) solver. The LSFEM formalism is applied in space, discrete ordinates is applied in angle, and the Multigroup method is applied in energy. The final linear system of equations produced is tightly coupled in space and angle. Our code written in CUDA-C was benchmarked on an Nvidia GeForce TITAN-X GPU against an Intel i7-6700K CPU. A spatial mesh of 30,950 tetrahedral elements was used with an S4 angular approximation. Results: To avoid repeating a full computationally intensive finite element matrix assembly at each Multigroup energy, a novel mapping algorithm was developed which minimized the operations required at each energy. Additionally, a parallelized memory mapping for the kronecker product between the sparse spatial and angular matrices, including Dirichlet boundary conditions, was created. Atomicity is preserved by graph-coloring overlapping nodes into separate kernel launches. The one-time mapping calculations for matrix assembly, kronecker product, and boundary condition application took 452±1ms on GPU. Matrix assembly for 16 energy groups took 556±3s on CPU, and 358±2ms on GPU using the mappings developed. The CG solver took 93±1s on CPU, and 468±2ms on GPU. Conclusion: Three computationally intensive subroutines in deterministically solving the LBTE have been formulated on GPU, resulting in two orders of magnitude speedup. Funding support from Natural Sciences and Engineering Research Council and Alberta Innovates Health Solutions. Dr. Fallone is a co-founder and CEO of MagnetTx Oncology Solutions (under discussions to license Alberta bi-planar linac MR for commercialization).

60 APPLIED LIFE SCIENCES↗

Dynamics of correlation-frozen antinodal quasiparticles in superconducting cuprates

Many puzzling properties of high–critical temperature ( T c ) superconducting (HTSC) copper oxides have deep roots in the nature of the antinodal quasiparticles, the elementary excitations with wave vector parallel to the Cu–O bonds. These electronic states are most affected by the onset of antiferromagnetic correlations and charge instabilities, and they host the maximum of the anisotropic superconducting gap and pseudogap. We use time-resolved extremeultraviolet photoemission with proper photon energy (18 eV) and time resolution (50 fs) to disclose the ultrafast dynamics of the antinodal states in a prototypical HTSC cuprate. After photoinducing a nonthermal charge redistribution within the Cu and O orbitals, we reveal a dramatic momentum-space differentiation of the transient electron dynamics. Whereas the nodal quasiparticle distribution is heated up as in a conventional metal, new quasiparticle states transiently emerge at the antinodes, similarly to what is expected for a photoexcited Mott insulator, where the frozen charges can be released by an impulsive excitation. This transient antinodal metallicity is mapped into the dynamics of the O-2p bands, thus directly demonstrating the intertwining between the lowand high-energy scales that is typical of correlated materials. Our results suggest that the correlation-driven freezing of the electrons moving along the Cu–O bonds, analogous to the Mott localization mechanism, constitutes the starting point for any model of high- T c superconductivity and other exotic phases of HTSC cuprates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

PETSc TSAdjoint: A Discrete Adjoint ODE Solver for First-Order and Second-Order Sensitivity Analysis

Here, we present a new software system PETSc TSAdjoint for first-order and second order adjoint sensitivity analysis of time-dependent nonlinear differential equations. The derivative calculation in PETSc TSAdjoint is essentially a high-level algorithmic differentiation process. The adjoint models are derived by differentiating the timestepping algorithms and implementing them based on the parallel infrastructure in PETSc. Full differentiation of the library code, including MPI routines, is avoided, and users do not need to derive their own adjoint models for their specific applications. PETSc TSAdjoint can compute the first-order derivative, that is, the gradient of a scalar functional, and the Hessian-vector product, which carries second-order derivative information, while requiring minimal input (a few callbacks) from the users. The adjoint model employs optimal checkpointing schemes in a manner that is transparent to users. Finally, usability, efficiency, and scalability are demonstrated through examples from a variety of applications.

79 ASTRONOMY AND ASTROPHYSICS↗

Human Factors for Advanced Reactors

Existing light water reactors in the U.S. are primarily large baseload electricity generating facilities. The concept of operations for these plants remains largely unchanged since the advent of commercial nuclear power—the main control room serves as the hub of plant activities and is staffed with multiple licensed operators who work in tandem under the shift supervisor, and staff such as field workers support the control room remotely. While newer plants have brought the advent of digital human-machine interfaces to replace earlier analog and mechanical instrumentation and controls, much of the control process remains unchanged and manual. It is simply a newer version of legacy concepts. Advanced reactors potentially bring considerable changes to the size, fuel type, automation, and staffing of nuclear power plants, necessitating a fundamental shift not just from analog to digital, but further from human to automation, from onsite to remote, from control to monitoring, and from many to few operators. Despite this multitude of parallel evolutions in reactor designs, many of the vendors developing the next generation of reactors represent smaller research and development enterprises. It is therefore not feasible to address all aspects of plant design at the same time. In particular, the competing design aspects of new reactors present a significant challenge to the development of robust and human factored systems at the plant. As vendors develop new reactor designs, much of the early focus is naturally on the fuel and reactor system technology. Looming behind these early advances is the daunting prospect of first-of-a-kind control concepts that have not yet been developed or validated. A failure to address the human element of reactor design early will lead to missed opportunities. The quickest development process is the replication of existing concepts of operations at legacy plants, even when such systems were long ago surpassed by better human-machine technologies outside the nuclear industry. Conversely, attempting to undertake novel concepts of operations late in the design life cycle of a plant could result in protracted development efforts and delays in licensing and deployment. This does not have to happen, and it is imperative that human factors be considered now, early in the design of new reactors.

99 GENERAL AND MISCELLANEOUS↗

A study of Channeling, Volume Reflection and Volume Capture of 3.35 - 14.0 GeV Electrons in a bent Silicon Crystal

We present the experimental data and analysis of experiments conducted at SLAC National Accelerator Laboratory investigating the processes of channeling, volume-reflection and volume-capture along the (111) plane in a strongly bent quasi-mosaic silicon crystal. Additionally, these phenomena were investigated at 5 energies: 3.35, 4.2, 6.3, 10.5 and 14.0 GeV with a crystal with bending radius of 0.15m, corresponding to curvatures of 0.070, 0.088, 0.13, 0.22 and 0.29 times the critical curvature respectively. We have extracted important parameters describing the channeling process such as the dechanneling length, the angle of volume reflection, the surface transmission and the widths of the distribution of channeled particles parallel and orthogonal to the plane.

ACCPHY↗

Microgrid energy scheduling under uncertain extreme weather: Adaptation from parallelized reinforcement learning agents

Microgrids are useful solutions for integrating renewable energy resources and providing seamless green electricity to minimize carbon footprint. In recent years, extreme weather events happened often worldwide and caused significant economic and societal losses. Such events bring uncertainties to the microgrid energy scheduling problems and increase the challenges of microgrid operation. Traditional optimization approaches suffer from the inaccuracy of the uncertain microgrid model and the unseen events. Existing reinforcement learning (RL) - based approaches are also hampered by the limited generalization and the increasing computational burden when stochastic formulations are required to accommodate the uncertainties. This paper proposes a new parallelized reinforcement learning (PRL) method based on the probabilistic events to handle the microgrid energy uncertainties. Specifically, several local learning agents are employed to interact with pertinent microgrid environments in a distributed manner and report outcomes to the global agent, which will optimize microgrid energy resources online during extreme events. The stochastic microgrid energy optimization problem is reformulated to include all possible scenarios with probabilities. The advantage estimate functions of learning agents are designed with a backward sweep to transfer the outcomes to the value function updating process. Two simulation studies, stochastic optimization and online testing, are performed to compare with several existing RL approaches. Results substantiate that the proposed PRL method can achieve up to 20% optimization performance improvement with 4 and 28 times less computation cost than Q-learning with experience replay and multi-agent Q-learning approaches, respectively.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Digital Time-Resolved Spot Diagnostic

The digital time-resolved spot diagnostic (D-TRSD) is intended to be a fast-response, time-dependent beam diagnostic for measuring the spot size emitted from radiographic x-ray sources. It will measure spot sizes from submillimeters up to several millimeters for photon energies in the range of 150 keV to 20 MeV with doses as low as 0.03 rad up to 500 rad. The D-TRSD performs fast data acquisition as radiographic x-ray sources typically produce pulses in the tens of nanosecond ranges. The system will be optimized to capture multiple individual pulses in time within the same data record. The data collected will be used to support beam tuning and source optimization with on-the-fly evaluation of radiographic source performance. The diagnostic is based on penumbral imaging techniques. A scintillating fiber array of parallel aligned fibers is positioned in the beam path with a high-Z rolled edge inserted between the source and the array, thereby creating the penumbral x-ray signal that is captured by the D-TRSD system. The light from the scintillation conversion is recorded by individual avalanche photodiode detectors integrated into custom digitizer electronics. The penumbra is used to calculate the spot of the beam on the axis orthogonal to the fiber array based on the magnification from the rolled edge, the resolution of the spaced fibers, and the rotation of the array normal to the beam path. Many arrays can be used to capture additional spot widths about different rotations. The D-TRSD system is modular, software-driven, and easier to use than legacy TRSD systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Parallel Spectral Acquisition with Orthogonal ICR Cells

FT-based high performance mass analyzers yield increased resolving power and mass measurement accuracy, yet require increased duration of signal acquisition that can limit many applications. The implementation of stronger magnetic fields, multiple detection electrodes for harmonic signal detection, and an array of multiple mass analyzers arranged along the magnetic field axis have been used to decrease required acquisition time. The results presented here show that multiple ion cyclotron resonance (ICR) mass analyzers can also be implemented orthogonal to the central magnetic field axis. The orthogonal ICR cell system presented here consisting of two cells (master and slave cells) was constructed with printed circuit boards and installed within a single superconducting magnet and vacuum system. A master cell was positioned, as is normally done with ICR cells, on the central magnetic field axis and a slave cell was located off this central axis, but directly adjacent and alongside the master cell. To achieve ion transfer between cells, ions that were initially trapped in the master cell were drifted across the magnetic field into the slave cell with application of a small DC field applied perpendicularly to the magnetic field axis. A subsequent population of ions was injected and accumulated in the master cell. Simultaneous excitation of cyclotron motion of ions in both cells was carried out; ICR signals from each cell were independently amplified and recorded in parallel. Presented here are the initial results of successful parallel spectral acquisition with this orthogonal dual ICR cell array. .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Predicting the effect of angular momentum on the dissociation dynamics of highly rotationally excited radical intermediates

We present a model which accurately predicts the net speed distributions of products resulting from the unimolecular decomposition of rotationally excited radicals. The radicals are produced photolytically from a halogenated precursor under collision-free conditions so they are not in a thermal distribution of rotational states. The accuracy relies on the radical dissociating with negligible energetic barrier beyond the endoergicity. We test the model predictions using previous velocity map imaging and crossed laser-molecular beam scattering experiments that photolytically generated rotationally excited CD{sub 2}CD{sub 2}OH and C{sub 3}H{sub 6}OH radicals from brominated precursors; some of those radicals then undergo further dissociation to CD{sub 2}CD{sub 2} + OH and C{sub 3}H{sub 6} + OH, respectively. We model the rotational trajectories of these radicals, with high vibrational and rotational energy, first near their equilibrium geometry, and then by projecting each point during the rotation to the transition state (continuing the rotational dynamics at that geometry). This allows us to accurately predict the recoil velocity imparted in the subsequent dissociation of the radical by calculating the tangential velocities of the CD{sub 2}CD{sub 2}/C{sub 3}H{sub 6} and OH fragments at the transition state. The model also gives a prediction for the distribution of angles between the dissociation fragments’ velocity vectors and the initial radical’s velocity vector. These results are used to generate fits to the previously measured time-of-flight distributions of the dissociation fragments; the fits are excellent. The results demonstrate the importance of considering the precession of the angular velocity vector for a rotating radical. We also show that if the initial angular momentum of the rotating radical lies nearly parallel to a principal axis, the very narrow range of tangential velocities predicted by this model must be convoluted with a J = 0 recoil velocity distribution to achieve a good result. The model relies on measuring the kinetic energy release when the halogenated precursor is photodissociated via a repulsive excited state but does not include any adjustable parameters. Even when different conformers of the photolytic precursor are populated, weighting the prediction by a thermal conformer population gives an accurate prediction for the relative velocity vectors of the fragments from the highly rotationally excited radical intermediates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ultrafast Time-Resolved Studies on Fluorescein for Recognition Strands Architecture in Amyloid Fibrils

Protein aggregation is associated with numerous devastating diseases such as Alzheimer’s, Parkinson’s, and prion diseases. Development of therapeutics would benefit from knowledge of the structural organization of protein molecules in these amyloid aggregates, particularly in their aqueous biological milieu. However, detailed structural studies to date have been mainly on the solid state and have required large quantities of purified aggregate. Moreover, these conventional methods require the aggregated assembly to remain structurally stable over days or weeks required to perform the experiment, whereas the pathologically relevant species of in vivo aggregates may be shorter lived. Here, we show the organization of protein chains in dissolved amyloid aggregates can be readily determined spectroscopically using minute quantities of fluorescein-labeled protein segments in a matter of minutes. Specifically, we investigated the possibility of using the ultrafast dynamics of fluorescein to distinguish among three categories of β-sheet geometry: (1) antiparallel in-register, (2) parallel in-register, or (3) antiparallel out-of-register. Fluorescein, the most commonly used staining dye in biology and medicine, was covalently attached to the N-termini of peptide sequences selected from a library of known amyloid crystal structures. We investigated the aggregates in solution using steady-state and time-resolved absorption and fluorescence spectroscopy. We found that the dynamics of fluorescein relaxation from the excited state revealed amyloid structure-specific information. Particularly, the nonfluorescent cation form of fluorescein showed remarkable sensitivity to local environments created during aggregation. We demonstrate that time-resolved absorption is capable of differentiating strand organization in β-sheet aggregates when strong intermolecular coupling between chromophores occurs. Here, this approach can be useful in optical recognition of specific fibril architectures of amyloid aggregates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Infrared MBE-Grown HgCdTe Focal Plane Arrays and Cameras After High Energy Neutron Irradiation

HgCdTe is one of the most important materials for the fabrication of infrared detectors and focal plane arrays (FPAs) deployed in environments where high-energy particles, such as protons and neutrons, are present. We designed and fabricated HgCdTe-based FPAs that can be used in high neutron radiation environments and we measured their characteristics. The influence of the radiation on the infrared FPAs and cameras is investigated. HgCdTe material and devices are capable of maintaining high performances in a high energy neutron irradiation environment. For MWIR FPA directly facing a $2.59 × 10^8$ n/cm$^2$ s neutron flux beam (with the highest energy 66 MeV) for 1 h, the noise equivalent differential temperature (NEDT) increased ~8 times after irradiation. However, NEDT decreased to 33 mK (compared to the original value of 21 mK) after one warming-up (to room temperature) and cooling-down cycle. The NEDT for the MWIR FPAs mounted parallel to the beam did not degrade (16 mK and 28 mK before irradiation, changed to 18 mK and 26 mK after irradiation, respectively).

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗