Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational Physics”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Flux-pumped impedance-engineered broadband Josephson parametric amplifier

Broadband quantum-limited amplifiers play a critical role in the single-shot readout of superconducting qubits, but a popular implementation, the traveling wave parametric amplifier, involves a complex design and fabrication process. Here, we present a simple design for a Josephson parametric amplifier, using a lumped element resonator comprising a superconducting quantum interference device whose useful bandwidth is enhanced with an on-chip impedance-matching circuit. Additionally, we demonstrate a flux-coupling geometry that maximizes the coupling to the Josepson loop and minimizes spurious excitation of the amplifier resonant circuit. The amplifier, which operates in a flux-pumped mode, is demonstrated with a power gain of more than 20dB over a bandwidth of about 300MHz, where approximate noise measurements indicate quantum-limited pehrformance. A procedure is given for optimizing the bandwidth for this kind of amplifier, using a linearized circuit simulation while minimizing non-linearities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Design and simulation of a muon detector to characterize geological overburden

This study presents the design, construction, and simulation of a mobile muon detector tailored for geological overburden characterization. The detector employs plastic scintillator paddles with silicon photomultipliers (SiPMs) and a QuarkNet data acquisition system, offering a portable solution suitable for remote field deployment. The simulator’s modular aluminum frame allows for adjustable geometry and directional sensitivity, while its battery system supports over a week of autonomous operation. Preliminary experimental tests confirmed that its muon flux measurements were consistent with theoretical expectations. A comprehensive simulation framework using Geant4 and CORSIKA was developed to model detector response and overburden effects. Analytical and Monte Carlo methods were used to assess quadrant resolution and infer muon directionality. This work lays the foundation for future overburden mapping and supports the development of reconstruction algorithms for geological applications.

72 - PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Structure and assembly of the diiron cofactor in the heme-oxygenase–like domain of the N-nitrosourea–producing enzyme SznF

In biosynthesis of the pancreatic cancer drug streptozotocin, the tridomain nonheme-iron oxygenase SznF hydroxylates N δ and N ω ' of N ω -methyl- L -arginine before oxidatively rearranging the triply modified guanidine to the N-methyl-N-nitrosourea pharmacophore. A previously published structure visualized the monoiron cofactor in the enzyme’s C-terminal cupin domain, which promotes the final rearrangement, but exhibited disorder and minimal metal occupancy in the site of the proposed diiron cofactor in the N-hydroxylating heme-oxygenase–like (HO-like) central domain. We leveraged our recent observation that the N-oxygenating µ-peroxodiiron(III/III) intermediate can form in the HO-like domain after the apo protein self-assembles its diiron(II/II) cofactor to solve structures of SznF with both of its iron cofactors bound. These structures of a biochemically validated member of the emerging heme-oxygenase–like diiron oxidase and oxygenase (HDO) superfamily with intact diiron cofactor reveal both the large-scale conformational change required to assemble the O 2 -reactive Fe 2 (II/II) complex and the structural basis for cofactor instability—a trait shared by the other validated HDOs. During cofactor (dis)assembly, a ligand-harboring core helix dynamically (un)folds. The diiron cofactor also coordinates an unanticipated Glu ligand contributed by an auxiliary helix implicated in substrate binding by docking and molecular dynamics simulations. The additional carboxylate ligand is conserved in another N-oxygenating HDO but not in two HDOs that cleave carbon–hydrogen and carbon–carbon bonds to install olefins. Among ~9,600 sequences identified bioinformatically as members of the emerging HDO superfamily, ~25% conserve this additional carboxylate residue and are thus tentatively assigned as N-oxygenases.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Polyubiquitin ligand-induced phase transitions are optimized by spacing between ubiquitin units

Biomolecular condensates form via multivalent interactions among key macromolecules and are regulated through ligand binding and/or posttranslational modifications. One such modification is ubiquitination, the covalent addition of ubiquitin (Ub) or polyubiquitin chains to target macromolecules. Specific interactions between polyubiquitin chains and partner proteins, including hHR23B, NEMO, and UBQLN2, regulate condensate assembly or disassembly. Here, we used a library of designed polyubiquitin hubs and UBQLN2 as model systems for determining the driving forces of ligand-mediated phase transitions. Perturbations to either the UBQLN2-binding surface of Ub or the spacing between Ub units reduce the ability of hubs to modulate UBQLN2 phase behavior. By developing an analytical model based on polyphasic linkage principles that accurately described the effects of different hubs on UBQLN2 phase separation, we determined that introduction of Ub to UBQLN2 condensates incurs a significant inclusion energetic penalty. This penalty antagonizes the ability of polyUb hubs to scaffold multiple UBQLN2 molecules and cooperatively amplify phase separation. The extent to which polyubiquitin hubs promote UBQLN2 phase separation is encoded in the spacings between Ub units. This spacing is modulated by chains of different linkages and designed chains of different architectures, thus illustrating how the ubiquitin code regulates functionality via the emergent properties of the condensate. The spacing in naturally occurring linear polyubiquitin chains is already optimized to promote phase separation with UBQLN2. We expect our findings to extend to other condensates, emphasizing the importance of ligand properties, including concentration, valency, affinity, and spacing between binding sites in studies and designs of condensates.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Where did the tumor start? An inverse solver with sparse localization for tumor growth models

In this work, we present a numerical scheme for solving an inverse problem for parameter estimation in tumor growth models for glioblastomas, a form of aggressive primary brain tumor. The growth model is a reaction–diffusion partial differential equation (PDE) for the tumor concentration. We use a PDE-constrained optimization formulation for the inverse problem. The unknown parameters are the reaction coefficient (proliferation), the diffusion coefficient (infiltration), and the initial condition field for the tumor PDE. Segmentation of magnetic resonance imaging (MRI) scans drive the inverse problem where segmented tumor regions serve as partial observations of the tumor concentration. Like most cases in clinical practice, we use data from a single time snapshot. Moreover, the precise time relative to the initiation of the tumor is unknown, which poses an additional difficulty for inversion. We perform a frozen-coefficient spectral analysis and show that the inverse problem is severely ill-posed. We introduce a biophysically motivated regularization on the structure and magnitude of the tumor initial condition. In particular, we assume that the tumor starts at a few locations (enforced with a sparsity constraint on the initial condition of the tumor) and that the initial condition magnitude in the maximum norm is equal to one. We solve the resulting optimization problem using an inexact quasi-Newton method combined with a compressive sampling algorithm for the sparsity constraint. Our implementation uses PETSc and AccFFT libraries. We conduct numerical experiments on synthetic and clinical images to highlight the improved performance of our solver over a previously existing solver that uses standard two-norm regularization for the calibration parameters. The existing solver is unable to localize the initial condition. Our new solver can localize the initial condition and recover infiltration and proliferation. In clinical datasets (for which the ground truth is unknown), our solver results in qualitatively different solutions compared to the two-norm regularized solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A domain decomposition preconditioning for an inverse volume scattering problem

In this work, we propose domain decomposition preconditioners for the solution of an integral equation formulation of the acoustic forward and inverse scattering problems. We study both forward and inverse volume problems and propose preconditioning techniques to accelerate the iterative solvers. For the forward scattering problem, we extend the domain decomposition based preconditioning techniques presented for partial differential equations in Cai and Sarkis (1999 SIAM J. Sci. Comput. 21 792–797), to integral equations. We combine this domain decomposition preconditioner with a low-rank correction, which is easy to construct, forming a new preconditioner. For the inverse scattering problem, we use the forward problem preconditioner as a building block for constructing a preconditioner for the Gauss–Newton Hessian. We present numerical results that demonstrate the performance of both preconditioning strategies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Gibbs mixing of partially distinguishable photons with a polarising beamsplitter membrane

For a thought experiment concerning the mixing of two classical gases, Gibbs concluded that the work that can be extracted from mixing is determined by whether or not the gases can be distinguished by a semi-permeable membrane; that is, the mixing work is a discontinuous function of how similar the gases are. Here we describe an optomechanical setup that generalises Gibbs’ thought experiment to partially distinguishable quantum gases. Specifically, we model the interaction between a polarisation dependent beamsplitter, that plays the role of a semi-permeable membrane, and two photon gases of non-orthogonal polarisation. We find that the work arising from the mixing of the gases is related to the potential energy associated with the displacement of the microscopic membrane, and we derive a general quantum mixing work expression, valid for any two photon gases with the same number distribution. The quantum mixing work is found to change continuously with the distinguishability of the two polarised gases. In addition, fluctuations of the work on the microscopic membrane become important, which we calculate for Fock and thermal states of the photon gases. Our findings generalise Gibbs’ mixing to the quantum regime and open the door for new quantum thermodynamic (thought) experiments with quantum gases with non-orthogonal polarisations and microscopic pistons that can distinguish orthogonal polarisations.

distinguishability↗
DOE OSTI2024

Transport and losses of energetic particles in tokamaks in the presence of Alfvén activity using the new full orbit TAPaS code coupled to FAR3d

Recent developments and tools integrated into the TAPaS code are presented, enabling realistic scenario simulations of particle dynamics within experimental tokamak magnetic equilibria. In particular, the enhanced capabilities of TAPaS enable seamless coupling with external simulations, provided the metric and equilibrium magnetic field of the external code are known. Coupling TAPaS with the gyro-fluid code FAR3d, the transport and losses of energetic particles in the presence Alfvén eigenmodes (AEs) in DIII-D plasma discharge #159243 were investigated. Detailed analyses of prompt losses with and without collisions were performed. Then, further analysis was performed in the presence of electromagnetic perturbations resulting from AEs activity. The results indicate that, for the energies and the initial conditions considered here, the presence of AEs enhances the particle losses.

Alfven eigenmodes↗
DOE OSTI2025

LTAU-FF: Loss Trajectory Analysis for Uncertainty in atomistic Force Fields

Model ensembles are effective tools for estimating prediction uncertainty in deep learning atomistic force fields. However, their widespread adoption is hindered by high computational costs and overconfident error estimates. In this work, we address these challenges by leveraging distributions of per-sample errors obtained during training and employing a distance-based similarity search in the model latent space. Our method, which we call LTAU (Loss Trajectory Analysis for Uncertainty), efficiently estimates the full probability distribution function of errors for any test point using the logged training errors, achieving speeds that are 2–3 orders of magnitudes faster than typical ensemble methods and allowing it to be used for tasks where training or evaluating multiple models would be infeasible. We apply LTAU towards estimating parametric uncertainty in atomistic force fields (LTAU-FF), demonstrating that it produces well-calibrated confidence intervals and predicts errors that correlate strongly with the true errors for data near the training domain. Furthermore, we show that the errors predicted by LTAU-FF can be used in practical applications for detecting out-of-domain data, tuning model performance, and predicting failure during simulations. We believe that LTAU will be a valuable tool for uncertainty quantification in atomistic force fields and is a promising method that should be further explored in other domains of machine learning.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Thresholded quantum sensing with a frustrated Kitaev trimer

We investigate the response of a Ramsey interferometric quantum sensor based on a frustrated, three-spin system (a Kitaev trimer) to a classical time-dependent field (signal). The system eigenspectrum is symmetric about a critical point, |𝒃| = 0, with four of the spectral components varying approximately linearly with the magnetic field and four exhibiting a nonlinear dependence. Under the adiabatic approximation and for appropriate initial states, we show that the sensor's response to a zero-mean signal is such that below a threshold, |𝒃| <⁢ 𝑏 th , the sensor does not respond to the signal, whereas above the threshold, the sensor acts as a detector that the signal has occurred. This thresholded response is approximately omnidirectional. Moreover, when deployed in an entangled multisensor configuration, the sensor achieves sensitivity at the Heisenberg limit. Such detectors could be useful both as stand-alone units for signal detection above a noise threshold and in two- or three-dimensional arrays, analogous to a quantum bubble chamber, for applications such as particle track detection and long-baseline telescopy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Improving qubit readout with hidden Markov models

We demonstrate the application of pattern recognition algorithms via hidden Markov models (HMM) for qubit readout. This scheme provides a state-path trajectory approach capable of detecting qubit-state transitions and makes for a robust classification scheme with higher starting-state assignment fidelity than when compared to a multivariate Gaussian or a support vector machine scheme. Therefore, the method also eliminates the qubit-dependent readout time optimization requirement in current schemes. Using a HMM state discriminator we estimate fidelities reaching the ideal limit. Unsupervised learning gives access to transition matrix, priors, and IQ distributions, providing a toolbox for studying qubit-state dynamics during strong projective readout.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Extending quantum probabilistic error cancellation by noise scaling

Here, we propose a general framework for quantum error mitigation that combines and generalizes two techniques: probabilistic error cancellation (PEC) and zero-noise extrapolation (ZNE). Similar to PEC, the proposed method represents ideal operations as linear combinations of noisy operations that are implementable on hardware. However, instead of assuming a fixed level of hardware noise, we extend the set of implementable operations by noise scaling. By construction, this method encompasses both PEC and ZNE as particular cases and allows us to investigate a larger set of hybrid techniques. For example, gate extrapolation can be used to implement PEC without requiring knowledge of the device’s noise model, e.g., avoiding gate-set tomography. Alternatively, probabilistic error reduction can be used to estimate expectation values at intermediate virtual noise strengths (below the hardware level), leading to partially mitigated results at a lower sampling cost. Moreover, multiple results obtained with different noise-reduction factors can be further postprocessed with ZNE to better approximate the zero-noise limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Possibilistic simulation of quantum circuits by classical circuits

In a breakthrough work, Bravyi, Gosset, and König (BGK) [Science 362, 308 (2018)] unconditionally proved that constant-depth quantum circuits are more powerful than their classical counterparts. Here, their result is equivalent to saying that a particular family of constant-depth quantum circuits takes classical circuits at least Ω⁡(log⁡n) depth to “simulate,” in a certain sense. In our paper, we formalize their sense of simulation, which we call “possibilistic simulation” or “p-simulation,” and construct explicit classical circuits that can p-simulate any depth-d quantum circuit with Clifford and t T-gates in depth O(d+t). Our classical circuits use {NOT,AND,OR} gates of fan-in ≤ 2.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Perturbative model of noisy quantum signal processing

Recent progress in quantum signal processing (QSP) and its generalization, quantum singular value transformation, has led to a grand unification of quantum algorithms. However, inherent experimental noise in quantum devices severely limits the length of realizable QSP sequences. Here, we consider a model of QSP with generic perturbative noise in the signal processing basis and present a diagrammatic notation useful for analyzing such errors. To demonstrate our technique, we study a specific coherent error, that of under- or overrotation of the signal processing operator parametrized by ε<<1. For this coherent error model, it is shown that while Pauli Z errors are not recoverable without additional resources, Pauli X and Y errors can be arbitrarily suppressed by coherently appending a noisy recovery QSP without the use of additional resources or ancillas. Furthermore, through a careful accounting of errors using our diagrammatic tools, we provide an upper and lower bound on the length of this recovery QSP operator. We anticipate that the perturbative technique and the diagrammatic notation proposed here will facilitate future study of generic noise in QSP and quantum algorithms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Exponential Escape Rate of Filamentary Incubation in Mott Spiking Neurons

Mott materials such as vanadium oxides, when subject to a strong applied voltage, present an inhomogeneous insulator-to-metal transition with formation of metallic filaments within the insulating bulk. This property is enabling the development of compact and power-efficient neuromorphic devices known as Mott neurons. However, the nature of the transition has not been fully understood yet, as it may be attributed to different effects, including Joule self-heating and hot-carrier injection. Moreover, the experimental determination of the threshold voltage needed to induce the transition has proven to be challenging, as the transition becomes increasingly unpredictable when the threshold is approached. The physical understanding of these issues would not only deepen our understanding of Mott insulators, but would also be an important step toward the realization of neuromorphic devices based on such materials. In this work we use numerical simulations based on the Mott resistor network model to study the nature of the filament incubation and formation process. Here, we show that both electronic and thermal effects, in the form of current density focusing and Joule self-heating, respectively, contribute to the filamentary incubation and growth. Remarkably, we find that the percolation of the metallic filaments near the threshold is intrinsically stochastic, qualitatively similar to the familiar Arrhenius activated behavior and to the stochastic firing of biological neurons. More precisely, we characterize the filament percolation as a Poisson point process, which has the same probability distribution as mathematical models of neuronal firing with an exponential escape rate. Finally, we support the numerical simulation results by performing experiments in VO 2 that are in agreement with the exponential escape rate behavior. Thus, we establish a functionality of Mott insulators that opens a path toward implementing neuromorphic hardware with quantum materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Efficient Entanglement of Spin Qubits Mediated by a Hot Mechanical Oscillator

Localized electronic and nuclear spin qubits in the solid state constitute a promising platform for storage and manipulation of quantum information, even at room temperature. However, the development of scalable systems requires the ability to entangle distant spins, which remains a challenge today. We propose and analyze an efficient, heralded scheme that employs a parity measurement in a decoherence free subspace to enable fast and robust entanglement generation between distant spin qubits mediated by a hot mechanical oscillator. We find that high-fidelity entanglement at cryogenic and even ambient temperatures is feasible with realistic parameters and show that the entangled pair can be subsequently leveraged for deterministic controlled-NOT operations between nuclear spins. Our results open the door for novel quantum processing architectures for a wide variety of solid-state spin qubits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗