Engineering Papers⌕ Search

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.

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

At least 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

The principal component analysis method used with polynomial Chaos expansion to propagate uncertainties through critical transport problems

A critical neutron transport problem with random material properties is introduced. The total cross section and the average neutron multiplicity are assumed to be uncertain, characterized by the mean and variance with a log-normal distribution. The average neutron multiplicity and the total cross section are assumed to be uncorrected and the material properties for differing materials are also assumed to be uncorrected. The principal component analysis method is used to decompose the covariance matrix into eigenvalues and eigenvectors and then 'realizations' of the material properties can be computed. A simple Monte Carlo brute force sampling of the decomposed covariance matrix is employed to obtain a benchmark result for each test problem. In order to save computational time and to characterize the moments and probability density function of the multiplication factor the polynomial chaos expansion method is employed along with the stochastic collocation method. A Gauss-Hermite quadrature set is convolved into a multidimensional tensor product quadrature set and is successfully used to compute the polynomial chaos expansion coefficients of the multiplication factor. Finally, for a particular critical fuel pin assembly the appropriate number of random variables and polynomial expansion order are investigated. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2013

Study of second phase in bioabsorbable magnesium alloys: Phase stability evaluation via Dmol{sup 3} calculation

Thermodynamical stabilities of four conventional second phases as well as magnesium matrix in bioabsorbable magnesium alloys were investigated theoretically via computer calculation method. Model of individual phase and systems including phase and four water molecular (phase-4H{sub 2}O) were established to simulate the in vitro and in vivo environment. Local orbital density functional theory approach was applied to calculate the total energy for the individual phase and phase-4H{sub 2}O system. The results demonstrated that all the second phases possessed higher phase stability compared with magnesium matrix, but the phase stability was quite different for different types of second phases or second phase-4H{sub 2}O systems. Furthermore, a schematic process of inflammation reaction caused by magnesium alloy implants was proposed for the further evaluation on biocompatibility of different second phases.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Selfsimilar and fractal analysis of n-type delta-doped quasiregular GaAs quantum wells

We study the electronic structure of n-type delta-doped quantum wells in GaAs in which the multiple well system is built according to the Fibonacci sequence. The building blocks A and B correspond to delta-doped wells with impurities densities n{sub 2DA} and n{sub 2DB}, and the same well width. The Thomas-Fermi approximation, the semi-empirical sp{sub 3}s* tight-binding model including spin, the Surface Green Function Matching method and the Transfer Matrix approach were implemented to obtain the confining potential, the electronic structure and the selfsimilarity of the spectrum. The fragmentation of the electronic spectra is observed whenever the building blocks A and B interact and it increases as the difference of impurities density between A and B increases as well. The wave function of the first sate of the fragmented bands presents critical characteristics, this is, it is not a localized state nor a extended one as well as it has selfsimilar features. So, the quasiregular characteristics are preserved irrespective of the complexity of the system and can affect the performance of devices based on these structures.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2014

First-order nonadiabatic coupling matrix elements between excited states: A Lagrangian formulation at the CIS, RPA, TD-HF, and TD-DFT levels

Analytic expressions for the first-order nonadiabatic coupling matrix elements between electronically excited states are first formulated exactly via both time-independent equation of motion and time-dependent response theory, and are then approximated at the configuration interaction singles, particle-hole/particle-particle random phase approximation, and time-dependent density functional theory/Hartree-Fock levels of theory. Note that, to get the Pulay terms arising from the derivatives of basis functions, the standard response theory designed for electronic perturbations has to be extended to nuclear derivatives. The results are further recast into a Lagrangian form that is similar to that for excited-state energy gradients and allows to use atomic orbital based direct algorithms for large molecules.

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

Effects of refractive index mismatch in optical CT imaging of polymer gel dosimeters

Purpose: Proposing an image reconstruction technique, algebraic reconstruction technique-refraction correction (ART-rc). The proposed method takes care of refractive index mismatches present in gel dosimeter scanner at the boundary, and also corrects for the interior ray refraction. Polymer gel dosimeters with high dose regions have higher refractive index and optical density compared to the background medium, these changes in refractive index at high dose results in interior ray bending. Methods: The inclusion of the effects of refraction is an important step in reconstruction of optical density in gel dosimeters. The proposed ray tracing algorithm models the interior multiple refraction at the inhomogeneities. Jacob’s ray tracing algorithm has been modified to calculate the pathlengths of the ray that traverses through the higher dose regions. The algorithm computes the length of the ray in each pixel along its path and is used as the weight matrix. Algebraic reconstruction technique and pixel based reconstruction algorithms are used for solving the reconstruction problem. The proposed method is tested with numerical phantoms for various noise levels. The experimental dosimetric results are also presented. Results: The results show that the proposed scheme ART-rc is able to reconstruct optical density inside the dosimeter better than the results obtained using filtered backprojection and conventional algebraic reconstruction approaches. The quantitative improvement using ART-rc is evaluated using gamma-index. The refraction errors due to regions of different refractive indices are discussed. The effects of modeling of interior refraction in the dose region are presented. Conclusions: The errors propagated due to multiple refraction effects have been modeled and the improvements in reconstruction using proposed model is presented. The refractive index of the dosimeter has a mismatch with the surrounding medium (for dry air or water scanning). The algorithm reconstructs the dose profiles by estimating refractive indices of multiple inhomogeneities having different refractive indices and optical densities embedded in the dosimeter. This is achieved by tracking the path of the ray that traverses through the dosimeter. Extensive simulation studies have been carried out and results are found to be matching that of experimental results.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

Multichannel filter application of a magnetized cold plasma defect in periodic structure of ZnS/TiO{sub 2} materials

In this article, transmission spectra of multilayer structure of zinc sulfide (ZnS) and titanium dioxide (TiO{sub 2}) with defect of magnetized cold plasma (MCP) have been calculated theoretically by using characteristics transfer matrix method. The transmission spectra versus frequency have been analyzed for symmetric and asymmetric structure. Due to have a multiple band gaps in the asymmetric structure, we have taken the asymmetric structure for analysis with the parameters of MCP and the thicknesses of the dielectric materials. It is well known that magnetized cold plasma is an abnormal material because the permittivity changes with applied external magnetic field, and shows tunable photonic band gap for permittivity of the MCP. Now, the transmission spectra of the asymmetric structure versus frequency by inserting one or two magnetized cold plasma layer were analyzed with the variation of incident angles, the parameters of MCP, and the thicknesses of ZnS and TiO{sub 2} material. The analyzed transmission spectra of asymmetric structure with the variation of electron density, thicknesses of dielectric materials were obtained better results with defect of two plasma layers as compare to the defect of one layer of plasma. These obtained results reveal that the defect of two magnetized cold plasma layers of one dimensional asymmetric structure of ZnS, TiO{sub 2} may be used for the tunable multichannel filter at microwave region.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Formation of Composite SiC/SiC Joints by Embedded Wire Chemical Vapor Deposition

The joining of ceramic monoliths or composites to date has primarily been limited to the formation of brittle monolithic joints using heterogeneous (dissimilar) materials, similar to brazing in metals. The development of a damage-tolerant joint layer by SiC fiber reinforcements is demonstrated here. Tube workpieces made of SiC fiber-SiC matrix composite are joined using a nonwoven SiC fiber mat densified by embedded wire chemical vapor deposition (EWCVD), creating a fiber-reinforced weld-like joint by homogeneous joining. EWCVD uses a localized heating method to target deposition and growth to the joint region specifically, while minimizing thermal damage to the surrounding composite tube material. X-ray computed tomography (XCT) is used to nondestructively characterize as-made joints for relative density, adhesion, and composition. In situ XCT analysis during mechanical testing revealed crack deflections in the bonding layer, which indicates a toughening mechanism typical of ceramic matrix composite phase. Gas permeation testing of these proof-of-concept composite joints identified relatively high leak rates in comparison to fully coated SiC/SiC composite tube workpieces. In conclusion, the novelty of the composite joining method and current technology challenges, including gas permeability, are discussed in comparison with traditional ceramic joints and materials.

SiC↗
DOE OSTI2019

O(N) Quantum Molecular Dynamics (QMD) Algorithms, Code Requirements, Release Notes, and Library Deployment

The Exascale Computing Project (ECP) Co-design Center for Particle Applications (CoPA) is addressing the challenges for particle applications to run on upcoming exascale computing architectures. This includes short- and long-range molecular dynamics (MD), particle-in-cell (PIC), and quantum molecular dynamics (QMD). This document addresses $\mathcal{O}(N)$ QMD. The goal of this CoPA team is to develop $\mathcal{O}(N)$ linear scaling methods for QMD optimized for new and next generation architectures. This will provide the capability to simulate larger systems and run for longer timeframes than is possible today. Focus is initially on the density matrix construction. This requires efficient matrix formats and linear algebra operations for distributed architectures. A significant reduction in the wallclock time for an MD timestep is necessary to achieve practical large scale QMD. Equally challenging is to develop scalable distributed implementations. This requires the exploration of efficient data formats (dense, sparse, etc.), data decomposition and distribution across compute nodes using block- and graph-based methods. Traditional programming models such as MPI/OpenMP, hybrid MPI/OpenMP/Accelerator, as well as asynchronous task-based models will be explored to achieve maximum parallelism. Parallelism at the level of a matrix operation, solver, and/or application needs to be explored to gain the most benefit on new and emerging architectures. This report serves as a “living document” of the additions and enhancements throughout the ECP CoPA project for $\mathcal{O}(N)$ QMD. The following sections address the applications we are working with, the quantum chemistry libraries (BML and PROGRESS), and our co-design exploration through the use of a proxy application.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Towards exact finite temperature electronic structure in solids and molecules (Final Technical Report)

This report describes the University of Iowa portion of a project that is now continuing at Michigan State University. We are developing novel methods, algorithms, and software to enable simulations of molecules and materials at high temperature. This is a key challenge in chemistry and materials science. By refining an approach called Density Matrix Quantum Monte Carlo (DMQMC), we developed faster and more accurate ways to conduct these simulations. These advances will help us understand how temperature affects the behavior of electrons, chemical bonds, and phase transitions in solids and molecules. These breakthroughs are especially important for applications where light and heat drive chemical reactions, superconductivity, and materials used in energy and sensing. In addition, this project involved the development of the open-source HANDE-QMC software package, supporting the broader community in benchmarking and developing finite-temperature electronic structure methods.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Controlled expansion for pairing in a polarized band with strong repulsion

Can strong repulsive interactions be shown to give rise to pairing in a controlled way? We find that for a single flavor polarized band, there is a small expansion parameter in the low density limit, once the Bloch wavefunction form factor is taken into account. A perturbative expansion is possible, even if the interaction is much stronger than the Fermi energy ε F . As a matter of principle, our work shows analytically how strong pairing can emerge from strong repulsion. Here, we illustrate our method with two examples: a 2D Dirac model and a 1D tight binding model with two orbitals. In the latter case, using density matrix renormalization group, we show that the analytical theory indeed guided us to discover the parameter regime where p-wave pairing with order-1 strength is dominant.

Ferromagnetic superconductors↗
DOE OSTI2012

A Practical Framework for Sharing and Rendering Real-World Bidirectional Scattering Distribution Functions

The utilization of real-world materials has been hindered by a lack of standards for sharing and interpreting measured data. This paper presents an XML representation and an Open Source C library to support bidirectional scattering distribution functions (BSDFs) in data-driven lighting simulation and rendering applications.The library provides for the efficient representation, query, and Monte Carlo sampling of arbitrary BSDFs in amodel-free framework. Currently, we support two BSDF data representations: one using a fixed subdivision of thehemisphere, and one with adaptive density. The fixed type has advantages for certain matrix operations, while theadaptive type can more accurately represent highly peaked data. We discuss advanced methods for data-drivenBSDF rendering for both types, including the proxy of detailed geometry to enhance appearance and accuracy.We also present an advanced interpolation method to reduce measured data into these standard representations.We end with our plan for future extensions and sharing of BSDF data.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

General Bayesian Inference over the Stiefel Manifold via the Givens Transform

We introduce here an approach based on the Givens representation for posterior inference in statistical models with orthogonal matrix parameters, such as factor models and probabilistic principal component analysis (PPCA). We show how the Givens representation can be used to develop practical methods for transforming densities over the Stiefel manifold into densities over subsets of Euclidean space. We demonstrate how to deal with issues arising from the topology of the Stiefel manifold and how to inexpensively compute the change-of-measure terms. We introduce an auxiliary parameter approach that limits the impact of topological issues. We provide both analysis of our methods and numerical examples demonstrating the effectiveness of the approach. We also discuss how our Givens representation can be used to define general classes of distributions over the space of orthogonal matrices. We then give demonstrations on several examples showing how the Givens approach performs in practice in comparison with other methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Finite and infinite matrix product states for Gutzwiller projected mean-field wave functions

Matrix product states (MPS) and “dressed” ground states of quadratic mean fields (e.g., Gutzwiller projected Slater determinants) are both important classes of variational wave functions. This latter class has played important roles in understanding superconductivity and quantum spin liquids. We present a method to obtain both the finite and infinite MPS (iMPS) representation of the ground state of an arbitrary fermionic quadratic mean-field Hamiltonian (which in the simplest case is a Slater determinant and in the most general case is a Pfaffian). We also show how to represent products of such states (e.g., determinants times Pfaffians). From this representation one can project to single occupancy and evaluate the entanglement spectra after Gutzwiller projection. We then obtain the MPS and iMPS representation of Gutzwiller projected mean-field states that arise from the variational slave-fermion approach to the S = 1 bilinear-biquadratic quantum spin chain. To accomplish this, we develop an approach to orthogonalize degenerate iMPS to find all the states in the degenerate ground-state manifold. We find the energies of the MPS and iMPS states match the variational energies closely, indicating the method is accurate and there is minimal loss due to truncation error. We then present an exploration of the entanglement spectra of projected slave-fermion states, exploring their qualitative features and finding good qualitative agreement with the respective exact ground-state spectra found from density matrix renormalization group.

1-dimensional spin chains↗
DOE OSTI2012

GRAFEC: A New Spanish Program to Investigate Waste Management Options for Radioactive Graphite - 12399

Spain has to manage about 3700 tons of irradiated graphite from the reactor Vandellos I as radioactive waste. 2700 tons are the stack of the reactor and are still in the reactor core waiting for retrieval. The rest of the quantities, 1000 tons, are the graphite sleeves which have been already retrieved from the reactor. During operation the graphite sleeves were stored in a silo and during the dismantling stage a retrieval process was carried out separating the wires from the graphite, which were crushed and introduced into 220 cubic containers of 6 m{sup 3} each and placed in interim storage. The graphite is an intermediate level radioactive waste but it contains long lived radionuclides like {sup 14}C which disqualifies disposal at the low level waste repository of El Cabril. Therefore, a new project has been started in order to investigate two new options for the management of this waste type. The first one is based on a selective decontamination of {sup 14}C by thermal methods. This method is based on results obtained at the Research Centre Juelich (FZJ) in the Frame of the EC programs 'Raphael' and 'Carbowaste'. The process developed at FZJ is based on a preferential oxidation of {sup 14}C in comparison to the bulk {sup 12}C. Explanations for this effect are the inhomogeneous distribution and a weaker bounding of {sup 14}C which is not incorporated in the graphite lattice. However these investigations have only been performed with graphite from the high temperature reactor Arbeitsgemeinschaft Versuchsreaktor Juelich AVR which has been operated in a non-oxidising condition or research reactor graphite operated at room temperature. The reactor Vandellos I has been operated with CO{sub 2} as coolant and significant amounts of graphite have been already oxidised. The aim of the project is to validate whether a {sup 14}C decontamination can also been achieved with graphite from Vandellos I. A second possibility under investigation is the encapsulation of the graphite in a long term stable glass matrix. The principal applicability has been already proved by FNAG. Crushed graphite mixed with a suitable glass powder has been pressed at elevated temperature under vacuum. The vacuum is required to avoid gas enclosures in the obtained product. The obtained products, named IGM for 'Impermeable Graphite Matrix', have densities above 99% of theoretical density. The amount of glass has been chosen with respect to the pore volume of the former graphite parts. The method allows the production of encapsulated graphite without increasing the disposal volume. This paper will give a short overview of characterisation results of different irradiated graphite materials obtained at CIEMAT and in the Carbowaste project as well as the proposed methods and the actual status of the program including first results about leaching of non-radioactive IGM samples and hopefully first tendencies concerning the C-14 separation from graphite of Vandellos I by thermal treatment. Both processes, the thermal treatment as well as the IGM, have the potential to solve problems related to the management of irradiated graphite in Spain. However the methods have only been tested with different types of i-graphite and virgin graphite, respectively. Only investigations with real i-graphite from Spain will reveal whether the described methods are applicable to graphite from Vandellos I. However all partners are convinced that one of these new methods or a combination of them will lead to a feasible option to manage i-graphite in Spain on an industrial scale. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

Accurate determination of the valence band edge in hard x-ray photoemission spectra using GW theory

We introduce a new method for determining accurate values of the valence-band maximum in x-ray photoemission spectra. Specifically, we align the sharpest peak in the valence-band region of the experimental spectrum with the corresponding feature of a theoretical valence-band density of states curve from ab initio GW theory calculations. This method is particularly useful for soft and hard x-ray photoemission studies of materials with a mixture of valence-band characters, where strong matrix element effects can render standard methods for extracting the valence-band maximum unreliable. We apply our method to hydrogen-terminated boron-doped diamond, which is a promising substrate material for novel solar cell devices. By carrying out photoemission experiments with variable light polarizations, we verify the accuracy of our analysis and the general validity of the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Thermophysical properties of sintered aluminum-silicon

The effect of silicon addition to aluminum on the thermophysical properties of Al was experimentally investigated. Thermal diffusivity, heat capacity, density, and thermal expansion of sintered Al-Si samples with Si content ranging 0-8 wt% were measured. For U-Mo/Al dispersion fuel for research reactors, an addition of Si (up to 8 wt% Si) to the Al matrix has been adopted as a method to reduce the formation of interaction layer between the U-Mo and the Al. In this study, the Si-added Al matrix was simulated by using a powder metallurgical method, in which an Al powder and a Si powder were mixed and sintered. All measured properties decreased with the addition of Si. The effect was most pronounced in thermal diffusivity while the decreases in heat capacity and density were relatively insignificant. The measured thermal conductivity from the present study was higher than those measured data for alloy samples available in the literature to the differences in the fabrication method. All measured data including this study, however, were lower than the prediction by the Bruggeman model. It is because the model assumes an ideal heat-transfer conditions such as spherical particles, homogeneous particle distribution, perfect bonding between Al and Si. (C) 2018 Elsevier B.V. All rights reserved.

Cho, Tae Won↗
NASA NTRS1988

Quantum mechanical algebraic variational methods for inelastic and reactive molecular collisions

The quantum mechanical problem of reactive or nonreactive scattering of atoms and molecules is formulated in terms of square-integrable basis sets with variational expressions for the reactance matrix. Several formulations involving expansions of the wave function (the Schwinger variational principle) or amplitude density (a generalization of the Newton variational principle), single-channel or multichannel distortion potentials, and primitive or contracted basis functions are presented and tested. The test results, for inelastic and reactive atom-diatom collisions, suggest that the methods may be useful for a variety of collision calculations and may allow the accurate quantal treatment of systems for which other available methods would be prohibitively expensive.

Schwenke, David W.↗
DOE OSTI2014

Rashba-Zeeman-effect-induced spin filtering energy windows in a quantum wire

We perform a numerical study on the spin-resolved transport in a quantum wire (QW) under the modulation of both Rashba spin-orbit coupling (SOC) and a perpendicular magnetic field by using the developed Usuki transfer-matrix method in combination with the Landauer-Büttiker formalism. Wide spin filtering energy windows can be achieved in this system for unpolarized spin injection. In addition, both the width of energy window and the magnitude of spin conductance within these energy windows can be tuned by varying Rashba SOC strength, which can be apprehended by analyzing the energy dispersions and spin-polarized density distributions inside the QW, respectively. Further study also demonstrates that these Rashba-SOC-controlled spin filtering energy windows show a strong robustness against disorders. These findings may not only benefit to further understand the spin-dependent transport properties of a QW in the presence of external fields but also provide a theoretical instruction to design a spin filter device.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗