Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Recursion”

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 649 records · Page 36

An efficient method for energy levels calculation using full symmetry and exact kinetic energy operator: Tetrahedral molecules

A simultaneous use of the full molecular symmetry and of an exact kinetic energy operator (KEO) is of key importance for accurate predictions of vibrational levels at a high energy range from a potential energy surface (PES). An efficient method that permits a fast convergence of variational calculations would allow iterative optimization of the PES parameters using experimental data. In this work, we propose such a method applied to tetrahedral AB{sub 4} molecules for which a use of high symmetry is crucial for vibrational calculations. A symmetry-adapted contracted angular basis set for six redundant angles is introduced. Simple formulas using this basis set for explicit calculation of the angular matrix elements of KEO and PES are reported. The symmetric form (six redundant angles) of vibrational KEO without the sin(q){sup −2} type singularity is derived. The efficient recursive algorithm based on the tensorial formalism is used for the calculation of vibrational matrix elements. A good basis set convergence for the calculations of vibrational levels of the CH{sub 4} molecule is demonstrated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Efficient block processing of long duration biotelemetric brain data for health care monitoring

In real time clinical environment, the brain signals which doctor need to analyze are usually very long. Such a scenario can be made simple by partitioning the input signal into several blocks and applying signal conditioning. This paper presents various block based adaptive filter structures for obtaining high resolution electroencephalogram (EEG) signals, which estimate the deterministic components of the EEG signal by removing noise. To process these long duration signals, we propose Time domain Block Least Mean Square (TDBLMS) algorithm for brain signal enhancement. In order to improve filtering capability, we introduce normalization in the weight update recursion of TDBLMS, which results TD-B-normalized-least mean square (LMS). To increase accuracy and resolution in the proposed noise cancelers, we implement the time domain cancelers in frequency domain which results frequency domain TDBLMS and FD-B-Normalized-LMS. Finally, we have applied these algorithms on real EEG signals obtained from human using Emotive Epoc EEG recorder and compared their performance with the conventional LMS algorithm. The results show that the performance of the block based algorithms is superior to the LMS counter-parts in terms of signal to noise ratio, convergence rate, excess mean square error, misadjustment, and coherence.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Calculation of the molecular integrals with the range-separated correlation factor

Explicitly correlated quantum chemical calculations require calculations of five types of two-electron integrals beyond the standard electron repulsion integrals. We present a novel scheme, which utilises general ideas of the McMurchie-Davidson technique, to compute these integrals when the so-called “range-separated” correlation factor is used. This correlation factor combines the well-known short range behaviour resulting from the electronic cusp condition, with the exact long-range asymptotics derived for the helium atom [Lesiuk, Jeziorski, and Moszynski, J. Chem. Phys. 139, 134102 (2013)]. Almost all steps of the presented procedure are formulated recursively, so that an efficient implementation and control of the precision are possible. Additionally, the present formulation is very flexible and general, and it allows for use of an arbitrary correlation factor in the electronic structure calculations with minor or no changes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of collisions on conservation laws in gyrokinetic field theory

Effects of collisions on conservation laws for toroidal plasmas are investigated based on the gyrokinetic field theory. Associating the collisional system with a corresponding collisionless system at a given time such that the two systems have the same distribution functions and electromagnetic fields instantaneously, it is shown how the collisionless conservation laws derived from Noether's theorem are modified by the collision term. Effects of the external source term added into the gyrokinetic equation can be formulated similarly with the collisional effects. Particle, energy, and toroidal momentum balance equations including collisional and turbulent transport fluxes are systematically derived using a novel gyrokinetic collision operator, by which the collisional change rates of energy and canonical toroidal angular momentum per unit volume in the gyrocenter space can be given in the conservative forms. The ensemble-averaged transport equations of particles, energy, and toroidal momentum given in the present work are shown to include classical, neoclassical, and turbulent transport fluxes which agree with those derived from conventional recursive formulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A simple extrapolation of thermodynamic perturbation theory to infinite order

Recent analyses of the third and fourth order perturbation contributions to the equations of state for square well spheres and Lennard-Jones chains show trends that persist across orders and molecular models. In particular, the ratio between orders (e.g., A{sub 3}/A{sub 2}, where A{sub i} is the ith order perturbation contribution) exhibits a peak when plotted with respect to density. The trend resembles a Gaussian curve with the peak near the critical density. This observation can form the basis for a simple recursion and extrapolation from the highest available order to infinite order. The resulting extrapolation is analytic and therefore cannot fully characterize the critical region, but it remarkably improves accuracy, especially for the binodal curve. Whereas a second order theory is typically accurate for the binodal at temperatures within 90% of the critical temperature, the extrapolated result is accurate to within 99% of the critical temperature. In addition to square well spheres and Lennard-Jones chains, we demonstrate how the method can be applied semi-empirically to the Perturbed Chain - Statistical Associating Fluid Theory (PC-SAFT)

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fault tolerant filtering and fault detection for quantum systems driven by fields in single photon states

The purpose of this paper is to solve the fault tolerant filtering and fault detection problem for a class of open quantum systems driven by a continuous-mode bosonic input field in single photon states when the systems are subject to stochastic faults. Optimal estimates of both the system observables and the fault process are simultaneously calculated and characterized by a set of coupled recursive quantum stochastic differential equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On regularizing the ML-MCTDH equations of motion

In a recent paper [H.-D. Meyer and H. Wang, J. Chem. Phys. 148, 124105 (2018)], we have examined the regularization of the equations of motion (EOMs) of the multiconfiguration time-dependent Hartree (MCTDH) approach. We could show that the standard regularization scheme used by almost all researchers in the field is not optimal. The improved regularization allows for larger values of the regularization parameter ϵ, is less sensitive to the actual choice of ϵ, and performs the rotation of initially unoccupied single-particle functions into the “correct” direction in Hilbert space much faster than the old scheme. Here, the latter point increases both the accuracy and efficiency of time propagation for challenging problems. For simple problems, the new scheme requires some additional numerical work as compared with the old scheme, ranging from negligible to almost doubling the total numerical labor. For demanding problems, on the other hand, the additional numerical work of the new scheme is often overcompensated by less steps taken by the integrator. In the present paper, we generalize the new regularization scheme to the multi-layer (ML) extension of MCTDH. Although the principle idea of the new regularization scheme remains unaltered, it was not obvious how the new scheme should be implemented into ML-MCTDH. The ML-MCTDH EOMs are much more complicated than the MCTDH ones, and for optimal numerical performance it was necessary to derive a recursive algorithm for implementing the new regularization scheme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Adaptive inverse ray-tracing for accurate and efficient modeling of cross beam energy transfer in hydrodynamics simulations

Integrated hydrodynamics simulations of inertial confinement fusion rely on reduced physics models. To reproduce experimental trends, these models often feature tuning parameters, but this comes with a risk: the over-tuning of one model can hide physics inadequacies in another. The ray-based models of cross-beam-energy transfer (CBET) represent this risk. In this paper, we present an accurate and efficient model of CBET suitable for inline implementation in 3D hydrodynamics simulations. Inverse Ray Tracing (IRT) is used to compute the ray field in a 3D permittivity profile described on an unstructured tetrahedral mesh using the Inline Field Reconstruction and Interaction using Inverse Tracing framework. CBET is accounted for through perturbations to the permittivity associated with ion acoustic waves driven by the overlapped fields. Large gradients in the permittivity are resolved by coupling the IRT to a recursive Adaptive Mesh Refinement (AMR) algorithm. The use of AMR also allows for the resolution of caustics, with accurate field reconstruction performed using the Etalon integral method. Comparisons of the model with wave-based solutions from the Laser Plasma Simulation Environment demonstrate its ability to control energy conservation and gain convergence through the AMR depth only, without the use of ad hoc physical models or artificial tuning parameters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Connected three-body terms in single-reference unitary many-body theories: Iterative and perturbative approximations

This work introduces various approaches to include connected three-body terms in unitary many-body theories, focusing on the driven similarity renormalization group (DSRG). Starting from the least approximate method—the linearized DSRG truncated to one-, two-, and three-body operators [LDSRG(3)]—we develop several approximate LDSRG(3) models with reduced computational cost. Through a perturbative analysis, we motivate a family of iterative LDSRG(3)-n and -n' (n = 1, 2, 3, 4) methods that contain a subset of the LDSRG(3) diagrams. Among these variants, the LDSRG(3)-2 scheme has the same computational complexity of coupled cluster theory with singles, doubles, and triples (CCSDT), but it outperforms CCSDT in the accuracy of the predicted correlation energies. We also propose and implement two perturbative triples corrections based on the linearized DSRG truncated to one- and two-body operators augmented with recursive semi-quadratic commutators [qDSRG(2)]. Overall, the resulting qDSRG(2)+(T) approach matches the accuracy of the “gold-standard” coupled cluster theory with singles, doubles, and perturbative triples model on the energetics of twenty-eight closed-shell atoms and small molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Finite temperature auxiliary field quantum Monte Carlo in the canonical ensemble

Finite temperature auxiliary field-based quantum Monte Carlo methods, including determinant quantum Monte Carlo and Auxiliary Field Quantum Monte Carlo (AFQMC), have historically assumed pivotal roles in the investigation of the finite temperature phase diagrams of a wide variety of multidimensional lattice models and materials. Despite their utility, however, these techniques are typically formulated in the grand canonical ensemble, which makes them difficult to apply to condensates such as superfluids and difficult to benchmark against alternative methods that are formulated in the canonical ensemble. Working in the grand canonical ensemble is furthermore accompanied by the increased overhead associated with having to determine the chemical potentials that produce desired fillings. Given this backdrop, in this work, we present a new recursive approach for performing AFQMC simulations in the canonical ensemble that does not require knowledge of chemical potentials. To derive this approach, we exploit the convenient fact that AFQMC solves the many-body problem by decoupling many-body propagators into integrals over one-body problems to which non-interacting theories can be applied. We benchmark the accuracy of our technique on illustrative Bose and Fermi–Hubbard models and demonstrate that it can converge more quickly to the ground state than grand canonical AFQMC simulations. We believe that our novel use of HS-transformed operators to implement algorithms originally derived for non-interacting systems will motivate the development of a variety of other methods and anticipate that our technique will enable direct performance comparisons against other many-body approaches formulated in the canonical ensemble.

Shen, Tong (ORCID:0000000206292468)↗

Efficient evaluation of the Breit operator in the Pauli spinor basis

The frequency-independent Coulomb–Breit operator gives rise to the most accurate treatment of two-electron interaction in the non-quantum-electrodynamics regime. The Breit interaction in the Coulomb gauge consists of magnetic and gauge contributions. The high computational cost of the gauge term limits the application of the Breit interaction in relativistic molecular calculations. In this work, we apply the Pauli component integral–density matrix contraction scheme for gauge interaction with a maximum spin- and component separation scheme. We also present two different computational algorithms for evaluating gauge integrals. One is the generalized Obara–Saika algorithm, where the Laplace transformation is used to transform the gauge operator into Gaussian functions and the Obara–Saika recursion is used for reducing the angular momentum. The other algorithm is the second derivative of Coulomb interaction evaluated with Rys-quadrature. This work improves the efficiency of performing Dirac–Hartree–Fock with the variational treatment of Breit interaction for molecular systems. In conclusion, we use this formalism to examine relativistic trends in the Periodic Table and analyze the relativistic two-electron interaction contributions in heavy-element complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Drift-ordered fluid vorticity equation with energy consistency

Although drift-ordered fluid models are widely applied in tokamak edge turbulence simulations, the models used are acknowledged not to conserve energy or even electrical charge. The present paper aims to remove many of the existing pitfalls in drift-fluid models, however, with the objective of finding a solution simple enough to be implemented in numerical applications. Our main result is an improved version of the drift-Braginskii equations involving a generalized vorticity function. In the new drift-Braginskii system, the quasi-neutrality condition translates into a transport equation for a generalized vorticity, expressed in conservation form, and related to the total mass-weighted circulation. It is found that kinetic energy conservation can be achieved if the polarization flow is defined recursively. In conclusion, the resulting model conserves the kinetic energy associated with E×B and diamagnetic flows and retains the associated perpendicular kinetic energy flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hybrid programming-model strategies for GPU offloading of electronic structure calculation kernels

To address the challenge of performance portability and facilitate the implementation of electronic structure solvers, we developed the basic matrix library (BML) and Parallel, Rapid O(N), and Graph-based Recursive Electronic Structure Solver (PROGRESS) library. The BML implements linear algebra operations necessary for electronic structure kernels using a unified user interface for various matrix formats (dense and sparse) and architectures (CPUs and GPUs). Focusing on density functional theory and tight-binding models, PROGRESS implements several solvers for computing the single-particle density matrix and relies on BML. In this paper, we describe the general strategies used for these implementations on various computer architectures, using OpenMP target functionalities on GPUs, in conjunction with third-party libraries to handle performance critical numerical kernels. In this study, we demonstrate the portability of this approach and its performance in benchmark problems.

36 MATERIALS SCIENCE↗

Polariton spectra under the collective coupling regime. I. Efficient simulation of linear spectra and quantum dynamics

We outline two general theoretical techniques to simulate polariton quantum dynamics and optical spectra under the collective coupling regimes described by a Holstein–Tavis–Cummings (HTC) model Hamiltonian. The first one takes advantage of sparsity of the HTC Hamiltonian, which allows one to reduce the cost of acting polariton Hamiltonian onto a state vector to the linear order of the number of states, instead of the quadratic order. The second one is applying the well-known Chebyshev series expansion approach for quantum dynamics propagation and to simulate the polariton dynamics in the HTC system; this approach allows us to use a much larger time step for propagation and only requires a few recursive operations of the polariton Hamiltonian acting on state vectors. These two theoretical approaches are general and can be applied to any trajectory-based non-adiabatic quantum dynamics methods. We apply these two techniques with our previously developed Lindblad-partially linearized density matrix approach to simulate the linear absorption spectra of the HTC model system, with both inhomogeneous site energy disorders and dipolar orientational disorders. Our numerical results agree well with the previous analytic and numerical work.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An algebraic convolution formulation for multiple-scattering correction in small-angle neutron scattering

Multiple scattering in small-angle neutron scattering (SANS) redistributes spectral weight and distorts structural interpretation, particularly for thick or strongly scattering samples. We develop a finite-dimensional spectral desmearing framework that corrects multiple scattering without resorting to integral transforms or model-dependent extrapolation. The primary intensity is expanded in an orthonormal basis adapted to the isotropic transverse-momentum measure, under which convolution reduces to a recursive tensor contraction, allowing the Poisson-weighted multiple-scattering series to be evaluated directly in a finite-dimensional basis representation. This formulation yields a stable forward–inverse mapping between apparent and primary spectra. Numerical tests demonstrate convergence under repeated convolution and accurate recovery of the single-scattering intensity. Application to SANS measurements collected at multiple neutron facilities, including the Spallation Neutron Source, the High Flux Isotope Reactor, and the Institut Laue-Langevin, shows the quantitative reconstruction of the underlying primary spectrum across a wide range of transmission conditions, including strongly attenuating samples. Here, the method provides a stable, model-agnostic framework for multiple-scattering correction in SANS and enables consistent structural interpretation across instruments and scattering regimes.

Tung, Chi-Huan [Oak Ridge National Laboratory (ORN↗

Asymptotic method for constructing a dynamic frequency response of a laser gyro

We consider an equation describing the dynamics of the phase of a beat signal at the output of a ring laser gyro in the presence of an alternating-sign frequency biasing in the form of a superposition of two meanders with different amplitudes and periods. It is assumed that the amplitude of one of the meanders is much larger than that of the other as well as than the half-widths of the static lock-in and the measured difference of eigenfrequencies of a ring resonator. A method is proposed for approximately calculating the beat frequency as a function of the measured difference between the eigenrequencies of the resonator. The result of applying this method represents analytically a recursive algorithm. Its implementation on a computer has made it possible to construct a dynamic frequency (output) response of a laser gyro with a biasing of given form. (laser applications and other topics in quantum electronics)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Tailored fatty acid synthesis via dynamic control of fatty acid elongation

Medium-chain fatty acids (MCFAs, 4-12 carbons) are valuable as precursors to industrial chemicals and biofuels, but are not canonical products of microbial fatty acid synthesis. We engineered microbial production of the full range of even-and odd-chain-length MCFAs and found that MCFA production is limited by rapid, irreversible elongation of their acyl-ACP precursors. To address this limitation, we programmed an essential ketoacyl synthase to degrade in response to a chemical inducer, thereby slowing acyl-ACP elongation and redirecting flux from phospholipid synthesis to MCFA production. Our results show that induced protein degradation can be used to dynamically alter metabolic flux, and thereby increase the yield of a desired compound. The strategy reported herein should be widely useful in a range of metabolic engineering applications in which essential enzymes divert flux away from a desired product, as well as in the production of polyketides, bioplastics, and other recursively synthesized hydrocarbons for which chain-length control is desired.

Torella, JP↗

The all configuration mean energy multiconfiguration self-consistent-field method. I. Equal configuration weights

The All Configuration Mean Energy (ACME) conditions are a special case of state averaging for Multiconfigurational Self-Consistent-Field (MCSCF) orbital optimisation. The method is formulated using the Graphical Unitary Group Approach (GUGA) in which the Configuration State Function (CSF) basis is represented as walks within a Shavitt graph. This graphical formulation leads to efficient recursive algorithms for the energy and reduced density matrices (RDM) that are independent of the CSF dimension and that scale only as O(n2) where n is the number of occupied orbitals. The Hamiltonian matrix diagonalization step is obviated and the CSF expansion coefficients are neither referenced nor required. This allows MCSCF orbital optimisation to be performed for essentially unlimited numbers of active orbitals and arbitrarily large CSF expansions. The discussion includes various types of CSF expansion spaces, the partitioning of the essential and redundant orbital optimisation parameters, the computation of the spin-density, and the formulation of state-specific analytic gradients and nonadiabatic coupling for high-level electronic structure methods that use the ACME MCSCF orbitals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗