Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Splitting”

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 1,261 records · Page 70

Importance of dynamic lattice effects for crystal field excitations in the quantum spin ice candidate Pr 2 Zr 2 O 7

We explore dynamic interactions between the crystal lattice and magnetic degrees of freedom in a frustrated magnetic system using the example of a pyrochlore quantum spin-ice candidate Pr 2 Zr 2 O 7 . Using Raman scattering spectroscopy we demonstrate that crystal electric field excitations of Pr3+, which define the magnetic properties of Pr 2 Zr 2 O 7 , cannot be understood within a model of a static lattice. We identify vibronic interactions with a phonon which lead to a splitting of a doublet crystal field excitation at around 55 meV. We also observe an unconventional behavior of a splitting of the non-Kramers ground state doublet of Pr 3+ , revealed by observing excitations to the first excited singlet state $E$$^{0}_{g}$ → A 1 g at around 10 meV. The splitting has a strong temperature dependence, where the doublet structure is most prominent between 50 and 100 K, and the weight of one of the components strongly decreases on cooling contrary to simple thermal population tendency. Furthermore, we suggest a static or dynamic deviation of Pr 3+ from the position in the ideal crystal structure can be the origin of the effect, with the deviation strongly decreasing at low temperatures.

36 MATERIALS SCIENCE↗

Trigonal symmetry breaking and its electronic effects in the two-dimensional dihalides MX 2 and trihalides MX 3

We study the consequences of the approximately trigonal (D 3d ) point symmetry of the transition metal (M) site in two-dimensional van der Waals MX 2 dihalides and MX 3 trihalides. The trigonal symmetry leads to a 2-2-1 orbital splitting of the transition metal d shell, which is best represented by d orbitals that describe the trigonal rather than Oh symmetry. The ligand-ligand bond length differences (rather than metal-ligand) take the role of a Jahn-Teller-like mode and in combination with interlayer distances and dimensionality effects tune the crystal field splittings and bandwidths. These effects, in turn, are amplified by electronic correlation effects—leading to crystal field splittings of the order of 0.1–1 eV between the singlet and the lowest orbital doublet—as opposed to the often assumed degenerate 3 t 2g states. Our calculations explain why most of the materials in this family are insulating, and why these considerations have to be taken into account in realistic models of them. Further, orbital order coupled to various lower symmetry lattice modes may lift the remaining orbital degeneracies, and we explain how these may support unique electronic states using ZrI 2 and CuCl 2 as examples, and offer a brief overview of possible electronic configurations in this class of materials. By building and analyzing Wannier models adapted to the appropriate symmetry we examine how the interplay among trigonal symmetry, electronic correlation effects, and – d orbital charge transfer leads to insulating, orbitally polarized magnetic and/or orbital-selective Mott states, and we provide a simple framework and numerical tool for others. Our paper establishes a rigorous framework to understand, control, and tune the electronic states in low-dimensional correlated halides. Furthermore, our analysis shows that trigonal symmetry and its breaking is a key feature of the two-dimensional halides that needs to be accounted for in search of novel electronic states in materials ranging from CrI 3 to α–RuCl 3 .

2-dimensional systems↗

Migration of Ga vacancies and interstitials in β – Ga 2 O 3

Pathways and energy barriers for the migration of Ga vacancies (V Ga ) and Ga interstitials (Ga i ) in β–Ga 2 O 3 are explored using hybrid functional calculations and the nudged elastic band method. Considering β–Ga 2 O 3 as primarily being an n-type semiconductor, we focus on defect charge states relevant under such conditions: $V^{3 –}_{Ga}$, Ga$^{3+}_{i}$, and Ga$^{+}_{i}$. Notably, we describe a mechanism by which V Ga can transform between its different split configurations. In all cases, the intermediate state consists of a vacancy split between three Ga sites—a three-split vacancy—which enables passage over a significantly lower energy barrier. This is because it avoids the unfavorable simple vacancy at the tetrahedral Ga site. The proposed mechanism lowers the overall barrier for $V^{3–}_{Ga}$ diffusion along the [001] crystal direction from 1.73 to 0.97 eV, whereas the 2.08 eV barrier for the [100] and [010] directions is unaffected. For Ga$^{3+}_{i}$, we obtain similar overall migration barriers of 0.72, 0.80, and 1.02 eV for the [010], [001], and [100] directions, respectively. Ga$^{+}_{i}$ exhibits a strong preference for diffusion within the large eight-sided channel; the overall migration barrier is 0.92 eV for the [010] direction, and 2.16 eV for the [001] and [100] directions. The limiting step for the two latter directions involves ionization of Ga$^{+}_{i}$ followed by a jump to an adjacent large eight-sided channel as Ga$^{3+}_{i}$. Finally, our results are discussed in light of experimental observations of thermally activated recovery processes in irradiated material.

36 MATERIALS SCIENCE↗

Local structure, bonding, and asymmetry of $\mathrm{((NH_2)_2CH)PbBr_3}$, $\mathrm{CsPbBr_3}$, and $\mathrm{(CH_3MH_3)PbBr_3}$

Here we report local structure measurements for CsPbBr 3 and ((NH 2 ) 2 CH)PbBr 3 (FAPbBr 3 ) and compare them with recent results for (CH 3 NH 3 )PbBr 3 (MAPbBr 3 ). The Pb-Br bonding is similar for all three systems; the effective spring constants, κ, are comparable (ranging from 1.20 to 1.95 eV/Å 2 ), but small in magnitude indicating very soft materials. However, there are also important differences between the three systems. Static disorder is very small for CsPbBr3 but increases somewhat with the size of the organic molecular ions MA + and FA + . At room temperature, dynamic disorder dominates in all compounds. The thermal disorder of the Pb-Br pair distribution function (PDF), i.e., the Debye-Waller factor σ 2 follows a correlated Debye or Einstein model up to 300 K in CsPbBr 3 (orthorhombic phase), but for FAPbBr 3 and MAPbBr 3 , there is a break in the σ 2 (T) curve at the orthorhombic-tetragonal transition (o-t) near 150 K, indicating a small change in the spring constant κ. κ increases for MAPbBr 3 but decreases in FAPbBr 3 at this transition. These changes are attributed to changes in the H-bonding between Br - and MA + or FA + at this transition, as a result of librations or rotations of these molecular cations. In addition, the Pb-Br PDF becomes asymmetric at a relatively low temperature for FAPbBr 3 and MAPbBr 3 , while this effect is significantly smaller for CsPbBr 3 . Finally, we address the question of a model to explain the asymmetric PDF. Two main models are discussed in the literature, an anharmonic pair potential and a split-pair distribution, possibly driven by the presence of a lone pair on the Pb ion. We show that the fourth cumulant C 4 can differentiate between these two models and other possible models. Experimentally C 4 is positive at 250 K and above, for all three systems and that is inconsistent with a split-peak model, for which C 4 is negative for splittings larger than 0.12 Å.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Comparison of interband related optical transitions and excitons in ZnGeN 2 and GaN

The optical dielectric function of ZnGeN 2 is calculated from the interband transitions using the energy bands calculated in the quasiparticle self-consistent (QS)GW^ method using two different levels of approximation: the independent particle approximation (IPA) and the Bethe-Salpeter Equation (BSE) approach. The first allows us to relate peaks in ε 2 to specific bands and k-points but does not include electron-hole interaction effects. The second includes electron-hole interaction or excitonic effects. The corresponding changes in the shape of ε 2 (ω) are found to be similar to those in GaN. The screened Coulomb interaction W^ is here calculated including electron-hole interactions in the polarization function and gives a band structure already going beyond the random phase approximation. The static dielectric constants including only electronic screening, commonly referred to as ε∞ , were calculated separately by extrapolating the wave vector dependent macroscopic dielectric function in the long-wavelength limit. Below the quasiparticle gap, we find three bound excitons optically active for different polarization. The convergence of these bound excitons with respect to the density of the k-mesh used in the BSE is studied and found to require a fine mesh. It is also found that these bound excitons originate from only the lowest conduction band and the top three valence bands. To incorporate the lattice screening, we include a scaling factor (ε∞/ε 0 ) 2 , which allows us to obtain exciton binding energies of the correct order of magnitude similar to those in GaN. Here, the excitons are related to each of the three fold split valence bands and the splittings of the latter are also studied as function of strain. Finally, a relation between the anisotropic effective masses and the valence band splitting is pointed out and explained.

36 MATERIALS SCIENCE↗

Distinctive orbital anisotropy observed in the nematic state of a FeSe thin film

Nematic state, where the system is translationally invariant but breaks the rotational symmetry, has drawn great attentions recently due to experimental observations of such a state in both cuprates and iron-based superconductors. The mechanism of nematicity that is likely tied to the pairing mechanism of high-T c , however, still remains controversial. Here, we studied the electronic structure of multilayer FeSe film by angle-resolved photoemission spectroscopy (ARPES). We found that the FeSe film enters the nematic state around 125 K, while the electronic signature of long range magnetic order has not been observed down to 20K indicating the non-magnetic origin of the nematicity. The band reconstruction in the nematic state is characterized by the splitting of the d xz and d yz bands. More intriguingly, such energy splitting is strong momentum dependent with the largest band splitting of ~80 meV at the zone corner. The simple on-site ferro-orbital ordering is insufficient to reproduce the nontrivial momentum dependence of the band reconstruction. Instead, our results suggest that the nearest-neighbor hopping of d xz and d yz is highly anisotropic in the nematic state, the origin of which holds the key in understanding the nematicity in iron-based superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dressed photon-orbital states in a quantum dot: Intervalley spin resonance

Because of the symmetry in silicon quantum wells, silicon quantum dots have an extra degree of freedom leading to a small energy splitting called the valley splitting. This degree of freedom has been viewed alternately as a hazard, especially when the lowest valley-orbit splitting is small compared to the thermal energy, or as an asset, most prominently in proposals to use the valley degree of freedom itself as a qubit. Here we present experiments in which microwave electric field driving induces transitions between both valley-orbit and spin states. We show that this system is highly nonlinear and can be understood through the use of dressed photon-orbital states, enabling a unified understanding of six resonance lines we observe in these experiments.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Strong Landau-quantization effects in high-magnetic-field superconductivity of a two-dimensional multiple-band metal near the Lifshitz transition

We investigate the onset of superconductivity in a magnetic field for a clean two-dimensional multiple-band superconductor in the vicinity of the Lifshitz transition when one of the bands is very shallow. Due to the small number of carriers in this band, the quasiclassical Werthamer-Helfand approximation breaks down and Landau quantization has to be taken into account. We found that the transition temperature T C2 ( H ) has giant oscillations and is resonantly enhanced at the magnetic fields corresponding to the matching of the chemical potential with the Landau levels in the shallow band. This enhancement is especially pronounced for the lowest Landau level. As a consequence, the reentrant superconducting regions in the temperature-field phase diagram emerge at low temperatures near the magnetic fields at which the shallow-band Landau levels cross the chemical potential. The specific behavior depends on the relative strength of the intraband and interband pairing interactions and the reentrance is most pronounced in the purely interband coupling scenario. The reentrant behavior is suppressed by the Zeeman spin splitting in the shallow band; the separated regions disappear already for very small spin-splitting factors. On the other hand, the reentrance is restored in the resonance cases when the spin-splitting energy exactly matches the separation between the Landau levels. As a result, the predicted behavior may be realized in the gate-tuned FeSe monolayer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Vibrational anomalies in A Fe 2 As 2 ( A =Ca, Sr, and Ba) single crystals

In this study, the detailed behavior of the in-plane infrared-active vibrational modes has been determined in A Fe 2 As 2 ( A =Ca, Sr, and Ba) above and below the structural and magnetic transitions at T N = 172, 195 and 138 K, respectively. Above T N , two infrared-active E u modes are observed. In all three compounds, below T N the low-frequency E u mode is observed to split into upper and lower branches; with the exception of the Ba material, the oscillator strength across the transition is conserved. In the Ca and Sr materials, the high-frequency E u mode splits into an upper and a lower branch; however, the oscillator strengths are quite different. Surprisingly, in both the Sr and Ba materials, below T N the upper branch appears to be either very weak or totally absent, while the lower branch displays an anomalous increase in strength. The frequencies and atomic characters of the lattice modes at the center of the Brillouin zone have been calculated for the high-temperature phase for each of these materials. The high-frequency E u mode does not change in position or character across this series of compounds. Below T N , the E u modes are predicted to split into features of roughly equal strength. Finally, we discuss the possibility that the anomalous increase in the strength of the lower branch of the high-frequency mode below T N in the Sr and Ba compounds, and the weak (silent) upper branch, may be related to the orbital ordering and a change in the bonding between the Fe and As atoms in the magnetically-ordered state.

36 MATERIALS SCIENCE↗

Calculating n-point charge correlations in evolving systems

Background: In dynamic systems such as heavy-ion collisions, charge susceptibilities and local charge correlations change with time. These changes are accompanied by nonlocal correlations, which spread diffusively with time and are constrained by local charge conservation. Such correlations have been measured at the Relativistic Heavy-Ion Collider (RHIC) and have been linked to the chemical evolution of the matter and possible phase separation or critical phenomena. The evolution of such correlations has been modeled for two-point correlators superimposed onto hydrodynamics, but not for three- or four-point correlators. These higher-order correlators represent an essential basis for calculating the kurtosis and skewness of charge distributions measured at RHIC. For short-range correlations, such as those between charges on the same particle, and the associated charge-balance correlations a theoretical formalism was lacking for handling higher-order correlations. In particular, one needed to understand how higher-order charge fluctuations was split onto lower numbers of particles. For example, correlations of order three have contributions with charge all on one particle, split onto two particle, or split onto three. By assuming local chemical equilibrium of short-range correlations, the sources of the various correlators were found to be uniquely determined. The relevant theoretical foundation is presented here, including a diagrammatic technique to correctly account for all the possible terms to higher-order correlations. Here, we found a consistent theoretical treatment and its applicability, viability, and tractability are discussed. The formalism derived here enables realistic and quantitative modeling of the three- and four-point charge correlations necessary for understanding measurements of the kurtosis and skewness of charge distributions at RHIC. These observables have been promoted as signals of phase transitions or measures of chemical evolution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Gluon contribution to open heavy-meson production in heavy-ion collisions

A sizable contribution to heavy-quark production in high-energy hadronic and nuclear collisions comes from heavy quark-antiquark pair production from gluon splitting during the parton shower evolution. Here, we investigate the effect of gluon-medium interaction on open heavy flavor spectra in ultrarelativistic heavy-ion collisions. The interaction of hard gluons and heavy quarks with the hot QCD medium is simulated by utilizing a Langevin transport model that simultaneously incorporates contributions from collisional and radiative processes. It is found that, while the gluon splitting channel has quite an important contribution to the single D-meson production cross section, its influence on the final heavy-meson nuclear modification turns out to be quite modest because the average lifetime of hard gluons is short before splitting into heavy-quark-antiquark pairs during the evolution and propagation of the parton shower.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Ab initio calculation of the potential bubble nucleus Si 34

Background: The possibility that an unconventional depletion (referred to as a “bubble”) happens in the middle of the charge density distribution of certain nuclei due to a purely quantum mechanical effect has attracted theoretical and experimental attention in recent years. Based on a mean-field rationale, a correlation between the occurrence of such a semibubble and an anomalously weak splitting between low angular-momentum spin-orbit partners has been further conjectured. Energy density functional and valence-space shell model calculations have been performed to identify and characterize the best candidates, among which Si 34 appears as a particularly interesting case. While the experimental determination of the charge density distribution of the unstable Si 34 is currently out of reach, ($d, p$) experiments on this nucleus have been performed recently to test the correlation between the presence of a bubble and an anomalously weak 1/2 – –3/2 – splitting in the spectrum of Si 35 as compared to S 37 . Purpose: We study the potential bubble structure of Si 34 on the basis of the state-of-the-art ab initio self-consistent Green's function many-body method. Methods: We perform the first ab initio calculations of Si 34 and S 36 . In addition to binding energies, the first observables of interest are the charge density distribution and the charge root-mean-square radius for which experimental data exist in S 36 . The next observable of interest is the low-lying spectroscopy of Si 35 and S 37 obtained from ($d, p$) experiments along with the spectroscopy of Al 33 and P 35 obtained from knock-out experiments. The interpretation in terms of the evolution of the underlying shell structure is also provided. The study is repeated using several chiral effective field theory Hamiltonians as a way to test the robustness of the results with respect to input internucleon interactions. The convergence of the results with respect to the truncation of the many-body expansion, i.e., with respect to the many-body correlations included in the calculation, is studied in detail. We eventually compare our predictions to state-of-the-art multireference energy density functional and shell model calculations. Results: The prediction regarding the (non)existence of the bubble structure in Si 34 varies significantly with the nuclear Hamiltonian used. Yet, demanding that the experimental charge density distribution and the root-mean-square radius of S 36 be well reproduced, along with Si 34 and S 36 binding energies, only leaves the NNLO sat Hamiltonian as a serious candidate to perform this prediction. In this context, a bubble structure, whose fingerprint should be visible in an electron scattering experiment of Si 34 , is predicted. Furthermore, a clear correlation is established between the occurrence of the bubble structure and the weakening of the 1/2 – –3/2 – splitting in the spectrum of Si 35 as compared to S 37 . Conclusions: The occurrence of a bubble structure in the charge distribution of Si 34 is convincingly established on the basis of state-of-the-art ab initio calculations. This prediction will have to be reexamined in the future when improved chiral nuclear Hamiltonians are constructed. On the experimental side, present results act as a strong motivation to measure the charge density distribution of Si 34 in future electron scattering experiments on unstable nuclei. In the meantime, it is of interest to perform one-neutron removal on Si 34 and S 36 in order to further test our theoretical spectral strength distributions over a wide energy range.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Kleinian black holes

We prove there is a unique vacuum solution in split-signature spacetimes with Kleinian S O ( 2 , 1 ) spherical symmetry. We extend our analysis to accommodate a positive or negative cosmological constant and we prove the Kleinian spherically symmetric solutions to Einstein’s equation are locally isomorphic to the split-signature analogs of Schwarzschild–(anti)-de Sitter or Nariai spacetimes. Our analysis provides a Kleinian extension of Birkhoff’s theorem to metrics with split signature. Axisymmetric vacuum solutions are also considered, including (2,2) signature formulations of the Kerr and Taub-Newman-Unti-Tamburino metrics. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Photon topology

The topology of photons in vacuum is interesting because there are no photons with k = 0, creating a hole in momentum space. We show that while the set of all photons forms a trivial vector bundle $γ$ over this momentum space, the R and L photons form topologically nontrivial subbundles $γ±$ with first Chern numbers ∓2. In contrast, $γ$ has no linearly polarized subbundles, and there is no Chern number associated with linear polarizations. It is a known difficulty that the standard version of Wigner’s little group method produces singular representations of the Poincaré group for massless particles. By considering representations of the Poincaré group on vector bundles we obtain a version of Wigner’s little group method for massless particles which avoids these singularities. Here we show that any massless bundle representation of the Poincaré group can be canonically decomposed into irreducible bundle representations labeled by helicity, which in turn can be associated to smooth irreducible Hilbert space representations. This proves that the R and L photons are globally well defined as particles and that the photon wave function can be uniquely split into R and L components. This formalism offers a method of quantizing the electromagnetic field without invoking discontinuous polarization vectors as in the traditional scheme. We also demonstrate that the spin-Chern number of photons is not a purely topological quantity. Lastly, there has been an extended debate on whether photon angular momentum can be split into spin and orbital parts. Our work explains the precise issues that prevent this splitting. Photons do not admit a spin operator; instead, the angular momentum associated with photons’ internal degree of freedom is described by a helicity-induced subalgebra corresponding to the translational symmetry of $γ$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Automated Parton-Shower Variations in Pythia 8

In the era of precision physics measurements at the LHC, efficient and exhaustive estimations of theoretical uncertainties play an increasingly crucial role. In the context of Monte Carlo (MC) event generators, the estimation of such uncertainties traditionally requires independent MC runs for each variation, for a linear increase in total run time. In this work, we report on an automated evaluation of the dominant (renormalization-scale and nonsingular) perturbative uncertainties in the pythia 8 event generator, with only a modest computational overhead. Each generated event is accompanied by a vector of alternative weights (one for each uncertainty variation), with each set separately preserving the total cross section. Explicit scale-compensating terms can be included, reflecting known coefficients of higher-order splitting terms and reducing the effect of the variations. The formalism also allows for the enhancement of rare partonic splittings, such as g→bb¯ and q→qγ, to obtain weighted samples enriched in these splittings while preserving the correct physical Sudakov factors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Measurement of Suppression of Large-Radius Jets and Its Dependence on Substructure in Pb+Pb Collisions at $\sqrt{s_{NN}}$ = 5.02 TeV with the ATLAS Detector

This letter presents a measurement of the nuclear modification factor of large-radius jets in $\sqrt{s_{NN}}$ = 5.02 TeV Pb + Pb collisions by the ATLAS experiment. The measurement is performed using 1.72 nb –1 and 257 pb –1 of Pb + Pb and pp data, respectively. The large-radius jets are reconstructed with the anti-k t algorithm using a radius parameter of R = 1.0, by reclustering anti-k t R = 0.2 jets, and are measured over the transverse momentum (p T ) kinematic range of 158 < p T <1000 GeV and absolute pseudorapidity |y| < 2.0. The large-radius jet constituents are further reclustered using the k t algorithm in order to obtain the splitting parameters, $\sqrt{d_{12}}$ and Δ⁢R 12 , which characterize the transverse momentum scale and angular separation for the hardest splitting in the jet, respectively. The nuclear modification factor, R AA , obtained by comparing the Pb + Pb jet yields to those in pp collisions, is measured as a function of jet transverse momentum (p T ) and $\sqrt{d_{12}}$ or Δ⁢R 12 . A significant difference in the quenching of large-radius jets having single subjet and those with more complex substructure is observed. Systematic comparison of jet suppression in terms of R AA for different jet definitions is also provided. Presented results support the hypothesis that jets with hard internal splittings lose more energy through quenching and provide a new perspective for understanding the role of jet structure in jet suppression.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Asymmetric equilibrium core structures of pyramidal-II < c + a > dislocations in ten hexagonal-close-packed metals

The structures of pyramidal-II < c + a > dislocations, one of the most important defects in structural hexagonal-close-packed (HCP) metals, have not been fully characterized for many of the HCP metals in use today. Here, we employ ab initio informed phase-field dislocation dynamics to determine the minimum energy structure of pyramidal {1¯1¯22} < 11¯23 > dislocations in ten HCP metals, including Be, Co, Mg, Re, Ti, Zn, Cd, Hf, Y, and Zr. As input for the simulations, we calculate, using first-principles density functional theory, the {1¯1¯22} generalized stacking fault energy (GSFE) curves for all ten metals. From these calculations, it is found that magnetism in Co is necessary for achieving a local minimum in the GSFE curve. We observe in simulations that edge and screw character dislocations split into two partials separated by a low-energy intrinsic stacking fault. The splitting distance is shown to scale inversely with the local minimum energy normalized by the product of its shear modulus and Burgers vector. Interestingly, some HCP metals exhibit an asymmetric structure, with either unequal partial Burgers vectors or widths, in contrast to the symmetric configuration expected from linear elastic dislocation theory. We explain these structures by properties of the local maxima in their GSFE curves. Metals with larger degrees of elastic anisotropy result in dislocations with larger splitting distances than would be expected under the commonly used assumption of elastic isotropy. Furthermore, these findings on the sizes and asymmetry in the structures of pyramidal-II < c + a > dislocations are fundamental to understanding how these dislocations glide and interact or react with other defects when these metals are mechanically strained.

36 MATERIALS SCIENCE↗

Quasiaperiodic grain boundary phases of Σ⁢5 tilt grain boundaries in refractory metals

We report ground-state structures and phase transitions in Σ⁢5⁢[001] tilt grain boundaries (GBs) in body-centered-cubic (bcc) refractory metals Nb, Ta, Mo, and W. Σ⁢5 tilt GBs have been extensively investigated over the past several decades, with their ground-state structure—composed of kite-shaped structural units—previously thought to be well understood. By performing a rigorous GB structure search that optimizes the number of atoms in the boundary core, we predict different quasiaperiodic “split kite” phases analogous to those previously found in GBs in face-centered-cubic metals. Furthermore, our results suggest that complex aperiodic phases of GBs appear to be a general phenomenon, as validated through density functional theory calculations. Moreover, the atoms in the split kite phase demonstrate distinct collective diffusion dynamics. Phase-contrast image simulations of split kites show better agreement with experimental observations, offering an alternative explanation for previous microscopy results and motivating future atomically resolved imaging of the GB structure.

Body-centered cubic↗