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 289 records · Page 16

In situ high-temperature 3D imaging of the damage evolution in a SiC nuclear fuel cladding material

Silicon carbide (SiC)-based nuclear fission fuel rod cladding has been considered as one of the possible designs for accident tolerant fuels. It is in the form of a SiC fibre reinforced SiC matrix composite tube (SiC f -SiC m ) with monolithic SiC outer and/or inner coating layers. This study focuses on the deformation and fracture processes in this material using in situ X-ray micro-computed tomography (XCT) at room temperature (RT) and 1200 °C in an inert gas environment in a C-ring compression loading configuration. Prior to testing, local properties and residual stresses were characterised using nanoindentation and Raman spectroscopy since they can impact the mechanical behaviour of the material. The 3D strain distribution, crack formation and propagation processes including the toughening mechanisms (e.g., crack deflection, micro-cracking, crack bridging and bifurcation) are investigated in the coating and underlying composites at RT and 1200 °C. There is no particular sequence which toughening mechanism occurs first – this is very different from the conventional toughening theory in ceramic-matrix composites under uniaxial tension loading. Indeed, no evidence of fibre pull-out or fibre fracture was observed in this SiC f -SiC m nuclear cladding material in the current C-ring compression configuration. The correlation between local measurements and bulk mechanical behaviour are discussed.

3D imaging↗

Electronic structure, magnetic properties, and pairing tendencies of the copper-based honeycomb lattice Na 2 Cu 2 TeO 6

Spin-1/2 chains with alternating antiferromagnetic (AFM) and ferromagnetic (FM) couplings have attracted considerable interest due to the topological character of their spin excitations. Here, using density functional theory and density-matrix renormalization-group (DMRG) methods, we have systematically studied the dimerized chain system Na 2 Cu 2 TeO 6 with a d 9 electronic configuration. Near the Fermi level, in the nonmagnetic phase the dominant states are mainly contributed by the Cu 3d x 2 -y 2 orbitals highly hybridized with the O 2p orbitals, leading to an “effective” single-orbital low-energy model. By calculating the relevant hoping amplitudes, we explain the size and sign of the exchange interactions in Na 2 Cu 2 TeO 6 . In addition, a single-orbital Hubbard model is constructed for this dimerized chain system where the quantum fluctuations are taken into account. Both AFM and FM couplings (leading to an ↑|-|↓-|↓-↑| state) along the chain were found in our DMRG and Lanczos calculations, in agreement with density functional theory and neutron-scattering results. The hole pairing binding energy ΔE is predicted to be negative at Hubbard U~ 11eV, suggesting incipient pairing tendencies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Quantum Monte Carlo calculation of 𝛿 NS in 10 C using an effective field theory approach

Here, we compute radiative corrections to the superallowed 𝛽 decay of 10 C in an effective field theory approach using nuclear matrix elements obtained from quantum Monte Carlo calculations. These corrections are an important ingredient in the extraction of the Cabibbo-Kobayashi-Masakawa quark mixing matrix element 𝑉 𝑢⁢𝑑 , and the role of this work is to illuminate the uncertainties arising from nuclear structure. Our results provide good agreement with both the traditional extraction of 𝑉 𝑢⁢𝑑 , as well as with a more recent evaluation performed using the no-core shell model and a dispersion formalism. The dominant uncertainty in this approach is the presence of two unknown low-energy constants that enter into the relevant nuclear matrix elements. Future determinations of these low-energy constants—either from QCD or modeling them with two nucleon amplitudes—would improve the precision of the extraction in this formalism.

beta decay↗

Tritium β decay and proton-proton fusion in pionless effective field theory

The Gamow-Teller and Fermi matrix elements, 〈 GT 〉 and 〈 F 〉 , respectively, for tritium β decay are calculated to next-to-leading order (NLO) in pionless effective field theory in the absence of Coulomb interactions and isospin violation giving the leading order predictions 〈 GT 〉 0 = 0.9807 and 〈 F 〉 0 = 1 . Using an experimentally determined value for the tritium β decay GT matrix element, the two-body axial current low energy constant is fixed at NLO yielding L 1 , A = 6.01 ± 2.08 fm 3 at the renormalization scale of the physical pion mass, which agrees with predictions based on naive dimensional analysis. The impact of L 1 , A on proton-proton fusion is also discussed. Finally, the consequences of Wigner-SU(4) spin-isospin symmetry are considered for the Gamow-Teller matrix element. Published by the American Physical Society 2024

Physics↗

An Improved Fracture Mechanics-Informed Multiscale Thermomechanical Damage Model for Ceramic Matrix Composites

This paper extends recent work done by the authors in modeling length scale-dependent damage behavior of ceramic matrix composites (CMCs) to include effects of local anisotropy introduced by matrix cracking. This model captures scale-dependent damage initiation and propagation behavior of the brittle matrix by employing internal state variable (ISV) theory within a multiscale modeling framework to obtain damaged matrix stress/strain constitutive relationships at each length scale. The damage ISV captures the effects of matrix cracking and growth by using fracture mechanics and the self-consistent scheme to determine the reduced stiffness of the cracked matrix. Matrix cracks, which activate when stress intensity factors near manufacturing induced cavities exceed the fracture toughness of the material, are assumed to be transversely isotropic in the plane of the crack, and matrix anisotropy occurs when the damaged stiffness tensor is rotated from the crack plane to the global axes. The crack progression and temporal evolution of the damage ISV are governed by fracture mechanics and crack growth kinetics. The model effectively captures first matrix cracking, which is the first significant deviation from linear elasticity. The nonlinear predictive capabilities of the material model are demonstrated for monolithic silicon carbide (SiC) and a 2D woven five-harness satin (5HS) carbon fiber SiC matrix (C/SiC) CMC.

Skinner, Travis↗

Analysis of the progressive failure of brittle matrix composites

This report investigates two of the most common modes of localized failures, namely, periodic fiber-bridged matrix cracks and transverse matrix cracks. A modification of Daniels' bundle theory is combined with Weibull's weakest link theory to model the statistical distribution of the periodic matrix cracking strength for an individual layer. Results of the model predictions are compared with experimental data from the open literature. Extensions to the model are made to account for possible imperfections within the layer (i.e., nonuniform fiber lengths, irregular crack spacing, and degraded in-situ fiber properties), and the results of these studies are presented. A generalized shear-lag analysis is derived which is capable of modeling the development of transverse matrix cracks in material systems having a general multilayer configuration and under states of full in-plane load. A method for computing the effective elastic properties for the damaged layer at the global level is detailed based upon the solution for the effects of the damage at the local level. This methodology is general in nature and is therefore also applicable to (0(sub m)/90(sub n))(sub s) systems. The characteristic stress-strain response for more general cases is shown to be qualitatively correct (experimental data is not available for a quantitative evaluation), and the damage evolution is recorded in terms of the matrix crack density as a function of the applied strain. Probabilistic effects are introduced to account for the statistical nature of the material strengths, thus allowing cumulative distribution curves for the probability of failure to be generated for each of the example laminates. Additionally, Oh and Finney's classic work on fracture location in brittle materials is extended and combined with the shear-lag analysis. The result is an analytical form for predicting the probability density function for the location of the next transverse crack occurrence within a crack bounded region. The results of this study verified qualitatively the validity of assuming a uniform crack spacing (as was done in the shear-lag model).

Thomas, David J.↗

Reliability analysis of ceramic matrix composite laminates

At a macroscopic level, a composite lamina may be considered as a homogeneous orthotropic solid whose directional strengths are random variables. Incorporation of these random variable strengths into failure models, either interactive or non-interactive, allows for the evaluation of the lamina reliability under a given stress state. Using a non-interactive criterion for demonstration purposes, laminate reliabilities are calculated assuming previously established load sharing rules for the redistribution of load as the failure of laminae occur. The matrix cracking predicted by ACK theory is modeled to allow a loss of stiffness in the fiber direction. The subsequent failure in the fiber direction is controlled by a modified bundle theory. Results using this modified bundle model are compared with previous models which did not permit separate consideration of matrix cracking, as well as to results obtained from experimental data.

Thomas, David J.↗

Residual thermal and moisture influences on the strain energy release rate analysis of local delaminations from matrix cracks

An analysis utilizing laminated plate theory is developed to calculate the strain energy release rate associated with local delaminations originating at off-axis, single ply, matrix cracks in laminates subjected to uniaxial loads. The analysis includes the contribution of residual thermal and moisture stresses to the strain energy released. Examples are calculated for the strain energy release rate associated with local delaminations originating at 90 degrees and angle-ply (non-90 degrees) matrix ply cracks in glass epoxy and graphite epoxy laminates. The solution developed may be used to assess the relative contribution of mechanical, residual thermal, and moisture stresses on the strain energy release rate for local delamination for a variety of layups and materials.

Obrien, T. K.↗

One-Body Properties and Their Perturbative Accuracy with Aufbau Suppressed Coupled Cluster Theory

In this work we derived and implemented the calculation of the one-body reduced density matrix for Aufbau suppressed coupled cluster theory, from which excited state natural orbitals and one-body properties, like atomic populations and dipole moments, are obtained. We utilized the natural orbitals to refine the ASCC solution for simple valence and Rydberg systems, exploring the process of repeatedly solving the ASCC equations in successive natural orbital bases to achieve independence from the starting molecular orbitals. For dipole moments in small molecules where high-level comparison data is available, we find that the accuracy of ASCC essentially matches that of linear response and equation-of-motion coupled cluster as long as care is taken to preserve the response's perturbative completeness.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Lattice QCD Inputs for nuclear double beta decay

Second order β -decay processes with and without neutrinos in the final state are key probes of nuclear physics and of the nature of neutrinos. Neutrinoful double- β decay is the rarest Standard Model process that has been observed and provides a unique test of the understanding of weak nuclear interactions. Observation of neutrinoless double- β decay would reveal that neutrinos are Majorana fermions and that lepton number conservation is violated in nature. While significant progress has been made in phenomenological approaches to understanding these processes, establishing a connection between these processes and the physics of the Standard Model and beyond is a critical task as it will provide input into the design and interpretation of future experiments. The strong-interaction contributions to double- β decay processes are non-perturbative and can only be addressed systematically through a combination of lattice Quantum Chromoodynamics (LQCD) and nuclear many-body calculations. In this review, current efforts to establish the LQCD connection are discussed for both neutrinoful and neutrinoless double- β decay. LQCD calculations of the hadronic contributions to the neutrinoful process $nn → ppe^-e^-\bar{v}_e\bar{v}_e$ and to various neutrinoless pionic transitions are reviewed, and the connections of these calculations to the phenomenology of double- β decay through the use of effective field theory (EFTs) is highlighted. At present, LQCD calculations are limited to small nuclear systems, and to pionic subsystems, and require matching to appropriate EFTs to have direct phenomenological impact. However, these calculations have already revealed qualitatively that there are terms in the EFTs that can only be constrained from double- β decay processes themselves or using inputs from LQCD. Finally, future prospects for direct calculations in larger nuclei are also discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Probing exciton/exciton interactions with entangled photons: Theory

Quantum entangled photons provide a sensitive probe of many-body interactions and offer a unique experimental portal for quantifying many-body correlations in a material system. In this paper, we present a theoretical demonstration of how photon–photon entanglement can be generated via interactions between coupled qubits. Here, we develop a model for the scattering of an entangled pair of photons from a molecular dimer. We develop a diagrammatic theory for the scattering matrix and show that one can correlate the von Neumann entropy of the outgoing bi-photon wave function with exciton exchange and repulsion interactions. Finally, we conclude by discussing possible experimental scenarios for realizing these ideas.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Recent developments in the PySCF program package

PySCF is a Python-based general-purpose electronic structure platform that supports first-principles simulations of molecules and solids as well as accelerates the development of new methodology and complex computational workflows. Here, we explain the design and philosophy behind PySCF that enables it to meet these twin objectives. With several case studies, we show how users can easily implement their own methods using PySCF as a development environment. We then summarize the capabilities of PySCF for molecular and solid-state simulations. Finally, we describe the growing ecosystem of projects that use PySCF across the domains of quantum chemistry, materials science, machine learning, and quantum information science.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the use of stereodynamical effects to control cold chemical reactions: The H + D 2 ↔ D + HD case study

Here, quantum calculations are reported for the stereodynamic control of the H + D 2 ↔ D + HD chemical reaction in the energy range of 1–50 K. Stereodynamic control is achieved by a formalism similar to that reported by Perreault et al. [Nat. Chem. 10, 561 (2018)] in recent experimental works in which the alignment of the molecular bond axis relative to the incident relative velocity is controlled by selective preparations of the molecule in a specific or superposition of magnetic projection quantum numbers of the initial molecular rotational level. The approach presented here generalizes the experimental scheme of Perreault et al. and offers additional degree of control through various experimental preparations of the molecular alignment angle. Illustrative results presented for the H + D 2 and D + HD reactions show significant control with the possibility of turning the reaction completely on or off with the appropriate stereodynamic preparation of the molecular state. Various scenarios for maximizing and minimizing the reaction outcomes are identified with the selective preparation of molecular rotational states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ultrafast charge dynamics and photoluminescence in bilayer MoS 2

Our examination of the interplay of ultrafast charge dynamics and electron–phonon interaction in the AA' stacked bilayer MoS 2 provides a microscopic basis for understanding the features (two peaks) in the emission spectrum. We demonstrate that while the initial accumulation of excited charge occurs at and near the Q point of the two-dimensional Brillioun zone, emission takes place predominantly through two pathways: direct charge recombination at the K point and indirect phonon-assisted recombination of electrons at the K valley and holes at the Γ hill of the Brillouin zone. Analysis of the wave vector dependencies of the electron–phonon interaction traces the higher energy peak to phonon-assisted relaxation of the excited electrons from the Q to the K valley in the conduction band. Our results thus reveal the importance of ultrafast charge dynamics in understanding photoemissive properties of a few-layer transition-metal dichalcogenide. These calculations are based on time dependent density functional theory in the density matrix formulation.

Materials Science↗

Spatially anisotropic 𝑆=1 square-lattice antiferromagnet with single-ion anisotropy realized in a Ni(II) pyrazine-𝑛,𝑛′-dioxide coordination polymer

The Ni(NCS) 2 (pyzdo) 2 coordination polymer is found to be an S = 1 spatially anisotropic square lattice with easy-axis single-ion anisotropy. This conclusion is based upon considering in concert the experimental probes x-ray diffraction, magnetic susceptibility, magnetic-field-dependent heat capacity, muon-spin relaxation, neutron diffraction, neutron spectroscopy, and pulsed-field magnetization. Long-range antiferromagnetic (AFM) order develops at T N = 18.5 K. Although the samples are polycrystalline, there is an observable spin-flop transition and saturation of the magnetization at ≈80 T. Linear spin-wave theory yields spatially anisotropic exchanges within an AFM square lattice, J x =0.235 meV, J y =2.014 meV, and an easy-axis single-ion anisotropy D=−1.622 meV (after renormalization). The anisotropy of the exchanges is supported by density functional theory.

2-dimensional systems↗

Bypassing the Energy Functional in Density Functional Theory: Direct Calculation of Electronic Energies from Conditional Probability Densities

Density functional calculations can fail for want of an accurate exchange-correlation approximation. The energy can instead be extracted from a sequence of density functional calculations of conditional probabilities (CP DFT). Simple CP approximations yield usefully accurate results for two-electron ions, the hydrogen dimer, and the uniform gas at all temperatures. CP DFT has no self-interaction error for one electron, and correctly dissociates H2, both major challenges. For warm dense matter, classical CP DFT calculations can overcome the convergence problems of Kohn-Sham DFT.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Unified description of the coupled-channels and statistical Hauser-Feshbach nuclear reaction theories for low energy neutron incident reactions

We incorporate the coupled-channels optical model into the statistical Hauser-Feshbach nuclear reaction theory, where the scattering matrix is diagonalized by performing the Engelbrecht-Weidenmüller transformation. This technique has been implemented in the coupled-channels optical model code ECIS by J. Raynal, and we extend this method so that all the open channels in a nucleon-induced reaction on a deformed nucleus can be calculated consistently.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗