Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Quantum Networking”

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.

1,368 records · Page 76

Logical Majorana zero modes in a nanowire network

We present a scheme to use physical Majorana quasizero modes at each junction of a two-dimensional nanowire network to build a logical Majorana zero mode, the location of which is controllable through gate voltages. The wire network is a way to realize a recent proposal to imprint a Kekulé vortex pattern on a honeycomb lattice via gate voltages. We show that a specific type of junction, other than a naive Y or T junction, is needed to realize, without breaking time-reversal symmetry, an artificial “graphene” system with Majorana fermions instead of complex ones at each site. The junction we propose (i) traps exactly one physical Majorana (quasi)zero mode at each site of either a brick wall or honeycomb lattice and (ii) allows this mode to hybridize with all three neighboring sites. Using a lattice of these junctions and starting from an electronic, tight-binding model for the wires, we imprint the voltage patterns corresponding to the Kekulé vortex and observe the emergence of the logical Majorana zero mode at the vortex core. By suitable adiabatic modulations of gate voltages as a function of time, the positions of individual Kekulé vortices can be moved, and in this manner, non-Abelian braiding of the logical Majoranas bound to the vortex cores can be performed. We provide the range of parameters where these excitations could be realized experimentally.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Data-driven model for divertor plasma detachment prediction

We present a fast and accurate data-driven surrogate model for divertor plasma detachment prediction leveraging the latent feature space concept in machine learning research. Our approach involves constructing and training two neural networks: an autoencoder that finds a proper latent space representation (LSR) of plasma state by compressing the multi-modal diagnostic measurements and a forward model using multi-layer perception (MLP) that projects a set of plasma control parameters to its corresponding LSR. By combining the forward model and the decoder network from autoencoder, this new data-driven surrogate model is able to predict a consistent set of diagnostic measurements based on a few plasma control parameters. In order to ensure that the crucial detachment physics is correctly captured, highly efficient 1D UEDGE model is used to generate training and validation data in this study. The benchmark between the data-driven surrogate model and UEDGE simulations shows that our surrogate model is capable of providing accurate detachment prediction (usually within a few per cent relative error margin) but with at least four orders of magnitude speed-up, indicating that performance-wise, it has the potential to facilitate integrated tokamak design and plasma control. Comparing with the widely used two-point model and/or two-point model formatting, the new data-driven model features additional detachment front prediction and can be easily extended to incorporate richer physics. This study demonstrates that the complicated divertor and scrape-off-layer plasma state has a low-dimensional representation in latent space. Understanding plasma dynamics in latent space and utilising this knowledge could open a new path for plasma control in magnetic fusion energy research.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Initial Test Results for NASA Goddard’s Low-Cost Optical Terminal Adaptive Optics System

This presentation presents the testing approach and initial laboratory test results of the Adaptive Optics (AO) Subassembly (AOS) for the Low-Cost Optical Terminal (LCOT), a low-cost and flexible optical communications ground terminal. Leveraging commercial components and subassemblies wherever possible, LCOT is being developed by the NASA Goddard Space Flight Center (GSFC) Advanced Communications Capabilities for Exploration and Science Systems (ACCESS) project. An economical ground terminal allows NASA to create a global optical terminal network and enables the expansion of direct-to-Earth (DTE) optical communications. Located at the Goddard Geophysical and Astronomical Observatory (GGAO), LCOT is in the development phase and will receive first light during initial testing beginning in 2023. This new optical terminal includes a 700 mm Ritchey-Chretien commercial telescope, able to track Low Earth Orbit (LEO) spacecraft, receive instruments specified from 1500 - 1600 nm, and an uplink transmit optical system. Demodulation of optical coherent modulation formats requires downlink signal coupling into single-mode fiber. Due to atmospheric turbulence, signals experience wavefront distortions creating speckle patterns larger than the fiber core and high jitter requiring precise beam steering. AO is necessary to correct atmospheric turbulence effects and allow efficient coupling of received signals. In 2020, GSFC contracted with General Atomics to provide the LCOT AO system. This system has been delivered to GSFC in September 2021 and has been tested in laboratory to evaluate performance. This paper describes the laboratory testing approach, turbulence phase plate design, test results, and the AO field testing plan when installed on the LCOT telescope.

Adaptive Optics↗

A General Solution Route to Nanoporous Metal Oxide Films

Metal oxide semiconductors are of interest as efficient, stable, and low-cost photoanode materials for photoelectrochemical water splitting, but their characteristically poor charge transport properties remain a fundamental challenge to practical use. Nanoporous metal oxide films that feature open bicontinuous networks of nanocrystals and nanopores can overcome such inefficient charge transport to enable high-performance photoanodes. Here, this paper describes a general method to make high-quality nanoporous films of metal oxides by spin coating and calcining molecular inks that contain a porosity-generating block copolymer. Phase-pure nanoporous films of BiFeO 3 , FeWO 4 , WO 3 , Fe 2 O 3 , and TiO 2 serve to demonstrate the versatility of the approach. It is demonstrated that the crystallite size, film porosity, and film thickness can be independently tuned by adjusting the ink composition and film processing conditions. The method is extended to fabricate core–shell nanoporous films consisting of a nanoporous film coated in a thin shell of a second metal oxide, using WO 3 –BiVO 4 and BiFeO 3 –BiVO 4 as examples. Given its simplicity and flexibility, this solution-phase route should prove useful for making nanoporous films of many different materials for a variety of applications, including energy conversion and storage, catalysis, and chemical sensing.

coating materials↗

Massive all-atom analysis of 2D materials with quantum properties (Final report)

Improvements in microscopy have enabled the acquisition of data at a scale that is difficult to process manually, making automated machine learning approaches to analyzing experimental images essential. In this project, we developed and applied machine learning (ML) workflows for atomic resolution scanning transmission electron microscopy (STEM) images. This development included improving both methodology as well as generating user-friendly codes. We developed machine learning architectures which, after training, automatically identify the location and types of defects throughout a material. We used these data to produce class-averaged images of 2D atomic coordinates with up to 0.3 pm precision, uncovering the structure and oscillations of long-range strain fields around point defects in WSe 2-2x Te 2x . We also resolved a long-standing problem in this field in the training of ML models, a lack of labeled experimental data, by developing a cycle-GAN that transformed simulated-generated labeled data into labeled data indistinguishable from experiment and therefore suitable for training. This removed the remaining parts of the ML data processing workflow where human intervention was still critical and therefore a bottleneck to working at scale. Codes have been developed and released for this full machine learning workflow. ML approaches to partially automate STEM acquisition were also developed. Finally we applied ML and other advanced data processing methods to several materials science problems in two-dimensional materials, including studying the evolution of hyperuniformity with defect concentration in WSe2, understanding phase transformations in transition metal dichalcogenides during in-situ heating in the STEM, and exploring how 2D interfaces transform from twisted into aligned structures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Improved energies and local energies with weighted variational Monte Carlo

Neural network parametrizations have increasingly been used to represent the ground and excited states in variational Monte Carlo (VMC) with promising results. However, traditional VMC methods only optimize the wave function in regions of peak probability. The wave function is uncontrolled in the tails of the probability distribution, which can limit the accuracy of the trained wave function. To improve the approximation accuracy in the probability tails, this paper interprets VMC as a gradient flow in the space of wave functions, followed by a projection step. From this perspective, arbitrary probability distributions can be used in the projection step, allowing the user to prioritize accuracy in different regions of state space. Motivated by this theoretical perspective, the paper tests a weighted VMC method on the antiferromagnetic Heisenberg model for a periodic spin chain. Compared to traditional VMC, weighted VMC reduces the error in the ground state energy by a factor of 2, and it reduces the errors in the local energies away from the mode by large factors of 10 2 –10 4 .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Studies of surface adsorbate electronic structure and femtochemistry at the fundamental length and time scales. Final report

The electronic structure and ultrafast (10-15 s-femtosecond timescale) electron dynamics were investigated for clean and atom/molecule covered metal surfaces. The studies were performed by scanning tunneling microscopy (STM) to measure the structure of adsorbed atoms and molecules on metal surfaces, and to investigate their electronic properties. The electronic structure of the observed molecular networks was calculated by electronic structure theory in collaboration with Prof. Jin Zhao, who is a long-time collaborator, a Professor at the University of Science and Technology of China, and holds an Adjunct Professorship at the University of Pittsburgh. We also investigated the electronic properties of C60 molecules when they are templated by corrugated black phosphorous surfaces. We found unexpected charge delocalization that is enabled by the templating. This research was done in collaboration with Professor Min Feng at the Wuhan University, and who also holds an Adjunct Professorship at the University of Pittsburgh. Moreover, the electronic structure and electron dynamics in metal surfaces were investigated by time-resolved photoemission electron spectroscopy. The focus of ultrafast spectroscopy has been on the plasmonic response of silver surfaces. One direction has been to develop multidimensional (energy, momentum, and time) photoelectron spectroscopy of the coherent response of solid surfaces. This method was applied to study the collective electron excitations known generally as plasmons, which screen optical fields from penetration into metals. Although this collective response has been known for more than 60 years and is used extensively to deposit optical energy into metals, how this happens is poorly known. We investigated the plasmonic response of silver at the point where the dielectric response passes through zero and bulk plasmon is excited by light. We discovered that the plasmon excitation decays by exciting electrons from the Fermi level of a metal, which is contrary to what is believed in the plasmonic science community. This research has been performed in collaboration with Dr. Marcel Reutzel, who was a postdoctoral fellow working on this research at the University of Pittsburgh, and now has a faculty position at the University of Göttingen in Germany. Prof. Branko Gumhalter from the Institute of Physics in Zagreb contributed on the theory of plasmonic decay processes. Furthermore, we investigated the Floquet engineering of electronic bands in metals leading to multiphoton photoemission and above threshold photoemission. Finally, we demonstrated that it is possible to change the electronic structure of metals by application of optical fields. Our studies indicated that this happens on subfemtosecond time scale and could potentially be used in ultrafast information processing and quantum computation. Related document information

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Performance of the rigorous renormalization group for first-order phase transitions and topological phases

Expanding and improving the repertoire of numerical methods for studying quantum lattice models is an ongoing focus in many-body physics. While the density matrix renormalization group (DMRG) has been established as a practically useful algorithm for finding the ground state in one-dimensional systems, a provably efficient and accurate algorithm remained elusive until the introduction of the rigorous renormalization group (RRG) by Landau [Nat. Phys. 11, 566 (2015)1745-247310.1038/nphys3345]. In this paper, we study the accuracy and performance of a numerical implementation of RRG at first-order phase transitions and in symmetry-protected topological phases. Our study is motivated by the question of when RRG might provide a useful complement to the more established DMRG technique. In particular, despite its general utility, DMRG can give unreliable results near first-order phase transitions and in topological phases, since its local update procedure can fail to adequately explore (near-)degenerate manifolds. Furthermore, the rigorous theoretical underpinnings of RRG, meanwhile, suggest that it should not suffer from the same difficulties. We show this optimism is justified, and that RRG indeed determines well-ordered, accurate energies even when DMRG does not. Moreover, our performance analysis indicates that in certain circumstances seeding DMRG with states determined by coarse runs of RRG may provide an advantage over simply performing DMRG.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Comparing the performance of density functionals in describing the adsorption of atoms and small molecules on Ni(111)

Reaction networks calculated using density functional theory (DFT) play a crucial role in understanding and predicting heterogeneous catalyst performance. However, results obtained using different density functionals can vary significantly. Accordingly, it is important to assess the accuracy of available functionals in capturing adsorption energetics derived from single-crystal adsorption microcalorimetry experiments. Here, we evaluate how different density functionals describe the molecular and dissociative adsorption of CH 3 I, the adsorption of CH 3 , I, and H, and the dissociative adsorption of CH 4 on Ni(111). First, we identify energetically-preferred adsorbed states for each adsorbate at three surface coverages and calculate the corresponding enthalpies of adsorption at 160 K. Then, we compare the results to experimental measurements published in the literature. For each adsorption process we find at least one functional that is quantitatively accurate within experimental error, but no functional is accurate (within experimental error) for all adsorption processes studied. Even when assuming an additional DFT error of ±20 kJ/mol on top of the experimental error, only the PBE-D3 and RPBE-D3 functionals are accurate for all considered adsorption systems. Overall, these results indicate that quantitative agreement between density functional theory calculations and experimental measurements is both system- and functional-dependent. Importantly, single-crystal adsorption microcalorimetry experiments will continue to play an important role for benchmarking density functionals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Solution-phase p-type doping of highly enriched semiconducting single-walled carbon nanotubes for thermoelectric thin films

Single-walled carbon nanotubes (SWCNTs) are attractive materials for next-generation energy-harvesting technologies, including thermoelectric generators, due to their tunable opto-electronic properties and high charge carrier mobilities. Controlling the Fermi level within these unique 1D nanomaterials is often afforded by charge transfer interactions between SWCNTs and electron or hole accepting species. Conventional methods to dope SWCNT networks typically involve the diffusion of molecular redox dopant species into solid-state thin films, but solution-phase doping could potentially provide routes and/or benefits for charge carrier transport, scalability, and stability. Here, we develop a methodology for solution-phase doping of polymer-wrapped, highly enriched semiconducting SWCNTs using a p-type charge transfer dopant, F4TCNQ. This allows doped SWCNT inks to be cast into thin films without the need for additional post-deposition doping treatments. We demonstrate that the introduction of the dopant at varying stages of the SWCNT dispersion process impacts the ultimate thermoelectric performance and observe that the dopant alters the polymer selectivity for semiconducting vs metallic SWCNTs. In contrast to dense semiconducting polymer films, where solution-phase doping typically leads to disrupted morphologies and poorer TE performance than solid-state doping, thin films of solution-doped s-SWCNTs perform similarly to their solid-state doped counterparts. Interestingly, our results also suggest that solution-phase F4TCNQ doping leads to fully ionized and dimerized F4TCNQ anions in solid-state films that are not observed in films doped with F4TCNQ after deposition. Our results provide a framework for the application of solution-phase doping to a broad array of high-performance SWCNT-based thermoelectric materials and devices that may require high-throughput deposition techniques.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Effects of pore connectivity and tortuosity on the dynamics of fluids confined in sub-nanometer pores

Dynamical behavior of fluids under nano-pore confinement is studied extensively as it has important implications for several industrial as well as geological processes. Pore network in many porous materials exhibits a varied degree of inter connections. The extent of this pore connectivity may affect the structural and dynamical behavior of the confined fluid. However, studies of fluid confinement addressing these effects systematically are lacking. In this paper we report molecular dynamics simulation studies addressing the effects of pore connectivity on the dynamics of two representative fluids – CO 2 and ethane in silicalite by systematically varying the degree of pore connectivity through selectively blocking some pore space with immobile methane molecules. By selectively turning off the pore spaces in the shape of straight, or tortuous zigzag channels, we also probe the effects of pore tortuosity. In general, pore connectivity is found to facilitate both the translational as well as rotational dynamics of both fluids, while the intermolecular modes of vibration in both fluids remain largely unaffected. The effects of providing connections between a set of straight or zigzag channel-like pores are however more nuanced. Pore tortuosity facilitates the rotational motion, but suppresses the translational motion of CO 2 , while its effects on the rotational and translational motion of ethane are less pronounced. The intermolecular vibrational modes of both fluids shift to higher energies with an increase in the number of tortuous pores. The results reported here provide a detailed molecular level understanding of the effects of pore connectivity on the dynamics of fluids and thus have implications for applications like fluid separation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Navigating Transition-Metal Chemical Space: Artificial Intelligence for First-Principles Design

Conspectus The variability of chemical bonding in open-shell transition-metal complexes not only motivates their study as functional materials and catalysts but also challenges conventional computational modeling tools. Here, tailoring ligand chemistry can alter preferred spin or oxidation states as well as electronic structure properties and reactivity, creating vast regions of chemical space to explore when designing new materials atom by atom. Although first-principles density functional theory (DFT) remains the workhorse of computational chemistry in mechanism deduction and property prediction, it is of limited use here. DFT is both far too computationally costly for widespread exploration of transition-metal chemical space and also prone to inaccuracies that limit its predictive performance for localized d electrons in transition-metal complexes. These challenges starkly contrast with the well-trodden regions of small-organic-molecule chemical space, where the analytical forms of molecular mechanics force fields and semiempirical theories have for decades accelerated the discovery of new molecules, accurate DFT functional performance has been demonstrated, and gold-standard methods from correlated wavefunction theory can predict experimental results to chemical accuracy. The combined promise of transition-metal chemical space exploration and lack of established tools has mandated a distinct approach. In this Account, we outline the path we charted in exploration of transition-metal chemical space starting from the first machine learning (ML) models (i.e., artificial neural network and kernel ridge regression) and representations for the prediction of open-shell transition-metal complex properties. The distinct importance of the immediate coordination environment of the metal center as well as the lack of low-level methods to accurately predict structural properties in this coordination environment first motivated and then benefited from these ML models and representations. Once developed, the recipe for prediction of geometric, spin state, and redox potential properties was straightforwardly extended to a diverse range of other properties, including in catalysis, computational “feasibility”, and the gas separation properties of periodic metal–organic frameworks. Interpretation of selected features most important for model prediction revealed new ways to encapsulate design rules and confirmed that models were robustly mapping essential structure–property relationships. Encountering the special challenge of ensuring that good model performance could generalize to new discovery targets motivated investigation of how to best carry out model uncertainty quantification. Distance-based approaches, whether in model latent space or in carefully engineered feature space, provided intuitive measures of the domain of applicability. With all of these pieces together, ML can be harnessed as an engine to tackle the large-scale exploration of transition-metal chemical space needed to satisfy multiple objectives using efficient global optimization methods. In practical terms, bringing these artificial intelligence tools to bear on the problems of transition-metal chemical space exploration has resulted in ML-model assessments of large, multimillion compound spaces in minutes and validated new design leads in weeks instead of decades.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Physics-Constrained Deep Learning Model for Simulating Multiphase Flow in 3D Heterogeneous Porous Media

Physics-based simulators for multiphase flow in porous media emulate nonlinear processes with coupled physics, and usually require extensive computational resources for software development, maintenance and simulation execution. As a result, a huge demand exists for fast modeling of coupled processes in a wide range of subsurface applications including geological sequestration, hydrocarbon recovery and geothermal energy extraction. In this work, an efficient physics-constrained deep learning model is developed for solving multiphase flow in 3-Dimensional (3D) heterogeneous porous media. The model fully leverages the spatial topology predictive capability of convolutional neural networks, specifically U-Net with successive contracting and expansive steps, and is coupled with an efficient continuity-based smoother to predict flow responses that need spatial continuity. Furthermore, the transient regions are penalized to steer the training process such that the model can accurately capture flow in these regions. The model takes inputs including properties of porous media, fluid properties and well controls, and predicts the temporal-spatial evolution of the state variables (pressure and saturation). While maintaining the continuity of fluid flow, the 3D spatial domain is decomposed into 2D images for reducing training cost, and the decomposition results in an increased number of training data samples and better training efficiency. Additionally, a surrogate model is separately constructed as a postprocessor to calculate well flow rate based on the predictions of state variables from the deep learning model. We use the example of CO 2 injection into saline aquifers, and apply the physics-constrained deep learning model that is trained from physics-based simulation data and emulates the physics process. The model performs prediction with a speedup of ~ 1400 times compared to physics-based simulations, and the average temporal errors of predicted pressure and saturation plumes are 0.27% and 0.099% respectively. Furthermore, water production rate is efficiently predicted by a surrogate model for well flow rate, with a mean error less than 5%. Therefore, with its unique scheme to cope with the fidelity in fluid flow in porous media, the physics-constrained deep learning model can become an efficient predictive model for computationally demanding inverse problems or other coupled processes.

58 GEOSCIENCES↗

VLBI Data Processing on Coronal Radio-sounding Experiments of Mars Express

The ESA’s Mars Express solar corona experiments were performed at two solar conjunctions in the years 2015 and 2017 by a number of radio telescopes in the European VLBI Network. This paper presents the methods to measure the frequency and phase fluctuations of the spacecraft radio signal, and the applications to study the characteristics of the plasma turbulence effects on the signal at a single station and at multiple stations via cross correlation. The power spectra of the frequency fluctuations observed between 4.9 and 76.3 R {sub s} have a power-law shape close to a Kolmogorov spectrum over the frequency interval ν {sub lo} < ν < ν {sub up}, where the nominal value of ν {sub lo} is set to 3 mHz and ν {sub up} is in the range of 0.03–0.15 Hz. The rms of the frequency fluctuations is presented as a function of the heliocentric distance. Furthermore, we analyze the variations of the electron column density fluctuations at solar offsets 4.9 R {sub s} and 9.9 R {sub s} and the cross-correlation products between the VLBI stations. The power density of the differential fluctuations between different stations decreases at ν < 0.01 Hz. Finally, the fast flow speeds of solar wind >700 km s{sup -1} are derived from the cross correlation of frequency fluctuations at ν < 0.01 Hz. The fast flow speeds of solar wind correspond to the high heliolatitude of the coronal region that the radio rays passed. The VLBI observations and analysis methods can be used to study the electron column density fluctuations and the turbulence at multiple spatial points in the inner solar wind by providing multiple lines of sight between the Earth and the spacecraft.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Final Reports of the 2021 Los Alamos National Laboratory Computational Physics Student Summer Workshop

Since 2011, the Los Alamos National Laboratory Computational Physics Student Summer Workshop has been bringing together a highly talented and diverse group of students every summer. Students work in teams of two, alongside typically two mentors, on research projects reflecting a broad range of topics within computational physics. In addition, students attend a series of lectures on topics within computational physics, facility tours, and networking events. The program lasts ten weeks, with this year’s workshop running from June 7 to August 13. At the end of the summer, students give a final presentation, along with a written report. Those reports are what make up the remaining sections of this document. Admission to the workshop is by a competitive process, with the mentors forming the selection committee. One of the important accomplishments of the workshop has been to create a student pipeline from diverse schools that sometimes are not normally tapped by LANL recruiting. Many workshop students maintain a continuing relationship with LANL, returning as student interns, post-doctoral researchers, and staff members. Additionally, workshop alumni act as ambassadors for LANL. The result is a wider awareness both of LANL as a potential employer, and of the technical work that happens at LANL. This year, the workshop was once again in an off-site, virtual format. Students worked on LANL virtual desktop systems remotely, also accessing LANL HPC resources. In order to facilitate communication, student were given accounts on both Webex, a video teleconferencing platform, and Mattermost, an online team collaboration and chat platform, similar to Slack. Daily communication between students and mentors was primarily on Mattermost, with Webex conferencing as needed. The lectures were all done on Webex. Given the difficulty of the virtual format, and a concern that students might have video teleconferencing burn-out, all lectures were optional this year. In spite of this, the attendance was generally high. Lecturers were asked to try to move to a more high-level, ”What is it?,” format. Once again, the students did a tremendous job. Over the course of ten weeks, they did important research across a staggering array of disciplines. The following pages contain the final report for each team’s research efforts. Enjoy!

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Final Reports of the 2020 Los Alamos National Laboratory Computational Physics Student Summer Workshop

For the past ten years, the workshop has been bringing a highly talented and diverse group of student every summer. Students work in teams of two, alongside typically two mentors, on research projects reflecting a broad range of topics within computational physics. In addition, students attend a series of lectures on topics within computational physics, facility tours, and networking events. The program lasts ten weeks, with this year’s workshop running from June 8 to August 14. At the end of the summer, students give a final presentation, along with a written report. Those reports are what make up the remaining sections of this document. Admission to the workshop is by a competitive process, with the mentors forming the selection committee. One of the important accomplishments of the workshop has been to create a student pipeline from diverse schools that sometimes are not normally tapped by LANL recruiting. Many workshop students maintain a continuing relationship with LANL, returning as students interns, post-doctoral researchers, and staff members. Additionally, workshop alumni act as ambassadors for LANL. The result is a wider awareness both of LANL as a potential employer, and of the technical work that happens at LANL. This year, the workshop format was changed in several ways, in order to accommodate the off-site, virtual format. Students worked on LANL virtual desktop systems remotely, also accessing LANL HPC resources. In order to facilitate communication, student were given accounts on both Webex, a video teleconferencing platform, and Mattermost, an online team collaboration and chat platform, similar to Slack. Daily communication between students and mentors was primarily on Mattermost, with Webex conferencing as needed. The lectures were all done on Webex. Given the difficulty of the virtual format, and a concern that students might have video teleconferencing burn-out after an academic semester largely moved to that format, all lectures were optional this year. In spite of this, the attendance was generally high. Lecturers were asked to try to move to a more high-level, ”What is it?,” format. Once again, the students did a tremendous job. Over the course of ten weeks, they did important research across a staggering array of disciplines. The following pages contain the final report for each team’s research efforts. We hope you will find reading them as exciting as it was for us to produce them.

36 MATERIALS SCIENCE↗

Multi-fidelity learning for interatomic potentials: low-level forces and high-level energies are all you need

The promise of machine learning interatomic potentials (MLIPs) has led to an abundance of public quantum mechanical (QM) training datasets. The quality of an MLIP is directly limited by the accuracy of the energies and atomic forces in the training dataset. Unfortunately, most of these datasets are computed with relatively low-accuracy QM methods, e.g. density functional theory with a moderate basis set. Due to the increased computational cost of more accurate QM methods, e.g. coupled-cluster theory with a complete basis set (CBS) extrapolation, most high-accuracy datasets are much smaller and often do not contain atomic forces. The lack of high-accuracy atomic forces is quite troubling, as training with force data greatly improves the stability and quality of the MLIP compared to training to energy alone. Because most datasets are computed with a unique level of theory, traditional single-fidelity (SF) learning is not capable of leveraging the vast amounts of published QM data. In this study, we apply multi-fidelity learning (MFL) to train an MLIP to multiple QM datasets of different levels of accuracy, i.e. levels of fidelity. Specifically, we perform three test cases to demonstrate that MFL with both low-level forces and high-level energies yields an extremely accurate MLIP—far more accurate than a SF MLIP trained solely to high-level energies and almost as accurate as a SF MLIP trained directly to high-level energies and forces. Therefore, MFL greatly alleviates the need for generating large and expensive datasets containing high-accuracy atomic forces and allows for more effective training to existing high-accuracy energy-only datasets. Indeed, low-accuracy atomic forces and high-accuracy energies are all that are needed to achieve a high-accuracy MLIP with MFL.

36 MATERIALS SCIENCE↗

HydraGNN_Predictive_GFM_2024 - Ensemble of predictive graph foundation models for ground state atomistic materials modeling

We provide the ensemble of fifteen pre-trained graph foundation models (GFMs) for atomistic materials modeling applications. Each one of the fifteen GFMs has been trained on five open-source datasets that (once aggregated) amount to over 154 million atomistic structures, which cover over two-thirds of the natural elements of the periodic table and that comprises a broad set of organic and inorganic compounds. This vast set of atomistic structures comprises ground state configurations that are dynamically stable (i.e., equilibrated structures with atomic forces approximately close to zero values) as well as dynamically unstable structures (i.e., non-equilibrium structures with non-negligible non-zero values of atomic forces). The ensemble of datasets aggregated does NOT include excited states. The datasets have been curated to remove atomistic structures with spectral norm of the force tensor above 100 eV/angstrom. Moreover, a linear term of the energy was computed for each dataset using a linear regression model that uses the chemical concentration of each natural element as regressor. The linear term predicted by the linear regression model has been subtracted from each original energy value to perform a re-alignment of the energy values across different electronic structures approximation theories performed to generate the diverse multi-source, multi-fidelity datasets. The folder "ADIOS_files" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "ADIOS_files" directory contains 6 sub-directories named as follows: - ANI1x-v3.bp - MPTrj-v3.bp - OC2020-20M-v3.bp - OC2020-v3.bp - OC2022-v3.bp - qm7x-v3.bp Each sub-directory contains the pre-processed datasets converted in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used to the development, training, and performance testing of the ensemble go predictive graph foundation models. Each GFM was developed using HydraGNN (https://github.com/ORNL/HydraGNN) as underlying graph neural network (GNN) architecture. The multi-task learning (MTL) capability of HydraGNN was used to simultaneously train the GFMs on labeled values for direct predictions of energy (a total system property of an atomistic structure that measures the chemical stability) and atomic forces (an atomic level property of an atomistic structure that measures the dynamical stability). The hyper parameters of the GFM have been tuned using scalable hyperparameter optimization (HPO) algorithms implemented in the software DeepHyper (https://github.com/deephyper/deephyper). The pre-training of each HPO trial was performed using distributed data parallelism (DDP) to scale the training across 128 compute nodes of the exascale OLCF supercomputer Frontier. Each HPO trial was trained only for 10 epochs and an early stopping was performed to avoid wasting significant computational resources on GNN architectures that were clearly underperforming. For each HPO trial, the 'omnistat' tool developed by (AMD Research - Advanced Micro Device) was used to measure the total energy consumption in kWh. The ensemble of GFMs was obtained by selecting the fifteen best performing HPO trials. Four models have been selected for their clear advantage in accuracy, and these are the GFMs with IDs 229, 156, 147, 260. Additional eleven models have been selected based on judicious balance between accuracy and energy consumption needed for training, and these are the GFMs with IDs 165, 78, 137, 1, 175, 171, 181, 67, 179, 167, 351. Each selected GFM of the ensemble was continued to cumulate a total of at most 30 epochs. In some cases, the total number of epochs actually performed was les than 30 due to two combined factors: (1) the size of the GFM (i.e., the number of model parameters to train) and (2) the total wall-clock time for which the computational resources could be allocated on OLCF-Frontier. The "Ensemble_of_models" directory contains 15 sub-directories named as follows: - gfm_0.229 - gfm_0.156 - gfm_0.147 - gfm_0.260 - gfm_0.165 - gfm_0.78 - gfm_0.137 - gfm_0.1 - gfm_0.175 - gfm_0.171 - gfm_0.181 - gfm_0.67 - gfm_0.179 - gfm_0.167 - gfm_0.351 Each one of these sub-directories refers to one of the fifteen HPO trials that have been selected to continue the pre-training with at most 30 epochs. With each sub-directory associated with a specific HPO trial, the following files can be found: - config.json: file for argument parsing to develop and train an HydraGNN architecture - gfm_0.ID_epoch_N.pk: file with model parameters for HPO ID trial after N epochs of training The ensemble of fifteen GFM architectures was used for (1) ensemble averaging to stabilize the predictions of energy and atomic forces after pre-training for post-processing analysis and (2) ensemble uncertainty quantification (UQ). The code used to develop, pre-train, and load the pre-trained models for post-processing analysis is available on the ORNL-GitHub at the following link: https://github.com/ORNL/HydraGNN/tree/Predictive_GFM_2024

36 MATERIALS SCIENCE↗