Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density matrix method”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
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↗
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.↗
NASA NTRS2000

Optical Control of Intersubband Absorption in a Multiple Quantum Well-Embedded Semiconductor Microcravity

Optical intersubband response of a multiple quantum well (MQW)-embedded microcavity driven by a coherent pump field is studied theoretically. The n-type doped MQW structure with three subbands in the conduction band is sandwiched between a semi-infinite medium and a distributed Bragg reflector (DBR). A strong pump field couples the two upper subbands and a weak field probes the two lower subbands. To describe the optical response of the MQW-embedded microcavity, we adopt a semi-classical nonlocal response theory. Taking into account the pump-probe interaction, we derive the probe-induced current density associated with intersubband transitions from the single-particle density-matrix formalism. By incorporating the current density into the Maxwell equation, we solve the probe local field exactly by means of Green's function technique and the transfer-matrix method. We obtain an exact expression for the probe absorption coefficient of the microcavity. For a GaAs/Al(sub x)Ga(sub 1-x)As MQW structure sandwiched between a GaAs/AlAs DBR and vacuum, we performed numerical calculations of the probe absorption spectra for different parameters such as pump intensity, pump detuning, and cavity length. We find that the probe spectrum is strongly dependent on these parameters. In particular, we find that the combination of the cavity effect and the Autler-Townes effect results in a triplet in the optical spectrum of the MQW system. The optical absorption peak value and its location can be feasibly controlled by varying the pump intensity and detuning.

Liu, Ansheng↗
DOE OSTI2025

Random Phase Approximation Correlation Energy Using Real-Space Density Functional Perturbation Theory

We present a real-space method for computing the random phase approximation (RPA) correlation energy within Kohn–Sham density functional theory, leveraging the low-rank nature of the frequency-dependent density response operator. In particular, we employ a cubic-scaling formalism based on density functional perturbation theory that circumvents the calculation of the response function matrix, instead relying on the ability to compute its product with a vector through the solution of the associated Sternheimer linear systems. We develop a large-scale parallel implementation of this formalism using the subspace iteration method in conjunction with the spectral quadrature method while employing the Kronecker product-based method for the application of the Coulomb operator and the conjugate orthogonal conjugate gradient method for the solution of the linear systems. We demonstrate convergence with respect to key parameters and verify the method’s accuracy by comparing with plane-wave results. We show that the framework achieves good strong scaling to many thousands of processors, reducing the time to solution for a lithium hydride system with 128 electrons to around 150 s on 4608 processors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2021

A Sulfur Dioxide Covariance-Based Retrieval Algorithm (COBRA): Application to TROPOMI Reveals New Emission Sources

Sensitive and accurate detection of sulfur dioxide (SO 2 ) from space is important for monitoring and estimating global sulfur emissions. Inspired by detection methods applied in the thermal infrared, we present here a new scheme to retrieve SO 2 columns from satellite observations of ultraviolet back-scattered radiances. The retrieval is based on a measurement error covariance matrix to fully represent the SO 2 -free radiance variability, so that the SO 2 slant column density is the only retrieved parameter of the algorithm. We demonstrate this approach, named COBRA, on measurements from the TROPOspheric Monitoring Instrument (TROPOMI) aboard the Sentinel-5 Precursor (S-5P) satellite. We show that the method reduces significantly both the noise and biases present in the current TROPOMI operational DOAS SO 2 retrievals. The performance of this technique is also benchmarked against that of the Principal Component Algorithm (PCA) approach. We find that the quality of the data is similar and even slightly better with the proposed COBRA approach. The ability of the algorithm to retrieve SO 2 accurately is also further supported by comparison with ground-based observations. We illustrate the great sensitivity of the method with a high-resolution global SO 2 map, considering two and a half years of TROPOMI data. In addition to the known sources, we detect many new SO 2 emission hotspots worldwide. For the largest sources, we use the COBRA data to estimate SO 2 emission rates. Results are comparable to other recently published TROPOMI-based SO 2 emissions estimates, but the associated uncertainties are significantly lower than with the operational data. Next, for a limited number of weak sources, we demonstrate the potential of our data for quantifying SO 2 emissions with a detection limit of about 8 kt yr -1 , a factor of 4 better than the emissions derived from the Ozone Monitoring Instrument (OMI). We anticipate that the systematic use of our TROPOMI COBRA SO 2 column data set at a global scale will allow identifying and quantifying missing sources, and help improving SO 2 emission inventories.

SO2↗
DOE OSTI2024

Entanglement and Bell inequalities with boosted t t ¯

The Large Hadron Collider provides a unique opportunity to study quantum entanglement and violation of Bell inequalities at the highest energy available today. In this paper, we will investigate these quantum correlations with top quark pair production, which represents a system of two-qubits. The semileptonic top pair channel has a factor of 6 increase in statistics and easier reconstruction with respect to the dileptonic channel. Although measuring the spin polarization of the hadronic top quark is known to be challenging, our study indicates that it is feasible to reconstruct the spin density matrix of the two-qubit system using an optimal hadronic polarimeter. This is achieved with the aid of jet substructure techniques and NN-inspired reconstruction methods, which improve the mapping between subjets and quarks. We find that entanglement can already be observed at more than the 5 σ level with existing data, and violation of Bell inequalities may be probed above the 4 σ level at the HL-LHC with 3 ab − 1 of data. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
NASA NTRS1974

Reaction formulation for radiation and scattering from plates, corner reflectors and dielectric-coated cylinders

The reaction concept is employed to formulate an integral equation for radiation and scattering from plates, corner reflectors, and dielectric-coated conducting cylinders. The surface-current density on the conducting surface is expanded with subsectional bases. The dielectric layer is modeled with polarization currents radiating in free space. Maxwell's equation and the boundary conditions are employed to express the polarization-current distribution in terms of the surface-current density on the conducting surface. By enforcing reaction tests with an array of electric test sources, the moment method is employed to reduce the integral equation to a matrix equation. Inversion of the matrix equation yields the current distribution, and the scattered field is then obtained by integrating the current distribution. The theory, computer program and numerical results are presented for radiation and scattering from plates, corner reflectors, and dielectric-coated conducting cylinders.

Wang, N. N.↗
NASA NTRS2003

Calculation of Free-Free Opacities

Free-free absorption is an important contribution to the opacity for radiation transport through hot materials Temperatures can be as high as several keV, such that it becomes a computational challenge to solve the Schrodinger equation efficiently for rapidly oscillating continuum functions for high angular momenta. Several groups\footnots, including ours, have studied the phase amplitude solution (PAS) of the Schrodinger equation, in which one solves equations for the wave function amplitude and phase, which are: smooth functions of the electron energy. It is also important to have an accurate Schroudinger benchmark for the development of the PAS method. We present results for dipole matrix elements, Gaunt factors, and cross sections for the absorption of radiation at various energies for Cs XIX at temperature=100 eV and density=0.187 g/cc for our newly developed PAS and Schrodinger benchmark.

Bhatia, A. K.↗
NASA NTRS2012

Mapping the Milky Way Galaxy with LISA

Gravitational wave detectors in the mHz band (such as the Laser Interferometer Space Antenna, or LISA) will observe thousands of compact binaries in the galaxy which can be used to better understand the structure of the Milky Way. To test the effectiveness of LISA to measure the distribution of the galaxy, we simulated the Close White Dwarf Binary (CWDB) gravitational wave sky using different models for the Milky Way. To do so, we have developed a galaxy density distribution modeling code based on the Markov Chain Monte Carlo method. The code uses different distributions to construct realizations of the galaxy. We then use the Fisher Information Matrix to estimate the variance and covariance of the recovered parameters for each detected CWDB. This is the first step toward characterizing the capabilities of space-based gravitational wave detectors to constrain models for galactic structure, such as the size and orientation of the bar in the center of the Milky Way

McKinnon, Jose A.↗
NASA NTRS2000

Equivalent Skin Analysis of Wing Structures Using Neural Networks

An efficient method of modeling trapezoidal built-up wing structures is developed by coupling. in an indirect way, an Equivalent Plate Analysis (EPA) with Neural Networks (NN). Being assumed to behave like a Mindlin-plate, the wing is solved using the Ritz method with Legendre polynomials employed as the trial functions. This analysis method can be made more efficient by avoiding most of the computational effort spent on calculating contributions to the stiffness and mass matrices from each spar and rib. This is accomplished by replacing the wing inner-structure with an "equivalent" material that combines to the skin and whose properties are simulated by neural networks. The constitutive matrix, which relates the stress vector to the strain vector, and the density of the equivalent material are obtained by enforcing mass and stiffness matrix equities with rec,ard to the EPA in a least-square sense. Neural networks for the material properties are trained in terms of the design variables of the wing structure. Examples show that the present method, which can be called an Equivalent Skin Analysis (ESA) of the wing structure, is more efficient than the EPA and still fairly good results can be obtained. The present ESA is very promising to be used at the early stages of wing structure design.

Liu, Youhua↗
NASA NTRS1990

Gaussian windows: A tool for exploring multivariate data

Presented here is a method for interactively exploring a large set of quantitative multivariate data, in order to estimate the shape of the underlying density function. It is assumed that the density function is more or less smooth, but no other specific assumptions are made concerning its structure. The local structure of the data in a given region may be examined by viewing the data through a Gaussian window, whose location and shape are chosen by the user. A Gaussian window is defined by giving each data point a weight based on a multivariate Gaussian function. The weighted sample mean and sample covariance matrix are then computed, using the weights attached to the data points. These quantities are used to compute an estimate of the shape of the density function in the window region. The local structure of the data is described by a method similar to the method of principal components. By taking many such local views of the data, we can form an idea of the structure of the data set. The method is applicable in any number of dimensions. The method can be used to find and describe simple structural features such as peaks, valleys, and saddle points in the density function, and also extended structures in higher dimensions. With some practice, we can apply our geometrical intuition to these structural features in any number of dimensions, so that we can think about and describe the structure of the data. Since the computations involved are relatively simple, the method can easily be implemented on a small computer.

Jaeckel, Louis A.↗
DOE OSTI2025

Exploiting a Shortcoming of Coupled-Cluster Theory: The Extent of Non-Hermiticity as a Diagnostic Indicator of Computational Accuracy

The fundamental non-Hermitian nature of the forms of the coupled-cluster (CC) theory widely used in quantum chemistry has usually been viewed as a negative, but the present paper shows how this can be used to an advantage. Specifically, the non-symmetric nature of the reduced one-particle density matrix (in the molecular orbital basis) is advocated as a diagnostic indicator of computational quality. In the limit of the full coupled-cluster theory [which is equivalent to full configuration interaction (FCI)], the electronic wave function and correlation energy are exact within a given one-particle basis set, and the symmetric character of the exact density matrix is recovered. The extent of the density matrix asymmetry is shown to provide a measure of “how difficult the problem is” (like the well-known T 1 diagnostic), but its variation with the level of theory also gives information about “how well this particular method works”, irrespective of the difficulty of the problem at hand. The proposed diagnostic is described and applied to a select group of small molecules, and an example of its overall utility for the practicing quantum chemist is illustrated through its application to the beryllium dimer (Be 2 ). Future application of this idea to excited states, open-shell systems, and symmetry-breaking problems and an extension of the method to the two-particle density are then proposed.

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

Ground-state-based model reduction with unitary circuits

Here, we present a method to numerically obtain low-energy effective models based on a unitary transformation of the ground state. The algorithm finds a unitary circuit that transforms the ground state of the original model to a projected wavefunction with only the low-energy degrees of freedom. The effective model can then be derived using the unitary transformation encoded in the circuit. We test our method on the one-dimensional and two-dimensional square-lattice Hubbard model at half-filling, and obtain more accurate effective spin models than the standard perturbative approach.

Hubbard model↗
DOE OSTI2026

Probing the Mechanism of Cadmium Sulfide Cluster Nucleation and Growth via Sequential Infiltration Synthesis

Few-atom metal chalcogenide clusters may be realized through the sequential infiltration of metal–organic precursors in polymer films. However, the underlying nucleation and growth mechanisms that allow for cluster synthesis with near atomic-scale precision are not fully resolved. The kinetics of the sequential infiltration synthesis (SIS) method that control the nucleation and growth mechanisms of primarily Cd 4 S 4 -core clusters within a poly­(4-vinylpyridine) (P4VP) matrix are probed with in situ UV–visible absorbance spectroscopy. Density functional theory (DFT) calculations that allow simulation of the optical properties of cluster fragments further reveal the thermodynamics that guide cluster growth within the P4VP matrix. Here, we conclude that a reactive capture mechanism for cluster nucleation and growth is dominant, although transient dimethyl cadmium adduction to the polymer backbone may contribute to cluster nucleation under shorter metal–organic purge process conditions. Grazing incidence X-ray diffraction (GI-XRD) and X-ray absorption spectroscopy (XAS) analyses further corroborate the cluster size and atom connectivity throughout the stepwise synthesis.

Jayaweera, Nuwanthaka P. [Argonne National Laborat↗
NASA NTRS2016

Influence of Constituents on Creep Properties of SiC/SiC Composites

SiC-SiC composites are being considered as potential candidate materials for next generation turbine components such as combustor liners, nozzle vanes and blades because of their low density, high temperature capability, and tailorable mechanical properties. These composites are essentially fabricated by infiltrating matrix into a stacked array of fibers or fiber preform by one or a combination of manufacturing methods such as, Melt Infiltration (MI) of molten silicon metal, Chemical Vapor Infiltration (CVI), Polymer Infiltration and Pyrolysis (PIP). To understand the influence of constituents, the SiC-SiC composites fabricated by MI, CVI, and PIP methods were creep tested in air between 12000 and 14500 degrees Centigrade for up to 500 hours. The failed specimens were analyzed under a scanning electron microscope to assess damage mechanisms. Also, knowing the creep deformation parameters of the fiber and the matrix under the testing conditions, the creep behavior of the composites was modeled and compared with the measured data. The implications of the results on the long term durability of these composites will be discussed.

modeling↗
NASA NTRS1988

Electron density-sensitive line ratios for Fe XII in the solar corona

The collision strengths for electron-impact excitation of fine-structure transitions from the ground 3s23p3 configuration to the excited 3s3p4 configuration in Fe XII are presented. These are calculated by transforming the K-matrices obtained in a nonrelativistic R-matrix method to K-matrices in intermediate coupling using term-coupling coefficients which take into account the effect of breakdown of target LS coupling. Extensive configuration-interaction wave functions are used to represent the target states. The effective collision strengths are determined assuming a Maxwellian distribution of electron energies. Rydberg series of resonances converging to the excited state thresholds are explicitly included in the calculation. The relative level populations and emission-line strengths for Fe XII are calculated using the new atomic data presented in this paper and that published by Tayal, Henry, and Pradhan (1987) and Tayal and Henry (1986) over a wide range of electron densities and temperatures (1-2 x 10 to the 6th K). Present results for collision strengths and emission-line strengths are compared with the previous calculation of Flower (1977). Some significant differences are noted.

Tayal, S. S.↗
NASA NTRS1972

Lumped parameter analysis of a stringer reinforced plate excited by band limited noise

The maximum root-mean-square response of a square clamped plate, subjected to a stationary random excitation, was determined both theoretically and empirically. For the tests, a 40 x 40 x 1/4 aluminum sheet was milled out to provide a thin waffle plate with 9 panels. The plate was acoustically excited, with the frequency range limited between 25 Hz and 500 Hz. The root-mean-square power of the random excitations, 149 decibels, was nearly constant for all natural frequencies of the plate system. Strain gage readings at 14 points were plotted as power spectral densities and root-mean-square displacements. A modal damping ratio matrix was determined from the response curves. Damping ratios were found to be a function of frequency, but not of location on the plate. For comparisons with the theoretical response, the plate was subdivided into 14 modes and anlayzed as a damped, lumped parameter system by use of an approximate normal-mode method. This method gave the best predictions of the power spectral densities for the lower frequencies. That prediction errors increased for the higher frequencies is attributed to improper assumptions for mass distribution and the insufficient number of lumped mass points.

Bilyeu, D. J.↗