Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flat”

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,675 records · Page 93

First-principles calculation of gate-tunable ferromagnetism in magic-angle twisted bilayer graphene under pressure

Magic-angle twisted bilayer graphene (MATBG) is notable as a highly tunable platform for investigating strongly correlated phenomena such as unconventional superconductivity and quantum spin liquids, due to easy control of doping level through gating and sensitive dependence of the magic angle on hydrostatic pressure. Experimental observations of correlated insulating states, unconventional superconductivity and ferromagnetism in MATBG indicate that this system exhibits rich exotic phases. In this work, using density functional theory calculations in conjunction with the effective screening medium method, we find the MATBG under pressure at a twisting angle of 2.88 °and simulate how its electronic states evolve when doping level and electric field perpendicular to plane are tuned by gating. Our calculations show that, at doping levels between two electrons and four holes per moiré unit cell, a ferromagnetic (FM) solution with spin density localized at AA stacking sites is lower in energy than the nonmagnetic solution. The magnetic moment of this FM state decreases with both electron and hole doping and vanishes at four electrons/holes doped per moiré unit cell. Hybridization between the flat bands at the Fermi level and the surrounding dispersive bands can take place at finite doping. On increasing the out-of-plane electric field at zero doping, a transition from the FM state to the nonmagnetic one is seen. An investigation of impurity effects shows that both absorption of O 2 molecules and occurrence of Stone–Wales impurities suppress the FM state, and the mechanisms are understood from our calculations. We also analyze the interlayer bonding character due to flat bands via Wannier functions. Finally, we report trivial band topology of the flat bands in the FM state at a certain doping level.

Physics↗

DESI 2024 VI: cosmological constraints from the measurements of baryon acoustic oscillations

We present cosmological results from the measurement of baryon acoustic oscillations (BAO) in galaxy, quasar and Lyman-α forest tracers from the first year of observations from the Dark Energy Spectroscopic Instrument (DESI), to be released in the DESI Data Release 1. DESI BAO provide robust measurements of the transverse comoving distance and Hubble rate, or their combination, relative to the sound horizon, in seven redshift bins from over 6 million extragalactic objects in the redshift range 0.1 < z < 4.2. To mitigate confirmation bias, a blind analysis was implemented to measure the BAO scales. DESI BAO data alone are consistent with the standard flat ΛCDM cosmological model with a matter density Ω m =0.295±0.015. Paired with a baryon density prior from Big Bang Nucleosynthesis and the robustly measured acoustic angular scale from the cosmic microwave background (CMB), DESI requires H 0 =(68.52±0.62) km s -1 Mpc -1 . In conjunction with CMB anisotropies from Planck and CMB lensing data from Planck and ACT, we find Ω m =0.307± 0.005 and H 0 =(67.97±0.38) km s -1 Mpc -1 . Extending the baseline model with a constant dark energy equation of state parameter w, DESI BAO alone requirew=-0.99 +0.15 -0.13 . In models with a time-varying dark energy equation of state parametrised by w 0 and w a , combinations of DESI with CMB or with type Ia supernovae (SN Ia) individually prefer w 0 > -1 and w a < 0. This preference is 2.6σ for the DESI+CMB combination, and persists or grows when SN Ia are added in, giving results discrepant with the ΛCDM model at the 2.5σ, 3.5σ or 3.9σ levels for the addition of the Pantheon+, Union3, or DES-SN5YR supernova datasets respectively. For the flat ΛCDM model with the sum of neutrino mass ∑ m ν free, combining the DESI and CMB data yields an upper limit ∑ m ν < 0.072 (0.113) eV at 95% confidence for a ∑ m ν > 0 (∑ m ν > 0.059) eV prior. These neutrino-mass constraints are substantially relaxed if the background dynamics are allowed to deviate from flat ΛCDM.

79 ASTRONOMY AND ASTROPHYSICS↗

The radial phase variation of reversed-shear and toroidicity-induced Alfvén eigenmodes in DIII-D

The eigenfunction of an instability contains information about energy flow in the wave. Here, the amplitude and phase of electron cyclotron emission radiometer data from hundreds of DIII-D reversed shear Alfvén eigenmodes (RSAE) and toroidicity-induced Alfvén eigenmodes (TAE) are analyzed along the outboard horizontal midplane. The radial phase profile can be flat, linearly rising or falling, convex or concave; in other words, a wide variety of shapes is observed. For a particular mode, often the radial phase profile remains approximately constant as the mode evolves in time but sometimes it changes rapidly. Many TAEs and some RSAEs have phase profiles that are rather flat where the mode amplitude is largest but rise steadily by ~2π at large major radius. Rapid phase changes are observed when the frequencies of an RSAE and TAE overlap and the modes couple. The phase profile depends weakly on the fast-ion gradient that would appear in the absence of wave-induced transport. Linear and quadratic fits to the phase profiles, together with many plasma parameters, are assembled into RSAE and TAE databases. In both cases, large variability is observed. For RSAEs, the strongest phase dependencies are on electron temperature T e , RSAE mode frequency, and the density of carbon impurities. For TAEs, the strongest dependencies are on beam power and major radius of the mode. In general, the average RSAE radial phase profile is essentially flat but the TAE profile has non-zero slope and curvature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Time-dependent scenario modeling for the ST-E1 fusion power plant

ST-E1 is a low aspect ratio fusion power plant being designed by Tokamak Energy targeting 1.5 GW of fusion power. Characterization of the ST-E1 flat-top scenario is described elsewhere McNamara et al (2026 Nucl. Fusion 66 086008); here we focus on addressing the question of how to ramp-up the ST-E1 plasma from an initial state following breakdown and flux-surface formation to the target flat-top state. Being low-aspect ratio, the available solenoid flux of ST-E1 is limited. Therefore, particular consideration is placed on developing ramp-up scenarios that predominantly use inductive flux provided by external vertical field coils. Through time-dependent modeling with METIS, we show that this is possible when the ramp-up is performed at relatively high plasma density: although auxiliary current drive efficiency is reduced, this is significantly outweighed by (1) higher electron-ion collisional equilibration, (2) higher fusion power once ions become sufficiently hot, (3) higher poloidal beta for increased vertical-field flux, and (4) potentially favorable exhaust compatibilities. Ultimately, we show the target ST-E1 flat-top performance can be reached after a ramp-up period lasting 150 s using less than 40 Vs of solenoid flux (with vertical field providing ∼ 90 Vs of flux). The sensitivity to model assumptions are presented, with the general observation that deleterious effects can be mitigated through minor alterations of the auxiliary power temporal waveform and/or total auxiliary power level. The impact of a solenoid and the auxiliary power mix (electron cyclotron heating only versus electron and ion cyclotron heating) on the ST-E1 ramp-up success are also discussed in appendices. On this latter topic, we show that the effect of direct-ion heating during ramp-up is obscured by the uncertainty in the pedestal dynamics, identifying a clear line of future work required to make a definite decision on the ST-E1 auxiliary power mix.

current ramp-up↗

Characterizing the uppermost 100 m structure of the San Jacinto fault zone southeast of Anza, California, through joint analysis of geological, topographic, seismic and resistivity data

We present results from complementary geological, topographic, seismic and electrical resistivity surveys at the Sagebrush Flat (SGB) site along the Clark fault (CF) strand of the San Jacinto fault zone trifurcation area southeast of Anza, California. Joint interpretation of these data sets, each with unique spatiotemporal sensitivities, allow us to better characterize the shallow (<100 m) fault zone at this structurally complex site. Geological mapping at the surface shows the CF has three main subparallel strands within a <100 m zone with varying degrees of rock damage. These strands intersect units of banded gneiss and tonalite, and various sedimentary units. Near the surface, the weathered but more intact tonalite and gneiss to the southwest have relatively high VP. The low-lying flat sedimentary basins around the two southwestern-most CF strands and elevated damaged gneiss to the northeast have lowest VP <500 m s–1. The high relief of the northeast gneiss unit may in part be explained by its extensive damage and inferred increased relative rock uplift. Resistivity imaging shows the unconsolidated dry basin sediments (maximum >1300 Ohm.m) contrasted against the compacted fine-grained (potentially wet) materials within the CF core and the Bautista Formation (minimum <40 Ohm.m), which is slightly elevated above the flat basins. The inverse relationship between VP (increases) and resistivity (decreases) in the uppermost ~15 m can be characterized as log–log linear with slopes of –2.6 to –4. At depths >30 m, the velocity heterogeneity near the surface merges into larger-scale structures that are generally slower on the northeast side of the CF core compared to the southwest side (as much as ~40 per cent reduction in average VP). A previous study revealed a 20–37 per cent variability in peak ground velocities across the SGB site from local earthquakes. The upper end of that range is associated with the near-surface unconsolidated sedimentary basins and northeast damaged gneiss unit. Preliminary analysis of time-dependent topography mostly shows effects of changing vegetation and anthropogenic activity.

Geochemistry & Geophysics↗

Cosmological constraints from higher redshift gamma-ray burst, H ii starburst galaxy, and quasar (and other) data

ABSTRACT We use higher redshift gamma-ray burst (GRB), H ii starburst galaxy (H iiG), and quasar angular size (QSO-AS) measurements to constrain six spatially flat and non-flat cosmological models. These three sets of cosmological constraints are mutually consistent. Cosmological constraints from a joint analysis of these data sets are largely consistent with currently accelerating cosmological expansion and with cosmological constraints derived from a combined analysis of Hubble parameter (H(z)) and baryon acoustic oscillation (BAO, with Planck-determined baryonic matter density) measurements. A joint analysis of the H(z) + BAO + QSO-AS + H iiG + GRB data provides fairly model-independent determinations of the non-relativistic matter density parameter $\Omega _{\rm m_0}=0.313\pm 0.013$ and the Hubble constant $H_0=69.3\pm 1.2\, \rm {km \, s^{-1} \, Mpc^{-1}}$. These data are consistent with the dark energy being a cosmological constant and with spatial hypersurfaces being flat, but they do not rule out mild dark energy dynamics or a little spatial curvature. We also investigate the effect of including quasar flux measurements in the mix and find no novel conclusions.

79 ASTRONOMY AND ASTROPHYSICS↗

Using Pantheon and DES supernova, baryon acoustic oscillation, and Hubble parameter data to constrain the Hubble constant, dark energy dynamics, and spatial curvature

ABSTRACT We use Pantheon Type Ia supernova (SN Ia) apparent magnitude, DES-3 yr binned SN Ia apparent magnitude, Hubble parameter, and baryon acoustic oscillation measurements to constrain six spatially flat and non-flat cosmological models. These sets of data provide mutually consistent cosmological constraints in the six cosmological models we study. A joint analysis of these data sets provides model-independent estimates of the Hubble constant, $H_0=68.8\pm 1.8\ \rm {km \, s^{-1} \ Mpc^{-1}}$, and the non-relativistic matter density parameter, $\Omega _{\rm m_0}=0.294\pm 0.020$. Although the joint constraints prefer mild dark energy dynamics and a little spatial curvature, they do not rule out dark energy being a cosmological constant and flat spatial hypersurfaces. We also add quasar angular size and H ii starburst galaxy measurements to the combined data set and find more restrictive constraints.

79 ASTRONOMY AND ASTROPHYSICS↗

Sub-percentage measure of distances to redshift of 0.1 by a new cosmic ruler

ABSTRACT Distance-redshift diagrams probe expansion history of the Universe. We show that the stellar mass-binding energy (massE) relation of galaxies proposed in our previous study offers a new distance ruler at cosmic scales. By using elliptical galaxies in the main galaxy sample of the Sloan Digital Sky Survey Data Release 7, we construct a distance-redshift diagram over the redshift range from 0.05 to 0.2 with the massE ruler. The best-fit dark energy density is 0.675 ± 0.079 for flat Λ-cold dark matter (ΛCDM) model, consistent with those by other probes. At the median redshift of 0.11, the median distance is estimated to have a fractional error of 0.34 per cent, much lower than those by supernova (SN) Ia and baryonic acoustic oscillation (BAO) and even exceeding their future capability at this redshift. The above low-$\mathit{ z}$ measurement is useful for probing dark energy that dominates at the late Universe. For a flat dark energy equation of state model (flat wCDM), the massE alone constrains w to an error that is only a factor of 2.2, 1.7, and 1.3 times larger than those by BAO, SN Ia, and cosmic microwave background (CMB), respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

The profiles of bars in spiral galaxies

ABSTRACT We present an analysis of major-axis surface-brightness profiles of bars in a volume-limited sample of 182 barred spiral galaxies, using Spitzer 3.6 $\mu$m images. Unlike most previous studies, we use the entire bar profile, and we classify profiles into four categories. These are ‘Peak+Shoulders’ (P+Sh) – updating the classic ‘flat bar’ profile – and three subtypes of the classic ‘exponential’ profile: (true) Exponential, ‘Two-Slope’ (shallow inner slope + steeper outer slope), and ‘Flat-Top’ (constant inner region, steep outer slope). P+Sh profiles are preferentially found in galaxies with high stellar masses, early Hubble types, red colours, and low gas fractions; the most significant factor is stellar mass, and previous correlations with Hubble type can be explained by the tendency of higher-mass galaxies to have earlier Hubble types. The most common type of non-P+Sh profile is Exponential, followed by Flat-Top profiles; all non-P+Sh profiles appear to have similar distributions of stellar mass, Hubble type, colour, and gas fraction. We also morphologically classify the bars of an inclined subsample into those with and without boxy/peanut-shaped (B/P) bulges; as previously reported, the presence of a B/P bulge is very strong function of stellar mass. Essentially all bars with B/P bulges have P+Sh profiles; we associate the profile shoulders with the outer, vertically thin part of the bar. We find a small number of P+Sh profiles in bars without clear B/P bulges, which may indicate that P+Sh formation precedes the formation of B/P bulges.

Astronomy & Astrophysics↗

The Dark Energy Survey supernova program: a reanalysis of cosmology results and evidence for evolving dark energy with an updated Type Ia supernova calibration

We present improved cosmological constraints from a re-analysis of the Dark Energy Survey (DES) 5-year sample of Type Ia supernovae (DES-SN5YR). This re-analysis includes an improved photometric cross-calibration, recent white dwarf observations to cross-calibrate between DES and low-redshift surveys, retraining the salt3 light-curve model and fixing a numerical approximation in the host-galaxy colour law. Our fully recalibrated sample, which we call DES-Dovekie, comprises ~1600 likely Type Ia SNe from DES and ~200 low-redshift SNe from other surveys. With DES-Dovekie, we obtain Ω m = 0.330 ± 0.015 in flat Lambda-cold dark matter (⁠ΛCDM) which changes Ω m by –0.022 compared to DES-SN5YR. Combining DES-Dovekie with cosmic microwave background data from Planck, Atacama Cosmology Telescope, and South Pole Telescope and the DESI DR2 measurements in a flat CDM cosmology, we find ω 0 = –0.803 ± 0.054 and ω a = –0.72 ± 0.21⁠. Our results hold a significance of 3.2σ, reduced from 4.2σ for DES-SN5YR, to reject the null hypothesis that the data are compatible with the cosmological constant. This significance is equivalent to a Bayesian model preference odds of approximately 5:1 in favour of the flat ω 0 ω a CDM model. Using generally accepted thresholds for model preference, our updated data exhibits only a weak preference for evolving dark energy.

dark energy↗

Unconventional Fractional Phases in Multiband Vortexable Systems

We study topological flat bands with distinct features that deviate from conventional Landau level behavior. We show that even in the ideal quantum geometry limit, moiré flat band systems can exhibit physical phenomena fundamentally different from Landau levels without lattices. In particular, we find new fractional quantum Hall states emerging from multiband vortexable systems, where multiple exactly flat bands appear at the Fermi energy. While the set of bands as a whole exhibits ideal quantum geometry, individual bands separately lose vortexability, and thus making them very different from a stack of Landau levels. At certain filling fractions, we find fractional states whose Hall conductivity deviates from the filling factor. Through careful numerical and analytical studies, we rule out all known mechanisms—such as fractional quantum Hall crystals or separate filling of trivial and topological bands—as possible explanations. Leveraging the exact solvability of vortexable systems, we use analytic Bloch wave functions to uncover the origin of these new fractional states, which arises from the commensurability between the moiré unit cell and the magnetic unit cell of an emergent effective magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Twisted symmetric trilayer graphene. II. Projected Hartree-Fock study

The Hamiltonian of the magic-angle twisted symmetric trilayer graphene (TSTG) can be decomposed into a twisted-bilayer-graphene- (TBG-) like flat band Hamiltonian and a high-velocity Dirac fermion Hamiltonian. We use Hartree-Fock mean field approach to study the projected Coulomb interacting Hamiltonian of TSTG developed in Calugaru et al. [Phys. Rev. B 103, 195411 (2021)] at integer fillings ν = -3, -2, -1, and 0 measured from charge neutrality. We study the phase diagram with w 0 /w 1 , the ratio of AA and AB interlayer hoppings, and the displacement field, which introduces an interlayer potential U and hybridizes the TBG-like bands with the Dirac bands. At small U, we find the ground states at all fillings ν are in the same phases as the tensor products of a Dirac semimetal with the filling ν TBG insulator ground states, which are spin-valley polarized at ν = -3, and fully (partially) intervalley coherent at ν = -2, 0(ν = -1) in the flat bands. An exception is ν = -3 with w 0 /w 1 ≳ 0.7 , which possibly becomes a metal with competing orders at small U due to charge transfers between the Dirac and flat bands. At strong U where the bandwidths exceed interactions, all the fillings ν enter a metal phase with small or zero valley polarization and intervalley coherence. Lastly, at intermediate U, semimetal or insulator phases with zero intervalley coherence may arise for ν = -2, -1, 0. Finally, our results provide a simple picture for the electron interactions in TSTG systems, and reveal the connection between the TSTG and TBG ground states.

2-dimensional systems↗

Hyper Suprime-Cam Year 3 results: Cosmology from cosmic shear two-point correlation functions

Here, we perform a blinded cosmology analysis with cosmic shear two-point correlation functions measured from more than 25 million galaxies in the Hyper Suprime-Cam three-year shear catalog in four tomographic redshift bins ranging from 0.3 to 1.5. After conservative masking and galaxy selection, the survey covers 416 deg 2 of the northern sky with an effective galaxy number density of 15 arcmin –2 over the four redshift bins. The 2PCFs adopted for cosmology analysis are measured in the angular range; 7.1 <θ/arcmin < 56.6 for ξ + and 31.2 < θ/arcmin < 248 for ξ – , with a total signal-to-noise ratio of 26.6. We apply a conservative, wide, flat prior on the photometric redshift errors on the last two tomographic bins, and the relative magnitudes of the cosmic shear amplitude across four redshift bins allow us to calibrate the photometric redshift errors. With this flat prior on redshift errors, we find Ω m = 0.25⁢6$^{+0.056}_{–0.044}$ and S 8 ≡ σ 8 ⁢$\sqrt{Ω_{m}/0.3}$ = 0.76⁢9$^{+0.031}_{–0.034}$ (both 68% C.I.) for a flat Λ cold dark matter cosmology. We find, after unblinding, that our constraint on S 8 is consistent with the Fourier space cosmic shear and the 3 × 2 pt analyses on the same HSC dataset. We carefully study the potential systematics from astrophysical and systematic model uncertainties in our fiducial analysis using synthetic data, and report no biases (including projection bias in the posterior space) greater than 0.5σ in the estimation of S 8 . Our analysis hints that the mean redshifts of the two highest tomographic bins are higher than initially estimated. In addition, a number of consistency tests are conducted to assess the robustness of our analysis. Comparing our result with Planck-2018 cosmic microwave background observations, we find a ~2σ tension for the Λ ⁢CDM model.

79 ASTRONOMY AND ASTROPHYSICS↗

Nature of the Correlated Insulator States in Twisted Bilayer Graphene

We use self-consistent Hartree-Fock calculations performed in the full $\pi$-band Hilbert space to assess the nature of the recently discovered correlated insulator states in magic-angle twisted bilayer graphene (TBG). In this work, we find that gaps between the flat conduction and valence bands open at neutrality over a wide range of twist angles, sometimes without breaking the system’s valley projected $\mathcal{C}_2\mathcal{T}$ symmetry. Broken spin-valley flavor symmetries then enable gapped states to form not only at neutrality, but also at total moiré band filling $\textit{n}$ = ± $\textit{p}/4$ with integer $\textit{p}$ = 1, 2, 3, when the twist angle is close to the magic value at which the flat bands are most narrow. Because the magic-angle flat band quasiparticles are isolated from remote band quasiparticles only for effective dielectric constants larger than ~20, the gapped states do not necessarily break $\mathcal{C}_2\mathcal{T}$ symmetry and as a consequence the insulating states at $\textit{n}$ = ± 1/ 4 and $\textit{n}$ = ± 3/4 need not exhibit a quantized anomalous Hall effect.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Enhanced Coherence in Superconducting Circuits via Band Engineering

In superconducting circuits interrupted by Josephson junctions, the dependence of the energy spectrum on offset charges on different islands is 2e periodic through the Aharonov-Casher effect and resembles a crystal band structure that reflects the symmetries of the Josephson potential. We show that higher-harmonic Josephson elements described by a cos(2φ) energy-phase relation provide an increased freedom to tailor the shape of the Josephson potential and design spectra featuring multiplets of flat bands and Dirac points in the charge Brillouin zone. Flat bands provide noise-insensitive energy levels, and consequently, engineering band pairs with flat spectral gaps can help improve the coherence of the system. We discuss a modified version of a flux qubit that achieves, in principle, no decoherence from charge noise and introduce a flux qutrit that shows a spin-1 Dirac spectrum and is simultaneously quite robust to both charge and flux noise.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Localization and Criticality in Antiblockaded Two-Dimensional Rydberg Atom Arrays

Controllable Rydberg atom arrays have provided new insights into fundamental properties of quantum matter both in and out of equilibrium. In this work, we study the effect of experimentally relevant positional disorder on Rydberg atoms trapped in a 2D square lattice under antiblockade (facilitation) conditions. We show that the facilitation conditions lead the connectivity graph of a particular subspace of the full Hilbert space to form a 2D Lieb lattice, which features a singular flat band. Remarkably, we find three distinct regimes as the disorder strength is varied: a critical regime, a delocalized but nonergodic regime, and a regime with a disorder-induced flat band. The critical regime’s existence depends crucially upon the singular flat band in our model, and is absent in any 1D array or ladder system. We propose to use quench dynamics to probe the three different regimes experimentally.

2-dimensional systems↗

Twisted bilayers of thin film magnetic topological insulators

Twisted bilayer graphene (TBG) near "magic angles" has emerged as a rich platform for strongly correlated states of two-dimensional (2D) Dirac semimetals. Here we show that twisted bilayers of thin film magnetic topological insulators (MTIs) with large in-plane magnetization can realize flat bands near 2D Dirac nodes. Using a simple model for thin films of MTIs, we derive a continuum model for two such MTIs, twisted by a small angle with respect to each other. When the magnetization is in plane, we show that interlayer tunneling terms act as effective SU (2) vector potentials, which are known to lead to flat bands in TBG. We show that by changing the in-plane magnetization, it is possible to tune the twisted bilayer MTI band dispersion to quadratic band touching or to flat bands, similar to the TBG. If realized, this system can be a highly tunable platform for strongly correlated phases of two-dimensional Dirac semimetals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Engineering 2D Square Lattice Hubbard Models in 90° Twisted GeX/SnX (X =S, Se) Moiré Superlattices

Because of the large-period superlattices emerging in moiré two-dimensional (2D) materials, electronic states in such systems exhibit low energy flat bands that can be used to simulate strongly correlated physics in a highly tunable setup. While many investigations have thus far focused on moiré flat bands and emergent correlated electron physics in triangular, honeycomb, and quasi-one-dimensional lattices, tunable moiré realizations of square lattices subject to strong correlations remain elusive. Here, in this work, we propose a feasible scheme to construct moiré square lattice systems by twisting two or more layers of 2D materials in a rectangular lattice by 90°. We demonstrate the concept with twisted GeX/SnX (X =S, Se) moiré superlattices and calculate their electronic structures from first principles. We show that the lowest conduction flat band in these systems can be described by a square lattice Hubbard model with parameters which can be controlled by varying the choice of host materials, number of layers, and external electric fields. In particular, twisted double bilayer GeSe realizes a square lattice Hubbard model with strong frustration due to the next-nearest-neighbor hopping that could host unconventional superconductivity, in close analogy to the Hubbard model for copper-oxygen planes of cuprate high-temperature superconductors. The presented scheme uses 90° twisted 2D materials with rectangular unit cells as a promising platform for realizing the physical phenomena of square lattice Hubbard models, establishing a new route for studying its rich phase diagram of magnetism, charge order, and unconventional superconductivity in a highly tunable setting.

2-dimensional systems↗