Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix theories”

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 37 records · Page 2

Theorems on symmetrics and flux conservation in radiative transfer using the matrix operator theory

The matrix operator approach to radiative transfer is shown to be a very powerful technique in establishing symmetry relations for multiple scattering in inhomogeneous atmospheres. Symmetries are derived for the reflection and transmission operators using only the symmetry of the phase function. These results will mean large savings in computer time and storage for performing calculations for realistic planetary atmospheres using this method. The results have also been extended to establish a condition on the reflection matrix of a boundary in order to preserve reciprocity. Finally energy conservation is rigorously proven for conservative scattering in inhomogeneous atmospheres.

Kattawar, G. W.↗

Theorems on symmetries and flux conservation in radiative transfer using the matrix operator theory.

The matrix operator approach to radiative transfer is shown to be a very powerful technique in establishing symmetry relations for multiple scattering in inhomogeneous atmospheres. Symmetries are derived for the reflection and transmission operators using only the symmetry of the phase function. These results will mean large savings in computer time and storage for performing calculations for realistic planetary atmospheres using this method. The results have also been extended to establish a condition on the reflection matrix of a boundary in order to preserve reciprocity. Finally energy conservation is rigorously proven for conservative scattering in inhomogeneous atmospheres.

Kattawar, G. W.↗

Pressure-induced 𝐵⁢1 to 𝐵⁢2 phase transition in CeN studied by 𝑎⁢𝑏 𝑖⁢𝑛⁢𝑖⁢𝑡⁢𝑖⁢𝑜 correlation matrix renormalization theory calculations

We apply correlation matrix renormalization theory (CMRT) to cerium nitride (CeN) under pressure. For 𝐵⁢1 (NaCl-type) phase, CMRT gives an equation of state consistent with ambient pressure experiments. It produces electronic density-of-state (DOS) characterized by a sharp 4⁢𝑓 quasiparticle resonance peak pinned at the Fermi level and two subbands formed by strong hybridization between the localized Ce-4⁢𝑓 electrons and the itinerant Ce-5⁢𝑑 and N-2⁢𝑝 electrons below the Fermi level, consistent with x-ray photoemission spectroscopy experiments. Upon compression, CMRT predicts a first-order 𝐵⁢1 → 𝐵⁢2 (CsCl-type) transition with ∼11% volume collapse in agreement with experiments. Across the transition, the 4⁢𝑓 spectrum broadens, the 4⁢𝑓 orbital occupancy increases, and the hybridization with conduction states enhances, signaling a crossover from partially localized to more itinerant 4⁢𝑓 behavior. Furthermore, these features are in excellent agreement with experimental observations, demonstrating that CMRT provides a parameter-free description and prediction of correlation-driven structural and electronic transitions in rare-earth compounds.

Ab initio calculations↗

A principle of maximum ignorance for semiclassical gravity

The principle of maximum ignorance posits that the coarse-grained description of a system is maximally agnostic about its underlying microscopic structure. We briefly review this principle for random matrix theory and for the eigenstate thermalization hypothesis. We then apply this principle in holography to construct ensembles of random mixed states. This leads to an ensemble of microstates which models our microscopic ignorance, and which on average reproduces the effective semiclassical physics of a given bulk state. We call this ensemble the state-averaging ansatz. The output of our model is a prediction for semiclassical contributions to variances and higher statistical moments over the ensemble of microstates. The statistical moments provide coarse-grained — yet gravitationally non-perturbative — information about the microstructure of the individual states of the ensemble. We show that these contributions exactly match the on-shell action of known wormhole configurations of the gravitational path integral. These results strengthen the view that wormholes simply parametrize the ignorance of the microstructure of a fundamental state, given a fixed semiclassical bulk description.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Finite-temperature density matrix embedding theory

We describe a formulation of the density matrix embedding theory at finite temperature. We present a generalization of the ground-state bath orbital construction that embeds a mean-field finite-temperature density matrix up to a given order in the Hamiltonian, or the Hamiltonian up to a given order in the density matrix. In this work, we assess the performance of the finite-temperature density matrix embedding on the one-dimensional Hubbard model both at half-filling and away from it, and the two-dimensional Hubbard model at half-filling, comparing to exact data where available, as well as results from finite-temperature density matrix renormalization group, dynamical mean-field theory, and dynamical cluster approximations. The accuracy of finite-temperature density matrix embedding appears comparable to that of the ground-state theory, with, at most, a modest increase in bath size, and competitive with that of cluster dynamical mean-field theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Computational Protocol to Evaluate Electron–Phonon Interactions Within Density Matrix Perturbation Theory

We present a computational protocol, based on density matrix perturbation theory, to obtain non-adiabatic, frequency-dependent electron–phonon self-energies for molecules and solids. Our approach enables the evaluation of electron–phonon interaction using hybrid functionals, for spin-polarized systems, and the computational overhead to include dynamical and non-adiabatic terms in the evaluation of electron–phonon self-energies is negligible. We discuss results for molecules, as well as pristine and defective solids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Matrix operator theory of radiative transfer. II - Scattering from maritime haze.

Matrix operator theory is used to calculate the reflected and transmitted radiance of photons that have interacted with plane-parallel maritime haze layers. The results are presented for three solar zenith angles, three values of the surface albedo, and a range of optical thicknesses from very thin to very thick. The diffuse flux at the lower boundary and the cloud albedo are tabulated. The forward peak and other features in the single-scattered phase function cause the radiance in many cases to be very different from that for Rayleigh scattering. In particular, the variation of the radiance with both the zenith or nadir angle and the azimuthal angle is more marked and the relative limb darkening under very thick layers is greater for haze M than for Rayleigh scattering. The downward diffuse flux at the lower boundary for A = 0 is always greater and the cloud albedo is always less for haze M than for Rayleigh layers.

Kattawar, G. W.↗

Ab initio description of bcc iron with correlation matrix renormalization theory

We applied the ab initio spin-polarized correlation matrix renormalization theory to the ferromagnetic state of the bulk bcc iron. We showed that it was capable of reproducing the equilibrium physical properties and the pressure-volume curve in good comparison with experiments. We then focused on the analysis of its local electronic correlations. By exploiting different local fluctuation-related physical quantities as measures of electronic correlation within target orbits, we elucidated the different roles of t 2g and e g states in both spin channels and presented compelling evidence to showcase this distinction in their electronic correlation.

3-dimensional systems↗

Real-time scattering in Ising field theory using matrix product states

We study scattering in Ising field theory (IFT) using matrix product states and the time-dependent variational principle. IFT is a one-parameter family of strongly coupled nonintegrable quantum field theories in 1+1 dimensions, interpolating between massive free fermion theory and Zamolodchikov's integrable massive 𝐸 8 theory. Particles in IFT may scatter either elastically or inelastically. In the postcollision wave function, particle tracks from all final-state channels occur in superposition; processes of interest can be isolated by projecting the wave function onto definite particle sectors, or by evaluating energy density correlation functions. Using numerical simulations we determine the time delay of elastic scattering and the probability of inelastic particle production as a function of collision energy. We also study the mass and width of the lightest resonance near the 𝐸 8 point in detail. Close to both the free fermion and 𝐸 8 theories, our results for both elastic and inelastic scattering are in good agreement with expectations from form-factor perturbation theory. Using numerical computations to go beyond the regime accessible by perturbation theory, we find that the high-energy behavior of the two-to-two particle scattering probability in IFT is consistent with a conjecture of Zamolodchikov. Our results demonstrate the efficacy of tensor-network methods for simulating the real-time dynamics of strongly coupled quantum field theories in 1+1 dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ab initio study on fcc Pr with correlation matrix renormalization theory

Here, we studied fcc praseodymium (Pr) with the ab initio correlation matrix renormalization theory (CMRT) explicitly calculating the Coulomb interactions among basis orbitals and without using any adjustable parameters to work for strongly correlated electron systems. We calculated its total energy in a paramagnetic ground state and studied the role of the correlated 4f electrons in the system. Good agreement was obtained between CMRT and experiments in the pressure volume dependence of the fcc phase. We also compared the CMRT results against other theoretical methods including local density approximation + dynamical mean-field theory and showed consistent results among them. Moreover, we found the normalized local charge fluctuation of the 4f electrons as the function of volume exhibits a clear slope change at the volume collapse region, indicating a switch in their correlation nature there.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Local Excitations of a Charged Nitrogen Vacancy in Diamond with Multireference Density Matrix Embedding Theory

Here, we investigate the negatively charged nitrogen-vacancy center in diamond using periodic density matrix embedding theory (pDMET). To describe the strongly correlated excited states of this system, the complete active space self-consistent field (CASSCF) followed by n-electron valence state second-order perturbation theory (NEVPT2) was used as the impurity solver. Since the NEVPT2-DMET energies show a linear dependence on the inverse of the size of the embedding subspace, we performed an extrapolation of the excitation energies to the nonembedding limit using a linear regression. The extrapolated NEVPT2-DMET first triplet–triplet excitation energy is 2.31 eV and that for the optically inactive singlet–singlet transition is 1.02 eV, both in agreement with the experimentally observed vertical excitation energies of ∼2.18 eV and ∼1.26 eV, respectively. This is the first application of pDMET to a charged periodic system and the first investigation of the NV – defect using NEVPT2 for periodic supercell models.

36 MATERIALS SCIENCE↗

Postbuckling analysis using a general purpose code

A new capability for solving postbuckling problems in shell structures is described. The matrix theory to adapt Newton's method to nonlinear finite element shell analysis is outlined first. The matrix theory is directed at writing consistent linear algebratic equations for problems where the tangent stiffness matrix is singular or nearly singular. The matrix theory suggests a change of variables as part of the usual iterative procedure in Newton's method. The change of variables is shown to be feasible for introduction into the algorithm programmed in general purpose codes for finite element analysis of structures. Numerical results from a new option that has been programmed in an existing general purpose code are presented. The analysis of shell structures for collapse and for branching at bifurcation loads is illustrated by the numerical examples.

Thurston, G. A.↗

Matrix operator theory of radiative transfer. I - Rayleigh scattering.

An entirely rigorous method for the solution of the equations for radiative transfer based on the matrix operator theory is reviewed. The advantages of the present method are: (1) all orders of the reflection and transmission matrices are calculated at once; (2) layers of any thickness may be combined, so that a realistic model of the atmosphere can be developed from any arbitrary number of layers, each with different properties and thicknesses; (3) calculations can readily be made for large optical depths and with highly anisotropic phase functions; (4) results are obtained for any desired value of the surface albedo including the value unity and for a large number of polar and azimuthal angles; (5) all fundamental equations can be interpreted immediately in terms of the physical interactions appropriate to the problem; and (6) both upward and downward radiance can be calculated at interior points from relatively simple expressions.

Plass, G. N.↗

Correlation matrix renormalization theory in multi-band lattice systems

An appropriate treatment of electronic correlation effects plays an important role in accurate descriptions of physical and chemical properties of real materials. The recently proposed correlation matrix renormalization theory with sum rule correction (CMR) for studying correlated-electron materials has shown good performance in molecular systems and a periodic hydrogen chain in comparison with various quantum chemistry and quantum Monte Carlo calculations. Additionally, this work gives a detailed formulation and computational code implementation of CMR in multi-band periodic lattice systems. This lattice CMR ab initio theory is highly efficient, has no material specific adjustable parameters, and has no double counting issues faced by the hybrid approaches like LDA + U, DFT + DMFT and DFT + GA type theories. Benchmark studies on materials with s and p orbitals in this study show that CMR in its current implementation consistently performs well for these systems as the electron correlation increases from the bonding region to the bond breaking region.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Role of Coulomb interaction in the phase formation of fcc Ce: Correlation matrix renormalization theory

The effect of electronic Coulomb interaction on the phase formation of fcc Ce lattice is investigated by full ab initio calculations without adjustable Coulomb U and J parameters using the Gutzwiller wavefunction-based correlation matrix renormalization theory (CMRT). Its total energy and pressure as a function of volume agree reasonably well with existing DFT+Gutzwiller calculations and experiments, indicating correct capture of electronic correlation and screening effects within the CMRT formalism. Here, a stable phase is found in line with the experimental α-Ce phase, and a lurking phase is identified supposedly linked with the experimental γ-Ce phase. A criterion based on the local 4f electron charge fluctuation is introduced to confirm the distinct electronic correlation natures of both phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Excited States of Crystalline Point Defects with Multireference Density Matrix Embedding Theory

Accurate and affordable methods to characterize the electronic structure of solids are important for targeted materials design. Embedding-based methods provide an appealing balance in the trade-off between cost and accuracy-particularly when studying localized phenomena. Here, we use the density matrix embedding theory (DMET) algorithm to study the electronic excitations in solid-state defects with a restricted open-shell Hartree-Fock (ROHF) bath and multireference impurity solvers, specifically, complete active space self-consistent field (CASSCF) and n-electron valence state second-order perturbation theory (NEVPT2). In this work, we apply the method to investigate the electronic excitations in an oxygen vacancy (OV) on a MgO(100) surface and find absolute deviations within 0.05 eV between DMET using the CASSCF/NEVPT2 solver, denoted as CAS-DMET/NEVPT2-DMET, and the nonembedded CASSCF/NEVPT2 approach. Next, we establish the practicality of DMET by extending it to larger supercells for the OV defect and a neutral silicon vacancy in diamond where the use of nonembedded CASSCF/NEVPT2 is extremely expensive.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamics of disordered mechanical systems with large connectivity, free probability theory, and quasi-Hermitian random matrices

Disordered mechanical systems with high connectivity represent a limit opposite to the more familiar case of disordered crystals. Individual ions in a crystal are subjected essentially to nearest-neighbor interactions. In contrast, the systems studied in this paper have all their degrees of freedom coupled to each other. Thus, the problem of linearized small oscillations of such systems involves two full positive-definite and non-commuting matrices, as opposed to the sparse matrices associated with disordered crystals. Consequently, the familiar methods for determining the averaged vibrational spectra of disordered crystals, introduced many years ago by Dyson and Schmidt, are inapplicable for highly connected disordered systems. In this paper we apply random matrix theory (RMT) to calculate the averaged vibrational spectra of such systems, in the limit of infinitely large system size. At the heart of our analysis lies a calculation of the average spectrum of the product of two positive definite random matrices by means of free probability theory techniques. We also show that this problem is intimately related with quasi-hermitian random matrix theory (QHRMT), which means that the ‘hamiltonian’ matrix is hermitian with respect to a non-trivial metric. This extends ordinary hermitian matrices, for which the metric is simply the unit matrix. The analytical results we obtain for the spectrum agree well with our numerical results. The latter also exhibit oscillations at the high-frequency band edge, which fit well the Airy kernel pattern. We also compute inverse participation ratios of the corresponding amplitude eigenvectors and demonstrate that they are all extended, in contrast with conventional disordered crystals. Finally, we compute the thermodynamic properties of the system from its spectrum of vibrations. In addition to matrix model analysis, we also study the vibrational spectra of various multi-segmented disordered pendula, as concrete realizations of highly connected mechanical systems. A universal feature of the density of vibration modes, common to both pendula and the matrix model, is that it tends to a non-zero constant at vanishing frequency.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗