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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Dynamical origins of weakly coupled relaxor behavior in Sn-doped (Ba,Ca)TiO 3 -BiScO 3

The peculiar characteristics of relaxors, viz., a frequency-dependent dielectric permittivity peak and good functional properties (dielectric, electromechanical, electrocaloric, etc.), are attributed to nanoscale regions with correlated dipoles, or polar nanoregions (PNRs). However, the exact nature of PNRs and their contribution to relaxor behavior remains debatable. In recent years, solid solutions of BaTiO3-BiMeO3 (where Me is a metal), have emerged as an interesting system with characteristics in between that of relaxors and dipole glasses. Here, we have examined the atomistic origins of weakly coupled relaxor behavior, specifically with regard to formation of PNRs, in Sn-doped (1-x)(Ba,Ca)TiO 3-x BiScO 3 using macroscopic polarization and neutron dynamic pair distribution function measurements. We show that the short-range atomic correlations observed within the PNRs dynamically fluctuate with frequencies of the order of THz. Furthermore the composition-dependent dielectric and polarization behaviors are critically influenced by the relative stability of the atomic correlations near ~1 THz, while the instantaneous atomic correlations are largely independent of x. The current results are discussed based on a model of intrinsic local modes distributed in a dielectrically soft matrix.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Ab initio calculation of muon capture on 24 Mg

In this work we study ordinary muon capture (OMC) on 24 Mg from a first-principles perspective. Starting from a particular two- and three-nucleon interaction derived from chiral effective field theory, we use the valence-space in-medium similarity renormalization group (VS-IMSRG) framework to construct effective Hamiltonians and muon-capture operators, which nonperturbatively account for many-body physics outside the valence space. Here, the obtained nuclear matrix elements are compared against those from the phenomenological shell model. The impact of including the correlations from the nuclear shell model (NSM) as well as including the induced twobody part is studied in detail. Furthermore, the effects of realistic bound-muon wave function on the operators is studied. Finally, predictions for capture rates to the lowest excited states in 24 Na are given and compared with available data. It is found that the spectroscopic properties of 24 Mg and its OMC daughter 24 Na are fairly well described by both the NSM and VS-IMSRG, and that the effect of the hadronic two-body currents significantly reduces the OMC rates. Both models have some difficulties in matching the measured OMC rates, especially for the 2 + final states. This calls for further studies in other light nuclei with available OMC data.

20 ≤ A ≤ 38↗
DOE OSTI2024

Correlations of conserved quantities at finite baryon density

Correlations involving the seven conserved quantities, namely energy, baryon number, electric charge, strangeness, and the three components of momentum, give rise to correlations in heavy-ion collisions. Through the utilization of a simple one-dimensional hydrodynamic model, we calculate the evolution of the entire 7×7 matrix of correlations as a function of relative spatial rapidity. This comprehensive analysis accounts for finite baryon density, which results in off-diagonal correlations between the charge-related quantities and the energy-momentum quantities. These correlations in coordinate space are subsequently transformed into correlations in momentum space using statistical weighting. As a result, the entire matrix of correlations is revealed to be highly sensitive to the equation of state, viscosity, and diffusivity.

Equations of state of nuclear matter↗
DOE OSTI2025

Bridging hadronic and vacuum structure by heavy quarkonia

We discuss the central and, mostly, spin-dependent potentials in heavy quarkonia $\overline{b}$$b$,$\overline{c}$$c$, with two goals in mind. The first is phenomenological: using the splitting between the 1S and 2S pairs, as well as the 1P and 2P quartet masses, we obtain pertinent matrix elements of all these potentials. The second is theoretical: using standard wave functions, we compute these matrix elements, quantifying the perturbative, nonperturbative and “string” contributions. The model for nonperturbative effects is a “dense instanton liquid model,” in which the QCD vacuum is made mostly of (strongly color correlated) instanton-anti-instanton pairs or “molecules.” Their contributions are calculated via standard Wilson lines, with or without extra powers of gauge fields. We find a reasonable description of all central, spin-spin, and spin-orbit forces at distances 𝑟 = 0–0.7 fm, relevant for $\overline{b}$$b$ and $\overline{c}$$c$ quarkonia.

Bound states↗
DOE OSTI2021

Origin of Insulating Ferromagnetism in Iron Oxychalcogenide Ce 2 O 2 FeSe 2

An insulating ferromagnetic (FM) phase exists in the quasi-one-dimensional iron oxychalcogenide Ce 2 O 2 FeSe 2 , but its origin is unknown. To understand the FM mechanism, here a systematic investigation of this material is provided, analyzing the competition between ferromagnetic and antiferromagnetic tendencies and the interplay of hoppings, Coulomb interactions, Hund’s coupling, and crystal-field splittings. Overall, our intuitive analysis based on second-order perturbation theory shows that large entanglements between doubly occupied and half filled orbitals play a key role in stabilizing the FM order in Ce 2 O 2 FeSe 2 . In addition, via many-body computational techniques applied to a multiorbital Hubbard model, the phase diagram confirms the proposed FM mechanism.

1-dimensional spin chains↗
DOE OSTI2024

Evidence of Dirac Quantum Spin Liquid in YbZn 2⁢ GaO 5

The emergence of a quantum spin liquid (QSL), a state of matter that can result when electron spins are highly correlated but do not become ordered, has been the subject of a considerable body of research in condensed matter physics [1,2]. Spin liquid states have been proposed as hosts for high-temperature superconductivity [3] and can host topological properties with potential applications in quantum information science [4]. The excitations of most quantum spin liquids are not conventional spin waves but rather quasiparticles known as spinons, whose existence is well established experimentally only in one-dimensional systems; the unambiguous experimental realization of QSL behavior in higher dimensions remains challenging. Here, in this work, we investigate the novel compound $YbZn_{2}GaO_{5}$, which hosts an ideal triangular lattice of effective spin-1/2 moments with no detectable inherent chemical disorder. Thermodynamic and inelastic neutron scattering measurements performed on high-quality single crystal samples of $YbZn_{2}GaO_{5}$ exclude the possibility of long-range magnetic ordering down to 0.06 K, demonstrate a quadratic power law for the specific heat and reveal a continuum of magnetic excitations in parts of the Brillouin zone. Both low-temperature thermodynamics and inelastic neutron scattering spectra suggest that $YbZn_{2}GaO_{5}$ is a U(1) Dirac QSL with spinon excitations concentrated at certain points in the Brillouin zone. We advanced these results by performing additional specific heat measurements under finite fields, further confirming the theoretical expectations for a Dirac QSL on the triangular lattice of $YbZn_{2}GaO_{5}$.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Controlled expansion for pairing in a polarized band with strong repulsion

Can strong repulsive interactions be shown to give rise to pairing in a controlled way? We find that for a single flavor polarized band, there is a small expansion parameter in the low density limit, once the Bloch wavefunction form factor is taken into account. A perturbative expansion is possible, even if the interaction is much stronger than the Fermi energy ε F . As a matter of principle, our work shows analytically how strong pairing can emerge from strong repulsion. Here, we illustrate our method with two examples: a 2D Dirac model and a 1D tight binding model with two orbitals. In the latter case, using density matrix renormalization group, we show that the analytical theory indeed guided us to discover the parameter regime where p-wave pairing with order-1 strength is dominant.

Ferromagnetic superconductors↗
DOE OSTI2021

Analyzing the Effects of Cyberattacks on Distribution System State Estimation

Key components of power systems—such as energy management systems, automatic generation control, and state estimation—are under serious vulnerability from cyberattacks. Cyber threats in electric grids have increased significantly because of the increased interconnectivity of supervisory control and data acquisition systems and public network infrastructure. As the penetration level of distributed energy resources increases, it is imperative to employ system-monitoring techniques such as state estimation for the reliable operation of distribution systems. Recently, multiple methods have been developed that exploit the low rank property of distribution system state matrix and are robust to bad data, such as matrix completion. This paper analyzes the impact of various realistic cyberattack scenarios on matrix completion. Realistic cyberattack scenarios are converted into data corruption models that are used in an extensive simulation of a custom IEEE 123-bus system.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2023

Bayesian High-Rank Hankel Matrix Completion for Nonlinear Synchrophasor Data Recovery

Phasor measurement units (PMUs) provide high temporal-resolution synchrophasor measurements for power system monitoring and control. The frequent data quality issues, such as missing and bad data, prevent the incorporation of synchrophasor data in real-time operations. Most existing data-driven data recovery methods assume the power system dynamics can be approximated by a linear dynamical system, and the recovery performance degrades significantly when the power system is experiencing nonlinear dynamics during significant events. Here, this paper proposes a data-driven Bayesian nonlinear synchrophasor data recovery method (Ba-NSDR) that can recover a consecutive time period of simultaneous data losses or errors across all channels, even when the underlying system is highly nonlinear. The idea is to lift the Hankel matrix of the spatial-temporal synchrophasor data to a higher dimension such that the lifted Hankel matrix is low-rank in that space and can be processed with the kernel trick. Our proposed Bayesian method then infers the probabilistic distributions of synchrophasor from the partial observations. Some distinctive features of Ba-NSDR include an uncertainty index to measure the accuracy of the recovery result and the robustness to parameter selections. Our method is verified on both synthetic and recorded event datasets.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Forced Oscillation Source Location of Bulk Power Systems Using Synchrosqueezing Wavelet Transform

Forced oscillation source location (FOSL) plays a significant role in mitigating forced oscillations (FOs), which threaten the power system stability. Here, this paper proposes a data-driven approach for FOSL in power systems using synchrosqueezing wavelet transform (SWT). The proposed approach conducts SWT on the measured system responses to obtain the SWT matrix. Then, the SWT-based dissipating energy flow (DEF) model in time-frequency domain and dissipating energy spectrum (DES) model in frequency domain are derived from the traditional DEF model. Further, the characteristics of SWT-based DEF and DES are revealed by referring to the traditional DEF, and the FOSL criteria of the SWT-based DEF and DES can be hereby obtained. Using the obtained FOSL criteria, the FO source can be located from the measured responses. The performance of the proposed FOSL method is evaluated by simulation data of the WECC 179-bus test system and field-measurement PMU data of the ISO New England. The results confirm the accuracy and efficiency of the proposed method in the FOSL.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Porewater Geochemical Assessment of Seismic Indications for Gas Hydrate Presence and Absence: Mahia Slope, East of New Zealand’s North Island

We compare sediment vertical methane flux off the Mahia Peninsula, on the Hikurangi Margin, east of New Zealand’s North Island, with a combination of geochemical, multichannel seismic and sub-bottom profiler data. Stable carbon isotope data provided an overview of methane contributions to shallow sediment carbon pools. Methane varied considerably in concentration and vertical flux across stations in close proximities. At two Mahia transects, methane profiles correlated well with integrated seismic and TOPAS data for predicting vertical methane migration rates from deep to shallow sediment. However, at our “control site”, where no seismic blanking or indications of vertical gas migration were observed, geochemical data were similar to the two Mahia transect lines. This apparent mismatch between seismic and geochemistry data suggests a potential to underestimate gas hydrate volumes based on standard seismic data interpretations. To accurately assess global gas hydrate deposits, multiple approaches for initial assessment, e.g., seismic data interpretation, heatflow profiling and controlled-source electromagnetics, should be compared to geochemical sediment and porewater profiles. A more thorough data matrix will provide better accuracy in gas hydrate volume for modeling climate change and potential available energy content.

03 NATURAL GAS↗
DOE OSTI2021

A Joint Measurement of $\nu_{\mu}$-Disappearance and $\nu_{e}$e-Appearance in the NuMI beam using the NOvA Experiment

The discovery of neutrino oscillations provides the first indication of a lepton flavor violating (LFV) process, one that isn't predicted by the Standard Model. As such, NOvA is part of a rich experimental program to constrain unknown parameters in the neutrino oscillation model, described for three neutrino flavors using the PMNS unitary matrix. It is a long baseline experiment utilizing two detectors, a Near Detector at Fermilab and a Far Detector in Ash River, Minnesota for a total baseline of 810 km. It receives a predominantly ${\nu_\mu}/\overline{\bar\nu_\mu}$ beam peaking at 1:8 GeV from the NuMI beam facility at Fermilab. There are four oscillation channels used in the analysis, ${\nu_\mu} \to {\nu_\mu}$, ${\nu_\mu} \to {\nu_e}$, ${\overline{\nu}_\mu} \to {\overline{\nu}_\mu}$ and ${\overline{\nu}_\mu} \to {\overline{\nu}_e }$. With a total exposure of 13:6_10$^{20}$ and 12:5_10$^{20}$ protons on target for the neutrino and anti-neutrino beam modes respectively, 82 candidates are seen in the ${\nu_\mu} \to {\nu_e}$ channel for a total predicted backgroundof 26:8 events. Similarly, 33 candidates are seen in the corresponding anti-neutrino channel for a total predicted background of 14:0 events. In the ${\nu_\mu} \to {\nu_\mu}$ $(\overline{\nu}_\mu \to \overline{\nu}_\mu}$ channel, 211 (105) candidates are seen with an expectation of 1156.1 (488.1) events at no oscillations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

A model for discrete fracture-clay rock interaction incorporating electrostatic effects on transport

Abstract A model based on the code CrunchClay is presented for a fracture-clay matrix system that takes electrostatic effects on transport into account. The electrostatic effects on transport include those associated with the development of a diffusion potential as captured by the Nernst-Planck equation, and the formation of a diffuse layer bordering negatively charged clay particles within which partial anion exclusion occurs. The model is based on a dual continuum formulation that accounts for diffuse layer and bulk water pore space, providing a more flexible framework than is found in the classical mean electrostatic potential models. The diffuse layer model is obtained by volume averaging ion concentrations in the Poisson-Boltzmann equation, but also includes the treatment of longitudinal transport within this continuum. The calculation of transport within the bulk and diffuse layer porosity is based on a new formulation for the Nernst-Planck equation that considers averaging of diffusion coefficients and accumulation factors at grid cell interfaces. Equations for function residuals and the associated Jacobian matrix are presented such that the system of nonlinear differential-algebraic equations can be solved with Newton’s method. As an example, we consider a 2D system with a single discrete fracture within which flow and advective transport occurs that is coupled to diffusion in the clay-rich matrix. The simulation results demonstrate the lack of retardation for anions (e.g., 36 Cl − ) of the contaminant plume within the fracture flow system because they are largely excluded from the charged clay rock, while the migration of cations (e.g., 90 Sr ++ ) is more strongly attenuated. The diffusive loss of divalent cations in particular from the fracture is accentuated by their accumulation in the diffuse layer within the clay-rich matrix.

58 GEOSCIENCES↗
DOE OSTI2021

Mitigating baryonic effects with a theoretical error covariance

ABSTRACT One of the primary sources of uncertainties in modelling the cosmic-shear power spectrum on small scales is the effect of baryonic physics. Accurate cosmology for stage-IV surveys requires knowledge of the matter power spectrum deep in the non-linear regime at the per cent level. Therefore, it is important to develop reliable mitigation techniques to take into account baryonic uncertainties if information from small scales is to be considered in the cosmological analysis. In this work, we develop a new mitigation method for dealing with baryonic physics for the case of the shear angular power spectrum. The method is based on an augmented covariance matrix that incorporates baryonic uncertainties informed by hydrodynamical simulations. We use the results from 13 hydrodynamical simulations and the residual errors arising from a fit to a ΛCDM model using the extended halo model code HMCode to account for baryonic physics. These residual errors are used to model a so-called theoretical error covariance matrix that is added to the original covariance matrix. In order to assess the performance of the method, we use the 2D tomographic shear from four hydrodynamical simulations that have different extremes of baryonic parameters as mock data and run a likelihood analysis comparing the residual bias on Ωm and σ8 of our method and the HMCode for an LSST-like survey. We use different modelling of the theoretical error covariance matrix to test the robustness of the method. We show that it is possible to reduce the bias in the determination of the tested cosmological parameters at the price of a modest decrease in the precision.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Matrix Diffusion Controls Mountain Hillslope Groundwater Ages and Inferred Storage Dynamics

Groundwater age distributions provide fundamental insights on coupled water and biogeochemical processes in mountain watersheds. Field-based studies have found mixtures of young and old-aged groundwater in mountain catchments underlain by bedrock; yet, the processes that dictate these groundwater age distributions are poorly understood. In this work, we use the coupled ParFlow-CLM integrated hydrologic and EcoSLIM particle tracking models to simulate groundwater age distributions on a lower montane hillslope in the East River Watershed, Colorado (USA). We develop a convolution-based approach to propagate fracture-matrix diffusion processes to the EcoSLIM advection-dominated age distributions. We compare observed 3 H and 4 He concentrations from two groundwater wells against model predictions that have varying advective transport times and matrix diffusion magnitudes. Based on a Monte Carlo analysis that considers uncertain matrix and fracture parameters, we find that matrix diffusion is needed to jointly predict 3 H and 4 He observations at both wells. The advection-dominated age distributions lack adequate mixing of young and old-aged water to capture the observed co-occurrence of 3 H and 4 He. The model scenario that best matches the 3 H, 4 He, and water level observations when considering both advective flowpath and matrix diffusion mixing processes has a dynamic bedrock groundwater reservoir that is susceptible to considerable storage losses during low-snow periods. This dynamic groundwater system amplifies the need to assimilate deeper bedrock groundwater into watershed hydro-biogeochemical predictions. This work further highlights the importance of considering matrix diffusion when interpreting environmental tracers in bedrock groundwater systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Quantum Monte Carlo Approaches to Na Intercalation on Bilayer Graphene

We have performed Quantum Monte Carlo (QMC) simulations on Na-intercalated bilayer graphene to study the evolution of electronic and optical properties upon Na intercalation into hard carbon layers. The objective was to model the optimal configuration of Na intercalation into a hard carbon matrix containing graphene regions. Our study showed that Na intercalation can be energetically stabilized at large interlayer distances (over 6 Å) in both AA- and AB-stacked bilayer graphene. In the QMC results, we found a significant band gap opening at the equilibrium interlayer distance of Na-intercalated bilayer graphene, while corresponding density functional theory (DFT) results showed no gap. This difference between DFT and QMC results indicates that the gap opening induced by Na intercalation into a hard carbon is underestimated within the DFT framework. In addition, a zigzag configuration of Na atoms was found to be energetically stable at interlayer distances up to 10 Å, leading us to predict the existence of a local minimum of Na intercalation at large interlayer distance. These computation and modeling results can provide guidance on how to synthesize and optimize hard carbon with bilayer graphene regions that permit a zigzag intercalation configuration that will maximize and stabilize sodium hosting.

Binding energy↗
DOE OSTI2020

Current and future neutrino oscillation constraints on leptonic unitarity

The unitarity of the lepton mixing matrix is a critical assumption underlying the standard neutrino-mixing paradigm. However, many models seeking to explain the as-yet-unknown origin of neutrino masses predict deviations from unitarity in the mixing of the active neutrino states. Motivated by the prospect that future experiments may provide a precise measurement of the lepton mixing matrix, we revisit current constraints on unitarity violation from oscillation measurements and project how next-generation experiments will improve our current knowledge. With the next-generation data, the normalizations of all rows and columns of the lepton mixing matrix will be constrained to ≲10% precision, with the e-row best measured at ≲1% and the τ-row worst measured at ~10% precision. The measurements of the mixing matrix elements themselves will be improved on average by a factor of 3. We highlight the complementarity of DUNE, T2HK, JUNO, and IceCube Upgrade for these improvements, as well as the importance of ντ appearance measurements and sterile neutrino searches for tests of leptonic unitarity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Plaquette valence bond state in the spin-$\frac{1}{2}$ J 1 –J 2 XY model on the square lattice

Here, we studied the ground-state phase diagram of spin-1/2 J 1 –J 2 XY model on the square lattice with first-neighbor J 1 and second-neighbor J 2 antiferromagnetic interactions using both infinite density matrix renormalization group (iDMRG) and DMRG approaches. We show that a plaquette valence bond phase is realized in an intermediate region 0.50 ≤ J 2 /J 1 < 0.54 between a Néel magnetic ordered phase at J 2 /J 1 < 0.50 and a stripy magnetic ordered phase at J 2 /J 1 ≥ 0.54 . The plaquette valence bond phase is characterized by finite dimer orders in both the horizontal and vertical directions. Contrary to the spin-1/2 J 1 –J 2 Heisenberg model, we do not find numerical evidence for a quantum spin liquid phase in the J 1 –J 2 XY model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗