Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix Models”

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,783 records · Page 99
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Measurements and interpretations of W ± Z production cross-sections in pp collisions at $\sqrt{s}=13$ TeV with the ATLAS detector

Measurements of integrated and differential cross-sections for W ± Z production in proton-proton collisions are presented. The data collected by the ATLAS detector at the Large Hadron Collider from 2015 to 2018 at a centre-of-mass energy of $\sqrt{s}=13$ TeV are used, corresponding to an integrated luminosity of 140 fb −1 . The W ± Z candidate events are reconstructed using leptonic decay modes of the gauge bosons into electrons or muons. The integrated cross-section per lepton flavour for the production of W ± Z is measured in the detector fiducial region with a relative precision of 4%. The measured value is compared with the Standard Model prediction at a precision of up to next-to-next-to-leading-order in QCD and next-to-leading-order in electroweak. Cross-sections for W + Z and W − Z production and their ratio are presented. The W ± Z production is also measured differentially as functions of various kinematic variables, including new observables sensitive to CP-violation effects. All measurements are compared with state-of-the-art Standard Model predictions from fixed-order calculations or Monte Carlo generators based on next-to-leading-order matrix elements interfaced with parton showers. An effective field theory interpretation of the measurements is performed, considering both CP-conserving and CP-violating dimension-6 operators modifying the W ± Z production. In the absence of observed deviations from the Standard Model, limits on CP-conserving Wilson coefficients are extracted using the transverse mass of the W ± Z system. For CP-violating coefficients a machine learning approach is used to construct an observable with enhanced sensitivity to CP-violation effects.

hadron-hadron scattering↗
DOE OSTI2020

Mixed quantum–classical approach to model non-adiabatic electron–nuclear dynamics: Detailed balance and improved surface hopping method

We develop a density matrix formalism to describe coupled electron–nuclear dynamics. To this end, we introduce an effective Hamiltonian formalism that describes electronic transitions and small (quantum) nuclear fluctuations along a classical trajectory of the nuclei. Using this Hamiltonian, we derive equations of motion for the electronic occupation numbers and for the nuclear coordinates and momenta. We show that, in the limit, when the number of nuclear degrees of freedom coupled to a given electronic transition is sufficiently high (i.e., the strong decoherence limit), the equations of motion for the electronic occupation numbers become Markovian. Furthermore, the transition rates in these (rate) equations are asymmetric with respect to the lower-to-higher energy transitions and vice versa. In thermal equilibrium, such asymmetry corresponds to the detailed balance condition. We also study the equations for the electronic occupations in the non-Markovian regime and develop a surface hopping algorithm based on our formalism. To treat the decoherence effects, we introduce additional “virtual” nuclear wave packets whose interference with the “real” (physical) wave packets leads to the reduction in coupling between the electronic states (i.e., decoherence) as well as to the phase shifts that improve the accuracy of the numerical approach. Remarkably, the same phase shifts lead to the detailed balance condition in the strong decoherence limit.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2013

Estimation of Sorption Behavior of Europium(III) Using Biotite Flakes - 13272

The interaction of biotite and Eu(III) (europium (III)) was examined by using secondary ion-microprobe mass spectrometer (SIMS), fluorescence emission spectrum and decay behavior of fluorescence emission spectrum in addition to the time-changes of Eu(III) and potassium ions concentrations in a solution, using the flake form samples. The results of SIMS showed that the intensity of Eu was gradually decreasing with depth, while the intensity of Eu in the case shaken for 30 days exceeded that in the case for 1 day. Furthermore, the spatial distribution of Eu(III) and potassium ions in the flake of biotite suggested that Eu ions diffuse mainly from the edges of biotite flake, while Eu ions can slightly diffuse through some small cracks existing on the flake surface far from the edges. Besides, the elution amount of potassium from the biotite flakes into a solution was proportional to the sorption amount of Eu(III). The changes nearly revealed ion exchange between these ions, while muscovite flake sample did not show such ion exchange reaction. In addition, from the time-change of Eu(III) concentration, an apparent diffusion coefficient was estimated to be 8.0x10{sup -12} m{sup 2}/s, by using two-dimensional diffusion model coupled with a film between the solid phase and the liquid phase. Furthermore, the fluorescent intensity decreased with the shaking (contacting) time. This means that Eu(III) gradually diffuses into the inside of biotite edges of the biotite flakes, after the sorption of Eu(III) in the edges. This tendency was observed also in the powder samples. The observed fluorescence decay (at 592 nm in wave length) showed almost similar curve in any samples, indicating a certain sorption form of Eu(III) onto the edges of the biotite flakes. These results mentioned above suggest that the diffusion processes through internal layer in biotite mainly control the sorption behavior of multivalent ions. Such diffusion processes affect the retardation-effects on fracture surfaces in the rock matrix, depending on the fluid flow velocity of groundwater. That is, a more reliable model considering the mass transfer in the internal layer of biotite may be required to estimate the sorption behavior of RNs with biotite which controls the whole sorption behavior of granite. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2012

Optical anisotropy of InAs/GaSb broken-gap quantum wells

We investigate in detail the optical anisotropy of absorption of linearly polarized light in InAs/GaSb quantum wells grown on GaSb along the [001] direction, which can be used as an active region of different laser structures. The energy level positions, the wave functions, the optical matrix elements, and the absorption coefficients are calculated using the eight-band k {center_dot} p model and the Burt-Foreman envelope function theory. In these calculations, the Schroedinger and Poisson equations are solved self-consistently taking the lattice-mismatched strain into account. We find that a realistic Hamiltonian, which has the C{sub 2v} symmetry, results in considerable anisotropy of optical matrix elements for different directions of light polarization and different directions of the initial-state in-plane wave vector, including low-symmetry directions. We trace how the optical matrix elements and absorption are modified when spin-orbit interaction and important symmetry breaking mechanisms are taken into account (structural inversion asymmetry, bulk inversion asymmetry, and interface Hamiltonian). These mechanisms result in an almost 100% anisotropy of the absorption coefficients as the light polarization vector rotates in the plane of the structure and in a plane normal to the interfaces.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Fully implicit crystal plasticity models representing orientations with modified Rodrigues parameters

Here, this work describes a crystal plasticity formulation combining several mathematical, numerical, and implementation choices to produce a highly efficient model. Specifically, the key choices in the implementation are (1) representing orientations with modified Rodrigues parameters, (2) implementing a fully coupled implicit time integration for the elastic stretch, the crystal orientations, and the model internal variables, (3) implementing the model in the NEML2 constitutive modeling framework, based on PyTorch, to vectorize the calculations and port the computation to GPUs and other hardware accelerators, and (4) an exact implementation of the consistent tangent matrix, even for arbitrary coupling to other field variables beyond the displacements, like temperature, neutron fluence, etc. The first two features of the model are, to our knowledge, novel. The paper considers each of these choices individually as well as the final model as a whole. This includes a full description of modified Rodrigues parameters, their advantages over other representations of orientations, the mathematical formulae and tools required to implement a model with modified Rodrigues parameters, and a detailed description of the geometry of the space of modified Rodrigues parameters (in an appendix). It also includes a description of a fully implicit time integration scheme for the orientations and the advantages in representing orientations with modified Rodrigues parameters in implementing such a model. The work then assess, via numerical examples, the advantages of fully coupled implicit time integration versus more common decoupled and explicit time integration schemes. These studies demonstrate the computational advantages of fully coupled integration versus other time integration algorithms, though the performance of the competing models depends on the complexity of the underlying single crystal model. The study concludes by demonstrating that the choice of time integration method affects the sharpness of the predicted texture, with explicit methods for integrating the orientations overestimating texture sharpness and implicit methods underestimating texture sharpness.

Crystal plasticity↗
DOE OSTI2019

Block excitonic condensate at n = 3.5 in a spin-orbit coupled t 2 g multiorbital Hubbard model

Theoretical studies recently predicted the condensation of spin-orbit excitons at momentum $q={\pi}$ in ${t}_{2g}^{4}$ spin-orbit coupled three-orbital Hubbard models at electronic density $n=4$. In parallel, experiments involving iridates with noninteger valence states for the Ir ions are starting to attract considerable attention. Using the density matrix renormalization group technique we present evidence for the existence of an excitonic condensate at $n=3.5$ in a one-dimensional Hubbard model with a degenerate ${t}_{2g}$ sector, when in the presence of spin-orbit coupling. At intermediate Hubbard $U$ and spin-orbit ${\lambda}$ couplings, we found an excitonic condensate at the unexpected momentum $q={\pi}/2$ involving ${j}_{\mathrm{eff}}=3/2,m=\pm{}1/2$, and ${j}_{\mathrm{eff}}=1/2,m=\pm{}1/2$ bands in the triplet channel, coexisting with an also unexpected block magnetic order. We also present the entire ${\lambda}$ vs $U$ phase diagram, at a fixed and robust Hund coupling. Interestingly, this new “block excitonic phase” is present even at large values of ${\lambda}$, unlike the $n=4$ excitonic phase discussed before. Our computational study helps to understand and predict the possible magnetic phases of materials with ${d}^{3.5}$ valence and robust spin-orbit coupling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

3D printed lignin/polymer composite with enhanced mechanical and anti-thermal-aging performance

Lignin is the most abundant natural aromatic polymer globally but is still underutilized as a renewable material, even with its versatile properties. One approach to using lignin is to incorporate it in polymer composites, but this application is often limited by the poor mechanical performance of the resultant composite due to the poor interfacial adhesion. Following the structure–property relationship, stronger interactions were designed to be enhanced by alternating the functional groups of lignin. Here, this study applied a demethylation method to hardwood lignin (Ori-Lig) to introduce more phenolic hydroxyl groups. Research studies with rheology behavior and molecular simulations demonstrated that an increased phenolic hydroxyl content could improve the adhesion between the modified lignin (OH-Lig) and the polymer matrix at the interface. The tensile performance of the samples from the fused depositional modeling (FDM) 3D printing technique was also improved due to the improved interfacial adhesion. Specifically, by adding 10 wt% of the modified lignin, the tensile strength of the 3D printed samples could reach 46.1 MPa (40.2 MPa of Polyamide 12, PA12), and Young’s Modulus could be improved to 1.73 GPa (1.48 GPa of PA12). OH-Lig also improved anti-aging performance so that the tensile strength of the OH-Lig composite after thermo-aging (100 h at 140 °C) remained ~48MPa. The enhanced mechanical performance with anti-aging properties indicated that the OH-Lig composite may replace PA12 to reduce the use of petrol-based materials. This work showed the method of purposefully designing and modifying lignin structures to fulfill the interaction requirements between lignin and polymer matrix. The as-prepared lignin could be used to prepare functional composites to achieve improved mechanical performance and targeted functions at the same time.

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

Manufacturing Induced Residual Stress Development and Deformation in PMDI-10 Structural Foam

In this report, we investigate how manufacturing conditions result in the warpage of moderate density PMDI polyurethane foam (12-50 lb/ft 3 ) when they are released from a mold. We have developed a multiphysics modeling framework to simulate the manufacturing process including resin injection, foaming and mold filling, gelation of the matrix, elevated cure, vitrification, cool down, and demolding. We have implemented this framework within the Sierra Mechanics Finite Element Code Suite. We couple Aria for flow, energy conservation, and foaming/curing kinetics with Adagio for the nonlinear viscoelastic solid response in a multi-staged simulation process flow. We calibrate a model for the PMDI-10S (10 lb/ft 3 free rise foam) through a suite of characterization data presented here to calibrate the solid cure behavior of the foam. The model is then used and compared to a benchmark experiment, the manufacturing and warpage over 1 year of a 10 cm by 10 cm by 2.5 cm foam "staple". This component features both slender and thick regions that warp considerably differently over time. Qualitative agreement between the model and the experiment is achieved but quantitative accuracy is not.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Interfaces in Oxides Formed on NiAlCr Doped with Y, Hf, Ti, and B

Abstract This study applies atom probe tomography (APT) to analyze the oxide scales formed on model NiAlCr alloys doped with Hf, Y, Ti, and B. Due to its ability to measure small amounts of alloying elements in the oxide matrix and its ability to quantify segregation, t he technique offers a possibility for detailed studies of the dopant’s fate during high-temperature oxidation. Three model NiAlCr alloys with different additions of Hf, Y, Ti, and B were prepared and oxidized in O 2 at 1,100°C for 100 h. All specimens showed an outer region consisting of different spinel oxides with relatively small grains and the protective Al 2 O 3 -oxide layer below. APT analyses focused mainly on this protective oxide layer. In all the investigated samples segregation of both Hf and Y to the oxide grain boundaries was observed and quantified. Neither B nor Ti were observed in the alumina grains or at the analyzed interfaces. The processes of formation of oxide scales and segregation of the alloying elements are discussed. The experimental challenges of the oxide analyses by APT are also addressed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Analyzer-free linear dichroic ptychography

Linear dichroism is an important tool to characterize the transmission matrix and determine the crystal or orbital orientation in a material. To achieve high-resolution mapping of transmission properties, we introduce the linear-dichroism scattering model in ptychographic imaging and develop an efficient two-stage reconstruction algorithm. Using the proposed algorithm on a uniaxial material, the dichroic transmission matrix can be recovered without an analyzer by using ptychography measurements with as few as three different polarization angles, with the help of an empty region to remove phase ambiguities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2011

Chemical order and local structure of the lead-free relaxor ferroelectric Na{sub 1/2}Bi{sub 1/2}TiO{sub 3}

The A-site mixed perovskite sodium bismuth titanate (Na{sub 1/2}Bi{sub 1/2})TiO{sub 3} (NBT) is investigated by means of first-principles calculations based on density functional theory. By studying different geometries with varying occupations of the A-site, the influence of chemical order on the thermodynamic stability and local structure is explored. We find that the hybridization of Bi 6sp with O 2p-states leads to stereochemically active Bi{sup 3+} lone pairs and increases the stability of structures with high Bi concentrations in {l_brace}001{r_brace}-planes. This goes along with displacive disorder on the oxygen sublattice, which up to now has been neglected in experimental studies. The calculated ordering energies are, however, small as compared to the thermal energy and therefore only short-range chemical order can be expected in experiments. Thus, it is conceivable that chemically ordered local areas can act as nucleation sites for polar nano-regions, which would explain the experimentally observed relaxor behavior of NBT. - Graphical abstract: First-principles calculations give relative stabilities of different chemically ordered structures. The results suggest a new model for the local structure of Na{sub 1/2}Bi{sub 1/2}TiO{sub 3} with 001-ordered nano-regions embedded in a chemically disordered matrix. Chemical order/disorder additionally induces displacive disorder within the oxygen sublattice. Highlights: > Lead-free relaxor ferroelectric Na{sub 1/2}Bi{sub 1/2}TiO{sub 3} is studied by ab initio-calculations. > Structural relaxations in the oxygen sublattice are decisive for relative stabilities. > Chemical environment of oxygen determines relaxation ability. > Bi 6s{sup 2} lone pair formation is the driving force for relaxation. > New structure model: Chemically 001-ordered nano-regions embedded in disordered matrix.

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

Global quantum phase diagram and non-Abelian chiral spin liquid in a spin- 3 2 square-lattice antiferromagnet

Since strong quantum fluctuations are essential for the emergence of quantum spin liquids, there have been extensive exploration and identification of spin liquid candidates in spin-$\frac{1}{2}$ systems, while such activities are rare in higher spin systems. Here we report an example of non-Abelian chiral spin liquid emerging in a spin-$\frac{3}{2}$ Heisenberg model on a square lattice. By tuning Heisenberg exchange interaction and scalar chirality interaction, we map out a quantum phase diagram enclosing three conventional magnetic orders and a chiral spin liquid based on density-matrix renormalization group studies. The nature of the spin liquid is identified as a long-sought bosonic version of the Read-Rezayi state that supports non-Abelian Fibonacci anyonic statistics, identified by the ground state entanglement spectrum. Significantly, we establish that the non-Abelian chiral spin liquid emerges through the enlarged local degrees of freedom and enhanced quantum fluctuations near the classical phase boundaries of competing magnetic orders. Finally, our numerical discovery of an exotic quantum spin liquid in a spin-$\frac{3}{2}$ system suggests a route for discovering fractionalized quantum phases in frustrated higher spin magnetic compounds.

2-dimensional systems↗
DOE OSTI2015

Underlying theory of a model for the Renner–Teller effect in tetra-atomic molecules: X{sup 2}Π{sub u} electronic state of C{sub 2}H{sub 2}{sup +}

In the present study, we prove the plausibility of a simple model for the Renner–Teller effect in tetra-atomic molecules with linear equilibrium geometry by ab initio calculations of the electronic energy surfaces and non-adiabatic matrix elements for the X{sup 2}Π{sub u} state of C{sub 2}H{sub 2}{sup +}. This phenomenon is considered as a combination of the usual Renner–Teller effect, appearing in triatomic species, and a kind of the Jahn–Teller effect, similar to the original one arising in highly symmetric molecules. Only four parameters (plus the spin–orbit constant, if the spin effects are taken into account), which can be extracted from ab initio calculations carried out at five appropriate (planar) molecular geometries, are sufficient for building up the Hamiltonian matrix whose diagonalization results in the complete low-energy (bending) vibronic spectrum. The main result of the present study is the proof that the diabatization scheme, hidden beneath the apparent simplicity of the model, can safely be carried out, at small-amplitude bending vibrations, without cumbersome computation of non-adiabatic matrix elements at large number of molecular geometries.

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

A dynamically bi-orthogonal method for time-dependent stochastic partial differential equations I: Derivation and algorithms

We propose a dynamically bi-orthogonal method (DyBO) to solve time dependent stochastic partial differential equations (SPDEs). The objective of our method is to exploit some intrinsic sparse structure in the stochastic solution by constructing the sparsest representation of the stochastic solution via a bi-orthogonal basis. It is well-known that the Karhunen–Loeve expansion (KLE) minimizes the total mean squared error and gives the sparsest representation of stochastic solutions. However, the computation of the KL expansion could be quite expensive since we need to form a covariance matrix and solve a large-scale eigenvalue problem. The main contribution of this paper is that we derive an equivalent system that governs the evolution of the spatial and stochastic basis in the KL expansion. Unlike other reduced model methods, our method constructs the reduced basis on-the-fly without the need to form the covariance matrix or to compute its eigendecomposition. In the first part of our paper, we introduce the derivation of the dynamically bi-orthogonal formulation for SPDEs, discuss several theoretical issues, such as the dynamic bi-orthogonality preservation and some preliminary error analysis of the DyBO method. We also give some numerical implementation details of the DyBO methods, including the representation of stochastic basis and techniques to deal with eigenvalue crossing. In the second part of our paper [11], we will present an adaptive strategy to dynamically remove or add modes, perform a detailed complexity analysis, and discuss various generalizations of this approach. An extensive range of numerical experiments will be provided in both parts to demonstrate the effectiveness of the DyBO method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Self-Consistent Two-Gap Description of MgB2 Superconductor

A self-consistent two-gap γ -model is used to quantitatively describe several thermodynamic properties of MgB 2 superconductor. The superconducting coupling matrix, ν i j , was obtained from the fitting of the superfluid density in the entire superconducting temperature range. Using this input, temperature-dependent superconducting gaps, specific heat, and upper critical fields were calculated with no adjustable parameters and compared with the experimental data as well as with the first-principles calculations. The observed agreement between fit and data shows that γ -model provides adequate quantitative description of the two-gap superconductivity in MgB 2 and may serve as a relatively simple and versatile self-consistent description of the thermodynamic quantities in multi-gap superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2012

Self-induced transparency in a dispersive medium

The evolution of the electromagnetic field in a two-level medium occurring in a matrix with finite response time has been studied. An integrable variant of the Maxwell-Bloch equations with allowance for nonlinear dispersion is derived and solved using the inverse scattering problem. It is shown that the nonlinear dispersion caused by the finite response time of the matrix yields a new possibility of controlling soliton parameters. A particular case of the constructed model can be used to describe field pulses in the parameter domain that occurs between the regions of applicability of the quasi-monochromatic approximation and the approximation of unidirectional propagation of pulses with durations on the order of the oscillation period.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

An inversion-relaxation approach for sampling stationary points of spin model Hamiltonians

Sampling the stationary points of a complicated potential energy landscape is a challenging problem. Here, we introduce a sampling method based on relaxation from stationary points of the highest index of the Hessian matrix. We illustrate how this approach can find all the stationary points for potentials or Hamiltonians bounded from above, which includes a large class of important spin models, and we show that it is far more efficient than previous methods. For potentials unbounded from above, the relaxation part of the method is still efficient in finding minima and transition states, which are usually the primary focus of attention for atomistic systems.

74 ATOMIC AND MOLECULAR PHYSICS↗