Search NASA⌕ Search

SEARCH · Search NASA

Results for “selected inversion”

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 19 records

PSelInv—A Distributed Memory Parallel Algorithm for Selected Inversion: The Symmetric Case

We describe an efficient parallel implementation of the selected inversion algorithm for distributed memory computer systems, which we call PSelInv. The PSelInv method computes selected elements of a general sparse matrix Athat can be decomposed as A= LU, where L is lower triangular and U is upper triangular. The implementation described in this article focuses on the case of sparse symmetric matrices. It contains an interface that is compatible with the distributed memory parallel sparse direct factorization SuperLU-DIST. Yet, the underlying data structure and design of PSelInv allows it to be easily combined with other factorization routines, such as PARDISO. We discuss general parallelization strategies such as data and task distribution schemes. Specifically, we describe how to exploit the concurrency exposed by the elimination tree associated with the LU factorization of A. We demonstrate the efficiency and accuracy of PSelInv by presenting several numerical experiments. In particular, we show that PSelInv can run efficiently on more than 4,000 cores for a modestly sized matrix. We also show how PSelInv can be used to accelerate large-scale electronic structure calculations.

97 MATHEMATICS AND COMPUTING↗

PSelInv – A distributed memory parallel algorithm for selected inversion: The non-symmetric case

This paper generalizes the parallel selected inversion algorithm called PSelInv to sparse non-symmetric matrices. We assume a general sparse matrix A has been decomposed as $PAQ$ = $LU$ on a distributed memory parallel machine, where L, U are lower and upper triangular matrices, and P, Q are permutation matrices, respectively. The PSelInvmethod computes selected elements of $A$ -1 . The selection is confined by the sparsity pattern of the matrix A T . Here our algorithm does not assume any symmetry properties of A, and our parallel implementation is memory efficient, in the sense that the computed elements of $A$ -T overwrites the sparse matrix $L + U$ in situ . PSelInv involves a large number of collective data communication activities within different processor groups of various sizes. In order to minimize idle time and improve load balancing, tree-based asynchronous communication is used to coordinate all such collective communication. Numerical results demonstrate that PSelInv can scale efficiently to 6,400 cores for a variety of matrices.

97 MATHEMATICS AND COMPUTING↗

Robust determination of the chemical potential in the pole expansion and selected inversion method for solving Kohn-Sham density functional theory

Fermi operator expansion (FOE) methods are powerful alternatives to diagonalization type methods for solving Kohn-Sham density functional theory (KSDFT). One example is the pole expansion and selected inversion (PEXSI) method, which approximates the Fermi operator by rational matrix functions and reduces the computational complexity to at most quadratic scaling for solving KSDFT. Unlike diagonalization type methods, the chemical potential often cannot be directly read off from the result of a single step of evaluation of the Fermi operator. Hence multiple evaluations are needed to be sequentially performed to compute the chemical potential to ensure the correct number of electrons within a given tolerance. This hinders the performance of FOE methods in practice. In this paper, we develop an efficient and robust strategy to determine the chemical potential in the context of the PEXSI method. The main idea of the new method is not to find the exact chemical potential at each self-consistent-field (SCF) iteration but to dynamically and rigorously update the upper and lower bounds for the true chemical potential, so that the chemical potential reaches its convergence along the SCF iteration. Instead of evaluating the Fermi operator for multiple times sequentially, our method uses a two-level strategy that evaluates the Fermi operators in parallel. In the regime of full parallelization, the wall clock time of each SCF iteration is always close to the time for one single evaluation of the Fermi operator, even when the initial guess is far away from the converged solution. We demonstrate the effectiveness of the new method using examples with metallic and insulating characters, as well as results from ab initio molecular dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Adaptive Pore Opening to Form Tailored Adsorption Sites in a Cooperatively Flexible Framework Enables Record Inverse Propane/Propylene Separation

A proposed low-energy alternative to the separation of alkanes from alkenes by energy-intensive cryogenic distillation is separation by porous adsorbents. Unfortunately, most adsorbents preferentially take up the desired, high-value major component alkene, requiring frequent regeneration. Adsorbents with inverse selectivity for the minor component alkane would enable the direct production of purified, reagent-grade alkene, greatly reducing global energy consumption. However, such materials are exceedingly rare, especially for propane/propylene separation. Here, we report that through adaptive and spontaneous pore size and shape adaptation to optimize an ensemble of weak noncovalent interactions, the structurally responsive metal-organic framework CdIF-13 (sod-Cd(benzimidazolate) 2 ) exhibits inverse selectivity for propane over propylene with record-setting separation performance under industrially relevant temperature, pressure, and mixture conditions. Powder synchrotron X-ray diffraction measurements combined with first-principles calculations yield atomic-scale insight and reveal the induced fit mechanism of adsorbate-specific pore adaptation and ensemble interactions between ligands and adsorbates. Dynamic column breakthrough measurements confirm that CdIF-13 displays selectivity under mixed-component conditions of varying ratios, with a record measured selectivity factor of α ≈ 3 at 95:5 propylene:propane at 298 K and 1 bar. When sequenced with a low-cost rigid adsorbent, we demonstrated the direct purification of propylene under ambient conditions. In conclusion, this combined atomic-level structural characterization and performance testing firmly establishes how cooperatively flexible materials can be capable of unprecedented separation factors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structural basis for receptor selectivity and inverse agonism in S1P5 receptors

This dataset contains a subset of raw diffraction images for the serial femtosecond crystallography data of S1P5-ONO-5430608 complex. The diffraction data contain a strong systematic background as described in Lyapina et al. No additional background subtraction was applied to the data. In order to avoid its effects, the background subtraction implemented in CrystFEL, was applied (--rings-grad subtraction). The final stream, supplemented with this dataset, was obtained with: indexamajig 0.8.0 Respectively, the hkl file was obtained with: partialator 0.8.0

Human sphingosine-1-phosphate receptor 5 (S1P5) in↗

Numerical solution of large scale Hartree–Fock–Bogoliubov equations

The Hartree–Fock–Bogoliubov (HFB) theory is the starting point for treating superconducting systems. However, the computational cost for solving large scale HFB equations can be much larger than that of the Hartree–Fock equations, particularly when the Hamiltonian matrix is sparse, and the number of electrons N is relatively small compared to the matrix size N b . We first provide a concise and relatively self-contained review of the HFB theory for general finite sized quantum systems, with special focus on the treatment of spin symmetries from a linear algebra perspective. We then demonstrate that the pole expansion and selected inversion (PEXSI) method can be particularly well suited for solving large scale HFB equations. For a Hubbard-type Hamiltonian, the cost of PEXSI is at most $\mathcal{O}$(N b 2 ) for both gapped and gapless systems, which can be significantly faster than the standard cubic scaling diagonalization methods. We show that PEXSI can solve a two-dimensional Hubbard-Hofstadter model with N b up to 2.88 × 10 6 , and the wall clock time is less than 100 s using 17 280 CPU cores. Finally, this enables the simulation of physical systems under experimentally realizable magnetic fields, which cannot be otherwise simulated with smaller systems.

97 MATHEMATICS AND COMPUTING↗

Conformational Selection as the Mechanism of Guest Binding in a Flexible Supramolecular Host

This work offers a detailed mechanistic investigation of host-guest encapsulation behavior in a new enzyme-mimetic metal-ligand host and provides the first observation of a conformational selection mechanism (as opposed to induced fit) in a supramolecular system. The Ga 4 L 4 host described features a C 3 -symmetric ligand motif with meta-substituted phenyl spacers, which enables the host to initially self-assemble into an S 4 -symmetric structure and then subsequently isomerize to a T -symmetric tetrahedron for better accommodation of a sufficiently large guest. Selective inversion recovery 1 H NMR studies provide structural insights into the self-exchange behaviors of the host and the guest individually in this dynamic system. Kinetic analysis of the encapsulation-isomerization event revealed that increasing the concentration of guest inhibits the rate of host-guest relaxation, a key distinguishing feature of conformational selection. A comprehensive study of this simple enzyme mimic provides insight into analogous behavior in biophysics and enzymology and aims to inform the design of efficient self-assembled microenvironment catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A molecular switch regulating transcriptional repression and activation of PPARγ

Nuclear receptor (NR) transcription factors use a conserved activation function-2 (AF-2) helix 12 mechanism for agonist-induced coactivator interaction and NR transcriptional activation. In contrast, ligand-induced corepressor-dependent NR repression appears to occur through structurally diverse mechanisms. We report two crystal structures of peroxisome proliferator-activated receptor gamma (PPARγ) in an inverse agonist/corepressor-bound transcriptionally repressive conformation. Helix 12 is displaced from the solvent-exposed active conformation and occupies the orthosteric ligand-binding pocket enabled by a conformational change that doubles the pocket volume. Paramagnetic relaxation enhancement (PRE) NMR and chemical crosslinking mass spectrometry confirm the repressive helix 12 conformation. PRE NMR also defines the mechanism of action of the corepressor-selective inverse agonist T0070907, and reveals that apo-helix 12 exchanges between transcriptionally active and repressive conformations—supporting a fundamental hypothesis in the NR field that helix 12 exchanges between transcriptionally active and repressive conformations.

59 BASIC BIOLOGICAL SCIENCES↗

Elucidation of puzzling questions regarding the CrO x /Al 2 O 3 catalyst I. X-ray absorption spectroscopy aided identification of the nature of the chromium oxide species in the CrO x /Al 2 O 3 dehydrogenation catalyst system

The CATOFIN© process produces propylene from propane using a CrO x /Al 2 O 3 catalyst. Despite its use, there are still uncertainties surrounding the nature of Cr species formed during synthesis, regeneration, and aging. It has been observed that ppm levels of Cr 2 O 3 can be oxidized to Cr(VI) in air; however, when Cr 2 O 3 is in direct contact with Al 2 O 3 , Cr(VI) formation increased more than 100-fold. Here, to confirm physical contact between Cr 2 O 3 and Al 2 O 3 results in increased Cr(VI) formation, Cr 2 O 3 and Cr 2 O 3 /Al 2 O 3 samples were synthesized, calcined, and characterized. Catalytic activity measurements confirm propylene selectivity inversely tracks with Cr(VI). Wet chemical titration and X-ray absorption spectroscopy (XAS) quantified and determined the structure of the Cr(VI) species. FEFF X-ray absorption near edge structure simulations determined the potential sensitivity of the XAS to the presence of Al 3+ /Cr 3+ neighboring atoms. XAS correlated to FEFF simulations suggested that Al 3+ x Cr 3+ 2–x Cr 6+ 3 O 12 comprises a significant fraction of the Cr(VI) species formed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

ELSI: A unified software interface for Kohn–Sham electronic structure solvers

Solving the electronic structure from a generalized or standard eigenproblem is often the bottleneck in large scale calculations based on Kohn-Sham density-functional theory. This problem must be addressed by essentially all current electronic structure codes, based on similar matrix expressions, and by high-performance computation. We here present a unified software interface, ELSI, to access different strategies that address the Kohn-Sham eigenvalue problem. Currently supported algorithms include the dense generalized eigensolver library ELPA, the orbital minimization method implemented in libOMM, and the pole expansion and selected inversion (PEXSI) approach with lower computational complexity for semilocal density functionals. The ELSI interface aims to simplify the implementation and optimal use of the different strategies, by offering (a) a unified software framework designed for the electronic structure solvers in Kohn-Sham density-functional theory; (b) reasonable default parameters for a chosen solver; (c) automatic conversion between input and internal working matrix formats, and in the future (d) recommendation of the optimal solver depending on the specific problem. As a result, comparative benchmarks are shown for system sizes up to 11,520 atoms (172,800 basis functions) on distributed memory supercomputing architectures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Impact of petrophysical uncertainty on Bayesian hydrogeophysical inversion and model selection

Quantitative hydrogeophysical studies rely heavily on petrophysical relationships that link geophysical properties to hydrogeological properties and state variables. Coupled inversion studies are frequently based on the questionable assumption that these relationships are perfect (i.e., no scatter). Using synthetic examples and crosshole ground-penetrating radar (GPR) data from the South Oyster Bacterial Transport Site in Virginia, USA, we investigate the impact of spatially-correlated petrophysical uncertainty on inferred posterior porosity and hydraulic conductivity distributions and on Bayes factors used in Bayesian model selection. Our study shows that accounting for petrophysical uncertainty in the inversion (I) decreases bias of the inferred variance of hydrogeological subsurface properties, (II) provides more realistic uncertainty assessment and (III) reduces the overconfidence in the ability of geophysical data to falsify conceptual hydrogeological models.

58 GEOSCIENCES↗

Inverse Bimetallic RuSn Catalyst for Selective Carboxylic Acid Reduction

Inverse bimetallic catalysts (IBCs), synthesized by sequential deposition of noble and oxophilic metals, offer potential reactivity enhancements to various reactions, including the reduction of carboxylic acids for renewable fuels and chemicals. Here, we demonstrate that an IBC comprising RuSn exhibits high selectivity for propionic acid reduction to 1-propanol, while Ru alone results in cracking. On RuSn, X-ray absorption spectroscopy identified Ru o nanoparticles with a near-surface bimetallic Ru o Sn o alloy and small SnO x domains. Corresponding model surfaces were examined with density functional theory to elucidate the observed selectivity difference. Only selective hydrogenation is predicted to be favorable on SnO x /Ru, with the SnO x clusters facilitating C-OH scission and Ru enabling hydrogen activation. Intrinsic barriers along nonselective pathways suggest that the RuSn alloy and SnO x resist cracking. SnO x /Ru hydrogenation activity was supported experimentally by inhibiting hydrogenation with phenylphosphonic acid, differentiating the system from fully alloyed RuSn metallic nanoparticles. Altogether, this work demonstrates a plausible mechanism for selective reduction of carboxylic acids and proposes a roadmap for rational design of IBCs.

09 BIOMASS FUELS↗

Full Waveform Inversion with Multisource Frequency Selection of Marine Streamer Data

The theory and practice of multisource full waveform inversion of marine supergathers are described with a frequency-selection strategy. The key enabling property of frequency selection is that it eliminates the crosstalk among sources, thus overcoming the aperture mismatch of marine multisource inversion. Tests on multisource full waveform inversion of synthetic marine data and Gulf of Mexico data show speedups of 4× and 8×, respectively, compared to conventional full waveform inversion.

58 GEOSCIENCES↗

Young inversion with multiple linked QTLs under selection in a hybrid zone

Fixed chromosomal inversions can reduce gene flow and promote speciation in two ways: by suppressing recombination and by carrying locally favoured alleles at multiple loci. However, it is unknown whether favoured mutations slowly accumulate on older inversions or if young inversions spread because they capture pre-existing adaptive quantitative trait loci (QTLs). By genetic mapping, chromosome painting and genome sequencing, we have identified a major inversion controlling ecologically important traits in Boechera stricta. The inversion arose since the last glaciation and subsequently reached local high frequency in a hybrid speciation zone. Furthermore, the inversion shows signs of positive directional selection. To test whether the inversion could have captured existing, linked QTLs, we crossed standard, collinear haplotypes from the hybrid zone and found multiple linked phenology QTLs within the inversion region. These findings provide the first direct evidence that linked, locally adapted QTLs may be captured by young inversions during incipient speciation.

59 BASIC BIOLOGICAL SCIENCES↗

Inverse chirality-induced spin selectivity effect in chiral assemblies of π -conjugated polymers

Coupling of spin and charge currents to structural chirality in non-magnetic materials, known as chirality-induced spin selectivity, is promising for application in spintronic devices at room temperature. Although the chirality-induced spin selectivity effect has been identified in various chiral materials, its Onsager reciprocal process, the inverse chirality-induced spin selectivity effect, remains unexplored. Here we report the observation of the inverse chirality-induced spin selectivity effect in chiral assemblies of π-conjugated polymers. Using spin-pumping techniques, the inverse chirality-induced spin selectivity effect enables quantification of the magnitude of the longitudinal spin-to-charge conversion driven by chirality-induced spin selectivity in different chiral polymers. By widely tuning conductivities and supramolecular chiral structures via a printing method, we found a very long spin relaxation time of up to several nanoseconds parallel to the chiral axis. Furthermore, our demonstration of the inverse chirality-induced spin selectivity effect suggests possibilities for elucidating the puzzling interplay between spin and chirality, and opens a route for spintronic applications using printable chiral assemblies.

36 MATERIALS SCIENCE↗

Delineating Facies Spatial Distribution by Integrating Ensemble Data Assimilation and Indicator Geostatistics with Level-Set Transformation

The facies-based approach has been widely adopted to delineate an aquifer into distinct geological units with unique distributions of hydraulic, physical, and/or chemical properties. The recent development in ensemble-based data assimilation methods allows both the direct and indirect data to be used to improve facies delineation. A major difficulty in those applications is to honor the spatial continuity and avoid overfitting after data assimilation. We introduce a new facies delineation framework to integrate ensemble data assimilation with traditional transition probability-based geostatistics. A level-set concept is used to parametrize discrete facies indicators and for updating facies shape. During the iterative data assimilation process, we impose spatial continuity by conditioning facies field generation on points selected adaptively based on their sensitivity to observation data. This reconditioning step is a key step to maintain spatial continuity and overcome over-fitting problems in inversion. We selected two examples to evaluate the performance of the new framework in estimating facies-based permeability field. The first example is a two-dimensional synthetic system with transient head data induced by pumping tests used for delineating two facies. The second example is a three-dimensional case with three facies, conceptualized from a field tracer experiment within the Columbia River corridor in Washington State, USA. Both examples demonstrate that the new method can adequately capture the spatial pattern of hydrofacies with reconditioning, which leads to the improved prediction of system behaviors.

Song, Xuehang↗