SEARCH · Engineering Papers

Results for “favorability”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Chapman–Jouguet deflagration criteria and compressibility dynamics of turbulent fast flames for turbulence-induced deflagration-to-detonation transition

This work characterizes the compressibility dynamics in turbulent fast flames for a range of turbulent flame speeds. These turbulent fast flames experience increased effects of compressibility through the formation of strong shocks and may develop a runaway acceleration combined with a pressure buildup that leads to turbulence induced deflagration-to-detonation transition (tDDT). Simultaneous high-speed particle image velocimetry, OH* chemiluminescence, schlieren, and pressure measurements are used to examine the reacting flow field and flame dynamics. We examine flames with turbulent flame speeds ranging from 100 to 600 m/s. At lower turbulent flame speeds, the flame is not able to produce favorable background conditions for deflagration-to-detonation transition (DDT) onset, and thus flame compressibility and turbulence amplification are less dominant, resulting in a weaker acoustic coupling between the flame and compressed region. As the turbulent burning velocities exceed the Chapman–Jouguet deflagration speed, favorable background conditions are produced, as we observe flame-generated shocks and flame-generated turbulence with higher turbulent velocities and larger turbulent scales. At this regime, the flame is categorized to be at the runaway transition regime that leads to tDDT.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Insights into the structure and dynamics of K + ions at the muscovite–water interface from machine learning potential simulations

The surfaces of many minerals are covered by naturally occurring cations that become partially hydrated and can be replaced by hydronium or other cations when the surface is exposed to water or an aqueous solution. These ion exchange processes are relevant to various chemical and transport phenomena, yet elucidating their microscopic details is challenging for both experiments and simulations. Here, in this work, we make a first step in this direction by investigating the behavior of the native K + ions at the interface between neat water and the muscovite mica (001) surface with ab-initio -based machine learning molecular dynamics and enhanced sampling simulations. Our results show that the desorption of the surface K + ions in pure ion-free water has a significant free energy barrier irrespective of their local surface arrangement. In contrast, facile K + diffusion between mica’s ditrigonal cavities characterized by different Al/Si orderings is observed. This behavior suggests that the K + ions may favor a dynamic disordered surface arrangement rather than complete desorption when exposed to deionized water.

Ab-initio molecular dynamics↗

DOE OSTI2023

Structure-guided approach to modulate small molecule binding to a promiscuous ligand-activated protein

Ligand-binding promiscuity in detoxification systems protects the body from toxicological harm but is a roadblock to drug development due to the difficulty in optimizing small molecules to both retain target potency and avoid metabolic events. Immense effort is invested in evaluating metabolism of molecules to develop safer, more effective treatments, but engineering specificity into or out of promiscuous proteins and their ligands is a challenging task. To better understand the promiscuous nature of detoxification networks, we have used X-ray crystallography to characterize a structural feature of pregnane X receptor (PXR), a nuclear receptor that is activated by diverse molecules (with different structures and sizes) to up-regulate transcription of drug metabolism genes. We found that large ligands expand PXR’s ligand-binding pocket, and the ligand-induced expansion occurs through a specific unfavorable compound-protein clash that likely contributes to reduced binding affinity. Removing the clash by compound modification resulted in more favorable binding modes with significantly enhanced binding affinity. We then engineered the unfavorable ligand-protein clash into a potent, small PXR ligand, resulting in marked reduction in PXR binding and activation. Structural analysis showed that PXR is remodeled, and the modified ligands reposition in the binding pocket to avoid clashes, but the conformational changes result in less favorable binding modes. Thus, ligand-induced binding pocket expansion increases ligand-binding potential of PXR but is an unfavorable event; therefore, drug candidates can be engineered to expand PXR’s ligand-binding pocket and reduce their safety liability due to PXR binding.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Precursor-mediated in situ growth of hierarchical N-doped graphene nanofibers confining nickel single atoms for CO 2 electroreduction

Despite the various strategies for achieving metal–nitrogen–carbon (M–N–C) single-atom catalysts (SACs) with different microenvironments for electrochemical carbon dioxide reduction reaction (CO 2 RR), the synthesis–structure–performance correlation remains elusive due to the lack of well-controlled synthetic approaches. Here, we employed Ni nanoparticles as starting materials for the direct synthesis of nickel (Ni) SACs in one spot through harvesting the interaction between metallic Ni and N atoms in the precursor during the chemical vapor deposition growth of hierarchical N-doped graphene fibers. By combining with first-principle calculations, we found that the Ni-N configuration is closely correlated to the N contents in the precursor, in which the acetonitrile with a high N/C ratio favors the formation of Ni-N 3 , while the pyridine with a low N/C ratio is more likely to promote the evolution of Ni-N 2 . Moreover, we revealed that the presence of N favors the formation of H-terminated edge of sp 2 carbon and consequently leads to the formation of graphene fibers consisting of vertically stacked graphene flakes, instead of the traditional growth of carbon nanotubes on Ni nanoparticles. With a high capability in balancing the *COOH formation and *CO desorption, the as-prepared hierarchical N-doped graphene nanofibers with Ni-N 3 sites exhibit a superior CO2RR performance compared to that with Ni-N 2 and Ni-N 4 ones.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Increased frequency of planetary wave resonance events over the past half-century

We demonstrate a tripling in the frequency of planetary wave resonance events over the past halfcentury, coinciding with the rise in persistent boreal summer weather extremes. This increase aligns with changes in the underlying climate conditions favoring these events, including amplified Arctic warming and land-sea thermal contrast. We also observe increased prevalence of resonant amplification events following the mature phase of strong El Niño events, suggesting that such events may precondition the mean state conditions in ways that favor large-scale quasi-stationary wave patterns and quasi-resonant wave amplification. Since the impact of anthropogenic warming on quasi-resonant amplification is not well captured by current-generation climate models, it is likely that models are underpredicting the potential increase, indicating even greater risk of persistent extreme summer weather events with ongoing warming.

Arctic amplification↗

DOE OSTI2021

Incorporation of Ti in epitaxial Fe2TiO4 thin films

The titanomagnetites (Fe2-xTixO4, x = 1) are a family of reducible spinel-structure oxides of interest for their favorable magnetic, catalytic, and electrical transport properties. To understand the stability of the system during low temperature deposition, epitaxial thin films of Fe2TiO4 were deposited by molecular beam epitaxy on MgO(001) at 250 – 375°C. The homogeneous incorporation of Ti, Fe valence, and film morphology were all found to be strongly dependent on the oxidation conditions at the low substrate temperatures employed. More oxidizing conditions led to phase separation into epitaxial, faceted Fe3O4 and rutile TiO2. Less oxidizing conditions resulted in polycrystalline films that exhibited Ti segregation to the film surface, as well as mixed Fe valence (Fe3+, Fe2+, Fe0). A narrow window of intermediate oxygen partial pressure during deposition yielded nearly homogeneous Ti incorporation and a large fraction of Fe2+. However, these films were poorly crystallized, and no occupation of tetrahedral sites in the spinel lattice by Fe2+ was detected by x-ray magnetic circular dichroism at the Fe L-edge. After vacuum annealing, a small fraction of Fe2+ was found to occupy tetrahedral sites. Comparison of these results with previous work suggests that the restriction of the deposition conditions to low temperature by the MgO substrate limits the incorporation of Ti into the spinel lattice. This work suggests a path towards obtaining stoichiometric, well-crystallized Fe2TiO4 by MBE by utilizing high substrate temperature and low oxygen partial pressure during deposition on temperature-stable substrates.

Kaspar, Tiffany C.↗

DOE OSTI2021

Numerical investigation of active control of tearing mode by magnetic coils and the role of Δ'

Magnetic feedback stabilization of the tearing mode (TM) is numerically investigated, utilizing the MARS-F code (Liu et al 2000 Phys. Plasmas 7 3681) for toroidal tokamak equilibria. With control coil configurations assumed in this study, magnetic feedback partially or fully stabilizes the TM, with either vanishing or finite equilibrium pressure. The best control is achieved by the combination of internal active coils and internal poloidal sensors. The internal and external tearing indices are evaluated for the close-loop system, based on the MARS-F computed mode eigenvalue and eigenfunction, respectively. In the absence of the favorable curvature effect, these two indices are real-valued and quantitatively agree well with each other. For the equilibrium with finite pressure gradient at the mode rational surface, the favorable average curvature effect becomes important and the close-loop tearing index also becomes complex-valued, partly due to interaction of the feedback system with the dissipative wall eddy current response. Isolating the inner layer and outer region response to magnetic feedback, with either proportional or proportional-derivative actions, allows to establish that feedback stabilization of the TM occurs mainly due to modification of the behavior of the external ideal solution, further confirming the analytic result reported in He et al 2021 Phys. Plasmas 28 012504.

Physics↗

DOE OSTI2022

Experimental study on the role of the target electron temperature as a key parameter linking recycling to plasma performance in JET-ILW*

Changes in global and edge plasma parameters (H98(y,2), dimensionless collisionality ν*, core density peaking, separatrix density ne,sep) with variations in the D2 fueling rate and divertor configuration are unified into a single trend when mapped to $\langle$Te,ot$\rangle$, the spatially averaged spectroscopically derived outer target electron temperature. Dedicated JET with the ITER-like wall (JET-ILW) experiments in combination with an extended JET-ILW database of unseeded low-triangularity H-mode plasmas spanning a wide range of D2 fueling rates, Ip, Bt and heating power have demonstrated the importance of $\langle$Te,ot$\rangle$ as a key physics parameter linking the recycling particle source and detachment with plasma performance. The remarkably robust H98(y,2) trend with $\langle$Te,ot$\rangle$ is connected to a strong inverse correlation between $\langle$Te,ot$\rangle$, ne,sep and ν*, thus directly linking changes in the divertor recycling moderated by $\langle$Te,ot$\rangle$ with the previously established relationship between ν*, core density peaking and core pressure resulting in a degradation in core plasma performance with decreasing $\langle$Te,ot$\rangle$ (increasing ν*). Additionally, a strong inverse correlation between the separatrix to pedestal density ratio, ne,sep/ne,ped, and $\langle$Te,ot$\rangle$ is also established, with the rise in ne,sep/ne,ped saturating at $\langle$Te,ot$\rangle$ > 10 eV. A strong reduction in H98(y,2) is observed as $\langle$Te,ot$\rangle$ is driven from 30 to 10 eV via additional D2 gas fueling, while the divertor remains attached. Consequently, the pronounced performance degradation in attached divertor conditions has implications for impurity seeding radiative divertor scenarios, in which H98(y,2) is already low (~0.7) before impurities are injected into the plasma since moderate gas fueling rates are required to promote high divertor neutral pressure. A favorable pedestal pressure, pe,ped, dependence on Ip has also been observed, with an overall increase in pe,ped at Ip = 3.4 MA as $\langle$Te,ot$\rangle$ is driven down from attached to high-recycling divertor conditions. In contrast, pe,ped is reduced with decreasing $\langle$Te,ot$\rangle$ in the lower Ip branches. Further work is needed to (i) clarify the potential role of edge opacity on the observed favorable pedestal pressure Ip scaling; as well as to (ii) project the global and edge plasma performance trends with $\langle$Te,ot$\rangle$ to reactor-scale devices to improve predictive capability of the coupling between recycling and confined plasma fueling in what are foreseen to be more opaque edge plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Isotope effects under the influence of global radial electric fields in a helical configuration

Isotope effects under the influence of a radial electric field are examined in a helical magnetic field configuration. We perform global gyrokinetic simulations with additional poloidal rotations to estimate quasi-linear heat flux due to ion temperature gradient mode under the mixing length model. In single-ion-species plasmas, the mass number dependency of heat flux agrees with gyro-Bohm scaling in the absence of a radial electric field. Favorable mass number dependencies violating gyro-Bohm scaling are observed in the presence of a global radial electric field or a heavy hydrogen component in multi-ion-species plasmas. The radial electric field and the heavy hydrogen component affect the heat flux through an increase of wavelength as well as mode stabilization. Poloidal Mach number characterizes the transition from unfavorable to favorable mass number dependency under radial electric fields. While the heat flux is independent of mass number for a given poloidal Mach number, the heat flux decreases for higher mass numbers in a given radial electric field. The heat flux is also independent of average mass number in multi-ion-species plasmas because the heavy hydrogen component effectively enhances the light hydrogen heat flux. The present results are potentially relevant to the violation of gyro-Bohm scaling observed in the recent deuterium experiments in the Large Helical Device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Feasibility of main thermal ion heating by ICRF waves using a top launcher in a tokamak with deuterium–tritium plasmas

A scenario of ion cyclotron range of frequency (ICRF) wave injection from a top launcher is proposed as an efficient and direct heating method for thermal deuterium ions in deuterium–tritium tokamak plasmas. Positioned between the tritium cyclotron layer and ion–ion hybrid layer, the top launcher allows effective wave penetration to the ion–ion hybrid layer and enables significant power transfer to thermal deuterium. This is achieved through favorable wave polarization for fundamental cyclotron damping. There is a Doppler broadening around the cyclotron resonance and this overlaps with the ion–ion hybrid layer. Low toroidal mode numbers and ion temperature in the range of 5–20 keV are favorable for enhancing the main ion damping relative to electron damping. In contrast to the neutral beam injection, which penetration strongly depends on machine size and plasma density, the proposed ICRF-based direct ion heating scenario is shown to be scalable and applicable to both larger and smaller tokamak devices within practical constraints.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Effect of outer divertor leg detachment on the high field side scrape-off layer in DIII-D and ASDEX Upgrade

In this work, evidence is presented that detachment of the outer divertor leg leads to a reduction of electron density and neutral pressure in the high-field side (HFS) scrape-off layer (SOL) of the ASDEX Upgrade (AUG) and DIII-D tokamaks with ion $B$ x $∇B$ drift directed toward the X-point (favorable configuration). These results are observed across multiple diagnostics and without the use of impurity seeding to reach detachment. In AUG, outer divertor leg detachment correlates with a decrease in electron density near the separatrix at the inner midplane, measured with HFS reflectometry. A concurrent reduction in inner divertor density and neutral pressure at the inner target is observed using divertor Thomson scattering and neutral pressure gauges. These effects are present in both L- and H-mode plasmas. In DIII-D, a similar reduction is detected through analysis of line-integrated hydrogenic emission measured by multiple diagnostics in the HFS SOL close to the separatrix. The consistent trends in both devices indicate that high electron density and strong hydrogenic emissivity in the HFS SOL are common features of H-mode plasmas in the favorable configuration, independent of wall material and divertor geometry. In L-mode plasmas, the reduction in electron density and neutral pressure is not observed in DIII-D, possibly due to differences in wall material. These results emphasize the importance of the divertor state in determining the two-dimensional neutral distribution and edge fueling.

detachment↗

DOE OSTI2020

DFT Plus DMFT Study of Spin-Charge-Lattice Coupling in Covalnt LaCoO{siub 3}.

We study energetics and the nature of both homogeneous and mixed spin (MS) states in LaCoO3 incorporating structural changes of the crystal volume expansion and the Co-O bond disproportionation (BD) during the spin-state transition using the density functional theory plus dynamical mean field theory (DFT + DMFT) method. DFT + DMFT predicts that energetics of both excited spin states are almost the same while DFT + U calculations of the same structures energetically favor the MS states and produce various metastable solutions whose energetics depend sensitively on final spin states. Within DFT + DMFT, the homogeneous spin state in the expanded crystal volume shows the multiconfigurational nature with non-negligible occupancy probabilities of both high spin (HS) and low spin (LS) states along with d(6) and d(7) charge configurations indicating the dynamically fluctuating nature of spin and charge states due to the Co-O covalency. The nature of the MS state under the BD structure reveals that Co sites with the long Co-O bonds develop a Mott insulating state and favor HS with a d(6) configuration, while more covalent Co sites with the short Co-O bonds occupy more LS states with a d(7) configuration and behave as a band insulator, as a result, charge ordering is induced in the BD structure from the spin-state ordering. We also find that both energetics and electronic structure sensitively depend on the Co-O covalency effect, which can be tuned by changing the double counting potential and the resulting d occupancy (N-d), and N-d close to 6.7 is consistent with the nature of the spin-state transition. Our results show that structural changes during the spin-state transition can play an important role in understanding energetics and electronic structure of LaCoO3.

Park, Hyowon↗

DOE OSTI2022

Pressure-induced creation and annihilation of Weyl points in T d - Mo 0.5 W 0.5 Te 2 and 1 T " - Mo 0.5 W 0.5 Te 2

By means of first-principles density-functional theory calculations, we investigate the role of hydrostatic pressure in the electronic structure of the T d (Pmn2 1 ) and 1T" (Pm) phases of the Weyl semimetal Mo 0.5 W 0.5 Te 2 , which is a promising material for phase-change memory technology and superconductivity. We particularly focus on changes occurring in the distribution of the gapless Weyl points (WPs) within the 0 to 45 GPa pressure range. Here we further investigate the structural phase transition and lattice dynamics of the T d and 1T" phases within the aforementioned pressure range. Our calculations suggest that both the T d and 1T" phases of Mo 0.5 W 0.5 Te 2 host four WPs in their full Brillouin zone at zero pressure. The total number of WPs increases to 44 (36) with increasing pressure via pair creation up to 20 (15) GPa for the T d (1T") phase, and beyond this pressure pair annihilation of WPs starts occurring, leaving only 16 WPs at 45 GPa in both phases. The enthalpy versus pressure data reveal that the 1T" phase is more favorable below the critical pressure of 7.5 GPa; however, beyond this critical pressure the T d phase becomes enthalpically favorable. We also provide the calculated x-ray diffraction spectra along with the calculated Raman- and infrared-active phonon frequencies to facilitate the experimental identification of the studied phases

36 MATERIALS SCIENCE↗

DOE OSTI2023

Measurement of the branching fractions for Cabibbo-suppressed decays D + → K + K - π + π 0 and D ( s ) + → K + π - π + π 0 at Belle

We present measurements of the branching fractions for the singly Cabibbo-suppressed decays and , and the doubly Cabibbo-suppressed decay , based on of data recorded by the Belle experiment at the KEKB e + ⁢e - collider. We measure these modes relative to the Cabibbo-favored modes D + → K - ⁢π + ⁢π + ⁢π 0 and D$^{+}_{s}$ → K + ⁢K -⁢ π + ⁢π 0 . Our results for the ratios of branching fractions are $\mathscr{B}$⁡(D + →K + ⁢K -⁢ π + ⁢π 0 )/$\mathscr{B}$⁡(D + → K - ⁢π + ⁢π + ⁢π 0 ) =(11.32±0.13±0.26)%, $\mathscr{B}$⁡(D + →K + π - ⁢π + ⁢π 0 )/$\mathscr{B}$(D + →K -⁢ π + ⁢π + ⁢π 0 ) =(1.68±0.11±0.03)%, and $\mathscr{B}$⁡(D$^{+}_{s}$ → K + ⁢π - π + ⁢π 0 )/$\mathscr{B}$⁡(D$^{+}_{s}$ → K + ⁢K - π + π 0 ) =(17.13±0.62±0.51)%, where the uncertainties are statistical and systematic, respectively. The second value corresponds to (5.83 ±0.42) ×tan 4 ⁡θ C , where θ C is the Cabibbo angle; this value is larger than other measured ratios of branching fractions for a doubly Cabibbo-suppressed charm decay to a Cabibbo-favored decay. Multiplying these results by world average values for $\mathscr{B}$⁡(D + → K - ⁢π + ⁢π + ⁢π 0 ) and $\mathscr{B}$⁡(D$^{+}_{s}$ → K + ⁢K - π + π 0 ) yields $\mathscr{B}$( + →K + ⁢K -⁢ π + ⁢π 0 ) = (7.08±0.08±0.16±0.20)×10 -3 , $\mathscr{B}$⁡(D + → K +⁢ π - π + ⁢π 0 ) =(1.05±0.07±0.02±0.03)×10 -3 , and $\mathscr{B}$⁡(D$^{+}_{s}$ → K + ⁢π - π + ⁢π 0 ) =(9.44±0.34±0.28±0.32)×10 -3 , where the third uncertainty is due to the branching fraction of the normalization mode. The first two results are consistent with, but more precise than, the current world averages. The last result is the first measurement of this branching fraction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Self-Interacting Dark Matter and the Origin of Ultradiffuse Galaxies NGC1052-DF2 and -DF4

Observations of ultradiffuse galaxies NGC 1052-DF2 and -DF4 show they may contain little dark matter, challenging our understanding of galaxy formation. Using controlled N-body simulations, we explore the possibility that their properties can be reproduced through tidal stripping from the elliptical galaxy NGC 1052, in both cold dark matter (CDM) and self-interacting dark matter (SIDM) scenarios. To explain the dark matter deficiency, we find that a CDM halo must have a very low concentration so that it can lose sufficient inner mass in the tidal field. In contrast, SIDM favors a higher and more reasonable concentration as core formation enhances tidal mass loss. Final stellar distributions in our SIDM benchmarks are more diffuse than the CDM one, and hence the former are in better agreement with the data. We further show that a cored CDM halo model modified by strong baryonic feedback is unlikely to reproduce the observations. Our results indicate that SIDM is more favorable for the formation of dark-matter-deficient galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Site-Specific Structure at Multiple Length Scales in Kagome Quantum Spin Liquid Candidates

Realizing a quantum spin liquid (QSL) ground state in a real material is a leading issue in condensed matter physics research. In this pursuit, it is crucial to fully characterize the structure and influence of defects, as these can significantly affect the fragile QSL physics. In this work, we perform a variety of cutting-edge synchrotron x-ray scattering and spectroscopy techniques, and we advance new methodologies for site-specific diffraction and L-edge Zn absorption spectroscopy. The experimental results along with our first-principles calculations address outstanding questions about the local and long-range structures of the two leading kagome QSL candidates, Zn-substituted barlowite Cu 3 Zn x Cu 1-x (OH) 6 FBr and herbertsmithite Cu 3 Zn(OH) 6 Cl 2 . On all length scales probed, there is no evidence that Zn substitutes onto the kagome layers, thereby preserving the QSL physics of the kagome lattice. Our calculations show that antisite disorder is not energetically favorable and is even less favorable in Zn-barlowite compared to herbertsmithite. Site-specific x-ray diffraction measurements of Zn-barlowite reveal that Cu 2+ and Zn 2+ selectively occupy distinct interlayer sites, in contrast to herbertsmithite. Using the first measured Zn L-edge inelastic x-ray absorption spectra combined with calculations, we discover a systematic correlation between the loss of inversion symmetry from pseudo-octahedral (herbertsmithite) to trigonal prismatic coordination (Zn-barlowite) with the emergence of a new peak. Overall, our measurements suggest that Zn-barlowite has structural advantages over herbertsmithite that make its magnetic properties closer to an ideal QSL candidate: its kagome layers are highly resistant to nonmagnetic defects while the interlayers can accommodate a higher amount of Zn substitution.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Dy adsorption on and intercalation under graphene on 6 H -SiC(0001) surface from first-principles calculations

Previous experimental observations motivate clarification of configuration stabilities and kinetic processes for intercalation of guest atoms into a layered van der Waals material such as a graphene-SiC system. From our first-principles density functional theory (DFT) calculations, we analyze Dy adsorption and intercalation for graphene on a 6H-SiC(0001) surface, where the system includes two single-atom-thick graphene layers: the top-layer graphene (TLG) and the underling buffer-layer graphene (BLG) above the terminal Si layer. Our chemical potential analysis shows that intercalation of a single Dy atom into the gallery between TLG and BLG is more favorable than adsorption on TLG but that intercalation into the gallery underneath BLG is highly unfavorable. We obtain diffusion barriers of ~0.45 and 0.54 eV for a Dy atom diffusing on and under TLG, respectively. We find that the direct penetration of a Dy atom from the graphene top into the gallery under TLG is almost inhibited below a temperature of ~1400 K due to a large global barrier of at least ~3.5 eV. Instead, we find that a single Dy atom on TLG can easily intercalate by crossing a TLG step (e.g., a zigzag step presaturated by a Dy chain or a reconstructed zigzag step zz57). Additionally, we also perform DFT calculations for different Dy coverages to demonstrate how the favorability of Dy intercalation, as well as the corresponding interlayer spacings, depend on the coverage. Consequently, we can provide general insight and guidance for extensively studied systems involving intercalation of foreign atoms into graphene on a SiC substrate.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Thermodynamics of boron incorporation in BGaN

We study the thermodynamics of boron (B) incorporation into gallium nitride (GaN) using first-principles calculations. In the dilute limit, we have calculated the formation energies of different configurations of the B impurity in GaN and found that substitution on the cation site is favored over substitution on the anion site. Furthermore, under p-type conditions, interstitial boron can become the more favorable configuration and will ultimately limit the p-type conductivity. At higher B concentrations we use the generalized quasi-chemical approximation to elucidate the thermodynamic stability of boron gallium nitride (BGaN) alloys. We also investigate the effects of strain, which will be present if BGaN alloys are grown pseudomorphically on a GaN substrate.Without strain, B incorporation at typical growth conditions is limited to about 1.4% at 800 °C. Pseudomorphic strain raises the limit to 3.0% at the same temperature, close to experimentally observed levels of B incorporation.

36 MATERIALS SCIENCE↗