Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “2d quantum materials”

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 559 records · Page 31

Nonlinear Nano‐Imaging of Interlayer Coupling in 2D Graphene‐Semiconductor Heterostructures

The emergent electronic, spin, and other quantum properties of 2D heterostructures of graphene and transition metal dichalcogenides are controlled by the underlying interlayer coupling and associated charge and energy transfer dynamics. However, these processes are sensitive to interlayer distance and crystallographic orientation, which are in turn affected by defects, grain boundaries, or other nanoscale heterogeneities. This obfuscates the distinction between interlayer charge and energy transfer. Here, nanoscale imaging in coherent four‐wave mixing (FWM) and incoherent two‐photon photoluminescence (2PPL) is combined with a tip distance‐dependent coupled rate equation model to resolve the underlying intra‐ and inter‐layer dynamics while avoiding the influence of structural heterogeneities in mono‐ to multi‐layer graphene/WSe 2 heterostructures. With selective insertion of hBN spacer layers, it is shown that energy, as opposed to charge transfer, dominates the interlayer‐coupled optical response. From the distinct nano‐FWM and ‐2PPL tip‐sample distance‐dependent modification of interlayer and intralayer relaxation by tip‐induced enhancement and quenching, an interlayer energy transfer time of τ ET ≈ ( 0 . 35 − 0.15 + 0.65 ) $\tau _{\rm ET} \approx (0.35^{+0.65}_{-0.15})$ ps consistent with recent reports is derived. As a local probe technique, this approach highlights the ability to determine intrinsic sample properties even in the presence of large sample heterogeneity.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Electronically Coupled 2D Polymer/MoS 2 Heterostructures

Emergent quantum phenomena in electronically coupled two-dimensional heterostructures are central to next-generation optical, electronic, and quantum information applications. Tailoring electronic band gaps in coupled heterostructures would permit control of such phenomena and is the subject of significant research interest. Two-dimensional polymers (2DPs) offer a compelling route to tailored band structures through the selection of molecular constituents. However, despite the promise of synthetic flexibility and electronic design, fabrication of 2DPs that form electronically coupled 2D heterostructures remains an outstanding challenge. Here, we report the rational design and optimized synthesis of electronically coupled semiconducting 2DP/2D transition metal dichalcogenide van der Waals heterostructures, demonstrate direct exfoliation of the highly crystalline and oriented 2DP films down to a few nanometers, and present the first thickness-dependent study of 2DP/MoS 2 heterostructures. In this work, control over the 2DP layers reveals enhancement of the 2DP photoluminescence by two orders of magnitude in ultrathin sheets and an unexpected thickness-dependent modulation of the ultrafast excited state dynamics in the 2DP/MoS 2 heterostructure. These results provide fundamental insight into the electronic structure of 2DPs and present a route to tune emergent quantum phenomena in 2DP hybrid van der Waals heterostructures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interfacial Charge Transfer in Lead Sulfide/Cadmium Sulfide Quantum Dot–Monolayer Molybdenum Disulfide Heterostructures

We investigate the ultrafast interfacial charge transfer dynamics in mixed-dimensional PbS/CdS quantum dot (QD)–monolayer MoS 2 heterostructures exhibiting a type-II band alignment. These 0D–2D heterostructures were assembled using oleic acid-capped PbS/CdS QDs with systematically varied core sizes deposited onto molecular beam epitaxy-grown monolayer MoS 2 . Steady-state and time-resolved optical spectroscopies revealed pronounced bidirectional charge transfer between the 0D and 2D components, involving both electron and hole exchange across the interface. The quenching of the monolayer MoS 2 A-exciton bleach and the modifications of its decay dynamics provide direct evidence for efficient hole transfer from photoexcited MoS 2 into the PbS QD core. Concurrently, attenuation of the QD ground-state bleach and the emergence of ultrafast decay components confirmed electron transfer from photoexcited QDs to monolayer MoS 2 . Both hole and electron transfer processes are dependent on the QD core size, with smaller cores enabling faster transfer rates and higher efficiencies due to enhanced band offsets and interfacial coupling. The complementary insights gained from transient absorption and time-resolved photoluminescence measurements clarify the mechanisms governing interfacial carrier exchange and highlight the importance of both hole and electron transfer as key components in the functioning of such 0D–2D heterostructures. This study underscores the potential of 0D–2D semiconductor heterostructures for advanced optoelectronic and light-energy conversion applications.

36 MATERIALS SCIENCE↗

cm 2 -Scale Synthesis of MoTe 2 Thin Films with Large Grains and Layer Control

Owing to the small energy differences between its polymorphs, MoTe 2 can access a full spectrum of electronic states from the 2H semiconducting state to the 1T' semimetallic state and from the T d Weyl semimetallic state to the superconducting state in the 1T' and T d phase at low temperature. Thus, it is a model system for phase transformation studies as well as quantum phenomena such as the quantum spin Hall effect and topological superconductivity. Careful studies of MoTe 2 and its potential applications require large-area MoTe 2 thin films with high crystallinity and thickness control. Here, we present cm 2 -scale synthesis of 2H-MoTe 2 thin films with layer control and large grains that span several microns. Layer control is achieved by controlling the initial thickness of the precursor MoO x thin films, which are deposited on sapphire substrates by atomic layer deposition and subsequently tellurized. Despite the van der Waals epitaxy, the precursor–substrate interface is found to critically determine the uniformity in thickness and grain size of the resulting MoTe 2 films: MoTe 2 grown on sapphire show uniform films while MoTe 2 grown on amorphous SiO 2 substrates form islands. Furthermore, this synthesis strategy decouples the layer control from the variabilities of growth conditions for robust growth results and is applicable to growing other transition-metal dichalcogenides with layer control.

2D materials↗

Angle- and polarization-resolved luminescence from suspended and hexagonal boron nitride encapsulated MoSe 2 monolayers

The polarized photoluminescence from atomically thin transition metal dichalcogenides is a frequently applied tool to scrutinize optical selection rules and valley physics, yet it is known to sensibly depend on a variety of internal and external material and sample properties. In this work, we apply combined angle- and polarization-resolved spectroscopy to explore the interplay of excitonic physics and phenomena arising from the commonly utilized encapsulation procedure on the optical properties of atomically thin M o S e 2 . We probe monolayers prepared in both suspended and encapsulated manners. We show that the hBN encapsulation significantly enhances the linear polarization of exciton photoluminescence emission at large emission angles. This degree of linear polarization of excitons can increase up to ∼ <#comment/> 17 % <#comment/> in the hBN encapsulated samples. As we confirm by finite-difference time-domain simulations, it can be directly connected to the optical anisotropy of the hBN layers. In comparison, the linear polarization at finite exciton momenta is significantly reduced in a suspended M o S e 2 monolayer, and becomes notable only in cryogenic conditions. This phenomenon strongly suggests that the effect is rooted in the k-dependent anisotropic exchange coupling inherent in 2D excitons. Our results have strong implications on further studies on valley contrasting selection rules and valley coherence phenomena using standard suspended and encapsulated samples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Axion insulator state in hundred-nanometer-thick magnetic topological insulator sandwich heterostructures

Abstract An axion insulator is a three-dimensional (3D) topological insulator (TI), in which the bulk maintains the time-reversal symmetry or inversion symmetry but the surface states are gapped by surface magnetization. The axion insulator state has been observed in molecular beam epitaxy (MBE)-grown magnetically doped TI sandwiches and exfoliated intrinsic magnetic TI MnBi 2 Te 4 flakes with an even number layer. All these samples have a thickness of ~ 10 nm, near the 2D-to-3D boundary. The coupling between the top and bottom surface states in thin samples may hinder the observation of quantized topological magnetoelectric response. Here, we employ MBE to synthesize magnetic TI sandwich heterostructures and find that the axion insulator state persists in a 3D sample with a thickness of ~ 106 nm. Our transport results show that the axion insulator state starts to emerge when the thickness of the middle undoped TI layer is greater than ~ 3 nm. The 3D hundred-nanometer-thick axion insulator provides a promising platform for the exploration of the topological magnetoelectric effect and other emergent magnetic topological states, such as the high-order TI phase.

36 MATERIALS SCIENCE↗

Four-Terminal Electrochemistry: A Back-Gate Controls the Electrochemical Potential of a 2D Working Electrode

We demonstrate that ultrathin semiconductor working electrodes integrated into metal–insulator–semiconductor (MIS) stacks are an enabling platform for understanding non-Faradaic semiconductor electrochemistry. Furthermore, 5 nm thick ZnO electrodes were deposited on 30 nm HfO 2 dielectric on a Pd “gate” electrode. Application of a bias V G between the Pd gate and the ZnO electrode causes electrons to accumulate in the ZnO layer as measured by recording the in-plane sheet conductance. By contacting the top surface of the ZnO layer with the electrolyte in a conventional three-electrode electrochemical cell, we show that the gate voltage V G modulates the electrochemical potential V ZnO of the ZnO film with respect to a reference electrode. Electrochemical potential changes ΔV ZnO up to –1 V vs Ag/Ag + are achieved for V G = +7 V. Furthermore, by measuring V ZnO vs V G , we extract the quantum capacitance CQ of the ZnO film as a function of the Fermi-level position, which provides a direct measure of the ZnO electronic density of states (DOS). Finally, we demonstrate that the gated ZnO working electrodes can disentangle the two principal components of electrochemical potential, namely, the Fermi-level shift Δδ and the double-layer charging energy eΔΦ EDL . This disentanglement hinges on a fundamental difference between back-gating and normal electrochemical control, namely, that electrochemical control requires double-layer charging, while back-gate control does not. Collectively, the results show that the backside gate electrode is an effective fourth terminal that enables measurements that are difficult to achieve in conventional three-terminal electrochemical setups.

36 MATERIALS SCIENCE↗

Giant Anisotropic Magnetoresistance in Few-Layer α-RuCl 3 Tunnel Junctions

The spin–orbit-assisted Mott insulator α-RuCl 3 is proximate to the coveted quantum spin liquid (QSL) predicted by the Kitaev model. In the search for the pure Kitaev QSL, reducing the dimensionality of this frustrated magnet by exfoliation has been proposed as a way to enhance magnetic fluctuations and Kitaev interactions. Here, in this work, we perform angle-dependent tunneling magnetoresistance (TMR) measurements on ultrathin α-RuCl 3 crystals with various layer numbers to probe their magnetic, electronic, and crystal structures. We observe a giant change in resistance, as large as ~2500%, when the magnetic field rotates either within or out of the α-RuCl 3 plane, a manifestation of the strongly anisotropic spin interactions in this material. In combination with scanning transmission electron microscopy, this tunneling anisotropic magnetoresistance (TAMR) reveals that few-layer α-RuCl 3 crystals remain in the high-temperature monoclinic phase at low temperatures. It also shows the presence of a zigzag antiferromagnetic order below the critical temperature T N ≃ 14 K, which is twice the one typically observed in bulk samples with rhombohedral stacking. Our work offers valuable insights into the relation between the stacking order and magnetic properties of this material, which helps lay the groundwork for creating and electrically probing exotic magnetic phases such as QSLs via van der Waals engineering.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Imaging two-dimensional generalized Wigner crystals

The Wigner crystal has fascinated condensed matter physicists for nearly 90 years. Signatures of two-dimensional (2D) Wigner crystals were first observed in 2D electron gases under high magnetic field, and recently reported in transition metal dichalcogenide moiré superlattices. Direct observation of the 2D Wigner crystal lattice in real space, however, has remained an outstanding challenge. Conventional scanning tunnelling microscopy (STM) has sufficient spatial resolution but induces perturbations that can potentially alter this fragile state. Here we demonstrate real-space imaging of 2D Wigner crystals in WSe 2 /WS 2 moiré heterostructures using a specially designed non-invasive STM spectroscopy technique. This employs a graphene sensing layer held close to the WSe 2 /WS 2 moiré superlattice. Local STM tunnel current into the graphene layer is modulated by the underlying Wigner crystal electron lattice in the WSe 2 /WS 2 heterostructure. Different Wigner crystal lattice configurations at fractional electron fillings of n = 1/3, 1/2 and 2/3, where n is the electron number per site, are directly visualized. The n = 1/3 and n = 2/3 Wigner crystals exhibit triangular and honeycomb lattices, respectively, to minimize nearest-neighbour occupations. The n = 1/2 state spontaneously breaks the original C3 symmetry and forms a stripe phase. Here our study lays a solid foundation for understanding Wigner crystal states in WSe 2 /WS 2 moiré heterostructures and provides an approach that is generally applicable for imaging novel correlated electron lattices in other systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Combining molecular beam epitaxy and low-energy electron microscopy with in situ magnetic susceptibility measurements within an integrated ultrahigh vacuum system

Quantum two-dimensional materials, including ultrathin superconducting films, are of great current research interest. These films are typically fabricated under ultra-high vacuum (UHV) conditions and are sensitive to the environment—prone to oxidation and contamination when exposed to the atmosphere. This hampers the study of their intrinsic properties by standard ex situ techniques. Here, we present a variable-temperature mutual inductance probe system integrated under UHV with molecular beam epitaxy (MBE) synthesis and low-energy electron microscopy, enabling nondestructive in situ characterization of superconducting thin films. The system employs a reflection-type configuration and reaches a low temperature (∼4 K) using a high-cooling-power, vibration-isolated cryocooler. In conclusion, we demonstrate the system performance by measuring the superconducting critical temperature in a copper-oxide thin film.

2D materials↗

Ionic Modulation at the LaAlO 3 /KTaO 3 Interface for Extreme High-Mobility Two-Dimensional Electron Gas

Due to the coexistence of many emergent phenomena, including 2D superconductivity and a large Rashba spin-orbit coupling, 5d transition metal oxides-based two-dimensional electron systems (2DESs) have been prospected as one of the potential intrants for modern electronics. However, despite the lighter electron mass, the mobility of carriers, a key requisite for high-performance devices, in 5d-oxides devices remains far behind their 3d-oxides analogs. The carriers’ mobility in these oxides is significantly hampered by the inevitable presence of defects. Here, in this study, very high mobility (≈22 650 cm 2 V –1 s –1 ) of 5d-2DES confined at the LaAlO 3 /KTaO 3 interface is reported. The high mobility, which is beyond the values observed in SrTiO 3 2DESs in the same carrier-density range, is achieved using the ionic-liquid gating at room temperature. The authors postulate that the ionic-liquid gating affects the oxygen vacancies and efficiently reduces any disorder at the interface. Investigating density and mobility in a broad range of back-gate voltage, the authors reveal that the mobility follows the power-law µ ∝ n 1.2 , indicating the very high quality of ionic-liquid-gated LaAlO 3 /KTaO 3 devices, consistent with the postulate. Furthermore, the analysis of the quantum oscillations confirms that the high-mobility electrons occupy the electronic sub-bands emerging from the Ta:5d orbitals of KTaO 3 .

36 MATERIALS SCIENCE↗

Non-perturbative cathodoluminescence microscopy of beam-sensitive materials

Cathodoluminescence microscopy is now a wellestablished and powerful tool for probing the photonic properties of nanoscale materials, but in many cases, nanophotonic materials are easily damaged by the electronbeam doses necessary to achieve reasonable cathodoluminescence signal-to-noise ratios. Two-dimensional materials have proven particularly susceptible to beam-induced modifications, yielding both obstacles to high spatial-resolution measurement and opportunities for beam-induced patterning of quantum photonic systems. Here pan-sharpening techniques are applied to cathodoluminescence microscopy in order to address these challenges and experimentally demonstrate the promise of pan-sharpening for minimally-perturbative high-spatial-resolution spectrum imaging of beam-sensitive materials.

2D materials↗

$t_1$-Noise eliminated dipolar heteronuclear multiple-quantum coherence solid-state NMR spectroscopy

Heteronuclear correlation (HETCOR) spectroscopy is one of the key tools in the arsenal of the solid-state NMR spectroscopist to probe chemical and spatial proximities between two different nuclei and enhance spectral resolution. Dipolar heteronuclear multiple-quantum coherence (D-HMQC) is a powerful technique that can be potentially utilized to obtain 1 H detected 2D HETCOR solid-state NMR spectra of any NMR active nucleus. A long-standing problem in 1 H detected D-HMQC solid-state NMR experiments is the presence of t 1 -noise which reduces sensitivity and impedes spectral interpretation. In this contribution, we describe novel pulse sequences, termed t 1 -noise eliminated (TONE) D-HMQC, that minimize t 1 -noise and can provide higher sensitivity and resolution than conventional D-HMQC. Monte-Carlo and numerical simulations confirm that t 1 -noise in conventional D-HMQC primarily occurs because random MAS frequency fluctuations cause variations in the NMR signal amplitude from scan to scan, leading to imperfect cancellation of uncorrelated signals by phase cycling. The TONE D-HMQC sequence uses 1 H π-pulses to refocus the evolution of 1 H CSA across each SR4 2 1 recoupling block, improving the stability of the pulse sequence to random MAS frequency fluctuations. The 1 H refocusing pulses also restore the orthogonality of in-phase and anti-phase magnetization for all crystallite orientations at the end of each recoupling block, enabling the use of 90° flip-back or LG spin-lock trim pulses to reduce the intensity of uncorrelated signals. We demonstrate the application of these methods to acquire 1 H detected 2D 1 H{ 35 Cl} and 1 H{ 13 C} HETCOR spectra of histidine·HCl·H 2 O with reduced t 1 -noise. To show generality, we also apply these methods to obtain 2D 1 H{ 17 O} spectra of 20%- 17 O fmoc-alanine and for the first time at natural abundance, 2D 1 H{ 25 Mg} HETCOR spectra of magnesium hydroxide. The TONE D-HMQC sequences are also used to probe 1 H– 25 Mg and 1 H– 27 Al proximities in Mg–Al layered double hydroxides and confirm the even mixing of Mg and Al in these materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tunable exciton valley-pseudospin orders in moiré superlattices

Excitons in two-dimensional (2D) semiconductors have offered an attractive platform for optoelectronic and valleytronic devices. Further realizations of correlated phases of excitons promise device concepts not possible in the single particle picture. Here we report tunable exciton “spin” orders in WSe 2 /WS 2 moiré superlattices. We find evidence of an in-plane (xy) order of exciton “spin”—here, valley pseudospin—around exciton filling v ex = 1, which strongly suppresses the out-of-plane “spin” polarization. Upon increasing v ex or applying a small magnetic field of ~10 mT, it transitions into an out-of-plane ferromagnetic (FM-z) spin order that spontaneously enhances the “spin” polarization, i.e., the circular helicity of emission light is higher than the excitation. The phase diagram is qualitatively captured by a spin-1/2 Bose–Hubbard model and is distinct from the fermion case. Our study paves the way for engineering exotic phases of matter from correlated spinor bosons, opening the door to a host of unconventional quantum devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Stabilized Synthesis of 2D Verbeekite: Monoclinic PdSe 2 Crystals with High Mobility and In-Plane Optical and Electrical Anisotropy

PdSe 2 has a layered structure with an unusual, puckered Cairo pentagonal tiling. Its atomic bond configuration features planar 4-fold-coordinated Pd atoms and intralayer Se–Se bonds that enable polymorphic phases with distinct electronic and quantum properties, especially when atomically thin. PdSe 2 is conventionally orthorhombic, and direct synthesis of its metastable polymorphic phases is still a challenge. Here, we report an ambient-pressure chemical vapor deposition approach to synthesize metastable monoclinic PdSe 2 . Monoclinic PdSe 2 is shown to be synthesized selectively under Se-deficient conditions that induce Se vacancies. These defects are shown by first-principles density functional theory calculations to reduce the free energy of the metastable monoclinic phase, thereby stabilizing it during synthesis. Further, the structure and composition of the monoclinic PdSe 2 crystals are identified and characterized by scanning transmission electron microscopy imaging, convergent beam electron diffraction, and electron energy loss spectroscopy. Polarized Raman spectroscopy of the monoclinic PdSe 2 flakes reveals their strong in-plane optical anisotropy. Electrical transport measurements show that the monoclinic PdSe 2 exhibits n-type charge carrier conduction with electron mobilities up to ~298 cm 2 V –1 s –1 and a strong in-plane electron mobility anisotropy of ~1.9. The defect-mediated growth pathway identified in this work is promising for phase-selective direct synthesis of other 2D transition metal dichalcogenides.

2D material↗

Origin of Large Effective Phonon Magnetic Moments in Monolayer MoS 2

Recent helicity-resolved magneto-Raman spectroscopy measurement demonstrates large effective phonon magnetic moments of ~2.5 μ B in monolayer MoS 2 , highlighting resonant excitation of bright excitons as a feasible route to activate Γ-point circularly polarized phonons in transition metal dichalcogenides. However, a microscopic picture of this intriguing phenomenon remains lacking. In this work, we show that an orbital transition between the split conduction bands (Δ 0 = 4 meV) of MoS 2 couples to the doubly degenerate E" phonon mode (Ω 0 = 33 meV), forming two hybridized states. Our phononic and electronic Raman scattering measurements capture these two states: (i) one with predominantly phonon contribution in the helicity-switched channels and (ii) one with primarily orbital contribution in the helicity-conserved channels. An orbital-phonon coupling model successfully reproduces the large effective magnetic moments of the circularly polarized phonons and explains their thermodynamic properties. Strikingly, the Raman mode from the orbital transition is superimposed on a strong quasi-elastic scattering background, indicating the presence of spin fluctuations. As a result, the electrons excited to the conduction bands through the exciton exhibit paramagnetic behavior although MoS 2 is generally considered as a nonmagnetic material. By depositing nanometer-thickness nickel thin films on monolayer MoS 2 , we tune the electronic structure so that the A exciton perfectly overlaps with the 633 nm laser. The optimization of resonance excitation leads to pronounced tunability of the orbital-phonon hybridized states. Our results generalize the orbital-phonon coupling model of effective phonon magnetic moments to material systems beyond the paramagnets and magnets.

2D transition metal dichalcogenide↗

Comparing Exciton-Polariton Couplings in Two-Dimensional MoS2 Cavities

Exciton-polaritons (EPs) are an emerging approach to achieve strong light-matter coupling and robust quantum phenomena, with application to areas such as photonic neuromorphic computing for instance reservoir computing. 2D transition metal dichalcogenides (TMDs) are promising candidates to host EPs due to their potential for strong coupling as evidenced by large Rabi splitting energies, which in the case of few-layer MoS2 has been demonstrated to reach 293 meV when coupled to the C exciton, placing it in the ultrastrong coupling regime, at room temperature. Here we examine how the Rabi splitting of CVD-grown monolayer MoS2 in a cavity can be optimized by tuning two variables: cavity geometry (thickness of silver and dielectric spacers) and the targeted exciton for coupling (A at about 1.9 eV, B at about 2.06 eV, or C at about 2.9 eV). With the cavity geometry, there is a tradeoff between optical field confinement and transmission. With the targeted exciton, there may be a tradeoff between Rabi splitting and lifetime. Transfer matrix simulations are used to inform how thick the cavity layers need to be, and at what angle measurements should be performed, to achieve Rabi splitting at a certain exciton energy. Experimentally, the Rabi splitting is measured by placing the cavity at a specified angle and measuring the optical reflectance and transmittance. By fitting the data, we find Rabi splittings of 124 meV, 194 meV, and 288 meV for the A, B, and C excitons, respectively, which places them in the strong, and on the edge of the ultrastrong, coupling regime.

exciton-polaritons↗

An atomic cluster expansion potential for twisted multilayer graphene

Twisted multilayer graphene, characterized by its moiré patterns arising from inter-layer rotational misalignment, serves as a rich platform for exploring quantum phenomena. Machine learning interatomic potentials (MLIPs) are a promising approach to model such systems. Our work develops a method to generate training and test datasets for fitting MLIPs that capture all possible misalignments but remain small-scale to facilitate efficient data generation and parameter estimation. To achieve this, we generate configurations with periodic boundary conditions suitable for density functional theory calculations, and then introduce an internal twist and shift within those supercell structures. Using this technique, supplemented with an active learning workflow, we fit an Atomic Cluster Expansion potential for simulating twisted multilayer graphene and test it for accuracy and robustness on a range of simulation tasks.

2D materials↗