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 793 records · Page 44

Non-robust Phase Transitions in the Generalized Clock Model on Trees

Pemantle and Steif provided a sharp threshold for the existence of a robust phase transition (RPT) for the continuous rotator model and the Potts model in terms of the branching number and the second eigenvalue of the transfer matrix whose kernel describes the nearest neighbor interaction along the edges of the tree. Here a RPT is said to occur if an arbitrarily weak coupling with symmetry-breaking boundary conditions suffices to induce symmetry breaking in the bulk. They further showed that for the Potts model RPT occurs at a different threshold than PT (phase transition in the sense of multiple Gibbs measures), and conjectured that RPT and PT should occur at the same threshold in the continuous rotator model. We consider the class of four- and five-state rotation-invariant spin models with reflection symmetry on general trees which contains the Potts model and the clock model with scalarproduct-interaction as limiting cases. The clock model can be viewed as a particular discretization which is obtained from the classical rotator model with state space S{sup 1}. We analyze the transition between PT=RPT and PT≠RPT, in terms of the eigenvalues of the transfer matrix of the model at the critical threshold value for the existence of RPT. The transition between the two regimes depends sensitively on the third largest eigenvalue.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

FEA-aided investigation of the effective thermal conductivity in a medium with embedded spheres

For multiple applications in nuclear energy, the ability to accurately represent material behavior with a simplified model is important to facilitate practical engineering-scale simulations. Here, we focus on the homogenized thermal response of a medium containing spherical inclusions, similar to a fuel form (compact or pebble) containing (TRISO) particles. A review on effective thermal conductivity (ETC) modeling is performed considering a random distribution of mono-sized spherical inclusions in a continuous matrix, with a primary focus on the analytical models. Finite element simulations are performed to evaluate each analytical model in varying material conditions. The model predictions are compared with the expected results obtained from the finite element predictions in addition to the Wiener and Hashin-Shtrikman bounds. Lastly, we included a practical discussion on the homogenization applied to a TRISO fuel pebble.

42 ENGINEERING↗

Mineral reaction kinetics constrain the length scale of rock matrix diffusion

Mass transport by aqueous fluids is a dynamic process in shallow crustal systems, redistributing nutrients as well as contaminants. Rock matrix diffusion into fractures (void space) within crystalline rock has been postulated to play an important role in the transient storage of solutes. The reacted volume of host rock involved, however, will be controlled by fluid-rock reactions. Here we present the results of a study which focusses on defining the length scale over which rock matrix diffusion operates within crystalline rock over timescales that are relevant to safety assessment of radioactive and other long-lived wastes. Through detailed chemical and structural analysis of natural specimens sampled at depth from an active system (Toki Granite, Japan), we show that, contrary to commonly proposed models, the length scale of rock matrix diffusion may be extremely small, on the order of centimetres, even over timescales of millions of years. This implies that in many cases the importance of rock matrix diffusion will be minimal. Additional analyses of a contrasting crystalline rock system (Carnmenellis Granite, UK) corroborate these results.

58 GEOSCIENCES↗

Quaternary AlInGaN/InGaN quantum well on vicinal c-plane substrate for high emission intensity of green wavelengths

Electronic and optical properties of non-trivial semipolar AlInGaN/InGaN quantum well (QW) structures are investigated by using the multiband effective-mass theory and non-Markovian optical model. On vicinal c-plane GaN substrate miscut by a small angle (θ < 40°) from c-plane, the AlInGaN/InGaN system is shown to have ∼3 times larger spontaneous emission peak intensity than the conventional InGaN/GaN system at green wavelength. It is attributed to much larger optical matrix element of the quaternary AlInGaN/InGaN system, derived from the reduction of internal electric field induced by polarizations. This effect exceeds the performance-degrading factor of smaller quasi-Fermi-level separation for the quaternary AlInGaN/InGaN system than that for the conventional InGaN/GaN system. Results indicate that the use of quaternary III-nitride QWs on vicinal substrates may be beneficial in improving the performance of optical devices emitting green light.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Transport properties of alkali metal doped fullerides

We have studied the intercage interactions between the adjacent C{sub 60} cages and expansion of lattice due to the intercalation of alkali atoms based on the spring model to estimate phonon frequencies from the dynamical matrix for the intermolecular alkali-C{sub 60} phonons. We considered a two-peak model for the phonon density of states to investigate the nature of electron pairing mechanism for superconducting state in fullerides. Coulomb repulsive parameter and the electron phonon coupling strength are obtained within the random phase approximation. Transition temperature, T{sub c}, is obtained in a situation when the free electrons in lowest molecular orbital are coupled with alkali-C{sub 60} phonons as 5 K, which is much lower as compared to reported T{sub c} (20 K). The superconducting pairing is mainly driven by the high frequency intramolecular phonons and their effects enhance it to 22 K. The importance of the present study, the pressure effect and normal state transport properties are calculated within the same model leading superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Stripe order enhanced superconductivity in the Hubbard model

Unidirectional (“stripe”) charge density wave order has now been established as a ubiquitous feature in the phase diagram of the cuprate high-temperature superconductors, where it generally competes with superconductivity. Nonetheless, on theoretical grounds it has been conjectured that stripe order (or other forms of “optimal” inhomogeneity) may play an essential positive role in the mechanism of high-temperature superconductivity. Here, we report density matrix renormalization group studies of the Hubbard model on long four- and six-leg cylinders, where the hopping matrix elements transverse to the long direction are periodically modulated—mimicking the effect of putative period 2 stripe order. We find that even modest amplitude modulations can enhance the long-distance superconducting correlations by many orders of magnitude and drive the system into a phase with a substantial spin gap and superconducting quasi–long-range order with a Luttinger exponent, K sc ~ 1.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Electronic specific heat capacities and entropies from density matrix quantum Monte Carlo using Gaussian process regression to find gradients of noisy data

In this work, we present a machine learning approach to calculating electronic specific heat capacities for a variety of benchmark molecular systems. Our models are based on data from density matrix quantum Monte Carlo, which is a stochastic method that can calculate the electronic energy at finite temperature. As these energies typically have noise, numerical derivatives of the energy can be challenging to find reliably. In order to circumvent this problem, we use Gaussian process regression to model the energy and use analytical derivatives to produce the specific heat capacity. From there, we also calculate the entropy by numerical integration. We compare our results to cubic splines and finite differences in a variety of molecules in which Hamiltonians can be diagonalized exactly with full configuration interaction. We finally apply this method to look at larger molecules where exact diagonalization is not possible and make comparisons with more approximate ways to calculate the specific heat capacity and entropy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Two-component quantum Hall effects in topological flat bands

Here in this paper, we study quantum Hall states for two-component particles (hardcore bosons and fermions) loading in topological lattice models. By tuning the interplay of interspecies and intraspecies interactions, we demonstrate that two-component fractional quantum Hall states emerge at certain fractional filling factors ν = 1/2 for fermions (ν = 2/3 for bosons) in the lowest Chern band, classified by features from ground states including the unique Chern number matrix (inverse of the K matrix), the fractional charge and spin pumpings, and two parallel propagating edge modes. Moreover, we also apply our strategy to two-component fermions at integer filling factor ν = 2 , where a possible topological Neel antiferromagnetic phase is under intense debate very recently. For the typical π -flux checkerboard lattice, by tuning the onsite Hubbard repulsion, we establish a first-order phase transition directly from a two-component fermionic ν = 2 quantum Hall state at weak interaction to a topologically trivial antiferromagnetic insulator at strong interaction, and therefore exclude the possibility of an intermediate topological phase for our system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Renormalization of the Standard Model Effective Field Theory from geometry

S-matrix elements are invariant under field redefinitions of the Lagrangian. They are determined by geometric quantities such as the curvature of the field-space manifold of scalar and gauge fields. We present a formalism where scalar and gauge fields are treated together, with a metric on the combined space of both types of fields. Scalar and gauge scattering amplitudes are given by the Riemann curvature Rijkl of this combined space, with indices i, j, k, l chosen to be scalar or gauge indices depending on the type of external particle. One-loop divergences can also be computed in terms of geometric invariants of the combined space, which greatly simplifies the computation of renormalization group equations. We apply our formalism to the Standard Model Effective Field Theory (SMEFT), and compute the renormalization group equations for even-parity bosonic operators to mass dimension eight.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Phase diagram of the one-dimensional t 1 - t 2 - J model: Ferromagnetism, triplet pairing, and charge and pair density waves

Here, we present a density matrix renormalization group (DMRG) study of an extended t-J model with hopping to the first and second neighbors -- the one dimensional t 1 -t 2 -J model. The full phase diagram as a function of the density n and exchange strength J, for both positive and negative values of t 2 , is obtained. For t 2 = -0.5 we observe that, in the strongly interacting region, Nagaoka Ferromagnetism (FM) is accompanied by a triplet pair density wave (PDW) upon doping. As the spin exchange J increases, a charge density wave (CDW) phase emerges and then gives way to singlet superconductivity (SC). This phase behaves as a singlet PDW with vanishing spacial average of the order parameter and a spin gap. When t 2 = 0.5, the physics is basically reminiscent of the conventional t-J model, undergoing a transition from a metallic to a SC phase as a function of J.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Magnetic properties and pairing tendencies of the iron-based superconducting ladder BaFe 2 S 3 : Combined ab initio and density matrix renormalization group study

The recent discovery of superconductivity under high pressure in the two-leg ladder compound BaFe 2 S 3 [H. Takahashi et al., Nat. Mater. 14, 1008 (2015)] opens a broad avenue of research, because it represents the first report of pairing tendencies in a quasi-one-dimensional iron-based high-critical-temperature superconductor. Similarly, as in the case of the cuprates, ladders and chains can be far more accurately studied using many-body techniques and model Hamiltonians than their layered counterparts, particularly if several orbitals are active. In this publication, we derive a two-orbital Hubbard model from first principles that describes individual ladders of BaFe 2 S 3 . The model is studied with the density matrix renormalization group. These first reported results are exciting for two reasons: (i) at half-filling, ferromagnetic order emerges as the dominant magnetic pattern along the rungs of the ladder, and antiferromagnetic order along the legs, in excellent agreement with neutron experiments; and (ii) with hole doping, pairs form in the strong coupling regime, as found by studying the binding energy of two holes doped on the half-filled system. In addition, orbital selective Mott phase characteristics develop with doping, with only oneWannier orbital receiving the hole carriers while the other remains half-filled. Lastly, these results suggest that the analysis of models for iron-based two-leg ladders could clarify the origin of pairing tendencies and other exotic properties of iron-based high-critical-temperature superconductors in general.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Full-dimensional diabatic potential energy surfaces including dissociation: The {sup 2}E{sup ″} state of NO{sub 3}

A scheme to produce accurate full-dimensional coupled diabatic potential energy surfaces including dissociative regions and suitable for dynamical calculations is proposed. The scheme is successfully applied to model the two-sheeted surface of the {sup 2}E{sup ″} state of the NO{sub 3} radical. An accurate potential energy surface for the NO{sub 3}{sup −} anion ground state is developed as well. Both surfaces are based on high-level ab initio calculations. The model consists of a diabatic potential matrix, which is expanded to higher order in terms of symmetry polynomials of symmetry coordinates. The choice of coordinates is key for the accuracy of the obtained potential energy surfaces and is discussed in detail. A second central aspect is the generation of reference data to fit the expansion coefficients of the model for which a stochastic approach is proposed. A third ingredient is a new and simple scheme to handle problematic regions of the potential energy surfaces, resulting from the massive undersampling by the reference data unavoidable for high-dimensional problems. The final analytical diabatic surfaces are used to compute the lowest vibrational levels of NO{sub 3}{sup −} and the photo-electron detachment spectrum of NO{sub 3}{sup −} leading to the neutral radical in the {sup 2}E{sup ″} state by full dimensional multi-surface wave-packet propagation for NO{sub 3} performed using the Multi-Configuration Time Dependent Hartree method. The achieved agreement of the simulations with available experimental data demonstrates the power of the proposed scheme and the high quality of the obtained potential energy surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Coreflood Tests to Evaluate Enhanced Oil Recovery Potential of Wettability-Altering Surfactants for Oil-Wet Heterogeneous Carbonate Reservoirs

Summary Oil-wetness and heterogeneity are two main factors that result in low oil recovery (OR) by waterflood in carbonate reservoirs. The injected water is likely to flow through high-permeability regions and bypass the oil in the low-permeability matrix. In this study, systematic coreflood tests were carried out in both “homogeneous” cores and “heterogeneous” cores with a wettability-altering surfactant. The homogeneous coreflood tests were conducted to evaluate surfactant retention, as well as to compare tertiary surfactant flooding with secondary surfactant flooding. The heterogeneous coreflood test was proposed to model bypassing in low-permeability matrix during waterfloods, and dynamic imbibition of surfactant into the low-permeability matrix. Surfactant retention results suggest that retention increases as initial oil saturation decreases. The retention of selected surfactant in the target reservoir cores was measured to be within a range of 0.07–0.12 mg/g-rock, which is economically acceptable. The results of homogeneous coreflood tests showed that both secondary waterflood and secondary surfactant flood can achieve high OR (>50%) from relatively homogeneous oil-wet cores. A shut-in phase after the surfactant injection resulted in a surge in oil production, which suggests that enough time should be given for wettability alteration by surfactants. The results of heterogeneous coreflood tests showed that more oil is bypassed in the tighter matrix by waterflood if the permeability is higher in the flooded layer and this bypassed oil is the target for the wettability-altering surfactant floods. Slow wettability-altering surfactant injection leads to imbibition into bypassed regions. When the oil-wet carbonate reservoirs have large unswept regions after waterflood, wettability-altering surfactants can significantly improve OR if enough time is given for imbibition.

Engineering↗

Inhibition of TGF-β pathway reverts extracellular matrix remodeling in T. cruzi-infected cardiac spheroids

Highlights: • A model of cardiac spheroids to study fibrosis in pathological processes is proposed. • The study evaluates a TGF-β inhibitor in restoring extracellular matrix balance disrupted by T. cruzi infection. • Inhibition of TGF-β signaling reduces fibrosis and hypertrophy in infected-cardiac spheroids. • The extracellular matrix remodeling is achieved through regulation of MMP-2 activity and TIMP-1 expression. Chagasic cardiomyopathy (CC) is the main manifestation of Chagas Disease (CD). CC is a progressive dysfunctional illness, in which transforming growth factor beta (TGF-β) plays a central role in fibrogenesis and hypertrophy. In the present study, we tested in a three-dimensional (3D) model of cardiac cells culture (named cardiac spheroids), capable of mimicking the aspects of fibrosis and hypertrophy observed in CC, the role of TGF-β pathway inhibition in restoring extracellular matrix (ECM) balance disrupted by T. cruzi infection. Treatment of T. cruzi-infected cardiac spheroids with SB 431542, a selective inhibitor of TGF-β type I receptor, resulted in a reduction in the size of spheroids, which was accompanied by a decrease in parasite load and in fibronectin expression. The inhibition of TGF-β pathway also promoted an increase in the activity of matrix metalloproteinase (MMP)-2 and a decrease in tissue inhibitor of matrix metalloproteinase (TIMP)-1 expression, which may be one of the mechanisms regulating extracellular matrix remodeling. Therefore, our study provides new insights into the molecular mechanisms by which inhibition of TGF-β signaling reverts fibrosis and hypertrophy generated by T. cruzi during CC and also highlights the use of cardiac spheroids as a valuable tool for the study of fibrogenesis and anti-fibrotic compounds.

60 APPLIED LIFE SCIENCES↗

The Influence of Shallow Subsurface Properties on Particle Motion in Acoustic-Seismic Coupling

Atmospheric acoustic waves transmit energy into the solid Earth through air-to-ground coupling. These waves are recorded by seismic sensors and provide insight into both atmospheric phenomena and subsurface properties. Interpreting these signals is often challenging because they are modulated by subsurface structure and the incidence angle of the acoustic wave. This study examines acoustic--seismic coupling generated by the 2012 Camp Minden Explosion, which was recorded by hundreds of seismoacoustic stations. We apply a novel technique to quantify the seismic particle motion, model coupled waves with a propagator matrix approach, and apply a Bayesian inversion to infer properties of the shallow subsurface. Our analysis reveals that prograde motion is widespread and focused in low shear-wave velocity regions, such as the Mississippi Embayment, and retrograde motion is more common in higher shear-velocity areas. Inversion results at some stations produce plausible subsurface models with strong waveform fits, while inversions at other sites are less successful. These results indicate prograde particle motion in air-to-ground coupled waves is more prevalent than previously recognized and may serve as a diagnostic for shallow velocity structure. Our comprehensive modeling and inversion framework provides a potential method to extract layered near-surface properties from acoustic-seismic coupling observations.

58 GEOSCIENCES↗

Production of heavy sterile neutrinos from vector boson decay at electroweak temperatures

Here, in the standard model extended with a seesaw mass matrix, we study the production of sterile neutrinos from the decay of vector bosons at temperatures near the masses of the electroweak bosons. We derive a general quantum kinetic equation for the production of sterile neutrinos and their effective mixing angles, which is applicable over a wide range of temperature, to all orders in interactions of the standard model and to leading order in a small mixing angle for the neutrinos. We emphasize the relation between the production rate and Landau damping at one-loop order and show that production rates and effective mixing angles depend sensitively upon the neutrino’s helicity. Sterile neutrinos with positive helicity interact more weakly with the medium than those with negative helicity, and their effective mixing angle is not modified significantly. Negative helicity states couple more strongly to the vector bosons, but their mixing angle is strongly suppressed by the medium. Consequently, if the mass of the sterile neutrino is ≲ 8.35 MeV , there are fewer states with negative helicity produced than those with positive helicity. There is an Mikheyev-Smirnov-Wolfenstein-type resonance in the absence of lepton asymmetry, but due to screening by the damping rate, the production rate is not enhanced. Sterile neutrinos with negative helicity freeze out at T$-\atop{f}$ ≃ 5 GeV , whereas positive helicity neutrinos freeze out at T$+\atop{f}$≃ 8 GeV , with both distributions far from thermal. As the temperature decreases, due to competition between a decreasing production rate and an increasing mixing angle, the distribution function for states with negative helicity is broader in momentum and hotter than that for those with positive helicity. Sterile neutrinos produced via vector boson decay do not satisfy the abundance, lifetime, and cosmological constraints to be the sole dark matter component in the Universe. Massive sterile neutrinos produced via vector boson decay might solve the 7 Li problem, albeit at the very edge of the possible parameter space. A heavy sterile neutrino with a mass of a few MeV could decay into light sterile neutrinos, of a few keV in mass, that contribute to warm dark matter. In conclusion, we argue that heavy sterile neutrinos with lifetime ≤1/H 0 reach local thermodynamic equilibrium.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

How does a probe inserted into the discharge influence the plasma structure?

Shielding the bias applied to the probe by the sheath formed around it and determination of parameters of unperturbed plasmas are in the basis of the probe diagnostics. The results from a two-dimensional model of a discharge with a probe inserted in it show that the probe influences the spatial distribution of the plasma parameters in the entire discharge. The increase (although slight) in the electron temperature, due to the increased losses of charged particles on the additional wall in the discharge (mainly the probe holder), leads to redistribution of the plasma density and plasma potential, as shown by the results obtained at the floating potential of the probe. The deviations due to the bias applied to the probe tip are stronger in the ion saturation region of the probe characteristics. The pattern of the spatial redistribution of the plasma parameters advances together with the movement of the probe deeper in the discharge. Although probe sheaths and probe characteristics resulting from the model are shown, the study does not aim at discussions on the theories for determination of the plasma density from the ion saturation current. Regardless of the modifications in the plasma behavior in the entire discharge, the deviations of the plasma parameters at the position of the probe tip and, respectively, the uncertainty which should be added as an error when the accuracy of the probe diagnostics is estimated do not exceed 10%. Consequently, the electron density and temperature obtained, respectively, at the position of the plasma potential on the probe characteristics and from its transition region are in reasonable agreement with the results from the model of the discharge without a probe. Being in the scope of research on a source of negative hydrogen ions with the design of a matrix of small radius inductive discharges, the model is specified for a low-pressure hydrogen discharge sustained in a small-radius tube.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗