SEARCH · Engineering Papers
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.
This presentation provides an overview of computational physics, the role LANL played in the history of computational physics, and how computational physics differs at LANL. Next, the presentation addresses computational physics as part of the LANL weapons program. Finally, the computational physics structure at LANL is addressed.
Performance of the Rigetti quantum computing platform was tested during the period of time from 12/16/2019 to 05/18/2020. The Grover's search (GS) algorithm was used as a testing tool, in particular, 3- and 4-level versions of the algorithm. As a result, a number of hardware issues were revealed, so the algorithm was split to smaller blocks and individual gates, and all of them were tested separately. The fidelity decay was found to be due to both decoherent processes and coherent errors of the native hardware gates. These errors were estimated from RX gate benchmarks and were in a form of extra rotation of the quantum state. As a consequence of that, performance of the algorithm was shown to be strongly dependent on the native gate decomposition of the program. Suggestions and possible improvements of the future Rigetti runs are made based on the obtained observations.
A large scientific and technological effort is underway to investigate the properties of two‐dimensional (2D) materials to become building blocks in integrated nano‐electronic and photonic circuits, composites, coatings, energy harvesting nano‐ systems, nano‐sensors, and nano‐electro‐mechanical systems (NEMS).While several experiments and calculations have revealed exciting novel phenomena in these nanostructures, many scientific and technological questions remain open. In particular, understanding and controlling the structure and mechanical properties of 2D materials at the interface with a solid surface is of key importance to enable the aforementioned nanotechnologies, as well as to fulfill the potential of strain engineering to tailor their properties. The vision of this DoE research program was to investigate novel mechanical properties and phenomena in 1D and 2D materials with the overarching goal of defining a new basic understanding of mechanical behavior in nano and quantum systems. The group of the PI has developed in the last years several studies on the mechanical properties of Carbon nanotubes and oxide nanobelts, more recently the PI has focused her attention on the properties of two‐dimensional materials, such as graphene and MoS 2 , which are a few‐atomic‐layer thick films and hold a great potential for technological and energy applications. The most studied 2D material is graphene, existing as a single layer of graphite or a few‐layer thick epitaxial graphene film. Graphene possesses a large in‐plane Young’s modulus as well as high intrinsic carrier mobility, and high in‐plane thermal conductivity. Besides graphene, also 2D films of graphene oxide (GO), hexagonal Boron Nitride (h‐BN), and transition metal dichalcogenides such as MoS 2 exhibit unique and excellent properties and hold great promise for nanotechnology applications. This research program was aimed at developing a new basic understanding of the mechanical properties of 1D and 2D materials at the interface with solid substrates, with a focus on the role of defects, materials structure, and substrate interaction. Ultimately, we aimed at developing the basic knowledge and experimental tools for manipulating 1D and 2D materials’ structure, stacking, substrate interaction, defects, and number of layers to realize an entirely new class of ultra‐hard, ultra‐thin, and ultra‐light active materials with ad‐ hoc mechanical properties for a variety of applications.
Experimentally measured reduced density matrices (RDMs) often violate constraints that ensure they represent N-electron states—known as N-representability conditions—because of statistical and hardware noise. In this work, we present a correlated purification framework based on semidefinite programming to restore the accuracy of a noisy, unphysical two-electron RDM (2-RDM). The method performs a bi-objective optimization that minimizes both the many-electron energy and the nuclear norm of the correction to the measured 2-RDM. The nuclear norm, often employed in matrix completion, promotes low-rank corrections, while the energy term acts as a regularization term that can improve the purity of the ground state. While the method is particularly effective for ground states, it can also be applied to excited and nonstationary states by decreasing the weight of the energy relative to the error norm. In an application to fermionic shadow tomography of large hydrogen chains, correlated purification yields substantial reductions in both energy and 2-RDM error, achieving chemical accuracy across dissociation curves. This framework provides a robust strategy for tomography in many-body quantum simulations.
Decays of the neutral and long-lived $\eta$ and $\eta'$ mesons provide a unique, flavor-conserving laboratory to test low-energy Quantum Chromodynamics and search for new physics beyond the Standard Model. They have drawn world-wide attention in recent years and have inspired broad experimental programs in different high-intensity-frontier centers. Additionally, new experimental data will offer critical inputs to precisely determine the light quark mass ratios, $\eta$-$\eta'$ mixing parameters, and hadronic contributions to the anomalous magnetic moment of the muon. At the same time, it will provide a sensitive probe to test potential new physics. This includes searches for hidden photons, light Higgs scalars, and axion-like particles that are complementary to worldwide efforts to detect new light particles below the GeV mass scale, as well as tests of discrete symmetry violation. In this review, we give an update on theoretical developments, discuss the experimental opportunities, and identify future research needed in this field.
Millikelvin cryogenics is an emergent area for Fermilab with broad applications in particle detection for the cosmic frontier and quantum information science. Fermilab Cryogenic Engineering currently supports multiple experiments, including the SuperCDMS SNOLAB and ADMX dark matter experiments, along with several millikelvin test facilities. More programs are expected to come online in future.
QSCOUT is the Quantum Scientific Computing Open User Testbed, a trapped-ion quantum computer testbed realized at Sandia National Laboratories on behalf of the Department of Energy's Office of Science and its Advanced Scientific Computing Research (ASCR) program.
Highly emissive organic molecular crystals find applications in several areas, such as organic electronics, solar cells, and sensors. Understanding the excited-state mechanisms underlying these applications is essential for optimizing and controlling them effectively. Exciton models coupled with nonadiabatic dynamics, particularly quantum dynamics, provide crucial insights into photochemical and photophysical processes in molecular crystals. Nevertheless, there remains a lack of general tools and automated workflows to facilitate such simulations. In this paper, we present a computational strategy to investigate the photoactivated dynamics of organic molecular crystals, bridging methodologies traditionally used for molecular systems and materials science, with a particular focus on the interplay between local excitations and charge transfer (CT) processes. We have implemented an interface between the fromage and Overdia programs, enabling the construction of vibronic Hamiltonians for molecular crystals within an excited-state ONIOM(QM:QM′) framework, incorporating long-range electrostatics through a RESP-based Ewald summation. Fragment-based diabatization provides a route to quantum dynamics simulations in weak-to-intermediate coupling regimes. The method was applied to the photophysics of dibenzo[g,p]chrysene (DBC) crystals using time-dependent DFT. The fromage/ Overdia interface was employed to compute the couplings of local excitations and CT states for 18 unique DBC dimers in the crystal and to quantify the influence of electrostatic embedding, which was found to be modest (10−20%). Simulations on π-stacked dimers reproduced the small red shift observed experimentally from solution to crystal, attributed to electronic interactions among fixed monomers rather than crystal electrostatics. Quantum dynamics simulations revealed ultrafast population transfer from bright local excitations to CT states. This approach establishes a robust framework linking molecular and solid-state excited-state dynamics, with potential applications for studying excitations, defects, and impurities in molecular crystals.
Black holes are chaotic quantum systems that are expected to exhibit random matrix statistics in their finite energy spectrum. Lin, Maldacena, Rozenberg and Shan (LMRS) have proposed a related characterization of chaos for the ground states of BPS black holes with finite area horizons. On a separate front, the “fuzzball program” has uncovered large families of horizon-free geometries that account for the entropy of holographic BPS systems, but only in situations with sufficient supersymmetry to exclude finite area horizons. The highly structured, non-random nature of these solutions seems in tension with strong chaos. We verify this intuition by performing analytic and numerical calculations of the LMRS diagnostic in the corresponding boundary quantum system. In particular we examine the 1/2 and 1/4-BPS sectors of \mathcal{N}=4 𝒩 = 4 SYM, and the two charge sector of the D1-D5 CFT. We find evidence that these systems are only weakly chaotic, with a Thouless time determining the onset of chaos that grows as a power of N N . In contrast, finite horizon area BPS black holes should be strongly chaotic, with a Thouless time of order one. In this case, finite energy chaotic states become BPS as N N is decreased through the recently discovered “fortuity” mechanism. Hence they can plausibly retain their strongly chaotic character.
We apply optimal control theory to drive a polar active fluid into new behaviors: relocating asters, reorienting waves, and on-demand switching between states. This study reveals general principles to program active matter for useful functions.
The project involved an integrated program of research which focused on the use of the extra-dimensional geometry of string theory as a tool in understanding a variety of phenomena in quantum field theory. This was used to construct and study new examples of quantum field theories with the eventual aim of also using these results to develop new particle physics and cosmology scenarios. The research involved the development of new methods for constructing such extra-dimensional geometries, and also the study of the resulting quantum field theories generated by this process.
We report the MechElastic Python package evaluates the mechanical and elastic properties of bulk and 2D materials using the elastic coefficient matrix ( C ij ) obtained from any ab-initio density-functional theory (DFT) code. The current version of this package reads the output of VASP, ABINIT, and Quantum Espresso codes (but it can be easily generalized to any other DFT code) and performs the appropriate post-processing of elastic constants as per the requirement of the user. This program can also detect the input structure's crystal symmetry and test the mechanical stability of all crystal classes using the Born-Huang criteria. Various useful material-specific properties such as elastic moduli, longitudinal and transverse elastic wave velocities, Debye temperature, elastic anisotropy, 2D layer modulus, hardness, Pugh's ratio, Cauchy's pressure, Kleinman's parameter, and Lame's coefficients, can be estimated using this program. Another existing feature of this program is to employ the ELATE package (2016) [29] and plot the spatial variation of several elastic properties such as Poisson's ratio, linear compressibility, shear modulus, and Young's modulus in three dimensions. Further, the MechElastic package can plot the equation of state (EOS) curves for energy and pressure for a variety of EOS models such as Murnaghan, Birch, Birch-Murnaghan, and Vinet, by reading the inputted energy/pressure versus volume data obtained via numerical calculations or experiments. This package is particularly useful for the high-throughput analysis of elastic and mechanical properties of materials.
Here, O 2 adsorption on MnO(100) precovered with sodium (Na) multilayers was investigated by X-ray photoelectron spectroscopy (XPS) and temperature programmed desorption (TPD). Deposition of Na multilayers leads to a first monolayer of oxidic Na followed by metallic Na island growth. XPS results for the oxidation of the metallic Na islands indicate an incomplete oxidation of Na at 125 K. Oxidation at 350 K completely oxidizes the metallic islands producing a mixture of Na 2 O and Na 2 O 2 . Thermal evolution of the oxidation products was examined. After oxidation at 350 K and flashing to 750 K, Na 2 O and Na 2 O 2 are the primary oxides on MnO(100). After flashing to 850 K, a solid state reaction of Na 2 O 2 /Na 2 O and the MnO(100) substrate forms a NaMnO 2 -like surface compound which decomposes above 850 K. Oxygen exchange between CO 2 and Na oxides is observed. The strong interaction between CO 2 and Na oxide islands forms Na 2 CO 3 on MnO(100). Heating the Na 2 CO 3 covered MnO(100) surface to 600 K for 10 min drives oxidic Na in the first monolayer into the MnO subsurface, and produces a surface exposing islands of Na 2 CO 3 on MnO(100).
Abstract Graph states are a broad family of entangled quantum states, each defined by a graph composed of edges representing the correlations between subsystems. Such states constitute versatile resources for quantum computation and quantum-enhanced measurement. Their generation and engineering require a high level of control over entanglement. Here we report on the generation of continuous-variable graph states of atomic spin ensembles, which form the nodes of the graph. We program the entanglement structure encoded in the graph edges by combining global photon-mediated interactions in an optical cavity with local spin rotations. By tuning the entanglement between two subsystems, we either localize correlations within each subsystem or enable Einstein–Podolsky–Rosen steering—a strong form of entanglement that enables the extraction of precise information from one subsystem through measurements on the other. We further engineer a four-mode square graph state, highlighting the flexibility of our approach. Our method is scalable to larger and more complex graphs, laying groundwork for measurement-based quantum computation and advanced protocols in quantum metrology.
Computational simulation of non-adiabatic molecular dynamics is an indispensable tool for understanding complex photoinduced processes such as internal conversion, energy transfer, charge separation, and spatial localization of excitons, to name a few. Here, we report an implementation of the fewest-switches surface hopping algorithm in the NWChem computational chemistry program. Here the surface hopping method is combined with linear-response time-dependent density functional theory calculations of adiabatic excited state potential energy surfaces. To treat quantum transitions between arbitrary electronic Born{Oppenheimer states, we have implemented both numerical and analytical differentiation schemes for derivative non-adiabatic couplings. A numerical approach for the time-derivative non-adiabatic couplings together with an analytical method for calculating non-adiabatic coupling vectors is an efficient combination for surface hopping approaches. Additionally, electronic decoherence schemes and a state reassigned unavoided crossings algorithm are also implemented to improve the accuracy of the simulated dynamics and to handle trivial unavoided crossings. We apply our code to study the ultrafast decay of photoexcited benzene, including a detailed analysis of the potential energy surface, population decay time scales, and vibrational coordinates coupled to the excitation dynamics. The development presented in this work is a baseline for future implementations of more sophisticated frameworks for simulating of non-adiabatic molecular dynamics in NWChem.
Approximate quantum error correction (AQEC) not only dictates the performance of discrete- and continuous-variable quantum error correction codes but also serves as a unifying framework across various physical disciplines. Identifying the optimal recovery channel to maximize the entanglement fidelity via standard semidefinite programming is computationally bottlenecked by the exponentially growing number of Kraus operators with system size, rendering large-scale optimization prohibitive. While analytical near-optimal maps exist, they typically work only when the Knill-Laflamme conditions are nearly satisfied. In this Letter, we establish an efficient framework by leveraging the duality between recovery and environment decoupling. This framework yields a tighter analytical lower bound on entanglement fidelity than the conventional limit set by the transpose channel. Furthermore, by exploiting the decayed weights of noise Kraus operators, we introduce a framework based on principal component analysis to reduce the dimension. In thermal loss channels where the weights decay exponentially, our approach achieves a 33-fold computational speedup while maintaining rigorous accuracy. Our approach enables high-precision optimization for AQEC codes that were previously intractable due to the curse of dimensionality.
Interactions govern the flow of information and the formation of correlations between constituents of many-body quantum systems, dictating phases of matter found in nature and forms of entanglement generated in the laboratory. Typical interactions decay with distance and thus produce a network of connectivity governed by geometry—such as the crystalline structure of a material or the trapping sites of atoms in a quantum simulator. However, many envisioned applications in quantum simulation and computation require more complex coupling graphs including non-local interactions, which feature in models of information scrambling in black hole, and mappings of hard optimization problems onto frustrated classical magnets. Here we describe the realization of programmable non-local interactions in an array of atomic ensembles within an optical cavity, in which photons carry information between atomic spins. By programming the distance dependence of the interactions, we access effective geometries for which the dimensionality, topology and metric are entirely distinct from the physical geometry of the array. As examples, we engineer an antiferromagnetic triangular ladder, a Möbius strip with sign-changing interactions and a treelike geometry inspired by concepts of quantum gravity. The tree graph constitutes a toy model of holographic duality, in which the quantum system lies on the boundary of a higher-dimensional geometry that emerges from measured correlations. Overall, our work provides broader prospects for simulating frustrated magnets and topological phases, investigating quantum optimization paradigms and engineering entangled resource states for sensing and computation.
Solving problems related to planning and operations of large-scale power systems is challenging on classical computers due to their inherent nature as mixed-integer and nonlinear problems. Quantum computing provides new avenues to approach these problems. We develop a hybrid quantum-classical algorithm for the Unit Commitment (UC) problem in power systems which aims at minimizing the total cost while optimally allocating generating units to meet the hourly demand of the power loads. The hybrid algorithm combines a variational quantum algorithm (VQA) with a classical Benders-type heuristic. The resulting algorithm computes approximate solutions to UC in three stages: i) a collection of UC vectors capable meeting the power demand with lowest possible operating costs is generated based on VQA; ii) a classical sequential least squares programming (SLSQP) routine is leveraged to find the optimal power level corresponding to a predetermined number of candidate vectors; iii) in the last stage, the approximate solution of UC along with generating units power level combination is given. To demonstrate the effectiveness of the presented method, three different systems with 3 generating units, 10 generating units, and 26 generating units were tested for different time periods. In addition, convergence of the hybrid quantum-classical algorithm for select time periods is proven out on IonQ's Forte system.