Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Real numbers”

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 559 records · Page 31

Accelerating network layouts using graph neural networks

Graph layout algorithms used in network visualization represent the first and the most widely used tool to unveil the inner structure and the behavior of complex networks. Current network visualization software relies on the force-directed layout (FDL) algorithm, whose high computational complexity makes the visualization of large real networks computationally prohibitive and traps large graphs into high energy configurations, resulting in hard-to-interpret “hairball” layouts. Here we use Graph Neural Networks (GNN) to accelerate FDL, showing that deep learning can address both limitations of FDL: it offers a 10 to 100 fold improvement in speed while also yielding layouts which are more informative. We analytically derive the speedup offered by GNN, relating it to the number of outliers in the eigenspectrum of the adjacency matrix, predicting that GNNs are particularly effective for networks with communities and local regularities. Finally, we use GNN to generate a three-dimensional layout of the Internet, and introduce additional measures to assess the layout quality and its interpretability, exploring the algorithm’s ability to separate communities and the link-length distribution. The novel use of deep neural networks can help accelerate other network-based optimization problems as well, with applications from reaction-diffusion systems to epidemics.

97 MATHEMATICS AND COMPUTING↗

EV Hosting Capacity Analysis on Distribution Grids

The increasing trend in electric vehicle (EV) adoption can cause challenges to traditional electric grid operations if utilities are not equipped with tools and methods to effectively manage these fleets. Growing EV charging loads will alter the magnitude and duration of conventional peaks in demand profiles and even significantly shift them, potentially causing operational violations in the distribution grid. This paper presents the development and results of an EV hosting capacity tool to quantify the impacts of injecting large numbers of EV charging loads and to determine the available capacity of existing distribution feeders to continue providing reliable and affordable grid operations. Tools like the hosting capacity analysis would enable utilities to better prepare for grid operations in the near future while exploring the impact and effectiveness of strategies to manage these loads, such as peak pricing and smart charging. This paper evaluates the hosting capacity of some real-world feeders to accommodate EV charging loads, including extreme fast-charging options.

distribution grid↗

CONDENSATION-EVAPORATION MODELING OF DRIFT FROM A MECHANICAL COOLING TOWER

Numerical studyfor investigatingthe phase change behavior of a drift from a mechanical draft cooling tower is presented. A k-ε turbulence model, discrete phase model, species transport equations, and a heat-mass exchange model are solved simultaneouslyto capture the heat and mass transfer effects of a drift condensing and evaporatingeffects. The simulation result is compared to an available real scale test data and its limitation is also discussed. A parametric study is carried out by varying the ambient dry bulb temperature, ambient relative humidity, wind speed, liquid droplet diameter, and the drift amount. It is found that by combining the evaporation and condensation models, a drift phase change phenomena is accurately depicted. A parametric study reveals that wind speed impact is affected by the ambient temperature and the number of liquid droplets and their size play a critical role in the phase change phen

Han, Kai↗

Navigating collinear superspace

We introduce a new set of effective field theory rules for constructing Lagrangians with \( \mathcal{N} \) = 1 supersymmetry in collinear superspace. In the standard superspace treatment, superfields are functions of the coordinates \( \left({x}^{\mu },{\theta}^{\alpha },{\theta}^{\dagger \overset{\cdot }{\alpha }}\right) \) , and supersymmetry preservation is manifest at the Lagrangian level in part due to the inclusion of auxiliary F - and D -term components. By contrast, collinear superspace depends on a smaller set of coordinates ( x μ , η , η † ), where η is a complex Grassmann number without a spinor index. This provides a formulation of supersymmetric theories that depends exclusively on propagating degrees of freedom, at the expense of obscuring Lorentz invariance and introducing inverse momentum scales. After establishing the general framework, we construct collinear superspace Lagrangians for free chiral matter and non-Abelian gauge fields. For the latter construction, an important ingredient is a superfield representation that is simultaneously chiral, anti-chiral, and real; this novel object encodes residual gauge transformations on the light cone. Additionally, we discuss a fundamental obstruction to constructing inter- acting theories with chiral matter; overcoming these issues is the subject of our companion paper, where we introduce a larger set of superfields to realize the full range of interactions compatible with \( \mathcal{N} \) = 1. Along the way, we provide a novel framing of reparametrization invariance using a spinor decomposition, which provides insight into this important light-cone symmetry.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

High-Fidelity CFD Simulation of Mixed Convection and Forced Convection in a Pebble Bed Test Reactor Core

The Hermes low-power [35-MW(thermal)] reactor will be built and operated by Kairos Power LLC (KP) to demonstrate its fluoride salt-cooled high-temperature reactor (FHR) technology. In the KP FHR, the reactor core is composed of randomly packed pebbles with TRISO fuel particles inside with FLiBe flow upward through the core acting as a coolant. Previous numerical and experimental studies have been limited to either a small-size bed or to a lack of detailed measurements for heat transfer. Here, to address the lack of high-fidelity heat transfer data in a real-size FHR core, in this study, we simulated a pebble bed core with 34 374 pebbles randomly packed, similar to the Hermes reactor's size. The core radius was 14 times that of the pebble diameter, while the core height was 45 times. In this work, we were particularly interested in a mixed convection regime, where buoyancy is important. Therefore, we performed several large-eddy simulations at different Reynolds numbers (160 to 1000) with gravitational force included. The spectral element computational fluid dynamics code NekRS with graphics processing unit acceleration was used for this study. The low-Mach number approximation was applied to address property changes in the FLiBe and to account for buoyancy. A pure hexahedral mesh with 60 million elements was generated by the Voronoi cell method. At the polynomial order of 5, the total degrees of freedom was 7.5 billion. The developed case in this work is the first of its kind in terms of size and complexity. The local numerical data across the domain were obtained and compared with empirical correlations. After examining the data, we found the following conclusions. For pressure drop, the Reger correlation predicted less than a 5% error. On the other hand, for heat transfer, the Wakao correlation outperformed the others. Based on our findings, we recommend the use of the Wakao correlation for the Nusselt number calculation, and for pressure drop, the KTA (Kerntechnischer Ausschuss) correclation, among the available experimental correlations. In conclusion, the Reger direct numerical simulation-driven correlation for pressure drops should also be considered, given its best agreement with our calculations.

Mixed Convection↗

Numerical evidence against advantage with quantum fidelity kernels on classical data

Quantum machine learning techniques are commonly considered one of the most promising candidates for demonstrating practical quantum advantage. In particular, quantum kernel methods have been demonstrated to be able to learn certain classically intractable functions efficiently if the kernel is well aligned with the target function. In the more general case, quantum kernels are known to suffer from exponential “flattening” of the spectrum as the number of qubits grows, preventing generalization and necessitating the control of the inductive bias by hyperparameters. We show that the general-purpose hyperparameter-tuning techniques proposed to improve the generalization of quantum kernels lead to the kernel becoming well approximated by a classical kernel, removing the possibility of quantum advantage. We provide extensive numerical evidence for this phenomenon utilizing multiple previously studied quantum feature maps and both synthetic and real data. Our results show that unless novel techniques are developed to control the inductive bias of quantum kernels, they are unlikely to provide a quantum advantage on classical data that lacks special structure.

Quantum Information, Science & Technology↗

Inverse modeling of circular lattices via orbit response measurements in the presence of degeneracy

The number and location of beam position monitors (BPMs) and steerers with respect to the quadrupoles in a circular lattice can lead to degeneracy in the context of fitting linear optics and extracting lattice information from measured closed orbits. Furthermore, the measurement uncertainties due to the imperfection of BPMs and steerers can be propagated by the fitting process in ways that prohibit the successful extraction of discrepancies between lattice elements in the real machine and their description in the corresponding model. We systematically studied the influence of the placement of BPMs and steerers on the reconstruction of linear optics and corresponding lattice information. The derivative of orbit response coefficients with respect to the quadrupole strengths, the Jacobian, is derived as an analytical formula. This analytical version of the Jacobian is used to further derive the theoretical limitations of fitting linear optics from closed orbits in terms of the placement of BPMs and steerers. It is further demonstrated that when evaluating the Jacobian during the fitting procedure, the analytical version can be used in place of the conventional finite-difference computation. This allows for greatly improved efficiency when computing the Jacobian during each iteration of the fitting procedure. The approach is tested with large-scale simulations and the findings are verified by measurement data taken on SIS18 synchrotron at GSI Helmholtz Centre for Heavy Ion Research. The presented methods are of general nature and can be applied to other accelerator lattices as well. The fitting procedure by using the analytical Jacobian is tested in conjunction with various methods for mitigating quasidegeneracy and the results agree with those obtained by using the conventional Jacobian via finite-difference approximation.

47 OTHER INSTRUMENTATION↗

Paired Hydro-Battery Hybrid System Operations Using MIQP Based Multi-Objective Optimization

Hybridizing a hydropower plant with a battery energy storage system is often a very expensive decision. So to make the case for feasibility of hybridization the cost benefit analysis must be comprehensive. Most literature found on this subject only uses one out of many different available value streams to carry out a feasibility analysis. To that end, in this study we propose a multi-objective optimization formulation and a mixed-integer quadratic programming optimization engine that considers multiple different value streams and optimizes hydro-battery hybrid system paired operations to maximize revenue generated from energy arbitrage while simultaneously minimizing the total hydro turbine mileage thus reducing turbine starts and stops and improving turbines life all while following all environmental constraints and not violating any limitations either regulatory or preferential (i.e., to support recreational activities like white water rafting, etc.). In this study we also implemented the developed optimization methodology to a real-world case study of Bagnell dam hydropower facility (8 units totaling 240 MW rated power capacity) which is owned and operated by our industry partner Ameren Energy Inc. The case study outcomes shows that by hybridizing the Bagnell dam hydropower facility with a 60 MW x 2-hr battery energy storage system, the annual benefits can be increased over \$6 million while reducing the annual mileage averaged per turbine and number of start/stops by over 98\% and 85\% respectively.

Chalishazar, Vishvas H.↗

Accelerating Bilevel Optimization With Hierarchical Many-Threaded Parallel Differential Evolution

Bilevel optimization is encountered in many relevant real-world applications. The main feature of this type of problem is that an upper-level optimization problem is constrained by a nested lower-level optimization problem. Because of this nested structure, bilevel problems (BLPs) are usually computationally expensive to solve. Differential evolution (DE) has demonstrated promising results in solving BLPs of relatively small scales. As the problem scale increases, the decision space becomes intrinsically larger, requiring a growing number of function evaluations for the method to work properly. In this context, heavy parallelization and high-performance computing techniques are indispensable to enable the resolution of more complex and challenging optimization problems. Hence, we propose a hierarchical many-threaded parallel DE approach for BLPs, where both levels are parallelized. The computational experiments demonstrate that the parallel implementation achieved runtime speeds ranging from 44 to 2559 times faster than the sequential version on a well-known scalable SMD benchmark test problem when executed on an NVIDIA A100 GPU. The findings indicate that the algorithm’s convergence is strongly influenced by the number of both upper- and lower-level generations. Moreover, the success of experiments with large-scale problems is closely linked to the choice of small population sizes.

Dufek, Amanda S↗

Distributed Real-time Plume Monitoring for Deep Sea Mineral Extraction​

In the emerging industry of deep-sea mining for minerals and deposits (e.g. polymetallic nodules for nickel, cobalt, copper, and manganese), more data is required to understand the effects of sediment plume generation and predict the distribution of disturbed sediment. There are two main sources of plume generation, the first being at the active mining site where the “collector” directly removes the top layer of the sea floor. The other is the “midwater plume” consisting of unwanted sediment that was collected during extraction that is pumped back into the aphotic zone. The vast majority of plume generation is caused by the collector, causing detrimental and long-lasting impacts on seafloor ecosystems due to the lack of wave activity or strong currents at the sea floor. Therefore, it is crucial to invest in the infrastructure to support the study and constant monitoring over a large area of the sea floor where plume generation is present. Due to the limited number of usable channels and power requirements, current subsea wireless communications technologies are not well suited to instrumenting the large areas of the sea floor needed to monitor plume migration. The scope of this effort is to transition experimental demonstrations of high-bandwidth, full-duplex scalable underwater laser communications to the seafloor in an open ocean environment. Specifically tackling challenges associated with the dynamic nature of the subsea world, including but not limited to, deployment logistics, sustainability, and range. The goal is to enable the internet of underwater things for deep sea industries by broadening the capabilities of subsea communications. By using high-precision laser transmitters, many of the challenges current subsea optical systems face can be circumvented, such as power consumption, interference, and bandwidth limitations. This approach lends itself to wireless interlinking multi-node networks, in series or parallel, facilitating the implementation of a wide array of sensor types. This interlinking allows all the data gathered from the network to be processed through a single hardline uplink to the surface, lowering the complexity required for near real-time data processing. Additionally, the laser control systems produce metadata that can be used to help characterize the water column between the nodes. Combining data from various sensors such as turbidity, temperature, current velocity with metadata such as beam attenuation and deflection can produce a high-resolution model of sea floor conditions around an active mining zone. The resulting near real-time model can be used to optimize location and flow rate of the mining operation to minimize and quantify the environmental impact.

Mons, Ishan↗

Author Correction: Proton irradiation-decelerated intergranular corrosion of Ni-Cr alloys in molten salt

In the original version of this article, numbers for the beam current densities were incorrectly given as 1.5, 2.0, and 2.5 mA cm -2 in various locations, instead of the correct values 0.3, 0.4, and 0.5 mA cm -2 . This was owing to a measurement error coming from the indirect correspondence between the Faraday cup and the beam profile monitor (BPM) on the accelerator as the proton beam traveling along the beamline resulted in mismatching between calculated values and real beam currents, which were obtained by ex-situ calibration. The following changes have been made in the correct version. The eighth sentence of the Results, and the figure legend of Fig. 2j state '0.3, 0.4, 0.5 mA cm -2 ' in place of '1.5, 2.0, 2.5 mA cm -2 '. The figure legend of Fig. 1j–l states '0.5, 0.4, 0.3 mA cm -2 ' in place of '2.5, 2.0, 1.5 mA cm -2 '. Figure 1e states '0.5 mA cm -2 ' in place of '2.5 mA cm -2 ', Fig. 1h states '0.4 mA cm -2 ' in place of '2.0 mA cm -2 ', and Fig. 1k states ‘0.3 mA cm -2 ' in place of ‘1.5 mA cm -2 '. Figure 2c states '0.5 mA cm -2 ' in place of '2.5 mA cm -2 ', Fig. 2d states '0.4 mA cm -2 ' in place of '2.0 mA cm -2 ', and Fig. 2e states '0.3 mA cm -2 ' in place of '1.5 mA cm -2 '. Figure 2j and k states ‘0.3 mA cm -2 ', '0.4 mA cm -2 ', '0.5 mA cm -2 ' in place of '1.5 mA cm -2 ', '2.0 mA cm -2 ', and '2.5 mA cm -2 ', respectively. The figure legend of Fig. 3b states '0.4 mA cm -2 ' in place of '2.0 mA cm -2 '. The figure legend of Supplementary Fig. 1 states '0.5 mA cm -2 ' in place of '2.5 mA cm -2 '. This has been corrected in the PDF and HTML versions of the Article.

36 MATERIALS SCIENCE↗

Extensive analysis of reconstruction algorithms for DESI 2024 baryon acoustic oscillations

Reconstruction of the baryon acoustic oscillation (BAO) signal has been a standard procedure in BAO analyses over the past decade and has helped to improve the BAO parameter precision by a factor of ∼2 on average. The Dark Energy Spectroscopic Instrument (DESI) BAO analysis for the first year (DR1) data uses the “standard” reconstruction framework, in which the displacement field is estimated from the observed density field by solving the linearized continuity equation in redshift space, and galaxy and random positions are shifted in order to partially remove non-linearities. There are several approaches to solving for the displacement field in real survey data, including the multigrid (MG), iterative Fast Fourier Transform (iFFT), and iterative Fast Fourier Transform particle (iFFTP) algorithms. In this work, we analyze these algorithms and compare them with various metrics including two-point statistics and the displacement itself using realistic DESI mocks. We focus on three representative DESI samples, the emission line galaxies (ELG), quasars (QSO), and the bright galaxy sample (BGS), which cover the extreme redshifts and number densities, and potential wide-angle effects. We conclude that the MG and iFFT algorithms agree within 0.4% in post-reconstruction power spectrum on BAO scales with the RecSym convention, which does not remove large-scale redshift space distortions (RSDs), in all three tracers. The RecSym convention appears to be less sensitive to displacement errors than the RecIso convention, which attempts to remove large-scale RSDs. However, iFFTP deviates from the first two; thus, we recommend against using iFFTP without further development. In addition, we provide the optimal settings for reconstruction for five years of DESI observation. The analyses presented in this work pave the way for DESI DR1 analysis as well as future BAO analyses.

79 ASTRONOMY AND ASTROPHYSICS↗

Opinion: A critical evaluation of the evidence for aerosol invigoration of deep convection

Deep convective updraft invigoration via indirect effects of increased aerosol number concentration on cloud microphysics is frequently cited as a driver of correlations between aerosol and deep convection properties. Here, we critically evaluate the theoretical, modeling, and observational evidence for warm- and cold-phase invigoration pathways. Though warm-phase invigoration is plausible and theoretically supported via lowering of the supersaturation with increased cloud droplet concentration in polluted conditions, the significance of this effect depends on substantial supersaturation changes in real-world convective clouds that have not been observed. Much of the theoretical support for cold-phase invigoration depends on unrealistic assumptions of instantaneous freezing and unloading of condensate in growing, isolated updrafts. When applying more realistic assumptions, impacts on buoyancy from enhanced latent heating via fusion in polluted conditions are largely canceled by greater condensate loading. Many foundational observational studies supporting invigoration have several fundamental methodological flaws that render their findings incorrect or highly questionable. Thus, much of the evidence for invigoration has come from numerical modeling, but different models and setups have produced a vast range of results. Furthermore, modeled aerosol impacts on deep convection are rarely tested for robustness, and microphysical biases relative to observations persist, rendering many results unreliable for application to the real world. Without clear theoretical, modeling, or observational support, and given that enervation rather than invigoration may occur for some deep convective regimes and environments, it is entirely possible that the overall impact of cold-phase invigoration is negligible. Substantial mesoscale variability of dominant thermodynamic controls on convective updraft strength coupled with substantial updraft and aerosol variability in any given event are poorly quantified by observations and present further challenges to isolating aerosol effects. Observational isolation and quantification of convective invigoration by aerosols is also complicated by limitations of available cloud condensation nuclei and updraft speed proxies, aerosol correlations with meteorological conditions, and cloud impacts on aerosols. Furthermore, many cloud processes, such as entrainment and condensate fallout, modulate updraft strength and aerosol–cloud interactions, varying with cloud life cycle and organization, but these processes remain poorly characterized. Considering these challenges, recommendations for future observational and modeling research related to aerosol invigoration of deep convection are provided.

54 ENVIRONMENTAL SCIENCES↗

Emulation of Synaptic Plasticity in WO 3 ‐Based Ion‐Gated Transistors

Neuromorphic systems, inspired by the human brain, promise significant advancements in computational efficiency and power consumption by integrating processing and memory functions, thereby addressing the von Neumann bottleneck. This paper explores the synaptic plasticity of a WO3-based ion-gated transistor (IGT) in [EMIM][TFSI] and a 0.1 mol L −1 LiTFSI in [EMIM][TFSI] for neuromorphic computing applications. Cyclic voltammetry (CV), transistor characteristics, and atomic force microscopy (AFM) force–distance (FD) profiling analyses reveal that Li + brings about ion intercalation, together with higher mobility and conductance, and slower response time (τ). WO 3 IGTs exhibit spike amplitude-dependent plasticity (SADP), spike number-dependent plasticity (SNDP), spike duration-dependent plasticity (SDDP), frequency-dependent plasticity (FDP), and paired-pulse facilitation (PPF), which are all crucial for mimicking biological synaptic functions and understanding how to achieve different types of plasticity in the same IGT. The findings underscore the importance of selecting the appropriate ionic medium to optimize the performance of synaptic transistors, enabling the development of neuromorphic systems capable of adaptive learning and real-time processing, which are essential for applications in artificial intelligence (AI).

36 MATERIALS SCIENCE↗

High-throughput computation of electric polarization in solids via Berry flux diagonalization

Electric polarization in the absence of an externally applied electric field is a key property of polar materials, but the standard interpolation-based ab initio approach to compute polarization differences within the modern theory of polarization presents challenges for automated high-throughput calculations. Berry flux diagonalization [J. Bonini et al., Phys. Rev. B 102, 045141 (2020)] has been proposed as an efficient and reliable alternative, though it has yet to be widely deployed. Here, we assess Berry flux diagonalization using ab initio calculations of a large set of materials, introducing and validating heuristics that ensure branch alignment with a minimal number of intermediate interpolated structures. Our automated implementation of Berry flux diagonalization succeeds in cases where prior interpolation-based workflows fail due to band-gap closures or branch ambiguities. Benchmarking with ab initio calculations of 176 candidate ferroelectrics, we demonstrate the efficacy of the approach on a broad range of insulating materials and obtain accurate effective polarization values with fewer interpolated structures than prior automated interpolation-based workflows. Our real-space heuristics that can predict gauge stability a priori from ionic displacements enable a general automated framework for reliable polarization calculations and efficient high-throughput screening of chemically and structurally diverse polar insulators. These results establish Berry flux diagonalization as a robust and efficient method to compute the effective polarization of solids and to accelerate the data-driven discovery of functional polar materials.

Poteshman, Abigail N. [University of Chicago, IL (↗

Measurement of 2D density profiles using a second-harmonic, dispersion interferometer

A second-harmonic, dispersion interferometer is used to image large-area (≃5 cm 2 ) plasma-jet and gas-jet density profiles. Achromatic telescopes magnify the diameters of the primary-laser beam (1064 nm) and its second-harmonic (532 nm) before probing the sample and de-magnify the beam diameters after the sample, where the primary beam transfers its phase change to a second, second-harmonic beam, allowing the sample’s dispersive-phase change to be measured between two, orthogonally polarized second harmonic beams. The telescopes produce an azimuthally symmetric, dispersive-phase shift in the sample + background phase-change image and in the background phase-change image, which is removed by digital subtraction. The interferometer’s performance was verified using standard-optical components as dispersive elements (BK7 lenses and wedge plates), resolving a minimum, phase-change sensitivity of ΔΦ min ≳ 15 mrad and spatial resolution of Δx res ≃ 100 μm. The phase change produced by unknown-density objects (a pulsed-plasma-jet and a pulsed-gas-jet) was measured, and their data were used to recover the original, 2D density profiles using an inverse Abel transform: peak-number density, N gas ≃ 6 × 10 20 cm –3 and N e ≃ 5 × 10 16 cm –3 ; line-integrated density, ∫N gas dl ~ 2 × 10 19 cm –2 and ∫N e dl ~ 1 × 10 16 cm –2 . The techniques and methods developed here are scalable to even larger probe-beam diameters and frame-capture rates, leading to a diagnostic capability that is well-suited for applications involving the real-time measurement of density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Trapped-ion quantum simulation of collective neutrino oscillations

It is well known that the neutrino flavor in extreme astrophysical environments changes under the effect of three contributions: the vacuum oscillation, the interaction with the surrounding matter, and the collective oscillations due to interactions between different neutrinos. The latter adds a nonlinear contribution to the equations of motion, making the description of their dynamics complex. In this work we study various strategies to simulate the coherent collective oscillations of a system of N neutrinos in the two-flavor approximation using quantum computation. This was achieved by using a pair-neutrino decomposition designed to account for the fact that the flavor Hamiltonian, in the presence of the neutrino-neutrino term, presents an all-to-all interaction that makes the implementation of the evolution dependent on the qubit topology. We analyze the Trotter error caused by the decomposition demonstrating that the complexity of the implementation of time evolution scales polynomially with the number of neutrinos and that the noise from near-term quantum device simulation can be reduced by optimizing the quantum circuit decomposition and exploiting a full-qubit connectivity. We find that the gate complexity using second order Trotter-Suzuki formulas scales better with system size than with other decomposition methods such as quantum signal processing. In conclusion, we finally present the application and the results of our algorithm on a real quantum device based on trapped-ion qubits.

79 ASTRONOMY AND ASTROPHYSICS↗

Gate-based quantum computing for protein design

Protein design is a technique to engineer proteins by permuting amino acids in the sequence to obtain novel functionalities. However, exploring all possible combinations of amino acids is generally impossible due to the exponential growth of possibilities with the number of designable sites. The present work introduces circuits implementing a pure quantum approach, Grover’s algorithm, to solve protein design problems. Our algorithms can adjust to implement any custom pair-wise energy tables and protein structure models. Moreover, the algorithm’s oracle is designed to consist of only adder functions. Quantum computer simulators validate the practicality of our circuits, containing up to 234 qubits. However, a smaller circuit is implemented on real quantum devices. Our results show that using iterations, the circuits find the correct results among all N possibilities, providing the expected quadratic speed up of Grover’s algorithm over classical methods (i.e.,).

59 BASIC BIOLOGICAL SCIENCES↗