Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ORTHOGONAL FUNCTION”

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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

SU-G-BRA-12: Development of An Intra-Fractional Motion Tracking and Dose Reconstruction System for Adaptive Stereotactic Body Radiation Therapy in High-Risk Prostate Cancer

Purpose: A clinical trial on stereotactic body radiation therapy (SBRT) for high-risk prostate cancer is undergoing at our institution. In addition to escalating dose to the prostate, we have increased dose to intra-prostatic lesions. Intra-fractional prostate motion deteriorates well planned radiation dose, especially for the small intra-prostatic lesions. To solve this problem, we have developed a motion tracking and 4D dose-reconstruction system to facilitate adaptive re-planning. Methods: Patients in the clinical trial were treated with VMAT using four arcs and 10 FFF beam. KV triggered x-ray projections were taken every 3 sec during delivery to acquire 2D projections of 3D anatomy at the direction orthogonal to the therapeutic beam. Each patient had three implanted prostate markers. Our developed system first determined 2D projection locations of these markers and then 3D prostate translation and rotation via 2D/3D registration of the markers. Using delivery log files, our GPU-based Monte Carlo tool (goMC) reconstructed dose corresponding to each triggered image. The calculated 4D dose distributions were further aggregated to yield the delivered dose. Results: We first tested each module in our system. MC dose engine were commissioned to our treatment planning system with dose difference of <0.5%. For motion tracking, 1789 kV projections from 7 patients were acquired. The 2D marker location error was <1 mm. For 3D motion tracking, root mean square (RMS) errors along LR, AP, and CC directions were 0.26mm, 0.36mm, and 0.01mm respectively in simulation studies and 1.99mm, 1.37mm, and 0.22mm in phantom studies. We also tested the entire system workflow. Our system was able to reconstruct delivered dose. Conclusion: We have developed a functional intra-fractional motion tracking and 4D dose re-construction system to support our clinical trial on adaptive high-risk prostate cancer SBRT. Comprehensive evaluations have shown the capability and accuracy of our system.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Reconfiguring DNA Nanotube Architectures via Selective Regulation of Terminating Structures

Molecular assemblies inside cells often undergo structural reconfiguration in response to stimuli to alter their function. Adaptive reconfiguration of cytoskeletal networks, for example, enables cellular shape change, movement, and cargo transport and plays a key role in driving complex processes such as division and differentiation. The cellular cytoskeleton is a self-assembling polymer network composed of simple filaments, so reconfiguration often occurs through the rearrangement of its component filaments’ connectivities. DNA nanotubes have emerged as promising building blocks for constructing programmable synthetic analogs of cytoskeletal networks. Nucleating seeds can control when and where nanotubes grow and capping structures can bind nanotube ends to stop growth. Such seeding and capping structure, collectively called termini, can organize nanotubes into larger architectures. However, these structures cannot be selectively activated or inactivated in response to specific stimuli to rearrange nanotube architectures, a key property of cytoskeletal networks. Here we demonstrate how selective regulation of the binding affinity of DNA nanotube termini for DNA nanotube monomers or nanotube ends can direct the reconfiguration of nanotube architectures. Using DNA hybridization and strand displacement reactions that specifically activate or inactivate four orthogonal nanotube termini, we demonstrate that nanotube architectures can be reconfigured by selective addition or removal of unique termini. Lastly, we show how terminus activation could be a sensitive detector and amplifier of a DNA sequence signal. These results could enable the development of adaptive and multifunctional materials or diagnostic tools.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

SnSe 1-x S x Alloys: Anisotropic Van der Waals Semiconductors with Tunable Bandgaps

Alloying is one of the main tools of bandgap engineering, allowing the tuning of crystal structure, lattice parameters, and electronic structure of 3D and 2D/layered semiconductors. Among the latter, it can play a key role in tailoring the properties of tin monochalcogenides, a class of van der Waals semiconductors of interest for optoelectronics, thermoelectrics, ferroelectrics, and valleytronics. Here, the study investigates the synthesis and properties of large flakes of the anion substitution alloys SnSe 1-x S x . Alloy flakes across a wide range of compositions are obtained systematically by repeated growth from the same mixed (SnS, SnSe) powder precursor. Combined experiment and theory show full miscibility for all compositions, along with tunable lattice constants, bandgaps, and vibrational modes. Atomic resolution imaging demonstrates the accumulation of S and Se in alternating layers in the SnSe 1-x S x unit cell, attributed to growth kinetics. Polarized Raman spectroscopy confirms anisotropic vibrational modes; the calculated and measured band structure shows systematic changes in the band edge energies and anisotropic electronic structure due to the anisotropic in-plane lattice of the monochalcogenides. Cathodoluminescence, finally, indicates that a unique configuration of two non-degenerate, direct valleys along orthogonal k-space directions persists all the way from SnS to SnSe, making SnSe 1-x S x alloys interesting for valleytronics.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Multi-layer Lanczos iteration approach to calculations of vibrational energies and dipole transition intensities for polyatomic molecules

We report a rigorous full dimensional quantum dynamics algorithm, the multi-layer Lanczos method, for computing vibrational energies and dipole transition intensities of polyatomic molecules without any dynamics approximation. The multi-layer Lanczos method is developed by using a few advanced techniques including the guided spectral transform Lanczos method, multi-layer Lanczos iteration approach, recursive residue generation method, and dipole-wavefunction contraction. The quantum molecular Hamiltonian at the total angular momentum J = 0 is represented in a set of orthogonal polyspherical coordinates so that the large amplitude motions of vibrations are naturally described. In particular, the algorithm is general and problem-independent. An application is illustrated by calculating the infrared vibrational dipole transition spectrum of CH{sub 4} based on the ab initio T8 potential energy surface of Schwenke and Partridge [Spectrochimica Acta, Part A 57, 887 (2001)] and the low-order truncated ab initio dipole moment surfaces of Yurchenko et al. [J. Mol. Spectrosc. 291, 69 (2013)]. A comparison with experiments is made. The algorithm is also applicable for Raman polarizability active spectra.

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

4D tumor centroid tracking using orthogonal 2D dynamic MRI: Implications for radiotherapy planning

Purpose: Current pretreatment, 4D imaging techniques are suboptimal in that they sample breathing motion over a very limited “snapshot” in time. Heretofore, long-duration, 4D motion characterization for radiotherapy planning, margin optimization, and validation have been impractical for safety reasons, requiring invasive markers imaged under x-ray fluoroscopy. To characterize 3D tumor motion and associated variability over durations more consistent with treatments, the authors have developed a practical dynamic MRI (dMRI) technique employing two orthogonal planes acquired in a continuous, interleaved fashion.Methods: 2D balanced steady-state free precession MRI was acquired continuously over 9–14 min at approximately 4 Hz in three healthy volunteers using a commercial 1.5 T system; alternating orthogonal imaging planes (sagittal, coronal, sagittal, etc.) were employed. The 2D in-plane pixel resolution was 2 × 2 mm{sup 2} with a 5 mm slice profile. Simultaneous with image acquisition, the authors monitored a 1D surrogate respiratory signal using a device available with the MRI system. 2D template matching-based anatomic feature registration, or tracking, was performed independently in each orientation. 4D feature tracking at the raw frame rate was derived using spline interpolation.Results: Tracking vascular features in the lung for two volunteers and pancreatic features in one volunteer, the authors have successfully demonstrated this method. Registration error, defined here as the difference between the sagittal and coronal tracking result in the SI direction, ranged from 0.7 to 1.6 mm (1σ) which was less than the acquired image resolution. Although the healthy volunteers were instructed to relax and breathe normally, significantly variable respiration was observed. To demonstrate potential applications of this technique, the authors subsequently explored the intrafraction stability of hypothetical tumoral internal target volumes and 3D spatial probability distribution functions. The surrogate respiratory information allowed the authors to show how this technique can be used to study correlations between internal and external (surrogate) information over these prolonged durations. However, compared against the gold standard of the time stamps in the dMRI frames, the temporal synchronization of the surrogate 1D respiratory information was shown to be likely unreliable.Conclusions: The authors have established viability of a novel and practical pretreatment, 4D tumor centroid tracking method employing a commercially available dynamic MRI sequence. Further developments from the vendor are likely needed to provide a reliably synchronized surrogate 1D respiratory signal, which will likely broaden the utility of this method in the pretreatment radiotherapy planning context.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Non-orthogonal configuration interaction for the calculation of multielectron excited states

We apply Non-orthogonal Configuration Interaction (NOCI) to molecular systems where multielectron excitations, in this case double excitations, play a substantial role: the linear polyenes and β-carotene. We demonstrate that NOCI when applied to systems with extended conjugation, provides a qualitatively correct wavefunction at a fraction of the cost of many other multireference treatments. We also present a new extension to this method allowing for purification of higher-order spin states by utilizing Generalized Hartree-Fock Slater determinants and the details for computing 〈S{sup 2}〉 for the ground and excited states.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Three-dimensional alignment of asymmetric-top molecules induced by polarization-shaped optimal laser pulses

Optimal control simulation is applied to numerically design nonresonant laser pulses that maximize the degrees of three-dimensional (3D) alignment of SO 2 using the lowest-order induced-dipole interaction. In our trials, combinations of more than two mutually orthogonal, linearly polarized subpulses are always obtained as the optimal solutions. Each subpulse in the optimal pulses impulsively excites the rotational wave packet. The optimal pulses effectively cooperate with the rotational dynamics up to only a few partial revival timings owing to the rotational dephasing that determines the effective control periods. The control mechanisms are interpreted in terms of the time derivatives of the expectation values of the squares of the direction cosines, that of the rotational energy, and the interplay between them. We find a special and important role of the last subpulses as they align the molecular axes using the interaction through the two smallest polarizability components, while the other subpulses excite the rotational wave packet mainly through the largest polarizability component. Lastly, the control pulses composed of the specified number of subpulses are also numerically optimized by actively utilizing the instantaneous penalty to systematically show the superiority of the use of more than two subpulses over that of two subpulses, the latter of which leads to the saturation of the degree of 3D alignment as a function of total fluence.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Transverse momentum-dependent heavy-quark fragmentation at next-to-leading order

The transverse momentum-dependent fragmentation functions (TMD FFs) of heavy (bottom and charm) quarks, which we recently introduced, are universal building blocks that enter predictions for a large number of observables involving final-state heavy quarks or hadrons. They enable the extension of fixed-order subtraction schemes to quasi-collinear limits, and are of particular interest in their own right as probes of the nonperturbative dynamics of hadronization. In this paper we calculate all TMD FFs involving heavy quarks and the associated TMD matrix element in heavy-quark effective theory (HQET) to next-to-leading order in the strong interaction. Our results confirm the renormalization properties, large-mass, and small-mass consistency relations predicted in our earlier work. We also derive and confirm a prediction for the large-z behavior of the heavy-quark TMD FF by extending, for the first time, the formalism of joint resummation to capture quark mass effects in heavy-quark fragmentation. Our final results in position space agree with those of a recent calculation by another group that used a highly orthogonal organization of singularities in the intermediate momentum-space steps, providing a strong independent cross check. As an immediate application, we present the complete quark mass dependence of the energy-energy correlator (EEC) in the back-to-back limit at $\mathcal{O}\left({\alpha}_s\right)$.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Multichannel Topological Kondo Effect

A Coulomb blockaded M-Majorana island coupled to normal metal leads realizes a novel type of Kondo effect where the effective impurity “spin” transforms under the orthogonal group SO⁡(M). The impurity spin stems from the nonlocal topological ground state degeneracy of the island and thus the effect is known as the topological Kondo effect. We introduce a physically motivated N-channel generalization of the topological Kondo model. Starting from the simplest case N =2, we conjecture a stable intermediate coupling fixed point and evaluate the resulting low-temperature impurity entropy. The impurity entropy indicates that an emergent Fibonacci anyon can be realized in the N =2 model. We also map the case N =2, M =4 to the conventional four-channel Kondo model and find the conductance at the intermediate fixed point. By using the perturbative renormalization group, we also analyze the large-N limit, where the fixed point moves to weak coupling. In the isotropic limit, we find an intermediate stable fixed point, which is stable to “exchange” coupling anisotropies, but unstable to channel anisotropy. We evaluate the fixed point impurity entropy and conductance to obtain experimentally observable signatures of our results. Here, in the large-N limit, we evaluate the full crossover function describing the temperature-dependent conductance.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2010

A Transportable Gravity Gradiometer Based on Atom Interferometry

A transportable atom interferometer-based gravity gradiometer has been developed at JPL to carry out measurements of Earth's gravity field at ever finer spatial resolutions, and to facilitate high-resolution monitoring of temporal variations in the gravity field from ground- and flight-based platforms. Existing satellite-based gravity missions such as CHAMP and GRACE measure the gravity field via precise monitoring of the motion of the satellites; i.e. the satellites themselves function as test masses. JPL's quantum gravity gradiometer employs a quantum phase measurement technique, similar to that employed in atomic clocks, made possible by recent advances in laser cooling and manipulation of atoms. This measurement technique is based on atomwave interferometry, and individual laser-cooled atoms are used as drag-free test masses. The quantum gravity gradiometer employs two identical atom interferometers as precision accelerometers to measure the difference in gravitational acceleration between two points (Figure 1). By using the same lasers for the manipulation of atoms in both interferometers, the accelerometers have a common reference frame and non-inertial accelerations are effectively rejected as common mode noise in the differential measurement of the gravity gradient. As a result, the dual atom interferometer-based gravity gradiometer allows gravity measurements on a moving platform, while achieving the same long-term stability of the best atomic clocks. In the laboratory-based prototype (Figure 2), the cesium atoms used in each atom interferometer are initially collected and cooled in two separate magneto-optic traps (MOTs). Each MOT, consisting of three orthogonal pairs of counter-propagating laser beams centered on a quadrupole magnetic field, collects up to 10(exp 9) atoms. These atoms are then launched vertically as in an atom fountain by switching off the magnetic field and introducing a slight frequency shift between pairs of lasers to create a moving rest frame for the trapped atoms. While still in this moving-frame molasses, the laser frequencies are further detuned from the atomic resonance (while maintaining this relative frequency shift) to cool the atom cloud's temperature to 2 K or below, corresponding to an rms velocity of less than 2 cm/s. After launch, the cold atoms undergo further state and velocity selection to prepare for atom interferometry. The atom interferometers are then realized using laser-induced stimulated Raman transitions to perform the necessary manipulations of each atom, and the resulting interferometer phase is measured using laser-induced fluorescence for state-normalized detection. More than 20 laser beams with independent controls of frequency, phase, and intensity are required for this measurement sequence. This instrument can facilitate the study of Earth's gravitational field from surface and air vehicles, as well as from space by allowing gravity mapping from a low-cost, single spacecraft mission. In addition, the operation of atom interferometer-based instruments in space offers greater sensitivity than is possible in terrestrial instruments due to the much longer interrogation times available in the microgravity environment. A space-based quantum gravity gradiometer has the potential to achieve sensitivities similar to the GRACE mission at long spatial wavelengths, and will also have resolution similar to GOCE for measurement at shorter length scales.

Yu, Nan↗
DOE OSTI2016

A novel software and conceptual design of the hardware platform for intensity modulated radiation therapy

Purpose: To deliver high quality intensity modulated radiotherapy (IMRT) using a novel generalized sparse orthogonal collimators (SOCs), the authors introduce a novel direct aperture optimization (DAO) approach based on discrete rectangular representation. Methods: A total of seven patients—two glioblastoma multiforme, three head & neck (including one with three prescription doses), and two lung—were included. 20 noncoplanar beams were selected using a column generation and pricing optimization method. The SOC is a generalized conventional orthogonal collimators with N leaves in each collimator bank, where N = 1, 2, or 4. SOC degenerates to conventional jaws when N = 1. For SOC-based IMRT, rectangular aperture optimization (RAO) was performed to optimize the fluence maps using rectangular representation, producing fluence maps that can be directly converted into a set of deliverable rectangular apertures. In order to optimize the dose distribution and minimize the number of apertures used, the overall objective was formulated to incorporate an L2 penalty reflecting the difference between the prescription and the projected doses, and an L1 sparsity regularization term to encourage a low number of nonzero rectangular basis coefficients. The optimization problem was solved using the Chambolle–Pock algorithm, a first-order primal–dual algorithm. Performance of RAO was compared to conventional two-step IMRT optimization including fluence map optimization and direct stratification for multileaf collimator (MLC) segmentation (DMS) using the same number of segments. For the RAO plans, segment travel time for SOC delivery was evaluated for the N = 1, N = 2, and N = 4 SOC designs to characterize the improvement in delivery efficiency as a function of N. Results: Comparable PTV dose homogeneity and coverage were observed between the RAO and the DMS plans. The RAO plans were slightly superior to the DMS plans in sparing critical structures. On average, the maximum and mean critical organ doses were reduced by 1.94% and 1.44% of the prescription dose. The average number of delivery segments was 12.68 segments per beam for both the RAO and DMS plans. The N = 2 and N = 4 SOC designs were, on average, 1.56 and 1.80 times more efficient than the N = 1 SOC design to deliver. The mean aperture size produced by the RAO plans was 3.9 times larger than that of the DMS plans. Conclusions: The DAO and dose domain optimization approach enabled high quality IMRT plans using a low-complexity collimator setup. The dosimetric quality is comparable or slightly superior to conventional MLC-based IMRT plans using the same number of delivery segments. The SOC IMRT delivery efficiency can be significantly improved by increasing the leaf numbers, but the number is still significantly lower than the number of leaves in a typical MLC.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

Reduced Order Modeling for Prediction and Control of Large-Scale Systems.

This report describes work performed from June 2012 through May 2014 as a part of a Sandia Early Career Laboratory Directed Research and Development (LDRD) project led by the first author. The objective of the project is to investigate methods for building stable and efficient proper orthogonal decomposition (POD)/Galerkin reduced order models (ROMs): models derived from a sequence of high-fidelity simulations but having a much lower computational cost. Since they are, by construction, small and fast, ROMs can enable real-time simulations of complex systems for onthe- spot analysis, control and decision-making in the presence of uncertainty. Of particular interest to Sandia is the use of ROMs for the quantification of the compressible captive-carry environment, simulated for the design and qualification of nuclear weapons systems. It is an unfortunate reality that many ROM techniques are computationally intractable or lack an a priori stability guarantee for compressible flows. For this reason, this LDRD project focuses on the development of techniques for building provably stable projection-based ROMs. Model reduction approaches based on continuous as well as discrete projection are considered. In the first part of this report, an approach for building energy-stable Galerkin ROMs for linear hyperbolic or incompletely parabolic systems of partial differential equations (PDEs) using continuous projection is developed. The key idea is to apply a transformation induced by the Lyapunov function for the system, and to build the ROM in the transformed variables. It is shown that, for many PDE systems including the linearized compressible Euler and linearized compressible Navier-Stokes equations, the desired transformation is induced by a special inner product, termed the “symmetry inner product”. Attention is then turned to nonlinear conservation laws. A new transformation and corresponding energy-based inner product for the full nonlinear compressible Navier-Stokes equations is derived, and it is demonstrated that if a Galerkin ROM is constructed in this inner product, the ROM system energy will be bounded in a way that is consistent with the behavior of the exact solution to these PDEs, i.e., the ROM will be energy-stable. The viability of the linear as well as nonlinear continuous projection model reduction approaches developed as a part of this project is evaluated on several test cases, including the cavity configuration of interest in the targeted application area. In the second part of this report, some POD/Galerkin approaches for building stable ROMs using discrete projection are explored. It is shown that, for generic linear time-invariant (LTI) systems, a discrete counterpart of the continuous symmetry inner product is a weighted L2 inner product obtained by solving a Lyapunov equation. This inner product was first proposed by Rowley et al., and is termed herein the “Lyapunov inner product“. Comparisons between the symmetry inner product and the Lyapunov inner product are made, and the performance of ROMs constructed using these inner products is evaluated on several benchmark test cases. Also in the second part of this report, a new ROM stabilization approach, termed “ROM stabilization via optimization-based eigenvalue reassignment“, is developed for generic LTI systems. At the heart of this method is a constrained nonlinear least-squares optimization problem that is formulated and solved numerically to ensure accuracy of the stabilized ROM. Numerical studies reveal that the optimization problem is computationally inexpensive to solve, and that the new stabilization approach delivers ROMs that are stable as well as accurate. Summaries of “lessons learned“ and perspectives for future work motivated by this LDRD project are provided at the end of each of the two main chapters.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

SU-E-T-603: PBS Prostate Plan Robustness: A Tool for Patient Specific Setup Tolerance

Purpose: Fiducial markers are commonly used for setup of prostate patients using orthogonal radiographs. After aligned with the markers, the displacement of the bony anatomy relative to the planned DRR can be up to 10 mm. Such offset can potentially have significant dosimetric effects because it changes the radiological path length of protons in differing amounts of bone. It is imperative to develop a method to evaluate its impact on target coverage and hence establish patient specific setup tolerance for prostate proton PBS treatment. Methods: Prostate patients were planned in RayStation according to the PCG protocol with bi-lateral beams. The primary planning objectives are: (1) 100% of CTV receives full prescription dose; (2) 98% of the prescription dose covers at least 98% of the PTV; (3) OARs meet criteria per protocol. For each patient 108 dose perturbations were automatically generated using an in-house script, which considered the isocenter shifting in S-I and A-P directions (up to ±15 mm with an interval of 6mm) as well as the range uncertainty (±3.5%). The target coverage was evaluated on the contour shifted along with prostate to mimic the daily treatment. Results: The minimum CTV coverage as a function of offsets in S-I and A-P directions is presented in a 2D contour map. The offsets along A-P direction generally have greater impact than along S-I direction. Both the CTV D98%>98% or CTV V98%>98% are achievable for most patients if the offset is <10 mm in either direction despite of range uncertainties. Conclusion: We developed a method to evaluate the plan robustness for proton PBS prostate treatment. It can provide patient specific setup tolerance of bony structure offset. For our current planning approach, a 1 cm displacement is acceptable. This approach can be generalized to other target structures that move relative to bony anatomy.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2015

Compressed-sensing-based content-driven hierarchical reconstruction: Theory and application to C-arm cone-beam tomography

Purpose: This paper addresses the reconstruction of x-ray cone-beam computed tomography (CBCT) for interventional C-arm systems. Subsampling of CBCT is a significant issue with C-arms due to their slow rotation and to the low frame rate of their flat panel x-ray detectors. The aim of this work is to propose a novel method able to handle the subsampling artifacts generally observed with analytical reconstruction, through a content-driven hierarchical reconstruction based on compressed sensing. Methods: The central idea is to proceed with a hierarchical method where the most salient features (high intensities or gradients) are reconstructed first to reduce the artifacts these features induce. These artifacts are addressed first because their presence contaminates less salient features. Several hierarchical schemes aiming at streak artifacts reduction are introduced for C-arm CBCT: the empirical orthogonal matching pursuit approach with the ℓ{sub 0} pseudonorm for reconstructing sparse vessels; a convex variant using homotopy with the ℓ{sub 1}-norm constraint of compressed sensing, for reconstructing sparse vessels over a nonsparse background; homotopy with total variation (TV); and a novel empirical extension to nonlinear diffusion (NLD). Such principles are implemented with penalized iterative filtered backprojection algorithms. For soft-tissue imaging, the authors compare the use of TV and NLD filters as sparsity constraints, both optimized with the alternating direction method of multipliers, using a threshold for TV and a nonlinear weighting for NLD. Results: The authors show on simulated data that their approach provides fast convergence to good approximations of the solution of the TV-constrained minimization problem introduced by the compressed sensing theory. Using C-arm CBCT clinical data, the authors show that both TV and NLD can deliver improved image quality by reducing streaks. Conclusions: A flexible compressed-sensing-based algorithmic approach is proposed that is able to accommodate for a wide range of constraints. It is successfully applied to C-arm CBCT images that may not be so well approximated by piecewise constant functions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

ELSI — An open infrastructure for electronic structure solvers

Routine applications of electronic structure theory to molecules and periodic systems need to compute the electron density from given Hamiltonian and, in case of non-orthogonal basis sets, overlap matrices. System sizes can range from few to thousands or, in some examples, millions of atoms. Different discretization schemes (basis sets) and different system geometries (finite non-periodic vs. infinite periodic boundary conditions) yield matrices with different structures. The ELectronic Structure Infrastructure (ELSI) project provides an open-source software interface to facilitate the implementation and optimal use of high-performance solver libraries covering cubic scaling eigensolvers, linear scaling density-matrix-based algorithms, and other reduced scaling methods in between. In this paper, we present recent improvements and developments inside ELSI, mainly covering (1) new solvers connected to the interface, (2) matrix layout and communication adapted for parallel calculations of periodic and/or spin-polarized systems, (3) routines for density matrix extrapolation in geometry optimization and molecular dynamics calculations, and (4) general utilities such as parallel matrix I/O and JSON output. The ELSI interface has been integrated into four electronic structure code projects (DFTB+, DGDFT, FHI-aims, SIESTA), allowing us to rigorously benchmark the performance of the solvers on an equal footing. Based on results of a systematic set of large-scale benchmarks performed with Kohn–Sham density-functional theory and density-functional tight-binding theory, we identify factors that strongly affect the efficiency of the solvers, and propose a decision layer that assists with the solver selection process. As a result, we describe a reverse communication interface encoding matrix-free iterative solver strategies that are amenable, e.g., for use with planewave basis sets.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2009

Target-Tracking Camera for a Metrology System

An analog electronic camera that is part of a metrology system measures the varying direction to a light-emitting diode that serves as a bright point target. In the original application for which the camera was developed, the metrological system is used to determine the varying relative positions of radiating elements of an airborne synthetic aperture-radar (SAR) antenna as the airplane flexes during flight; precise knowledge of the relative positions as a function of time is needed for processing SAR readings. It has been common metrology system practice to measure the varying direction to a bright target by use of an electronic camera of the charge-coupled-device or active-pixel-sensor type. A major disadvantage of this practice arises from the necessity of reading out and digitizing the outputs from a large number of pixels and processing the resulting digital values in a computer to determine the centroid of a target: Because of the time taken by the readout, digitization, and computation, the update rate is limited to tens of hertz. In contrast, the analog nature of the present camera makes it possible to achieve an update rate of hundreds of hertz, and no computer is needed to determine the centroid. The camera is based on a position-sensitive detector (PSD), which is a rectangular photodiode with output contacts at opposite ends. PSDs are usually used in triangulation for measuring small distances. PSDs are manufactured in both one- and two-dimensional versions. Because it is very difficult to calibrate two-dimensional PSDs accurately, the focal-plane sensors used in this camera are two orthogonally mounted one-dimensional PSDs.

Liebe, Carl↗
DOE OSTI2013

Extending Romanovski polynomials in quantum mechanics

Some extensions of the (third-class) Romanovski polynomials (also called Romanovski/pseudo-Jacobi polynomials), which appear in bound-state wavefunctions of rationally extended Scarf II and Rosen-Morse I potentials, are considered. For the former potentials, the generalized polynomials satisfy a finite orthogonality relation, while for the latter an infinite set of relations among polynomials with degree-dependent parameters is obtained. Both types of relations are counterparts of those known for conventional polynomials. In the absence of any direct information on the zeros of the Romanovski polynomials present in denominators, the regularity of the constructed potentials is checked by taking advantage of the disconjugacy properties of second-order differential equations of Schrödinger type. It is also shown that on going from Scarf I to Scarf II or from Rosen-Morse II to Rosen-Morse I potentials, the variety of rational extensions is narrowed down from types I, II, and III to type III only.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Accurate quantum dynamics calculations using symmetrized Gaussians on a doubly dense Von Neumann lattice

In a series of earlier articles [B. Poirier, J. Theor. Comput. Chem. 2, 65 (2003); B. Poirier and A. Salam, J. Chem. Phys. 121, 1690 (2004); and ibid. 121, 1704 (2004)], a new method was introduced for performing exact quantum dynamics calculations. The method uses a 'weylet' basis set (orthogonalized Weyl-Heisenberg wavelets) combined with phase space truncation, to defeat the exponential scaling of CPU effort with system dimensionality-the first method ever able to achieve this long-standing goal. Here, we develop another such method, which uses a much more convenient basis of momentum-symmetrized Gaussians. Despite being non-orthogonal, symmetrized Gaussians are collectively local, allowing for effective phase space truncation. A dimension-independent code for computing energy eigenstates of both coupled and uncoupled systems has been created, exploiting massively parallel algorithms. Results are presented for model isotropic uncoupled harmonic oscillators and coupled anharmonic oscillators up to 27 dimensions. These are compared with the previous weylet calculations (uncoupled harmonic oscillators up to 15 dimensions), and found to be essentially just as efficient. Coupled system results are also compared to corresponding exact results obtained using a harmonic oscillator basis, and also to approximate results obtained using first-order perturbation theory up to the maximum dimensionality for which the latter may be feasibly obtained (four dimensions).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗