Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “differential mode”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

On the relationship between eastern China aerosols and western North Pacific tropical cyclone activity

How aerosols affect the weather and climate system has received increasing attention. This study finds a significant negative correlation between March–May eastern China aerosol optical depth (AOD) and July–November western North Pacific (WNP) tropical cyclone (TC) frequency during 2003–2020. This time period spans when several aerosol reanalyses are available and both Terra and Aqua Moderate Resolution Imaging Spectroradiometer AOD retrievals are assimilated therein. Composite analyses and budget analyses of dynamical genesis potential indices indicate the importance of large-scale environmental factors, especially vertical velocity and vertical wind shear, associated with changes in AOD that in turn modulate changes in WNP TC frequency. Increased eastern China AOD may facilitate negative Pacific meridional mode development via modulation of the westerly jet, which then forces an anticyclonic circulation over the WNP basin. Increased AOD can also directly decrease the inter-hemispheric temperature differential and increases the intra-hemispheric temperature differential between the equator and the mid-latitudes, thus weakening ascending motion and enhancing vertical wind shear over the WNP, especially the southeastern portion of the basin. Further, all of these large-scale environment changes induced by increased eastern China AOD tend to suppress WNP TCs. This study highlights the potential influence of eastern China aerosol loadings on WNP TCs, thus improving our understanding of TC climate variability over the WNP.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

On the mechanism of ionization oscillations in Hall thrusters

Low-frequency ionization oscillations involving plasma and neutral density (breathing modes) are the most violent perturbations in Hall thrusters for electric propulsion. Because of its simplicity, the zero-dimensional (0D) predator–prey model of two nonlinearly coupled ordinary differential equations for plasma and neutral density has often been used for the characterization of such oscillations and scaling estimates. In this work, we investigate the properties of its continuum analog, the one-dimensional (1D) system of two nonlinearly coupled equations in partial derivatives (PDEs) for plasma and neutral density. This is a more general model, of which the standard 0D predator–prey model is a special limit case. We show that the 1D model is stable and does not show any oscillations for the boundary conditions relevant to Hall thrusters and the uniform ion velocity. We then propose a reduced 1D model based on two coupled PDEs for plasma and neutral densities that is unstable and exhibit oscillations if the ion velocity profile with the near-the-anode back-flow (toward the anode) region is used. Comparisons of the reduced model with the predictions of the full model that takes into account the self-consistent plasma response show that the main properties of the breathing mode are well captured. In particular, it is shown that the frequency of the breathing mode oscillations is weakly dependent on the final ion velocity but shows a strong correlation with the width of the ion back-flow region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Analytic techniques for solving the transport equations in electroweak baryogenesis

We develop an efficient method for solving transport equations, particularly in the context of electroweak baryogenesis. It provides fully-analytical results under mild approximations and can also test semi-analytical results, which are applicable in more general cases. Key elements of our method include the reduction of the second-order differential equations to first order, representing the set of coupled equations as a block matrix of the particle densities and their derivatives, identification of zero modes, and block decomposition of the matrix. We apply our method to calculate the baryon asymmetry of the Universe (BAU) in a Standard Model effective field theory framework of complex Yukawa couplings to determine the sensitivity of the resulting BAU to modifications of various model parameters and rates, and to estimate the effect of the commonly-used thin-wall approximation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Canonical and noncanonical Hamiltonian operator inference

Here, a method for the nonintrusive and structure-preserving model reduction of canonical and noncanonical Hamiltonian systems is presented. Based on the idea of operator inference, this technique is provably convergent and reduces to a straightforward linear solve given snapshot data and gray-box knowledge of the system Hamiltonian. Examples involving several hyperbolic partial differential equations show that the proposed method yields reduced models which, in addition to being accurate and stable with respect to the addition of basis modes, preserve conserved quantities well outside the range of their training data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A Theory of Ultrafast Charge Transfer Relaxation with Non-Innocent Solvent Molecules

In this work we revisit the photodynamics of tetracyanoethylene-hexamethylbenzene (TCNE- HMB), the molecular complex studied by Hochstrasser et al. [J. Chem. Phys. 100, 4797–4810, 1994] that has long challenged the applicability of Marcus’s theory of elec- tron transfer for predicting photochemical reactions. Using a novel black-box electronic structure algorithm (time-dependent density functional theory with one double, TD- DFT-1D) to efficiently run molecular dynamics that can treat charge recombination, we run ab initio surface hopping molecular dynamics and confirm that, for a polar solvent, charge recombination rates can be incredibly fast (indeed faster than the sol- vent relaxation time); for non-polar solvents, the rate is much slower. We demonstrate that, although Marcus theory cannot be directly applied, these nonequilibrium (and sometimes incredibly fast) photoexcited dynamics can be effectively explained within a two-state model without any evidence of a transition through a conical intersection. Most importantly, for this paradigmatic model system, we are able to identify two nuclear coordinates of interest (rather than the single coordinate predicted by Marcus or a full set of internal quantum modes studied by Bixon and Jortner): the solvent relaxation in the first shell (that strongly modulates the energies of the charge trans- fer state and differentiates time scales for relaxation) and a nuclear displacement in the TCNE-HMB complex arising from a handful of vibrations that induces non-Born Oppenheimer motion and eventually facilitates an abrupt electronic transition to the ground state. Altogether, these findings suggest a tractable generalization of Marcus theory for future simulations of photochemistry with non-innocent solvent environ- ments in the spirit of a Hamiltonian suggested by Stuchebrukhov (J. Chem. Phys. 107, 3821, 1997).

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

Instability and turbulent relaxation in a stochastic magnetic field

An analysis of instability dynamics in a stochastic magnetic field is presented for the tractable case of the resistive interchange. Externally prescribed static magnetic perturbations convert the eigenmode problem to a stochastic differential equation, which is solved by the method of averaging. The dynamics are rendered multi-scale, due to the size disparity between the test mode and magnetic perturbations. Maintaining quasi-neutrality at all orders requires that small-scale convective cell turbulence be driven by disparate scale interaction. Here, the cells in turn produce turbulent mixing of vorticity and pressure, which is calculated by fluctuation-dissipation type analyses, and are relevant to pump-out phenomena. The development of correlation between the ambient magnetic perturbations and the cells is demonstrated, showing that turbulence will 'lock on' to ambient stochasticity. Magnetic perturbations are shown to produce a magnetic braking effect on vorticity generation at large scale. Detailed testable predictions are presented. The relations of these findings to the results of available simulations and recent experiments are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Experimental investigation of beam-ion losses induced by magnetic perturbations using the light ion beam probe technique in the ASDEX Upgrade tokamak

The impact of externally applied magnetic perturbations (MPs) on fast-ion losses has been investigated by means of the light ion beam probe (LIBP) technique in the ASDEX Upgrade tokamak. The LIBP technique allows to experimentally infer the fast-ion orbit displacement induced by MPs via first-orbit losses using scintillator based fast-ion loss detector (FILD) measurements. The fast-ion orbit displacement against different applied MP spectra has been studied. These shots were conducted in ELM mitigated H-mode plasmas. A rigid rotation of the MP coils was applied with a frequency of 1 Hz, with an n = 2 configuration and changing the differential phase between the upper and lower set of coils (ΔΦ ul ) on a shot-to-shot basis. Beam sources Q7 (tangential) and Q8 (radial) were used to probe different fast-ion orbits with FILD1. The measured fast-ion orbit displacement ranges from 3 to 20 mm approximately, and no qualitative difference is observed between ions from beam sources Q7 and Q8. The minimum is found for a ΔΦ ul ~ 50°, which is shifted with respect to the minimum of the plasma boundary displacement, found at ΔΦ ul ~ 0°. A first attempt to validate the orbit following code ASCOT—including the plasma response calculated with the MARS-F code—against these experimental measurements is performed. Finally, while the dependence of the first-orbit fast-ion displacement with ΔΦ ul does not match the experimental measurements, these simulations do capture other features such as the order of magnitude of the orbit displacement and the importance of the toroidal spectrum of the applied perturbation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Optimizing the differential connection schemes for detecting 3D magnetic perturbations in DIII-D

Optimizing the differential pair connection scheme (i.e. the set of pairs) of a toroidal array of magnetic sensors dedicated to measuring slowly rotating asymmetric fields can enhance the mode number detection capability and failure-resilience. In this work the condition number obtained from singular value decomposition (SVD) of the design matrix is used as a metric to evaluate the quality of a connection scheme. A large number of possible pair connections are usually available, so evaluating all of them may require extensive use of computational resources and can be very time consuming. Alternative methods to reduce the number of pairs evaluated without losing the capabilities of toroidal mode detection are presented in this paper. Here, three examples of the applications of such analysis for the 3D magnetic diagnostic system of DIII-D are also presented: the addition of two new toroidal arrays with n > 3 detection capabilities, the modification of an existing toroidal array in the low field side of the machine to accommodate the addition of a helicon antenna, and the design of changes in several toroidal arrays in the high field side to accommodate the addition of a lower hybrid current drive antenna on the center post.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Magnetic Field Strength from Turbulence Theory. I. Using Differential Measure Approach

The mean plane-of-sky magnetic field strength is traditionally obtained from the combination of polarization and spectroscopic data using the Davis–Chandrasekhar–Fermi (DCF) technique. However, we identify the major problem of the DCF technique to be its disregard of the anisotropic character of MHD turbulence. On the basis of the modern MHD turbulence theory we introduce a new way of obtaining magnetic field strength from observations. Unlike the DCF technique, the new technique uses not the dispersion of the polarization angle and line-of-sight velocities, but increments of these quantities given by the structure functions. To address the variety of astrophysical conditions for which our technique can be applied, we consider turbulence in both media with magnetic pressure higher than the gas pressure, corresponding, e.g., to molecular clouds, and media with gas pressure higher than the magnetic pressure, corresponding to the warm neutral medium. We provide general expressions for arbitrary admixtures of Alfvén, slow, and fast modes in these media and consider in detail particular cases relevant to diffuse media and molecular clouds. We successfully test our results using synthetic observations obtained from MHD turbulence simulations. We demonstrate that our differential measure approach, unlike the DCF technique, can be used to measure the distribution of magnetic field strengths, can provide magnetic field measurements with limited data, and is much more stable in the presence of induced large-scale variations of nonturbulent nature. Furthermore, our study uncovers the deficiencies of earlier DCF research.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Double-differential Measurements of Mesonless Charged-current Muon Neutrino Interactions on Argon with Final-State Protons using the MicroBooNE Detector

This note reports the first double-differential measurements of charged-current $ν_µ$ scattering on argon leading to final states containing zero mesons and one or more protons. This event topology (hereafter abbreviated as CC0$πNp$) is the most common at the neutrino energies typically produced by the Fermilab Booster Neutrino Beam. A detailed understanding of neutrino-argon scattering in the CC0$πNp$ channel is therefore crucial for the success of the precision neutrino oscillation analyses planned for the Short-Baseline Neutrino (SBN) program. This remains true for the upcoming Deep Underground Neutrino Experiment (DUNE), but the higher mean neutrino energy used there will ensure that more inelastic reaction modes, such as single pion production, will also play a major role. The analysis described in this note builds on a previous MicroBooNE study of the CC0$πNp$ channel which obtained the first single-differential cross-section measurements on an argon target. Since that foundational work, significant improvements have been made to MicroBooNE’s simulation software, event reconstruction algorithms, and procedure for calculating systematic uncertainties. When combined with a larger dataset (corresponding to a beam exposure of 6.79 × 10 20 protons-on-target versus 1.60 × 10 20 in Ref. [3]), these enhancements allow the important CC0$πNp$ channel to be studied in more detail. This note begins with a description of the data and simulation samples used as input to the analysis. Section 3 then defines the CC0$πNp$ signal event topology, and Section 4 describes a set of selection criteria designed to identify these events in MicroBooNE data. Binning schemes are then defined in Section 5 for two double-differential measurements of event rates. The first of these considers the momentum and scattering cosine of the outgoing muon, while the second reports the same observables for the leading proton, i.e., the final-state proton with the largest momentum. After a discussion of systematic uncertainties in Section 6, the note concludes by comparing the predictions of MicroBooNE Monte Carlo (MC) simulations to the measured double-differential distributions. These results will form the basis for a future extraction of flux-averaged double-differential CC0$πNp$ cross sections that will be immediately comparable to the theoretical predictions of multiple neutrino event generators. The selection described herein may also be used to study various other observables in CC0$πNp$ events, including those which are sensitive to correlations between leptonic and hadronic kinematics in the final state.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Calorimetric study of skutterudite (CoAs2.92) and heazlewoodite (Ni3S2)

Abstract Nickel and cobalt arsenides, sulfarsenides, and sulfides occur in many hydrothermal ore deposits, but their thermodynamic properties are not well known, in some cases not known at all. In this work, we determined a full set of thermodynamic properties for heazlewoodite and skutterudite. Both phases were synthesized in evacuated silica tubes at elevated temperatures, and electron microprobe analyses gave their compositions as Ni3S2 and CoAs2.92, respectively. Enthalpies of formation were measured by high-temperature oxide-melt solution calorimetry. The reference phases were pure elements, thus eliminating any systematic errors related to such phases. The enthalpies of formation at T = 298.15 K and P = 105 Pa are –216.0 ± 8.4(2σ) and –88.2 ± 6.1 kJ·mol−1 for Ni3S2 and CoAs2.92, respectively. Entropies were calculated from low-temperature heat capacity (CP) data from relaxation (PPMS) calorimetry and are 133.8 ± 1.6 and 106.4 ± 1.3 J·mol–1·K–1, respectively. The calculated Gibbs free energies of formation are –210.0 ± 8.4 and –79.9 ± 6.2 kJ·mol−1 for Ni3S2 and CoAs2.92, respectively. The PPMS CP data, together with a set of differential scanning calorimetry measurements, were used to derive CP polynomials up to 700 K with the Kieffer model based on previously published frequencies of acoustic and optic modes. Equilibrium constants for selected reactions with an aqueous phase were calculated up to 700 K. Geochemical modeling in these systems, however, should await until more reliable data for other phases from the system Co-Ni-As-S are available.

Geochemistry & Geophysics↗
DOE OSTI2021

Undersizing of aged African biomass burning aerosol by an ultra-high-sensitivity aerosol spectrometer

The ultra-high-sensitivity aerosol spectrometer (UHSAS) differs from most other optical particle spectrometers by using a high-power infrared (IR) laser to detect small particles and reduce the sizing ambiguity due to the non-monotonicity of scattering with particle size. During the NASA ORACLES project (ObseRvations of Aerosols above CLouds and their intEractionS) over the southeast Atlantic Ocean, the UHSAS clearly undersized particles in the biomass burning plume extending from southern Africa. Since the horizontal and vertical extent of the plume was vast, the NASA P-3B research aircraft often flew through a fairly uniform biomass burning plume for periods exceeding 30 min, sufficient time to explore the details of the UHSAS response by selecting single particle sizes with a differential mobility analyzer (DMA) and passing them to the UHSAS. This was essentially an in-flight calibration of the UHSAS using the particles of interest. Two modes of responses appeared. Most particles were undersized by moderate amounts, ranging from not at all for 70 nm aerosols to 15 % for 280 nm particles. Mie scattering calculations show that composition-dependent refractive index of the particles cannot explain the pattern. Heating of brown carbon or tarballs in the beam causing evaporation and shrinking of the particles is the most plausible explanation, though mis-sizing due to non-sphericity cannot be ruled out. A small fraction (10 %–30 %) of the particles were undersized by 25 % to 35 %. Those were apparently the particles containing refractory black carbon. Laboratory calibrations confirm that black carbon is drastically undersized by the UHSAS, because particles heat to their vaporization point and shrink. A simple empirical correction equation was implemented that dramatically improves agreement with DMA distributions between 100 and 500 nm. It raised the median particle diameter by 18 nm, from 163 to 181 nm, during the August 2017 deployment and by smaller amounts during deployments with less intense pollution. Calculated scattering from UHSAS size distributions increased by about 130 %, dramatically improving agreement with scattering measured by nephelometers. The correction is only valid in polluted instances; clean marine boundary layer and free troposphere aerosols behaved more like the calibration spheres. We were unable to directly test the correction between 500 and 1000 nm, though aerodynamic particle sizer (APS) data appear to show that the correction is poor at the largest diameters, which is no surprise as the composition of those particles is likely to be quite different than that of the accumulation mode. This adds to the evidence that UHSAS data must be treated cautiously whenever the aerosol may absorb infrared light. Similar corrections may be required whenever brown carbon aerosol is present.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

First measurement of polarized spin-density matrix elements and differential cross sections dσ/dt in ω photoproduction off the proton for 2.7 < Eγ < 5.2 GeV using CLAS at Jefferson Lab

We report on the differential cross sections dσ/dt, the unpolarized spin-density matrix elements ρ 00 0 , ρ 1 − 1 0 , Re ρ 10 0 , and the first extraction of the polarized elements Im ρ 10 3 , Im ρ 1 − 1 3 for the reaction γp → pω using the CLAS spectrometer at Jefferson Laboratory. The ω mesons were detected in their dominant charged decay mode, ω → π + π − π 0 , and all t-dependent results are presented in a fine binning for incident photon energies between 2.73 and 5.16 GeV (corresponding to the center-of-mass energy range W ∈ [ 2.45, 3.25 ] GeV). All matrix elements are first measurements for − t > 0.6 GeV2. Moreover, differential cross sections dσ/d(cos Θ c . m . ω ) and the corresponding angle-dependent unpolarized spin-density matrix elements in the Adair frame are presented for the incident photon energy range 1.56–3.80 GeV (corresponding to W ∈ [ 1.95, 2.83 ] GeV). These new ω photoproduction data are consistent with earlier CLAS results but extend the energy range well beyond the nucleon resonance region into the Regge regime. The comparison with Regge-theory-based model predictions shows that the new data impose more stringent constraints on our understanding of ω photoproduction.

Hu, T.↗
DOE OSTI2024

Correlation Between Microscopic Current Fluctuations Observed at Ultra-Microelectrodes and Macroscopic Bulk Electrolysis Performance in Redox-Active Microemulsions

Microemulsions (μEs) have been proposed as redox flow battery (RFB) electrolytes that maximize ionic conductivity and charge capacity by synergizing two immiscible phases. However, charge transfer during electrolysis in μEs is poorly understood. Here, we show that ultramicroelectrode electrolysis of ferrocene-loaded μEs –20%, 60%, and 90% water - reveals stochastic current fluctuations. These are differentiated in the scanning electrochemical microscopy (SECM) geometry, where power spectral density analysis showed distinct changes in the frequency contributions. SECM in the substrate generation-tip collection mode showed that fluctuations arise under mass-transfer control. Significant differences in the diffusion coefficient of ferrocene species were deducted from SECM approach curves, suggesting phase transfer behavior. Using bulk electrolysis, we calculated the charge accessibility and cycling behavior in the μEs. A decrease in the stochastic behavior of the μEs seems to correlate to a higher accessibility and cycling performance, with the 90% water μE displaying the best reversibility and the 60% the lowest. Altogether, these results suggest that Marangoni-type convection driven by concentration gradients and/or μE restructuring during charge transfer play a role in the electrochemical performance of μEs. This presents opportunities for screening and diagnosing the performance of these emerging RFB electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2024

Tearing Stable Stationary ITER Baseline Operation in DIII-D

We present a tearing stable, stationary and reproducible operating solution for DIII-D plasmas characterized by the ITER Baseline Scenario normalized parameter set and shape. In these plasmas low differential rotation (Δf1,2) between the core and edge is identified as the single root cause of the onset of the primary limiting instability, the 2,1 Neoclassical Tearing Mode (NTM), while a group of other examined parameters have a much weaker impact on the stability. Explanation is offered by drift kinetic simulations, which show a six-fold reduction in the NTM onset threshold due to diminished stabilizing ion polarization currents when the seed island drift frequency ceases in the local plasma frame. Based on the experimental observation and the explanation provided by the theory, a new method is developed to sustain Δf1,2 by modifying the edge neoclassical potential through neutral gas fueling. This method enables stationary ITER baseline operation free of disruptive 2,1 tearing modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Influence of overhead HVAC and aerosol control strategies on coarse mode particle dispersion and exposure in a full-scale room experiment

Coarse mode respiratory aerosols can carry viral loads over long distances and have very different dynamics than submicron particles, but experimental studies under realistic conditions remain limited. Here, to study the differential impacts on exposure under different mixing conditions, we co-released 7–10 µm particles and carbon dioxide (CO 2 )—which served as an indicator of gas and submicron particle dynamics—in a 158 m 3 room at LBNL's FLEXLAB facility with an overhead heating, ventilation, and air conditioning (HVAC) system. The room was arranged as a distanced meeting then a classroom with eight heated manikins and a researcher. Spatial variability was measured using 16 particle counters and 22–26 CO 2 sensors throughout the space. Conditions included: HVAC off or supply air at 1000-1060 m 3 h -1 at neutral, cooling, or heating temperatures; with and without 20% outdoor air; and added HVAC filtration, portable air cleaners (PACs), or a physical barrier between the speaker and occupants. We found that good mixing via neutral or cooling supply air or use of PACs under heating lowered coarse particle exposure at some locations, but increased exposure for one-quarter to two-thirds of manikins compared to poor mixing under heating. A physical barrier reduced direct transfer of coarse particles during heating, but less during cooling. High spatial variability shows that a single measurement cannot represent occupant exposure. Instantaneous air mixing assumptions overstate the effectiveness of ventilation, HVAC filtration, and upper-room germicidal ultraviolet disinfection for coarse particles, as relatively few particles reach the return grille or upper room under most conditions.

Air mixing↗
DOE OSTI2022

Concealed Cathode Degradation in Lithium-Ion Cells with a Ni-Rich Oxide

Difficulties with sourcing cobalt and the interest in increasing cell energy have motivated the development of Ni-rich oxide materials for lithium-ion battery cathodes. Despite the intense research on the topic, there is limited information available on the long-term performance of novel cathode formulations. Here, we evaluate the stability of LiNi 0.9 Mn 0.05 Co 0.05 O 2 in full-cells tested for over five months, assessing how cycling, voltage and electrolyte additives impact cathode aging. We use differential voltage analysis to extract insights about the cathode from the full-cell data and identify the isolation of cathode particles in the delithiated (charged) state as a relevant mode of aging. Importantly, this particular mechanism of cathode aging does not cause immediate cell capacity fade, causing the simple analysis of cell capacity retention to overestimate the cathode stability under the investigated conditions. So, our observations serve as cautionary tale indicating that careful analysis of data from extended testing may be required for assessing the performance of Ni-rich cathodes and for evaluating how these materials are affected by electrolyte additives.

25 ENERGY STORAGE↗
DOE OSTI2025

Solving the Bernstein-Vazirani problem using Majorana-based topological quantum algorithms

Executing quantum algorithms using Majorana zero modes—a major milestone for the field of topological quantum computing—requires a platform that can be scaled to large quantum registers, can be controlled in real time and space, and a braiding protocol that uses the unique properties of these exotic particles. Here, we demonstrate the first successful simulation of a Majorana-based, fault-tolerant quantum algorithm to solve the Bernstein-Vazirani problem in two-dimensional magnet-superconductor hybrid structures from initialization to read-out of the final many-body state. Utilizing the Majorana zero modes’ topological properties, we introduce an optimized braiding protocol for the algorithm and a scalable architecture for its implementation with an arbitrary number of qubits. We visualize the algorithm protocol in real time and space by computing the non-equilibrium density of states, which is proportional to the time-dependent differential conductance, and the non-equilibrium charge density, which assigns a unique signature to each final state of the algorithm.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗