Abstract: Ripples in pristine freestanding graphene naturally orient themselves in an array that is alternately curved-up and curved-down; maintaining an average height of zero. Using scanning tunneling microscopy (STM) to apply a local force, the graphene sheet will reversibly rise and fall in height until the height reaches 60%-70% of its maximum at which point a sudden, permanent jump occurs. We successfully model the ripples as a spin-half Ising magnetic system, where the height of the graphene plays the role of the spin. The permanent jump in height, controlled by the tunneling current, is found to be equivalent to an antiferromagnetic-to-ferromagnetic phase transition. The thermal load underneath the STM tip alters the local tension and is identified as the responsible mechanism for the phase transition. Four universal critical exponents are measured from our STM data, and the model provides insight into the statistical role of graphene's unusual negative thermal expansion coefficient.

Abstract: The appearance of massless Dirac fermions in graphene requires two equivalent carbon sublattices of trigonal shape. While the generation of an effective mass and a band gap at the Dirac point remains an unresolved problem for freestanding extended graphene, it is well established by breaking translational symmetry by confinement and by breaking sublattice symmetry by interaction with a substrate. One of the strongest sublattice-symmetry-breaking interactions with predicted and measured band gaps ranging from 400 meV to more than 3 eV has been attributed to the interfaces of graphene with Ni and Co, which are also promising spin-filter interfaces. Here, we apply angle-resolved photoemission to epitaxial graphene on Ni (111) and Co(0001) to show the presence of intact Dirac cones 2.8 eV below the Fermi level. Our results challenge the common belief that the breaking of sublattice symmetry by a substrate and the opening of the band gap at the Dirac energy are in a straightforward relation. A simple effective model of a biased bilayer structure composed of graphene and a sublattice-symmetry-broken layer, corroborated by density-functional-theory calculations, demonstrates the general validity of our conclusions.

Abstract: We compare three first-principles methods of calculating the Curie temperature in two-dimensional (2D) ferromagnetic materials (FM), modeled using the Heisenberg model, and propose a simple formula for estimating the Curie temperature with high accuracy that works for all common 2D lattice types. First, we study the effect of exchange anisotropy on the Curie temperature calculated using the Monte Carlo (MC), the Green's function, and the renormalized spin-wave (RNSW) methods. We find that the Green's function method overestimates the Curie temperature in high-anisotropy regimes compared to the MC method, whereas the RNSW method underestimates the Curie temperature compared to the MC and the Green's function methods. Next, we propose a closed-form formula for calculating the Curie temperature of 2D FMs, which provides an estimate of the Curie temperature that is greatly improved over the mean-field expression for magnetic material screening. We apply the closed-form formula to predict the Curie temperature 2D magnets screened from the C2DB database and discover several high Curie temperature FMs, with Fe2F2 and MoI2 emerging as the most promising 2D ferromagnets. Finally, by comparing to experimental results for CrI3, CrCl3, and CrBr3, we conclude that for small effective anisotropies, the Green's-function-based equations are preferable, while for larger anisotropies, MC-based results are more predictive.

Abstract: The in-plane dc conductivity of twisted bilayer graphene is calculated using an expansion of the real-space Kubo-Bastin conductivity in terms of Chebyshev polynomials. We investigate within a tight-binding approach the transport properties as a function of rotation angle, applied perpendicular electric field, and vacancy disorder. We find that for high-angle twists, the two layers are effectively decoupled, and the minimum conductivity at the Dirac point corresponds to double the value observed in monolayer graphene. This remains valid even in the presence of vacancies, hinting that chiral symmetry is still preserved. On the contrary, for low twist angles, the conductivity at the Dirac point depends on the twist angle and is not protected in the presence of disorder. Furthermore, for low angles and in the presence of an applied electric field, we find that the chiral boundary states emerging between AB and BA regions contribute to the dc conductivity, despite the appearance of localized states in the AA regions. The results agree qualitatively with recent transport experiments in low-angle twisted bilayer graphene.

Abstract: We investigate the diffusive electron-transport properties of charge-doped graphene ribbons and nanoribbons with imperfect edges. We consider different regimes of edge scattering, ranging from wide graphene ribbons with (partially) diffusive edge scattering to ribbons with large width variations and nanoribbons with atomistic edge roughness. For the latter, we introduce an approach based on pseudopotentials, allowing for an atomistic treatment of the band structure and the scattering potential, on the self-consistent solution of the Boltzmann transport equation within the relaxation-time approximation and taking into account the edge-roughness properties and statistics. The resulting resistivity depends strongly on the ribbon orientation, with zigzag (armchair) ribbons showing the smallest (largest) resistivity and intermediate ribbon orientations exhibiting intermediate resistivity values. The results also show clear resistivity peaks, corresponding to peaks in the density of states due to the confinement-induced subband quantization, except for armchair-edge ribbons that show a very strong width dependence because of their claromatic behavior. Furthermore, we identify a strong interplay between the relative position of the two valleys of graphene along the transport direction, the correlation profile of the atomistic edge roughness, and the chiral valley modes, leading to a peculiar strongly suppressed resistivity regime, most pronounced for the zigzag orientation.

Abstract: Motivated by the potential of atomically thin magnets towards achieving tunable high-frequency magnonics, we detail the spin-wave dispersion of bilayer CrI3. We demonstrate that the magnonic behavior of the bilayer strongly depends on its stacking configuration and the interlayer magnetic ordering, where a topological band gap opens in the dispersion caused by the Dzyaloshinskii-Moriya and Kitaev interactions, classifying bilayer CrI3 as a topological magnon insulator. We further reveal that both the size and the topology of the band gap in a CrI3 bilayer with an antiferromagnetic interlayer ordering are tunable by an external magnetic field.

Abstract: The wide measured range of thermal conductivities (k) for monolayer MoS2 and the corresponding incongruent calculated values in the literature all suggest that extrinsic defect thermal resistance is significant and varied in synthesized samples of this material. Here we present defect-mediated thermal transport calculations of MoS2 using interatomic forces derived from density functional theory combined with Green's function methods to describe phonon-point-defect interactions and a Peierls-Boltzmann formalism for transport. Conductivity calculations for bulk and monolayer MoS2 using different density functional formalisms are compared. Nonperturbative first-principles methods are used to describe defect-mediated spectral functions, scattering rates, and phonon k, particularly from sulfur vacancies (VS), and in the context of the plethora of measured and calculated literature values. We find that k of monolayer MoS2 is sensitive to phonon-VS scattering in the range of experimentally observed densities, and that first-principles k calculations using these densities can explain the range of measured values found in the literature. Furthermore, measured k values for bulk MoS2 are more consistent because VS defects are not as prevalent.

Abstract: Two-dimensional (2D) honeycomb ferromagnets, such as monolayer chromium trihalides, are predicted to behave as topological magnon insulators, characterized by an insulating bulk and topologically protected edge states, giving rise to a thermal magnon Hall effect. Here we report the behavior of the topological magnons in monolayer CrI3 under external stimuli, including biaxial and uniaxial strain, electric gating, as well as in-plane and out-of-plane magnetic field, revealing that one can thereby tailor the magnetic states as well as the size and the topology of the magnonic bandgap. These findings broaden the perspective of using 2D magnetic materials to design topological magnonic devices.