Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Random Matrix Theory”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Model dynamics for quantum computing

A model master equation suitable for quantum computing dynamics is presented. In an ideal quantum computer (QC), a system of qubits evolves in time unitarily and, by virtue of their entanglement, interfere quantum mechanically to solve otherwise intractable problems. In the real situation, a QC is subject to decoherence and attenuation effects due to interaction with an environment and with possible short-term random disturbances and gate deficiencies. The stability of a QC under such attacks is a key issue for the development of realistic devices. We assume that the influence of the environment can be incorporated by a master equation that includes unitary evolution with gates, supplemented by a Lindblad term. Lindblad operators of various types are explored; namely, steady, pulsed, gate friction, and measurement operators. In the master equation, we use the Lindblad term to describe short time intrusions by random Lindblad pulses. The phenomenological master equation is then extended to include a nonlinear Beretta term that describes the evolution of a closed system with increasing entropy. An external Bath environment is stipulated by a fixed temperature in two different ways. Here we explore the case of a simple one-qubit system in preparation for generalization to multi-qubit, qutrit and hybrid qubit–qutrit systems. This model master equation can be used to test the stability of memory and the efficacy of quantum gates. The properties of such hybrid master equations are explored, with emphasis on the role of thermal equilibrium and entropy constraints. Several significant properties of time-dependent qubit evolution are revealed by this simple study. - Highlights: • General realistic density matrix evolution model shown for one qubit. • Unitary pulses & bias yield quantum gates without harmful phase accumulations. • Gate friction, dissipation included by random Lindblad noise pulses. • Bath & entropy terms yield equilibrium; test algorithms, gates act on Larmor grid. • Measurement modeled by strong, fast Lindblad pulses.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Impact of the stacking sequence on the bandgap and luminescence properties of bulk, bilayer, and monolayer hexagonal boron nitride

We examine the effects of stacking sequence and number of layers on the electronic and luminescence properties of hexagonal boron nitride (h-BN) structures with first-principles calculations based on density functional and many-body perturbation theory. We explored the variations of the magnitude and character (direct or indirect) of the quasiparticle bandgap and interband optical matrix elements for bulk, bilayer, and monolayer stacking polytypes. Although the fundamental gap for most structures is indirect, phonon-assisted transitions are strong (typically 600 times stronger than bulk Si) and enable efficient deep-ultraviolet (UV) luminescence. The polarization of the emitted light is transverse electric, which facilitates light extraction perpendicularly to the h-BN basal plane. Random stacking in turbostratic BN breaks the crystal symmetry and enables optical transitions across the quasi-direct bandgap, albeit with a weak matrix element. Our results demonstrate that h-BN is a promising material for efficient deep-UV light emitters.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Radiative and Nonradiative Recombinations in Organic Radical Emitters: The Effect of Guest–Host Interactions

Radical-carrying organic molecules have received significant attention to bypass the issue related to harvesting triplet excitons in current light-emitting materials. While the computational efforts conducted so far have treated these radical emitters as isolated entities, in actual devices, they are embedded in a host matrix and subject to emitter–host interactions. In this study, by combining molecular dynamics simulations and density functional theory calculations, the impact of the host matrix on the optoelectronic performance of radical emitters is evaluated, taking as a representative example the (4-ncarbazolyl-2,6-dichlorophenyl)bis(2,4,6-trichlorophenyl)-methyl (TTM-3NCz) radical emitter dispersed in a 4,4-bis(carbazol-9-yl)biphenyl (CBP) host. A morphological analysis shows that steric effects around the radical centers, carried by the TTM electron-poor moieties of the emitters, disfavor π–π interactions with the host molecules, which leads to random intermolecular orientations around the TTM moieties. The 3NCz electron-rich moieties of the emitters, however, have much lesser spatial hindrance for intermolecular π–π stacking, which modulates the structural and electronic properties of the emitters in the host matrix. The influence of dynamic and static disorders on the radiative and nonradiative recombination processes is also investigated and it is found that the rates of nonradiative recombination are small, which opens the way to 100% internal quantum efficiency for the doublet-based emission process.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Phase-field based Multiscale Modeling of Heterogeneous Solid Electrolytes: Applications to Nanoporous Li 3 PS 4

Modeling the effective ion conductivities of heterogeneous solid electrolytes typically involves the use of a computer-generated microstructure consisting of randomly or uniformly oriented fillers in a matrix. But, the structural features of the filler/matrix interface, which critically determine the interface ion conductivity and the microstructure morphology, have not been considered during the microstructure generation. In using nanoporous β-Li 3 PS 4 electrolyte as an example, we develop a phase-field model that enables generating nanoporous microstructures of different porosities and connectivity patterns based on the depth and the energy of the surface (pore/electrolyte interface), both of which are predicted through density functional theory (DFT) calculations. Room-temperature effective ion conductivities of the generated microstructures are then calculated numerically, using DFT-estimated surface Li-ion conductivity (3.14×10 -3 S/cm) and experimentally measured bulk Li-ion conductivity (8.93×10 -7 S/cm) of β-Li 3 PS 4 as the inputs. We also use the generated microstructures to inform effective medium theories to rapidly predict the effective ion conductivity via analytical calculations. Furthemore, when porosity approaches the percolation threshold, both the numerical and analytical methods predict a significantly enhanced Li-ion conductivity (1.74×10 -4 S/cm) that is in good agreement with experimental data (1.64×10 -4 S/cm). The present phase-field based multiscale model is generally applicable to predict both the microstructure patterns and the effective properties of heterogeneous solid electrolytes.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Development of the self-consistent approximation and its application to the problem of magnetoelastic resonance in an inhomogeneous medium

Our previously proposed approximation involving both the first and second terms of the expansion of the vertex function is generalized to the system of two interacting wavefields of different physical nature. A system of self-consistent equations for the matrix Green’s function and matrix vertex function is derived. On the basis of this matrix generalization of the new self-consistent approximation, a theory of magnetoelastic resonance is developed for a ferromagnetic model, where the magnetoelastic coupling parameter ε(x) is inhomogeneous. Equations for magnetoelastic resonance are analyzed for one-dimensional inhomogeneities of the coupling parameter. The diagonal and off-diagonal elements of the matrix Green’s function of the system of coupled spin and elastic waves are calculated with the change in the ratio between the average value ε and rms fluctuation Δε of the coupling parameter between waves from the homogeneous case (ε ≠ 0, Δε = 0) to the extremely randomized case (ε = 0, Δε ≠ 0) at various correlation wavenumbers of inhomogeneities k{sub c}. For the limiting case of infinite correlation radius (k{sub c} = 0), in addition to approximate expressions, exact analytical expressions corresponding to the summation of all diagrams of elements of the matrix Green’s function are obtained. The results calculated for an arbitrary k{sub c} value in the new self-consistent approximation are compared to the results obtained in the standard self-consistent approximation, where only the first term of the expansion of the vertex function is taken into account. It is shown that the new approximation corrects disadvantages of the Green’s functions calculated in the standard approximation such as the dome shape of resonances and bends on the sides of resonance peaks. The appearance of a fine structure of the spectrum in the form of a narrow resonance on the Green’s function of spin waves and a narrow antiresonance on the Green’s function of elastic waves, which was previously predicted in the standard self-consistent approximation, is confirmed. With an increase in the parameter k{sub c}, the Green’s functions calculated in the standard and new approximations approach each other and almost coincide with each other at k{sub c}/k ≥ 0.5. At the same time, the results of this work indicate that the new self-consistent approximation has a certain advantage for studying the problems of stochastic radiophysics in media with long-wavelength inhomogeneities (small k{sub c} values), because it describes both the shape and width of peaks much better than the standard approximation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

QuTree: A tree tensor network package

Here we present QuTree, a C++ library for tree tensor network approaches. QuTree provides class structures for tensors, tensor trees, and related linear algebra functions that facilitate the fast development of tree tensor network approaches such as the multilayer multiconfigurational time-dependent Hartree approach or the density matrix renormalization group approach and its various extensions. We investigate the efficiency of relevant tensor and tensor network operations and show that the overhead for managing the network structure is negligible, even in cases with a million leaves and small tensors. QuTree focuses on providing simple, high-level routines while retaining easy access to the backend to facilitate novel developments. We demonstrate the capabilities of the package by computing the eigenstates of coupled harmonic oscillator Hamiltonians and performing random circuit simulations on a virtual quantum computer.

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

Quantum theory of a two-mode open-cavity laser

We develop the quantum theory of an open-cavity laser assuming that only two modes compete for gain. We show that the modes interact to build up a composite mode that becomes the lasing mode when pumping exceeds a threshold. This composite mode exhibits all the features of a typical laser mode, whereas its precise behavior depends explicitly on the openness of the cavity. We approach the problem by using the density-matrix formalism and derive the master equation for the light field. Our results are of particular interest in the context of random laser systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Progress on Complex Langevin simulations of a finite density matrix model for QCD

We study the Stephanov model, which is an RMT model for QCD at finite density, using the Complex Langevin algorithm. Naive implementation of the algorithm shows convergence towards the phase quenched or quenched theory rather than to intended theory with dynamical quarks. A detailed analysis of this issue and a potential resolution of the failure of this algorithm are discussed. We study the effect of gauge cooling on the Dirac eigenvalue distribution and time evolution of the norm for various cooling norms, which were specifically designed to remove the pathologies of the complex Langevin evolution. The cooling is further supplemented with a shifted representation for the random matrices. Unfortunately, none of these modifications generate a substantial improvement on the complex Langevin evolution and the final results still do not agree with the analytical predictions.

Bloch, Jacques↗
DOE OSTI2019

Negative specific heat from non-planar interactions and small black holes in AdS/CFT

The gravity side of the gauge/gravity duality predicts the existence of small black holes with negative specific heat. A free theory of strings has a Hagedorn behavior, but it does not lead to negative specific heat. To understand such states one needs to consider a theory of interacting strings. In the dual gauge theory, the string interactions are related to non-planar diagrams. In this paper the simplest gauge matrix model of two free matrices, that has Hagedorn behavior is analyzed in detail. A simple double trace deformation of the Hamiltonian, proportional to square of the free Hamiltonian square with a negative sign that mimics a gravitational attraction is enough to produce states with negative specific heat perturbatively and one can still compute the equation of state relating the entropy and the energy. A more general argument based on non-planar interactions that are random and that grow faster in strength than the energy suggests that states with negative specific heat appear generically.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

SU-E-J-83: Ion Imaging to Better Estimate In-Vivo Relative Stopping Powers Using X-Ray CT Prior-Knowledge Information

Purpose: To reduce uncertainties in relative stopping power (RSP) estimates for ions (alpha and carbon) by using Ion radiographic-imaging and X-ray CT prior-knowledge. Methods: A 36×36 phantom matrix composed of 9 materials with different thicknesses and randomly placed is generated. Theoretical RSPs are calculated using stopping power (SP) data from three references (Janni, ICRU49 and Bischel). We introduced an artificial systematic error (1.5%, 2.5% or 3.5%) and a random error (<0.5%) to the SP to simulated patient ion-range errors present in clinic environment. Carbon/alpha final energy for each RSPs set (theoretical and from CT images) is obtained with a ray-tracing algorithm. A gradient descent (GD) method is used to minimize the difference in exit particle energy, between theory and X-ray CT RSP maps, by iteratively correcting the RSP map from X-ray CT. Once a new set of RSPs is obtained for a direction a new optimization is done over other direction using the RSPs from the previous optimization. Theoretical RSPs are compared with experimental RSPs obtained with Gammex Phantom. Results: Preliminary results show that optimized RSP values can be obtained with smaller uncertainties (<1%) than clinical RSPs (1.5% to 3.5%). Theoretical values from three different references show uncertainties, up to 3% from experimental values. Further investigation will consider prior-knowledge from RSP obtained with CT images and ion radiographies from Monte Carlo Simulations. Conclusion: GD and ray-tracing methods have been implemented to reduce RSP uncertainties from values obtained for clinical treatment. Experimental RSPs will be obtained using carbon/alpha beams to consider the existence of material dependent systematic errors. Based on the results, it is hoped to show that using ray-tracing optimization with ion radiography and prior knowledge on RPSs, treatment planning accuracy and cost-effectiveness can be improved.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Pseudorandomness from Subset States

We show it is possible to obtain quantum pseudorandomness and pseudoentanglement from random subset states -- i.e. quantum states which are equal superpositions over (pseudo)random subsets of strings. This answers an open question of Aaronson et al. [arXiv:2211.00747], who devised a similar construction augmented by pseudorandom phases. Our result follows from a direct calculation of the trace distance between t copies of random subset states and the Haar measure, via the representation theory of the symmetric group. We show that the trace distance is negligibly small, as long as the subsets are of an appropriate size which is neither too big nor too small. In particular, we analyze the action of basis permutations on the symmetric subspace, and show that the largest component is described by the Johnson scheme: the double-cosets of the symmetric group S N by the subgroup S t x S N-t . The Gelfand pair property of this setting implies that the matrix eigenbasis coincides with the symmetric group irreducible blocks, with the largest eigenblock asymptotically approaching the Haar average. An immediate corollary of our result is that quantum pseudorandom and pseudoentangled state ensembles do not require relative phases.

Computational Complexity (cs.CC)↗
DOE OSTI2016

First-principles thermodynamics and defect kinetics guidelines for engineering a tailored RRAM device

Resistive Random Access Memories are among the most promising candidates for the next generation of non-volatile memory. Transition metal oxides such as HfOx and TaOx attracted a lot of attention due to their CMOS compatibility. Furthermore, these materials do not require the inclusion of extrinsic conducting defects since their operation is based on intrinsic ones (oxygen vacancies). Using Density Functional Theory, we evaluated the thermodynamics of the defects formation and the kinetics of diffusion of the conducting species active in transition metal oxide RRAM materials. The gained insights based on the thermodynamics in the Top Electrode, Insulating Matrix and Bottom Electrode and at the interfaces are used to design a proper defect reservoir, which is needed for a low-energy reliable switching device. The defect reservoir has also a direct impact on the retention of the Low Resistance State due to the resulting thermodynamic driving forces. The kinetics of the diffusing conducting defects in the Insulating Matrix determine the switching dynamics and resistance retention. The interface at the Bottom Electrode has a significant impact on the low-current operation and long endurance of the memory cell. Our first-principles findings are confirmed by experimental measurements on fabricated RRAM devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Quickest Detection for Changes in Maximal kNN Coherence of Random Matrices

This paper addresses the problem of quickest detection of a change in the maximal coherence between columns of a n p random matrix based on a sequence of matrix observations having a single unknown change point. The random matrix is assumed to have identically distributed rows and the maximal coherence is defined as the largest of the ( 2 p) correlation coefficients associated with any row. Likewise the κ nearest neighbor (kNN) coherence is defined as the κ-th largest of these correlation coefficients. The forms of the preand postchange distributions of the observed matrices are assumed to belong to the family of elliptically contoured densities with sparse dispersion matrices but are otherwise unknown. A nonparametric stopping rule is proposed that is based on the maximal kNN sample coherence between columns of each observed random matrix. This is a summary statistic that is related to a test of existence of a hub vertex in a sample correlation graph having degree at least κ. Performance bounds on the delay and false alarm performance of the proposed stopping rule are obtained in the purely high-dimensional regime where p → ∞ and n is fixed. When the prechange dispersion matrix is diagonal it is shown that, among all functions of the proposed summary statistic, the proposed stopping rule is asymptotically optimal under a minimax quickest change detection model as the stopping threshold approaches infinity. Furthermore, the theory developed also applies to sequential hypothesis testing and fixed sample size tests.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2021

Modeling and simulation of vascular tumors embedded in evolving capillary networks

Here in this work, we present a coupled 3D–1D model of solid tumor growth within a dynamically changing vascular network to facilitate realistic simulations of angiogenesis. Additionally, the model includes erosion of the extracellular matrix, interstitial flow, and coupled flow in blood vessels and tissue. We employ continuum mixture theory with stochastic Cahn–Hilliard type phase-field models of tumor growth. The interstitial flow is governed by a mesoscale version of Darcy’s law. The flow in the blood vessels is controlled by Poiseuille flow, and Starling’s law is applied to model the mass transfer in and out of blood vessels. The evolution of the network of blood vessels is orchestrated by the concentration of the tumor angiogenesis factors (TAFs); blood vessels grow towards the increasing TAFs concentrations. This process is not deterministic, allowing random growth of blood vessels and, therefore, due to the coupling of nutrients in tissue and vessels, makes the growth of tumors stochastic. We demonstrate the performance of the model by applying it to a variety of scenarios. Numerical experiments illustrate the flexibility of the model and its ability to generate satellite tumors. Simulations of the effects of angiogenesis on tumor growth are presented as well as sample-independent features of cancer.

3D–1D coupled blood flow models↗
DOE OSTI2016

Purcell effect at the percolation transition

Here, we investigate the spontaneous emission rate of a two-level quantum emitter next to a composite medium made of randomly distributed metallic inclusions embedded in a dielectric host matrix. In the near field, the Purcell factor can be enhanced by two orders of magnitude relative to the case of a homogeneous metallic medium and reaches its maximum precisely at the insulator-metal transition. By unveiling the role of the decay pathways in the emitter's lifetime, we demonstrate that, close to the percolation threshold, the radiation emission process is dictated by electromagnetic absorption in the heterogeneous medium. We show that our findings are robust against change in material properties and shape of inclusions and apply for different effective-medium theories as well as for a wide range of transition frequencies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Statistical Hauser-Feshbach theory with width-fluctuation correction including direct reaction channels for neutron-induced reactions at low energies

Here, a model to calculate particle-induced reaction cross sections with statistical Hauser-Feshbach theory including direct reactions is given. The energy average of the scattering matrix from the coupled-channels optical model is diagonalized by the transformation proposed by Engelbrecht and Weidenmüller [C. A. Engelbrecht and H. A. Weidenmüller, Phys. Rev. C 8, 859 (1973)]. The ensemble average of S -matrix elements in the diagonalized channel space is approximated by a model of Moldauer [P. A. Moldauer, Phys. Rev. C 12, 744 (1975)] using the newly parametrized channel degree-of-freedom ν a to better describe the Gaussian orthogonal ensemble (GOE) reference calculations. The Moldauer approximation is confirmed by a Monte Carlo study using a randomly generated S matrix, as well as the GOE threefold integration formula. Lastly, the method proposed is applied to the 238U(n,n') cross-section calculation in the fast-energy range, showing an enhancement in the inelastic scattering cross sections.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Uncertainty quantification in (α,n) neutron source calculations for an oxide matrix

Here we present a methodology to propagate nuclear data covariance information in neutron source calculations from (α,n) reactions. The approach is applied to estimate the uncertainty in the neutron generation rates for uranium oxide fuel types due to uncertainties on 1) 17,18 O( α,n) reaction cross sections and 2) uranium and oxygen stopping power cross sections. The procedure to generate reaction cross section covariance information is based on the Bayesian fitting method implemented in the R-matrix SAMMY code. The evaluation methodology uses the Reich-Moore approximation to fit the 17,18 O(α,n) reaction cross-sections in order to derive a set of resonance parameters and a related covariance matrix that is then used to calculate the energydependent cross section covariance matrix. The stopping power cross sections and related covariance information for uranium and oxygen were obtained by the fit of stopping power data in the -energy range of 1 keV up to 12 MeV. Cross section perturbation factors based on the covariance information relative to the evaluated 17,18 O( α,n) reaction cross sections, as well as uranium and oxygen stopping power cross sections, were used to generate a varied set of nuclear data libraries used in SOURCES4C and ORIGEN for inventory and source term calculations. The set of randomly perturbed output (α,n) source responses, provide the mean values and standard deviations of the calculated responses reflecting the uncertainties in nuclear data used in the calculations. Lastly, the results and related uncertainties are compared with experiment thick target (α,n) yields for uranium oxide.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Matrix Diffusion in Fractured Media: New Insights Into Power Law Scaling of Breakthrough Curves

Here, we develop a theoretical model for power law tailing behavior of transport in fractured rock based on the relative dominance of the decay rate of the advective travel time distribution, modeled using a Pareto distribution (with tail decaying as ~ time -(1+α) ), versus matrix diffusion, modeled using a Lévy distribution. The theory predicts that when the advective travel time distribution decays sufficiently slowly (α<1), the late-time decay rate of the breakthrough curve is -(1+α/2) rather than the classical -3/2. Yet, if α>1, the -3/2 decay rate is recovered. For weak matrix diffusion or short advective first breakthrough times, we identify an early-time regime where the breakthrough curve follows the Pareto distribution, before transitioning to the late-time decay rate. The theoretical predictions are validated against particle tracking simulations in the three-dimensional discrete fracture network simulator dfnWorks, where matrix diffusion is incorporated using a time domain random walk.

58 GEOSCIENCES↗