Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantum mechanical systems and processes”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Interligand Coupling Drives Fast Triplet Energy Transfer Routes in PbS/Tetracene Quantum Dot Hybrids

The binding of photoactive organic ligands to inorganic quantum dots (QDs) creates a versatile hybrid architecture that allows access to photophysical processes such as efficient triplet exciton generation with near-infrared radiation. Here we report the subnanosecond generation of a hybrid triplet state with mixed ligand-QD character by replacing native oleate ligands on small PbS QDs with 5,12-tetracenepropiolic acid, a bifunctional ligand with two carboxylic acids that tends to lie face-on with the QD surface at low loadings. The face-on geometry engenders a regime of strong electronic coupling that is evident in steady-state absorption and hastens triplet energy flow by several orders of magnitude compared with more typical tetracene-based ligands exhibiting weak coupling. We further determined via Fourier transform infrared (FTIR) and supported by density functional theory (DFT)-based geometry optimizations that high ligand loading causes a shift in QD-ligand mutual disposition toward an edge-on geometry that instigates the formation of intermolecular excited states characterized by triplet excimer-like features in photoluminescence and transient absorption. Our results demonstrate the ability to control strongly coupled ligand-QD systems toward ultrafast generation of photophysically relevant species such as triplets that are valuable for photon upconversion and catalysis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

First-Principles Simulation of Beam-Induced Processes Underlying Atomic Manipulation in Electron Microscopes

The development of experimental methods and apparatuses capable of promoting atomically precise material manipulations holds great promise for realizing the ultimate limit of feature miniaturization in materials and devices. The ability to modify materials atom by atom is anticipated to usher in new technologies in areas as diverse as separation science, medicine, and quantum information science. Historically, scanning probe-based techniques have been the most prominent approaches in this space. However, these methods are best suited for the manipulation of surface-exposed regions of materials, as the strong perturbations required for bond scission are delivered most effectively to atoms in the near-proximity to the scanning probe. In contrast, convergent electron beams with energies tuned slightly below the threshold for inducing irreversible knock-on damage have recently been employed (within scanning transmission electron microscopy) to promote atomic-scale bond rearrangements in various beam-stable solids. Currently, however, the efficiency and selectivity of beam-induced atomic manipulation processes with focused electron beams are such that long irradiation times are required to induce a desired atomic rearrangement. With a better understanding of the underlying physics dictating the outcome of a given irradiation event, methods can be devised to improve the efficiency of these techniques so that their promise can be fully realized through widespread adoption.To this end, this Account details our recent efforts to develop and apply tractable first-principles simulation approaches for studying the response of materials to electric beam-like external electric potentials applied in real space. We briefly review the concepts and capabilities in the area of atomically precise materials manipulation and review the early demonstrations of accomplishments in this area, focusing on studies using scanned convergent electron beam probes in particular. We expound upon the depth of the challenge and identify critical shortcomings of theoretical methods that have previously been employed in the simulation of beam-induced processes. We then describe the computational methods that we have generalized from the concepts and tools most commonly applied to the study of molecular photochemistry and how our adaptations of these methods can be employed to capture the relevant dynamical phenomena for beam-induced processes ranging from the initial electron scattering to the ensuing multistate reactions. Here, we contextualize these methods within the current state of the art in this area, which has historically focused primarily on the simulation of inelastic image formation in the electron microscope for the purpose of interpreting the results of quantitative electron microscopy experiments. We demonstrate that the spatial distribution of state-specific excitation rates due to the presence of an external (probe) electric charge is inhomogeneous, such that irradiation at particular locations in materials can favor specific electronic transitions (and disallow others). In addition to the potential for excited-state reaction pathways to be accessed through the initial inelastic scattering of the tightly focused electron beam from the targeted atoms, we also identify favorable conditions for the electronically nonadiabatic evolution of the highly vibrationally excited system to open complex multistate reaction pathways. Implications of the early results for understanding the mechanisms and potential routes to improved efficiency and selectivity in beam-induced reactions are discussed. We conclude with a summary of the current state of theory and modeling capabilities in this area and provide our perspective on future directions for theoretical and experimental developments that we view as crucial to advancing the use of convergent electron beams in mode-specific, atomically precise platforms for direct-write materials modifications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Charge Transfer and Spin Dynamics in a Zinc Porphyrin Donor Covalently Linked to One or Two Naphthalenediimide Acceptors

We report quantum coherence effects on charge transfer and spin dynamics in a system having two degenerate electron acceptors are studied using a zinc 5,10,15-tri(n-pentyI)-20-phenylporphyrin (ZnP) electron donor covalently linked to either one or two naphthalene-1,8:4,5-bis(dicarboximide) (NDI) electron acceptors using an anthracene (An) spacer, ZnP-An-NDI (1) and ZnP-An-NDI 2 (2), respectively. Following photoexcitation of 1 and 2 in toluene at 295 K, femtosecond transient absorption spectroscopy shows that the electron transfer (ET) rate constant for 2 is about three times larger than that of 1, which can be accounted for by the statistical nature of incoherent ET as well as the electron couplings for the charge separation reactions. In contrast, the rate constant for charge recombination (CR) of 1 is about 25% faster than that of 2. Using femtosecond transient infrared spectroscopy and theoretical analysis, we find that the electron on NDI 2 ∙- in 2 localizes onto one of the two NDIs prior to CR, thus precluding electronically coherent CR from NDI 2 ∙- . Conversely, CR in both 1 and 2 is spin coherent as indicated by the observation of a resonance in the 3* ZnP yield following CR as a function of applied magnetic field, giving spin-spin exchange interaction energies of 2J = 210 and 236 mT, respectively, where the line width of the resonance for 2 is greater than 1. These data show that while CR is a spin-coherent process, incoherent hopping of the electron between the two NDIs in 2, consistent with the lack of delocalization noted above, results in greater spin decoherence in 2 relative to 1.

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

Metal-Dependent Photodissociation of Hydrazone Photoswitches from Rare-Earth Complexes

Rare-earth element separation processes often rely on a small decrease in ionic radii along the series of elements. Separation processes based on the distinct optical properties of REs remain less explored, although photochemical methods may offer a viable alternative. Accurate selection of the synthetic precursors of a photoswitchable acylhydrazonic ligand led to a system that could quantitatively isomerize (E−to−Z) upon irradiation with commercial LED lights. Coordination of the photoswitch with RE III nitrates (RE = La−Lu except Pm and Y) resulted in the retention of the photoswitching properties observed in solution. The lower binding affinity of the generated Z−isomer with RE III ions yielded the dissociation of the complexes upon irradiation (photodissociation) with the release of RE−nitrates in solution. The rate of the reaction was found to be dependent on the optical properties of the RE III ions, with nonemissive complexes (no 4f excited states) dissociating faster than emissive ones (having accessible 4f excited states). A thorough solid-state characterization of the complexes was performed by using crystallographic and photochemical methods. Ultimately, the accessibility of the 4f excited states of the metals following light irradiation and excitation of the ligand led to a decrease in the rate of the reaction due to quenching of the ligand excited state. These results demonstrate that the direct modulation of the metal coordination environment, combined with the metal-dependent reaction rate, could provide a strategy for the development of RE−separation processes based on differences in their optical properties.

Ions↗
DOE OSTI2024

Circularly Polarized Attosecond Pulses Enabled by an Azimuthal Phase and Polarization Grating

High-harmonic generation (HHG) is an extreme nonlinear optical process that can map the properties of an infrared driving laser beam onto short wavelength attosecond pulse trains. However, current techniques for generating circularly polarized high harmonics for probing magnetic materials and chiral systems have limitations: two-color collinear counter-rotating driving lasers result in a low cutoff photon energy, while single-color non-collinear counter-rotating schemes suffer from low conversion efficiency. In this work, we generate circularly polarized attosecond pulse trains by using a structured laser driver which has a rotating polarization and phase grating along the azimuthal coordinate. Furthermore, our experimental and numerical results demonstrate the production of left and right circularly polarized harmonics, which naturally separate upon propagation. Our approach uses a single laser color in a collinear geometry, that can be scaled for high efficiency. Simulations show this scheme can extend into the soft x-ray region when driven by mid-infrared driving lasers, while preserving the same high phase-matching cutoff photon energy as for linearly-polarized high harmonics.

attosecond science↗
DOE OSTI2022

Ultrafast carrier dynamics in terahertz photoconductors and photomixers: beyond short-carrier-lifetime semiconductors

Abstract Efficient terahertz generation and detection are a key prerequisite for high performance terahertz systems. Major advancements in realizing efficient terahertz emitters and detectors were enabled through photonics-driven semiconductor devices, thanks to the extremely wide bandwidth available at optical frequencies. Through the efficient generation and ultrafast transport of charge carriers within a photo-absorbing semiconductor material, terahertz frequency components are created from the mixing products of the optical frequency components that drive the terahertz device – a process usually referred to as photomixing. The created terahertz frequency components, which are in the physical form of oscillating carrier concentrations, can feed a terahertz antenna and get radiated in case of a terahertz emitter, or mix with an incoming terahertz wave to down-convert to DC or to a low frequency photocurrent in case of a terahertz detector. Realizing terahertz photoconductors typically relies on short-carrier-lifetime semiconductors as the photo-absorbing material, where photocarriers are quickly trapped within one picosecond or less after generation, leading to ultrafast carrier dynamics that facilitates high-frequency device operation. However, while enabling broadband operation, a sub-picosecond lifetime of the photocarriers results in a substantial loss of photoconductive gain and optical responsivity. In addition, growth of short-carrier-lifetime semiconductors in many cases relies on the use of rare elements and non-standard processes with limited accessibility. Therefore, there is a strong motivation to explore and develop alternative techniques for realizing terahertz photomixers that do not rely on these defect-introduced short-carrier-lifetime semiconductors. This review will provide an overview of several promising approaches to realize terahertz emitters and detectors without short-carrier-lifetime semiconductors. These novel approaches utilize p-i-n diode junctions, plasmonic nanostructures, ultrafast spintronics, and low-dimensional materials to offer ultrafast carrier response. These innovative directions have great potentials for extending the applicability and accessibility of the terahertz spectrum for a wide range of applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

ORNL_AISD_NiPt_108atoms

This dataset describes the nickel-platinum (NiPt) solid solution binary alloy, where the two constituent elements nickel (Ni) and platinum (Pt) are randomly placed on the face centered cubic (FCC) crystal structure, with the lattice constant of 3.840 angstroms. The dataset comprises data for crystal structures with 108 atoms with 1,900 configurations. The data set was generated for concentrations ranging from 0at% of Pt to 100at% of Pt in the NiPt binary system, with increasing the concentration of Pt in the system every 5at%. For each one of the chemical compositions, 100 random configurations were generated, each with a different random seed. Each of the output files contains the mass, type, atomic coordinates, energy per atom, and forces in x, y, and z directions respectively. For each atomic configuration, the output was collected every 150 steps during the minimization stage and every 1000 steps during the replica exchange stage. Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) [1], which is a molecular dynamics code, was used to generate data for NiPt alloy. The simulation used the interatomic potential for NiPt binary system 'MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001' [3] from the OpenKIM library (Open Knowledgebase of Interatomic Models) [2]. This potential was developed based on the second nearest-neighbor modified embedded-atom method (2NN MEAM). The simulation process begins with the generation of the random NiPt structure and follows with the short minimization and replica exchange simulation. The minimization procedure adjusts atomic coordinates and performs energy minimization, which typically leads to a local potential energy minimum. The method used for the minimization was the conjugate gradient algorithm. A short replica exchange (parallel tempering) simulation involves four replicas (ensembles) of a system and follows the minimization stage. Multiple snapshots of the configuration were collected during the minimization and replica exchange stages. NiPt alloy is interesting due to its magnetic and charge transfer properties [4]. The data is provided in a compressed zipped folders atoms108.zip. The zipped folder contains the data structured in the following way: - Ni_ground_state.cfg --> atomic configuration for the pure nickel - Pt_ground_state.cfg --> atomic configuration for the pure platinum - Pt#_filtered --> folders containing atomic configurations for #at% concentration of platinum. The folder contains 100 atomic configurations, each saved in a subfolder - Each subfolder named config* is associated with a specific atomic configuration. Each of these subfolders contains files with .cfg format, corresponding to outputs for each atomic configuration The total number of atomic configurations contained in atoms108.zip is 66,132. This dataset is an extension to the dataset ORNL_AISD_NiPt [5] that has been previously released with crystal structures of 256 atoms, 864 atoms, and 2,048 atoms, with the same methodology for data collection. References [1] https://www.lammps.org/ [2] https://openkim.org/ [3] https://openkim.org/id/MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001 [4] El-Gendy, Ahmed A. and Hampel, Silke and Büchner, Bernd and Klingeler, Rüdiger, Tuneable magnetic properties of carbon-shielded NiPt-nanoalloys, RSC Adv., volume 6, issue 57, pages 52427-52433, 2016, The Royal Society of Chemistry, doi:10.1039/C6RA05910D [5] M. Karabin, M. Lupo Pasini, and M. Eisenbach. ORNL_AISD_NiPt. United States: N. p., 2023. Web. doi:10.13139/OLCF/1958172.

36 MATERIALS SCIENCE↗
DOE OSTI2023

ORNL_AISD_NiPt

This dataset describes the nickel-platinum (NiPt) solid solution binary alloy, where the two constituent elements nickel (Ni) and platinum (Pt) are randomly placed on the face centered cubic (FCC) crystal structure, with the lattice constant of 3.840 angstroms. The dataset comprises data for three different sizes of the crystal structure: 256 atoms, 864 atoms, and 2,048 atoms, each of which contains 1900 configurations. For each size of the crystal structure, the data set was generated for concentrations ranging from 0at% of Pt to 100at% of Pt in the NiPt binary system, with increasing the concentration of Pt in the system every 5at%. For each one of the chemical compositions, 100 random configurations were generated, each with a different random seed. Each of the output files contains the mass, type, atomic coordinates, energy per atom, and forces in x, y, and z directions respectively. For each atomic configuration, the output was collected every 150 steps during the minimization stage and every 1000 steps during the replica exchange stage. Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) [1], which is a molecular dynamics code, was used to generate data for NiPt alloy. The simulation used the interatomic potential for NiPt binary system MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001 [3] from the OpenKIM library (Open Knowledgebase of Interatomic Models) [2]. This potential was developed based on the second nearest-neighbor modified embedded-atom method (2NN MEAM). The simulation process begins with the generation of the random NiPt structure and follows with the short minimization and replica exchange simulation. The minimization procedure adjusts atomic coordinates and performs energy minimization, which typically leads to a local potential energy minimum. The method used for the minimization was the conjugate gradient algorithm. A short replica exchange (parallel tempering) simulation involves four replicas (ensembles) of a system and follows the minimization stage. Multiple snapshots of the configuration were collected during the minimization and replica exchange stages. NiPt alloy is interesting due to its magnetic and charge transfer properties [4]. The data is provided in three compressed zipped folders: atoms256.zip, atoms864.zip, atoms2048.zip Each zipped folder contains the data that describes crystals of size 256 atoms, 864 atoms, and 2,048 atoms respectively. Each one of the three zipped folders contains the data structured in the following way: -Ni_ground_state.cfg --> atomic configuration for the pure nickel -Pt_ground_state.cfg --> atomic configuration for the pure platinum -Pt#_filtered --> folders containing atomic configurations for #at% concentration of platinum. The folder contains 100 atomic configurations, each saved in a subfolder. Each subfolder named config* is associated with a specific atomic configuration. Each of these subfolders contains files with .cfg format, corresponding to outputs for each atomic configuration The total number of atomic configurations contained in atoms256.zip is 65,046. The total number of atomic configurations contained in atoms864.zip is 63,936. The total number of atomic configurations contained in atoms2048.zip is 61,997. The total number of atomic configurations spanned by the entire dataset is 190,979. References [1] https://www.lammps.org/ [2] https://openkim.org/ [3] https://openkim.org/id/MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001 [4] El-Gendy, Ahmed A. and Hampel, Silke and Büccchner, Bernd and Klingeler, Rüdiger, Tuneable magnetic properties of carbon-shielded NiPt-nanoalloys, RSC Adv., volume 6, issue 57, pages 52427-52433, 2016, The Royal Society of Chemistry, doi:10.1039/C6RA05910D

36 MATERIALS SCIENCE↗
DOE OSTI2021

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↗
DOE OSTI2024

Project DarkStar: Vision for LLNL in 2030

DarkStar was a Strategic Initiative (FY2021-FY2024) to investigate applications of Artificial Intelligence (AI) and Machine Learning (ML) to scientific problems of complex hydrodynamics, shockwave physics and energetic materials. The research focused on physics and engineering design as a process that can be tremendously accelerated through merging AI with advanced physics simulation on exascale-class platforms, and to experimentally validate this revolutionary new approach through dynamic materials campaigns. A central thread of scientific inquiry was in the application of AI to enable human understanding of how to control hydrodynamic instability (which has impacts to areas such as inertial confinement fusion) via engineering features and time-dependent sources. Motivated by an unfinished line of research started by Dr. Johnny von Neumann, AI-enabled simulation approaches were developed that allowed DarkStar researchers to uncover several ground-breaking discoveries regarding hydrodynamic instability, including how to completely suppress Richtmyer-Meshkov instability (RMI). These S&T discoveries, along with other advances, have shown the way for an entirely new approach to time-dependent problems known as inverse design – the idea that complex systems can be developed directly from a final state that is to be achieved and resolve the initial design via satisfying several constraints simultaneously via AI/ML. Through experimental campaigns conducted across a wide range of facilities in the NNSA complex (the High Explosive Application Facility at LLNL, the Dynamic Compression Sector/Advanced Photon Source at Argonne National Lab, and Special Technologies Laboratory at MSTS) the radical new AI/ML approach to engineering complex material dynamics was verified, establishing a new field of study within the realm of shock physics. As advanced manufacturing capabilities continue to develop, the great importance of inverse design as a means to apply that technology effectively for NNSA missions will feature prominently over this decade. DarkStar has positioned NNSA as a world-leader in this newly emerging cross-disciplinary area of AI methods for advanced physics simulation and pioneered multiple novel approaches that have enabled the broader scientific community. By allowing us to see past the horizon, to 2030 and beyond, DarkStar has illuminated the vast potential of AI/ML to impact a wide range of new national security missions and, consequently, multiple areas of further research have already emerged across the NNSA and DOD complex.

42 ENGINEERING↗
DOE OSTI2020

NWChem: Past, present, and future

Specialized computational chemistry packages have permanently reshaped the landscape of chemical sciences by providing tools to support and guide the experimental effort and for prediction of chemical and materials properties. In this regard, a special role has been played by electronic structure packages where complex chemical and materials processes can be modeled using first-principle-driven methodologies. Over the last few decades, the rapid development of computing technologies and a tremendous increase in computational power has offered a unique chance to study complex chemical transformations using sophisticated and predictive many-body techniques to describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, a critical role has been played by novel parallel algorithms capable of taking advantage of computational resources to address polynomial scaling of electronic structure methods. NWChem was among the first electronic structure codes that focused on delivering scalable performance for electronic structure simulations. Herein, we briefly review the NWChem suite of computational codes including its history, design principles, parallel tools, current capabilities, outreach and outlook.

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

Importance of Electron Correlation on the Geometry and Electronic Structure of [2Fe–2S] Systems: A Benchmark Study of the [Fe 2 S 2 (SCH 3 ) 4 ] 2–,3–,4– , [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S- p -tol) 4 ] 2– , and [Fe 2 S 2 (S- o -xyl) 4 ] 2– Complexes

Iron–sulfur clusters are crucial for biological electron transport and catalysis. Obtaining accurate geometries, energetics, manifolds of their excited electronic states, and reduction energies is important to understand their role in these processes. Here, using a [2Fe–2S] model complex with Fe II and Fe III oxidation states, which leads to different charges, i.e., [Fe 2 S 2 (SMe) 4 ] 2–,3–,4– , we benchmarked a variety of computational methodologies ranging from density functional theory (DFT) to post-Hartree–Fock methods, including complete active space self-consistent field (CASSCF), multireference configuration interaction, the second-order N-electron valence state perturbation theory (NEVPT2), and the linearized integrand approximation of adiabatic connection (AC0) approaches. Additionally, we studied three experimentally well-characterized complexes, [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S-o-tol) 4 ] 2– , and [Fe 2 S 2 (S-o-xyl) 4 ] 2– , via DFT methods. We conclude that the dynamic electron correlation is important for accurately predicting the geometry of these complexes. Broken symmetry (BS) DFT correctly predicts experimental geometries of low-spin multiplicity, while CASSCF does not. However, BS-DFT significantly overestimates the difference between the low- and high-spin electronic states for a given oxidation state. At the same time, CASSCF underestimates it but provides relative energies closer to the reference NEVPT2 results. Finally, AC0 provides energetics of NEVPT2 quality with the additional advantage of being able to use large CASSCF sizes. NEVPT2 gives the best estimates of the Fe III /Fe III → Fe II /Fe III (4.27 eV) and Fe II /F III → Fe II /F II (7.72 eV) reduction energies. The results provide insight into the electronic structure of these complexes and assist in the understanding of their physical properties.

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

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↗
DOE OSTI2025

Schrödinger cat states of a nuclear spin qudit in silicon

High-dimensional quantum systems are a valuable resource for quantum information processing. They can be used to encode error-correctable logical qubits, which has been demonstrated using continuous-variable states in microwave cavities or the motional modes of trapped ions. For example, high-dimensional systems can be used to realize ‘Schrödinger cat’ states, which are superpositions of widely displaced coherent states that can be used to illustrate quantum effects at large scales. Recent proposals have suggested encoding qubits in high-spin atomic nuclei, which are finite-dimensional systems that can host hardware-efficient versions of continuous-variable codes. Here, in this study, we demonstrate the creation and manipulation of Schrödinger cat states using the spin-7/2 nucleus of an antimony atom embedded in a silicon nanoelectronic device. We use a multi-frequency control scheme to produce spin rotations that preserve the symmetry of the qudit, and we constitute logical Pauli operations for qubits encoded in the Schrödinger cat states. Our work demonstrates the ability to prepare and control non-classical resource states, which is a prerequisite for applications in quantum information processing and quantum error correction, using our scalable, manufacturable semiconductor platform.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Characterizing non-Markovian and coherent errors in quantum simulation

Quantum simulation of many-body systems, particularly using ultracold atoms and trapped ions, presents a unique form of quantum control—it is a direct implementation of a multi-qubit gate generated by the Hamiltonian. As a consequence, it also faces a unique challenge in terms of benchmarking, because the well-established gate benchmarking techniques are unsuitable for this form of quantum control. Here we show that the symmetries of the target many-body Hamiltonian can be used not only to benchmark but to characterize experimental errors in the quantum simulation. We use our results to develop protocols to characterize these errors, which can be implemented using state-of-the-art technology. We consider two forms of errors: (i) unitary errors arising out of systematic errors in the applied Hamiltonian and (ii) canonical non-Markovian errors arising out of random shot-to-shot fluctuations in the applied Hamiltonian. We show that the dynamics of the expectation value of the target Hamiltonian itself, which is ideally constant in time, can be used to characterize these errors. In the presence of errors, the expectation value of the target Hamiltonian shows a characteristic thermalization dynamics, when it satisfies the operator thermalization hypothesis (OTH). That is, an oscillation in the short time followed by relaxation to a steady-state value in the long time limit. We show that while the steady-state value can be used to characterize the coherent errors, the amplitude of the oscillations can be used to estimate the non-Markovian errors. We prove a sandwich theorem to establish a linear relation between the amplitude of the oscillations and the magnitude of the non-Markovian errors. Moreover, by varying the initial state, we show that the steady state values can be used to completely construct the generator of the coherent errors. Using these results, we develop two experimental protocols to characterize the unitary errors based on these results, one of which requires single-qubit addressing and the other one doesn't. We also develop a protocol to characterize non-Markovian errors. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Progress report on analytic and numerical studies of x-ray-induced impulse

In a variety of high-energy-density (HED) systems, x-rays of a given energy are used to generate shockwaves, bulk motion, and impulse in materials. As radiation energy is deposited within the material, the material is heated. The number of photons of a certain wavelength absorbed is determined by the spectral intensity of the radiation and the material’s opacity evaluated at that wavelength. This heating results in a pressure increase dictated by the material’s equation of state. Depending on the intensity of the absorbed radiation, it may also cause the material to change phase into a liquid, gas, or plasma. The increased pressure drives the heated surface layer to blow off, imparting impulse and sending a compressive wave into the bulk of the material. Additionally, the compression wave interacts with the solid boundary of the material, resulting in a tensile wave that may cause the material to spall. The impulse generated by the deposition of x-ray energy within the sample can be modeled using purely analytical methods, e.g. the Bethe, Bade, Averell, and Yost (BBAY) model. However, the blow-off process is rather complicated, and proper modeling efforts must account for material ejected by spallation, vaporization, jetting, and plasma ablation. For this reason, analytical models have an unclosed term describing the final energy of the blown-off material Ef(z). Prior modeling efforts have arbitrarily fixed this at some value or modeled it using a limiting set of thermodynamic assumptions. The work we are currently performing uses validated simulations using sophisticated photon transport, equation of state, and strength models/data to provide a fit for Ef(z) that is useful for predictive calculation of impulse. We will apply our methodology and show results for different materials and x-ray sources. This work is particularly useful for the design of experiments studying x-ray impulse generation. The present report is outlined as follows. Section 1.2 describes a series of HED experiments investigating x-ray-generated impulse in materials. Section 1.3 describes the computational approach we employ in this study, and presents validation results against the aforementioned experiments. Section 1.4 introduces analytical models for impulse generation, as well as our method for utilizing impulse from simulations to close the models. We discuss the concept of impulse-spectrum sensitivity, it’s application to uncertainty quantification, and derive a very useful analytical expression for it in 1.5. We then summarize recent progress in this project and discuss future work in section 1.6.

36 MATERIALS SCIENCE↗
DOE OSTI2024

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↗
DOE OSTI2025

Direct experimental constraints on the spatial extent of a neutrino wavepacket

Despite their high relative abundance in our Universe, neutrinos are the least understood fundamental particles of nature. In fact, the quantum properties of neutrinos emitted in experimentally relevant sources are theoretically contested and the spatial extent of the neutrino wavepacket is only loosely constrained by reactor neutrino oscillation data with a spread of 13 orders of magnitude. Here we present a method to directly access this quantity by precisely measuring the energy width of the recoil daughter nucleus emitted in the radioactive decay of beryllium-7. The final state in the decay process contains a recoiling lithium-7 nucleus, which is entangled with an electron neutrino at creation. The lithium-7 energy spectrum is measured to high precision by directly embedding beryllium-7 radioisotopes into a high-resolution superconducting tunnel junction that is operated as a cryogenic sensor. Under this approach, we set a lower limit on the Heisenberg spatial uncertainty of the recoil daughter of 6.2 pm, which implies that the final-state system is localized at a scale more than a thousand times larger than the nucleus itself. From this measurement, the first, to our knowledge, direct lower limit on the spatial extent of a neutrino wavepacket is extracted. These results may have implications in several areas including the theoretical understanding of neutrino properties, the nature of localization in weak nuclear decays and the interpretation of neutrino physics data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗