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

An anisotropic viscoplasticity model for shale based on layered microstructure homogenization

Viscoplastic deformation of shale is frequently observed in many subsurface applications. Many studies have suggested that this viscoplastic behavior is anisotropic—specifically, transversely isotropic—and closely linked to the layered composite structure at the microscale. In this work, we develop a two-scale constitutive model for shale in which anisotropic viscoplastic behavior naturally emerges from semianalytical homogenization of a bilayer microstructure. The microstructure is modeled as a composite of soft layers, representing a ductile matrix formed by clay and organics, and hard layers, corresponding to a brittle matrix composed of stiff minerals. This layered microstructure renders the macroscopic behavior anisotropic, even when the individual layers are modeled with isotropic constitutive laws. Using a common correlation between clay and organic content and magnitude of creep, we apply a viscoplastic modified Cam-Clay plasticity model to the soft layers, while treating the hard layers as a linear elastic material to minimize the number of calibration parameters. We then describe the implementation of the proposed model in a standard material update subroutine. The model is validated with laboratory creep data on samples from three gas shale formations. We also demonstrate the computational behavior of the proposed model through simulation of time-dependent borehole closure in a shale formation with different bedding plane directions.

58 GEOSCIENCES↗
DOE OSTI2025

Classical Benchmarks for Variational Quantum Eigensolver Simulations of the Hubbard Model

Simulating the Hubbard model is of great interest to a wide range of applications within condensed matter physics, however its solution on classical computers remains challenging in dimensions larger than one. The relative simplicity of this model, embodied by the sparseness of the Hamiltonian matrix, allows for its efficient implementation on quantum computers, and for its approximate solution using variational algorithms such as the variational quantum eigensolver. While these algorithms have been shown to reproduce the qualitative features of the Hubbard model, their quantitative accuracy in terms of producing true ground state energies and other properties, and the dependence of this accuracy on the system size and interaction strength, the choice of variational ansatz, and the degree of spatial inhomogeneity in the model, remains unknown. Here we present a rigorous classical benchmarking study, demonstrating the potential impact of these factors on the accuracy of the variational solution of the Hubbard model on quantum hardware, for systems with up to 32 qubits. We find that even when using the most accurate wavefunction ansätze for the Hubbard model, the error in its ground state energy and wavefunction plateaus for larger lattices, while stronger electronic correlations magnify this issue. Concurrently, spatially inhomogeneous parameters and the presence of off-site Coulomb interactions only have a small effect on the accuracy of the computed ground state energies. Our study highlights the capabilities and limitations of current approaches for solving the Hubbard model on quantum hardware, and we discuss potential future avenues of research.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Block Mott insulating state induced by next-nearest-neighbor hopping in the S = 3 2 zigzag chain BaCoTe 2 O 7

Quasi-one-dimensional correlated electronic multiorbital systems with either ladder or chain geometries continue attracting considerable interest due to their complex electronic phases arising from the interplay of the hopping matrix, the crystal-field splitting, the electronic correlations (Hubbard repulsion U and Hund coupling J H ), and strong quantum fluctuations. Recently, the intriguing cobalt zigzag chain system BaCoTe 2 ⁢ O 7 , with electronic density n=7, was prepared experimentally. Here, we systematically study the electronic and magnetic properties of this quasi-one-dimensional compound from the theoretical perspective. Based on first-principles density functional theory calculations, strongly anisotropic one-dimensional electronic Co 3⁢d bands were found near the Fermi level. By evaluating the relevant hopping amplitudes, we provide the magnitude and origin of the nearest-neighbor (NN) and next-nearest-neighbor (NNN) hopping matrices in BaCoTe 2 ⁢ O 7 . With this information, we constructed a three-orbital electronic Hubbard model for this zigzag chain system, and studied two cases: with only a NN hopping matrix, and with NN plus NNN hopping matrices. Introducing the Hubbard and Hund couplings and studying the model via the density matrix renormalization group method, we constructed the ground-state phase diagram. A robust staggered ↑ - ↓ - ↑ - ↓ antiferromagnetic (AFM) region was found when only the NN hopping matrix in the chain direction was employed. However, for the realistic case where the NNN hopping matrix is also included, the dominant state becomes instead a block AFM ↑ - ↑ - ↓ - ↓ order, in agreement with experiments. The system displays Mott insulator characteristics with three half-filled orbitals, when the block AFM order is stable. In conclusion, our results for BaCoTe 2 ⁢O 7 provide guidance to experimentalists and theorists working on this zigzag one-dimensional chain and related materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

A nonlinear Bloch model for Coulomb interaction in quantum dots

In this paper, we first derive a Coulomb Hamiltonian for electron–electron interaction in quantum dots in the Heisenberg picture. Then we use this Hamiltonian to enhance a Bloch model, which happens to be nonlinear in the density matrix. The coupling with Maxwell equations in case of interaction with an electromagnetic field is also considered from the Cauchy problem point of view. The study is completed by numerical results and a discussion about the advisability of neglecting intra-band coherences, as is done in part of the literature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Search for New Physics via EFT in collisions containing a top quark and a boosted Z or Higgs boson

We present a generator-level sensitivity study of physics beyond the Standard Model within the framework of Effective Field Theory (EFT). The targeted signal processes are single top quark production in association with a Lorentz-boosted Higgs or Z boson (tHq and tZq) in proton-proton collisions at a center-of-mass energy of 13.6 TeV. Event selection requires one electron or muon, a high-transverse-momentum, large-radius (AK8) jet reconstructing hadronic H → b̄b or Z → b̄b decays, at least two distinct, non-overlapping small-radius (AK4) jets with one originating from a b quark, and a missing transverse energy greater than 30 GeV. Signal and background samples are generated via Monte Carlo simulation using MadGraph with MadJax to perform EFT matrix-element calculations. Deviations from the Standard Model are parametrized through dimension-six EFT operators that probe modified top quark interactions. This generator-level study identifies kinematic observables with enhanced sensitivity to EFT operators and motivates future reconstruction-level analyses of tHq and tZq channels.

Meshramkar, Daniel [Baylor U., Waco]↗
DOE OSTI2026

Search for New Physics via EFT in collisions containing a top quark and a boosted Z or Higgs boson

We present a generator-level sensitivity study of physics beyond the Standard Model within the framework of Effective Field Theory (EFT). The targeted signal processes are single top quark production in association with a Lorentz-boosted Higgs or Z boson (tHq and tZq) in proton-proton collisions at a center-of-mass energy of 13.6 TeV. Event selection requires one electron or muon, a high-transverse-momentum, large-radius (AK8) jet reconstructing hadronic H → b̄b or Z → b̄b decays, at least two distinct, non-overlapping small-radius (AK4) jets with one originating from a b quark, and a missing transverse energy greater than 30 GeV. Signal and background samples are generated via Monte Carlo simulation using MadGraph with MadJax to perform EFT matrix-element calculations. Deviations from the Standard Model are parametrized through dimension-six EFT operators that probe modified top quark interactions. This generator-level study identifies kinematic observables with enhanced sensitivity to EFT operators and motivates future reconstruction-level analyses of tHq and tZq channels.

Meshramkar, Daniel [Baylor U., Waco]↗
DOE OSTI2021

Seasonal Disorder in Urban Traffic Patterns: A Low Rank Analysis

This article proposes several advances to sparse nonnegative matrix factorization (SNMF) as a way to identify large-scale patterns in urban traffic data. The input to our model is traffic counts organized by time and location. Nonnegative matrix factorization additively decomposes this information, organized as a matrix, into a linear sum of temporal signatures. Penalty terms encourage this factorization to concentrate on only a few temporal signatures, with weights which are not too large. Our interest here is to quantify and compare the regularity of traffic behavior, particularly across different broad temporal windows. In addition to the rank and error, we adapt a measure introduced by Hoyer to quantify sparsity in the representation. Combining these, we construct several curves which quantify error as a function of rank (the number of possible signatures) and sparsity; as rank goes up and sparsity goes down, the approximation can be better and the error should decreases. Plots of several such curves corresponding to different time windows leads to a way to compare disorder/order at different time scalewindows. In this paper, we apply our algorithms and procedures to study a taxi traffic dataset from New York City. In this dataset, we find weekly periodicity in the signatures, which allows us an extra framework for identifying outliers as significant deviations from weekly medians. We then apply our seasonal disorder analysis to the New York City traffic data and seasonal (spring, summer, winter, fall) time windows. We do find seasonal differences in traffic order.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Markov state models from short non-equilibrium simulations—Analysis and correction of estimation bias

To begin, many state-of-the-art methods for the thermodynamic and kinetic characterization of large and complex biomolecular systems by simulation rely on ensemble approaches, where data from large numbers of relatively short trajectories are integrated. In this context, Markov state models (MSMs) are extremely popular because they can be used to compute stationary quantities and long-time kinetics from ensembles of short simulations, provided that these short simulations are in “local equilibrium” within the MSM states. However, over the last 15 years since the inception of MSMs, it has been controversially discussed and not yet been answered how deviations from local equilibrium can be detected, whether these deviations induce a practical bias in MSM estimation, and how to correct for them. In this paper, we address these issues: We systematically analyze the estimation of MSMs from short non-equilibrium simulations, and we provide an expression for the error between unbiased transition probabilities and the expected estimate from many short simulations. We show that the unbiased MSM estimate can be obtained even from relatively short non-equilibrium simulations in the limit of long lag times and good discretization. Further, we exploit observable operator model (OOM) theory to derive an unbiased estimator for the MSM transition matrix that corrects for the effect of starting out of equilibrium, even when short lag times are used. Finally, we show how the OOM framework can be used to estimate the exact eigenvalues or relaxation time scales of the system without estimating an MSM transition matrix, which allows us to practically assess the discretization quality of the MSM. Applications to model systems and molecular dynamics simulation data of alanine dipeptide are included for illustration. The improved MSM estimator is implemented in PyEMMA of version 2.3.

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

Optimizing the Performance of Sparse-Matrix Vector Products on Next-Generation Processors

Matrix-vector products are ubiquitous in high-performance scientific applications and have a growing set of occurrences in advanced data analysis activities. Achieving high performance for these kernels is therefore paramount, in part, because these operations can consume vast amounts of application execution time. In this report we document the development of several sparse-matrix vector product kernel implementations using a variety of programming models and approaches. Each kernel is run on a broad set of matrices selected to demonstrate the wide variety of matrix structure and sparsity that is possible with a single, generic kernel. For benchmarking and performance analysis, we utilize leading computing architectures for the NNSA/ASC program including Intel's Knights Landing processor and IBM's POWER8.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Preliminary Interpretations from Fall 2018 Field Test at Ne’ot Hovav

Presentation conclusions: Conservative tracer transport was reasonably consistent with earlier tests between RH11c and RH11a; No attempt to model bentonite colloids because of minimal elevation of Al concentrations over background; Modeling confirms what we knew already for Cs and Ce; Cs transports as solute that sorbs in matrix (after diffusion into matrix); Ce transports as colloid that is rapidly filtered; Sr not modeled, but qualitatively appears consistent with combination of colloid and solute transport

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

The Influence of Temperature on the Cyclic Properties of the Transversely Isotropic Nanocomposite System Under Kinematic Harmonic Loading

A micromechanical model is developed to determine the effective inelastic properties of nanocomposites under monoharmonic deformation by taking into account the detailed microstructural geometries and constitutive models of the constituents. By using the correspondence principle in viscoelasticity and the modified Mori–Tanaka method, the effects of interface between the inclusion and the matrix are taken into account. By applying the developed model, we perform the numerical analysis aimed at the determination of complex moduli for a polymeric nanocomposite reinforced with nanofibers composed from carbon nanotubes under the isothermal conditions. The dependences of complex moduli on temperature and the amplitude of strain intensity are analyzed. The composites reinforced with unidirectionally aligned nanofibers are considered. The accumulated results demonstrate a strong dependence of the storage and loss moduli on temperature within a broad temperature range. The storage and loss moduli are found to monotonically increase with the volume fraction of nanofibers. At the same time, they decrease as temperature increases. The obtained results show that the strength of material decreases as temperature increases in the elastic and inelastic regions and the inelastic behavior occurs for lower strain amplitude as temperature increases.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Effect of Nucleation Heterogeneity on Mineral Precipitation in Confined Environments

Abstract The formation of new mineral phases in confined environments, especially in porous media, is crucial for various geological processes like mineralization and diagenesis. The nucleation and precipitation of minerals are initiated at the microscale through fluid‐rock interaction, where dissolution of primary phases leads to supersaturated conditions and nucleation and growth of secondary ones. Previous research has focused primarily on either precipitation or nucleation, without fully exploring their combined impact. Our study introduces a computational framework that integrates classical nucleation theory with the micro‐continuum method. We validated our model by comparing with experiments, and discovered that different surface nucleation rate changes the mode of precipitation from a preferential to uniform precipitate textures. Furthermore, our study uncovered that the conventional deterministic precipitation method tends to underestimate the permeability of the porous matrix. In contrast, the new framework significantly improves model accuracy by incorporating preferential precipitation and heterogeneous nucleation.

58 GEOSCIENCES↗
DOE OSTI2021

Approaching the self-dual point of the sinh-Gordon model

One of the most striking but mysterious properties of the sinh-Gordon model (ShG) is the b → 1/b self-duality of its S-matrix, of which there is no trace in its Lagrangian formulation. Here b is the coupling appearing in the model’s eponymous hyperbolic cosine present in its Lagrangian, cosh(bΦ). In this paper we develop truncated spectrum methods (TSMs) for studying the sinh-Gordon model at a finite volume as we vary the coupling constant. We obtain the expected results for b << 1 and intermediate values of b, but as the self-dual point b = 1 is approached, the basic application of the TSM to the ShG breaks down. We find that the TSM gives results with a strong cutoff Ec dependence, which disappears according only to a very slow power law in Ec. Standard renormalization group strategies — whether they be numerical or analytic — also fail to improve upon matters here. We thus explore three strategies to address the basic limitations of the TSM in the vicinity of b = 1. In the first, we focus on the small-volume spectrum. We attempt to understand how much of the physics of the ShG is encoded in the zero mode part of its Hamiltonian, in essence how ‘quantum mechanical’ vs ‘quantum field theoretic’ the problem is. In the second, we identify the divergencies present in perturbation theory and perform their resummation using a supra-Borel approximate. In the third approach, we use the exact form factors of the model to treat the ShG at one value of b as a perturbation of a ShG at a different coupling. In the light of this work, we argue that the strong coupling phase b > 1 of the Lagrangian formulation of model may be different from what is naïvely inferred from its S-matrix. In particular, we present an argument that the theory is massless for b > 1.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Bayesian chain graph models to characterize microbe-environment dynamics

Microbiome data require statistical models that can simultaneously decode microbes' reaction to the environment and interactions among microbes. While a multiresponse linear regression model seems like a straight-forward solution, we argue that treating it as a graphical model is problematic given that the regression coefficient matrix does not encode the conditional dependence structure between response and predictor nodes. This observation is especially important in biological settings when we have prior knowledge on the edges from specific experimental interventions that can only be properly encoded under a conditional dependence model. Here, we propose a chain graph model with two sets of nodes (predictors and responses) whose solution yields a graph with edges that indeed represent conditional dependence, thus agreeing with the experimenter's intuition on the average behavior of nodes under treatment. The solution to our model is sparse via the Bayesian linear regression (LASSO). In addition, we propose an adaptive extension so that different shrinkages can be applied to different edges to incorporate edge-specific prior knowledge. Our model is computationally inexpensive through an efficient Gibbs sampling algorithm and can account for binary, counting, and compositional responses via an appropriate hierarchical structure. We test the performance of our model in a variety of simulated datasets, thereby showing superior performance to state-of-the-art approaches. We further apply our model to human gut and soil microbial compositional datasets, and we highlight that CG-LASSO can estimate biologically meaningful network structures in the data.

compositional data↗
DOE OSTI2016

Theory of optical transitions in π-conjugated macrocycles

We describe a theoretical and computational investigation of the optical properties of π-conjugated macrocycles. Since the low-energy excitations of these systems are Frenkel excitons that couple to high-frequency dispersionless phonons, we employ the quantized Frenkel-Holstein model and solve it via the density matrix renormalization group (DMRG) method. First we consider optical emission from perfectly circular systems. Owing to optical selection rules, such systems radiate via two mechanisms: (i) within the Condon approximation, by thermally induced emission from the optically allowed j = ± 1 states and (ii) beyond the Condon approximation, by emission from the j = 0 state via coupling with a totally non-symmetric phonon (namely, the Herzberg-Teller effect). Using perturbation theory, we derive an expression for the Herzberg-Teller correction and show via DMRG calculations that this expression soon fails as ħ ω/J and the size of the macrocycle increase. Next, we consider the role of broken symmetry caused by torsional disorder. In this case the quantum number j no longer labels eigenstates of angular momentum, but instead labels localized local exciton groundstates (LEGSs) or quasi-extended states (QEESs). As for linear polymers, LEGSs define chromophores, with the higher energy QEESs being extended over numerous LEGSs. Within the Condon approximation (i.e., neglecting the Herzberg-Teller correction) we show that increased disorder increases the emissive optical intensity, because all the LEGSs are optically active. We next consider the combined role of broken symmetry and curvature, by explicitly evaluating the Herzberg-Teller correction in disordered systems via the DMRG method. The Herzberg-Teller correction is most evident in the emission intensity ratio, I{sub 00}/I{sub 01}. In the Condon approximation I{sub 00}/I{sub 01} is a constant function of curvature, whereas in practice it vanishes for closed rings and only approaches a constant in the limit of vanishing curvature. We calculate the optical spectra of a model system, cyclo-poly(para-phenylene ethynylene), for different amounts of torsional disorder within and beyond the Condon approximation. We show how broken symmetry and the Herzberg-Teller effect explain the spectral features. The Herzberg-Teller correction to the 0–1 emission vibronic peak is always significant. Finally, we note the qualitative similarities between the optical properties of conformationally disordered linear polymers and macrocycles in the limit of sufficiently large disorder, because in both cases they are determined by the optical properties of curved chromophores.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Fano Resonance by Means of the One-Dimensional Superconductor Photonic Crystals

In the present work, we introduce a study about the reflectance properties of a one-dimensional superconductor photonic crystal in a different manner. Our structure is designed from alternately layers of high-temperature superconducting material (BSCCO) and a conventional superconducting material (Nb) terminated by the dielectric cap layer. The investigated numerical results are essentially based on the basics of the transfer matrix method and the two-fluid model. The numerical results showed the appearance of two sharp peaks, which are referred to Fano and electromagnetic-induced reflectance resonances. Also, we demonstrate the effects of many parameters on the reflectance properties of Fano resonance and the electromagnetic-induced reflectance. The effects of the constituent layer thicknesses, the number of periods, the operating temperature, the angle of incidence, and the refractive index of the dielectric cap layer are all considered. Finally, the effects of hydrostatic pressure on the Fano and EIR resonances are studied and calculated. This structure can play an important role in optical switching devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Optimal system size for complex dynamics in random neural networks near criticality

In this article, we consider a model of dynamical agents coupled through a random connectivity matrix, as introduced by Sompolinsky et al. [Phys. Rev. Lett. 61(3), 259–262 (1988)] in the context of random neural networks. When system size is infinite, it is known that increasing the disorder parameter induces a phase transition leading to chaotic dynamics. We observe and investigate here a novel phenomenon in the sub-critical regime for finite size systems: the probability of observing complex dynamics is maximal for an intermediate system size when the disorder is close enough to criticality. We give a more general explanation of this type of system size resonance in the framework of extreme values theory for eigenvalues of random matrices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗