SEARCH · Engineering Papers

Results for “Quantum control theory”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Convergence of Eigenvector Continuation

Eigenvector continuation is a computational method that finds the extremal eigenvalues and eigenvectors of a Hamiltonian matrix with one or more control parameters. It does this by projection onto a subspace of eigenvectors corresponding to selected training values of the control parameters. The method has proven to be very efficient and accurate for interpolating and extrapolating eigenvectors. However, almost nothing is known about how the method converges, and its rapid convergence properties have remained mysterious. In this Letter, we present the first study of the convergence of eigenvector continuation. In order to perform the mathematical analysis, we introduce a new variant of eigenvector continuation that we call vector continuation. We first prove that eigenvector continuation and vector continuation have identical convergence properties and then analyze the convergence of vector continuation. Our analysis shows that, in general, eigenvector continuation converges more rapidly than perturbation theory. The faster convergence is achieved by eliminating a phenomenon that we call differential folding, the interference between nonorthogonal vectors appearing at different orders in perturbation theory. From our analysis we can predict how eigenvector continuation converges both inside and outside the radius of convergence of perturbation theory. Further, while eigenvector continuation is a nonperturbative method, we show that its rate of convergence can be deduced from power series expansions of the eigenvectors. Our results also yield new insights into the nature of divergences in perturbation theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Alkali Metal Control of Triplet-Mediated C–H Activation in Iron-Mediated Coupling of Dinitrogen and Benzene

Direct coupling of N 2 with abundant feedstocks like benzene to form N-containing organic compounds is a promising strategy for N 2 fixation pathways. The challenge of coupling N 2 activation and C–H bond oxidative addition was recently solved by introducing a reversible benzene C–H bond activation process mediated by a low-valent Fe(0) complex, which gave an organometallic product that could couple with partially reduced N 2 . Interestingly, the energetics of the C–H oxidative addition/reductive elimination step depends on the choice of alkali metal. However, the reason why the alkali metal influences the C–H bond activation remained elusive. In this work, we present a comprehensive study on this Fe(0)-mediated reversible C–H activation. Through density functional theory combined with high-level coupled cluster calculations, we discovered that the intermediate-spin triplet (S = 1) controls the energy of the transition state for C–H cleavage, while the high-spin quintet (S = 2) controls the position of the equilibrium. Na + drives the equilibrium toward oxidative addition due to an electrostatic effect, while K + and Rb + are dominated by a steric effect that favors the iron(0) species. The key role played by nonbonding interactions in the Fe-mediated C–H activation provides a conceptual model for alkali control over organometallic transformations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Tree tensor network hierarchical equations of motion based on time-dependent variational principle for efficient open quantum dynamics in structured thermal environments

In this work, we introduce an efficient method, TTN-HEOM, for exactly calculating the open quantum dynamics for driven quantum systems interacting with highly structured bosonic baths by combining the tree tensor network (TTN) decomposition scheme with the bexcitonic generalization of the numerically exact hierarchical equations of motion (HEOM). The method yields a series of quantum master equations for all core tensors in the TTN that efficiently and accurately capture the open quantum dynamics for non-Markovian environments to all orders in the system–bath interaction. These master equations are constructed based on the time-dependent Dirac–Frenkel variational principle, which isolates the optimal dynamics for the core tensors given the TTN ansatz. The dynamics converges to the HEOM when increasing the rank of the core tensors, a limit in which the TTN ansatz becomes exact. We introduce TENSO, tensor equations for non-Markovian structured open systems, as a general-purpose Python code to propagate the TTN-HEOM dynamics. We implement three general propagators for the coupled master equations: two fixed-rank methods that require a constant memory footprint during the dynamics and one adaptive-rank method with a variable memory footprint controlled by the target level of computational error. We exemplify the utility of these methods by simulating a two-level system coupled to a structured bath containing one Drude–Lorentz component and eight Brownian oscillators, which is beyond what can presently be computed using the standard HEOM. Our results show that the TTN-HEOM is capable of simulating both dephasing and relaxation dynamics of driven quantum systems interacting with structured baths, even those of chemical complexity, with an affordable computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Transferable dispersion-aware machine learning interatomic potentials for multilayer transition metal dichalcogenide heterostructures

Stacking atomically thin transition metal dichalcogenides (TMDs) into heterostructures enables exploration of exotic quantum phases, particularly through twist-angle-controlled moiré superlattices. These structures exhibit novel electronic and optical behaviors driven by atomic-scale structural reconstruction. However, studying such systems with DFT is computationally demanding due to their large unit cells and van der Waals (vdW) interactions between layers. To address this, we develop a transferable neural network potential (NNP) that includes long-range vdW corrections up to 12Å with minimal overhead. Trained on vdW-corrected DFT data for Mo- and W-based TMDs with S, Se, and Te, the NNP accurately models monolayers, bilayers, heterostructures, and their interaction with h-BN substrates. It reproduces equilibrium structures, energy landscapes, phonon dispersions, and matches experimental atomic reconstructions in twisted WS2 and MoS2/WSe2 systems. We demonstrate that our NNP achieves DFT-level accuracy and high computational efficiency, enabling large-scale simulations of TMD-based moiré superlattices both with and without substrates.

materials for energy and catalysis↗

DOE OSTI2021

De Sitter decays to infinity

Bubbles of nothing are a class of vacuum decay processes present in some theories with compactified extra dimensions. We investigate the existence and properties of bubbles of nothing in models where the scalar pseudomoduli controlling the size of the extra dimensions are stabilized at positive vacuum energy, which is a necessary feature of any realistic model. We map the construction of bubbles of nothing to a four-dimensional Coleman-De Luccia problem and establish necessary conditions on the asymptotic behavior of the scalar potential for the existence of suitable solutions. We perform detailed analyses in the context of five-dimensional theories with metastable dS 4 x S 1 vacua, using analytic approximations and numerical methods to calculate the decay rate. We find that bubbles of nothing sometimes exist in potentials with no ordinary Coleman-De Luccia decay process, and that in the examples we study, when both processes exist, the bubble of nothing decay rate is typically faster. Our methods can be generalized to other stabilizing potentials and internal manifolds.

classical theories of gravity↗

DOE OSTI2022

Quantum advantages for Pauli channel estimation

We show that entangled measurements provide an exponential advantage in sample complexity for Pauli channel estimation, which is both a fundamental problem and a practically important subroutine for benchmarking near-term quantum devices. The specific task we consider is to simultaneously learn all the eigenvalues of an n-qubit Pauli channel to ±ε precision. We give an estimation protocol with an n-qubit ancilla that succeeds with high probability using only O(n/ε 2 ) copies of the Pauli channel, while proving that any ancilla-free protocol (possibly with adaptive control and channel concatenation) would need at least Ω(2 n/3 ) rounds of measurement. We further study the advantages provided by a small number of ancillas. For the case that a k-qubit ancilla (k≤n) is available, we obtain a sample complexity lower bound of Ω(2 (n-k)/3 ) for any nonconcatenating protocol, and a stronger lower bound of Ω(n 2n-k ) for any nonadaptive, nonconcatenating protocol, which is shown to be tight. We also show how to apply the ancilla-assisted estimation protocol to a practical quantum benchmarking task in a noise-resilient and sample-efficient manner, given reasonable noise assumptions. Our results provide a practically interesting example for quantum advantages in learning and also bring insights for quantum benchmarking.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Casimir-Polder potential on an excited atom near an atomic array

We develop a microscopic description of the fluctuation-mediated Casimir-Polder (CP) shifts on a 'test' two-level atom placed near a two-dimensional atomic array of two-level atoms. We derive the resonant and off-resonant CP potentials experienced by the excited test atom using fourth-order perturbation theory, under the assumption that the test atom resonance is far detuned from those of the array atoms. The total potential on the test atom can be described as the sum of the pairwise resonant and off-resonant potentials resulting from its interaction with the individual atoms of the array. We analyze the asymptotic scaling of CP shifts as a function of the test atom-array separation, and its dependence on various system parameters: array spacing and size, and dipole orientation of the array atoms. Our results bridge the description of CP potential across two distinct regimes: (i) from a single-atom limit where we recover the well-known two-atom Van der Waals potential, (ii) to a macroscopic boundary limit, where we demonstrate new asymptotic scaling laws. We demonstrate that these scaling laws can be tuned via the microscopic parameters of the atomic array, establishing atomically-controlled arrays as a versatile platform for tailoring fluctuation-induced QED phenomena.

FOS: Physical sciences↗

DOE OSTI2020

Photoinduced Dirac semimetal in ZrTe5

Abstract Novel phases of matter with unique properties that emerge from quantum and topological protection present an important thrust of modern research. Of particular interest is to engineer these phases on demand using ultrafast external stimuli, such as photoexcitation, which offers prospects of their integration into future devices compatible with optical communication and information technology. Here, we use MeV Ultrafast Electron Diffraction (UED) to show how a transient three-dimensional (3D) Dirac semimetal state can be induced by a femtosecond laser pulse in a topological insulator ZrTe 5 . We observe marked changes in Bragg diffraction, which are characteristic of bond distortions in the photoinduced state. Using the atomic positions refined from the UED, we perform density functional theory (DFT) analysis of the electronic band structure. Our results reveal that the equilibrium state of ZrTe 5 is a topological insulator with a small band gap of ~ 25 meV, consistent with angle-resolved photoemission (ARPES) experiments. However, the gap is closed in the presence of strong spin-orbit coupling (SOC) in the photoinduced transient state, where massless Dirac fermions emerge in the chiral band structure. The time scale of the relaxation dynamics to the transient Dirac semimetal state is remarkably long, τ ~ 160 ps, which is two orders of magnitude longer than the conventional phonon-driven structural relaxation. The long relaxation is consistent with the vanishing density of states in Dirac spectrum and slow spin-repolarization of the SOC-controlled band structure accompanying the emergence of Dirac fermions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

SIERRA Multimechanics Module: Aria Thermal Theory Manual (V.4.56)

Aria is a Galerkin finite element based program for solving coupled-physics problems described by systems of PDEs and is capable of solving nonlinear, implicit, transient and direct-to-steady state problems in two and three dimensions on parallel architectures. The suite of physics currently supported by Aria includes thermal energy transport, species transport, and electrostatics as well as generalized scalar, vector and tensor transport equations. Additionally, Aria includes support for manufacturing process flows via the incompressible Navier-Stokes equations specialized to a low Reynolds number ( Re < 1) regime. Enhanced modeling support of manufacturing processing is made possible through use of either arbitrary Lagrangian-Eulerian (ALE) and level set based free and moving boundary tracking in conjunction with quasi-static nonlinear elastic solid mechanics for mesh control. Coupled physics problems are solved in several ways including fully-coupled Newton’s method with analytic or numerical sensitivities, fully-coupled Newton-Krylov methods and a loosely-coupled nonlinear iteration about subsets of the system that are solved using combinations of the aforementioned methods. Error estimation, uniform and dynamic h -adaptivity and dynamic load balancing are some of Aria’s more advanced capabilities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Comprehensive view of microscopic interactions between DNA-coated colloids

Abstract The self-assembly of DNA-coated colloids into highly-ordered structures offers great promise for advanced optical materials. However, control of disorder, defects, melting, and crystal growth is hindered by the lack of a microscopic understanding of DNA-mediated colloidal interactions. Here we use total internal reflection microscopy to measure in situ the interaction potential between DNA-coated colloids with nanometer resolution and the macroscopic melting behavior. The range and strength of the interaction are measured and linked to key material design parameters, including DNA sequence, polymer length, grafting density, and complementary fraction. We present a first-principles model that screens and combines existing theories into one coherent framework and quantitatively reproduces our experimental data without fitting parameters over a wide range of DNA ligand designs. Our theory identifies a subtle competition between DNA binding and steric repulsion and accurately predicts adhesion and melting at a molecular level. Combining experimental and theoretical results, our work provides a quantitative and predictive approach for guiding material design with DNA-nanotechnology and can be further extended to a diversity of colloidal and biological systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

The cosmological bootstrap: inflationary correlators from symmetries and singularities

Scattering amplitudes at weak coupling are highly constrained by Lorentz invariance, locality and unitarity, and depend on model details only through coupling constants and the particle content of the theory. For example, four-particle amplitudes are analytic for contact interactions and have simple poles with appropriately positive residues for tree-level exchange. In this paper, we develop an understanding of inflationary correlators which parallels that of flat-space scattering amplitudes. Specifically, we study slow-roll inflation with weak couplings to extra massive particles, for which all correlation functions are controlled by an approximate conformal symmetry on the boundary of the spacetime. After systematically classifying all possible contact terms in de Sitter space, we derive an analytic expression for the four-point function of conformally coupled scalars mediated by the tree-level exchange of massive scalars. Conformal symmetry implies that the correlator satisfies a pair of differential equations with respect to spatial momenta, encoding bulk time evolution in purely boundary terms. The absence of unphysical singularities (and the correct normalization of physical ones) completely fixes this correlator. Moreover, a “spin-raising” operator relates it to the correlators associated with the exchange of particles with spin, while “weight-shifting” operators map it to the four-point function of massless scalars. We explain how these de Sitter four-point functions can be perturbed to obtain inflationary three-point functions. Using our formalism, we reproduce many classic results in the literature, such as the three-point function of slow-roll inflation, and provide a complete classification of all inflationary three- and four-point functions arising from weakly broken conformal symmetry. Remarkably, the inflationary bispectrum associated with the exchange of particles with arbitrary spin is completely characterized by the soft limit of the simplest scalar-exchange four-point function of conformally coupled scalars and a series of contact terms. Finally, we demonstrate that the inflationary correlators contain flat-space scattering amplitudes via a suitable analytic continuation of the external momenta, which can also be directly connected with the signals for particle production seen in the squeezed limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

SIERRA Multimechanics Module: Aria Thermal Theory Manual - Version 5.8

Aria is a Galerkin finite element based program for solving coupled-physics problems described by systems of PDEs and is capable of solving nonlinear, implicit, transient and direct-to-steady state problems in two and three dimensions on parallel architectures. The suite of physics currently supported by Aria includes thermal energy transport, species transport, and electrostatics as well as generalized scalar, vector and tensor transport equations. Additionally, Aria includes support for manufacturing process ows via the incompressible Navier-Stokes equations specialized to a low Reynolds number (Re < 1) regime. Enhanced modeling support of manufacturing processing is made possible through use of either arbitrary Lagrangian-Eulerian (ALE) and level set based free and moving boundary tracking in conjunction with quasi-static nonlinear elastic solid mechanics for mesh control. Coupled physics problems are solved in several ways including fully-coupled Newton’s method with analytic or numerical sensitivities, fully-coupled Newton-Krylov methods and a loosely-coupled nonlinear iteration about subsets of the system that are solved using combinations of the aforementioned methods. Error estimation, uniform and dynamic ℎ-adaptivity and dynamic load balancing are some of Aria’s more advanced capabilities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Constitutive Model of Radiation Aging Effects in Filled Silicone Elastomers under Strain

Filled silicone elastomers, an essential component in many technological applications, are often subjected to controlled or unintended radiation for a variety of reasons. Radiation exposure can lead to permanent mechanical and structural changes in the material, which is manifested as altered mechanical response, and in some cases, a permanent set. For unfilled elastomers, network theories developed and refined over decades can explain these effects in terms of chain-scission and cross-link formation and a hypothesis involving independent networks formed at different strain levels of the material. In this work, we expose a filled silicone rubber to gamma radiation while being under finite elongational strain and show that the observed mechanical and structural changes can be quantitatively modeled within the same theoretical framework developed for unfilled elastomers as long as nuances associated with the Mullins effect are accounted for in a consistent manner. In this work, we employ Ogden’s incompressible hyperelastic model within the framework of Tobolsky’s two-network scheme to describe the observed permanent set and mechanical modulus changes as a function of radiation dosage. In the process, we conclude that gamma radiation induces both direct cross-linking at chain crossings (H-links) and main-chain-scission followed by cross-linking (Y-links). We provide an estimate of the ratio of chain-scission to cross-linking rates, which is in reasonable agreement with previous experimental estimate from Charlesby–Pinner analysis. We use density functional theory (DFT)-based quantum mechanical calculations to explore the stability of –Si and –SiO radicals that form upon a radiation-induced chain-scission event, which sheds light on the relative rates of Y-linking and H-linking processes.

36 MATERIALS SCIENCE↗

DOE OSTI2020

The concept of spin ice graphs and a field theory for their charges

Originally detected in rare earth pyrochlores, spin ice physics is now being artificially extended to a variety of geometries that control collective behavior and exotic properties, making graph theory their proper framework. We relate spin ice notions, such as ice rule, ice manifold, Coulomb phases, charges, and monopoles, to graph-theoretical notions, such as balance, in/out-degrees, and Euler paths. We then propose a field-theoretical treatment in which topological charges and monopoles are the degrees of freedom, while the binary spins are subsumed in an entropic interaction among charges. We show that for a spin ice on a graph in a Gaussian approximation, the kernel of the entropic interaction is the inverse of the graph Laplacian, and we compute screening functions from the graph spectra as Green operators for the screened Poisson problem on a graph. We then apply the treatment to star graphs, tournaments, cycles, and regular spin ice in different dimensions. Our aim is twofold: to set spin ice physics in a proper graph setting, where only topological rather than geometrical notions hold, and to invite graph theorists to contribute their powerful tools to the field of spin ice.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Performance of Heterostructural TaC/AlGaN Schottky Diodes Based on First Principles Electronic Structure Properties

Advances in ultra-wide bandgap materials, such as high Al-content AlxGa1-xN (AlGaN), are essential for next generation power electronics, but the requirement for lattice matched substrates is currently a significant obstacle. Recently, conductive TaC has emerged as a promising virtual substrate for AlGaN heteroepitaxy, with wurtzite (0001) AlxGa1-xN lattice-matched to rocksalt (111) TaC at x ~ 0.5. Thus, understanding and controlling the electronic properties of the TaC/AlGaN interface is key for developing technological applications based on TaC/AlGaN devices. Using density functional theory and electronic structure calculations, we here investigate TaC/Al0.5Ga0.5N interfaces, where we include explicit alloy models in the slab calculations. We predict the Schottky barrier height and the electric field discontinuity resulting from interface charges. Considering all possible combinations of (Ta, C) substrate termination, (Al/Ga, N) nucleation, and (Al/Ga, N) polarity, we construct a chemical potential phase diagram to identify the stable interfaces that can be accessed through variation of the synthesis conditions. The predicted interface electronic properties are implemented in device performance simulations to demonstrate a practical design for a strain-free, high-efficiency TaC/AlGaN Schottky diode with a low barrier height and without interface charges, underscoring the potential of TaC as a substrate for ultra-wide bandgap devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Scattering neutrinos, spin models, and permutations

We consider a class of Heisenberg all-to-all coupled spin models inspired by neutrino interactions in a supernova with N degrees of freedom. These models are characterized by a coupling matrix that is relatively simple in the sense that there are only a few, relative to N , nontrivial eigenvalues, in distinction to the classic Heisenberg spin-glass models, leading to distinct behavior in both the high-temperature and low-temperature regimes. When the momenta of the neutrinos are uniform and random in directions, we can calculate the large- N partition function for the quantum Heisenberg model. In particular, the high-temperature partition function predicts a non-Gaussian density of states, providing interesting counterexamples showing the limits of general theorems on the density of states for quantum spin models. We can repeat the same argument for classical Heisenberg models, also known as rotor models, and we find the high-temperature expansion is completely controlled by the eigenvalues of the coupling matrix, and again predicts non-Gaussian behavior for the density of states as long as the number of eigenvalues does not scale linearly with N . Indeed, we derive the amusing fact that these partition functions are essentially the generating function for counting permutations in the high-temperature regime. Finally, for the case relevant to neutrinos in a supernova, we identify the low-temperature phase as a unique state with the direction of the momenta of the neutrino dictating its coherent state in flavor-space, a state we dub the “flavor-momentum-locked” state. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

In situ electric-field control of ferromagnetic resonance in the low-loss organic-based ferrimagnet V[TCNE] x ∼2

We demonstrate indirect electric-field control of ferromagnetic resonance (FMR) in devices that integrate the low-loss, molecule-based, room-temperature ferrimagnet vanadium tetracyanoethylene (V[TCNE]x∼2) mechanically coupled to PMN-PT piezoelectric transducers. Upon straining the V[TCNE]x films, the FMR frequency is tuned by more than 6 times the resonant linewidth with no change in Gilbert damping for samples with α = 6.5 × 10−5. We show this tuning effect is due to a strain-dependent magnetic anisotropy in the films and find the magnetoelastic coefficient |λs| ∼ (1–4.4) ppm, backed by theoretical predictions from density-functional theory calculations and magnetoelastic theory. Noting the rapidly expanding application space for strain-tuned FMR, we define a new metric for magnetostrictive materials, magnetostrictive agility, given by the ratio of the magnetoelastic coefficient to the FMR linewidth. This agility allows for a direct comparison between magnetostrictive materials in terms of their comparative efficacy for magnetoelectric applications requiring ultra-low loss magnetic resonance modulated by strain. With this metric, we show V[TCNE]x is competitive with other magnetostrictive materials, including YIG and Terfenol-D. This combination of ultra-narrow linewidth and magnetostriction, in a system that can be directly integrated into functional devices without requiring heterogeneous integration in a thin film geometry, promises unprecedented functionality for electric-field tuned microwave devices ranging from low-power, compact filters and circulators to emerging applications in quantum information science and technology.

Materials Science↗

DOE OSTI2021

Cold ion beam in a storage ring as a platform for large-scale quantum computers and simulators: Challenges and directions for research and development

The purpose of this paper is to evaluate the possibility of constructing a large-scale storage-ring-type ion-trap system capable of storing, cooling, and controlling a large number of ions as a platform for scalable quantum computing (QC) and quantum simulations. In such a trap, the ions form a crystalline beam moving along a circular path with a constant velocity determined by the frequency and intensity of the cooling lasers. In this paper, we consider a large leap forward in terms of the number of ions that serve as qubits in QC, from fewer than 100 available in state of the art linear ion-trap devices today to an order of 10 5 crystallized ions in the storage-ring setup. This new trap design unifies two different concepts: the storage rings of charged particles and the linear ion traps used for QC and mass spectrometry. In this paper, we use the language of particle accelerators to discuss the ion state and dynamics. We outline the differences between the above concepts, analyze challenges of the large ring with a revolving chain of ions, and propose goals for the research and development required to enable future quantum computers with 1000 times more qubits than available today. The challenge of creating such a large-scale quantum system while maintaining the necessary coherence of the qubits and the high fidelity of quantum logic operations is significant. Performing analog quantum simulations may be an achievable initial goal for such a device. Quantum calculations and simulations of complex quantum systems will move forward both the fundamental science and the applied research. Nuclear and particle physics, many-body quantum systems, lattice gauge theories, and nuclear structure calculations are just a few examples in which a large-scale quantum simulation system will become a very powerful tool to move forward our understanding of nature.

36 MATERIALS SCIENCE↗