Quantum Computing at Los Alamos
Explore the source record for details and available documents.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Understanding how isolated quantum systems thermalize is central to both fundamental physics and the development of quantum technologies. In this Perspective, we present recent and new results on thermalization in nonabelian gauge theory from exact real-time simulations of two-dimensional lattice models. We discuss tests of the eigenstate thermalization hypothesis, entanglement entropy dynamics, the absence of persistent non-thermal states, spectral signatures, and measures of quantum complexity. Our new results include a two-step thermalization process in localized regions and the role of higher gauge field representations. These findings suggest that thermalization in gauge theories may offer a test case for quantum advantage.
A search for a pseudoscalar a produced in association with a top-quark pair, or in association with a single top quark plus a W boson, with the pseudoscalar decaying into b-quarks $(a → b\bar{b})$, is performed using the full Run 2 data sample using a dileptonic decay mode signature. The search covers pseudoscalar boson masses between 12 and 100 GeV and involves both the kinematic regime where the decay products of the pseudoscalar are reconstructed as two standard b-tagged small-radius jets, or merged into a large-radius jet due to its Lorentz boost. No significant excess relative to expectations is observed. Assuming a branching ratio BR$(a → b\bar{b})$ = 100%, the range of pseudoscalar masses between 50 and 80 GeV is excluded at 95% confidence level for a coupling of the pseudoscalar to the top quark of 0.5, while a coupling of 1.0 is excluded at 95% confidence level for the masses considered, with the coupling defined as the strength modifier of the Standard Model Yukawa coupling.
A search for pair-production of vector-like leptons is presented, considering their decays into a third-generation Standard Model (SM) quark and a vector leptoquark ( U 1 ) as predicted by an ultraviolet-complete extension of the SM, referred to as the ‘4321’ model. Given the assumed decay of U 1 into third-generation SM fermions, the final state can contain multiple τ-leptons and b-quarks. This search is based on a dataset of pp collisions at $\sqrt{s}=13$ TeV recorded with the ATLAS detector during Run 2 of the Large Hadron Collider, corresponding to an integrated luminosity of up to $140~\textrm{fb}^{-1}$. No significant excess above the SM background prediction is observed, and 95% confidence level limits on the cross-section times branching ratio are derived as a function of the vector-like lepton mass. A lower observed (expected) limit of 910 GeV (970 GeV) is set on the vector-like lepton mass. Additionally, the results are interpreted for a supersymmetric model with an R-parity violating coupling to the third-generation quarks and leptons. Lower observed (expected) limits are obtained on the higgsino mass at 880 GeV (940 GeV) and on the wino mass at 1170 GeV (1170 GeV).
Magnetic skyrmions are chiral spin textures with rich physics and great potential for unconventional computing. Typically, skyrmions form in bulk crystals with reduced symmetry or ultrathin film multilayers involving heavy metals. Here, the formation of fractional Bloch skyrmion tubes at room temperature is demonstrated by 3D printing ferromagnetic double‐helix nanowires with two regions of opposite chirality. Using X‐ray microscopy and micromagnetic simulations, it is shown that the coexistence of vortex and anti‐parallel spin states induces the formation of fractional skyrmion tubes at zero magnetic fields, minimizing the energy cost of breaking the coupling between geometric and magnetic chirality. Control over zero‐field states is also demonstrated, including pure vortex, or mixed skyrmion‐vortex states, highlighting the magnetic reconfigurability of these 3D nanowires. This work shows how interfacing chiral geometries at the nanoscale can enable advanced forms of topological spintronics.
A search for a light charged Higgs boson produced in decays of the top quark, $t \rightarrow H^{\pm } b$ with $H^{\pm } \rightarrow cs$, is presented. This search targets the production of top-quark pairs $t\bar{t} \rightarrow WbH^{\pm } b$, with $W \rightarrow ℓv(ℓ = e, μ)$, resulting in a lepton-plus-jets final state characterised by an isolated electron or muon and at least four jets. The search exploits b-quark and c-quark identification techniques as well as multivariate methods to suppress the dominant $t\bar{t}$ background. The data analysed correspond to 140 fb -1 of $pp$ collisions at $\sqrt{s}$ = 13 TeV recorded with the ATLAS detector at the LHC between 2015 and 2018. Observed (expected) 95% confidence-level upper limits on the branching fraction $\mathscr{B}(t \rightarrow H^{\pm } b)$, assuming $\mathscr{B}(t \rightarrow Wb) + \mathscr{B}(t \rightarrow H^{\pm }(\rightarrow cs)b$, are set between 0.066% (0.077%) and 3.6% (2.3%) for a charged Higgs boson with a mass between 60 and 168 GeV.
Abstract Ranking sectors and countries within global value chains is of paramount importance to estimate risks and forecast growth in large economies. However, this task is often non-trivial due to the lack of complete and accurate information on the flows of money and goods between sectors and countries, which are encoded in input–output (I–O) tables. In this work, we show that an accurate estimation of the role played by sectors and countries in supply chain networks can be achieved without full knowledge of the I–O tables, but only relying on local and aggregate information, e.g., the total intermediate demand per sector. Our method, based on a rank-1 approximation to the I–O table, shows consistently good performance in reconstructing rankings (i.e., upstreamness and downstreamness measures for countries and sectors) when tested on empirical data from the world input–output database. Moreover, we connect the accuracy of our approximate framework with the spectral properties of the I–O tables, which ordinarily exhibit relatively large spectral gaps. Our approach provides a fast and analytically tractable framework to rank constituents of a complex economy without the need of matrix inversions and the knowledge of finer intersectorial details.
Double shell targets are an alternative ignition platform for inertial confinement fusion. One design consideration for double shell targets is the choice of inner shell material to help trap radiation emitted by the hot fuel to aid ignition. Materials such as molybdenum and tungsten are of interest for the inner shell layer of the targets. While molybdenum has a lower density that could inhibit instability growth and allow for radiography and code benchmarking, tungsten has a higher density that could provide more compression and confinement. These tradeoffs have been explored using optimized designs for each material. Our previous work [Vazirani et al., “Coupling 1D xRAGE simulations with machine learning for graded inner shell design optimization in double shell capsules,” Phys. Plasmas 28, 122709 (2021); Vazirani et al., “Coupling multi-fidelity xRAGE with machine learning for graded inner shell design optimization in double shell capsules,” Phys. Plasmas 30, 062704 (2023); and Vazirani et al., “Bayesian batch optimization for molybdenum versus tungsten inertial confinement fusion double shell target design,” Stat. Anal. Data Min. 17, e11698 (2024)] resulted in a multi-fidelity Bayesian optimization framework to find yield-optimized double shell target geometries. By leveraging simulations of varying fidelities (one-dimensional and two-dimensional) to inform one another, the multi-fidelity optimization was able to optimize a design in the highest fidelity with significantly fewer simulations than would be used in a systematic parameter scan. In this work, we apply the multi-fidelity Bayesian optimization to explore the optimized designs of double shell targets with molybdenum and tungsten inner shells as well as the physics producing the high performing implosions. A physics exploration of all the simulations used in this study shows trends in designs that contribute to high yields, ion temperatures, and fuel areal densities. Comparison of molybdenum and tungsten simulations shows that they can produce similar implosion conditions with different geometries, which would be important to study in experiments. Graded density layers produce varying performances with the two materials but continue to be of interest for future studies along with studies of doped inner shell materials and applied surface roughness.
This report represents the first status update from the Creation of Next-generation Nuclear Energy Computational Technology (CONNECT) effort within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program and is intended to satisfy the reporting requirements of the following milestones: • M3MS-24OR0101321: Assess needs, requirements and opportunities for high-fidelity neutronics and transport as used in DOE-NE program. • M3MS-24AN0101301: Generate a community report on the needs, requirements, and opportunities for high-fidelity neutronics and transport as used in DOE-NE programs and industry. • M3MS-24OR0202433: Implement and assess numerical strategies leveraging Monte Carlo neutron transport on GPUs for production analysis. In particular, Sections 1, 2, and 3 detail the assessment of the current state of practice and the landscape of the broader advanced computing world to meet the objectives of M3MS-24OR0101321 and M3MS-24AN0101301, while Sections 4 and 5 describe new work during FY24 to extend current GPU Monte Carlo capabilities in areas relevant to active NEAMS development to satisfy M3MS-24OR0202433. Section 6 offers some thoughts about the potential future impact of Monte Carlo neutronics on NEAMS and the nuclear industry.
We investigate identifying differential equations in the frequency domain. Fourier analysis is an important tool in theoretical analysis and numerical solvers of differential equations, yet there is limited work in exploring this connection in the identification of differential equations. This paper aims to identify the underlying differential equation in the frequency domain, from a given single realization of the differential equation perturbed by noise. Such setting imposes difficulties which are different from other identification methods where computation is carried out in the physical domain. We propose several ways to mitigate the challenges arising from noise in data and large differences in the magnitudes of frequency responses. The main takeaways are that identifying differential equations solely in the frequency domain is challenging, the method we propose is based on a form of domain partitions in the frequency domain, and this method shows benefits for complex data even with high level of noise. We introduce a Fourier feature denoising, and define the meaningful data region and the core regions of features to reduce the effect of noise in the frequency domain and to enhance the accuracy in coefficient identification. The proposed method is tested on various differential equations with linear, nonlinear, and high-order derivative feature terms, and shows advantages on complex data with many frequency modes, even under high level of noise.
Microgrids rely on communication networks for reliable operation, which makes them inherently vulnerable to cyberattacks. Such attacks can destabilise system dynamics and drive states away from their nominal operating trajectories. Although several physics-informed and machine learning-based strategies have been developed to counter these threats, the rapidly evolving cyber landscape enables adversaries to bypass static defences or rules-based mitigation approaches. This paper proposes a dynamic, online-trained and fully decentralised reinforcement learning (RL)-based cyberattack response framework to protect microgrids from evolving cyberattacks. The proposed framework deploys multiple deep Q-networks (DQNs), each associated with a distributed energy resource (DER), to enable localised and adaptive attack mitigation. In this framework, each DQN processes local voltage and frequency measurements—combined with intrusion detection system (IDS) alerts—as observations and rewards to guide decision-making. Extensive simulation studies demonstrate the robustness of the proposed framework under diverse attack scenarios and varying IDS-induced detection delays. Comparative analysis highlights its superiority over existing static or preexisting rules-based mitigation approaches. Finally, we present an analysis that shows the framework's scalability to real-life microgrids with more interacting agents.
Counterdiabatic (CD) driving has the potential to speed up adiabatic quantum state preparation by suppressing unwanted excitations. However, existing approaches either require intractable classical computations or are based on approximations that do not have performance guarantees. We propose and analyze a nonvariational, system-agnostic CD expansion method and analytically show that it converges exponentially quickly in the expansion order. In finite systems, the required resources scale inversely with the spectral gap, which we argue is asymptotically optimal. To extend our method to the thermodynamic limit and suppress errors stemming from high-frequency transitions, we leverage finite-time adiabatic protocols. In particular, we show that a time determined by the quantum speed limit is sufficient to prepare the desired ground state, without the need to optimize the adiabatic trajectory. Numerical tests of our method on the quantum Ising chain show that our method can outperform state-of-the-art variational CD approaches.
This paper presents a search for massive, charged, long-lived particles with the ATLAS detector at the Large Hadron Collider using an integrated luminosity of 140 fb −1 of proton-proton collisions at $\sqrt{s}$ = 13 TeV. These particles are expected to move significantly slower than the speed of light. In this paper, two signal regions provide complementary sensitivity. In one region, events are selected with at least one charged-particle track with high transverse momentum, large specific ionisation measured in the pixel detector, and time of flight to the hadronic calorimeter inconsistent with the speed of light. In the other region, events are selected with at least two tracks of opposite charge which both have a high transverse momentum and an anomalously large specific ionisation. The search is sensitive to particles with lifetimes greater than about 3 ns with masses ranging from 200 GeV to 3 TeV. The results are interpreted to set constraints on the supersymmetric pair production of long-lived R-hadrons, charginos and staus, with mass limits extending beyond those from previous searches in broad ranges of lifetime.
Subspace methods are powerful, noise-resilient methods that can effectively prepare ground states on quantum computers. The challenge is to get a subspace with a small condition number that spans the states of interest using minimal quantum resources. In this work, we will use eigenvector continuation to build a subspace from the low-lying states of a set of Hamiltonians. The basis vectors are prepared using truncated versions of standard state preparation methods such as imaginary-time evolution (ITE), adiabatic state preparation (ASP), and variational quantum eigensolver. By using these truncated methods combined with eigenvector continuation, we can directly improve upon them, obtaining more accurate ground-state energies at a reduced cost. We use several spin systems to demonstrate convergence even when methods like ITE and ASP fail, such as ASP in the presence of level crossings and ITE with vanishing energy gaps. We also showcase the noise resilience of this approach beyond the gains already made by having shallower quantum circuits. Furthermore, our findings suggest that eigenvector continuation can be used to improve existing state preparation methods in the near term.
This paper presents a search for supersymmetric particles in models with highly compressed mass spectra, in events consistent with being produced through vector boson fusion. The search uses 140 fb −1 of proton-proton collision data at $\sqrt{s}$ = 13 TeV collected by the ATLAS experiment at the Large Hadron Collider. Events containing at least two jets with a large gap in pseudorapidity, large missing transverse momentum, and no reconstructed leptons are selected. A boosted decision tree is used to separate events consistent with the production of supersymmetric particles from those due to Standard Model backgrounds. The data are found to be consistent with Standard Model predictions. The results are interpreted using simplified models of R-parity-conserving supersymmetry in which the lightest supersymmetric partner is a bino-like neutralino with a mass similar to that of the lightest chargino and second-to-lightest neutralino, both of which are wino-like. Lower limits at 95% confidence level on the masses of next-to-lightest supersymmetric partners in this simplified model are established between 117 and 120 GeV when the lightest supersymmetric partners are within 1 GeV in mass.
In this CRADA project, Lawrence Livermore National Laboratory, Stanford University, SLAC National Laboratory, and TotalEnergies collaboratively executed a multidisciplinary investigation of electrochemical reduction of CO2 to produce sustainable fuels and chemicals. Overall, the project led to an increased understanding of the fundamental processes involved in CO2 electrolysis, from the atomistic scale to the full electrolyzer device scale, ultimately leading to design guidelines for CO2 electrolyzers that will help in their future commercialization. As the model systems, Ag- and Cu-based catalysts were investigated in various forms depending on the electrochemical platform that was utilized to study the activity, selectivity, and durability towards electrochemical CO2 reduction. By employing experimental, theoretical, and computational techniques, the project team experimentally validated multi-physics models, evaluated the experimental levers that lead to increased electrolyzer reaction selectivity and energy efficiency, and used computational optimization to design higher performance electrodes. The learnings of this project were extensively documented in publicly available peer-reviewed journal publications and conference presentations, which serve as a foundation for further work to build from.
Many extensions of the Standard Model, including those with dark matter particles, propose new mediator particles that decay into hadrons. This paper presents a search for such low mass narrow resonances decaying into hadrons using 140 fb −1 of proton-proton collision data recorded with the ATLAS detector at a centre-of-mass energy of 13 TeV. The resonances are searched for in the invariant mass spectrum of large-radius jets with two-pronged substructure that are recoiling against an energetic photon from initial state radiation, which is used as a trigger to circumvent limitations on the maximum data recording rate. This technique enables the search for boosted hadronically decaying resonances in the mass range 20–100 GeV hitherto unprobed by the ATLAS Collaboration. The observed data are found to agree with Standard Model predictions and 95% confidence level upper limits are set on the coupling of a hypothetical new spin-1 Z′ resonance with Standard Model quarks as a function of the assumed Z′-boson mass in the range between 20 and 200 GeV.
A search for T and Y vector-like quarks produced in proton-proton collisions at a centre-of-mass energy of 13 TeV and decaying into Wb in the fully hadronic final state is presented. The search uses 139 fb −1 of data collected by the ATLAS detector at the LHC from 2015 to 2018. The final state is characterised by a hadronically decaying W boson with large Lorentz boost and a b-tagged jet, which are used to reconstruct the invariant mass of the vector-like quark candidate. The main background is QCD multijet production, which is estimated using a data-driven method. Upon finding no significant excess in data, mass limits at 95% confidence level are obtained as a function of the global coupling parameter, κ. The observed lower limits on the masses of Y quarks with κ = 0.5 and κ = 0.7 are 2.0 TeV and 2.4 TeV, respectively. For T quarks, the observed mass limits are 1.4 TeV for κ = 0.5 and 1.9 TeV for κ = 0.7.