Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SLOPE”

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

Self-consistent core-pedestal ITER scenario modeling

The purpose of these integrated ITER simulations is to identify dependencies that can impact the performance of ITER. The dependence of fusion power production, temperature and density pedestals, and core profiles on varying magnetic-q, edge density fueling strength, neoclassical transport, magnetic field strength, and alpha heating, are examined. It is noted that ion modes usually dominate the interior, while electron modes dominate the barrier. The electron mode without flow shear is substantially stronger. The q increase tends to cause more kinetic ballooning modes at the edge and the q decrease means that the peeling is getting worse, but the ballooning is getting better. The large edge q tends to provide hollow density while small edge q gives normal density profile. This is consistent with the observation that the density of particles is inversely proportional to the density of the plasma current. The rise in the source of the edge particle increases the density of the edge and the reaction rate close to the edge increases temperature fluxes, resulting in weak pedestal barriers. When neoclassical transport is turned off or the strength of B-field is increased, the pronounced edge barriers are identified. In addition, a distinct effect on the magnetohydrodynamic (MHD) modes is observed due to the increased B-field. The interesting aspect is that the barrier collapses about q = 2 as the stability of MHD has decreased to replace the barrier with a continuous slope but with roughly the same fusion power. This would actually be good for a reactor. The slope of the H-mode pedestal is found to be reduced due to the alpha heating. It is expected that the problem of edge-localized modes causing damage to the first wall will be eased. The effects of circular, elongated, and general geometry on plasma profiles are also compared and contrasted. The circular geometry has a hollow density profile and no edge or internal thermal transport barrier (ETB or ITB). However, both ETBs and ITBs can be found in the elongated geometry. In general geometry, there is a weaker ETB, comparable ITBs, but a higher edge particle transport barrier than in elongated geometry. Higher density near the edge, on the other hand, causes more wall erosion, so elongated geometry might be the best option.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Range shifts in a foundation sedge potentially induce large Arctic ecosystem carbon losses and gains

Foundation species have disproportionately large impacts on ecosystem structure and function. As a result, future changes to their distribution may be important determinants of ecosystem carbon (C) cycling in a warmer world. We assessed the role of a foundation tussock sedge (Eriophorum vaginatum) as a climatically vulnerable C stock using field data, a machine learning ecological niche model, and an ensemble of terrestrial biosphere models (TBMs). Field data indicated that tussock density has decreased by ~0.97 tussocks per m2 over the past ~38 years on Alaska's North Slope from ~1981 to 2019. This declining trend is concerning because tussocks are a large Arctic C stock, which enhances soil organic layer C stocks by 6.9% on average and represents 745 Tg C across our study area. By 2100, we project that changes in tussock density may decrease the tussock C stock by 41% in regions where tussocks are currently abundant (e.g. -0.8 tussocks per m2 and -85 Tg C on the North Slope) and may increase the tussock C stock by 46% in regions where tussocks are currently scarce (e.g. +0.9 tussocks per m2 and +81 Tg C on Victoria Island). These climate-induced changes to the tussock C stock were comparable to, but sometimes opposite in sign, to vegetation C stock changes predicted by an ensemble of TBMs. Our results illustrate the important role of tussocks as a foundation species in determining future Arctic C stocks and highlight the need for better representation of this species in TBMs..

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

The ionomer as an oxygen evolution reaction promoter: piperidinium’s impact on mechanistic pathways on NiO, IrO 2 , and Fe–NiO

The commercial viability of anion exchange membrane (AEM) electrolysis requires optimization of various stack components, with specific catalyst-ionomer combinations often yielding higher current densities, lowered Tafel slopes, and improved mass activity. In this joint theoretical-experimental study, theoretical calculations detail the impact of Versogen’s piperidinium functional group on the complex, kinetically limiting oxygen evolution reaction, finding that the functional group can act as a promoter of specific steps (O*/O 2 * formation; H 2 O/O 2 desorption with reaction enthalpies ranging between 0.2–0.6 eV at higher coverages of O x H y intermediates) on NiO and NiFeO x catalysts. In particular, Fe sites on the NiFeO x catalyst facilitate concerted mechanisms of O*/O 2 * formation and H 2 O desorption with a low enthalpy of 0.5 eV; O 2 desorption alone requires only 0.3 eV. In contrast, Versogen-IrO 2 results in stronger Ir–O bonds, where the enthalpies for bond breaking (Ir–OH 2 and Ir–O 2 ) are considerably higher (1.4 eV and 1.6 eV, respectively). Rotating disk electrode studies utilized commercially available NiO and IrO 2 and synthesized 7.5 wt % Fe in NiFeO x catalysts in combination with Versogen, a common AEM ionomer, and Nafion, an alternative binder. Electrochemical testing validated the impact of these mechanistic changes on ionomer-catalyst combinations, finding that Versogen particularly activates NiO and NiFeO x compared to IrO 2 . Following a 13.5 h hold at 1.8 V, mass activities and Tafel slopes improved to 34 ± 13 A g −1 and 79 ± 2 mV dec −1 (NiO) and 82 ± 4.9 A g −1 and 72 ± 2 mV dec −1 (NiFeO x ). In contrast, Versogen-IrO 2 only reached 17 ± 2.9 A g −1 and 81 ± 3 mV dec −1 . Optimization of the ionomer-catalyst can yield significant increases in performance from initial activity and after an electrochemical conditioning procedure: this enhancement to the mass activity resulted in a 200.9 ± 106.1% improvement for Versogen-NiFeO x and 1284.2 ± 260.5% for Versogen-NiO. In contrast, Nafion-NiFeO x and -NiO offered moderate improvements of 39.1 ± 30.5% and 120.9 ± 59.1%, respectively.

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

Stellar splashback: the edge of the intracluster light

We examine the outskirts of galaxy clusters in the C-EAGLE simulations to quantify the ‘edges’ of the stellar and dark matter distribution. The radius of the steepest slope in the dark matter, commonly used as a proxy for the splashback radius, is located at $\sim \, r_{200 \rm m}$; the strength and location of this feature depends on the recent mass accretion rate, in good agreement with previous work. Interestingly, the stellar distribution (or intracluster light, ICL) also has a well-defined edge, which is directly related to the splashback radius of the halo. Thus, detecting the edge of the ICL can provide an independent measure of the physical boundary of the halo, and the recent mass accretion rate. We show that these caustics can also be seen in the projected density profiles, but care must be taken to account for the influence of substructures and other non-diffuse material, which can bias and/or weaken the signal of the steepest slope. This is particularly important for the stellar material, which has a higher fraction bound in subhaloes than the dark matter. Finally, we show that the ‘stellar splashback’ feature is located beyond current observational constraints on the ICL, but these large projected distances (>> 1 Mpc) and low surface brightnesses (μ>>32 mag arcsec -2 ) can be reached with upcoming observational facilities such as the Vera C. Rubin Observatory, the Nancy Grace Roman Space Telescope, and Euclid.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

CODEX weak lensing mass catalogue and implications on the mass–richness relation

The COnstrain Dark Energy with X-ray clusters (CODEX) sample contains the largest flux limited sample of X-ray clusters at 0.35 < z < 0.65. It was selected from ROSAT data in the 10 000 square degrees of overlap with BOSS, mapping a total number of 2770 high-z galaxy clusters. We present here the full results of the CFHT CODEX programme on cluster mass measurement, including a reanalysis of CFHTLS Wide data, with 25 individual lensing-constrained cluster masses. We employ lensfit shape measurement and perform a conservative colour–space selection and weighting of background galaxies. Using the combination of shape noise and an analytic covariance for intrinsic variations of cluster profiles at fixed mass due to large-scale structure, miscentring, and variations in concentration and ellipticity, we determine the likelihood of the observed shear signal as a function of true mass for each cluster. We combine 25 individual cluster mass likelihoods in a Bayesian hierarchical scheme with the inclusion of optical and X-ray selection functions to derive constraints on the slope α, normalization β, and scatter σln λ|μ of our richness–mass scaling relation model in log-space: ${\langle {\rm In}\,\, \lambda\!\!\mid\!\!\mu\rangle = \alpha\mu + \beta,} $ with μ = ln (M 200c /M piv ), and M piv = 10 14.81 M ⊙ . We find a slope $\alpha = 0.49^{+0.20}_{-0.15}$, normalization $\exp (\beta) = 84.0^{+9.2}_{-14.8}$, and $\sigma _{\ln \lambda | \mu } = 0.17^{+0.13}_{-0.09}$ using CFHT richness estimates. In comparison to other weak lensing richness–mass relations, we find the normalization of the richness statistically agreeing with the normalization of other scaling relations from a broad redshift range (0.0 < z < 0.65) and with different cluster selection (X-ray, Sunyaev–Zeldovich, and optical).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

X-ray properties of the X-CLASS-redMaPPer galaxy cluster sample: the luminosity–temperature relation

ABSTRACT This article presents the results of a spectroscopic analysis of the X-CLASS-redMaPPer (XC1-RM) galaxy cluster sample. X-CLASS is a serendipitous search for clusters in X-ray wavebands based on the XMM–Newton archive, whereas redMaPPer is an optical cluster catalogue derived from the Sloan Digital Sky Survey (SDSS). The present sample comprises 92 X-ray extended sources identified in optical images within 1 arcmin separation. The area covered by the cluster sample is ∼ 27 deg2. The clusters span a wide redshift range (0.05 < z < 0.6) and 88 clusters benefit from spectrosopically confirmed redshifts using data from SDSS Data Release 14. We present an automated pipeline to derive the X-ray properties of the clusters in three distinct apertures: R500 (at fixed mass overdensity), Rfit (at fixed signal-to-noise ratio) and ${R}_{300\, {\rm kpc}}$ (fixed physical radius). The sample extends over wide temperature and luminosity ranges: from 1–10 keV and from 6 × 1042 to 11 × 1044 erg s−1, respectively. We investigate the luminosity–temperature (L–T) relation of the XC1-RM sample and find a slope equal to 3.03 ± 0.26. It is steeper than predicted by self-similar assumptions, in agreement with independent studies. A simplified approach is developed to estimate the amount and impact of selection biases that might be affecting our recoveredL–Tparameters. The result of this simulation process suggests that the measuredL–Trelation is biased to a steeper slope and higher normalization.

Molham, Mona↗
DOE OSTI2020

Blazar variability power spectra from radio up to TeV photon energies: Mrk 421 and PKS 2155−304

ABSTRACT We present the results of the power spectral density (PSD) analysis for the blazars Mrk 421 and PKS 2155−304, using good-quality, densely sampled light curves at multiple frequencies, covering 17 decades of the electromagnetic spectrum, and variability time-scales from weeks up to a decade. The data were collected from publicly available archives of observatories at radio from Owens Valley Radio Observatory, optical and infrared (B, V, R, I, J, H, and Kbands), X-rays from the Swift and the Rossi X-ray Timing Explorer, high and very high energy (VHE) γ-rays from the Fermi and Very Energetic Radiation Imaging Telescope Array System as well as the High Energy Stereoscopic System. Our results are: (1) the power-law form of the variability power spectra at radio, infrared, and optical frequencies have slopes ∼1.8, indicative of random-walk-type noise processes; (2) the power-law form of the variability power spectra at higher frequencies, from X-rays to VHE γ-rays, however, have slopes ∼1.2, suggesting a flicker noise-type process; and (3) there is significantly more variability power at X-rays, high and VHE γ-rays on time-scales ≲ 100 d, as compared to lower energies. Our results do not easily fit into a simple model, in which a single compact emission zone is dominating the radiative output of the blazars across all the time-scales probed in our analysis. Instead, we argue that the frequency-dependent shape of the variability power spectra points out a more complex picture, with highly inhomogeneous outflow producing non-thermal emission over an extended, stratified volume.

Goyal, Arti↗
DOE OSTI2021

The rapid transition from star formation to AGN-dominated rest-frame ultraviolet light at z ≃ 4

ABSTRACT With the advent of deep optical-to-near-infrared extragalactic imaging on the degree scale, samples of high-redshift sources are being selected that contain both bright star-forming (SF) galaxies and faint active galactic nuclei (AGN). In this study, we investigate the transition between SF- and AGN-dominated systems at z ≃ 4 in the rest-frame ultraviolet (UV). We find a rapid transition to AGN-dominated sources brightward of MUV ≃ −23.2. The effect is observed in the rest-frame UV morphology and size–luminosity relation, where extended clumpy systems become point-source-dominated, and also in the available spectra for the sample. These results allow us to derive the rest-frame UV luminosity function (LF) for the SF- and AGN-dominated subsamples. We find the SF-dominated LF is best fit with a double power law, with a lensed Schechter function being unable to explain the existence of extremely luminous SF galaxies at MUV ≃ −23.5. If we identify AGN-dominated sources according to a point-source morphology criterion, we recover the relatively flat faint-end slope of the AGN LF determined in previous studies. If we instead separate the LF according to the current spectroscopic AGN fraction, we find a steeper faint-end slope of α = −1.83 ± 0.11. Using a simple model to predict the rest-frame AGN LF from the z = 4 galaxy LF, we find that the increasing impact of host galaxy light on the measured morphology of faint AGN can explain our observations.

Bowler, R. A. A.↗
DOE OSTI2021

SDSS-IV MaNGA: Stellar M/L gradients and the M/L-colour relation in galaxies

The stellar mass-to-light ratio gradient in SDSS r-band ∇(M * /L r ) of a galaxy depends on its mass assembly history, which is imprinted in its morphology and gradients of age, metallicity, and stellar initial mass function (IMF). Taking a MaNGA sample of 2051 galaxies with stellar masses ranging from 10 9 to 10 12 M ⊙ released in SDSS DR15, we focus on face-on galaxies, without merger and bar signatures, and investigate the dependence of the 2D ∇(M * /L r ) on other galaxy properties, including M * /L r -colour relationships by assuming a fixed Salpeter IMF as the mass normalization reference. The median gradient is ∇M * /L r ~ –0.1 (i.e. the M * /L r is larger at the centre) for massive galaxies, becomes flat around M * ~ 10 10 M ⊙ and change sign to ∇M * /L r ~ 0.1 at the lowest masses. The M * /L r inside a half-light radius increases with increasing galaxy stellar mass; in each mass bin, early-type galaxies have the highest value, while pure-disc late-type galaxies have the smallest. Correlation analyses suggest that the mass-weighted stellar age is the dominant parameter influencing the M * /L r profile, since a luminosity-weighted age is easily affected by star formation when the specific star formation rate (sSFR) inside the half-light radius is higher than 10 –3 Gyr –1 . With increased sSFR gradient, one can obtain a steeper negative ∇(M * /L r ). The scatter in the slopes of M * /L-colour relations increases with increasing sSFR, for example, the slope for post-starburst galaxies can be flattened to 0.45 from the global value 0.87 in the M * /L versus g – r diagram. Hence converting galaxy colours to M * /L should be done carefully, especially for those galaxies with young luminosity-weighted stellar ages, which can have quite different star formation histories.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Rotation curves and scaling relations of extremely massive spiral galaxies

We study the kinematics and scaling relations of a sample of 43 giant spiral galaxies that have stellar masses exceeding $10^{11} \, {\rm M}_\odot$ and optical discs up to 80 kpc in radius. We use a hybrid 3D–1D approach to fit 3D kinematic models to long-slit observations of the H α-$\rm{[N\, \small {II}]}$ emission lines and we obtain robust rotation curves of these massive systems. We find that all galaxies in our sample seem to reach a flat part of the rotation curve within the outermost optical radius. We use the derived kinematics to study the high-mass end of the two most important scaling relations for spiral galaxies: the stellar/baryonic mass Tully–Fisher relation and the Fall (mass-angular momentum) relation. All galaxies in our sample, with the possible exception of the two fastest rotators, lie comfortably on both these scaling relations determined at lower masses, without any evident break or bend at the high-mass regime. When we combine our high-mass sample with low-mass data from the Spitzer Photometry & Accurate Rotation Curves catalogue, we find a slope of α = 4.25 ± 0.19 for the stellar Tully–Fisher relation and a slope of γ = 0.64 ± 0.11 for the Fall relation. Our results indicate that most, if not all, of these rare, giant spiral galaxies are scaled up versions of less massive discs and that spiral galaxies are a self-similar population of objects up to the very high-mass end.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

OzDES Reverberation Mapping Programme: Mg ii lags and R − L relation

ABSTRACT The correlation between the broad line region radius and continuum luminosity (R–L relation) of active galactic nuclei (AGNs) is critical for single-epoch mass estimates of supermassive black holes (SMBHs). At z ∼ 1–2, where AGN activity peaks, the R–L relation is constrained by the reverberation mapping (RM) lags of the Mg ii line. We present 25 Mg ii lags from the Australian Dark Energy Survey RM project based on 6 yr of monitoring. We define quantitative criteria to select good lag measurements and verify their reliability with simulations based on both the damped random walk stochastic model and the rescaled, resampled versions of the observed light curves of local, well-measured AGN. Our sample significantly increases the number of Mg ii lags and extends the R–L relation to higher redshifts and luminosities. The relative iron line strength $\mathcal {R}_{\rm Fe}$ has little impact on the R–L relation. The best-fitting Mg iiR–L relation has a slope α = 0.39 ± 0.08 with an intrinsic scatter $\sigma _{\rm rl} = 0.15^{+0.03}_{-0.02}$ . The slope is consistent with previous measurements and shallower than the H β R–L relation. The intrinsic scatter of the new R–L relation is substantially smaller than previous studies and comparable to the intrinsic scatter of the H β R–L relation. Our new R–L relation will enable more precise single-epoch mass estimates and SMBH demographic studies at cosmic noon.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

MaNGA DynPop – I. Quality-assessed stellar dynamical modelling from integral-field spectroscopy of 10K nearby galaxies: a catalogue of masses, mass-to-light ratios, density profiles, and dark matter

ABSTRACT This is the first paper in our series on the combined analysis of the Dynamics and stellar Population (DynPop) for the MaNGA survey in the final SDSS Data Release 17 (DR17). Here, we present a catalogue of dynamically determined quantities for over 10 000 nearby galaxies based on integral-field stellar kinematics from the MaNGA survey. The dynamical properties are extracted using the axisymmetric Jeans Anisotropic Modelling (JAM) method, which was previously shown to be the most accurate for this kind of study. We assess systematic uncertainties using eight dynamical models with different assumptions. We use two orientations of the velocity ellipsoid: either cylindrically aligned JAMcyl or spherically aligned JAMsph. We also make four assumptions for the models’ dark versus luminous matter distributions: (1) mass-follows-light, (2) free NFW dark halo, (3) cosmologically constrained NFW halo, (4) generalized NFW dark halo, i.e. with free inner slope. In this catalogue, we provide the quantities related to the mass distributions (e.g. the density slopes and enclosed mass within a sphere of a given radius for total mass, stellar mass, and dark matter mass components). We also provide the complete models which can be used to compute the full luminous and mass distribution of each galaxy. Additionally, we visually assess the qualities of the models to help with model selections. We estimate the observed scatter in the measured quantities which decreases as expected with improvements in quality. For the best data quality, we find a remarkable consistency of measured quantities between different models, highlighting the robustness of the results.

Astronomy & Astrophysics↗
DOE OSTI2023

Cosmic variance and the inhomogeneous UV luminosity function of galaxies during reionization

When the first galaxies formed and starlight escaped into the intergalactic medium to reionize it, galaxy formation and reionization were both highly inhomogeneous in time and space, and fully coupled by mutual feedback. To show how this imprinted the UV luminosity function (UVLF) of reionization-era galaxies, we use our large-scale, radiation-hydrodynamics simulation CoDa II to derive the time- and space-varying halo mass function and UVLF, from z ≃ 6–15. That UVLF correlates strongly with local reionization redshift: earlier-reionizing regions have UVLFs that are higher, more extended to brighter magnitudes, and flatter at the faint end than later-reionizing regions observed at the same z. In general, as a region reionizes, the faint-end slope of its local UVLF flattens, and, by z = 6 (when reionization ended), the global UVLF, too, exhibits a flattened faint-end slope, ‘rolling-over’ at M UV ≳ -17. CoDa II’s UVLF is broadly consistent with cluster-lensed galaxy observations of the Hubble Frontier Fields at z = 6–8, including the faint end, except for the faintest data point at z = 6, based on one galaxy at M UV = -12.5. According to CoDa II, the probability of observing the latter is $\sim 5~{{\ \rm per\ cent}}$. However, the effective volume searched at this magnitude is very small, and is thus subject to significant cosmic variance. Here, we find that previous methods adopted to calculate the uncertainty due to cosmic variance underestimated it on such small scales by a factor of 2–4, primarily by underestimating the variance in halo abundance when the sample volume is small.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Rapidity-dependent spin decomposition of the nucleon

We revisit the two-dimensional Fourier transform of generalized parton distributions (GPDs) at nonzero skewness. At 𝜂 = 0 it reduces to the standard impact-parameter density, while at 𝜂 ≠ 0 it is an off-forward amplitude that we interpret as a genuine parton–nucleon correlation. Its overall strength (the transverse-plane integral of the density) is fixed by the GPD at the kinematic point 𝑡 =−𝑐 𝜂 =−4⁢𝜂 2 ⁢𝑚$^{2}_{𝑁}$/(1 − 𝜂 2 ) and decreases monotonically with the rapidity gap Δ⁢𝑦 = ln⁡[(1+𝜂)/(1−𝜂)] = 2 artanh⁡(𝜂). This rapidity dependence implies rapidity-modified Ji identities that connect helicity, orbital, and total angular momenta of the correlation in closed form. To quantify these effects, we construct leading-twist quark and gluon GPDs in a string-based conformal framework: conformal moments are parametrized by linear open- and closed-string Regge trajectories with slopes constrained by parton distribution functions (PDFs), hadron/glueball spectroscopy, and form-factor data, and GPDs are reconstructed over the full (𝑥,𝜂,𝑡) domain by Mellin-Barnes inversion with next-to-leading order evolution. We find qualitative agreement (and fair quantitative agreement within quoted uncertainties) for several moments and selected nonsinglet 𝑥-space channels at 𝜇 = 2 GeV when compared with lattice QCD, while we also identify channels with visible tension and discuss likely sources (PDF priors and 𝑡-slope systematics).

Gauge-gravity dualities↗
DOE OSTI2026

Importance of finite-size corrections for accurate ab initio modeling of carrier capture at semiconductor defects: A case study of substitutional C N in GaN

In ab initio studies of carrier-capture processes in defective semiconductor materials, the single-effective-mode formalism and the static-coupling approximation have become the predominant theoretical approaches for determining carrier-capture coefficients. The single-mode formalism relies on accurate nonequilibrium defect energies obtained from density-functional theory (DFT), where required inputs are a series of configurationally displaced, defect-containing supercells obtained using an interpolative ansatz, and where the DFT outputs are corresponding total energies that have traditionally been postprocessed using a long-established ground-state formulation of finite-size corrections and defect-formation energies. This formulation remains commonly used even though the defects that form a configuration-coordinate (CC) diagram typically exist as structures that are displaced from the ground state. To remedy this inconsistency, Kumagai has recently proposed novel methods for implementing finite-size corrections specifically intended for DFT calculations of the defect energies used to construct CC diagrams and implement the single-mode formalism [Y. Kumagai, Phys. Rev. B 107, L220101 (2023)]. Kumagai's approach builds on the latest finite-size-correction methods introduced to describe vertical charge-state transitions for charge-localizing point defects in semiconductors and insulators [T. Gake et al., Phys. Rev. B 101, 020102 (2020); S. Falletta et al., Phys. Rev. B 102, 041115 (2020)]. The newly identified finite-size artifact treated in these studies is the polarization charge induced on a configurationally frozen defect and its subsequent interaction with a vertical transition in charge state. In this work, we evaluate Kumagai's proposed methodology by applying it in a high-precision DFT study of carrier capture by substitutional C N in GaN, a well-characterized and technologically relevant defect and material. We have rigorously calculated C N defect energies across various supercell sizes for each defect configuration and charge state on the hole-capture CC diagram of C N (𝑞=−1), enabling a direct comparison of the slopes of the defect energies versus inverse cell size with those predicted by Kumagai. The most consequential prediction of Kumagai's method is that these slopes distinctly vary as the square of the linear-interpolation parameter used to construct the nonequilibrium defect configurations. Our results quantitatively support this prediction. Moreover, with these new finite-size corrections and multiple-cell-size DFT calculations in place, we find that the classical energy barrier for hole capture by C N (𝑞=−1) in GaN decreases to 0.092–0.127 eV. This finding confirms the recent ≈ 0.1 eV prediction of Reshchikov based on the weak temperature dependence for hole capture observed in photoluminescence experiments [M. A. Reshchikov, J. Appl. Phys. 129, 121101 (2021)]. These results stand in stark contrast to previously calculated barriers of 0.486 and 0.73 eV, which also used the single-mode formalism but were obtained by instead using ground-state-based finite-size corrections. Our reduced classical barrier for capture increases the temperature-dependent hole-capture coefficient of a C N (𝑞=−1) defect by more than two to four orders of magnitude for temperatures of 100–600 K, compared to the previous 0.486 eV results. While other defects may not be as dramatically affected as here, we suggest that incorporating proper finite-size corrections for the vertical-transition-like states embedded within CC diagrams is an essential, yet previously unrecognized, component of accurate modeling of carrier-capture when using the single-effective-mode formalism.

dielectric properties↗
DOE OSTI2021

Shock-compressed silicon: Hugoniot and sound speed up to 2100 GPa

In this work, the high-pressure equation of state and isentropic sound speed of fluid silicon were studied using laser-driven shock waves. Principal Hugoniot measurements to 2100 GPa were performed using impedance matching to an α-quartz reference. Sound speed was determined by time correlating imposed shock-velocity perturbations in both the sample (Si) and reference material (α-quartz). A change in shock velocity versus particle velocity (u s –u p ) slope on the fluid silicon principal Hugoniot was detected at 200 GPa. Density functional theory-based quantum molecular dynamics simulations suggest that an increase in ionic coordination and 50% increase in average ionization is coincident with the experimentally observed change in slope.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Magnetic excitations, nonclassicality, and quantum wake spin dynamics in the Hubbard chain

Recent work has demonstrated that quantum Fisher information (QFI), a witness of multipartite entanglement, and magnetic Van Hove correlations G(r,t), a probe of local real-space real-time spin dynamics, can be successfully extracted from inelastic neutron scattering on spin systems through accurate measurements of the dynamical spin structure factor S(k,ω). In this work, we apply theoretically these ideas to the half-filled Hubbard chain with nearest-neighbor hopping, away from the strong-coupling limit. This model has nontrivial redistribution of spectral weight in S(k,ω) going from the noninteracting limit (U = 0) to strong coupling (U→∞), where it reduces to the Heisenberg quantum spin chain. We use the density matrix renormalization group to find S(k,ω), from which QFI is then calculated. We find that QFI grows with U. With realistic energy resolution it becomes capable of witnessing bipartite entanglement above U=2.5 (in units of the hopping), where it also changes slope. This point is also proximate to slope changes of the bandwidth W(U) and the half-chain von Neumann entanglement entropy. We compute G(r,t) by Fourier transforming S(k,ω). The results indicate a crossover in the short-time short-distance dynamics at low U characterized by ferromagnetic light-cone wavefronts, to a Heisenberg-type behavior at large U featuring antiferromagnetic light cones and spatially period-doubled antiferromagnetism. We find this crossover has largely been completed by U=3. Our results thus provide evidence that, in several aspects, the strong-coupling limit of the Hubbard chain is reached qualitatively already at a relatively modest interaction strength. We discuss experimental candidates for observing the G(r,t) dynamics found at low U.

1-dimensional systems↗
DOE OSTI2023

Anomalous polarization reversal in strained thin films of CuIn⁢P 2 ⁢S 6

Strain-induced transitions of polarization reversal in thin films of a ferrielectric CuIn⁢P 2 ⁢S 6 (CIPS) with ideally conductive electrodes are explored using the Landau-Ginzburg-Devonshire approach with an eighth-order free energy expansion in polarization powers. Due to multiple potential wells, the height and position of which are temperature and strain dependent, the energy profiles of CIPS can flatten in the vicinity of the nonzero polarization states. Thereby, we reveal an unusually strong effect of a mismatch strain on the out-of-plane polarization reversal, hysteresis loop shape, dielectric susceptibility, and piezoelectric response of CIPS films. By varying the sign of the mismatch strain and its magnitude in a fairly narrow range, quasistatic hysteresisless paraelectric curves can transform into double, triple, and other types of pinched and single hysteresis loops. The strain effect on the polarization reversal is opposite, i.e., anomalous, in comparison with many other ferroelectric films, for which the out-of-plane remanent polarization and coercive field increase strongly for tensile strains and decrease or vanish for compressive strains. For definite values of temperature and mismatch strain, the low-frequency hysteresis loops of polarization may exhibit negative slope in a relatively narrow range of external field amplitude and frequency. The low-frequency susceptibility hysteresis loops, which correspond to the negative slope of polarization loops, contain only positive values, which can be giant in the entire range of field changes. In conclusion, the corresponding piezoresponse also reaches giant values, being maximal near coercive fields.

Dielectric properties↗