SEARCH · Engineering Papers

Results for “Density matrix methods”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 631 records · Page 35

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Stimulated X-ray Raman and Absorption Spectroscopy of Iron–Sulfur Dimers

Iron–sulfur complexes play an important role in biological processes such as metabolic electron transport. A detailed understanding of the mechanism of long-range electron transfer requires knowledge of the electronic structure of the complexes, which has traditionally been challenging to obtain, either by theory or by experiment, but the situation has begun to change with advances in quantum chemical methods and intense free electron laser light sources. We compute the spectra for stimulated X-ray Raman spectroscopy (SXRS) and absorption spectroscopy of homovalent and mixed-valence [2Fe–2S] complexes, using the ab initio density matrix renormalization group algorithm. The simulated spectra show clear signatures of the theoretically predicted dense low-lying excited states within the d–d manifold. Finally, the difference in spectral intensity between the absorption-active and Raman-active states provides a potential mechanism to selectively excite states by a proper tuning of the excitation pump, to access the electronic dynamics within this manifold.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2002

Electronic Absorption Spectra of Neutral Perylene (C20H12), Terrylene (C30H16), and Quaterrylene (C40H20) and their Positive and Negative Ions: Ne Matrix-Isolation Spectroscopy and Time Dependent Density Functional Theory Calculations

We present a full experimental and theoretical study of an interesting series of polycyclic aromatic hydrocarbons, the oligorylenes. The absorption spectra of perylene, terrylene and quaterrylene in neutral, cationic and anionic charge states are obtained by matrix-isolation spectroscopy in Ne. The experimental spectra are dominated by a bright state that red shifts with growing molecular size. Excitation energies and state symmetry assignments are supported by calculations using time dependent density functional theory methods. These calculations also provide new insight into the observed trends in oscillator strength and excitation energy for the bright states: the oscillator strength per unit mass of carbon increases along the series.

Halasinski, Thomas M.↗

NASA NTRS2004

Signal Processing Variables for Optimization of Flaw Detection in Composites Using Ultrasonic Guided Wave Scanning

This study analyzes the effect of signal processing variables on the ability of the ultrasonic guided wave scan method at NASA Glenn Research Center to distinguish various flaw conditions in ceramic matrix composites samples. In the ultrasonic guided wave scan method, several time- and frequency-domain parameters are calculated from the ultrasonic guided wave signal at each scan location to form images. The parameters include power spectral density, centroid mean time, total energy (zeroth moment), centroid frequency, and ultrasonic decay rate. A number of signal processing variables are available to the user when calculating these parameters. These signal processing variables include 1) the time portion of the time-domain waveform processed, 2) integration type for the properties requiring integrations, 3) bounded versus unbounded integrations, 4) power spectral density window type, 5) and the number of time segments chosen if using the short-time fourier transform to calculate ultrasonic decay rate. Flaw conditions examined included delamination, cracking, and density variation.

Roth, Don J.↗

NASA NTRS2005

Multi-Scale Modeling of Cross-Linked Nanotube Materials

The effect of cross-linking single-walled carbon nanotubes on the Young's modulus of a nanotube-reinforced composite is modeled with a multi-scale method. The Young's modulus is predicted as a function of nanotube volume fraction and cross-link density. In this method, the constitutive properties of molecular representative volume elements are determined using molecular dynamics simulation and equivalent-continuum modeling. The Young's modulus is subsequently calculated for cross-linked nanotubes in a matrix which consists of the unreacted cross-linking agent. Two different cross-linking agents are used in this study, one that is short and rigid (Molecule A), and one that is long and flexible (Molecule B). Direct comparisons between the predicted elastic constants are made for the models in which the nanotubes are either covalently bonded or not chemically bonded to the cross-linking agent. At a nanotube volume fraction of 10%, the Young's modulus of Material A is not affected by nanotube crosslinking, while the Young's modulus of Material B is reduced by 64% when the nanotubes are cross-linked relative to the non-cross-linked material with the same matrix.

Frankland, S. J. V.↗

DOE OSTI2017

Efficient block preconditioned eigensolvers for linear response time-dependent density functional theory

We present two efficient iterative algorithms for solving the linear response eigen- value problem arising from the time dependent density functional theory. Although the matrix to be diagonalized is nonsymmetric, it has a special structure that can be exploited to save both memory and floating point operations. In particular, the nonsymmetric eigenvalue problem can be transformed into a product eigenvalue problem that is self-adjoint with respect to a K-inner product. This product eigenvalue problem can be solved efficiently by a modified Davidson algorithm and a modified locally optimal block preconditioned conjugate gradient (LOBPCG) algorithm that make use of the K-inner product. The solution of the product eigenvalue problem yields one component of the eigenvector associated with the original eigenvalue problem. However, the other component of the eigenvector can be easily recovered in a postprocessing procedure. Therefore, the algorithms we present here are more efficient than existing algorithms that try to approximate both components of the eigenvectors simultaneously. The efficiency of the new algorithms is demonstrated by numerical examples.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2017

Monte Carlo/CFD Coupling for Accurate Modeling of the Delayed Neutron Precursors and Compressibility Effects in Molten Salt Reactors

Multiphysics reactor analysis has become one of the leading activities of most nuclear engineering research groups. Coupled neutronics and thermal/hydraulics (T/H) simulations have traditionally been performed by adopting deterministic transport tools and reactor system T/H codes. The main reason for avoiding higher fidelity numerical methods like continuous-energy Monte Carlo analysis has been the high computational cost of coupled, full-scale simulations. Nowadays, high fidelity multiphysics simulations represent a hot research topic, with several applications involving advanced nuclear systems. Liquid fuel Molten Salt Reactors (MSRs) feature peculiar physical phenomena that are not present in solid fuel reactors. Some of these phenomena are not correctly resolved by legacy reactor physics tools, and require specific treatments. This work focuses in particular on the analysis of the effect of the delayed neutron precursors drift and the impact of liquid fluid thermal expansion and compressibility in postulated criticality excursion transients in fast-spectrum MSRs. A new modeling tool based on a Monte Carlo approach for neutronics (Serpent code) coupled to a Computational Fluid Dynamic (CFD) analysis (OpenFOAM) was developed to obtain accurate modeling and simulation of these phenomena. The development of the Serpent Monte Carlo code [1] at VTT in the past few years has made this code particularly attractive for performing multiphysics studies. The availability of a flexible multiphysics interface, on-the-fly Doppler broadening and density-based collision rejection, and Delta-tracking, proved of paramount importance for the computationally efficient performance of coupled Monte Carlo/CFD calculations [2]. OpenFOAM is a C++ multiphysics toolkit featured by automatic matrix construction and solution capabilities for scalar and vector partial differential equations (PDEs), adopting several possible differencing and interpolation schemes, based on the Finite-Volume method. The ease of implementing the constitutive equations of different physical phenomena makes OpenFOAM particularly suitable for the modeling of complex coupled phenomena. In the past, several research activities proposed the adoption of both Serpent and OpenFOAM for the modeling and simulation of MSRs and the associated fuel cycles. In this work, these tools have been internally coupled for the steady-state and transient analysis of fast-spectrum systems, accounting for the effect of delayed neutron precursors transport and the salt thermal expansion and compressibility.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2014

Electron correlation in solids via density embedding theory

Density matrix embedding theory [G. Knizia and G. K.-L. Chan, Phys. Rev. Lett. 109, 186404 (2012)] and density embedding theory [I. W. Bulik, G. E. Scuseria, and J. Dukelsky, Phys. Rev. B 89, 035140 (2014)] have recently been introduced for model lattice Hamiltonians and molecular systems. In the present work, the formalism is extended to the ab initio description of infinite systems. An appropriate definition of the impurity Hamiltonian for such systems is presented and demonstrated in cases of 1, 2, and 3 dimensions, using coupled cluster theory as the impurity solver. Additionally, we discuss the challenges related to disentanglement of fragment and bath states. The current approach yields results comparable to coupled cluster calculations of infinite systems even when using a single unit cell as the fragment. The theory is formulated in the basis of Wannier functions but it does not require separate localization of unoccupied bands. The embedding scheme presented here is a promising way of employing highly accurate electronic structure methods for extended systems at a fraction of their original computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Polymer-fiber-reinforced polymers with enhanced interfacial bonding between polypropylene fiber and polyethylene matrix

Self-reinforced composites (SRCs) consist of reinforcing fibers and a base matrix made of the same thermoplastic polymer, offering lightweight, recyclability, and sustainability benefits. However, limited research exists on composites where the reinforcing thermoplastic polymer fibers differ from the base thermoplastic matrix. Here, this study focuses on investigating the mechanical behavior of such composites and exploring different surface modification methods to enhance the fiber/matrix interfacial bonding using polypropylene fibers and a polyethylene matrix as an example. It is shown that surface treatment with a commercial adhesion promoter containing n-butyl acetate significantly improves the interfacial shear strength between polypropylene fibers and the polyethylene matrix, increasing it by 145% compared to other methods investigated. Additionally, increasing the length of the embedded polymer fiber in the matrix leads to a notable increase in specific interfacial energy. Consequently, the thermoplastic polymer-fiber-reinforced polymers (PFRPs) using surface-treated woven polypropylene fabrics and a polyethylene matrix exhibit a 20% higher tensile strength and a 65% higher toughness compared to non-treated PFRPs. This study also shows that specific mechanical properties (normalized by the composite density) of the investigated woven PFRPs are similar to those of non-treated SRCs under uni-axial tension. Particularly, their ductility outperforms carbon-/glass-/aramid-fiber-reinforced polymers by at least 6 times at a same fiber volume fraction. The investigation of such composites and the exploration of surface modification methods present important progress in the field of thermoplastic PFRPs, which serve as a solution for addressing concerns related to recyclability and sustainability.

Fiber pull-out↗

DOE OSTI2019

The HANDE-QMC Project: Open-Source Stochastic Quantum Chemistry from the Ground State Up

Building on the success of Quantum Monte Carlo techniques such as diffusion Monte Carlo, alternative stochastic approaches to solve electronic structure problems have emerged over the past decade. The full configuration interaction quantum Monte Carlo (FCIQMC) method allows one to systematically approach the exact solution of such problems, for cases where very high accuracy is desired. The introduction of FCIQMC has subsequently led to the development of coupled cluster Monte Carlo (CCMC) and density matrix quantum Monte Carlo (DMQMC), allowing stochastic sampling of the coupled cluster wave function and the exact thermal density matrix, respectively. In this Article, we describe the HANDE-QMC code, an open-source implementation of FCIQMC, CCMC and DMQMC, including initiator and semistochastic adaptations. We describe our code and demonstrate its use on three example systems; a molecule (nitric oxide), a model solid (the uniform electron gas), and a real solid (diamond). An illustrative tutorial is also included.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

A novel assay to measure low-density lipoproteins binding to proteoglycans

Background: The binding of low-density lipoprotein (LDL) to proteoglycans (PGs) in the extracellular matrix (ECM) of the arterial intima is a key initial step in the development of atherosclerosis. Although many techniques have been developed to assess this binding, most of the methods are labor-intensive and technically challenging to standardize across research laboratories. Thus, sensitive, and reproducible assay to detect LDL binding to PGs is needed to screen clinical populations for atherosclerosis risk. Objectives: The aim of this study was to develop a quantitative, and reproducible assay to evaluate the affinity of LDL towards PGs and to replicate previously published results on LDL-PG binding. Methods: Immunofluorescence microscopy was performed to visualize the binding of LDL to PGs using mouse vascular smooth muscle (MOVAS) cells. An in-cell ELISA (ICE) was also developed and optimized to quantitatively measure LDL-PG binding using fixed MOVAS cells cultured in a 96-well format. Results: We used the ICE assay to show that, despite equal APOB concentrations, LDL isolated from adults with cardiovascular disease bound to PG to a greater extent than LDL isolated from adults without cardiovascular disease (p<0.05). Conclusion: We have developed an LDL-PG binding assay that is capable of detecting differences in PG binding affinities despite equal APOB concentrations. Future work will focus on candidate apolipoproteins that enhance or diminish this interaction.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1990

Slow plastic deformation of extruded NiAl-10TiB2 particulate composites at 1200 and 1300 K

A dispersion of 1-micron TiB2 particles in the B2 crystal structure NiAl intermetallic can effectively increase its elevated temperature strength, in association with increasing deformation resistance with TiB2 volume fraction. Attention is presently given to alternative densification methods, which may increase the initial as-fabricated dislocation density and lead to enhanced elevated-temperature strength. The 'XD' extrusion method was used to produce NiAl with 10 vol pct TiB2. Although apparent extrusion defects were occasionally found, neither grain-boundary cracking nor particle-matrix separation occurred.

Whittenberger, J. D.↗

DOE OSTI2023

A fast, dense Chebyshev solver for electronic structure on GPUs

Matrix diagonalization is almost always involved in computing the density matrix needed in quantum chemistry calculations. In the case of modest matrix sizes (≲4000), performance of traditional dense diagonalization algorithms on modern GPUs is underwhelming compared to the peak performance of these devices. This motivates the exploration of alternative algorithms better suited to these types of architectures. We newly derive, and present in detail, an existing Chebyshev expansion algorithm whose number of required matrix multiplications scales with the square root of the number of terms in the expansion. Focusing on dense matrices of modest size, our implementation on GPUs results in large speed ups when compared to diagonalization. Additionally, we improve upon this existing method by capitalizing on the inherent task parallelism and concurrency in the algorithm. Furthermore, this improvement is implemented on GPUs by using CUDA and HIP streams via the MAGMA library and leads to a significant speed up over the serial-only approach for smaller (≲1000) matrix sizes. Finally, we apply our technique to a model system with a high density of states around the Fermi level, which typically presents significant challenges.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2017

Towards optimal extraction of cosmological information from nonlinear data

One of the main unsolved problems of cosmology is how to maximize the extraction of information from nonlinear data. If the data are nonlinear the usual approach is to employ a sequence of statistics (N-point statistics, counting statistics of clusters, density peaks or voids etc.), along with the corresponding covariance matrices. However, this approach is computationally prohibitive and has not been shown to be exhaustive in terms of information content. Here we instead develop a hierarchical Bayesian approach, expanding the likelihood around the maximum posterior of linear modes, which we solve for using optimization methods. By integrating out the modes using perturbative expansion of the likelihood we construct an initial power spectrum estimator, which for a fixed forward model contains all the cosmological information if the initial modes are gaussian distributed. We develop a method to construct the window and covariance matrix such that the estimator is explicitly unbiased and nearly optimal. We then generalize the method to include the forward model parameters, including cosmological and nuisance parameters, and primordial non-gaussianity. We apply the method in the simplified context of nonlinear structure formation, using either simplified 2-LPT dynamics or N-body simulations as the nonlinear mapping between linear and nonlinear density, and 2-LPT dynamics in the optimization steps used to reconstruct the initial density modes. We demonstrate that the method gives an unbiased estimator of the initial power spectrum, providing among other a near optimal reconstruction of linear baryonic acoustic oscillations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

Optical Properties of New Type of Superconductor-Semiconductor Metamaterial Photonic Crystals

Based on the fundamentals of the characteristic matrix method, we theoretically investigate the optical properties of onedimensional superconductor metamaterial photonic crystals. The photonic crystals are composed of a superconductor layer and two semiconductor metamaterial layers of In{sub 0.53}Ga{sub 0.47}As with different doping densities. The numerical results show negative values in permittivity of the metamaterial layer along a broad band of the incident radiation. The negative values show a significant effect on the optical properties of the present structure. Moreover, the optical properties of our design can be controlled by different parameters such as thicknesses of the layers, the operating temperature and the doping density. Our results may be suitable for different applications in the optical community.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2015

The Coulomb, exchange, and correlation components of the electron-electron repulsion in harmonium atoms

Highly accurate Coulomb, exchange, and correlation components of the electron-electron repulsion energies of the three-electron harmonium atoms in the {sup 2}P{sub −} and {sup 4}P{sub +} states are obtained for 19 values of the confinement strength ω ranging from 10{sup −3} to 10{sup 3}. The computed data are consistent with their ω → 0 and ω → ∞ asymptotics that are given by closed-form algebraic expressions. Robust approximants that accurately reproduce the actual values of the energy components while strictly conforming to these limits are constructed, opening an avenue to stringent tests capable of predicting the performance of electronic structure methods for systems with varying extents of the dynamical and nondynamical electron correlation. The values of the correlation components, paired with the computed 1-matrices are expected to be particularly useful in the context of benchmarking of approximate density matrix functionals.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS2007

A Simple Stochastic Model for Generating Broken Cloud Optical Depth and Top Height Fields

A simple and fast algorithm for generating two correlated stochastic twodimensional (2D) cloud fields is described. The algorithm is illustrated with two broken cumulus cloud fields: cloud optical depth and cloud top height retrieved from Moderate Resolution Imaging Spectrometer (MODIS). Only two 2D fields are required as an input. The algorithm output is statistical realizations of these two fields with approximately the same correlation and joint distribution functions as the original ones. The major assumption of the algorithm is statistical isotropy of the fields. In contrast to fractals and the Fourier filtering methods frequently used for stochastic cloud modeling, the proposed method is based on spectral models of homogeneous random fields. For keeping the same probability density function as the (first) original field, the method of inverse distribution function is used. When the spatial distribution of the first field has been generated, a realization of the correlated second field is simulated using a conditional distribution matrix. This paper is served as a theoretical justification to the publicly available software that has been recently released by the authors and can be freely downloaded from http://i3rc.gsfc.nasa.gov/Public codes clouds.htm. Though 2D rather than full 3D, stochastic realizations of two correlated cloud fields that mimic statistics of given fields have proved to be very useful to study 3D radiative transfer features of broken cumulus clouds for better understanding of shortwave radiation and interpretation of the remote sensing retrievals.

Prigarin, Sergei M.↗

DOE OSTI2019

First-principles calculations of solute transport in zirconium: Vacancy-mediated diffusion with metastable states and interstitial diffusion

Zirconium alloys are the most widely used nuclear fuel cladding materials for light water power reactors where irradiation damage causes solute redistribution, leading to degradation of alloy properties such as corrosion resistance. Designing radiation-tolerant zirconium alloys requires a thorough understanding of the atomic-scale transport behavior of the alloying elements in Zr. In this work, we perform density function theory calculations to investigate the diffusion of Sn, Cr, Fe, Be, Al, and Ni in the hexagonal close-packed (HCP) Zr matrix. We develop a methodology to accurately model the metastable vacancy states along the basal migration path, known to occur in group IV metals. We compute the vacancy-mediated solute diffusion coefficients and drag ratios using the kinetic Monte Carlo method and an analytic Green's function method—the agreement between the two validates our methodology. The computed diffusion coefficients of Sn and Al show good agreement with the experimental data and we expect these solutes to diffuse via the vacancy-mediated mechanism. We use a Green's function approach, parameterized with data from density functional theory calculations, to compute the interstitial diffusion coefficients of Cr, Fe, Be, and Ni in the HCP Zr lattice. The computed diffusion coefficients of Cr, Ni, and Be agree with the experimental measurements within one order of magnitude, while those of Fe are within two orders of magnitude of the experimental measurements. The drag ratios for Cr, Fe, Be, and Ni are positive up to 1235 K, which suggests that nonequilibrium vacancy fluxes could drag these solutes toward sinks such as dislocation loops and grain boundaries. We also compute the transport coefficients without including the metastable states, and using the eight- and thirteen-frequency model. Our results show significant differences in drag ratio for the eight- and thirteen-frequency model predictions compared with the Green's function methodology, but smaller errors in the solute diffusivity. Combining interstitial and vacancy-mediated diffusivities, we predict the unusual result that increased vacancy concentration slows down solute diffusivity, while a sufficiently high vacancy concentration can change the dominant mechanism to an accelerated vacancy-mediated diffusion.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Quantum simulation of molecules without fermionic encoding of the wave function

Abstract Molecular simulations generally require fermionic encoding in which fermion statistics are encoded into the qubit representation of the wave function. Recent calculations suggest that fermionic encoding of the wave function can be bypassed, leading to more efficient quantum computations. Here we show that the two-electron reduced density matrix (2-RDM) can be expressed as a unique functional of the unencoded N -qubit-particle wave function without approximation, and hence, the energy can be expressed as a functional of the 2-RDM without fermionic encoding of the wave function. In contrast to current hardware-efficient methods, the derived functional has a unique, one-to-one (and onto) mapping between the qubit-particle wave functions and 2-RDMs, which avoids the over-parametrization that can lead to optimization difficulties such as barren plateaus. An application to computing the ground-state energy and 2-RDM of H 4 is presented.

74 ATOMIC AND MOLECULAR PHYSICS↗