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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Theoretical Treatment of the Coriolis Effect Using Hyperspherical Coordinates, with Application to the Ro-Vibrational Spectrum of Ozone

Here, several alternative methods for the description of the interaction between rotation and vibration are compared and contrasted using hyperspherical coordinates for a triatomic molecule. These methods differ by the choice of the z -axis and by the assumption of a prolate or oblate rotor shape of the molecule. For each case, a block-structure of the rotational–vibrational Hamiltonian matrix is derived and analyzed, and the advantages and disadvantages of each method are made explicit. This theory is then employed to compute ro-vibrational spectra of singly substituted ozone; roughly, 600 vibrational states of 16 O 18 O 16 O and 16 O 16 O 18 O isomers combined, with rotational excitations up to J = 5 and both inversion parities (21600 coupled ro-vibrational states in total). Splittings between the states of different parities, so-called K -doublings, are calculated and analyzed. The roles of the asymmetric-top rotor term and the Coriolis coupling term are determined individually, and it is found that they both affect these splittings, but in the opposite directions. Thus, the two effects partially cancel out, and the residual splittings are relatively small. Energies of the ro-vibrational states reported in this work for 16 O 18 O 16 O and 16 O 16 O 18 O are in excellent agreement with literature (available for low-vibrational excitation). New data obtained here for the highly excited vibrational states enable the first systematic study of the Coriolis effect in symmetric and asymmetric isotopomers of ozone.

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

Second harmonic generation in topological insulators under quantizing magnetic fields

We theoretically investigate the second harmonic generation (SHG) of topological insulator surface states in a perpendicular magnetic field. Our theory is based on the microscopic expression of the second-order magneto-optical conductivity developed from the density matrix formalism, taking into account hexagonal warping effects on the surface states' band structure. Using numerically exact Landau-level energies and wave functions including hexagonal warping, we calculate the spectrum of SHG conductivities under normal incidence for different values of magnetic field and chemical potential. The imaginary parts of the SHG conductivities show prominent resonant peaks corresponding to one-photon and two-photon inter-Landau-level transitions. Treating the hexagonal warping term perturbatively, these transitions are clarified analytically within a perturbation theory from which approximate selection rules for the allowable optical transitions for SHG are determined. Furthermore, our results show extremely high SHG susceptibility that is easily tunable by magnetic field and doping level for topological surface states in the far-infrared regime, exceeding that of many conventional nonlinear materials. This work highlights the key role of hexagonal warping effects in generating second-order optical responses and provides insights into the nonlinear magneto-optical properties of the topological insulators.

Landau levels↗
DOE OSTI2024

Renormalization of beta decay at three loops and beyond

The anomalous dimension for heavy-heavy-light effective theory operators describing nuclear beta decay is computed through three-loop order in the static limit. The result at order Z 2 α 3 corrects a previous result in the literature. An all-orders symmetry is shown to relate the anomalous dimensions at leading and subleading powers of Z at a given order of α . The first unknown coefficient for the anomalous dimension now appears at O ( Z 2 α 4 ) . Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Accelerating the density-functional tight-binding method using graphical processing units

Acceleration of the density-functional tight-binding (DFTB) method on single and multiple graphical processing units (GPUs) was accomplished using the MAGMA linear algebra library. Herein two major computational bottlenecks of DFTB ground-state calculations were addressed in our implementation: the Hamiltonian matrix diagonalization and the density matrix construction. The code was implemented and benchmarked on two different computer systems: (1) the SUMMIT IBM Power9 supercomputer at the Oak Ridge National Laboratory Leadership Computing Facility with 1–6 NVIDIA Volta V100 GPUs per computer node and (2) an in-house Intel Xeon computer with 1–2 NVIDIA Tesla P100 GPUs. The performance and parallel scalability were measured for three molecular models of 1-, 2-, and 3-dimensional chemical systems, represented by carbon nanotubes, covalent organic frameworks, and water clusters.

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

Symmetry-resolved entanglement entropy, spectra & boundary conformal field theory

We perform a comprehensive analysis of the symmetry-resolved (SR) entanglement entropy (EE) for one single interval in the ground state of a 1 + 1D conformal field theory (CFT), that is invariant under an arbitrary finite or compact Lie group, G. We utilize the boundary CFT approach to study the total EE, which enables us to find the universal leading order behavior of the SREE and its first correction, which explicitly depends on the irreducible representation under consideration and breaks the equipartition of entanglement. We present two distinct schemes to carry out these computations. The first relies on the evaluation of the charged moments of the reduced density matrix. This involves studying the action of the defect-line, that generates the symmetry, on the boundary states of the theory. This perspective also paves the way for discussing the infeasibility of studying symmetry resolution when an anomalous symmetry is present. The second scheme draws a parallel between the SREE and the partition function of an orbifold CFT. This approach allows for the direct computation of the SREE without the need to use charged moments. From this standpoint, the infeasibility of defining the symmetry-resolved EE for an anomalous symmetry arises from the obstruction to gauging. Finally, we derive the symmetry-resolved entanglement spectra for a CFT invariant under a finite symmetry group. We revisit a similar problem for CFT with compact Lie group, explicitly deriving an improved formula for U(1) resolved entanglement spectra. Using the Tauberian formalism, we can estimate the aforementioned EE spectra rigorously by proving an optimal lower and upper bound on the same. In the abelian case, we perform numerical checks on the bound and find perfect agreement.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1986

Remote sensing of Earth Terrain

The Mueller matrix and polarization covariance matrix for polarimetric radar systems were studied. The clutter is modelled by a layer of random permittivity, described by a three dimensional correlation function, with variance, and horizontal and vertical correlation lengths. This model is applied, using the wave theory with Born approximations carried to the second order, to find the backscattering elements of the polarimetric matrices. It is found that 8 out of 16 elements of the Mueller matrix are identically zero, corresponding to a covariance matrix with four zero elements. Theoretical predictions are matched with experimental data for vegetation fields. A general mixing formula is derived for discrete scatterers immersed in a host medium. The inclusion particles are assumed to be ellipsoidal. The electric field inside the scatterers is determined by quasistatic analysis, assuming the diameter of the inclusion particles to be much smaller than the wavelength. The results are applicable to general multiphase mixtures.

Kong, J. A.↗
DOE OSTI2021

Designing high-performance superconductors with nanoparticle inclusions: Comparisons to strong pinning theory

One of the most promising routes for achieving high critical currents in superconductors is to incorporate dispersed, non-superconducting nanoparticles to control the dissipative motion of vortices. However, these inclusions reduce the overall superconducting volume and can strain the interlaying superconducting matrix, which can detrimentally reduce T c . Consequently, an optimal balance must be achieved between the nanoparticle density n p and size d. Determining this balance requires garnering a better understanding of vortex–nanoparticle interactions, described by strong pinning theory. Here, we map the dependence of the critical current on nanoparticle size and density in (Y 0.77 , Gd 0.23 )Ba 2 Cu 3 O 7-δ films in magnetic fields of up to 35 T and compare the trends to recent results from time-dependent Ginzburg–Landau simulations. We identify consistency between the field-dependent critical current J c (B) and expectations from strong pinning theory. Specifically, we find that J c ∝ B -α , where α decreases from 0.66 to 0.2 with increasing density of nanoparticles and increases roughly linearly with nanoparticle size d/ξ (normalized to the coherence length). At high fields, the critical current decays faster (~B -1 ), suggesting that each nanoparticle has captured a vortex. When nanoparticles capture more than one vortex, a small, high-field peak is expected in J c (B). Due to a spread in defect sizes, this novel peak effect remains unresolved here. Finally, we reveal that the dependence of the vortex creep rate S on nanoparticle size and density roughly mirrors that of α, and we compare our results to low-T nonlinearities in S(T) that are predicted by strong pinning theory.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Metal–Nitrogen–Carbon Cluster–Decorated Titanium Carbide is a Durable and Inexpensive Oxygen Reduction Reaction Electrocatalyst

Clusters of nitrogen- and carbon-coordinated transition metals dispersed in a carbon matrix (e. g., Fe-N-C) have emerged as an inexpensive class of electrocatalysts for the oxygen reduction reaction (ORR). In this work, it was shown that optimizing the interaction between the nitrogen-coordinated transition metal clusters embedded in a more stable and corrosion-resistant carbide matrix yielded an ORR electrocatalyst with enhanced activity and stability compared to Fe-N-C catalysts. Utilizing first-principles calculations, an electrostatics-based descriptor of catalytic activity was identified, and nitrogen-coordinated iron (FeN 4 ) clusters embedded in a TiC matrix were predicted to be an efficient platinum-group metal (PGM)-free ORR electrocatalyst. Guided by theory, selected catalyst formulations were synthesized, and it was demonstrated that the experimentally observed trends in activity fell exactly in line with the descriptor-derived theoretical predictions. The Fe-N-TiC catalyst exhibited enhanced activity (20 %) and durability (3.5-fold improvement) compared to a traditional Fe-N-C catalyst. It was posited that the electrostatics-based descriptor provides a powerful platform for the design of active and stable PGM-free electrocatalysts and heterogenous single-atom catalysts for other electrochemical reactions.

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

Notes on resonances and unitarity from celestial amplitudes

We study the celestial description of the O(N) sigma model in the large N limit as introduced by Coleman, Jackiw and Politzer. Focusing on three dimensions, we analyze the implications of a UV complete, all-loop order 4-point amplitude of pions in terms of correlation functions defined on the celestial circle. We find these retain many key features from the previously studied tree-level case, such as their relation to Generalized Free Field theories and crossing-symmetry, but also incorporate new properties such as IR/UV softness and S-matrix metastable states. In particular, to understand unitarity, we propose a form of the optical theorem that controls the imaginary part of the correlator based solely on the presence of these resonances. We also explicitly analyze the conformal block expansions and factorization of four-point functions into three-point functions. We find that summing over resonances is key for these factorization properties to hold. We end with some topics for future study.

1/N expansion↗
DOE OSTI2022

Accurate localization of Kosterlitz-Thouless-type quantum phase transitions for one-dimensional spinless fermions

We investigate the charge-density wave (CDW) transition for one-dimensional spinless fermions at half band filling with nearest-neighbor electron transfer amplitude t and interaction V. The model is equivalent to the anisotropic XXZ Heisenberg model for which the Bethe Ansatz provides an exact solution. For V>V c =2t, the CDW order parameter and the single-particle gap are finite but exponentially small, as is characteristic for a Kosterlitz-Thouless transition. It is notoriously difficult to locate such infinite-order phase transitions in the phase diagram using approximate analytical and numerical approaches. Second-order Hartree-Fock theory is qualitatively applicable for all interaction strengths, and predicts the CDW transition to occur at V$^{(2)}_{c,2}$≈1.5t. Second-order Hartree Fock theory is almost variational because the density of quasiparticle excitations is small. We apply the density-matrix renormalization group (DMRG) for periodic boundary conditions for system sizes up to 514 sites, which permits a reliable extrapolation of all physical quantities to the thermodynamic limit, apart from the critical region. We investigate the ground-state energy, the gap, the order parameter, the momentum distribution, the quasiparticle density, and the density-density correlation function to locate V c from the DMRG data. In conclusion, tracing the breakdown of the Luttinger liquid and the peak in the quasiparticle density at the band edge permits us to reproduce V c with an accuracy of one percent.

1-dimensional spin chains↗
DOE OSTI2021

Impact of three-body forces on elastic nucleon-nucleus scattering observables

In a previous series of papers we investigated the domain of applicability of chiral potentials to the construction of a microscopic optical potential (OP) for elastic nucleon-nucleus scattering. The OP was derived at the first order of the spectator expansion of the Watson multiple scattering theory and its final expression was a folding integral between the nucleon-nucleon ( N N ) t matrix and the nuclear density of the target. In the calculations N N and three-nucleon ( 3 N ) chiral interactions were used for the target density and only the N N interaction for the N N t matrix. The purpose of this study is to achieve another step towards the calculation of a more consistent OP introducing the 3 N force also in the dynamic part of the OP. The full treatment of the 3 N interaction is beyond our present capabilities. Thus, here it is approximated with a density dependent N N interaction obtained after the averaging over the Fermi sphere. In practice, in our model the 3 N force acts as a medium correction of the bare N N interaction used to calculate the t matrix. Even if the 3 N force is treated in an approximate way, this method naturally extends our previous model of the OP and allows a direct comparison of our present and previous results. We consider as case studies the elastic scattering of nucleons off C 12 and O 16 . We present results for the differential cross section and the spin observables for different values of the projectile energy. From the comparison with the experimental data and with the results of our previous model we assess the importance of the 3 N interaction in the dynamic part of the OP. Conclusions: Our analysis indicates that the contribution of the 3 N force in the t matrix is small for the differential cross section and it is sizable for the spin observables, in particular, for the analyzing power. We find that the two-pion exchange term is the major contributor to the 3 N force. A chiral expansion order-by-order analysis of the scattering observables confirms the convergence of our results at the next-to-next-to-next-to-leading-order, as already established in our previous work.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Quantum Melting of Generalized Wigner Crystals in Transition Metal Dichalcogenide Moiré Systems

The generalized Wigner crystal (GWC) is a novel quantum phase of matter driven by further-range interaction at fractional fillings of a lattice. The role of further-range interaction as the driver for the incompressible state is akin to the Wigner crystal. On the other hand, the significant role of commensurate filling is akin to the Mott insulator. Recent progress in simulator platforms presents unprecedented opportunities to investigate quantum melting in the strongly interacting regime through synergy between theory and experiments. However, the earlier theory literature presents diverging predictions. We study the quantum freezing of GWC through large-scale density matrix renormalization group simulations of a triangular lattice extended Hubbard model. We find a single first-order phase transition between the Fermi liquid and the √3×√3 GWC state. Furthermore, the GWC state shows long-range antiferromagnetic 120° Néel order. Our results present the simplest answers to the question of the quantum phase transition into the GWC phase and the properties of the GWC phase.

Density matrix renormalization group↗
DOE OSTI2021

A Modified Network Alteration Model to Predict Quasi-Static Behavior of the Cross-Linked Polymers During Hydrolytic Aging

Abstract The network alteration theory is utilized to capture the effects of hydrolytic aging on the mechanical integrity of the polymer matrix. Hydrolytic degradation generally considered as a result of two simultaneous chemical phenomena; (i) reduction of the cross-links, and (ii) active chain loss. Assuming the alteration of the chain density is identical to the change of the number average molecular weight along the aging trajectory, the strain energy of polymer matrix can be rewritten as a function of deformation, deformation history, storage time, and aging temperature. The proposed model includes only five physically inspired material parameters. Thus, while it is computationally efficient, it shows good agreement with own experimental data, which performed on various range of accelerated aging temperatures and times. With respect to its computational efficiency, simplicity, accuracy, and interpret-ability, the model is a good choice for advanced implementations in FE programs.

Bahrololoumi, Amir↗
NASA NTRS1973

Minimum impulse three-body trajectories.

A rapid and accurate method of calculating optimal impulsive transfers in the restricted problem of three bodies has been developed. The technique combines a multi-conic method of trajectory integration with primer vector theory and an accelerated gradient method of trajectory optimization. A unique feature is that the state transition matrix and the primer vector are found analytical without additional integrations or differentiations. The method has been applied to the determination of optimal two and three impulse transfers between the L2 libration point and circular orbits about both the earth and the moon.

D'Amario, L.↗
NASA NTRS1980

Predicting the time-temperature dependent axial failure of B/A1 composites

Experimental and theoretical studies were conducted in order to understand and predict the effects of time, temperature, and stress on the axial failure modes of boron fibers and B/A1 composites. Due to the anelastic nature of boron fiber deformation, it was possible to determine simple creep functions which can be employed to accurately describe creep and fracture stress of as-produced fibers. Analysis of damping and strength data for B/6061 A1 composites indicates that fiber creep effects of creep on fiber fracture are measurably reduced by the composite fabrication process. The creep function appropriate for fibers with B/Al composites was also determined. A fracture theory is presented for predicting the time-temperature dependence of the axial tensile strength for metal matrix composites in general and B/A1 composites in particular.

Dicarlo, J. A.↗
NASA NTRS1987

Recursive Robot-Arm Dynamics via Filtering and Smoothing

Forward and inverse dynamics solved using Kalman filtering and Bryson-Frazier smoothing. Dynamics of serial-link robot arm solved by using recursive techniques from linear filtering and smoothing theory. Solutions of dynamical equations give forces, moments, and accelerations at joints between links, and multilink inertia matrix and its inverse. Theoretical developments lay foundation for use of filtering and smoothing techniques in design of robot controls.

Rodriguez, Guillermo↗
NASA NTRS1990

Stabilization of Taylor-Couette flow due to time-periodic outer cylinder oscillation

The linear stability of circular Couette flow between concentric infinite cylinders is considered for the case when the inner cylinder is rotated at a constant angular velocity and the outer cylinder is driven sinusoidally in time with zero mean rotation. This configuration was studied experimentally by Walsh and Donnelly. The critical Reynolds numbers calculated from linear stability theory agree well with the experimental values, except at large modulation amplitudes and small frequencies. The theoretical values are obtained using Floquet theory implemented in two distinct approaches: a truncated Fourier series representation in time, and a fundamental solution matrix based on a Chebyshev pseudospectral representation in space. For large amplitude, low frequency modulation, the linear eigenfunctions are temporally complex, consisting of a quiescent region followed by rapid change in the perturbed flow velocities.

Murray, B. T.↗