Stimulus induced strain in spin transition heterostructures
Explore the source record for details and available documents.
99 GENERAL AND MISCELLANEOUS↗
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.
Band-mixing effects in one-dimensional charge transfer insulators
The low-energy properties of transition metal oxides (TMOs) are governed by the electrons occupying strongly correlated $d$-orbitals that are hybridized with surrounding ligand oxygen $p$ orbitals to varying degrees. Their physics is thus established by a complex interplay between the transition-metal (TM)-ligand hopping $t$, charge transfer energy $\Delta_\mathrm{CT}$, and on-site TM Hubbard repulsion $U$. Here, we study the spectral properties of a one-dimensional (1D) analog of such a $pd$ system, with alternating TM $d$ and ligand anion $p$ orbitals situated along a chain. Using the density matrix renormalization group method, we study the model's single-particle spectral function, x-ray absorption spectrum, and dynamical spin structure factor as a function of $\Delta_\mathrm{CT}$ and $U$. In particular, we present results spanning from the Mott insulating ($\Delta_\mathrm{CT} > U$) to negative charge transfer regime $\Delta_\mathrm{CT} < 0$ to gain a better understanding of the ground and momentum-resolved excited state properties of these different regimes. Our results can guide new studies on TMOs that seek to situate them within the Mott-Hubbard/charge transfer insulator classification scheme.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Revisiting the three-kaon interaction and its relation with 𝐾(1460)
Here, we test the hypothesis of 𝐾(1460) being a hadronic 𝐽 𝑃 = 0 − molecule through nonperturbative 𝑆-wave $KK\bar{K}$, 𝐾𝜋𝜋, 𝐾𝜋𝜂 coupled-channel dynamics in a three-body unitary isobar approach. The scalar two-body coupled-channel resonances 𝑓 0 (500), 𝑓 0 (980), 𝑎 0 (980), and 𝐾$^*_0$(700) are generated in the subsytems in amplitudes that match phase shifts from experiment. In a first step, previous results in the limit of mass-degenerate, stable 𝑓 0 and 𝑎 0 isobars are reproduced. Once the full seven-coupled channel model is switched on, other 𝑆-matrix effects obscure and modify the resonance signal, including complex thresholds, a two-body cusp, and a triangle singularity at threshold.
Döring, M. [George Washington Univ., Washington, D↗
Third moments of nucleon unpolarized, polarized, and transversity parton distribution functions from physical-point lattice QCD
Using forward matrix elements of local leading-twist operators, we present a determination of the isovector third Mellin moments ⟨𝑥 2 ⟩ of nucleon unpolarized, polarized, and transversity parton distribution functions. Two lattice QCD ensembles at the physical pion mass are used, which were generated using a tree-level Symanzik-improved gauge action and 2+1 flavor tree-level improved Wilson Clover fermions coupling via 2-level HEX-smearing. Leveraging a wide set of operators, two extraction methods for the matrix elements, and the automatic inclusion of model uncertainties via bootstrapped model averages, we extract values of the third Mellin moments. Furthermore, this is the first direct calculation of these observables performed at the physical pion mass.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Ultrastrong coupling of vibrationally dressed organic Frenkel excitons with Bloch surface waves in a one-sided all-dielectric structure
We demonstrate a transition from weak, to strong, to ultrastrong coupling of Frenkel molecular excitons and Bloch surface wave photons at room temperature using a one-sided, all-dielectric optical structure. The all-dielectric structure comprises an organic semiconductor thin film on the surface of a distributed Bragg reflector. We investigate the evolution of multiple vibronic polariton branches and their dominant absorption peaks as a function of coupling and in-plane momentum, which is absent in previous ultrastrong coupling systems. Measurements are interpreted using both the transfer matrix method and a coupled-oscillator model without the rotating wave approximation. The dependence of Rabi splitting on the number of excitons and electrical field amplitude is also modeled showing a transition to ultrastrong coupling at film thicknesses ≥ 50 nm. Furthermore, this low-loss polaritonic structure enables us to study phenomena such as organic exciton-polariton dynamics, ultralong range polariton propagation, and high efficiency energy transport in the ultrastrong coupling regime at room temperature.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Quadrupolar interactions between acceptor pairs in $\mathcal{p}$-doped semiconductors
Here, we consider the interaction between acceptor pairs in doped semiconductors in the limit of large interacceptor separation relevant for low doping densities. Modeling individual acceptors via the spherical model of Baldereschi and Lipari, we calculate matrix elements of the quadrupole tensor between the four degenerate ground states and show that the acceptor has a nonzero quadrupole moment. As a result, the dominant contribution to the large-separation acceptor-acceptor interaction comes from direct (charge-density) terms rather than exchange terms. The quadrupole is the leading nonzero moment, so the electric quadrupole-quadrupole interaction dominates for large separation. We calculate the matrix elements of the quadrupole-quadrupole interaction Hamiltonian in a product-state basis and diagonalize, obtaining a closed-form expression for the energies and degeneracies of the sixteen-state energy spectrum. All dependence on material parameters enters via an overall prefactor, resulting in surprisingly simple and universal results. This simplicity is due, in part, to a mathematical happenstance, the nontrivial vanishing of a particular Wigner 6-$\mathcal{j}$ symbol, $^\{\frac{^2_3}{\small2} \frac{^2_3}{\small2} \frac{^2_3}{\small2}$$^\}$ = 0. Results are relevant to the control of two-qubit interactions in quantum computing implementations based on acceptor spins, as well as calculations of the thermodynamic properties of insulating $\mathcal{p}$-type semiconductors.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Coupled dimer and bond-order-wave order in the quarter-filled one-dimensional Kondo lattice model
Motivated by experiments on the organic compound (Per) 2 [Pt(mnt) 2 ], we study the ground state of the onedimensional Kondo lattice model at quarter filling with the density matrix renormalization group method. Here, we show a coupled dimer and bond-order-wave (BOW) state in the weak-coupling regime for the localized spins and itinerant electrons, respectively. The quantum phase transitions for the dimer and the BOW orders occur at the same critical coupling parameter J c , with the opening of a charge gap. The emergence of the combination of dimer and BOW order agrees with the experimental findings of the simultaneous Peierls and spin-Peierls transitions at low temperatures, which provides a theoretical understanding of such a phase transition. We also show that the localized spins in this insulating state have quasi-long-ranged spin correlations with collinear configurations, which resemble the classical dimer order in the absence of a magnetic order.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Unified ab initio description of Fröhlich electron-phonon interactions in two-dimensional and three-dimensional materials
Ab initio calculations of electron-phonon interactions including the polar Fröhlich coupling have advanced considerably in recent years. The Fröhlich electron-phonon matrix element is by now well understood in the case of bulk three-dimensional (3D) materials. In the case of two-dimensional (2D) materials, the standard procedure to include Fröhlich coupling is to employ Coulomb truncation, so as to eliminate artificial interactions between periodic images of the 2D layer. While these techniques are well established, the transition of the Fröhlich coupling from three to two dimensions has not been investigated. Furthermore, it remains unclear what error one makes when describing 2D systems using the standard bulk formalism in a periodic supercell geometry. In this work, we generalize previous work on the ab initio Fröhlich electron-phonon matrix element in bulk materials by investigating the electrostatic potential of atomic dipoles in a periodic supercell consisting of a 2D material and a continuum dielectric slab. We obtain a unified expression for the matrix element, which reduces to the existing formulas for three-dimensional and 2D systems when the interlayer separation tends to zero or infinity, respectively. This expression enables an accurate description of the Fröhlich matrix element in 2D systems without resorting to Coulomb truncation. We validate our approach by direct ab initio density-functional perturbation theory calculations for monolayer BN and MoS 2 , and we provide a simple expression for the 2D Fröhlich matrix element that can be used in model Hamiltonian approaches. The formalism outlined in this work may find applications in calculations of polarons, quasiparticle renormalization, transport coefficients, and superconductivity, in 2D and quasi-2D materials.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Topological and magnetic properties of the interacting Bernevig-Hughes-Zhang model
We investigate the effects of electronic correlations on the Bernevig-Hughes-Zhang model using the real-space density matrix renormalization group (DMRG) algorithm. We introduce a method to probe topological phase transitions in systems with strong correlations using DMRG, substantiated by an unsupervised machine learning methodology that analyzes the orbital structure of the real-space edges. Including the full multi-orbital Hubbard interaction term, we construct a phase diagram as a function of a gap parameter (m) and the Hubbard interaction strength (U) via exact DMRG simulations on N×4 cylinders. Our analysis confirms that the topological phase persists in the presence of interactions, consistent with previous studies, but it also reveals an intriguing phase transition from a paramagnetic to a stripey antiferromagnetic topological insulator. The combination of the magnetic structure factor, strength of magnetic moments, and the orbitally resolved density, provides real-space information on both topology and magnetism in a strongly correlated system.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Entanglement of elastic and inelastic scattering
The entanglement properties of systems in which elastic and inelastic reactions occur in projectile-target interactions is studied. A new measure of entanglement, the scattering entropy, based on the unitarity of the S matrix (probability conservation), is suggested. Using simple models for both low- and high-energy interactions, the amount of entanglement is found to track with the strength of the inelastic interaction. Here, the familiar example of the classical “black disk,” or total absorption, model is found to correspond to maximum entanglement. An analysis of high-energy pp scattering data shows that entanglement is near maximum for laboratory energies greater than about 1 GeV, showing that the total absorption model is a reasonable starting point for understanding the data.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Hadronic vacuum polarization for the muon 𝑔 − 2 from lattice QCD: Complete short and intermediate windows
We present complete results for the hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment 𝑎 𝜇 in the short- and intermediate-distance window regions, which account for roughly 10% and 35% of the total HVP contribution to 𝑎 𝜇 , respectively. In particular, we perform lattice-QCD calculations for the isospin-symmetric connected and disconnected contributions, as well as corrections due to strong-isospin breaking. For the short-distance window observables, we investigate the so-called log-enhancement effects as well as the significant oscillations associated with staggered quarks in this region. For the dominant, isospin-symmetric light-quark-connected contribution, we obtain 𝑎$^{𝑙𝑙,SD}_{𝜇}$(conn) = 48.139(11) stat (91) syst [92] total × 10 −10 and 𝑎$^{𝑙𝑙,W}_{𝜇}$(conn) = 206.90(14) stat (61) syst [63] total × 10 −10 . We use Bayesian model averaging to fully estimate the covariance matrix between the individual contributions. Our determinations of the complete window contributions are 𝑎$^{SD}_{𝜇}$ = 69.05(1) stat (21) syst [21] total × 10 −10 and 𝑎$^{W}_{𝜇}$ = 236.45(17) stat (83) syst [85] total × 10 −10 . This work is part of our ongoing effort to compute all contributions to HVP with an overall uncertainty at the few-permille level.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Power-Law Entanglement and Hilbert Space Fragmentation in Nonreciprocal Quantum Circuits
Quantum circuits utilizing measurement to evolve a quantum wave function offer a new and rich playground to engineer unconventional entanglement dynamics. Here, in this work, we introduce a hybrid, nonreciprocal setup featuring a quantum circuit, whose updates are conditioned on the state of a classical dynamical agent. In our example the circuit is represented by a Majorana quantum chain controlled by a classical N-state Potts chain undergoing pair flips. The local orientation of the classical spins controls whether randomly drawn local measurements on the quantum chain are allowed or not. This imposes a dynamical kinetic constraint on the entanglement growth, described by the transfer matrix of an N-colored loop model. It yields an equivalent description of the circuit by an SU(N)-symmetric Temperley-Lieb Hamiltonian or by a kinetically constrained surface growth model for an N-component height field. For N = 2, we find a diffusive growth of the half-chain entanglement toward a stationary profile S(L) ~ L 1/2 for L sites. For N ≥ 3, the kinetic constraints impose Hilbert space fragmentation, yielding subdiffusive growth toward S(L) ~ L 0.57 . This showcases how the control by a classical dynamical agent can enrich the entanglement dynamics in quantum circuits, paving a route toward novel entanglement dynamics in nonreciprocal hybrid circuit architectures.
1-dimensional spin chains↗
First-principles calculations of solute transport in zirconium: Vacancy-mediated diffusion with metastable states and interstitial diffusion
Zirconium alloys are the most widely used nuclear fuel cladding materials for light water power reactors where irradiation damage causes solute redistribution, leading to degradation of alloy properties such as corrosion resistance. Designing radiation-tolerant zirconium alloys requires a thorough understanding of the atomic-scale transport behavior of the alloying elements in Zr. In this work, we perform density function theory calculations to investigate the diffusion of Sn, Cr, Fe, Be, Al, and Ni in the hexagonal close-packed (HCP) Zr matrix. We develop a methodology to accurately model the metastable vacancy states along the basal migration path, known to occur in group IV metals. We compute the vacancy-mediated solute diffusion coefficients and drag ratios using the kinetic Monte Carlo method and an analytic Green's function method—the agreement between the two validates our methodology. The computed diffusion coefficients of Sn and Al show good agreement with the experimental data and we expect these solutes to diffuse via the vacancy-mediated mechanism. We use a Green's function approach, parameterized with data from density functional theory calculations, to compute the interstitial diffusion coefficients of Cr, Fe, Be, and Ni in the HCP Zr lattice. The computed diffusion coefficients of Cr, Ni, and Be agree with the experimental measurements within one order of magnitude, while those of Fe are within two orders of magnitude of the experimental measurements. The drag ratios for Cr, Fe, Be, and Ni are positive up to 1235 K, which suggests that nonequilibrium vacancy fluxes could drag these solutes toward sinks such as dislocation loops and grain boundaries. We also compute the transport coefficients without including the metastable states, and using the eight- and thirteen-frequency model. Our results show significant differences in drag ratio for the eight- and thirteen-frequency model predictions compared with the Green's function methodology, but smaller errors in the solute diffusivity. Combining interstitial and vacancy-mediated diffusivities, we predict the unusual result that increased vacancy concentration slows down solute diffusivity, while a sufficiently high vacancy concentration can change the dominant mechanism to an accelerated vacancy-mediated diffusion.
36 MATERIALS SCIENCE↗
In the grass species Brachypodium distachyon , the production of mixed-linkage (1,3;1,4)-β-glucan (MLG) occurs in the Golgi apparatus
Mixed-inkage (1,3;1,4)-β-glucan (MLG) is a glucose polymer with beneficial effects on human health and high potential for the agricultural industry. MLG is present predominantly in the cell wall of grasses and is synthesized by cellulose synthase-like F or H families of proteins, with CSLF6 being the best-characterized MLG synthase. Although the function of this enzyme in MLG production has been established, the site of MLG synthesis in the cell is debated. It has been proposed that MLG is synthesized at the plasma membrane, as occurs for cellulose and callose; in contrast, it has also been proposed that MLG is synthesized in the Golgi apparatus, as occurs for other matrix polysaccharides of the cell wall. Testing these conflicting possibilities is fundamentally important in the general understanding of the biosynthesis of the plant cell wall. Using immuno-localization analyses with MLG-specific antibody in Brachypodium and in barley, we found MLG present in the Golgi, in post-Golgi structures and in the cell wall. Accordingly, analyses of a functional fluorescent protein fusion of CSLF6 stably expressed in Brachypodium demonstrated that the enzyme is localized in the Golgi. We also established that overproduction of MLG causes developmental and growth defects in Brachypodium as also occur in barley. Our results indicated that MLG production occurs in the Golgi similarly to other cell wall matrix polysaccharides, and supports the broadly applicable model in grasses that tight mechanisms control optimal MLG accumulation in the cell wall during development and growth. This work addresses the fundamental question of where mixed linkage (1,3;1,4)-β-glucan (MLG) is synthesized in plant cells. By analyzing the subcellular localization of MLG and MLG synthase in an endogenous system, we demonstrated that MLG synthesis occurs at the Golgi in Brachypodium and barley. In conclusion, a growth inhibition due to overproduced MLG in Brachypodium supports the general applicability of the model that a tight control of the cell wall polysaccharides accumulation is needed to maintain growth homeostasis during development.
59 BASIC BIOLOGICAL SCIENCES↗
Characterization of scatter in digital mammography from physical measurements
Purpose: That scattered radiation negatively impacts the quality of medical radiographic imaging is well known. In mammography, even slight amounts of scatter reduce the high contrast required for subtle soft-tissue imaging. In current clinical mammography, image contrast is partially improved by use of an antiscatter grid. This form of scatter rejection comes with a sizeable dose penalty related to the concomitant elimination of valuable primary radiation. Digital mammography allows the use of image processing as a method of scatter correction that might avoid effects that negatively impact primary radiation, while potentially providing more contrast improvement than is currently possible with a grid. For this approach to be feasible, a detailed characterization of the scatter is needed. Previous research has modeled scatter as a constant background that serves as a DC bias across the imaging surface. The goal of this study was to provide a more substantive data set for characterizing the spatially-variant features of scatter radiation at the image detector of modern mammography units. Methods: This data set was acquired from a model of the radiation beam as a matrix of very narrow rays or pencil beams. As each pencil beam penetrates tissue, the pencil widens in a predictable manner due to the production of scatter. The resultant spreading of the pencil beam at the detector surface can be characterized by two parameters: mean radial extent (MRE) and scatter fraction (SF). The SF and MRE were calculated from measurements obtained using the beam stop method. Two digital mammography units were utilized, and the SF and MRE were found as functions of target, filter, tube potential, phantom thickness, and presence or absence of a grid. These values were then used to generate general equations allowing the SF and MRE to be calculated for any combination of the above parameters. Results: With a grid, the SF ranged from a minimum of about 0.05 to a maximum of about 0.16, and the MRE ranged from about 3 to 13 mm. Without a grid, the SF ranged from a minimum of 0.25 to a maximum of 0.52, and the MRE ranged from about 20 to 45 mm. The SF with a grid demonstrated a mild dependence on target/filter combination and kV, whereas the SF without a grid was independent of these factors. The MRE demonstrated a complex relationship as a function of kV, with notable difference among target/filter combinations. The primary source of change in both the SF and MRE was phantom thickness. Conclusions: Because breast tissue varies spatially in physical density and elemental content, the effective thickness of breast tissue varies spatially across the imaging field, resulting in a spatially-variant scatter distribution in the imaging field. The data generated in this study can be used to characterize the scatter contribution on a point-by-point basis, for a variety of different techniques.
62 RADIOLOGY AND NUCLEAR MEDICINE↗
Determination of the Disperse Composition of a PbO Suspension Containing Aggregates of Particles of Lamellar Shape by the Laser-Polarimetry Method
The results of measurements of the scattering matrix at a wavelength of 0.63 μm in the range of scattering angles of 10°–155° are presented for an aqueous suspension of lead oxide containing particles of plate form and their aggregates of monomers with dimensions of ~5 nm. The results of the measurements are compared with the results of calculations for axially symmetric scatterers (ellipsoids of rotation, cylinders). It is shown that the presence of aggregates affects the scattering properties of such a medium. The results of reconstructing the distribution of particles of a disperse medium in sizes from the measurements data of the scattering matrix are presented. The reconstruction of the distributions was carried out by solving the problem of optimizing the sum of the squared deviations of the experimental and calculated values of matrix elements in the framework of the model of axially symmetric scatterers. It is shown that the distribution of particles by sizes is more accurately reconstructed by minimizing the sum of the squares of the deviations for the sum of the diagonal elements. The obtained distribution is compared with the distribution measured by the method of dynamic light scattering.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Report on the Implementation of Homogeneous Nucleation Scheme in MARMOT-based Phase Field Simulation
In this report, we summarized our effort in developing mesoscale phase field models for predicting precipitation kinetics in alloys during thermal aging and/or under irradiation in nuclear reactors. The first part focused on developing a method to predict the thermodynamic properties of critical nuclei such as the sizes and concentration profiles of critical nuclei, and nucleation barrier. These properties are crucial for quantitative simulations of precipitate evolution kinetics with phase field models. Fe-Cr alloy was chosen as a model alloy because it has valid thermodynamic and kinetic data as well as it is an important structural material in nuclear reactors. A constrained shrinking dimer dynamics (CSDD) method was developed to search for the energy minimum path during nucleation. With the method we are able to predict the concentration profiles of the critical nuclei of Cr-rich precipitates and nucleation energy barriers. Simulations showed that Cr concentration distribution in the critical nucleus strongly depends on the overall Cr concentration as well as temperature. The Cr concentration inside the critical nucleus is much smaller than the equilibrium concentration calculated by the equilibrium phase diagram. This implies that a non-classical nucleation theory should be used to deal with the nucleation of Cr precipitates in Fe-Cr alloys. The growth kinetics of both classical and non-classical nuclei was investigated by the phase field approach. A number of interesting phenomena were observed from the simulations: 1) a critical classical nucleus first shrinks toward its non-classical nucleus and then grows; 2) a non-classical nucleus has much slower growth kinetics at its earlier growth stage compared to the diffusion-controlled growth kinetics. 3) a critical classical nucleus grows faster at the earlier growth stage than the non-classical nucleus. All of these results demonstrated that it is critical to introduce the correct critical nuclei into phase field modeling in order to correctly capture the kinetics of precipitation. In most alloys the matrix phase and precipitate phase have different concentrations as well as different crystal structures. For example, Cu precipitates in FeCu alloys have fcc crystal structure while the matrix Fe-Cu solid solution has bcc structure at low temperature. The WBM model and KimS model, where both concentrations and order parameters are chosen to describe the microstructures, are commonly used to model precipitations in such alloys. The WBM and KimS models have not been implemented into Marmot yet. In the second part of this report, we focused on implementing the WBM and KimS models into Marmot. The Fe-Cu alloys, which are important structure materials in nuclear reactors, was taken as the model alloys to test the models.
Fe Alloys↗
Feasibility of Geophysical Monitoring of Carbon-Sequestrated Deep Saline Aquifers
As carbon dioxide (CO{sub 2}) is sequestered from the bottom of a brine reservoir and allowed to migrate upward, the effects of the relative permeability hysteresis due to capillary trapping and buoyancy driven migration tend to make the reservoir patchy saturated with different fluid phases over time. Seismically, such a patchy saturated reservoir induces an effective anisotropic behavior whose properties are primarily dictated by the nature of the saturation of different fluid phases in the pores and the elastic properties of the rock matrix. By combining reservoir flow simulation and modeling with seismic modeling, it is possible to derive these effective anisotropic properties, which, in turn, could be related to the saturation of CO{sub 2} within the reservoir volume any time during the post-injection scenario. Therefore, if time-lapse seismic data are available and could be inverted for the effective anisotropic properties of the reservoir, they, in combination with reservoir simulation could potentially predict the CO{sub 2} saturation directly from the time-lapse seismic data. It is therefore concluded that the time-lapse seismic data could be used to monitor the carbon sequestrated saline reservoirs. But for its successful implementation, seismic modeling and inversion methods must be integrated with the reservoir simulations. In addition, because CO{sub 2} sequestration induces an effective anisotropy in the sequestered reservoir and anisotropy is best detected using multicomponent seismic data compared to single component (P-wave) data, acquisition, processing, and analysis is multicomponent seismic data is recommended for these time-lapse studies. Finally, a successful implementation of using time-lapse seismic data for monitoring the carbon sequestrated saline reservoirs will require development of a robust methodology for inverting multicomponent seismic data for subsurface anisotropic properties.
54 ENVIRONMENTAL SCIENCES↗