Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “second order phase transitions”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Flare Statistics for Young Stars from a Convolutional Neural Network Analysis of TESS Data

All-sky photometric time-series missions have allowed for the monitoring of thousands of young (t age < 800 Myr) stars in order to understand the evolution of stellar activity. In this work, we developed a convolutional neural network (CNN), stella, specifically trained to find flares in Transiting Exoplanet Survey Satellite (TESS) short-cadence data. We applied the network to 3200 young stars in order to evaluate flare rates as a function of age and spectral type. The CNN takes a few seconds to identify flares on a single light curve. We also measured rotation periods for 1500 of our targets and find that flares of all amplitudes are present across all spot phases, suggesting high spot coverage across the entire surface. Additionally, flare rates and amplitudes decrease for stars t age > 50 Myr across all temperatures T eff ≥ 4000 K, while stars from 2300 ≤ T eff < 4000 K show no evolution across 800 Myr. Stars of T eff ≤ 4000 K also show higher flare rates and amplitudes across all ages. We investigate the effects of high flare rates on photoevaporative atmospheric mass loss for young planets. In the presence of flares, planets lose 4%–7% more atmosphere over the first 1 Gyr. stella is an open-source Python toolkit hosted on GitHub and PyPI.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2020

Flare Statistics for Young Stars from a Convolutional Neural Network Analysis of TESS Data

All-sky photometric time-series missions have allowed for the monitoring of thousands of young (t(age) < 800 Myr) stars in order to understand the evolution of stellar activity. Here, we developed a convolutional neural network (CNN), stella, specifically trained to find flares in Transiting Exoplanet Survey Satellite (TESS) short-cadence data. We applied the network to 3200 young stars in order to evaluate flare rates as a function of age and spectral type. The CNN takes a few seconds to identify flares on a single light curve. We also measured rotation periods for 1500 of our targets and find that flares of all amplitudes are present across all spot phases, suggesting high spot coverage across the entire surface. Additionally, flare rates and amplitudes decrease for stars t(age) > 50 Myr across all temperatures T(eff) ≥ 4000 K, while stars from 2300 ≤ T(eff) < 4000 K show no evolution across 800 Myr. Stars of T(eff) ≤ 4000 K also show higher flare rates and amplitudes across all ages. We investigate the effects of high flare rates on photoevaporative atmospheric mass loss for young planets. In the presence of flares, planets lose 4%–7% more atmosphere over the first 1 Gyr. stella is an open-source Python toolkit hosted on GitHub and PyPI.

Adina D. Feinstein↗
DOE OSTI2025

Strange quark matter as dark matter: 40 yr later, a reappraisal

ABSTRACT Forty years ago Witten suggested that dark matter could be composed of macroscopic clusters of strange quark matter. This idea was very popular for several years, but it dropped out of fashion once lattice quantum chromodynamics calculations indicated that the confinement/deconfinement transition, at small baryonic chemical potential, is not first order, which seemed to be a crucial requirement in order to produce large clusters of quarks. Here, we revisit the conditions under which strangelets can be produced in the Early Universe. We discuss the impact of an instability in the hadronic phase separating a low density, positive-strange-charge phase from a high-density phase with a negative strange charge. This second phase can rapidly stabilize by forming colour-superconducting gaps. The strangelets then undergo partial evaporation. In this way, we obtain distributions of their sizes in agreement with the observational constraints and we discuss the many astrophysical and cosmological implications of these objects. Finally, we examine the most promising techniques to detect this type of strangelets. We also show that strangelets can exist with masses $\lesssim $1017 g, while primordial black holes are ruled out in that mass range, allowing us to distinguish between these two dark matter candidates.

Clemente, Francesco Di (ORCID:0000000282573819)↗
DOE OSTI2025

Bayesian quantification of observability and equation of state of twin stars

The possibility of discovering twin stars, two neutron stars (NSs) with the same mass but different radii, is usually studied in forward modelings by using a restricted number of NS matter equations of state (EOSs) encapsulating a first-order phase transition from hadronic to quark matter (QM). Informing our likelihood function with the NS radius data from GW170817 and using a metamodel with nine parameters capable of mimicking most NS EOSs available in the literature, we conduct a Bayesian quantification of the observability and underlying EOSs of twin stars. Of the accepted EOSs, between 12 and 18% yield twin stars, depending on the restrictions we place on the second branch. The possibility of twin stars remains robust even under recent observational constraints. Here, we show that many of these twin star scenarios are observable with currently available levels of accuracy in measuring NS radii. We also present the marginalized posterior probability density functions (PDFs) of every EOS parameter for each of four mass-radius correlation topologies. We find that the inferred EOS depends sensitively on not only whether twin stars are present, but also the category of twin stars, indicating that the observation of twin stars would provide a strong constraint on the underlying EOS. In particular, for two coexisting hybrid stars having QM cores at different densities, the PDF for QM speed of sound squared 𝑐$^2_{qm}$ has two peaks, one below and another above the conformal limit 𝑐$^2_{qm}$ = 1/3 predicted by perturbative QCD.

QCD in nuclear reactions↗
DOE OSTI2025

Nontrivial critical behavior at magnetic transitions: A case study of Sm 7 ⁢Pd 3

We present a comprehensive analysis of the critical behavior of Sm 7 ⁢Pd 3 in the vicinity of its second-order magnetoelastic transition at 𝑇 c =173 K. The critical exponents (CEs) 𝛽 and 𝛾, determined using both the standard convergence procedure and the average normalized slope (ANS) method, diverge at 𝑇 c –-a characteristic typically associated with first-order transitions. Notably, none of the established universality classes satisfactorily describe the critical behavior of Sm 7 ⁢Pd 3 , and we discuss the possible origins of this deviation in the context of the strong spin-lattice coupling intrinsic to the sample. We emphasize the importance of accurately selecting the critical temperature and magnetic field ranges to ensure robust critical behavior analysis, and propose a quantitative approach to assess the reliability of the extracted CEs. Additionally, we demonstrate that in the ANS method, the critical exponents 𝛽 and 𝛾 should be calculated separately using data for 𝑇 ⩽ 𝑇 c and 𝑇 ⩾ 𝑇 c , respectively. In conclusion, our findings underscore the need for a revised theoretical framework to accurately describe second-order magnetoelastic transitions.

Ferrimagnets↗
DOE OSTI2025

Full complex amplitude control of second-harmonic generation via electrically tunable intersubband polaritonic metasurfaces

Nonlinear intersubband polaritonic metasurfaces based on coupling of the intersubband nonlinear optical response of quantum-engineered semiconductor heterostructures and electromagnetic modes of nanoresonators provide efficient frequency mixing with moderate pump intensities. The resonant nonlinear optical response, represented as a complex function, can be modulated via Stark tuning of intersubband transition energies under applied voltages. However, achieving full complex amplitude control (both phase and magnitude) remains challenging. In this work, we present and experimentally validate electrically tunable nonlinear intersubband polaritonic metasurfaces that achieve complete complex amplitude control for second-harmonic generation (SHG). Through a design featuring two in-plane flipped meta-atoms per unit cell, we achieve complete electrical control of both the amplitude and phase of the metasurface second-order nonlinear susceptibility, with a tuning range of 0 to 30 nm V −1 for the magnitude and 0-2π for the phase of the nonlinear optical response. Using these properties, we achieve complete on-off SHG modulation and beam diffraction tuning through electrically controlled amplitude and phase gratings.

Science & Technology - Other Topics↗
DOE OSTI2022

Wetting behavior of polyelectrolyte complex coacervates on solid surfaces

The wetting behavior of complex coacervates underpins their use in many emerging applications of surface science, particularly wet adhesives and coatings. Many factors dictate if a coacervate phase will condense on a solid surface, including solution conditions, the nature of the polymer–substrate interaction, and the underlying supernatant–coacervate bulk phase behavior. In this work, we use a simple inhomogeneous mean-field theory to study the wetting behavior of complex coacervates on solid surfaces both off-coexistence (wetting transitions) and on-coexistence (contact angles). Here we focus on the effects of salt concentration, the polycation/polyanion surface affinity, and the applied electrostatic potential on the wettability. We find that the coacervate generally wets the surface via a first order wetting transition with second order transitions possible above a surface critical point. Applying an electrostatic potential to a solid surface always improves the surface wettability when the polycation/polyanion–substrate interaction is symmetric. For asymmetric surface affinity, the wettability has a nonmonotonic dependence with the applied potential. We use simple scaling and thermodynamic arguments to explain our results.

36 MATERIALS SCIENCE↗
NASA NTRS1974

A second-order all-digital phase-locked loop

A simple second-order digital phase-locked loop has been designed to synchronize itself to a square-wave subcarrier. Analysis and experimental performance are given for both acquisition behavior and steady-state phase error performance. In addition, the damping factor and the noise bandwidth are derived analytically. Although all the data are given for the square-wave subcarrier case, the results are applicable to arbitrary subcarriers that are odd symmetric about their transition region.

Holmes, J. K.↗
DOE OSTI2021

Phase transitions and critical behavior in hadronic transport with a relativistic density functional equation of state

We develop a flexible, relativistically covariant parametrization of the dense nuclear matter equation of state suited for inclusion in computationally demanding hadronic transport simulations. Within an implementation in the hadronic transport code smash, we show that effects due to bulk thermodynamic behavior are reproduced in dynamic hadronic systems, demonstrating that hadronic transport can be used to study critical behavior in dense nuclear matter, both at and away from equilibrium. We also show that two-particle correlations calculated from hadronic transport simulation data follow theoretical expectations based on the second-order cumulant ratio, and constitute a clear signature of the crossover region above the critical point.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Microchannel-based Membrane-less Extraction of Li from Unconventional Lithium Sources & the Separation of REE

This final report provides an overview of the Project's entire duration, covering July 1, 2021 to December 31, 2023. It primarily focuses on the achievements, technological developments, and unique challenges the team faced while working on separating and extracting Lithium from produced waters. The project's primary aim was to create an integrated, high-throughput, membrane-less, and modular microfluidic platform that could extract Lithium from unconventional sources. We have successfully met all goals and milestones envisioned in the SOPO document. The most critical primary milestones, including the Go-No-Go milestone (refer to the Gantt chart in the Appendices), were successfully accomplished. We demonstrated phase separation (>90%) and extraction (>85%) performance in the MPSE using synthetic, and representative produced water composition feed at 50 ml/min total flow through MPSE 36. We have also performed a parametric study of the MPSE operations, beyond the scope of SOPO, exploring operating conditions of current and broader interest. The extended investigation of operational parameters is concurrent with our efforts to seek further development of the MPSE technology beyond the scope of the Project. Along these lines of development, we have made efforts to be responsive to DOE calls for technological developments of other types of resources (beyond PW) for the recovery of Critical Materials and higher TRL development (beyond TRL 4). During the work on this Project, we developed and implemented three innovative technical approaches that emerged from our efforts to successfully meet the Project milestones. The innovative & original technical approaches developed and implemented in this Project are now the contributions to process engineering that could be clearly credited to the Project. First, Convergent Design Approach is a comprehensive feedforward & feedback loop of four design phases: i) design for functionality, ii) design for manufacturing, iii) design for sustainability, and iv) design for market. Next was Process Intensification. A major aim of this Project was to create an innovative phase separation & extraction microscale-based technology for Li separation – thus the words microchannel-based in the Project title. A microscale-based technology is intrinsically in the center of the Process Intensification domain as defined by its unique principles. Therefore, Process Intensification was implicitly envisioned in the Project’s SOPO. Lastly, Time Scale Analysis is a novel tool for discovering the needs and directions of Process Intensification implementations in any process technology. This Project is fully credited for developing and implementing the three novel technical approaches mentioned above. These are general contributions to process engineering that emerged from this Project. Beyond the original SOPO scope, the OSU-U.Pitt research group utilized a Convergent Design methodology, integrating first-principles mathematical modeling with experimental validation on the Minimum Development Vehicle. By creating these Digital Twins, the team rapidly assessed manufacturing iterations to support TEA analysis. This framework further enabled the development of advanced Surface Modification Techniques, where hydrophobic and oleophobic coating strategies were optimized via Digital Twin tools and validated through rigorous 100-hour longevity testing. TEA Analysis: The closing efforts of this Project were focused on the TEA analysis. TEA analysis had two primary functions: i) enabling critical assessments of design variations withing 10 the Concurrent Design Approach, thus enabling evolution of the MPSE design to reach faster- better-cheaper alternatives; and ii) to create a bridge between the accomplishments of this Project and future projects of higher TRL, beyond TRL 6 level. It is important to note that the TEA model created in the Project stirred the technological solutions for the recovery of critical materials toward a vision of a very profitable modular plant that has unique zero-waste water discharge signature. More importantly, thanks to our experimental performance data and conservative assumptions, the TEA model predicts minimal technological and investment risks. Low cost of a modular unit of a nominal capacity of [1000 tons of Li 2 CO 3 /year] positions the MPSE based technology within the reach of community investors, thus offering a paradigm shift in the development of critical technologies. The project successfully navigated two primary challenges: solvent selection and manufacturing adaptation. Restricted by the SOPO to existing literature for lithium recovery, the team identified a critical need for a "material excellence program" to develop next-generation solvents, eventually concluding with a preliminary investigation into promising Ionic Liquids (ILs). Simultaneously, COVID-19 supply chain disruptions forced a pivot from traditional manufacturing to advanced additive methods at ATAMI-OSU. By transitioning from stainless steel to 3D-printed polymer substrates, the team achieved a transformative three-order-of- magnitude reduction in manufacturing costs and compressed prototyping timelines from several months to just two days. The MPSE technology offers significant energy, environmental, and economic advantages by overcoming the traditional bottlenecks of phase-separation hardware and contactor size. Unlike conventional mixer-settlers or membrane-based systems, MPSE operates without moving parts or fouling-prone membranes, achieving robust performance even with challenging, viscous, or particulate-heavy feeds. Key performance metrics include an energy intensity reduction of 5–50x (3–40 kJ/m 3 ) compared to incumbent technologies and a dramatic reduction of processing time to under 60 seconds, which drastically reduces the physical plant footprint. These technical efficiencies translate into superior economic outcomes; for a 100 t/year Li 2 CO 3 facility, implementing MPSE is projected to nearly halve contactor CAPEX (from $\$$6.08M to $\$$3.01M) and significantly increase the project's Net Present Value (NPV), derisking new investment and enabling distributed critical-mineral processing configurations. The commercialization of MPSE technology is being spearheaded by Vigsur Dynamics Inc., which has adopted a structured, parallel approach to technical and business development since its formation in January 2026. Following extensive customer discovery and engagement with the Oregon State University accelerator, Vigsur Dynamics is working to establish a business model that transitions from pilot demonstrations to modular hardware sales, ultimately aiming for a "build-own-operate" service strategy. Current technical milestones—including 100 hours of continuous operation, superior energy efficiency, and successful 6-unit modular scale-up— provide a foundation for this transition. Backed by ongoing IP licensing and a growing network of industrial and venture advisors, the company is actively de-risking the platform to replace conventional mixer-settler systems in the critical minerals market.

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

Ferroelastic Domain-Induced Electronic Modulation in Halide Perovskites

Lead halide perovskites have emerged as promising materials for optoelectronic applications due to their exceptional properties. In the all-inorganic CsPbBr 3 perovskites, ferroelastic domains formed during phase transitions enhance bulk transport and emissive efficiency. However, the microscopic mechanisms governing carrier dynamics remain poorly understood. In this study, we employ cathodoluminescence (CL) and micro-Raman spectroscopy to image and investigate the electronic properties of the ferroelastic domain walls in CsPbBr 3 single crystals. CL measurements reveal a reduced emissive yield and a slight redshift in emission at the domain walls. Further, micro-Raman studies provide spatially resolved mapping of vibrational modes, exhibiting second-order phonon modes localized at the domain boundaries. In conclusion, our findings suggest that electron–phonon coupling at the twin domain walls plays a critical role in facilitating efficient charge separation, thereby improving the optoelectronic performance of the CsPbBr 3 perovskites.

Carrier dynamics↗
DOE OSTI2022

Nematicity and nematic fluctuations in iron-based superconductors

The reduction of rotational symmetry in a crystalline solid driven by an electronic mechanism is referred to as electronic nematicity. This phenomenon – initially thought to be rare – has by now been observed in an increasing number of systems. Here, the iron-based superconductors present an ideal material platform to study nematicity in crystalline solids. Their nematic transition is pronounced, it can be studied with a wide range of experimental techniques; it is easily tunable; and high-quality samples are widely available. As research on nematicity in iron-based materials is now in its second decade, it builds on tremendous progress in theoretical concepts and experimental techniques. Thus, a stage has been reached at which the nematic phase can be addressed with confidence in its full complexity, including momentum-, time- and material- dependence of the order parameter. Important open questions concern the mechanism by which nematicity affects superconducting pairing and normal-state properties, with a central role for the challenging issue of nematic quantum criticality.

Electronic properties and materials↗
DOE OSTI2021

Asymmetric Character of the Ferroelectric Phase Transition and Charged Domain Walls in a Hybrid Improper Ferroelectric

Abstract In improper ferroelectrics, the spontaneous ordering is typically driven by a structural distortion or a magnetic spin alignment. The induced electric polarization is only a secondary effect. This dependence is a rich source for unusual phenomena and ferroelectric domain configurations for proper, polarization‐driven ferroelectrics. This study focuses on the polar domain structure and the hysteretic behavior at the ferroelectric phase transition in Ca 3 Mn 1.9 Ti 0.1 O 7 as a representative of the recently discovered hybrid improper ferroelectric class of multiferroics. Combining optical second harmonic generation and Raman spectroscopy gives access to the spontaneous structural distortion and the resulting improper electric polarization. This study shows that hybrid improper ferroelectrics contrast proper and improper ferroelectrics in several ways. Most intriguingly, adjacent ferroelectric domains favor head‐to‐head and tail‐to‐tail domain walls over charge‐neutral configurations. Furthermore, the phase transition occurs in an asymmetric fashion. The regime of phase coexistence of the nonpolar and polar phases shows a clear and abrupt upper temperature limit. In contrast, the coexistence toward low temperatures is best described as a fade‐out process, where 100‐nm‐sized islands of the nonpolar phase expand deep into the polar phase.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Low-Temperature Structural Transition in Canfieldite, Ag 8 SnS 6 , Single Crystals

Canfieldite, Ag 8 SnS 6 , is a semiconducting mineral notable for its high ionic conductivity, photosensitivity, and low thermal conductivity. In this paper, we report the solution growth of large single crystals of Ag 8 SnS 6 of mass up to 1 g from a ternary Ag–Sn–S melt. On cooling from high temperature, Ag 8 SnS 6 undergoes a known cubic ($F\bar{4}3m$) to orthorhombic ($Pna2_1$) phase transition at ≈460 K. By studying the magnetization and thermal expansion between 5–300 K, we discover a second structural transition at ≈120 K. Single crystal X-ray diffraction reveals the low-temperature phase adopts a different orthorhombic structure with space group $Pmn2_1$ ($\textit{a}$ = 7.662 9(5) Å, $\textit{b}$ = 7.539 6(5) Å, $\textit{c}$ = 10.630 0(5) Å, Z = 2 at 90 K) that is isostructural to the room-temperature forms of the related Se-based compounds Ag 8 SnSe 6 and Ag 8 GeSe 6 . The 120 K transition is first-order and has a large thermal hysteresis. On the basis of the magnetization and thermal expansion data, the room-temperature polymorph can be kinetically arrested into a metastable state by rapidly cooling to temperatures below 40 K. We last compare the room- and low-temperature forms of Ag8SnS6 with its argyrodite analogues, Ag 8 TQ 6 ($\textit{T}$ = Si, Ge, Sn; $\textit{Q}$ = S, Se), and identify a trend relating the preferred structures to the unit cell volume, suggesting smaller phase volume favors the $Pna2_1$ arrangement. We support this picture by showing that the transition to the $Pmn2_1$ phase is avoided in Ge alloyed Ag 8 Sn 1–x Ge x S 6 samples as well as in pure Ag 8 GeS 6 .

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

Induced ferroelectric phases in SrTiO 3 by a nanocomposite approach

Inducing new phases in thick films via vertical lattice strain is one of the critical advantages of vertically aligned nanocomposites (VANs). In SrTiO 3 (STO), the ground state is ferroelastic, and the ferroelectricity in STO is suppressed by the orthorhombic transition. In this study, we explore whether vertical lattice strain in three-dimensional VANs can be used to induce new ferroelectric phases in SrTiO 3 :MgO (STO:MgO) VAN thin films. The STO:MgO system incorporates ordered, vertically aligned MgO nanopillars into a STO film matrix. Strong lattice coupling between STO and MgO imposes a large lattice strain in the STO film. We have investigated ferroelectricity in the STO phase, existing up to room temperature, using piezoresponse force microscopy, phase field simulation and second harmonic generation. We also serendipitously discovered the formation of metastable TiO nanocores in MgO nanopillars embedded in the STO film matrix. Our results emphasize the design of new phases via vertical epitaxial strain in VAN thin films.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Anisotropic magnetism of polymorphic ErAl 3

ErAl 3 can form in either a trigonal (α) or cubic (β) polymorph and this paper investigates the physical properties of these polymorphs through characterizations of single crystals grown in an aluminum flux. Here, we demonstrate that polymorph selection can be achieved based on the nominal composition of the crystal growth. Magnetic measurements confirm that both β-ErAl 3 and α-ErAl 3 order antiferromagnetically at low temperatures. β-ErAl 3 undergoes antiferromagnetic ordering at a Néel temperature T N = 5.1 K, and the transition is suppressed continually with applied field. α-ErAl 3 displays more complex behavior, with successive magnetic transitions at T N = 5.7 K and T 2 = 4.6 K for zero field, where heat capacity and dilatometry measurements evidence that these transitions are second and first order, respectively. Under magnetic field, strong anisotropy is revealed in α-ErAl 3 , with several steplike metamagnetic transitions observed below T 2 for H ∥ c. These transitions produce sequential magnetization plateaus near one-half of the apparent saturation magnetization. The electrical resistivity of α-ErAl 3 is strongly coupled to its magnetism. At T = 2 K, we observe a positive magnetoresistance reaching 60%, with distinct anomalies at the metamagnetic transitions. The results are summarized in H-T phase diagrams that demonstrate complex magnetic behavior for α-ErAl 3 , suggesting an important role of competing interactions in this metallic system that possesses characteristics of Ising physics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Twisted bilayer graphene. IV. Exact insulator ground states and phase diagram

Here, we derive the exact insulator ground states of the projected Hamiltonian of magic-angle twisted bilayer graphene (TBG) flat bands with Coulomb interactions in various limits, and study the perturbations away from these limits. We define the (first) chiral limit where the AA stacking hopping is zero, and a flat limit with exactly flat bands. In the chiral-flat limit, the TBG Hamiltonian has a U(4) × U(4) symmetry, and we find that the exact ground states at integer filling – 4 ≤ $\textit{ν}$ ≤ 4 relative to charge neutrality are Chern insulators of Chern numbers $ν_C$ = 4 – |$\textit{ν}$ |, 2 – |$\textit{ν}$ |, $\cdots$, |$\textit{ν}$ | – 4, all of which are degenerate. This confirms recent experiments where Chern insulators are found to be competitive low-energy states of TBG. When the chiral-flat limit is reduced to the nonchiral-flat limit which has a U(4) symmetry, we find $\textit{ν}$ = 0 , ± 2 has exact ground states of Chern number 0, while $\textit{ν}$ = ±1, ± 3 has perturbative ground states of Chern number $ν_C$ = ± 1, which are U(4) ferromagnetic. In the chiral-nonflat limit with a different U(4) symmetry, different Chern number states are degenerate up to second-order perturbations. In the realistic nonchiral-nonflat case, we find that the perturbative insulator states with Chern number $ν_C$ = 0 (0 < |$ν_C$| < 4 – |$\textit{ν}$ |) at integer fillings $\textit{ν}$ are fully (partially) intervalley coherent, while the insulator states with Chern number |$ν_C$| = 4 – |$\textit{ν}$ | are valley polarized. However, for 0 < |$ν_C$| ≤ 4 – |$\textit{ν}$ |, the fully intervalley coherent states are highly competitive (0.005 meV/electron higher). At nonzero magnetic field |$\textit{B}$ | > 0, a first-order phase transition for $\textit{ν}$ = ± 1, ± 2 from Chern number $ν_C$ = sgn ($\textit{νB}$)(2 – |$\textit{ν}$ |) to $ν_C$ = sgn ($\textit{νB}$)(4 – |$\textit{ν}$ |) is expected, which agrees with recent experimental observations. Lastly, the TBG Hamiltonian reduces into an extended Hubbard model in the stabilizer code limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗