Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “POLYNOMIAL”

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 865 records · Page 48

Jordan–Wigner Transformation for the Description of Strong Correlation in Fermionic Systems

Seniority is a useful way of organizing Hilbert space for strongly correlated systems. The exact zero-seniority wave function, doubly occupied configuration interaction (DOCI), provides accurate results (given the right orbitals) for many strongly correlated electronic systems but has a combinatorial computational cost. In many cases, pair coupled cluster doubles provide a polynomial-cost approximation that closely reproduces the energies of DOCI, but it breaks down in some cases and, as shown herein, it does not provide particularly good density matrices. In this article, we demonstrate that by using the Jordan–Wigner transformation to turn the seniority zero problem back into a Fermionic one, we can provide mean-field variational results of DOCI quality for the Hubbard model and a few small molecular dissociation examples, with polynomial cost, both for the energies and for density matrices, all while being protected from collapse. This success is rooted in the proof we provide, showing that the Hartree–Fock wave function on the Jordan–Wigner-transformed Hamiltonian transforms back to variational coupled cluster doubles in the seniority zero representation, but restricted to have determinant rather than permanent amplitude coefficients, without compromising its overall accuracy.

74 ATOMIC AND MOLECULAR PHYSICS↗

State Preparation of Antisymmetrized Geminal Power on a Quantum Computer without Number Projection

The antisymmetrized geminal power (AGP) is equivalent to the number projected Bardeen–Cooper–Schrieffer (PBCS) wave function. It is also an elementary symmetric polynomial (ESP) state. We generalize previous research on deterministically implementing the Dicke state to a state preparation algorithm for an ESP state, or equivalently AGP, on a quantum computer. Our method is deterministic and has polynomial cost, and it does not rely on number symmetry breaking and restoration. We also show that our circuit is equivalent to a disentangled unitary paired coupled cluster operator and a layer of unitary Jastrow operator acting on a single Slater determinant. Here, the method presented herein highlights the ability of disentangled unitary coupled cluster to capture nontrivial entanglement properties that are hardly accessible with traditional Hartree–Fock based electronic structure methods.

97 MATHEMATICS AND COMPUTING↗

Energy Minimization of Discrete Protein Titration State Models Using Graph Theory

There are several applications in computational biophysics which require the optimization of discrete interacting states; e.g., amino acid titration states, ligand oxidation states, or discrete rotamer angles. Such optimization can be very time-consuming as it scales exponentially in the number of sites to be optimized. In this paper, we describe a new polynomial-time algorithm for optimization of discrete states in macromolecular systems. This algorithm was adapted from image processing and uses techniques from discrete mathematics and graph theory to restate the optimization problem in terms of maximum flow-minimum cut graph analysis. The interaction energy graph, a graph in which vertices (amino acids) and edges (interactions) are weighted with their respective energies, is transformed into a flow network in which the value of the minimum cut in the network equals the minimum free energy of the protein, and the cut itself encodes the state that achieves the minimum free energy. Because of its deterministic nature and polynomial-time performance, this algorithm has the potential to allow for the ionization state of larger proteins to be discovered.

energy minimization↗

Distribution–Based Global Sensitivity Analysis in Hydrology

Global sensitivity analysis (GSA) is routinely used in academic setting to quantify the influence of input variability and uncertainty on predictions of a quantity of interest. Practical applications of GSA are hampered by its high computational cost, which arises from the need to run large (e.g., groundwater) models multiple times, and by its reliance on the analysis of variance, which formally requires input parameters to be uncorrelated. The former difficulty can be alleviated by replacing expensive models with inexpensive (e.g., polynomial) surrogates, while adoption of distribution-based (rather than variance-based) metrics can, in principle, overcome the latter but at significantly increased computational cost. To make use of distribution-based GSA feasible for regional-scale models with a large number of degrees of freedom, we supplement it with a surrogate model built with polynomial chaos expansions with analytically updated coefficients. Here, we demonstrate the computational efficiency of our algorithm on a case study dealing with evaluation of the effects of temperature variability on annual evapotranspiration at the regional scale.

42 ENGINEERING↗

Efficient multimode Wigner tomography

Abstract Advancements in quantum system lifetimes and control have enabled the creation of increasingly complex quantum states, such as those on multiple bosonic cavity modes. When characterizing these states, traditional tomography scales exponentially with the number of modes in both computational and experimental measurement requirement, which becomes prohibitive as the system size increases. Here, we implement a state reconstruction method whose sampling requirement instead scales polynomially with system size, and thus mode number, for states that can be represented within such a polynomial subspace. We demonstrate this improved scaling with Wigner tomography of multimode entangled W states of up to 4 modes on a 3D circuit quantum electrodynamics (cQED) system. This approach performs similarly in efficiency to existing matrix inversion methods for 2 modes, and demonstrates a noticeable improvement for 3 and 4 modes, with even greater theoretical gains at higher mode numbers.

Science & Technology - Other Topics↗

Quantum computational phase transition in combinatorial problems

Quantum Approximate Optimization algorithm (QAOA) aims to search for approximate solutions to discrete optimization problems with near-term quantum computers. As there are no algorithmic guarantee possible for QAOA to outperform classical computers, without a proof that bounded-error quantum polynomial time (BQP) ≠ nondeterministic polynomial time (NP), it is necessary to investigate the empirical advantages of QAOA. We identify a computational phase transition of QAOA when solving hard problems such as SAT—random instances are most difficult to train at a critical problem density. We connect the transition to the controllability and the complexity of QAOA circuits. Moreover, we find that the critical problem density in general deviates from the SAT-UNSAT phase transition, where the hardest instances for classical algorithms lies. Then, we show that the high problem density region, which limits QAOA’s performance in hard optimization problems (reachability deficits), is actually a good place to utilize QAOA: its approximation ratio has a much slower decay with the problem density, compared to classical approximate algorithms. Indeed, it is exactly in this region that quantum advantages of QAOA over classical approximate algorithms can be identified.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantum Zeno Monte Carlo for computing observables

The recent development of logical quantum processors marks a pivotal transition from the noisy intermediate-scale quantum (NISQ) era to the fault-tolerant quantum computing (FTQC) era. These devices have the potential to address classically challenging problems with polynomial computational time using quantum properties. However, they remain susceptible to noise, necessitating noise resilient algorithms. We introduce Quantum Zeno Monte Carlo (QZMC), a classical-quantum hybrid algorithm that demonstrates resilience to device noise and Trotter errors while showing polynomial computational cost for a gapped system. QZMC computes static and dynamic properties without requiring initial state overlap or variational parameters, offering reduced quantum circuit depth.

Han, Mancheon [Korea Institute for Advanced Study ↗

Rapidly convergent quantum Monte Carlo using a Chebyshev projector

The multireference coupled-cluster Monte Carlo (MR-CCMC) algorithm is a determinant-based quantum Monte Carlo (QMC) algorithm that is conceptually similar to Full Configuration Interaction QMC (FCIQMC). It has been shown to offer a balanced treatment of both static and dynamic correlation while retaining polynomial scaling, although application to large systems with significant strong correlation remained impractical. In this paper, we document recent algorithmic advances that enable rapid convergence and a more black-box approach to the multireference problem. These include a logarithmically scaling metric-tree-based excitation acceptance algorithm to search for determinants connected to the reference space at the desired excitation level and a symmetry-screening procedure for the reference space. We show that, for moderately sized reference spaces, the new search algorithm brings about an approximately 8-fold acceleration of one MR-CCMC iteration, while the symmetry screening procedure reduces the number of active reference space determinants with essentially no loss of accuracy. We also introduce a stochastic implementation of an approximate wall projector, which is the infinite imaginary time limit of the exponential projector, using a truncated expansion of the wall function in Chebyshev polynomials. Notably, this wall-Chebyshev projector can be used to accelerate any projector-based QMC algorithm. We show that it requires significantly fewer applications of the Hamiltonian to achieve the same statistical convergence. We benchmark these acceleration methods on the beryllium and carbon dimers, using initiator FCIQMC and MR-CCMC with basis sets up to cc-pVQZ quality.

Zhao, Zijun↗

Elucidating the local structure of Li 1+ x Al x Ti 2– x (PO 4 ) 3 and Li 3 Al x Ti 2– x (PO 4 ) 3 ( x = 0, 0.3) via total scattering

Li 1+x Al x Ti 2–x (PO 4 ) 3 (LATP) and Li 3 Al x Ti 2–x (PO 4 ) 3 (x = 0, 0.3) are promising candidates in all-solid-state batteries due to their high room temperature conductivity of 10 –3 S cm –1 and air- and moisture-stability. They also exhibit unusual thermal expansion properties, with Li 1+x Al x Ti 2–x (PO 4 ) 3 showing near-zero thermal expansion along the a axis while Li 3 Al x Ti 2–x (PO 4 ) 3 exhibits polynomial positive thermal expansion along the a axis and polynomial negative thermal expansion along the c axis. A crucial component to understanding these properties is understanding the local structure. Total scattering is a powerful analytical technique as it provides information on the long-range, average structure as well as the local structure. Here, we report the first X-ray and neutron total scattering experiments performed on Li 1+x Al x Ti 2–x (PO 4 ) 3 and Li 3 Al x Ti 2–x (PO 4 ) 3 (x = 0, 0.3). We show that the PO 4 and TiO 6 polyhedra experience very little expansion of the P/Ti–O bonds up to 800 °C, nor is there much expansion when the Li content increases significantly. The minor thermal expansion of the nearest-neighbor bonds of the polyhedra is revealed to be the reason behind the unusual thermal expansion properties, causing the near-zero thermal expansion along a in Li 1+x Al x Ti 2–x (PO 4 ) 3 and moving as whole units in Li 3 Al x Ti 2–x (PO 4 ) 3 . The structural robustness of the framework is also the reason for the increased conductivity as Li content increases, as the framework remains undistorted as Li content increases, permitting Li-ion mobility as the number of charge carriers increases. Finally, this suggests that phosphate-based framework materials beyond LATP would also be a good material space to explore for new Li-ion (and other ion-) conducting materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A mixed $\ell_1$ regularization approach for sparse simultaneous approximation of parameterized PDEs

We introduce and assess a novel sparse polynomial technique for the simultaneous approximation of parameterized partial differential equations (PDEs) with deterministic and stochastic inputs. Our approach treats the numerical solution as a jointly sparse reconstruction problem through the reformulation of the standard basis pursuit denoising, where the set of jointly sparse vectors is infinite. To achieve global reconstruction of sparse solutions to parameterized elliptic PDEs over both physical and parametric domains, we combine the standard measurement scheme developed for compressed sensing in the context of bounded orthonormal systems with a novel mixed-norm based $\ell_1$ regularization method that exploits both energy and sparsity. Moreover, we are able to prove that, with minimal sample complexity, error estimates comparable to the best $s$-term and quasi-optimal approximations are achievable, while requiring only {\em a priori} bounds on polynomial truncation error with respect to the energy norm.Finally, we perform extensive numerical experiments on several high-dimensional parameterized elliptic PDE models to demonstrate the superior recovery properties of the proposed approach.

97 MATHEMATICS AND COMPUTING↗

Proposed linear energy transfer areal detector for protons using radiochromic film

Radiation therapy depends on predictably and reliably delivering dose to tumors and sparing normal tissues. Protons with kinetic energy of a few hundred MeV can selectively deposit dose to deep seated tumors without an exit dose, unlike x-rays. The better dose distribution is attributed to a phenomenon known as the Bragg peak. The Bragg peak is due to relatively high energy deposition within a given distance or high Linear Energy Transfer (LET). In addition, biological response to radiation depends on the dose, dose rate, and localized energy deposition patterns or LET. At present, the LET can only be measured at a given fixed point and the LET spatial distribution can only be inferred from calculations. The goal of this study is to develop and test a method to measure LET over extended areas. Traditionally, radiochromic films are used to measure dose distribution but not for LET distribution. We report the first use of these films for measuring the spatial distribution of the LET deposited by protons. The radiochromic film sensitivity diminishes for large LET. A mathematical model correlating the film sensitivity and LET is presented to justify relating LET and radiochromic film relative sensitivity. Protons were directed parallel to radiochromic film sandwiched between solid water slabs. This study proposes the scaled-normalized difference (SND) between the Treatment Planning system (TPS) and measured dose as the metric describing the LET. The SND is correlated with a Monte Carlo (MC) calculation of the LET spatial distribution for a large range of SNDs. A polynomial fit between the SND and MC LET is generated for protons having a single range of 20 cm with narrow Bragg peak. Coefficients from these fitted polynomial fits were applied to measured proton dose distributions with a variety of ranges. An identical procedure was applied to the protons deposited from Spread Out Bragg Peak and modulated by 5 cm. Gamma analysis is a method for comparing the calculated LET with the LET measured using radiochromic film at the pixel level over extended areas. Failure rates using gamma analysis are calculated for areas in the dose distribution using parameters of 25% of MC LET and 3 mm. The processed dose distributions find 5%–10% failure rates for the narrow 12.5 and 15 cm proton ranges and 10%–15% for proton ranges of 15, 17.5, and 20 cm and modulated by 5 cm. It is found through gamma analysis that the measured proton energy deposition in radiochromic film and TPS can be used to determine LET. This modified film dosimetry provides an experimental areal LET measurement that can verify MC calculations, support LET point measurements, possibly enhance biologically based proton treatment planning, and determine the polymerization process within the radiochromic film.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Compact CH 4 sensor system based on a continuous-wave, low power consumption, room temperature interband cascade laser

A tunable diode laser absorption spectroscopy-based methane sensor, employing a dense-pattern multi-pass gas cell and a 3.3 µm, CW, DFB, room temperature interband cascade laser (ICL), is reported. The optical integration based on an advanced folded optical path design and an efficient ICL control system with appropriate electrical power management resulted in a CH 4 sensor with a small footprint (32 x 20 x 17 cm 3 ) and low-power consumption (6 W). Polynomial and least-squares fit algorithms are employed to remove the baseline of the spectral scan and retrieve CH 4 concentrations, respectively. An Allan-Werle deviation analysis shows that the measurement precision can reach 1.4 ppb for a 60 s averaging time. Continuous measurements covering a seven-day period were performed to demonstrate the stability and robustness of the reported CH 4 sensor system.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Communication: Fitting potential energy surfaces with fundamental invariant neural network

A more flexible neural network (NN) method using the fundamental invariants (FIs) as the input vector is proposed in the construction of potential energy surfaces for molecular systems involving identical atoms. Mathematically, FIs finitely generate the permutation invariant polynomial (PIP) ring. In combination with NN, fundamental invariant neural network (FI-NN) can approximate any function to arbitrary accuracy. Because FI-NN minimizes the size of input permutation invariant polynomials, it can efficiently reduce the evaluation time of potential energy, in particular for polyatomic systems. In this work, we provide the FIs for all possible molecular systems up to five atoms. Potential energy surfaces for OH{sub 3} and CH{sub 4} were constructed with FI-NN, with the accuracy confirmed by full-dimensional quantum dynamic scattering and bound state calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Prime factorization using magnonic holographic devices

Determining the prime factors of a given number $N$ is a problem that requires super-polynomial time for conventional digital computers. A polynomial-time algorithm was invented by Shor for quantum computers. In this paper, we present experimental data that demonstrate prime factorization using spin-wave interference but without quantum entanglement. Prime factorization includes three major steps. First, a general-type computer calculates the sequence of numbers m k mod($N$), where $N$ is the number to be factorized, m is a randomly chosen positive integer, and k = 1, 2, 3, 4, 5, 6…. Next, the period of the calculated sequence $r$ is determined by exploiting spin-wave interference. Finally, the general-type computer determines the primes based on the obtained $r$. The experiment for period finding was conducted on a six-terminal Y 3 Fe 2 (FeO 4 ) 3 device. We chose number 15 for testing and determined its primes using a sequence of measurements. The obtained experimental data for a micrometer-sized prototype aimed to demonstrate the benefits of using spin-wave devices to solve complex computational problems. Scalability is one of the major strengths inherent in this type of wave-based device, which may provide a route to nanometer-sized logic circuits. We discuss the physical and technological limitations of this approach, which define the maximum size of N and the computational speed. Although this classical approach cannot compete with the quantum algorithm in terms of efficiency, magnonic holographic devices can potentially be used as complementary logic units aimed at speeding up prime factorization for classical computers.

97 MATHEMATICS AND COMPUTING↗

Neural network based quasi-diabatic Hamiltonians with symmetry adaptation and a correct description of conical intersections

Previously, we demonstrated in a report that artificial neural networks (NNs) can be used to generate quasidiabatic Hamiltonians (H d ) that are capable of representing adiabatic energies, energy gradients, and derivative couplings. In this work, two additional issues are addressed. First, symmetry-adapted functions such as permutation invariant polynomials are introduced to account for complete nuclear permutation inversion symmetry. Second, a partially diagonalized representation is introduced to facilitate a better description of near degeneracy points. The diabatization of 1, 2 1 A states of NH 3 is used as an example. The NN fitting findings are compared to that of a previous fitting with symmetry adapted polynomials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Determination of acoustic nonlinearity parameters using thermal modulation of ultrasonic waves

This study presents a test method and its theoretical framework to determine the acoustic nonlinearity parameters (α, β, δ) of material using thermal modulation of ultrasonic waves. Temperature change-induced thermal strain excites the nonlinear response of the material and modulates the ultrasonic wave propagating in it. Experimental results showed a strong correlation between the relative wave velocity change and the temperature change. With a quadratic polynomial model, the acoustic nonlinearity parameters were obtained from the polynomial coefficients by curve fitting the experimental data. Their effects on thermal-induced velocity change were discussed. The parameters α, β, and δ govern the hysteretic gap, average slope, and curvature of the correlation curve, respectively. The proposed theory was validated on aluminum, steel, intact and damaged concrete samples. The obtained nonlinear parameters show reasonable agreement with values reported in the literature. Compared to other nonlinear acoustic methods using vibration or acoustic excitation, the thermal modulation method generates more uniform, slow changing, and larger strain field in the test sample. Employing the thermal effect as the driving force for nonlinearity instead of an undesired influencing factor, this method can measure the absolute values of α, β, and δ with good accuracy using a simple ultrasonic test setup.

Physics↗

Active learning for robust, high-complexity reactive atomistic simulations

Machine learned reactive force fields based on polynomial expansions have been shown to be highly effective for describing simulations involving reactive materials. Nevertheless, the highly flexible nature of these models can give rise to a large number of candidate parameters for complicated systems. In these cases, reliable parameterization requires a well-formed training set, which can be difficult to achieve through standard iterative fitting methods. In this paper, we present an active learning approach based on cluster analysis and inspired by Shannon information theory to enable semi-automated generation of informative training sets and robust machine learned force fields. The use of this tool is demonstrated for development of a model based on linear combinations of Chebyshev polynomials explicitly describing up to four-body interactions, for a chemically and structurally diverse system of C/O under extreme conditions. We show that this flexible training database management approach enables development of models exhibiting excellent agreement with Kohn–Sham density functional theory in terms of structure, dynamics, and speciation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Novel 2D velocity estimation method for large transient events in plasmas

Dynamics of fast transient events are challenging to be analyzed with high time resolution. Such events can occur in fusion plasmas such as the filaments during edge-localized modes (ELMs). Here, we present a robust method—the spatial displacement estimation—for estimating the displacements of structures with fast dynamics from high spatial and time resolution imaging diagnostics [e.g., gas-puff imaging (GPI)] with sampling time temporal resolution. First, a background suppression method is shown, which suppresses the slowly time-evolving and spatially non-uniform background in the signal. In the second step, a two-dimensional polynomial trend subtraction method is presented to tackle the remaining polynomial order trend in the signal. After performing these pre-processing steps, the spatial displacement of the propagating structure is estimated from the two-dimensional spatial cross-correlation coefficient function calculated between consecutive frames. The method is tested for its robustness and accuracy by simulated Gaussian events and spatially displaced random noise. An example application of the method is presented on propagating ELM filaments measured by the GPI system on the National Spherical Torus Experiment spherical tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗