Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multilayer graphene”

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 73 records · Page 4

Current-Driven Magnetization Reversal in Orbital Chern Insulators

Graphene multilayers with flat moiré minibands can exhibit the quantized anomalous Hall effect due to the combined influence of spontaneous valley polarization and topologically nontrivial valley-projected bands. The sign of the Hall effect in these Chern insulators can be reversed either by applying an external magnetic field, or by driving a transport current through the system. Here we propose a current-driven mechanism whereby reversal occurs along lines in the (current I , magnetic-field B ) control parameter space with slope d I / d B = ( e / h ) M A M ( 1 – γ 2 ) / γ , where M is the magnetization, A M is the moiré unit cell area, and γ < 1 is the ratio of the chemical potential difference between valleys along a domain wall to the electrical bias e V .

2-dimensional systems↗

Dimer description of the SU(4) antiferromagnet on the triangular lattice

In systems with many local degrees of freedom, high-symmetry points in the phase diagram can provide an important starting point for the investigation of their properties throughout the phase diagram. In systems with both spin and orbital (or valley) degrees of freedom such a starting point gives rise to SU(4)-symmetric models. Here we consider SU(4)-symmetric "spin" models, corresponding to Mott phases at half-filling, i.e. the six-dimensional representation of SU(4). This may be relevant to twisted multilayer graphene. In particular, we study the SU(4) antiferromagnetic "Heisenberg" model on the triangular lattice, both in the classical limit and in the quantum regime. Carrying out a numerical study using the density matrix renormalization group (DMRG), we argue that the ground state is non-magnetic. We then derive a dimer expansion of the SU(4) spin model. An exact diagonalization (ED) study of the effective dimer model suggests that the ground state breaks translation invariance, forming a valence bond solid (VBS) with a 12-site unit cell. Finally, we consider the effect of SU(4)-symmetry breaking interactions due to Hund's coupling, and argue for a possible phase transition between a VBS and a magnetically ordered state.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Current rectification based on noncentrosymmetric quantum materials

Rectification is a process that converts electromagnetic fields into direct current (DC). Such a process underlies a wide range of technologies, including wireless communication, wireless charging, energy harvesting, and infrared detection. Existing rectifiers are mostly based on semiconductor diodes, with limited applicability to small voltages or high frequency inputs. Here, we present an alternative approach to current rectification that uses the electronic properties of quantum crystals without semiconductor junctions. We identify a new mechanism for rectification from skew scattering due to the chirality of itinerant electrons in time-reversal-invariant but inversion-breaking materials. Our calculations reveal large, tunable rectification effects in graphene multilayers and transition metal dichalcogenides. These effects can be used in high-frequency rectifiers by rational material design and quantum wavefunction engineering.

Isobe, Hiroki↗

Enhanced van der Waals epitaxy of germanium by out-of-plane dipole moment induced from transferred graphene on TiN/AlN multilayers

We report recent advances in 3D/2D heterostructures have opened up tremendous opportunities in building highly flexible and durable optoelectronic devices. However, the inherit lack of interfacial bonding and low surface energy of van der Waals surfaces limit the nucleation and growth of 3D materials. Enhancing wettability by providing a porous buffer is effective in growing compound semiconductors on graphene while van der Waals epitaxy of Ge remains challenging. Here, the nucleation of Ge has been significantly improved from an islanded mode to granular modes by using a TiN/AlN multilayered buffer prior to Ge/graphene integration. Highly textured Ge growth with dominating (111), (220), and (311) peaks are identified by x-ray diffraction. The microstructure of the buffer TiN/AlN demonstrates a polycrystalline quality with clean interfaces between each interlayer and the substrate. Kelvin probe force microscopy measurements along the lateral TiN/AlN interface identify a potential drop corresponding to the AlN phase. This contact potential difference between TiN and AlN is the key in generating the out-of-plane dipole moment that modifies the surface energy of the monolayer graphene, resulting in enhanced wettability of the Ge adatoms nucleated on top. Surface dipole induced nucleation of 3D semiconductor thin films on 2D materials via the proper design of buffer layer is fundamentally important to enhance the 3D/2D growth toward flexible optoelectronic applications.

36 MATERIALS SCIENCE↗

The structural modification and magnetism of many-layer epitaxial graphene implanted with low-energy light ions

Modifying the properties of graphene has gained wide interest for a plethora of potential applications, including spintronics. One approach has demonstrated that proton irradiation can induce ferromagnetism in graphene as well as in graphite. However, little is known about how the protons interact with graphene, the mechanism that creates the ferromagnetism, or whether the protons remain in the graphene. Here we report an investigation, broadly relevant to graphitic carbon, using low-energy (360–2000 eV) ions of hydrogen, deuterium, and helium implanted into multilayer epitaxial graphene. Complementary x-ray and neutron reflectivity demonstrate that essentially all of the implanted hydrogen remains chemisorbed in graphene. In situ x-ray diffraction reveals significantly different rates of interlayer expansion of the multilayer graphene. Analysis of these data demonstrates that the interlayer expansion arises entirely from the interstitials created by the ions and not from hydrogen that remains in the graphene. The results also establish a quantitative measure of the layer expansion due to carbon interstitials. Magnetometry and x-ray diffraction studies show that the magnetic moment relates to the amount of interstitial carbon rather than the amount of hydrogen, demonstrating that the induced room-temperature ferromagnetism arises directly from the disrupted bonding of the carbon lattice.

36 MATERIALS SCIENCE↗

Magic-angle twisted symmetric trilayer graphene as a topological heavy-fermion problem

Recently, Song and Bernevig [Phys. Rev. Lett. 129, 047601 (2022)] reformulated magic-angle twisted bilayer graphene as a topological heavy fermion problem, and used this reformulation to provide a deeper understanding for the correlated phases at integer fillings. Here, in this work, we generalize this heavy-fermion paradigm to magic-angle twisted symmetric trilayer graphene, and propose a low-energy f–c–d model that reformulates magic-angle twisted symmetric trilayer graphene as heavy localized f modes coupled to itinerant topological semimetalic c modes and itinerant Dirac d modes. Our f–c–d model well reproduces the single-particle band structure of magic-angle twisted symmetric trilayer graphene at low energies for displacement field $\mathcal{E}$ ϵ [0,300]⁢ meV. By performing Hartree-Fock calculations with the f–c–d model for v = 0,–1,–2 electrons per Moiré unit cell, we reproduce all the correlated ground states obtained from the previous numerical Hartree-Fock calculations with the Bistritzer-MacDonald-type model, and we find additional new correlated ground states at high displacement field. Based on the numerical results, we propose a simple rule for the ground states at high displacement fields by using the f–c–d model, and provide analytical derivation for the rule at charge neutrality. We also provide analytical symmetry arguments for the (nearly) degenerate energies of the high-$\mathcal{E}$ ground states at all the integer fillings of interest, and make experimental predictions of which charge-neutral states are stabilized in magnetic fields. Our f–c–d model provides a new perspective for understanding the correlated phenomena in magic-angle twisted symmetric trilayer graphene, suggesting that the heavy fermion paradigm of Song and Bernevig [Phys. Rev. Lett. 129, 047601 (2022)] should be the generic underpinning of correlated physics in multilayer moire graphene structures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Frank-van der Merwe growth in bilayer graphene

Bilayer graphene has attracted significant interest because of its unique properties, including fascinating electrical behavior when one layer is slightly rotated relative to the other. However, the quality of large-area bilayer graphene is often limited by the layer-plus-island growth mode in which islands of thicker graphene present as unavoidable impurities. Here, we report the observation of the layer-by-layer, Frank-van der Merwe (FM) growth mode in bilayer graphene where multilayer impurities are suppressed. Instead of the conventional surface adhesive energy, we found it possible to tune interface adhesive energy with an oxidative pretreatment. Finally, the FM-grown bilayer graphene is of AB stacking or with small twisting angle (θ = 0°–5°), which is more mechanically robust compared with monolayer graphene, facilitating a free-standing wet transfer technology.

36 MATERIALS SCIENCE↗

Moiréless correlations in ABCA graphene

Atomically thin van der Waals materials stacked with an interlayer twist have proven to be an excellent platform toward achieving gate-tunable correlated phenomena linked to the formation of flat electronic bands. In this work we demonstrate the formation of emergent correlated phases in multilayer rhombohedral graphene––a simple material that also exhibits a flat electronic band edge but without the need of having a moiré superlattice induced by twisted van der Waals layers. We show that two layers of bilayer graphene that are twisted by an arbitrary tiny angle host large (micrometer-scale) regions of uniform rhombohedral four-layer (ABCA) graphene that can be independently studied. Scanning tunneling spectroscopy reveals that ABCA graphene hosts an unprecedentedly sharp van Hove singularity of 3–5-meV half-width. We demonstrate that when this van Hove singularity straddles the Fermi level, a correlated many-body gap emerges with peak-to-peak value of 9.5 meV at charge neutrality. Mean-field theoretical calculations for model with short-ranged interactions indicate that two primary candidates for the appearance of this broken symmetry state are a charge-transfer excitonic insulator and a ferrimagnet. Finally, we show that ABCA graphene hosts surface topological helical edge states at natural interfaces with ABAB graphene which can be turned on and off with gate voltage, implying that small-angle twisted double-bilayer graphene is an ideal programmable topological quantum material.

42 ENGINEERING↗

Mechanical, electronic, optical, piezoelectric and ferroic properties of strained graphene and other strained monolayers and multilayers: an update

Abstract This is an update of a previous review (Naumis et al 2017 Rep. Prog. Phys. 80 096501). Experimental and theoretical advances for straining graphene and other metallic, insulating, ferroelectric, ferroelastic, ferromagnetic and multiferroic 2D materials were considered. We surveyed (i) methods to induce valley and sublattice polarisation ( P ) in graphene, (ii) time-dependent strain and its impact on graphene’s electronic properties, (iii) the role of local and global strain on superconductivity and other highly correlated and/or topological phases of graphene, (iv) inducing polarisation P on hexagonal boron nitride monolayers via strain, (v) modifying the optoelectronic properties of transition metal dichalcogenide monolayers through strain, (vi) ferroic 2D materials with intrinsic elastic ( σ ), electric ( P ) and magnetic ( M ) polarisation under strain, as well as incipient 2D multiferroics and (vii) moiré bilayers exhibiting flat electronic bands and exotic quantum phase diagrams, and other bilayer or few-layer systems exhibiting ferroic orders tunable by rotations and shear strain. The update features the experimental realisations of a tunable two-dimensional Quantum Spin Hall effect in germanene, of elemental 2D ferroelectric bismuth, and 2D multiferroic NiI 2 . The document was structured for a discussion of effects taking place in monolayers first, followed by discussions concerning bilayers and few-layers, and it represents an up-to-date overview of exciting and newest developments on the fast-paced field of 2D materials.

Physics↗

Design and Characterization of Hybrid Multilayer Structures: Layer-by-Layer Growth of Polymer and Graphene Oxide Assemblies and Their Utility in Fuel Cell Applications

We report the layer-by-layer growth of poly- (ethylenimine) (PEI) that has been modified with sulfamate functionalities (S-PEI) using Zr 4+ -complexation interlayer linking chemistry. We have also deposited adlayers of graphene oxide (GO) that have been modified to possess sulfate functionalities (SGO), onto the S-PEI layers. The multilayer assemblies are formed with sulfamate/sulfate and sulfate/sulfate (S-PEI + S-GO) interlayer linking chemistry. In all cases the adlayer thickness is consistent with predictions based on van der Waals volume and/or molecular mechanics calculations. X-ray photoelectron spectroscopy (XPS) is used to characterize the Zr/S stoichiometry in the multilayer assembly. The utility of these hybrid multilayer structures is demonstrated in a Proton Exchange Membrane (PEM) fuel cell, where they are shown to reduce H 2 gas crossover by 15% with only a 12 nm thick layer.

fuel cell↗

Highly Asymmetric Graphene Layer Doping and Band Structure Manipulation in Rare Earth–Graphene Heterostructure by Targeted Bonding of the Intercalated Gadolinium

Heterostructures consisting of vertically stacked two-dimensional (2D) materials have recently gained large attention due to their highly controllable electronic properties and resulting quantum phases. In contrast to the mechanically stacked multilayered systems, which offer exceptional control over a stacking sequence or interlayer twist angles, the epitaxially grown 2D materials express unprecedented quality and stability over wafer-scale lengths. However, controlling the growth conditions remains a major obstacle toward the formation of complex, epitaxial heterostructures with well-defined electronic properties. Here, we synthesized a trilayer graphene heterostructure on the SiC(0001) substrate with two specific interlayer locations occupied by gadolinium. We applied multitechnique methodology based on low-temperature scanning tunneling microscopy/spectroscopy (STM/S) and angle-resolved photoelectron spectroscopy (ARPES) to determine the intercalant’s locations in the complex, epitaxial graphene heterostructure. Our approach relies on very high quality and large, micrometer-scale homogeneity of the synthesized system. The experimentally determined electronic structure is dominated by the two topmost graphene layers. Overall, our spectroscopic results show quantitative agreement between global ARPES, local STM/S, and density functional theory predictions. The characterized electronic properties primarily reflect highly anisotropic doping levels between the two corresponding graphene layers, which significantly affect the band structure topology. Two pairs of hybridized massive Dirac bands from our initial synthesis–the bilayer graphene on the SiC(0001) substrate–are transformed upon Gd intercalation into two pairs of massless Dirac bands with a new hybridization region in between. Our results open perspectives in the realization of exotic 2D quantum materials via atomically precise synthesis of epitaxial, multilayered graphene–rare earth heterostructures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A combined SECM and electrochemical AFM approach to probe interfacial processes affecting molecular reactivity at redox flow battery electrodes

Redox flow batteries are attractive technologies for grid energy storage since they use solutions of redox-active molecules that enable a superior scalability and the decoupling of power and energy density. However, the reaction mechanisms of the redox active components at RFB electrodes are complex, and there is currently a pressing need to understand how interfacial processes impact the kinetics and operational reversibility of RFB systems. Here, we developed a combined electrochemical imaging methodology rooted in scanning electrochemical microscopy (SECM) and atomic force microscopy (AFM) for exploring the impact of electrode structure and conditioning on the electron transfer properties of model redox-active dialkoxybenzene derivatives, 2,5-di-tert-butyl-1,4-bis(2-methoxyethoxy)benzene (C1) and 2,3-dimethyl-1,4-dialkoxybenzene (C7). Using AFM and secondary-ion mass spectrometry (SIMS), we observed the formation of interfacial films with distinct mechanical properties compared to those of cleaved graphitic surfaces, and exclusively during reduction of electrogenerated radical cations. These films had an impact on the median rate and distribution of the electron transfer rate constant at the basal plane of multilayer and single layer graphene electrodes, displaying kinetically-limited values that did not yield the activation expected per the Butler–Volmer model with a transfer coefficient ~0.5. These changes were dependent on redoxmer structure: furthermore, SECM showed strong attenuation of C7 kinetics by a surface layer on MLG and SLG, while C1 kinetics were only affected by SLG. SECM and AFM results together show that these limiting films operate exclusively on the basal plane of graphite, with the edge plane showing a relative insensitivity to cycling and operation potential. Ultimately, this integrated electrochemical imaging methodology creates new opportunities to understand the unique role of interfacial processes on the heterogeneous reactivity of redoxmers at electrodes for RFBs, with a future role in elucidating phenomena at high active concentrations and spatiotemporal variations in electrode dynamics.

25 ENERGY STORAGE↗

Density functionals combined with van der Waals corrections for graphene adsorbed on layered materials

Standard density functionals like Perdew-Burke-Ernzerhof (PBE) or Strongly Constrained and Appropriately Normed (SCAN) need a correction to account for long-range van der Waals (vdW) interaction. The damped Zaremba-Kohn model (dZK) starts from a formula for the vdW interaction of a distant atom with a solid surface, both with known dielectric properties, damps this formula at short range, and then treats an adsorbed molecule or atomic layer as a collection of renormalized atoms. Here, we extend this model to graphene adsorbed on semiconducting layered materials [bulk graphite and hBN, and multilayer transition metal dichalcogenides (TMD)] by including the C 4 asymptotic term and multiple electrostatic image effects due to the two surfaces of the substrate slabs in the vdW calculations. The resulting SCAN-vdW-dZK and PBE-vdW-dZK give approximately the same results for the systems considered here, in agreement with available reference values. The predicted binding energies are roughly 25% lower than those from SCAN+rVV10 (revised Vydrov and Van Voorhis 2010), and ~15% lower than those from PBE+rVV10. Since SCAN+rVV10 usually overbinds, the predicted binding energies by SCAN-vdW-dZK and PBE-vdW-dZK are expected to be closer to the true values. The predicted equilibrium binding distances from SCAN-vdW-dZK and PBE-vdW-dZK are slightly larger (~ 0.1 Å) than those from SCAN+rVV10, and close to those from SCAN. The binding energy depends upon the number of substrate layers more strongly in vdW-dZK than in rVV10. The C 4 -term contributions can be 40% of the total vdW interactions, and the C 5 term contributes about 10%. The effects of images and the back surfaces of slabs can contribute about 4–10%. The vdW interaction energy power laws from the vdW-dZK model for graphene adsorbed on multilayer MoS 2 show slowly varying decay to the pairwise exponent –4 with increasing separation $\textit{D}$, very similar to those obtained from the random-phase approximation and renormalization group approaches by Ambrosetti et al. Both the PBE-vdW-dZK and SCAN-vdW-dZK give a greater increase in interlayer binding energy when the TMD substrate changes from monolayer to four-layer for the graphene/TMD adsorption systems. This is consistent with the relevant electron-energy-loss spectroscopy experimental results showing increased dielectric response from the substrate with increasing substrate layer number.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

matsim-agents v1.0

matsim-agents is a multi-agent AI framework for atomistic materials simulation and discovery. It orchestrates large language models (LLMs), machine-learned interatomic potentials (MLIPs), and DFT codes into a single agentic loop running on laptops and DOE leadership-class supercomputers. MULTI-AGENT ORCHESTRATION A LangGraph state machine with three nodes: a Planner that converts a natural-language research objective into structured tasks; an Executor that dispatches atomistic tools and loops until the queue is empty; and an Analyst that summarizes results into a human-readable report. State is checkpointed after every step and human-in-the-loop gates can be inserted at any edge. HYPOTHESIS-DRIVEN DISCOVERY CHAT An interactive REPL (matsim-agents chat) that couples LLM dialogue with atomistic simulation. Chemical formulas are automatically detected in conversation turns and trigger a full crystal-phase exploration: structure generation → relaxation → stability scoring → result injection back into the conversation, creating a closed hypothesis-refinement loop. CRYSTAL PHASE ENUMERATION Given a composition, the phase explorer enumerates prototypes by stoichiometry: elemental (fcc/bcc/hcp/sc/diamond), binary 1:1 (rocksalt/CsCl/zincblende/ wurtzite/fluorite/rutile), ternary 1:1:3 (cubic perovskite), ternary 1:2:4 (perovskite + spinel), quaternary 1:1:2:6 (Fm-3m double perovskite). 2-D prototypes (graphene, h-BN, MoS2 2H/1T) and multilayer stacking are also supported via --include-2d and --num-layers. SUPERCELL GENERATION AND SITE DECORATION Auto-tiling to a minimum atom count (--min-atoms), explicit NxNxN tiling (--supercell), symmetry-distinct site decorations (--n-orderings), and isotropic lattice-scale sweeps (--lattice-scales) for volume bracketing. MLFF RELAXATION AND STABILITY SCORING HydraGNN (multi-headed GNN) drives structure relaxation via ASE with FIRE, BFGS, or BFGSLineSearch. Stability output: delta-E/atom ranking across phases and a max-residual-force dynamical-stability proxy. Other MLIPs (MACE, NequIP, Orb) can be plugged in through the same interface. DFT BACKENDS Quantum ESPRESSO pw.x and VASP 6.6 are first-class labellers. Both have validated GPU builds and SLURM/PBS launchers for three DOE platforms: Frontier (AMD MI250X, ROCm), Aurora (Intel PVC, oneAPI), Perlmutter (NVIDIA A100, CUDA). QE produces ~100 binaries (pw.x, ph.x, epw.x, ...). VASP supports scf, relax, vc-relax, and vc-relax-shape run types. ACTIVE-LEARNING LOOP matsim-agents al run CONFIG.yaml drives an iterative HydraGNN-DFT loop: MD generates candidates → ensemble/MC-dropout uncertainty selects the most informative → DFT labels them in parallel inside one allocation → dataset grows → HydraGNN retrains → repeat. DFT backend is a single YAML toggle (dft.backend: vasp | qe). LLM-generated seed structures are supported (no curated POSCAR library needed). Config uses ${VAR}, ${VAR:-default}, ${VAR:?msg} shell-style substitution for cross-user/cross-site portability. LLM BACKENDS Ollama (local, default), vLLM (HPC multi-GPU serving), OpenAI, Anthropic, HuggingFace Transformers+Accelerate. Selected at runtime via flag or env var with no code changes. HPC PORTABILITY Same Python entry points run on Frontier (ROCm 7.2), Aurora (oneAPI), and Perlmutter (CUDA 12). DFT and ML stacks are never co-loaded in the same shell; they couple through the scheduler and filesystem. Advanced multi-node launchers (serve, discovery-chat, single-relaxation, active-learning, QE warm-start) are provided for all three platforms. CODABENCH COMPETITION BUNDLE A self-contained benchmark: 159 atomistic test structures across 11 material classes, 5 tasks (formation energy, forces, ML relaxation, AI-DFT relaxation, phase stability ranking), public/private leaderboard split (30/70), and four ready-to-run baselines: MACE-MP-0, HydraGNN, UMA, AllScAIP.

Lupo Pasini, Massimiliano [Oak Ridge National Labo↗

Domain-Dependent Surface Adhesion in Twisted Few-Layer Graphene: Platform for Moiré-Assisted Chemistry

Twisted van der Waals multilayers are widely regarded as a rich platform to access novel electronic phases thanks to the multiple degrees of freedom available for controlling their electronic and chemical properties. Here, in this work, we propose that the stacking domains that form naturally due to the relative twist between successive layers act as an additional ”knob” for controlling the behavior of these systems and report the emergence and engineering of stacking domain-dependent surface chemistry in twisted few-layer graphene. Using mid-infrared near-field optical microscopy and atomic force microscopy, we observe a selective adhesion of metallic nanoparticles and liquid water at the domains with rhombohedral stacking configurations of minimally twisted double bi- and trilayer graphene. Furthermore, we demonstrate that the manipulation of nanoparticles located at certain stacking domains can locally reconfigure the moiré superlattice in their vicinity at the micrometer scale. Our findings establish a new approach to controlling moiré-assisted chemistry and nanoengineering.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Chemical Vapor Deposition-Grown Hexagonal Boron Nitride and Graphene for Tandem Dielectric Capacitive Polymer Devices

To address the low energy density of polymeric capacitors, this work explores how tandem combinations of antagonistic properties can enhance energy storage. We introduce multilayer architecture integrating chemical vapor deposition (CVD)-grown ≈2 nm hexagonal boron nitride (hBN) and graphene with ferroelectric polyvinylidene fluoride (PVDF) and linear polyetherimide (PEI) layers to suppress leakage paths, reduce loss, and promote Maxwell–Wagner–Sillars polarization. Raman spectroscopy, X-ray diffraction (XRD), atomic force microscopy (AFM), scanning electron microscopy (SEM), and plasma focused ion beam (PFIB) analyses confirm robust material integration. In conclusion, for the tandem device, E BD ≈600 V/μm, while individual layers show lower values (PVDF|Graphene (Gr)|PVDF ≈ 50 V/μm and PEI|hBN|PEI ≈ 275 V/μm), yielding an overall ≈6000% enhancement and demonstrating the effectiveness of the 2D multilayer design.

Chemical vapor deposition (CVD)↗

Hysteresis-free high mobility graphene encapsulated in tungsten disulfide

High mobility is a crucial requirement for a large variety of electronic device applications. The state of the art for high-quality graphene devices is based on heterostructures made with graphene encapsulated in >40 nm-thick flakes of hexagonal boron nitride (hBN). Unfortunately, scaling up multilayer hBN while precisely controlling the number of layers remains an outstanding challenge, resulting in a rough material unable to enhance the mobility of graphene. This leads to the pursuit of alternative, scalable materials, which can be used as substrates and encapsulants for graphene. Tungsten disulfide (WS2) is a transition metal dichalcogenide, which was grown in large (∼mm-size) multi-layers by chemical vapor deposition. However, the resistance vs gate voltage characteristics when gating graphene through WS2 exhibit largely hysteretic shifts of the charge neutrality point on the order of Δn∼ 3 × 1011 cm−2, hindering the use of WS2 as a reliable encapsulant. The hysteresis originates due to the charge traps from sulfur vacancies present in WS2. In this work, we report the use of WS2 as a substrate and overcome the hysteresis issues by chemically treating WS2 with a super-acid, which passivates these vacancies and strips the surface from contaminants. The hysteresis is significantly reduced by about two orders of magnitude, down to values as low as Δn∼ 2 × 109 cm−2, while the room-temperature mobility of WS2-encapsulated graphene is as high as ∼62 × 103 cm2 V−1 s−1 at a carrier density of n ∼ 1 ×1012 cm−2. Our results promote WS2 as a valid alternative to hBN as an encapsulant for high-performance graphene devices.

Physics↗