Search NASA⌕ Search

SEARCH · Search NASA

Results for “operator exponential”

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 253 records · Page 14

Statistical Correlation of Heliostat Pointing Deviation With Wind

This work was carried out as part of the Heliostat Consortium (HelioCon) Field Deployment subtask with the aim to develop a reduced order model framework for correlating wind speed and pointing deviation of a heliostat facet. There are only sparse field measurements of heliostat pointing deviations and accompanying wind conditions published in the literature. Heliostat test standards, such as IEC 62862-4-3, propose a suite of tests including laser pointing repeatability at wind speeds below 4 m/s, and provide technical requirements for heliostat slope and tracking deviations in coarse average wind speed bins of 4 m/s, 6 m/s, and 8 m/s. In addressing the gap of the variation of heliostat pointing deviation with wind speed, field measurements of laser pointing on a grid target and wind conditions were analyzed in this study at the Third-Party Metrology Platform at the National Laboratory of the Rockies (NLR) Flatirons Campus. Horizontal pointing deviations were found to follow a logarithmic relationship with peak wind speed, whereas vertical pointing deviations follow an exponential relationship with peak wind speed. Both horizontal and vertical pointing deviations also follow a second order polynomial relationship, as expected from the proportionality of elastic loads and deformations with the square of wind speed. The results indicate that heliostat facet pointing deviations in the vertical direction increase at a faster rate than in the horizontal direction with increasing wind speed over the tested range, however these are dependent on the heliostat structural design. Next steps are recommended for additional field measurements to confirm a linear relationship of pointing deviation with applied moment on a heliostat facet, and to distinguish between gravity-induced and wind-induced pointing deviations at different elevation angles. The derived correlations in the preliminary analysis in this report serve as a case study for heliostat developers and plant operators to estimate the wind-induced pointing deviations and their variation with peak gust wind speed. Next steps in future work would recommend higher resolution and longer duration datasets for different elevation angles and wind directions to reduce uncertainties and variance of collected laser beam spot data and their correlations with bin-averaged wind speed.

17 WIND ENERGY↗

Description of Sensor Assignment Optimization Method as Deployed on a Multi-Node Cluster

Data analytic methods are being developed to address the problem of how to assign a sensor set in a nuclear facility such that a requisite level of process monitoring capability is realized and that the sensor set is sufficiently rich to determine the status of the individual sensors with respect to need for calibration. There is an awareness in the nuclear industry that data analytics combined with rich sensor sets represent a means to improve operations and reduce costs. In the industry the calibration problem has been previously approached as an empirical data-driven problem with several methods having been developed. However, the experience of the utilities over the past ten years with these methods indicates that the absence of physics-based information renders the data-driven approach less reliable. Complicating factors such as the inherent variability of operation (both equipment alignment and operating condition) can confound a pure data-driven approach while there are no rigorous guidelines for determining what constitutes an adequate sensor set. The solution under development to overcome these shortcomings supplements the data analytic method with process information in a so-called process-constrained data-analytic approach. Simple balance equations are written for generic components (e.g., mechanical pump, valve, and heat exchanger). These do not require a priori knowledge of process parameters, such as heat transfer coefficients or friction factors. All that is needed on the part of the utility user is to identify the components and how they are connected. This report describes the development of a parallel computing capability for determining the optimal sensor set. The optimal sensor set problem suffers from the curse of dimensionality. Computation time increases exponentially as the size of the system grows. To overcome this difficulty a pre-conditioner algorithm is developed to find an approximate solution close the actual solution. This serves as a seed for the full-blown algorithm and acts to constrain the space that must searched. The optimization algorithms are described and the implementation on a parallel computing platform is described. The application of the method to a use case we are solving in collaboration with our utility partner served to illustrate how the default sensor set in a nuclear plant may not provide sufficient coverage to infer sensor calibration status.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The DESC stellarator code suite. Part 1. Quick and accurate equilibria computations

Three-dimensional equilibrium codes are vital for stellarator design and operation, and high-accuracy equilibria are also necessary for stability studies. This paper details comparisons of two three-dimensional equilibrium codes: VMEC, which uses a steepest-descent algorithm to reach a minimum-energy plasma state, and DESC, which minimizes the magnetohydrodynamic (MHD) force error in real space directly. Accuracy as measured by satisfaction of MHD force balance is presented for each code, along with the computation time. It is shown that DESC is able to achieve more accurate solutions, especially near axis. The importance of higher-accuracy equilibria is shown in DESC's better agreement of stability metrics with asymptotic formulae. DESC's global Fourier–Zernike basis also yields solutions with analytic derivatives explicitly everywhere in the plasma volume, provides improved accuracy in the radial direction versus conventional finite differences and allows for exponential convergence. Further, DESC can compute a solution with the same accuracy as VMEC in order-of-magnitude less time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Grid Integration of Small-Scale Photovoltaic Systems in Secondary Distribution Network—A Review

The relative share of renewable energy, specifically the solar photovoltaic (PV), is increasing exponentially in the world electric energy sector. This is a cumulative result of reduction in the cost of solar panels, improvement in the panel efficiency, and advancement in the associated power electronics. Among different types of PV plants, installation of small-scale rooftop PV is growing rapidly due to direct end-user benefits and lucrative governmental incentives. There are various standards developed in regards to grid integration of PVs and other distributed generations (DGs). Different power converter topologies are developed to interface the PV panel with the utility grid. To keep up with the stringent regulations imposed by the standards, various control strategies and grid synchronization methods have been developed. This review article amalgamates and summarizes all of the aforementioned aspects of a grid-integrated PV system including various standards, power stage architectures, grid synchronization methods, operation under extreme events, and control methodologies, pertaining to small-scale PV plants. This article will help freshman researchers to gain some familiarity with the topic and introduce them to some of the key issues encountered in this field.

14 SOLAR ENERGY↗

Advancing quantum simulations of the nuclear shell model with Gray-code–based resource-efficient protocols

Background: Some of the computational limitations in solving the nuclear many-body problem could be overcome by utilizing quantum computers. The nuclear shell-model calculations providing deeper insights into the properties of atomic nuclei are one such case with high demand for resources, as the size of the Hilbert space grows exponentially with the number of particles involved. Quantum algorithms are being developed to overcome these challenges and advance such calculations. Purpose: To develop quantum circuits for the nuclear shell-model, leveraging the capabilities of noisy intermediate-scale quantum (NISQ) devices. Here, we aim to minimize resource requirements (specifically in terms of qubits and gates) and strive to reduce the impact of noise by employing relevant mitigation techniques. Methods: We achieve noise resilience by designing an optimized Ansatz for the variational quantum eigensolver (VQE) based on Givens rotations and incorporating qubit-ADAPT-VQE in combination with variational quantum deflation (VQD) to compute ground and excited states, incorporating the zero-noise extrapolation mitigation technique. Furthermore, the qubit requirements are significantly reduced by mapping the basis states to qubits using Gray-code encoding and generalizing transformations of fermionic operators to efficiently represent many-body states. Results: By employing the resource-efficient protocols, we achieve the ground and excited state energy levels of 38 Ar and 6 Li with better accuracy. These energy levels are presented for noiseless simulations, noisy conditions, and after applying noise mitigation techniques. Results are compared for Jordan-Wigner and Gray-code encoding using VQE, qubit-ADAPT-VQE, and VQD. Conclusions: Our work highlights the potential of resource-efficient protocols to leverage the full potential of NISQ devices in scaling the nuclear shell model calculations, offering a pathway toward more complex quantum simulations in nuclear physics. This approach establishes a framework for studying other nuclear systems with improved quantum resource efficiency, marking a significant advancement in applying quantum computing to realistic nuclear physics applications.

Physics - Nuclear physics and radiation physics↗

Importance of gas heating in capacitively coupled radiofrequency plasma-assisted synthesis of carbon nanomaterials

In pursuit of diamond nanoparticles, a capacitively-coupled radio frequency flow-through plasma reactor was operated with methane-argon gas mixtures. Signatures of the final product obtained microscopically and spectroscopically indicated that the product was an amorphous form of graphite. This result was consistent irrespective of combinations of the macroscopic reactor settings. To explain the observed synthesis output, measurements of C 2 and gas properties were carried out by laser-induced fluorescence and optical emission spectroscopy. Strikingly, the results indicated a strong gas temperature gradient of 100 K per mm from the center of the reactor to the wall. Based on additional plasma imaging, a model of hot constricted region (filamentation region) was then formulated. It illustrated that, while the hot constricted region was present, the bulk of the gas was not hot enough to facilitate diamond sp 3 formation: characterized by much lower reaction rates, when compared to sp 2 , sp 3 formation kinetics are expected to become exponentially slow. This result was further confirmed by experiments under identical conditions but with a H 2 /CH 4 mixture, where no output material was detected: if graphitic sp 2 formation was expected as the main output material from the methane feedstock, atomic hydrogen would then be expected to etch it away in situ, such that the net production of that sp 2 -hybridized solid material is nearly a zero. Finally, the crucial importance of gas heating was corroborated by replacing RF with microwave source whereby facile sp 3 production was attained with H 2 /CH 4 gas mixture.

36 MATERIALS SCIENCE↗

Dynamic Cooling on Contemporary Quantum Computers

We study the problem of dynamic cooling whereby a target qubit is cooled at the expense of heating up N − 1 further identical qubits by means of a global unitary operation. A standard back-of-the-envelope high-temperature estimate establishes that the target qubit temperature can be dynamically cooled by at most a factor of 1 / N . Here we provide the exact expression for the minimum temperature to which the target qubit can be cooled and reveal that there is a crossover from the high initial temperature regime, where the scaling is 1 / N , to a low initial temperature regime, where a much faster scaling of 1 / N occurs. This slow, 1 / N scaling, which was relevant for early high-temperature NMR quantum computers, is the reason dynamic cooling was dismissed as ineffectual around 20 years ago; the fact that current low-temperature quantum computers fall in the fast, 1 / N scaling regime, reinstates the appeal of dynamic cooling today. We further show that the associated work cost of cooling is exponentially more advantageous in the low-temperature regime. We discuss the implementation of dynamic cooling in terms of quantum circuits and examine the effects of hardware noise. We successfully demonstrate dynamic cooling in a three-qubit system on a real quantum processor. Since the circuit size grows quickly with N , scaling dynamic cooling to larger systems on noisy devices poses a challenge. We therefore propose a suboptimal cooling algorithm, whereby relinquishing a small amount of cooling capability results in a drastically reduced circuit complexity, greatly facilitating the implementation of dynamic cooling on near-future quantum computers. Published by the American Physical Society2024

Physics↗

Controlling Symmetry Breaking Charge Transfer in BODIPY Pairs

Symmetry breaking charge transfer (SBCT) is a process in which a pair of identical chromophores absorb a photon and use its energy to transfer an electron from one chromophore to the other, breaking the symmetry of the chromophore pair. This excited state phenomenon is observed in photosynthetic organisms where it enables efficient formation of separated charges that ultimately catalyze biosynthesis. SBCT has also been proposed as a means for developing photovoltaics and photocatalytic systems that operate with minimal energy loss. It is known that SBCT in both biological and artificial systems is in part made possible by the local environment in which it occurs, which can move to stabilize the asymmetric SBCT state. However, how environmental degrees of freedom act in concert with steric and structural constraints placed on a chromophore pair to dictate its ability to generate long-lived charge pairs via SBCT remain open topics of investigation. In this work we compare a broad series of dipyrrin dimers that are linked by distinct bridging groups to discern how the spatial separation and mutual orientation of linked chromophores and the structural flexibility of their linker each impact SBCT efficiency. Across this material set, we observe a general trend that SBCT is accelerated as the spatial separation between dimer chromophores decreases, consistent with the expectation that the electronic coupling between these units varies exponentially with their separation. However, one key observation is that the rate of charge recombination following SBCT was found to slow with decreasing interchromophore separation, rather than speed up. This stems from an enhancement of the dimers’ structural rigidity due to increasing steric repulsion as the length of their linker shrinks. This rigidity further inhibits charge recombination in systems where symmetry has already enforced zero HOMO-LUMO overlap. Additionally, for the forward transfer the active torsion is shown to increase LUMO-LUMO coupling, allowing for faster SBCT within bridging groups. Further, by understanding trends for how rates of SBCT and charge recombination depend on a dimer’s internal structure and their environment, we identify design guidelines for creating artificial systems for driving sustained light-induced charge separation. Such systems can find application in solar energy technologies and photocatalytic applications but can serve as a model for light-induced charge separation in biological systems.

14 SOLAR ENERGY↗

Performance Evaluations of Noisy Approximate Quantum Fourier Arithmetic

The Quantum Fourier Transform (QFT) grants competitive advantages, especially in resource usage and circuit approximation, for performing arithmetic operations on quantum computers, and offers a potential route towards a numerical quantum-computational paradigm. In this paper, we utilize efficient techniques to implement QFT-based integer addition and multiplications. These operations are fundamental to various quantum applications including Shor’s algorithm, weighted sum optimization problems in data processing and machine learning and quantum algorithms requiring inner products. We carry out performance evaluations of these implementations based on IBM’s superconducting qubit architecture using different compatible noise models. We isolate the sensitivity of the component quantum circuits on both one-/two-qubit gate error rates, and the number of the arithmetic operands’ superposed integer states. We analyze performance, and identify the most effective approximation depths for quantum add and quantum multiply within the given context. We observe significant dependency of the optimal approximation depth on the degree of machine noise and the number of superposed states in certain performance regimes. Finally, we elaborate on the algorithmic challenges - relevant to signed, unsigned, modular and non-modular versions - that could also be applied to current implementations of QFT-based subtraction, division, exponentiation, and their potential tensor extensions. Here, we analyze performance trends in our results and speculate on possible future development within this computational paradigm.

97 MATHEMATICS AND COMPUTING↗

Score-Based Physics-Informed Neural Networks for High-Dimensional Fokker–Planck Equations

The Fokker-Planck (FP) equation is a foundational partial differential equation (PDE) in stochastic processes involving Brownian motions. However, the curse of dimensionality (CoD) poses a formidable challenge when dealing with high-dimensional FP equations. Although Monte Carlo simulation and (vanilla) Physics-Informed Neural Networks (PINNs) have shown the potential to tackle CoD, both methods exhibit significant numerical errors in high dimensions when dealing with the probability density function (PDF) associated with Brownian motion. The point-wise PDF values tend to decrease exponentially as dimensionality increases, surpassing the precision of numerical simulations and resulting in substantial errors. In addition, due to its massive sampling, Monte Carlo fails to offer fast sampling. Modeling the logarithm likelihood (LL) via vanilla PINNs transforms the FP equation into a notoriously difficult Hamilton-Jacobi-Bellman (HJB) equation, which is impractical for PINN learning, whose error grows rapidly with dimension. To this end, we propose a novel approach utilizing a score-based solver to fit the score function in stochastic differential equations (SDEs). The score function, defined as the gradient of the LL, plays a fundamental role in inferring LL and PDF and enables fast SDE sampling, offering an effective means to overcome the CoD. Three fitting methods, Score Matching (SM), Sliced Score Matching (SSM), and Score-PINN, are introduced, each contributing unique advantages in computational complexity, accuracy, and generality. The proposed score-based SDE solver operates in two stages: first, employing score matching or Score-PINN to acquire the score function; and second, solving the LL via an ordinary differential equation (ODE) using the obtained score function. Comparative evaluations across these methods showcase varying trade-offs. The proposed methodology is evaluated across diverse SDEs, including anisotropic Ornstein-Uhlenbeck processes, geometric Brownian motion, and Brownian motion with varying eigenspace. We also test various distributions, including Gaussian, Log-normal, Laplace, and Cauchy distributions. The numerical results demonstrate the score-based SDE solver’s stability, speed, and performance across different experimental settings, solidifying its potential as a solution to CoD for high-dimensional FP equations.

97 MATHEMATICS AND COMPUTING↗

Quantifying fission gas adsorption onto natural clinoptilolite in the presence of environmental air and water

Adsorption of noble gas fission products onto naturally occurring minerals is of interest for its potential to retain or retard emissions from nuclear fuel reprocessing operations or underground nuclear explosions. However, experimental studies of trace noble gas adsorption in the presence of air and water have largely focused on synthetic materials, such as activated carbon or metal-organic frameworks. Here, in this study, adsorption of Kr and Xe onto the naturally occurring zeolitic mineral clinoptilolite is studied in the presence of nitrogen and water. By varying the composition of the gas phase and monitoring the change in the combined adsorbate mass, the adsorbed concentration of noble gas is calculated gravimetrically. For dry clinoptilolite, the concentration of adsorbed Kr and Xe is linearly correlated with noble gas pressure and Henry's Law appears satisfactory, despite the presence of nitrogen at atmospheric pressures. However, the presence of water significantly reduces the adsorbed concentration of both Kr and Xe, which is typical in nanoporous sorbents. Here, an empirical bivariate model is presented, combining the Henry's Law adsorption model for a dry adsorbent with the exponential reduction in the presence of water, as reported by Lungu and Underhill in 1999. This model provides a means to estimate the adsorbate concentration at the trace partial pressures and higher water contents relevant to field-scale modeling of fission gas transport through the vadose zone.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

High Yield of B-Side Electron Transfer at 77 K in the Photosynthetic Reaction Center Protein from Rhodobacter sphaeroides

In this work, the primary electron transfer (ET) processes at 295 K and 77 K are compared for the Rhodobacter sphaeroides reaction center (RC) pigment-protein complex from thirteen mutants including a wild-type control. The engineered RCs bear mutations in the L and M polypeptides that largely inhibit ET from the excited state P* of the primary electron donor (P, a bacteriochlorophyll dimer) to the normally photoactive A-side cofactors and enhance ET to the C 2 -symmetry related, and normally photoinactive, B-side cofactors. P* decay is multi-exponential at both temperatures and modeled as arising from subpopulations that differ in contributions of two-step ET (e.g., P* → P + B B – → P + H B – ), one-step superexchange ET (e.g., P* → P + H B – ) and P* → ground state. [H B and B B are monomeric bacteriopheophytin and bacteriochlorophyll, respectively.] The relative abundances of the subpopulations and the inherent rate constants of the P* decay routes vary with temperature. Regardless, ET to produce P + H B – is generally faster at 77 K than at 295 K by about a factor of two. A key finding is that the yield of P + H B – , which ranges from ~5% to ~90% among the mutant RCs, is essentially the same at 77 K as at 295 K in each case. Overall, the results show that ET from P* to the B-side cofactors in these mutants does not require thermal activation and involves combinations of ET mechanisms analogous to those operative on the A side in the native RC.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Performance evaluations of signed and unsigned noisy approximate quantum Fourier arithmetic

The Quantum Fourier Transform (QFT) grants competitive advantages, especially in resource usage and circuit approximation, for performing arithmetic operations on quantum computers, and offers a potential route toward a numerical quantum-computational paradigm. In this paper, we utilize efficient techniques to implement QFT-based integer addition and multiplications. These operations are fundamental to various quantum applications including Shor’s algorithm, weighted-sum optimization problems in data processing and machine learning, and quantum algorithms requiring inner products. We carry out performance evaluations of these implementations based on IBM’s superconducting-qubit architecture using different compatible noise models. We isolate the sensitivity of the component quantum circuits on both one-/two-qubit gate error rates, and the number of the arithmetic operands’ superposed integer states. We analyze performance and identify the most effective approximation depths for unsigned quantum addition and quantum multiplication within the given context. We then perform a similar analysis of signed addition and compare to the unsigned results. We observe significant dependency of the optimal approximation depth on the degree of machine noise and the number of superposed states in certain performance regimes. Finally, we elaborate on the algorithmic challenges—relevant to signed, unsigned, modular and non-modular versions—that could also be applied to current implementations of QFT-based subtraction, division, exponentiation, and their potential tensor extensions. Here, we analyze the performance trends in our results and speculate on possible future developments within this computational paradigm.

Computational models↗

Wettability and CHF limits of Accident-Tolerant nuclear fuel cladding materials in light water reactor conditions

Here, we present the results of experimental investigations aimed at evaluating the thermal–hydraulic performance of chromium-coated zircaloy, i.e., one of the most promising accident tolerant fuel (ATF) cladding material for light water nuclear reactors. Precisely, we investigate the wettability and critical heat flux (CHF) limits of chromium-coated and conventional zircaloy surfaces in prototypical reactor conditions. For both surface types, we measure the contact angle in a vapor-saturated atmosphere from ambient pressure to the operating pressure of pressurized water reactors (PWRs), i.e., ~15 MPa. We measure the ambient-pressure steady-state flow boiling CHF with a spatially uniform heat flux. We measure the high-pressure steady-state flow boiling CHF with a cosine shape heat flux (up to 20 MPa) and with a uniform heat flux (up to 15 MPa), also exploring the effect of CRUD deposits on the chromium-coated surface. Our results reveal that the chromium surface and the bare zircaloy surface have similar wettability and both become super-hydrophilic in PWR conditions, and that there is practically no difference in the steady-state CHF limits, both at low-pressure and high-pressure conditions, also when the chromium-coated surface is covered by a CRUD deposit. However, while the chromium-coating does not improve the CHF compared to the bare zircaloy surface, it improves the post-CHF behavior. The chromium coating prevents the reaction between zircaloy and steam, which results in the formation of a brittle zirconium oxide through the surface of the cladding. We also measure the transient CHF under exponentially escalating heat flux inputs of a nano-smooth and a rough surface mimicking a chromium-coated zircaloy cladding. Interestingly, the results of the transient heat flux tests suggest that the CHF limit for very short periods (i.e., fast transients) is independent of the surface finish, being the same for a rough chromium surface or a nano-smooth surface.

36 MATERIALS SCIENCE↗

Biased degenerate ground-state sampling of small Ising models with converged quantum approximate optimization algorithm

The quantum alternating operator ansatz, a generalization of the quantum approximate optimization algorithm (QAOA), is a quantum algorithm used for approximately solving combinatorial optimization problems. QAOA typically uses the transverse field mixer as the driving Hamiltonian. One of the interesting properties of the transverse field driving Hamiltonian is that it results in nonuniform sampling of degenerate ground states of optimization problems. In this study, we numerically examine the fair sampling properties of the transverse field mixer QAOA, and Grover mixer QAOA (GM-QAOA), which provides theoretical guarantees of fair sampling of degenerate optimal solutions, up to a large enough p such that the mean expectation value converges to an optimal approximation ratio of 1. This comparison is performed with high-quality heuristically computed, but not necessarily optimal, QAOA angles, which give strictly monotonically improving solution quality as p increases. These angles are computed using the Julia based numerical simulation software JuliQAOA. Fair sampling of degenerate ground states is quantified using the Shannon entropy of the ground-state amplitudes distribution. The fair sampling properties are reported on several quantum signature Hamiltonians from previous quantum annealing fair sampling studies. Small random fully connected spin glasses are shown, which exhibit exponential suppression of some degenerate ground states with transverse field mixer QAOA. The transverse field mixer QAOA simulations show that some problem instances clearly saturate the Shannon entropy of 0 with a maximally biased distribution that occurs when the learning converges to an approximation ratio of 1 while other problem instances never deviate from a maximum Shannon entropy (uniform distribution) at any p step. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Classically verifiable quantum advantage from a computational Bell test

Existing experimental demonstrations of quantum computational advantage have had the limitation that verifying the correctness of the quantum device requires exponentially costly classical computations. Here we propose and analyse an interactive protocol for demonstrating quantum computational advantage, which is efficiently classically verifiable. Our protocol relies on a class of cryptographic tools called trapdoor claw-free functions. Although this type of function has been applied to quantum advantage protocols before, our protocol employs a surprising connection to Bell’s inequality to avoid the need for a demanding cryptographic property called the adaptive hardcore bit, while maintaining essentially no increase in the quantum circuit complexity and no extra assumptions. Leveraging the relaxed cryptographic requirements of the protocol, we present two trapdoor claw-free function constructions, based on Rabin’s function and the Diffie–Hellman problem, which have not been used in this context before. We also present two independent innovations that improve the efficiency of our implementation and can be applied to other quantum cryptographic protocols. First, we give a scheme to discard so-called garbage bits, removing the need for reversibility in the quantum circuits. Second, we show a natural way of performing postselection that reduces the fidelity needed to demonstrate quantum advantage. Combining these results, we describe a blueprint for implementing our protocol on Rydberg atom-based quantum devices, using hardware-native operations that have already been demonstrated experimentally.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Simulating scalar field theories on quantum computers with limited resources

Here we present a quantum algorithm for implementing φ 4 lattice scalar field theory on qubit computers. The field is represented in the discretized field amplitude basis. The number of qubits and elementary gates required by the implementation of the evolution operator is proportional to the lattice size. The algorithm allows efficient φ 4 state preparation for a large range of input parameters in both the normal and broken-symmetry phases. The states are prepared using a combination of variational and adiabatic evolution methods. First, the ground state of a local Hamiltonian, which includes the φ 4 self-interaction, is prepared using short variational circuits. Next, this state is evolved by switching on the coupling between the lattice sites adiabatically. The parameters defining the local Hamiltonian are adjustable and constitute the input of our algorithm. We present a method to optimize these parameters in order to reduce the adiabatic time required for state preparation. For preparing broken-symmetry states, the adiabatic evolution problems caused by crossing the phase transition critical line and by the degeneracy of the broken-symmetry ground state can be addressed using an auxiliary external field which gradually turns off during the adiabatic process. We show that the time dependence of the external field during the adiabatic evolution is important for addressing the broken-symmetry ground state degeneracy. The adiabatic time dependence on the inverse error tolerance can be reduced from quadratic to linear by using a field strength that decreases exponentially in time relative to one that decreases linearly.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Inferring the scrape-off layer heat flux width in a divertor with a low degree of axisymmetry

Plasma facing components (PFCs) in the next generation of tokamak devices will operate in challenging environments, with heat loads predicted to exceed 10 MWm -2 . The magnitude of these heat loads is set by the width of the channel, the ‘scrape-off layer’ (SOL), into which heat is exhausted, and can be characterised by an e-folding length scale for the decay of heat flux across the channel. It is expected this channel will narrow as tokamaks move towards reactor relevant conditions. Understanding the processes involved in setting the SOL heat flux width is imperative to be able to predict the heat loads PFCs must handle in future devices. Measurements of the SOL width are performed on the high-field spherical tokamak, ST40, using a newly commissioned infrared thermography system. With its high on-axis toroidal magnetic field (≥1.5 T) ST40 is uniquely positioned to investigate the influence of toroidal field on the heat flux width in spherical tokamaks, whilst also extending measurements of the SOL width in spherical tokamaks to increased poloidal field (≥0.3 T). Due to the divertor on ST40 having a low degree of axisymmetry, it is necessary for a set of radial measurements of the heat flux to be taken across the divertor, made possible using an automated toolchain that fully incorporates its 3D geometry. These radial profiles are combined with the magnetic geometry of the plasma to infer the width of the SOL, with both Eich and double exponential profiles of heat flux observed. A reduction in the heat flux is observed toroidally across part of the divertor, along with increased heat loads observed locally around the edges of the tiles. Future work in characterising the impact of tile misalignment and uncertainties in the reconstructed divertor magnetic geometry is required in order to further understand the observed heat flux patterns, as are additional investigations into the role potentially being played by an inhomogeneous sheath electric field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗