Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HoPS”

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 163 records · Page 9

π-Orbital Yin–Yang Kagome bands in anilato-based metal–organic frameworks

π-Orbital bonding plays an important role not only in traditional molecular science and solid-state chemistry but also in modern quantum physics and materials, such as the relativistic Dirac states formed by bonding and antibonding π-bands in graphene. In this work, we disclose an interesting manifestation of π-orbitals in forming the Yin–Yang Kagome bands, which host potentially a range of exotic quantum phenomena. Based on first-principles calculations and tight-binding orbital analyses, we show that the frontier π 2 - and π 3 -orbitals in anilato-based metal–organic frameworks form concurrently a conduction and valence set of Kagome bands, respectively, but with opposite signs of lattice hopping to constitute a pair of enantiomorphic Yin and Yang Kagome bands, as recently proposed in a diatomic Kagome lattice. The twisted configuration of neighboring benzene-derived organic ligands bridged by an octahedrally O-coordinated metal ion is found to play a critical role in creating the opposite sign of lattice hopping for the π 2 - versus π 3 -orbitals. Our finding affords a new material platform to study π-orbital originated quantum chemistry and physics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular engineering of fluoroether electrolytes for lithium metal batteries

Fluoroether solvents are promising electrolyte candidates for high-energy-density lithium metal batteries, where high ionic conductivity and oxidative stability are important metrics for design of new systems. Recent experiments have shown that these performance metrics, particularly stability, can be tuned by changing the fraction of ether and fluorine content. However, little is known about how different molecular architectures influence the underlying ion transport mechanisms and conductivity. Here, we use all-atom molecular dynamics simulations to elucidate the ion transport and solvation characteristics of fluoroether chains of varying length, and having different ether segment and fluorine terminal group contents. The design rules that emerge from this effort are that solvent size determines lithium-ion transport kinetics, solvation structure, and solvation energy. In particular, the mechanism for lithium-ion transport is found to shift from ion hopping between solvation sites located in different fluoroether chains in short-chain solvents, to ion–solvent co-diffusion in long-chain solvents, indicating that an optimum exists for molecules of intermediate length, where hopping is possible but solvent diffusion is fast. Consistent with these findings, our experimental measurements reveal a non-monotonic behavior of the effects of solvent size on lithium-ion conductivity, with a maximum occurring for medium-length solvent chains. Finally, a key design principle for achieving high ionic conductivity is that a trade-off is required between relying on shorter fluoroether chains having high self-diffusivity, and relying on longer chains that increase the stability of local solvation shells.

25 ENERGY STORAGE↗

Machine learning reveals factors that control ion mobility in anti-perovskite solid electrolytes

Solid-state batteries are projected to exhibit improved energy densities and safety compared to liquid-electrolyte-based systems. Consequently, the development of solid electrolytes that can conduct cations at rates comparable to that of liquid electrolytes is an active area of research. To facilitate this search it is desirable to understand the chemical and structural features that control ion transport through a crystalline lattice. Here, machine learning (ML) is used to identify features that influence ion mobility in solids, quantify the relative importance of those features, and demonstrate how the variation of individual features alters ion mobility. Adopting the anti-perovskite lattice as a model system, ML algorithms were trained to predict migration barriers from a data set containing hundreds of barriers calculated using density functional theory. In this study, more than 10 6 feature sets of non-redundant descriptors were evaluated for their ability to predict migration barriers. The present analysis correctly reiterates the importance of several features that are known to influence mobility – for example, anion polarizability – but also identifies others whose importance is not widely recognized. Lattice properties such as hopping distance and channel width exhibit the greatest influence on cation mobility. Individual conditional expectation analyses performed on these features shows that barriers decrease as hopping distance decreases and as channel width increases. The defect formation energy – which is commonly associated with the concentration of defects, but not with their mobility – is identified as an important feature for interstitial migration. In sum, this analysis aids in the design of optimal solid electrolytes by simplifying the multi-dimensional design space to a sub-set of properties that are the most important.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hydrogen in disordered titania: connecting local chemistry, structure, and stoichiometry through accelerated exploration

Hydrogen incorporation in native surface oxides of metal alloys often controls the onset of metal hydriding, with implications for materials corrosion and hydrogen storage. A key representative example is titania, which forms as a passivating layer on a variety of titanium alloys for structural and functional applications. These oxides tend to be structurally diverse, featuring polymorphic phases, grain boundaries, and amorphous regions that generate a disparate set of unique local environments for hydrogen. Here, we introduce a workflow that can efficiently and accurately navigate this complexity. First, a machine learning force field, trained on ab initio molecular dynamics simulations, was used to generate amorphous configurations. Density functional theory calculations were then performed on these structures to identify local oxygen environments, which were compared against experimental observations. Second, to classify subtle differences across the disordered configuration space, we employ a graph-based sampling procedure. Finally, local hydrogen binding energies and hopping kinetics are computed using exhaustive density functional theory calculations on representative configurations. Here, we leverage this methodology to show that hydrogen binding energetics are described by local oxygen coordination, which in turn is affected by stoichiometry, and form the basis of hopping kinetics and diffusion. Together these results imply that hydrogen incorporation and transport in TiO x can be tailored through compositional engineering, with implications for improving the performance and durability of titanium-derived alloys in hydrogen environments.

36 MATERIALS SCIENCE↗

Surface facet dependence of competing alloying mechanisms

Metal alloys are ubiquitous in many branches of heterogeneous catalysis, and it is now fairly well established that the local atomic structure of an alloy can have a profound influence on its chemical reactivity. While these effects can be difficult to probe in nanoparticle catalysts, model studies using well defined single crystal surfaces alloyed with dopants enable these structure–function correlations to be drawn. The first step in this approach involves understanding the alloying mechanism and the type of ensembles formed. In this study, we examined the atomic structure of RhCu single-atom alloys formed on Cu(111), Cu(100), and Cu(110) surfaces. Our results show a striking difference between Rh atoms alloying in Cu(111) vs the more open Cu(100) and Cu(110) surface facets. Unlike Cu(111) on which Rh atoms preferentially place-exchange with Cu atoms in the local regions above step edges leaving the majority of the Cu surface free of Rh, highly dispersed, homogeneous alloys are formed on the Cu(100) and (110) surfaces. These dramatically different alloying mechanisms are understood by quantifying the energetic barriers for atomic hopping, exchange, swapping, and vacancy filling events for Rh atoms on different Cu surfaces through theoretical calculations. Density functional theory results indicate that the observed differences in the alloying mechanism can be attributed to a faster hopping rate, relatively high atomic exchange barriers, and stronger binding of Rh atoms in the vicinity of step edges on Cu(111) compared to Cu(110) and Cu(100). Furthermore, these model systems will serve as useful platforms for examining structure sensitive chemistry on single-atom alloys.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Off-state magnetoresistance in long-channel germanium Schottky-barrier MOSFETs

An increasing magnetic field perpendicular to an undoped semiconductor surface at low temperature is known to strengthen the binding of localized electrons to stationary ions, as the wavefunction's tails evolve from exponential to Gaussian. It is also known that application of a high bias voltage to a depleted semiconductor can liberate bound charge and induce a large drop in electrical resistance. We connect these established results to experimental electrical transport measurements on off-state germanium Schottky-barrier metal–oxide–semiconductor field-effect transistor (MOSFETs) with an aluminum oxide insulating dielectric and platinum germanide contacts. We make measurements at the three distinct orientations of the magnetic field with respect to the substrate and the current. At 6 K, we observe sharp attenuation of current by more than 2 orders of magnitude, within 60 mT, at a crossover magnetic field perpendicular to the substrate. A 1 T magnetic field attenuates the current by more than 4 orders of magnitude. The strength of the attenuation and the value of the crossover field are controlled by both the gate–source and drain–source voltages. The attenuation is much weaker when the magnetic field is parallel to the current. Finally, we orient the magnetic field parallel to the substrate, but perpendicular to the current, allowing us to distinguish charge hopping at the oxide interface from charge hopping in the bulk. In conclusion, this large off-state magnetoresistance can be exploited for cryogenic magnetic- and photo-detection, and for high-bias, low-leakage MOSFETs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Fast lithium ion diffusion in brownmillerite Li x Sr 2 Co 2 O 5

Transition metal oxides not only exhibits novel magnetic properties but also provides outstanding ionic transports. Ionic conductors have great potential for interesting tunable physical properties via ionic liquid gating and novel energy storage applications such as all-solid-state lithium batteries. In particular, low migration barriers and high hopping attempt frequency are the keys to achieve fast ion diffusion in solids. Taking advantage of the oxygen-vacancy channel in Li x Sr 2 Co 2 O 5 , we show that migration barriers of lithium ion are as small as 0.28–0.17 eV depending on the lithium concentration rates. Our first-principles calculation also investigated hopping attempt frequency and concluded the room temperature ionic diffusivity and ion conductivity are high as 10 −7 –10 −6 cm 2 s −1 and 10 −3 –10 −2 Scm −1 , respectively, which outperform most of perovskite-type, garnet-type, and sulfide Li-ion solid-state electrolytes. This work proves Li x Sr 2 Co 2 O 5 as a promising super-ionic conductor.

Crystallographic defects↗

Single molecule tracking of bacterial cell surface cytochromes reveals dynamics that impact long-distance electron transport

Using a series of multiheme cytochromes, the metal-reducing bacterium Shewanella oneidensis MR-1 can perform extracellular electron transfer (EET) to respire redox-active surfaces, including minerals and electrodes outside the cell. While the role of multiheme cytochromes in transporting electrons across the cell wall is well established, these cytochromes were also recently found to facilitate long-distance (micrometer-scale) redox conduction along outer membranes and across multiple cells bridging electrodes. Recent studies proposed that long-distance conduction arises from the interplay of electron hopping and cytochrome diffusion, which allows collisions and electron exchange between cytochromes along membranes. However, the diffusive dynamics of the multiheme cytochromes have never been observed or quantified in vivo, making it difficult to assess their hypothesized contribution to the collision-exchange mechanism. Here, we use quantum dot labeling, total internal reflection fluorescence microscopy, and single-particle tracking to quantify the lateral diffusive dynamics of the outer membrane-associated decaheme cytochromes MtrC and OmcA, two key components of EET in S. oneidensis. We observe confined diffusion behavior for both quantum dot-labeled MtrC and OmcA along cell surfaces (diffusion coefficients D MtrC = 0.0192 ± 0.0018 µm 2 /s, D OmcA = 0.0125 ± 0.0024 µm 2 /s) and the membrane extensions thought to function as bacterial nanowires. We find that these dynamics can trace a path for electron transport via overlap of cytochrome trajectories, consistent with the long-distance conduction mechanism. The measured dynamics inform kinetic Monte Carlo simulations that combine direct electron hopping and redox molecule diffusion, revealing significant electron transport rates along cells and membrane nanowires.

59 BASIC BIOLOGICAL SCIENCES↗

Global phase diagram of D-wave superconductivity in the square-lattice t-J model

The Hubbard and closely related t-J models are exciting platforms for unconventional superconductivity (SC). Through state-of-the-art density matrix renormalization group calculations using the grand canonical ensemble, we address open issues regarding the ground-state phase diagram of the extended t-J model on a square lattice in the parameter regime relevant to cuprate superconductors. On large 8-leg cylinders, we demonstrate that the pure t-J model with only nearest-neighbor hoppings and superexchange interactions, for a wide range of doping ( δ = 0.1 to 0.2 ), hosts robust d-wave superconductivity possibly coexisting with weak unidirectional pair density wave. Furthermore, a small next nearest neighbor hopping t 2 suppresses pair and charge density waves, resulting in a uniform d-wave SC phase in both electron- and hole-doped cuprate model systems. Our work validates the t-J model as a proper minimum model for the emergence of superconductivity in cuprate superconductors.

Chen, Feng (ORCID:000000018196030X)↗

An efficient random-sampling method for calculating double occupancy of Gutzwiller wave function in single-band 1D and 2D lattices

In this paper, we report a random sampling method for computing the expectation value of physical quantities based on the Gutzwiler variational wave function. As the first application, we calculated the double occupancy, which is a critical quantity for under- standing the correlation effects in many-body systems, for single-band 1D and 2D lattices. We demonstrated that the random sampling scheme is more efficient than an existing Metropolis Monte-Carlo algorithm. For the 1D Hubbard model with only nearest-neighbor hopping, our results are almost identical to the exact analytic solution. We have also studied systems to which analytic solutions are not available, including the 1D lattices with next-nearest-neighbor hopping and 2D lattices. In addition, constraints on real-space con gurations can be easily implemented in the current scheme to further improve the Gutzwiller wave function. As an example, we calculated the double occupancy for 1D Hubbard model by applying the constraint that all double-occupied sites are paired with an empty site. With enhanced correlation between double-occupied and empty sites, the constraint results in much improved ground-state energy for 1D Hubbard model with strong on-site repulsion.

74 ATOMIC AND MOLECULAR PHYSICS↗

Estimation of biquadratic and bicubic Heisenberg effective couplings from multiorbital Hubbard models

Abstract We studied a multi-orbital Hubbard model at half-filling for two and three orbitals per site on a two-site cluster via full exact diagonalization, in a wide range for the onsite repulsion U , from weak to strong coupling, and multiple ratios of the Hund coupling J H to U . The hopping matrix elements among the orbitals were also varied extensively. At intermediate and large U , we mapped the results into a Heisenberg model. For two orbitals per site, the mapping is into a S = 1 Heisenberg model where by symmetry both nearest-neighbor ( S i ⋅ S j ) and ( S i ⋅ S j ) 2 are allowed, with respective couplings J 1 and J 2 . For the case of three orbitals per site, the mapping is into a S = 3/2 Heisenberg model with ( S i ⋅ S j ), ( S i ⋅ S j ) 2 , and ( S i ⋅ S j ) 3 terms, and respective couplings J 1 , J 2 , and J 3 . The strength of these coupling constants in the Heisenberg models depend on the U , J H , and hopping amplitudes of the underlying Hubbard model. Our study provides a first crude estimate to establish bounds on how large the ratios J 2 / J 1 and J 3 / J 1 can be. We show that those ratios appear rather limited and, as a qualitative guidance, we conclude that J 2 / J 1 is less than 0.4 and J 3 / J 1 is less than 0.2, establishing bounds on effective models for strongly correlated Hubbard systems. Moreover, the intermediate Hubbard U regime was found to be the most promising to enhance J 2 / J 1 and J 3 / J 1 .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Compressive strain turns 𝑠 ± - into 𝑑-wave pairing in a one-unit-cell La 3 ⁢Ni 2 ⁢O 7 thin film via substrate-induced hole doping

Motivated by recent reports of ambient-pressure superconductivity in La 3 ⁢Ni 2⁢ O 7 films grown on LaSrAlO 4 , we investigate the superconducting instability in a one-unit-cell (1UC) thin film using ab initio and random-phase approximation techniques. Compared to the high-pressure bulk system, the ratio of interlayer 𝑑 3⁢𝑧 2 −𝑟 2 hopping to intralayer 𝑑 𝑥 2 −𝑦 2 hopping is suppressed in the 1UC thin film, and the crystal-field splitting of the 𝑒 𝑔 orbitals is increased. Here, our calculation indicates that spin-fluctuation-driven pairing correlations are weak for the stoichiometric case at ambient pressure, but increase significantly under hole doping. The leading pairing symmetry is also found to change by hole doping. Specifically, we obtain a leading 𝑑 𝑥 2 −𝑦 2 pairing state at moderate hole doping, followed by a 𝑑 𝑥⁢𝑦 state at higher doping. These states are driven by intraband spin-fluctuation scattering within the 𝛾 hole pocket centered around the 𝑀 point, and arise primarily from states in the Ni layer farther from the substrate. These results strongly suggest that the thin-film superconducting samples are hole-doped and that pairing in this system predominantly arises in the layer, as opposed to the interlayer pairing in the pressurized bulk system.

Zhang, Yang [Oak Ridge National Laboratory (ORNL),↗

Four-Dimensional Quasielastic Neutron Scattering (4D-QENS) Analysis of Correlated Ionic Conduction in Sr⁢Cl2

Methods of elucidating the mechanisms of fast-ion conduction in solid-state materials are pivotal for advancements in energy technologies such as batteries, fuel cells, sensors, and supercapacitors. In this study, we examine the ionic conduction pathways in single crystal Sr⁢Cl2, which is a fast-ion conductor above 900 K, using four-dimensional quasielastic neutron scattering (4D-QENS). We explore both coherent and incoherent neutron scattering at temperatures above the transition temperature into the superionic phase to explore the correlated motion of hopping anions. Refinements of the incoherent QENS yield residence times and jump probabilities between lattice sites in good agreement with previous studies, confirming that ionic hopping along nearest-neighbor directions is the most probable conduction pathway. However, the coherent QENS reveals evidence of de Gennes narrowing, indicating the importance of ionic correlations in the conduction mechanism. This highlights the need for improvements both in the theory of ionic transport in fluorite compounds and the modeling of coherent 4D-QENS in single crystals.

Coles, Jared [Materials Science Division, Argonne ↗

General trends of superconducting pairing and magnetic correlations in the Ruddlesden-Popper nickelate 𝑚-layered superconductors La 𝑚+1 ⁢Ni 𝑚 ⁢O 3⁢𝑚+1

Here, we report a comprehensive theoretical analysis of the Ruddlesden-Popper layered nickelates La 𝑚+1 ⁢Ni 𝑚 ⁢O 3⁢𝑚+1 (𝑚 = 1 to 6) under pressure. These materials have recently received significant attention due to the discovery of superconductivity in some nickelates under pressure. Our results suggest that, while these Ruddlesden-Popper layered nickelates display many similarities, they also show noticeable differences. One of the common features of La 𝑚+1 ⁢Ni 𝑚 ⁢O 3⁢𝑚+1 is that the electronic states near the Fermi level are mainly contributed by Ni 3⁢𝑑 orbitals, slightly hybridized with O 2⁢𝑝 orbitals. The Ni 𝑑 3⁢𝑧 2 −𝑟 2 orbitals display bonding-antibonding, or bonding-antibonding-nonbonding, characteristic splittings, depending on the even or odd number of stacking layers 𝑚. In addition, the ratio of the in-plane interorbital hopping between 𝑑 3⁢𝑧 2 −𝑟 2 and 𝑑 𝑥 2 −𝑦 2 orbitals and in-plane intraorbital hopping between 𝑑 𝑥 2 −𝑦 2 orbitals was found to be large in La 𝑚+1 ⁢Ni 𝑚 ⁢O 3⁢𝑚+1 (𝑚 = 1 to 6), and this ratio increases from 𝑚 = 1 to 𝑚 = 6, suggesting that the in-plane hybridization will increase as the layer number 𝑚 increases. In contrast to the dominant 𝑠 ± -wave state driven by spin fluctuations in the bilayer La 3 ⁢Ni 2⁢ O 7 and trilayer La 4 ⁢Ni 3 ⁢O 10 , two nearly degenerate 𝑑 𝑥 2 −𝑦 2 -wave and 𝑠 ± -wave leading states were obtained in the four-layer stacking La 5⁢ Ni 4 ⁢O 13 and five-layer stacking La 6 ⁢Ni 5 ⁢O 16 . The leading 𝑠 ± -wave state was recovered in the six-layer material La 7 ⁢Ni 6 ⁢O 19 with slightly higher calculated pairing strength 𝜆 than that of the 𝑑 𝑥 2 −𝑦 2 -wave state. All this evidence suggests that both 𝑠 ± -wave and 𝑑 𝑥 2 −𝑦 2 -wave channels are strongly competing in the high-order niceklates based on our random-phase approximation calculations. In general, at the level of the random-phase approximation treatment, the superconducting transition temperature 𝑇 𝑐 decreases in stoichiometric bulk systems from the bilayer La 3 ⁢Ni 2 ⁢O 7 to the six-layer La 7 ⁢Ni 6 ⁢O 19 , despite the 𝑚-dependent dominant pairing. Both in-plane and out-of-plane magnetic correlations are found to be quite complex. Within the in-plane direction, we obtained the peak of the magnetic susceptibility at 𝐪 = (0.6⁢𝜋, 0.6⁢𝜋) for La 5 ⁢Ni 4 ⁢O 13 (𝑚 = 4) and La 7 ⁢Ni 6 ⁢O 19 (𝑚 = 6) and at 𝐪 = (0.7⁢𝜋, 0.7⁢𝜋) for La 6⁢ Ni 5⁢ O 16 (𝑚 = 5). Along the out-of-plane direction, four layers are coupled as ↓−↑−↑−↓ in La 5 ⁢Ni 4⁢ O 13 , five layers are coupled as ↑−↑−↓−↑−↑ in La 6 ⁢Ni 5 ⁢O 16 , and six layers are coupled as ↑−↓−↓−↑−↑−↓ in La 7 ⁢Ni 6 ⁢O 19 .

Zhang, Yang [Univ. of Tennessee, Knoxville, TN (Un↗

Extending the Gutzwiller approximation to intersite interactions

In this work, we develop an extension of the Gutzwiller approximation (GA) formalism that includes the effects of Coulomb interactions of arbitrary range (including density density, exchange, pair hopping, and Coulomb-assisted hopping terms). This formalism reduces to the ordinary GA formalism for the multiband Hubbard models in the presence of only local interactions. This is accomplished by combining the 1 / z expansion—where z is the coordination number, and only the leading-order terms contribute in the limit of infinite dimensions—with a P R † P R - I expansion, where P R is the Gutzwiller projector on the site R . Furthermore, the method is conveniently formulated in terms of a Gutzwiller Lagrange function. We apply our theory to the extended single-band Hubbard model. Similarly to the usual Brinkman-Rice mechanism, we find a Mott transition. A valence skipping transition is observed, where the occupation of the empty and doubly occupied states for the Gutzwiller wave function is enhanced with respect to the uncorrelated Slater determinant wave function.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dynamics of almost strong edge modes in spin chains away from integrability

Results are presented for the dynamics of an almost strong edge mode which is the quasistable Majorana edge mode occurring in nonintegrable spin chains. The dynamics of the edge mode is studied using exact diagonalization, and compared with time evolution with respect to an effective semi-infinite model in Krylov space obtained from the recursion method. The effective Krylov Hamiltonian is found to resemble a spatially inhomogeneous Su-Schrieffer-Heeger model where the hopping amplitude increases linearly with distance into the bulk, typical of thermalizing systems, but also has a staggered or dimerized structure superimposed on it. Here, the nonperturbatively long lifetime of the edge mode is shown to be due to this staggered structure which diminishes the effectiveness of the linearly growing hopping amplitude. On taking the continuum limit of the Krylov Hamiltonian, the edge mode is found to be equivalent to the quasistable mode of a Dirac Hamiltonian on a half line, with a mass which is nonzero over a finite distance, before terminating into a gapless metallic bulk. The analytic estimates are found to be in good agreement with the numerically obtained lifetimes of the edge mode.

1-dimensional systems↗

Antiferromagnetic ordering and possible lattice response to dynamic uranium valence in U 3 O 8

Determining the correct electronic structure of U 3 O 8 remains a formidable experimental and theoretical challenge. In the low-temperature phase, two crystallographic U sites are separated into a distinct 2U(V)+1U(VI) oxidation configuration. At low temperatures, the U(V) sites form a distorted honeycomb lattice, but the U(VI) sit on a triangular sublattice, suggesting potential for magnetic frustration effects. The spin configuration of the unpaired f electrons on the U(V) sites is likely antiferromagnetic (AFM) from susceptibility measurements, but this has not been confirmed. Here in this paper, we present a neutron scattering investigation of the structure and dynamics of U 3 O 8 from 1.7 to 600 K. We confirm static AFM ordering onset at between 22 and 25 K, which is present down to at least 1.7 K with AFM peaks corresponding to [0.5 1 1] and [0.5 2 2] in the orthorhombic phase. These measurements rule out static AFM order along the a axis of the Amm2 phase, a configuration previously suggested by theory. Above 100 K a quasielastic scattering channel opens that we speculate arises from a lattice relaxation response to thermally activated electron hopping. This term does not conform to a magnetic form factor, so it is not related to spin relaxations. If correct, this mechanism stabilizes a continuous valence transition from 2U(V)+1U(VI) in the low-temperature (T<600 K) orthorhombic phase to the hexagonal phase that contains only one degenerate U site, wherein the U valence can be dynamically stabilized between U(V)↔U(VI) by phonon-assisted electron hopping.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Twisted bilayer graphene. III. Interacting Hamiltonian and exact symmetries

We derive the explicit Hamiltonian of twisted bilayer graphene (TBG) with Coulomb interaction projected into the flat bands and study the symmetries of the Hamiltonian. First, we show that all projected TBG Hamiltonians can be written as positive semidefinite Hamiltonians, an example of which was found in work by Kang and Vafek [Phys. Rev. Lett. 122, 246401 (2019)]. We then prove that the interacting TBG Hamiltonian exhibits an exact U(4) symmetry in the exactly flat band (nonchiral-flat) limit. We further define, besides a first chiral limit where the AA stacking hopping is zero, a second chiral limit where the AB/BA stacking hopping is zero. In the first chiral-flat limit (or second chiral-flat limit) with exactly flat bands, the TBG is enhanced to have an exact U(4) × U(4) symmetry, whose generators are different between the two chiral limits. While in the first chiral limit and in the nonchiral case these symmetries have been found in work by Bultinck et al. [Phys. Rev. X 10, 031034 (2020)], for the eight lowest bands, we here prove that they are valid for projection into any 8 n max particle-hole symmetric TBG bands, with n max > 1 being the practical case for small twist angles < 1°. Furthermore, in the first or second chiral-nonflat limit without flat bands, an exact U(4) symmetry still remains. We also elucidate the link between the U(4) symmetry presented here and the similar but different U(4) of Kang and Vafek [Phys. Rev. Lett. 122, 246401 (2019)]. Additionally, we show that our projected Hamiltonian can be viewed as the normal-ordered Coulomb interaction plus a Hartree-Fock term from passive bands, and exhibits a many-body particle-hole symmetry which renders the physics symmetric around charge neutrality. We also provide an efficient parametrization of the interacting Hamiltonian. The existence of two chiral limits with an enlarged symmetry suggests a possible duality of the model yet undiscovered.

2-dimensional systems↗