Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Extension of frozen natural orbital approximation to open-shell references: Theory, implementation, and application to single-molecule magnets

Natural orbitals are often used to achieve a more compact representation of correlated wave-functions. Using natural orbitals computed as eigenstates of the virtual–virtual block of the state density matrix instead of the canonical Hartree–Fock orbitals results in smaller errors when the same fraction of virtual space is frozen. This strategy, termed frozen natural orbital (FNO) approach, is effective in reducing the cost of regular coupled-cluster (CC) calculations and some multistate methods, such as EOM-IP-CC (equation-of-motion CC for ionization potentials). This contribution extends the FNO approach to the EOM-SF-CC ansatz (EOM-CC with spin-flip). In contrast to EOM-IP-CCSD, EOM-SF-CCSD relies on high-spin open-shell references. Using FNOs computed for an open-shell reference leads to an erratic behavior of the EOM-SF-CC energies and properties due to an inconsistent truncation of the α and β orbital spaces. A general solution to problems arising in the EOM-CC calculations utilizing open-shell references, termed OSFNO (open-shell FNO), is proposed. Furthermore, by means of singular value decomposition (SVD) of the overlap matrix of the α and β orbitals, the OSFNO algorithm identifies the corresponding orbitals and determines virtual orbitals corresponding to the singly occupied space. This is followed by SVD of the singlet part of the state density matrix in the remaining virtual orbital subspace. The so-computed FNOs preserve the spin purity of the open-shell orbital subspace to the extent allowed by the original reference, thus facilitating a safe truncation of the virtual space. The performance of OSFNO is benchmarked for selected diradicals and triradicals.

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

Sparsity of the electron repulsion integral tensor using different localized virtual orbital representations in local second-order Møller–Plesset theory

Utilizing localized orbitals, local correlation theory can reduce the unphysically high system-size scaling of post-Hartree–Fock (post-HF) methods to linear scaling in insulating molecules. The sparsity of the four-index electron repulsion integral (ERI) tensor is central to achieving this reduction. For second-order Møller–Plesset theory (MP2), one of the simplest post-HF methods, only the (ia|jb) ERIs are needed, coupling occupied orbitals i, j and virtuals a, b. In this paper, we compare the numerical sparsity (called the “ragged list”) and two other approaches revealing the low-rank sparsity of the ERI. The ragged list requires only one set of (localized) virtual orbitals, and we find that the orthogonal valence virtual-hard virtual set of virtuals originally proposed by Subotnik et al. gives the sparsest ERI tensor. To further compress the ERI tensor, the pair natural orbital (PNO) type representation uses different sets of virtual orbitals for different occupied orbital pairs, while the occupied-specific virtual (OSV) approach uses different virtuals for each occupied orbital. Here, our results indicate that while the low-rank PNO representation achieves significant rank reduction, it also requires more memory than the ragged list. The OSV approach requires similar memory to that of the ragged list, but it involves greater algorithmic complexity. An approximation (called the “fixed sparsity pattern”) for solving the local MP2 equations using the numerically sparse ERI tensor is proposed and tested to be sufficiently accurate and to have highly controllable error. A low-scaling local MP2 algorithm based on the ragged list and the fixed sparsity pattern is therefore promising.

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

Effective rationality for local unitary invariants of mixed states of two qubits

Abstract We calculate the field of rational local unitary invariants for mixed states of two qubits, by employing methods from algebraic geometry. We prove that this field is rational (i.e. purely transcendental), and that it is generated by nine algebraically independent polynomial invariants. We do so by constructing a relative section, in the sense of invariant theory, whose Weyl group is a finite abelian group. From this construction, we are able to give explicit expressions for the generating invariants in terms of the Bloch matrix representation of mixed states of two qubits. We also prove similar rationality results for the local unitary invariants of symmetrically mixed states of two qubits. We also provide a sketch of how to generalize our results to the case of an arbitrary number of qubits. Our results apply to both complex-valued and real-valued invariants.

Physics↗
DOE OSTI2026

Variance reduction in lattice QCD observables via normalizing flows

Normalizing flows can be used to construct unbiased, reduced-variance estimators for lattice field theory observables that are defined by a derivative with respect to action parameters. This work implements the approach for observables involving gluonic operator insertions in the SU(3) Yang-Mills theory and two-flavor QCD in four space-time dimensions. Variance reduction by factors of 10–60 is achieved in glueball correlation functions and in gluonic matrix elements related to hadron structure, with demonstrated computational advantages. The observed variance reduction is found to be approximately independent of the lattice volume, so volume transfer can be utilized to minimize training costs.

Abbott, Ryan [Columbia U.; MIT, Cambridge, CTP; IA↗
DOE OSTI2020

Superdiffusion from Emergent Classical Solitons in Quantum Spin Chains

Finite-temperature spin transport in the quantum Heisenberg spin chain is known to be superdiffusive, and has been conjectured to lie in the Kardar-Parisi-Zhang (KPZ) universality class. Using a kinetic theory of transport, we compute the KPZ coupling strength for the Heisenberg chain as a function of temperature, directly from microscopics; the results agree well with density-matrix renormalization group simulations. Here, we establish a rigorous quantum-classical correspondence between the “giant quasiparticles” that govern superdiffusion and solitons in the classical continuous Landau-Lifshitz ferromagnet. We conclude that KPZ universality has the same origin in classical and quantum integrable isotropic magnets: a finite-temperature gas of low-energy classical solitons.

1-dimensional spin chains↗
DOE OSTI2022

Advances in Nuclear Data Evaluation Theory for NCSP [Slides]

The generalized form of the Bayes’ Theorem improves UQ for NCSP of differential and/or integral data evaluations and of applications. It is being implemented in an API for use by SCALE, SAMMY, etc. Extended R-matrix formalism improves evaluations of RRR including direct reactions, direct capture, and the Reich-Moore approx., doorway state reactions with applications to RPI data on lead isotopes. Extended R-Matrix formalism will be implemented into SAMMY post modernization.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1983

Ultrasonic wave propagation in two-phase media: Spherical inclusions

The scattering theory, recently developed via the extended method of equivalent inclusion, is used to study the propagation of time-harmonic waves in two-phase media of elastic matrix with randomly distributed elastic spherical inclusion materials. The elastic moduli and mass density of the composite medium are determined as functions of frequencies when given properties and concentration of the spheres and the matrix. Velocity and attenuation of ultrasonic waves in two-phase media are determined for cases of distributed spheres and localized damage. An averaging theorem that requires the equivalence of the strain energy and the kinetic energy between the effective medium and the original matrix with spherical inhomogeneities is employed to derive the effective moduli and mass density. The functional dependency of these quantities upon frequencies and concentration provides a method of data analysis in ultrasonic evaluation of material properties. Numerical results or moduli, velocity and/or attenuation as functions of concentration of inclusion material, or porosity, are graphically displayed.

Fu, L. S.↗
NASA NTRS1985

Thermal expansion of selected graphite reinforced polyimide-, epoxy-, and glass-matrix composite

The thermal expansion of three epoxy-matrix composites, a polyimide-matrix composite and a borosilicate glass-matrix composite, each reinforced with continuous carbon fibers, has been measured and compared. The expansion of a composite with a rubber toughened epoxy-matrix and P75S carbon fibers was very different from the expansion of two different single phase epoxy-matrix composites with P75S fibers although all three had the same stacking sequence. Reasonable agreement was obtained between measured thermal-expansion data and results from classical laminate theory. The thermal expansion of a material may change markedly as a result of thermal cycling. Microdamage, induced by 250 cycles between -156 C and 121 C in the graphite/polyimide laminate, caused a 53 percent decrease in the coefficient of thermal expansion. The thermal expansion of the graphite/glass laminate was not changed by 100 thermal cycles from -129 C to 38 C; however, a residual strain of about 10 x 10 to the minus 6 power was measured for the laminate tested.

Tompkins, S. S.↗
NASA NTRS1984

Ultrasonic wave propagation in two-phase media - Spherical inclusions

The scattering theory, recently developed via the extended method of equivalent inclusion, is used to study the propagation of time-harmonic waves in two-phase media of elastic matrix with randomly distributed elastic spherical inclusion materials. The elastic moduli and mass density of the composite medium are determined as functions of frequencies when given properties and concentration of the spheres and the matrix. Velocity and attenuation of ultrasonic waves in two-phase media are determined for cases of distributed spheres and localized damage. An averaging theorem that requires the equivalence of the strain energy and the kinetic energy between the effective medium and the original matrix with spherical inhomogeneities is employed to derive the effective moduli and mass density. The functional dependency of these quantities upon frequencies and concentration provides a method of data analysis in ultrasonic evaluation of material properties. Numerical results or moduli, velocity and/or attenuation as functions of concentration of inclusion material, or porosity, are graphically displayed.

Fu, L. S.↗
NASA NTRS1985

Accretionary rims on inclusions in the Allende meteorite

The origin and composition of the rim sequence on the refractory inclusion in the Allende meteorite are studied. The different textures, mineralogy, and mineral-chemistry of the four layers of the rim are described. The layers are composed of: pyroxene, needles, olivine, hedenbergite, and andradite. Tables of the element and chemical compositions of the layers are presented. The data reveals that: (1) the layers are highly porous masses of euhedral crystals with no intergrowth; (2) layers contain highly disequilibrium mineral assemblages; and (3) the thickness of the layers varies with the underlying topography. These results support the theory that rim structures are accretionary aggregates formed from accretion of independently grown particles onto the surface of inclusions. The formation of the grains in the layers and matrix from nebular condensates is studied.

Macpherson, G. J.↗
NASA NTRS1993

Frequency domain system identification methods - Matrix fraction description approach

This paper presents the use of matrix fraction descriptions for least-squares curve fitting of the frequency spectra to compute two matrix polynomials. The matrix polynomials are intermediate step to obtain a linearized representation of the experimental transfer function. Two approaches are presented: first, the matrix polynomials are identified using an estimated transfer function; second, the matrix polynomials are identified directly from the cross/auto spectra of the input and output signals. A set of Markov parameters are computed from the polynomials and subsequently realization theory is used to recover a minimum order state space model. Unevenly spaced frequency response functions may be used. Results from a simple numerical example and an experiment are discussed to highlight some of the important aspect of the algorithm.

Horta, Luca G.↗
NASA NTRS1995

Comparison of the Compressive Strengths for Stitched and Toughened Composite Systems

The compression strength of a stitched and a toughened matrix graphite/epoxy composite was determined and compared to a baseline unstitched untoughened composite. Two different layups with a variety of test lengths were tested under both ambient and hot/wet conditions. No significant difference in strength was seen for the different materials when the gage lengths of the specimens were long enough to lead to a buckling failure. For shorter specimens, a 30% reduction in strength from the baseline was seen due to stitching for both a 48-ply quasi-isotropic and a (0/45/0/-45/90/-45/0/45/0)s laminate. Analysis of the results suggested that the decrease in strength was due to increased fiber misalignment due to the stitches. An observed increasing strength with decreasing gage length, which was seen for all materials, was explained with a size effect model. The model assumed a random distribution of flaws (misaligned fibers). The toughened material showed a small increase in strength over the baseline material for both laminates presumably due to the compensating effects of a more compliant matrix and straighter fibers in the toughened material. The hot/wet strength of the stitched and baseline material fell 30% below their ambient strengths for shorter, non-buckling specimen, while the strength of the toughened matrix material only fell 20%. Video images of the failing specimen were recorded and showed local failures prior to global collapse of the specimen. These images support the theory of a random distribution of flaws controlling composite failure. Failed specimen appearance however, seems to be a misleading indication of the cause of failure.

Reeder, James R.↗
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↗
NASA NTRS2001

Tunable Optical Properties of Metal Nanoparticle Sol-Gel Composites

We demonstrate that the linear and non-linear optical properties of sol-gels containing metal nanoparticles are highly tunable with porosity. Moreover, we extend the technique of immersion spectroscopy to inhomogeneous hosts, such as aerogels, and determine rigorous bounds for the average fractional composition of each component, i.e., the porosity of the aerogel, or equivalently, for these materials, the catalytic dispersion. Sol-gels containing noble metal nanoparticles were fabricated and a significant blue-shift in the surface plasmon resonance (SPR) was observed upon formation of an aerogel, as a result of the decrease in the dielectric constant of the matrix upon supercritical extraction of the solvent. However, as a result of chemical interface damping and aggregation this blue-shift does not strictly obey standard effective medium theories. Mitigation of these complications is achieved by avoiding the use of alcohol and by annealing the samples in a reducing atmosphere.

Smith, David D.↗
NASA NTRS2002

Fracture Mechanisms For SiC Fibers And SiC/SiC Composites Under Stress-Rupture Conditions at High Temperatures

The successful application of SiC/SiC ceramic matrix composites as high-temperature structural materials depends strongly on maximizing the fracture or rupture life of the load-bearing fiber and matrix constituents. Using high-temperature data measured under stress-rupture test conditions, this study examines in a mechanistic manner the effects of various intrinsic and extrinsic factors on the creep and fracture behavior of a variety of SiC fiber types. It is shown that although some fiber types fracture during a large primary creep stage, the fiber creep rate just prior to fracture plays a key role in determining fiber rupture time (Monkman-Grant theory). If it is assumed that SiC matrices rupture in a similar manner as fibers with the same microstructures, one can develop simple mechanistic models to analyze and optimize the stress-rupture behavior of SiC/SiC composites for applied stresses that are initially below matrix cracking.

DiCarlo, James A.↗
NASA NTRS2001

Fully Coupled Micro/Macro Deformation, Damage, and Failure Prediction for SiC/Ti-15-3 Laminates

The deformation, failure, and low cycle fatigue life of SCS-6/Ti-15-3 composites are predicted using a coupled deformation and damage approach in the context of the analytical generalized method of cells (GMC) micromechanics model. The local effects of inelastic deformation, fiber breakage, fiber-matrix interfacial debonding, and fatigue damage are included as sub-models that operate on the micro scale for the individual composite phases. For the laminate analysis, lamination theory is employed as the global or structural scale model, while GMC is embedded to operate on the meso scale to simulate the behavior of the composite material within each laminate layer. While the analysis approach is quite complex and multifaceted, it is shown, through comparison with experimental data, to be quite accurate and realistic while remaining extremely efficient.

Bednarcyk, Brett A.↗
DOE OSTI2021

AdS bulk locality from sharp CFT bounds

It is a long-standing conjecture that any CFT with a large central charge and a large gap Δ gap in the spectrum of higher-spin single-trace operators must be dual to a local effective field theory in AdS. We prove a sharp form of this conjecture by deriving numerical bounds on bulk Wilson coefficients in terms of Δ gap using the conformal bootstrap. Our bounds exhibit the scaling in Δ gap expected from dimensional analysis in the bulk. Our main tools are dispersive sum rules that provide a dictionary between CFT dispersion relations and S-matrix dispersion relations in appropriate limits. This dictionary allows us to apply recently-developed flat-space methods to construct positive CFT functionals. We show how AdS 4 naturally resolves the infrared divergences present in 4D flat-space bounds. Our results imply the validity of twice-subtracted dispersion relations for any S-matrix arising from the flat-space limit of AdS/CFT.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗