Recently, symmetry-broken ground states, such as correlated insulating states, magnetic order and superconductivity, have been discovered in twisted bilayer graphene (tBLG) and twisted trilayer graphene (tTLG) near the so-called magic angles. Understanding the magnetic order in these systems is challenging, however, as atomistic methods become extremely expensive near the magic angle and continuum approaches fail to capture important atomistic details. In this work, we develop an approach to incorporate short-ranged Hubbard interactions self-consistently in a continuum model. In addition, we include long-ranged Coulomb interactions which are known to be important when doping the flat bands of tBLG and tTLG. Therefore, for the first time, magnetic order in moiré graphene multilayers is self-consistently explored in a continuum model with atomistic detail. With this approach, we perform a systematic analysis of the magnetic phase diagram of tBLG as a function of doping level and twist angle, near the magic angle. Our results are consistent with previous perturbative atomistic Hartree+Ucalculations. Furthermore, we investigated magnetic order of tTLG, which were found to be similar to those in tBLG. In the future, the developed continuum model can be utilized to investigate magnetic ordering tendencies from short-range exchange interactions in other moiré graphene multilayers as a function of doping, twist angle, screening environment, among other variables.

Moir & eacute; materials provide a versatile platform for engineering excitons, enabling next-generation optoelectronic applications. Continuum models are widely used to study moir & eacute; excitons due to their efficiency, but they often disagree with ab initio many-body approaches, as seen for intralayer excitons in WS2/WSe2 heterobilayers. Here, we resolve these discrepancies using an atomistic, quantum-mechanical framework based on the Bethe-Salpeter equation with Wannier functions as the electronic structure basis, showing that dielectric screening from hBN encapsulation is essential to reproduce experimentally observed exciton features. Our analysis reveals that exciton behavior emerges from a subtle competition between Wannier and charge-transfer characters, driven by stacking-dependent intralayer bandgap variations and environment-tuned electron-hole interactions. We show that the lowest-energy bright intralayer excitons are Wannier-like in WS2/WSe2 heterobilayers but charge-transfer-like in twisted WSe2 homobilayers, despite comparable moir & eacute; sizes. These results establish atomistic modeling as a powerful tool for understanding and controlling excitonic phenomena in moir & eacute; materials.

Niobium diselenide (NbSe_2) has garnered significant attention due to the coexistence of superconductivity and charge density waves (CDWs) down to the monolayer limit. However, realistic modeling of CDWs-capturing effects such as layer number, twist angle, and strain-remains challenging due to the high computational cost of first-principles methods. Here, we develop a physically informed workflow for training machine-learning interatomic potentials (MLIPs) based on the E(3)-equivariant Allegro architecture, tailored to capture the subtle structural and dynamical signatures of CDWs in mono- and bilayer NbSe_2.We find that while CDW lattice distortions are relatively easy to learn, modeling vibrational properties remains more challenging. It requires targeted dataset design and careful hyperparameter tuning, pushing the boundaries and testing the extensibility of current MLIP frameworks. Our MLIPs enable reliable simulations of commensurate and incommensurate CDW phases, including their sensitivity to dimensionality and stacking, as well as CDW dynamics, phonons, and transition temperatures estimated via the stochastic self-consistent harmonic approximation. This work opens new possibilities for studying and tuning CDWs in NbSe_2 and other two-dimensional systems, with implications for electron-phonon coupling, superconductivity, and advanced materials design.

Twisted bilayers of two-dimensional materials have emerged as a highly tunable platform to study and engineer properties of excitons. However, the atomistic description of these properties has remained a significant challenge as a consequence of the large unit cells of the emergent moiré superlattices. To address this problem, we introduce an efficient atomistic quantum-mechanical approach to solve the Bethe–Salpeter equation that exploits the localization of atomic Wannier functions. We then use this approach to study intra- and interlayer excitons in twisted WS2/WSe2 at a range of twist angles. In agreement with experiment, we find that the optical spectrum exhibits three low-energy peaks for twist angles smaller than 2∘. The energy splitting between the peaks is described accurately. We also find two low-energy interlayer excitons with weak oscillator strengths. Our approach opens up new opportunities for the atomistic design of light-matter interactions in ultrathin materials.

Atomically thin crystals hosting flat electronic bands have recently been identified as a rich playground for exploring and engineering strongly correlated phases. Yet, their variety remains limited, primarily to two-dimensional moiré superlattices. Here, we predict the formation of reversible, electrically induced ultraflat bands and π-electron magnetism in one-dimensional chevron graphene nanoribbons. Our ab initio calculations show that the application of a transverse electric field to these nanoribbons generates a pair of isolated, nearly perfectly flat bands with widths of approximately 1 meV around the Fermi level. Upon charge doping, these flat bands undergo a Stoner-like electronic instability, resulting in the spontaneous emergence of local magnetic moments at the edges of the otherwise nonmagnetic nanoribbon, akin to a one-dimensional spin-1/2 chain. Our findings expand the class of carbon-based nanostructures exhibiting flat bands and establish a novel route for inducing correlated electronic phases in chevron graphene nanoribbons.

We examine the electronic structure of recently fabricated in-plane heterojunctions of zigzag graphene nanoribbons embedded in hexagonal boron nitride. We focus on hitherto unexplored interface configurations in which both edges of the nanoribbon are bonded to the same chemical species, either boron or nitrogen atoms. Using ab initio and mean-field Hubbard model calculations, we reveal the emergence of one-dimensional magnetic conducting channels at these interfaces. These channels originate from the energy shift of the magnetic interface states that is induced by charge transfer between the nanoribbon and hexagonal boron nitride. We further address the response of these heterojunctions to external electric and magnetic fields, demonstrating the tunability of energy and spin splittings in the electronic structure. Our findings establish that zigzag graphene nanoribbon/hexagonal boron nitride heterojunctions are a suitable platform for exploring and engineering spin transport in the atomically thin limit, with potential applications in integrated spintronic devices.

Over the last two decades, following the early developments on maximally localized Wannier functions, an ecosystem of electronic-structure simulation techniques and software packages leveraging the Wannier representation has flourished. This environment includes codes to obtain Wannier functions and interfaces with first-principles simulation software, as well as an increasing number of related post-processing packages. Wannier functions can be obtained for isolated or extended systems (both crystalline and disordered), and can be used to understand chemical bonding, to characterize electric polarization, magnetization, and topology, or as an optimal basis set, providing very accurate interpolations in reciprocal space or large-scale Hamiltonians in real space. In this review, we summarize the current landscape of techniques, materials properties and simulation codes based on Wannier functions that have been made accessible to the research community, and that are now well integrated into what we term a Wannier function software ecosystem. First, we introduce the theory and practicalities of Wannier functions, starting from their broad domains of applicability to advanced minimization methods using alternative approaches beyond maximal localization. Then we define the concept of a Wannier ecosystem and its interactions and interoperability with many quantum simulations engines and post-processing packages. We focus on some of the key properties and capabilities that are empowered by such ecosystem—from band interpolations and large-scale simulations to electronic transport, Berryology, topology, electron-phonon couplings, dynamical mean-field theory, embedding, and Koopmans functionals—concluding with the current status of interoperability and automation. [...]

Work function is a fundamental property of metals and is related to many surface-related phenomena of metals. Theoretically, it can be calculated with a metal slab supercell in density functional theory (DFT) calculations. In this paper, we discuss how the commensurability of atomic structure with the underlying fast Fourier transform (FFT) grid affects the accuracy of work function obtained from plane-wave pseudopotential DFT calculations. We show that the macroscopic average potential, which is an important property in work function calculations under the 'bulk reference' method, is more numerically stable when it is calculated with commensurate FFT grids than with incommensurate FFT grids. Due to the stability of the macroscopic average potential, work function calculated with commensurate FFT grids shows better convergence with respect to basis set size, vacuum length and slab thickness of a slab supercell. After we control the FFT grid commensurability issue in our work function calculations, we obtain well-converged work functions for Al, Pd, Au and Pt of (100), (110) and (111) surface orientations. For all the metals considered, the ordering of our calculated work functions of the three surface orientations agrees with experiment. Our findings reveal the importance of the FFT grid commensurability issue, which is usually neglected in practice, in obtaining accurate metal work functions, and are also meaningful to other DFT calculations which can be affected by the FFT grid commensurability issue.

Twisted bilayers of two-dimensional materials have emerged as a highly tunable platform for studying broken symmetry phases. While most interest has been focused on emergent states in systems whose constituent monolayers do not feature broken symmetry states, assembling monolayers that exhibit ordered states into twisted bilayers can also give rise to interesting phenomena. Here, we use large-scale first-principles density-functional theory calculations to study the atomic structure of twisted bilayer $\mathrm{{N}b{S}e_2}$ whose constituent monolayers feature a charge density wave. We find that different charge density wave states coexist in the ground state of the twisted bilayer: monolayer-like $3\times 3$ triangular and hexagonal charge density waves are observed in low-energy stacking regions, while stripe charge density waves are found in the domain walls surrounding the low-energy stacking regions. These predictions, which can be tested by scanning tunneling microscopy experiments, highlight the potential to create complex charge density wave ground states in twisted bilayer systems and can serve as a starting point for understanding superconductivity occurring at low temperatures.

We present theoretical calculations of the optical spectrum of monolayer MoS$_2$ with a charged defect. In particular, we solve the Bethe-Salpeter equation based on an atomistic tight-binding model of the MoS$_2$ electronic structure which allows calculations for large supercells. The defect is modelled as a point charge whose potential is screened by the MoS$_2$ electrons. We find that the defect gives rise to new peaks in the optical spectrum approximately 100-200 meV below the first free exciton peak. These peaks arise from transitions involving in-gap bound states induced by the charged defect. Our findings are in good agreement with experimental measurements.

Armchair graphene nanoribbons are a highly promising class of semiconductors for all-carbon nanocircuitry. Here, we present a new perspective on their electronic structure from simple model Hamiltonians and ab initio calculations. We focus on a specific set of nanoribbons of width n=3p+2, where n is the number of carbon atoms across the nanoribbon axis and p is a positive integer. We demonstrate that the energy-gap opening in these nanoribbons originates from the breaking of a previously unidentified hidden symmetry by long-ranged hopping of π electrons and structural distortions occurring at the edges. This hidden symmetry can be restored or manipulated through the application of in-plane lattice strain, which enables continuous energy-gap tuning, the emergence of Dirac points at the Fermi level, and topological quantum phase transitions. Our work establishes an original interpretation of the semiconducting character of armchair graphene nanoribbons and offers guidelines for rationally designing their electronic structure.

We investigate the effect of thermal fluctuations on the atomic and electronic structure of a twisted MoSe2/WSe2 heterobilayer using a combination of classical molecular dynamics and ab initio density functional theory calculations. Our calculations reveal that thermally excited phason modes give rise to an almost rigid motion of the moiré lattice. Electrons and holes in low-energy states are localized in specific stacking regions of the moiré unit cell and follow the thermal motion of these regions. In other words, charge carriers surf phason waves that are excited at finite temperatures. We also show that such surfing survives in the presence of a substrate and frozen potential. This effect has potential implications for the design of charge and exciton transport devices based on moiré materials.

Half-metals have been envisioned as active components in spintronic devices by virtue of their completely spin-polarized electrical currents. Actual materials hosting half-metallic phases, however, remain scarce. Here, we predict that recently fabricated heterojunctions of zigzag nanoribbons embedded in two-dimensional hexagonal boron nitride are half-semimetallic, featuring fully spin-polarized Dirac points at the Fermi level. The half-semimetallicity originates from the transfer of charges from hexagonal boron nitride to the embedded graphene nanoribbon. These charges give rise to opposite energy shifts of the states residing at the two edges, while preserving their intrinsic antiferromagnetic exchange coupling. Upon doping, an antiferromagnetic-to-ferrimagnetic phase transition occurs in these heterojunctions, with the sign of the excess charge controlling the spatial localization of the net magnetic moments. Our findings demonstrate that such heterojunctions realize tunable one-dimensional conducting channels of spin-polarized Dirac fermions seamlessly integrated into a two-dimensional insulator, thus holding promise for the development of carbon-based spintronics.

In graphene, charged defects break the electron-hole symmetry and can even give rise to exotic collapse states when the defect charge exceeds a critical value which is proportional to the Fermi velocity. In this work, we investigate the electronic properties of twisted bilayer graphene (tBLG) with charged defects using tight-binding calculations. Like monolayer graphene, tBLG exhibits linear bands near the Fermi level but with a dramatically reduced Fermi velocity near the magic angle (approximately 1.1°). This suggests that the critical value of the defect charge in magic-angle tBLG should also be very small. We find that charged defects give rise to significant changes in the low-energy electronic structure of tBLG. Depending on the defect position in the moiré unit cell, it is possible to open a band gap or to induce an additional flattening of the low-energy valence and conduction bands. Our calculations suggest that the collapse states of the two monolayers hybridize in the twisted bilayer. However, their in-plane localization remains largely unaffected by the presence of the additional twisted layer because of the different length scales of the moiré lattice and the monolayer collapse state wavefunctions. These predictions can be tested in scanning tunneling spectroscopy experiments.

Suppression of superconductivity in favour of a striped phase, and its coincidence with a structural transition from a low-temperature orthorhombic (LTO) to a low-temperature tetragonal (LTT) phase, is a ubiquitous feature of hole-doped lanthanum cuprates. We study the effect of anisotropic strain on this transition using density-functional theory on both La$_2$CuO$_4$ and the recently-synthesised surrogate La$_2$MgO$_4$ to decouple electronic and structural effects. Strikingly, we find that compressive strain applied diagonally to the in-plane metal-oxygen bonds dramatically stabilises the LTO phase. Given the mutual exclusivity of 3D superconductivity and long-range static stripe order, we thereby suggest a structural mechanism for understanding experimentally-observed trends in the superconducting $T_{\mathrm{c}}$ under uniaxial pressure, and suggest principles for tuning it.

We discuss the effect of long-range interactions within the self-consistent Hartree-Fock (HF) approximation in comparison to short-range atomic Hubbard interactions on the band structure of twisted bilayer graphene (TBG) at charge neutrality for various twist angles. Starting from atomistic calculations, we determine the quasi-particle band structure of TBG with Hubbard interactions for various magnetic orderings: modulated anti-ferromagnetic (MAFM), nodal anti-ferromagnetic (NAFM) and hexagonal anti-ferromagnetic (HAFM). Then, we develop an approach to incorporate these magnetic orderings along with the HF potential in the continuum approximation. Away from the magic angle, we observe a drastic effect of the magnetic order on the band structure of TBG compared to the influence of the HF potential. Near the magic angle, however, the HF potential seems to play a major role on the band structure compared to the magnetic order. These findings suggest that the spin-valley degenerate broken symmetry state often found in HF calculations of charge neutral TBG near the magic angle should favour magnetic order, since the atomistic Hubbard interaction will break this symmetry in favour of spin polarization.

Charge injection at metal/polymer interfaces is a critical process in many technological devices, including high voltage capacitors and cables in which polyolefin materials, such as polyethylene (PE) and polypropylene (PP), are often used as insulation materials. We use simulations based on density-functional theory to study charge injection at aluminum/PE and aluminum/PP interfaces. Specifically, we investigate the influence of incorporating a variety of polar chemical impurities at the PE and PP chain ends on electron and hole injection barriers. Crucially, we account for the effect of thermal disorder by considering ensembles of thousands of interface structures obtained from ab initio molecular dynamics trajectories at 373 K. We show that the mean injection barrier can change by up to 1.1 eV for Al/PE and 0.6 eV for Al/PP, as compared to the pristine case, depending on which chemical impurity is introduced. We also show that the spread of injection barriers from thermal fluctuations also depends strongly on the chemistry of the impurity. The observed trends can be understood with a simple model based on thermal fluctuations of the dipole moment density associated with the chemical impurity at the interface. We further verify this model by considering larger interface models with lower impurity densities. Our results demonstrate that small chemical modifications, which may arise from oxidation, for example, have a significant influence on charge injection barriers in metal/polyolefin interfaces.

Magic-angle twisted trilayer graphene (MATTG) recently emerged as a highly tunable platform for studying correlated phases of matter, such as correlated insulators and superconductivity. Superconductivity occurs in a range of doping levels that is bounded by van Hove singularities, which stimulates the debate of the origin and nature of superconductivity in this material. In this work, we discuss the role of spin-fluctuations arising from atomic-scale correlations in MATTG for the superconducting state. We show that in a phase diagram as a function of doping ( ν ) and temperature, nematic superconducting regions are surrounded by ferromagnetic states and that a superconducting dome with T c ≈ 2 K appears between the integer fillings ν = −2 and ν = −3. Applying a perpendicular electric field enhances superconductivity on the electron-doped side which we relate to changes in the spin-fluctuation spectrum. We show that the nematic unconventional superconductivity leads to pronounced signatures in the local density of states detectable by scanning tunneling spectroscopy measurements.

We undertake a detailed high-resolution diffraction study of a novel plain band insulator, La$_2$MgO$_4$, which may be viewed as a structural surrogate system of the undoped end-member of the high-T$_c$ superconductors, La$_{2-x-y}$A$^{2+}_x$RE$^{3+}_y$CuO$_{4}$ (A = Ba, Sr, RE= Rare Earth). We find that La$_2$MgO$_4$ exhibits the infamous low-temperature orthorhombic (LTO) to low-temperature tetragonal (LTT) phase transition that has been linked to the suppression of superconductivity in a variety of underdoped cuprates, including the well known La$_{2-x}$Ba$_{x}$CuO$_4$ ($x=0.125$). Furthermore, we find that the LTO-to-LTT phase transition in La$_2$MgO$_4$ occurs for an octahedral tilt angle in the 4 $^{\circ}$ to 5 $^{\circ}$ range, similar to that which has previously been identified as a critical tipping point for superconductivity in these systems. We show that this phase transition, occurring in a system lacking spin correlations and competing electronic states such as charge-density waves and superconductivity, can be understood by simply navigating the density-functional theory ground-state energy landscape as a function of the order parameter amplitude. This result calls for a careful re-investigation of the origins of the phase transitions in high-T$_c$ superconductors based on the hole-doped, $n = 1$ Ruddelsden-Popper lanthanum cuprates.

We investigate the chirality of phonon modes in twisted bilayer WSe2 and demonstrate distinct chiral behavior of the K/K' valley phonons for twist angles close to 0 degrees and close to 60 degrees. In particular, multiple chiral nondegenerate K/K' valley phonons are found for twist angles near 60 degrees whereas no nondegenerate chiral modes are found for twist angles close to 0 degrees. Moreover, we discover two sets of emergent chiral valley modes that originate from an inversion symmetry breaking at the moire scale and find similar modes in moire patterns of strain-engineered bilayers WSe2 and MoSe2/WSe2 heterostructures. At the energy gap between acoustic and optical modes, the formation of flat phonon bands for a broad range of twist angles is observed in twisted bilayer WSe2. Our findings are relevant for understanding electron-phonon and exciton-phonon scattering in moire materials and also for the design of phononic analogues of flat band electrons.