Showing 1–50 of 80 results for author: Luskin, M

Tunable atomically enhanced moiré Berry curvatures in twisted triple bilayer graphene

Authors: Konstantin Davydov, Ziyan Zhu, Noah Friedman, Ethan Gramowski, Yaotian Li, Jack Tavakley, Kenji Watanabe, Takashi Taniguchi, Mitchell Luskin, Efthimios Kaxiras, Ke Wang

Abstract: We report a twisted triple bilayer graphene platform consisting of three units of Bernal bilayer graphene (BLG) consecutively twisted at 1.49° and 1.68°. We observe inter-moiré Hofstadter butterflies from two co-existing moiré superlattices and a Hofstadter butterfly from reconstructed moiré-of-moiré lattice, and show that their Brown-Zak (BZ) oscillations quantitatively agree with each other, bot… ▽ More We report a twisted triple bilayer graphene platform consisting of three units of Bernal bilayer graphene (BLG) consecutively twisted at 1.49° and 1.68°. We observe inter-moiré Hofstadter butterflies from two co-existing moiré superlattices and a Hofstadter butterfly from reconstructed moiré-of-moiré lattice, and show that their Brown-Zak (BZ) oscillations quantitatively agree with each other, both evidencing strong atomic reconstruction with a lattice constant of 18.1 nm. We further demonstrate such atomic reconstruction strongly enhances the Berry curvature of each moiré and moiré-of-moiré band-insulator state, characterized by measured strong non-local valley Hall effect (VHE) that sensitively depends on the inter-moiré competition strength, tunable by manipulating the out-of-the-plane carrier distribution which controls the magnitude of the valley currents. Our study sheds new light on the microscopic mechanism of atomic and electronic reconstruction in twisted-multilayer systems, by investigating novel emergent quantum phenomena of reconstructed quasi-crystalline moiré-of-moiré superlattice, including a new type of moiré-of-moiré band-insulator states and atomically enhanced moiré Berry curvature. We show that the reconstructed electronic band can be versatilely tuned by electrostatics, providing an approach towards engineering the band structure and its topology for a novel quantum material platform with designer electrical and optical properties. △ Less

Submitted 11 August, 2024; originally announced August 2024.

Learning the local density of states of a bilayer moiré material in one dimension

Authors: Diyi Liu, Alexander B. Watson, Stephen Carr, Mitchell Luskin

Abstract: Recent work of three of the authors showed that the operator which maps the local density of states of a one-dimensional untwisted bilayer material to the local density of states of the same bilayer material at non-zero twist, known as the twist operator, can be learned by a neural network. In this work, we first provide a mathematical formulation of that work, making the relevant models and opera… ▽ More Recent work of three of the authors showed that the operator which maps the local density of states of a one-dimensional untwisted bilayer material to the local density of states of the same bilayer material at non-zero twist, known as the twist operator, can be learned by a neural network. In this work, we first provide a mathematical formulation of that work, making the relevant models and operator learning problem precise. We then prove that the operator learning problem is well-posed for a family of one-dimensional models. To do this, we first prove existence and regularity of the twist operator by solving an inverse problem. We then invoke the universal approximation theorem for operators to prove existence of a neural network capable of approximating the twist operator. △ Less

Submitted 6 May, 2024; originally announced May 2024.

Comments: 28 pages, 4 figures

MSC Class: 81; 65; 35

Electron Collimation in Twisted Bilayer Graphene via Gate-defined Moiré Barriers

Authors: Wei Ren, Xi Zhang, Ziyan Zhu, Moosa Khan, Kenji Watanabe, Takashi Taniguchi, Efthimios Kaxiras, Mitchell Luskin, Ke Wang

Abstract: Electron collimation via a graphene pn-junction allows electrostatic control of ballistic electron trajectories akin to that of an optical circuit. Similar manipulation of novel correlated electronic phases in twisted-bilayer graphene (tBLG) can provide additional probes to the underlying physics and device components towards advanced quantum electronics. In this work, we demonstrate collimation o… ▽ More Electron collimation via a graphene pn-junction allows electrostatic control of ballistic electron trajectories akin to that of an optical circuit. Similar manipulation of novel correlated electronic phases in twisted-bilayer graphene (tBLG) can provide additional probes to the underlying physics and device components towards advanced quantum electronics. In this work, we demonstrate collimation of the electron flow via gate-defined moiré barriers in a tBLG device, utilizing the band-insulator gap of the moiré superlattice. A single junction can be tuned to host a chosen combination of conventional pseudo barrier and moiré tunnel barriers, from which we demonstrate improved collimation efficiency. By measuring transport through two consecutive moiré collimators separated by 1 um, we demonstrate evidence of electron collimation in tBLG in the presence of realistic twist-angle inhomogeneity. △ Less

Submitted 30 March, 2024; originally announced April 2024.

Tunable Inter-Moiré Physics in Consecutively-Twisted Trilayer Graphene

Authors: Wei Ren, Konstantin Davydov, Ziyan Zhu, Jaden Ma, Kenji Watanabe, Takashi Taniguchi, Efthimios Kaxiras, Mitchell Luskin, Ke Wang

Abstract: We fabricate a twisted trilayer graphene device with consecutive twist angles of 1.33 and 1.64 degrees, in which we electrostatically tune the electronic states from each of the two co-existing moiré superlattices and the interactions between them. When both moiré superlattices contribute equally to electrical transport, we report a new type of inter-moiré Hofstadter butterfly. Its Brown-Zak oscil… ▽ More We fabricate a twisted trilayer graphene device with consecutive twist angles of 1.33 and 1.64 degrees, in which we electrostatically tune the electronic states from each of the two co-existing moiré superlattices and the interactions between them. When both moiré superlattices contribute equally to electrical transport, we report a new type of inter-moiré Hofstadter butterfly. Its Brown-Zak oscillation corresponds to one of the intermediate quasicrystal length scales of the reconstructed moiré of moiré (MoM) superlattice, shedding new light on emergent physics from competing atomic orders. △ Less

Submitted 16 November, 2023; originally announced November 2023.

Modeling of electronic dynamics in twisted bilayer graphene

Authors: Tianyu Kong, Diyi Liu, Mitchell Luskin, Alexander B. Watson

Abstract: We consider the problem of numerically computing the quantum dynamics of an electron in twisted bilayer graphene. The challenge is that atomic-scale models of the dynamics are aperiodic for generic twist angles because of the incommensurability of the layers. The Bistritzer-MacDonald PDE model, which is periodic with respect to the bilayer's moiré pattern, has recently been shown to rigorously des… ▽ More We consider the problem of numerically computing the quantum dynamics of an electron in twisted bilayer graphene. The challenge is that atomic-scale models of the dynamics are aperiodic for generic twist angles because of the incommensurability of the layers. The Bistritzer-MacDonald PDE model, which is periodic with respect to the bilayer's moiré pattern, has recently been shown to rigorously describe these dynamics in a parameter regime. In this work, we first prove that the dynamics of the tight-binding model of incommensurate twisted bilayer graphene can be approximated by computations on finite domains. The main ingredient of this proof is a speed of propagation estimate proved using Combes-Thomas estimates. We then provide extensive numerical computations which clarify the range of validity of the Bistritzer-MacDonald model. △ Less

Submitted 25 February, 2024; v1 submitted 20 August, 2023; originally announced August 2023.

Comments: 27 pages, 8 figures

MSC Class: 81; 65; 35

From incommensurate bilayer heterostructures to Allen-Cahn: An exact thermodynamic limit

Authors: Michael Hott, Alexander B. Watson, Mitchell Luskin

Abstract: We give a complete and rigorous derivation of the mechanical energy for twisted 2D bilayer heterostructures without any approximation beyond the existence of an empirical many-body site energy. Our results apply to both the continuous and discontinuous continuum limit. Approximating the intralayer Cauchy-Born energy by linear elasticity theory and assuming an interlayer coupling via pair potential… ▽ More We give a complete and rigorous derivation of the mechanical energy for twisted 2D bilayer heterostructures without any approximation beyond the existence of an empirical many-body site energy. Our results apply to both the continuous and discontinuous continuum limit. Approximating the intralayer Cauchy-Born energy by linear elasticity theory and assuming an interlayer coupling via pair potentials, our model reduces to a modified Allen-Cahn functional. We rigorously control the error, and, in the case of sufficiently smooth lattice displacements, provide a rate of convergence for twist angles satisfying a Diophantine condition. △ Less

Submitted 27 June, 2024; v1 submitted 29 May, 2023; originally announced May 2023.

Comments: 67 pages, 9 figures, V4: Streamlined presentation

MSC Class: 82-10; 74-10; 70-10; 37A44; 37A46; 35Qxx; 35Q56

Mathematical aspects of the Kubo formula for electrical conductivity with dissipation

Authors: Alexander B. Watson, Dionisios Margetis, Mitchell Luskin

Abstract: In this expository article, we present a systematic formal derivation of the Kubo formula for the linear-response current due to a time-harmonic electric field applied to non-interacting, spinless charged particles in a finite volume in the quantum setting. We model dissipation in a transparent way by assuming a sequence of scattering events occurring at random-time intervals modeled by a Poisson… ▽ More In this expository article, we present a systematic formal derivation of the Kubo formula for the linear-response current due to a time-harmonic electric field applied to non-interacting, spinless charged particles in a finite volume in the quantum setting. We model dissipation in a transparent way by assuming a sequence of scattering events occurring at random-time intervals modeled by a Poisson distribution. By taking the large-volume limit, we derive special cases of the formula for free electrons, continuum and tight-binding periodic systems, and the nearest-neighbor tight-binding model of graphene. We present the analogous formalism with dissipation to derive the Drude conductivity of classical free particles. △ Less

Submitted 13 August, 2023; v1 submitted 9 April, 2023; originally announced April 2023.

Comments: 35 pages. Small corrections and improvements to exposition

MSC Class: 81; 78

Relaxation and domain wall structure of bilayer moire systems

Authors: Paul Cazeaux, Drake Clark, Rebecca Engelke, Philip Kim, Mitchell Luskin

Abstract: Moire patterns result from setting a 2D material such as graphene on another 2D material with a small twist angle or from the lattice mismatch of 2D heterostructures. We present a continuum model for the elastic energy of these bilayer moire structures that includes an intralayer elastic energy and an interlayer misfit energy that is minimized at two stackings (disregistries). We show by theory an… ▽ More Moire patterns result from setting a 2D material such as graphene on another 2D material with a small twist angle or from the lattice mismatch of 2D heterostructures. We present a continuum model for the elastic energy of these bilayer moire structures that includes an intralayer elastic energy and an interlayer misfit energy that is minimized at two stackings (disregistries). We show by theory and computation that the displacement field that minimizes the global elastic energy subject to a global boundary constraint gives large alternating regions of one of the two energy-minimizing stackings separated by domain walls. We derive a model for the domain wall structure from the continuum bilayer energy and give a rigorous asymptotic estimate for the structure. We also give an improved estimate for the L2-norm of the gradient on the moire unit cell for twisted bilayers that scales at most inversely linearly with the twist angle, a result which is consistent with the formation of one-dimensional domain walls with a fixed width around triangular domains at very small twist angles. △ Less

Submitted 26 March, 2023; v1 submitted 22 November, 2022; originally announced November 2022.

Comments: 21 pages, 14 figures

MSC Class: 74B99; 74G65; 74K99

On the Su-Schrieffer-Heeger model of electron transport: low-temperature optical conductivity by the Mellin transform

Authors: Dionisios Margetis, Alexander B. Watson, Mitchell Luskin

Abstract: We describe the low-temperature optical conductivity as a function of frequency for a quantum-mechanical system of electrons that hop along a polymer chain. To this end, we invoke the Su-Schrieffer-Heeger \emph{tight-binding} Hamiltonian for non-interacting spinless electrons on a one-dimensional (1D) lattice. Our goal is to show via asymptotics how the interband conductivity of this system behave… ▽ More We describe the low-temperature optical conductivity as a function of frequency for a quantum-mechanical system of electrons that hop along a polymer chain. To this end, we invoke the Su-Schrieffer-Heeger \emph{tight-binding} Hamiltonian for non-interacting spinless electrons on a one-dimensional (1D) lattice. Our goal is to show via asymptotics how the interband conductivity of this system behaves as the smallest energy bandgap tends to close. Our analytical approach includes: (i) the Kubo-type formulation for the optical conductivity with a nonzero damping due to microscopic collisions; (ii) reduction of this formulation to a 1D momentum integral over the Brillouin zone; and (iii) evaluation of this integral in terms of elementary functions via the three-dimensional Mellin transform with respect to key physical parameters and subsequent inversion in a region of the respective complex space. Our approach reveals an intimate connection of the behavior of the conductivity to particular singularities of its Mellin transform. The analytical results are found in good agreement with direct numerical computations. △ Less

Submitted 26 September, 2022; originally announced September 2022.

Comments: 35 pages, 10 figures

Bistritzer-MacDonald dynamics in twisted bilayer graphene

Authors: Alexander B. Watson, Tianyu Kong, Allan H. MacDonald, Mitchell Luskin

Abstract: The Bistritzer-MacDonald (BM) model, introduced in \cite{Bistritzer2011}, attempts to capture the electronic properties of twisted bilayer graphene (TBG), even at incommensurate twist angles, by an effective periodic model over the bilayer moiré pattern. Starting from a tight-binding model, we identify a regime where the BM model emerges as the effective dynamics for electrons modeled as wave-pack… ▽ More The Bistritzer-MacDonald (BM) model, introduced in \cite{Bistritzer2011}, attempts to capture the electronic properties of twisted bilayer graphene (TBG), even at incommensurate twist angles, by an effective periodic model over the bilayer moiré pattern. Starting from a tight-binding model, we identify a regime where the BM model emerges as the effective dynamics for electrons modeled as wave-packets spectrally concentrated at the monolayer Dirac points, up to error that can be rigorously estimated. Using measured values of relevant physical constants, we argue that this regime is realized in TBG at the first "magic" angle. △ Less

Submitted 13 October, 2022; v1 submitted 27 July, 2022; originally announced July 2022.

Comments: 50 pages, 9 figures. Fixed wave-packet scaling, bounded some higher-order terms in the residual, and added reference to Bal-Cazeaux-Massatt-Quinn. Bounding the higher-order terms required strengthening our assumption on the regularity of initial data

Seeing moiré: convolutional network learning applied to twistronics

Authors: Diyi Liu, Mitchell Luskin, Stephen Carr

Abstract: Moiré patterns made of two-dimensional (2D) materials represent highly tunable electronic Hamiltonians, allowing a wide range of quantum phases to emerge in a single material. Current modeling techniques for moiré electrons requires significant technical work specific to each material, impeding large-scale searches for useful moiré materials. In order to address this difficulty, we have developed… ▽ More Moiré patterns made of two-dimensional (2D) materials represent highly tunable electronic Hamiltonians, allowing a wide range of quantum phases to emerge in a single material. Current modeling techniques for moiré electrons requires significant technical work specific to each material, impeding large-scale searches for useful moiré materials. In order to address this difficulty, we have developed a material-agnostic machine learning approach and test it here on prototypical one-dimensional (1D) moiré tight-binding models. We utilize the stacking dependence of the local density of states (SD-LDOS) to convert information about electronic bandstructure into physically relevant images. We then train a neural network that successfully predicts moiré electronic structure from the easily computed SD-LDOS of aligned bilayers. This network can satisfactorily predict moiré electronic structures, even for materials that are not included in its training data. △ Less

Submitted 22 July, 2022; originally announced July 2022.

Comments: 11 pages, 12 figures

Non-Abelian topological defects and strain mapping in 2D moiré materials

Authors: Rebecca Engelke, Hyobin Yoo, Stephen Carr, Kevin Xu, Paul Cazeaux, Richard Allen, Andres Mier Valdivia, Mitchell Luskin, Efthimios Kaxiras, Minhyong Kim, Jung Hoon Han, Philip Kim

Abstract: We present a general method to analyze the topological nature of the domain boundary connectivity that appeared in relaxed moiré superlattice patterns at the interface of 2-dimensional (2D) van der Waals (vdW) materials. At large enough moiré lengths, all moiré systems relax into commensurated 2D domains separated by networks of dislocation lines. The nodes of the 2D dislocation line network can b… ▽ More We present a general method to analyze the topological nature of the domain boundary connectivity that appeared in relaxed moiré superlattice patterns at the interface of 2-dimensional (2D) van der Waals (vdW) materials. At large enough moiré lengths, all moiré systems relax into commensurated 2D domains separated by networks of dislocation lines. The nodes of the 2D dislocation line network can be considered as vortex-like topological defects. We find that a simple analogy to common topological systems with an $S^1$ order parameter, such as a superconductor or planar ferromagnet, cannot correctly capture the topological nature of these defects. For example, in twisted bilayer graphene, the order parameter space for the relaxed moiré system is homotopy equivalent to a punctured torus. Here, the nodes of the 2D dislocation network can be characterized as elements of the fundamental group of the punctured torus, the free group on two generators, endowing these network nodes with non-Abelian properties. Extending this analysis to consider moiré patterns generated from any relative strain, we find that antivortices occur in the presence of anisotropic heterostrain, such as shear or anisotropic expansion, while arrays of vortices appear under twist or isotropic expansion between vdW materials. Experimentally, utilizing the dark field imaging capability of transmission electron microscopy (TEM), we demonstrate the existence of vortex and antivortex pair formation in a moiré system, caused by competition between different types of heterostrains in the vdW interfaces. We also present a methodology for mapping the underlying heterostrain of a moiré structure from experimental TEM data, which provides a quantitative relation between the various components of heterostrain and vortex-antivortex density in moiré systems. △ Less

Submitted 16 July, 2022; v1 submitted 11 July, 2022; originally announced July 2022.

Comments: 15 pages with 11 figures

Electronic Observables for Relaxed Bilayer 2D Heterostructures in Momentum Space

Authors: Daniel Massatt, Stephen Carr, Mitchell Luskin

Abstract: Momentum space transformations for incommensurate 2D electronic structure calculations are fundamental for reducing computational cost and for representing the data in a more physically motivating format, as exemplified in the Bistritzer-MacDonald model. However, these transformations can be difficult to implement in more complex systems such as when mechanical relaxation patterns are present. In… ▽ More Momentum space transformations for incommensurate 2D electronic structure calculations are fundamental for reducing computational cost and for representing the data in a more physically motivating format, as exemplified in the Bistritzer-MacDonald model. However, these transformations can be difficult to implement in more complex systems such as when mechanical relaxation patterns are present. In this work, we aim for two objectives. Firstly, we strive to simplify the understanding and implementation of this transformation by rigorously writing the transformations between the four relevant spaces, which we denote real space, configuration space, momentum space, and reciprocal space. This provides a straight-forward algorithm for writing the complex momentum space model from the original real space model. Secondly, we implement this for twisted bilayer graphene with mechanical relaxation affects included. We also analyze the convergence rates of the approximations, and show the tight-binding coupling range increases for smaller relative twists between layers, demonstrating that the 3-nearest neighbor coupling of the Bistritzer-MacDonald model is insufficient when mechanical relaxation is included for very small angles. We quantify this and verify with numerical simulation. △ Less

Submitted 31 August, 2023; v1 submitted 30 September, 2021; originally announced September 2021.

Comments: 40 pages, 10 figures

MSC Class: 65 (primary); 81-10 (secondary)

Gate-tunable Veselago Interference in a Bipolar Graphene Microcavity

Authors: Xi Zhang*, Wei Ren*, Elliot Bell, Ziyan Zhu, Kenji Watanabe, Takashi Taniguchi, Efthimios Kaxiras, Mitchell Luskin, Ke Wang

Abstract: The relativistic charge carriers in monolayer graphene can be manipulated in manners akin to conventional optics (electron-optics): angle-dependent Klein tunneling collimates an electron beam (analogous to a laser), while a Veselago refraction process focuses it (analogous to an optical lens). Both processes have been previously investigated, but the collimation and focusing efficiency have been r… ▽ More The relativistic charge carriers in monolayer graphene can be manipulated in manners akin to conventional optics (electron-optics): angle-dependent Klein tunneling collimates an electron beam (analogous to a laser), while a Veselago refraction process focuses it (analogous to an optical lens). Both processes have been previously investigated, but the collimation and focusing efficiency have been reported to be relatively low even in state-of-the-art ballistic pn-junction devices. These limitations prevented the realization of more advanced quantum devices based on electron-optical interference, while understanding of the underlying physics remains elusive. Here, we present a novel device architecture of a graphene microcavity defined by carefully-engineered local strain and electrostatic fields. We create a controlled electron-optic interference process at zero magnetic field as a consequence of consecutive Veselago refractions in the microcavity and provide direct experimental evidence through low-temperature electrical transport measurements. The experimentally observed first-, second-, and third-order interference peaks agree quantitatively with the Veselago physics in a microcavity. In addition, we demonstrate decoherence of the interference by an external magnetic field, as the cyclotron radius becomes comparable to the interference length scale. For its application in electron-optics, we utilize Veselago interference to localize uncollimated electrons and characterize its contribution in further improving collimation efficiency. Our work sheds new light on relativistic single-particle physics and provides important technical improvements toward next-generation quantum devices based on the coherent manipulation of electron momentum and trajectory. △ Less

Submitted 17 June, 2021; originally announced June 2021.

Existence of the first magic angle for the chiral model of bilayer graphene

Authors: Alexander B. Watson, Mitchell Luskin

Abstract: We consider the chiral model of twisted bilayer graphene introduced by Tarnopolsky-Kruchkov-Vishwanath (TKV). TKV have proved that for inverse twist angles $α$ such that the effective Fermi velocity at the moiré $K$ point vanishes, the chiral model has a perfectly flat band at zero energy over the whole Brillouin zone. By a formal expansion, TKV found that the Fermi velocity vanishes at… ▽ More We consider the chiral model of twisted bilayer graphene introduced by Tarnopolsky-Kruchkov-Vishwanath (TKV). TKV have proved that for inverse twist angles $α$ such that the effective Fermi velocity at the moiré $K$ point vanishes, the chiral model has a perfectly flat band at zero energy over the whole Brillouin zone. By a formal expansion, TKV found that the Fermi velocity vanishes at $α\approx .586$. In this work, we give a proof that the Fermi velocity vanishes for at least one $α$ between $.57$ and $.61$ by rigorously justifying TKV's formal expansion of the Fermi velocity over a sufficiently large interval of $α$ values. The idea of the proof is to project the TKV Hamiltonian onto a finite dimensional subspace, and then expand the Fermi velocity in terms of explicitly computable linear combinations of modes in the subspace, while controlling the error. The proof relies on two propositions whose proofs are computer-assisted, i.e., numerical computation together with worst-case estimates on the accumulation of round-off error which show that round-off error cannot possibly change the conclusion of the computation. The propositions give a bound below on the spectral gap of the projected Hamiltonian, an Hermitian $80 \times 80$ matrix whose spectrum is symmetric about $0,$ and verify that two real 18th order polynomials, which approximate the numerator of the Fermi velocity, take values with definite sign when evaluated at specific values of $α$. Together with TKV's work our result proves existence of at least one perfectly flat band of the chiral model. △ Less

Submitted 10 August, 2021; v1 submitted 13 April, 2021; originally announced April 2021.

Comments: 61 pages (last 7 pages are supplementary material), 11 figures. Improved after review comments: (1) proved that the Fermi velocity perturbation series converges for $|α| < 1/3$ (2) added proof that round-off error cannot affect the conclusions of our numerical computations (3) more extensive discussion of related work

MSC Class: 35Q; 35P; 81V

Journal ref: J. Math. Phys. 62 9 (2021) 091502

Dipole excitation of collective modes in viscous two-dimensional electron systems

Authors: Vera Andreeva, Denis A. Bandurin, Mitchell Luskin, Dionisios Margetis

Abstract: We describe the structure of the time-harmonic electromagnetic field of a vertical Hertzian electric dipole source radiating over an infinite, translation invariant two-dimensional electron system. Our model for the electron flow takes into account the effects of shear and Hall viscosities as well as an external static magnetic field perpendicular to the sheet. We identify two wave modes, namely,… ▽ More We describe the structure of the time-harmonic electromagnetic field of a vertical Hertzian electric dipole source radiating over an infinite, translation invariant two-dimensional electron system. Our model for the electron flow takes into account the effects of shear and Hall viscosities as well as an external static magnetic field perpendicular to the sheet. We identify two wave modes, namely, a surface plasmon and a diffusive mode. In the presence of an external static magnetic field, the diffusive mode combines the features of both the conventional and Hall diffusion and may exhibit a negative group velocity. In our analysis, we solve exactly a boundary value problem for the time-harmonic Maxwell equations coupled with linearized hydrodynamic equations for the flat, two-dimensional material. By numerically evaluating the integrals for the electromagnetic field on the sheet, we find that the plasmon contribution dominates in the intermediate-field region of the dipole source. In contrast, the amplitude of the diffusive mode reaches its maximum value in the near-field region, and quickly decays with the distance from the source. We demonstrate that the diffusive mode can be distinguished from the plasmon in the presence of the static magnetic field, when the highly oscillatory plasmon is gapped and tends to disappear. △ Less

Submitted 2 August, 2020; originally announced August 2020.

Comments: 18 pages, 7 figures

Journal ref: Phys. Rev. B 102, 205411 (2020)

Twisted Trilayer Graphene: a Precisely Tunable Platform for Correlated Electrons

Authors: Ziyan Zhu, Stephen Carr, Daniel Massatt, Mitchell Luskin, Efthimios Kaxiras

Abstract: We introduce twisted trilayer graphene (tTLG) with two independent twist angles as an ideal system for the precise tuning of the electronic interlayer coupling to maximize the effect of correlated behaviors. As established by experiment and theory in the related twisted bilayer graphene system, van Hove singularities (VHS) in the density of states can be used as a proxy of the tendency for correla… ▽ More We introduce twisted trilayer graphene (tTLG) with two independent twist angles as an ideal system for the precise tuning of the electronic interlayer coupling to maximize the effect of correlated behaviors. As established by experiment and theory in the related twisted bilayer graphene system, van Hove singularities (VHS) in the density of states can be used as a proxy of the tendency for correlated behaviors. To explore the evolution of VHS in the twist-angle phase space of tTLG, we present a general low-energy electronic structure model for any pair of twist angles. We show that the basis of the model has infinite dimensions even at a finite energy cutoff and that no Brillouin zone exists even in the continuum limit. Using this model, we demonstrate that the tTLG system exhibits a wide range of magic angles at which VHS merge and the density of states has a sharp peak at the charge-neutrality point through two distinct mechanisms: the incommensurate perturbation of twisted bilayer graphene's flat bands or the equal hybridization between two bilayer moiré superlattices. △ Less

Submitted 8 October, 2020; v1 submitted 30 May, 2020; originally announced June 2020.

Journal ref: Phys. Rev. Lett. 125, 116404 (2020)

Finite-size effects in wave transmission through plasmonic crystals: A tale of two scales

Authors: Matthias Maier, Mitchell Luskin, Dionisios Margetis

Abstract: We study optical coefficients that characterize wave propagation through layered structures called plasmonic crystals. These consist of a finite number of stacked metallic sheets embedded in dielectric hosts with a subwavelength spacing. By adjustment of the frequency, spacing, number as well as geometry of the layers, these structures may exhibit appealing transmission properties in a range of fr… ▽ More We study optical coefficients that characterize wave propagation through layered structures called plasmonic crystals. These consist of a finite number of stacked metallic sheets embedded in dielectric hosts with a subwavelength spacing. By adjustment of the frequency, spacing, number as well as geometry of the layers, these structures may exhibit appealing transmission properties in a range of frequencies from the terahertz to the mid-infrared regime. Our approach uses a blend of analytical and numerical methods for the distinct geometries with infinite, translation invariant, flat sheets and nanoribbons. We describe the transmission of plane waves through a plasmonic crystal in comparison to an effective dielectric slab of equal total thickness that emerges from homogenization, in the limit of zero interlayer spacing. We demonstrate numerically that the replacement of the discrete plasmonic crystal by its homogenized counterpart can accurately capture a transmission coefficient akin to the extinction spectrum, even for a relatively small number of layers. We point out the role of a geometry-dependent corrector field, which expresses the effect of subwavelength surface plasmons. In particular, by use of the corrector we describe lateral resonances inherent to the nanoribbon geometry. △ Less

Submitted 5 August, 2020; v1 submitted 26 May, 2020; originally announced May 2020.

Journal ref: Phys. Rev. B 102, 075308 (2020)

Homogenization of hydrodynamic transport in Dirac fluids

Authors: Guillaume Bal, Andrew Lucas, Mitchell Luskin

Abstract: Large-scale electrical and thermal currents in ordinary metals are well approximated by effective medium theory: global transport properties are governed by the solution to homogenized coupled diffusion equations. In some metals, including the Dirac fluid of nearly charge neutral graphene, microscopic transport is not governed by diffusion, but by a more complicated set of linearized hydrodynamic… ▽ More Large-scale electrical and thermal currents in ordinary metals are well approximated by effective medium theory: global transport properties are governed by the solution to homogenized coupled diffusion equations. In some metals, including the Dirac fluid of nearly charge neutral graphene, microscopic transport is not governed by diffusion, but by a more complicated set of linearized hydrodynamic equations, which form a system of degenerate elliptic equations coupled with the Stokes equation for fluid velocity. In sufficiently inhomogeneous media, these hydrodynamic equations reduce to homogenized diffusion equations. We re-cast the hydrodynamic transport equations as the infimum of a functional over conserved currents, and present a functional framework to model and compute the homogenized diffusion tensor relating electrical and thermal currents to charge and temperature gradients. We generalize to this system two well-known results in homogenization theory: Tartar's proof of local convergence to the homogenized theory in periodic and highly oscillatory media, and sub-additivity of the above functional in random media with highly oscillatory, stationary and ergodic coefficients. △ Less

Submitted 22 April, 2020; originally announced April 2020.

Comments: 27 pages

Nonlinear eigenvalue problems for coupled Helmholtz equations modeling gradient-index graphene waveguides

Authors: Jung Heon Song, Matthias Maier, Mitchell Luskin

Abstract: We discuss a quartic eigenvalue problem arising in the context of an optical waveguiding problem involving atomically thick 2D materials. The waveguide configuration we consider consists of a gradient-index (spatially dependent) dielectric equipped with conducting interior interfaces. This leads to a quartic eigenvalue problem with mixed transverse electric and transverse magnetic modes, and stron… ▽ More We discuss a quartic eigenvalue problem arising in the context of an optical waveguiding problem involving atomically thick 2D materials. The waveguide configuration we consider consists of a gradient-index (spatially dependent) dielectric equipped with conducting interior interfaces. This leads to a quartic eigenvalue problem with mixed transverse electric and transverse magnetic modes, and strongly coupled electric and magnetic fields. We derive a weak formulation of the quartic eigenvalue problem and introduce a numerical solver based on a quadratification approach in which the quartic eigenvalue problem is transformed to a spectrally equivalent companion problem. We verify our numerical framework against analytical solutions for prototypical geometries. As a practical example, we demonstrate how an improved quality factor (defined by the ratio of the real and the imaginary part of the computed eigenvalues) can be obtained for a family of gradient-index host materials with internal conducting interfaces. We outline how this result lays the groundwork for solving related shape optimization problems. △ Less

Submitted 22 September, 2020; v1 submitted 13 March, 2020; originally announced March 2020.

MSC Class: 65N30; 78M10; 78M30; 35P30

Correlated Insulating States and Transport Signature of Superconductivity in Twisted Trilayer Graphene Moiré of Moiré Superlattices

Authors: Kan-Ting Tsai, Xi Zhang, Ziyan Zhu, Yujie Luo, Stephen Carr, Mitchell Luskin, Efthimios Kaxiras, Ke Wang

Abstract: Layers of two-dimensional materials stacked with a small twist-angle give rise to beating periodic patterns on a scale much larger than the original lattice, referred to as a moiré superlattice. When the stacking involves more than two layers with independent twist angles between adjacent layers, it generates moiré of moiré superlattices, with multiple length scales that control the system's behav… ▽ More Layers of two-dimensional materials stacked with a small twist-angle give rise to beating periodic patterns on a scale much larger than the original lattice, referred to as a moiré superlattice. When the stacking involves more than two layers with independent twist angles between adjacent layers, it generates moiré of moiré superlattices, with multiple length scales that control the system's behavior. Here we demonstrate these effects of a high-order moiré superlattice in twisted trilayer graphene with two consecutive small twist angles. We report correlated insulating states near the half filling of the moiré of moiré superlattice at an extremely low carrier density (~1010 cm-2), near which we also report a zero-resistance transport behavior typically expected in a 2D superconductor. Moreover, the temperature dependence of the measured resistances at full-occupancy (v = -4 and v = 4) states are semi-metallic, distinct from the insulating behavior of twisted bilayer systems, providing the first demonstration of emergent correlated transport behaviors from continuous, non-isolated higher-order moiré flat bands. Our findings shed new insights into the microscopic mechanisms of moiré correlated states and provide the impetus for future studies on this material platform, such as the demonstration of phase coherence and Meissner-like effect. △ Less

Submitted 1 December, 2020; v1 submitted 6 December, 2019; originally announced December 2019.

Moiré of Moiré: Modeling Mechanical Relaxation in Incommensurate Trilayer van der Waals Heterostructures

Authors: Ziyan Zhu, Paul Cazeaux, Mitchell Luskin, Efthimios Kaxiras

Abstract: The incommensurate stacking of multi-layered two-dimensional materials is a challenging problem from a theoretical perspective and an intriguing avenue for manipulating their physical properties. Here we present a multi-scale model to obtain the mechanical relaxation pattern of twisted trilayer van der Waals (vdW) heterostructures with two independent twist angles, a generally incommensurate syste… ▽ More The incommensurate stacking of multi-layered two-dimensional materials is a challenging problem from a theoretical perspective and an intriguing avenue for manipulating their physical properties. Here we present a multi-scale model to obtain the mechanical relaxation pattern of twisted trilayer van der Waals (vdW) heterostructures with two independent twist angles, a generally incommensurate system without a supercell description. We adopt the configuration space as a natural description of such incommensurate layered materials, based on the local environment of atomic positions, bypassing the need for commensurate approximations. To obtain the relaxation pattern, we perform energy minimization with respect to the relaxation displacement vectors. We use a continuum model in combination with the Generalized Stacking Fault energy to describe the interlayer coupling, obtained from first-principles calculations based on Density Functional Theory. We show that the relaxation patterns of twisted trilayer graphene and $\mathrm{WSe_2}$ are "moiré of moiré", as a result of the incommensurate coupling two bilayer moiré patterns. We also show that, in contrast to the symmetry-preserving in-plane relaxation in twisted bilayers, trilayer relaxation can break the two-fold rotational symmetry about the xy-plane when the two twist angles are equal. △ Less

Submitted 20 April, 2020; v1 submitted 13 November, 2019; originally announced November 2019.

Journal ref: Phys. Rev. B 101, 224107 (2020)

Efficient Computation of Kubo Conductivity for Incommensurate 2D Heterostructures

Authors: Daniel Massatt, Stephen Carr, Mitchell Luskin

Abstract: Here we introduce a numerical method for computing conductivity via the Kubo Formula for incommensurate 2D bilayer heterostructures using a tight-binding framework. We begin with deriving the momentum space formulation and Kubo Formula from the real space tight-binding model using the appropriate Bloch transformation operator. We further discuss the resulting algorithm along with its convergence r… ▽ More Here we introduce a numerical method for computing conductivity via the Kubo Formula for incommensurate 2D bilayer heterostructures using a tight-binding framework. We begin with deriving the momentum space formulation and Kubo Formula from the real space tight-binding model using the appropriate Bloch transformation operator. We further discuss the resulting algorithm along with its convergence rate and computation cost in terms of parameters such as relaxation time and temperature. In particular, we show that for low frequencies, low temperature, and long relaxation times conductivity can be computed very efficiently using momentum space for a wide class of materials. We then demonstrate our method by computing conductivity for twisted bilayer graphene (tBLG) for small twist angles. △ Less

Submitted 24 October, 2019; originally announced October 2019.

Comments: 16 pages, 3 figures

Nonretarded edge plasmon-polaritons in anisotropic two-dimensional materials

Authors: Dionisios Margetis, Matthias Maier, Tobias Stauber, Tony Low, Mitchell Luskin

Abstract: By an integral equation approach to the time-harmonic classical Maxwell equations, we describe the dispersion in the nonretarded frequency regime of the edge plasmon-polariton (EPP) on a semi-infinite flat sheet. The sheet has an arbitrary, physically admissible, tensor valued and spatially homogeneous conductivity, and serves as a model for a family of two-dimensional conducting materials. We for… ▽ More By an integral equation approach to the time-harmonic classical Maxwell equations, we describe the dispersion in the nonretarded frequency regime of the edge plasmon-polariton (EPP) on a semi-infinite flat sheet. The sheet has an arbitrary, physically admissible, tensor valued and spatially homogeneous conductivity, and serves as a model for a family of two-dimensional conducting materials. We formulate a system of integral equations for the electric field tangential to the sheet in a homogeneous and isotropic ambient medium. We show how this system is simplified via a length scale separation. This view entails the quasi-electrostatic approximation, by which the tangential electric field is replaced by the gradient of a scalar potential, $\varphi$. By the Wiener-Hopf method, we solve an integral equation for $\varphi$ in some generality. The EPP dispersion relation comes from the elimination of a divergent limiting Fourier integral for $\varphi$ at the edge. We connect the existence, or lack thereof, of the EPP dispersion relation to the index for Wiener-Hopf integral equations, an integer of topological character. We indicate that the values of this index may express an asymmetry due to the material anisotropy in the number of wave modes propagating on the sheet away from the edge with respect to the EPP direction of propagation. We discuss extensions such as the setting of two semi-infinite, coplanar sheets. Our theory forms a generalization of the treatment by Volkov and Mikhailov (1988 Sov. Phys. JETP 67 1639). △ Less

Submitted 10 October, 2019; originally announced October 2019.

Modeling and Computation of Kubo Conductivity for 2D Incommensurate Bilayers

Authors: Simon Etter, Daniel Massatt, Mitchell Luskin, Christoph Ortner

Abstract: This paper presents a unified approach to the modeling and computation of the Kubo conductivity of incommensurate bilayer heterostructures at finite temperature. Firstly, we derive an expression for the large-body limit of Kubo-Greenwood conductivity in terms of an integral of the conductivity function with respect to a current-current correlation measure. We then observe that the incommensurate s… ▽ More This paper presents a unified approach to the modeling and computation of the Kubo conductivity of incommensurate bilayer heterostructures at finite temperature. Firstly, we derive an expression for the large-body limit of Kubo-Greenwood conductivity in terms of an integral of the conductivity function with respect to a current-current correlation measure. We then observe that the incommensurate structure can be exploited to decompose the current-current correlation measure into local contribution and deduce an approximation scheme which is exponentially convergent in terms of domain size. Secondly, we analyze the cost of computing local conductivities via Chebyshev approximation. Our main finding is that if the inverse temperature $β$ is sufficiently small compared to the inverse relaxation time $η$, namely $β\lesssim η^{-1/2}$, then the dominant computational cost is $\mathcal{O}\bigl(η^{-3/2}\bigr)$ inner products for a suitably truncated Chebyshev series, which significantly improves on the $\mathcal{O}\bigl(η^{-2}\bigr)$ inner products required by a naive Chebyshev approximation. Thirdly, we propose a rational approximation scheme for the low temperature regime $η^{-1/2} \lesssim β$, where the cost of the polynomial method increases up to $\mathcal{O}\bigl(β^2\bigr),$ but the rational scheme scales much more mildly with respect to $β$. △ Less

Submitted 23 September, 2020; v1 submitted 2 July, 2019; originally announced July 2019.

Comments: 40 pages, 10 figures

Homogenization of plasmonic crystals: Seeking the epsilon-near-zero effect

Authors: Matthias Maier, Marios Mattheakis, Efthimios Kaxiras, Mitchell Luskin, Dionisios Margetis

Abstract: By using an asymptotic analysis and numerical simulations, we derive and investigate a system of homogenized Maxwell's equations for conducting material sheets that are periodically arranged and embedded in a heterogeneous and anisotropic dielectric host. This structure is motivated by the need to design plasmonic crystals that enable the propagation of electromagnetic waves with no phase delay (e… ▽ More By using an asymptotic analysis and numerical simulations, we derive and investigate a system of homogenized Maxwell's equations for conducting material sheets that are periodically arranged and embedded in a heterogeneous and anisotropic dielectric host. This structure is motivated by the need to design plasmonic crystals that enable the propagation of electromagnetic waves with no phase delay (epsilon-near-zero effect). Our microscopic model incorporates the surface conductivity of the two-dimensional (2D) material of each sheet and a corresponding line charge density through a line conductivity along possible edges of the sheets. Our analysis generalizes averaging principles inherent in previous Bloch-wave approaches. We investigate physical implications of our findings. In particular, we emphasize the role of the vector-valued corrector field, which expresses microscopic modes of surface waves on the 2D material. We demonstrate how our homogenization procedure may set the foundation for computational investigations of: effective optical responses of reasonably general geometries, and complicated design problems in the plasmonics of 2D materials. △ Less

Submitted 19 April, 2019; v1 submitted 21 September, 2018; originally announced September 2018.

MSC Class: 35B27; 35Q60; 74Q10

Adaptive finite element simulations of waveguide configurations involving parallel 2D material sheets

Authors: Jung Heon Song, Matthias Maier, Mitchell Luskin

Abstract: We discuss analytically and numerically the propagation and energy transmission of electromagnetic waves caused by the coupling of surface plasmon polaritons (SPPs) between two spatially separated layers of 2D materials, such as graphene, at subwavelength distances. We construct an adaptive finite-element method to compute the ratio of energy transmitted within these waveguide structures reliably… ▽ More We discuss analytically and numerically the propagation and energy transmission of electromagnetic waves caused by the coupling of surface plasmon polaritons (SPPs) between two spatially separated layers of 2D materials, such as graphene, at subwavelength distances. We construct an adaptive finite-element method to compute the ratio of energy transmitted within these waveguide structures reliably and efficiently. At its heart, the method is built upon a goal-oriented a posteriori error estimation with the dual-weighted residual method (DWR). Further, we derive analytic solutions of the two-layer system, compare those to (known) single-layer configurations, and compare and validate our numerical findings by comparing numerical and analytical values for optimal spacing of the two-layer configuration. Additional aspects of our numerical treatment, such as local grid refinement, and the utilization of perfectly matched layers (PMLs) are examined in detail. △ Less

Submitted 11 March, 2019; v1 submitted 17 September, 2018; originally announced September 2018.

MSC Class: 65N30; 78M10; 78M30; 78A45

Nonperturbative nonlinear effects in the dispersion relations for TE and TM plasmons on two-dimensional materials

Authors: Vera Andreeva, Mitchell Luskin, Dionisios Margetis

Abstract: We analytically obtain the dispersion relations for transverse-electric (TE) and transverse-magnetic (TM) surface plasmon-polaritons in a nonlinear two-dimensional (2D) conducting material lying between two Kerr-type dielectric media. To this end, we use Maxwell's equations within the quasi-electrostatic, weakly dissipative regime. We show that the wavelength and propagation distance of surface pl… ▽ More We analytically obtain the dispersion relations for transverse-electric (TE) and transverse-magnetic (TM) surface plasmon-polaritons in a nonlinear two-dimensional (2D) conducting material lying between two Kerr-type dielectric media. To this end, we use Maxwell's equations within the quasi-electrostatic, weakly dissipative regime. We show that the wavelength and propagation distance of surface plasmons decrease due to the nonlinearity of the surrounding dielectric. In contrast, the effect of the nonlinearity of the 2D material depends on the signs of the real and imaginary parts of the third-order conductivity. Notably, the dispersion relations obtained by naively replacing the permittivity of the dielectric medium by its nonlinear counterpart in the respective dispersion relations of the linear regime are not accurate. We apply our analysis to the case of doped graphene and make predictions for the TM-polarized surface plasmon wavelength and propagation distance. △ Less

Submitted 16 October, 2018; v1 submitted 4 August, 2018; originally announced August 2018.

Comments: 18 pages, 8 figures

Journal ref: Phys. Rev. B 98, 195407 (2018)

Switchable and unidirectional plasmonic beacons in hyperbolic 2D materials

Authors: Andrei Nemilentsau, Tobias Stauber, Guillermo Gómez-Santos, Mitchell Luskin, Tony Low

Abstract: In hyperbolic 2D materials, energy is channeled to their deep subwavelength polaritonic modes via four narrow beams. Here we consider the launching of surface polaritons in the hyperbolic 2D materials and demonstrate that efficient uni-directional excitation is possible with an elliptically polarized electric dipole, with the optimal choice of dipole ellipticity depending on the materials optical… ▽ More In hyperbolic 2D materials, energy is channeled to their deep subwavelength polaritonic modes via four narrow beams. Here we consider the launching of surface polaritons in the hyperbolic 2D materials and demonstrate that efficient uni-directional excitation is possible with an elliptically polarized electric dipole, with the optimal choice of dipole ellipticity depending on the materials optical constants. The selection rules afforded by the choice of dipole polarization allow turning off up to two beams, and even three if the dipole is placed close to an edge. This makes the dipole a directionally switchable beacon for the launching of sub-difractional polaritonic beams, a potential logical gate. We develop an analytical approximation of the excitation process which describes the results of the numerical simulations well and affords a simple physical interpretation. △ Less

Submitted 6 June, 2019; v1 submitted 2 August, 2018; originally announced August 2018.

Comments: 20 pages, 13 figures

Journal ref: Phys. Rev. B 99, 201405(R) (2019)

Stability and convergence of the string method for computing minimum energy paths

Authors: Brian Van Koten, Mitchell Luskin

Abstract: We analyze the convergence of the string method of E, Ren, and Vanden-Eijnden to a minimum energy path. Under some assumptions relating to the critical points on the minimum energy path, we show that the string method initialized in a neighborhood of the minimum energy path converges to an arbitrarily small neighborhood of the minimum energy path as the number of images is increased. We analyze the convergence of the string method of E, Ren, and Vanden-Eijnden to a minimum energy path. Under some assumptions relating to the critical points on the minimum energy path, we show that the string method initialized in a neighborhood of the minimum energy path converges to an arbitrarily small neighborhood of the minimum energy path as the number of images is increased. △ Less

Submitted 10 August, 2018; v1 submitted 16 July, 2018; originally announced July 2018.

Comments: 26 pages, 4 figures

MSC Class: 65

Energy minimization of 2D incommensurate heterostructures

Authors: Paul Cazeaux, Mitchell Luskin, Daniel Massatt

Abstract: We derive and analyze a novel approach for modeling and computing the mechanical relaxation of incommensurate 2D heterostructures. Our approach parametrizes the relaxation pattern by the compact local configuration space rather than real space, thus bypassing the need for the standard supercell approximation and giving a true aperiodic atomistic configuration. Our model extends the computationally… ▽ More We derive and analyze a novel approach for modeling and computing the mechanical relaxation of incommensurate 2D heterostructures. Our approach parametrizes the relaxation pattern by the compact local configuration space rather than real space, thus bypassing the need for the standard supercell approximation and giving a true aperiodic atomistic configuration. Our model extends the computationally accessible regime of weakly coupled bilayers with similar orientations or lattice spacing, for example materials with a small relative twist where the widely studied large-scale moire patterns arise. Our model also makes possible the simulation of multi-layers for which no interlayer empirical atomistic potential exists, such as those composed of MoS2 layers, and more generally makes possible the simulation of the relaxation of multi-layer heterostructures for which a planar moire pattern does not exist. △ Less

Submitted 3 August, 2019; v1 submitted 27 June, 2018; originally announced June 2018.

Comments: 30 pages, 4 figures

MSC Class: 74; 70

Relaxation and Domain Formation in Incommensurate 2D Heterostructures

Authors: Stephen Carr, Daniel Massatt, Steven B. Torrisi, Paul Cazeaux, Mitchell Luskin, Efthimios Kaxiras

Abstract: We introduce configuration space as a natural representation for calculating the mechanical relaxation patterns of incommensurate two-dimensional (2D) bilayers, bypassing supercell approximations to encompass aperiodic relaxation patterns. The approach can be applied to a wide variety of 2D materials through the use of a continuum model in combination with a generalized stacking fault energy for i… ▽ More We introduce configuration space as a natural representation for calculating the mechanical relaxation patterns of incommensurate two-dimensional (2D) bilayers, bypassing supercell approximations to encompass aperiodic relaxation patterns. The approach can be applied to a wide variety of 2D materials through the use of a continuum model in combination with a generalized stacking fault energy for interlayer interactions. We present computational results for small-angle twisted bilayer graphene and molybdenum disulfide (MoS$_2$), a representative material of the transition metal dichalcogenide (TMDC) family of 2D semiconductors. We calculate accurate relaxations for MoS$_2$ even at small twist-angle values, enabled by the fact that our approach does not rely on empirical atomistic potentials for interlayer coupling. The results demonstrate the efficiency of the configuration space method by computing relaxations with minimal computational cost for twist angles down to $0.05^\circ$, which is smaller than what can be explored by any available real space techniques. We also outline a general explanation of domain formation in 2D bilayers with nearly-aligned lattices, taking advantage of the relationship between real space and configuration space. △ Less

Submitted 17 May, 2018; originally announced May 2018.

Comments: 6 pages, 4 figures

Journal ref: Phys. Rev. B 98, 224102 (2018)

Atomic and electronic reconstruction at van der Waals interface in twisted bilayer graphene

Authors: Hyobin Yoo, Rebecca Engelke, Stephen Carr, Shiang Fang, Kuan Zhang, Paul Cazeaux, Suk Hyun Sung, Robert Hovden, Adam W. Tsen, Takashi Taniguchi, Kenji Watanabe, Gyu-Chul Yi, Miyoung Kim, Mitchell Luskin, Ellad B. Tadmor, Efthimios Kaxiras, Philip Kim

Abstract: Control of the interlayer twist angle in two-dimensional (2D) van der Waals (vdW) heterostructures enables one to engineer a quasiperiodic moiré superlattice of tunable length scale. In twisted bilayer graphene (TBG), the simple moiré superlattice band description suggests that the electronic band width can be tuned to be comparable to the vdW interlayer interaction at a 'magic angle', exhibiting… ▽ More Control of the interlayer twist angle in two-dimensional (2D) van der Waals (vdW) heterostructures enables one to engineer a quasiperiodic moiré superlattice of tunable length scale. In twisted bilayer graphene (TBG), the simple moiré superlattice band description suggests that the electronic band width can be tuned to be comparable to the vdW interlayer interaction at a 'magic angle', exhibiting strongly correlated behavior. However, the vdW interlayer interaction can also cause significant structural reconstruction at the interface by favoring interlayer commensurability, which competes with the intralayer lattice distortion. Here we report the atomic scale reconstruction in TBG and its effect on the electronic structure. We find a gradual transition from incommensurate moiré structure to an array of commensurate domain structures as we decrease the twist angle across the characteristic crossover angle, $θ_c$ ~1°. In the twist regime smaller than $θ_c$ where the atomic and electronic reconstruction become significant, a simple moiré band description breaks down. Upon applying a transverse electric field, we observe electronic transport along the network of one-dimensional (1D) topological channels that surround the alternating triangular gapped domains, providing a new pathway to engineer the system with continuous tunability. △ Less

Submitted 2 October, 2018; v1 submitted 11 April, 2018; originally announced April 2018.

Duality between atomic configurations and Bloch states in twistronic materials

Authors: Stephen Carr, Daniel Massatt, Mitchell Luskin, Efthimios Kaxiras

Abstract: The relative orientation (twist) of successive layers of stacked two-dimensional (2D) materials creates variations in the interlayer atomic registry. The variations often form a super lattice, called a moiré pattern, which can alter electronic properties. In this work we introduce a classification of the single-particle electronic structures that can occur in twisted stacks of 2D layers by charact… ▽ More The relative orientation (twist) of successive layers of stacked two-dimensional (2D) materials creates variations in the interlayer atomic registry. The variations often form a super lattice, called a moiré pattern, which can alter electronic properties. In this work we introduce a classification of the single-particle electronic structures that can occur in twisted stacks of 2D layers by characterizing them as "moiré molecules" or "moiré crystals". The molecules generate localized electronic states and moiré flat bands, while the crystals are sometimes unconventional and produce electronic banding in the configuration basis. The underpinning of this classification is the duality between interlayer configuration and monolayer Bloch momentum in moiré Hamiltonians. We apply this understanding to diagrams of local electron density in untwisted geometries to produce intuitive and quantitative predictions of twistronic properties. We provide a conceptual introduction to this framework through a one-dimensional model, and then apply it to 2D twisted bilayers of the semi-metal graphene, and of MoS$_2$, a representative material of the transition metal dichalcogenide (TMDC) family of semiconductors. This level of thorough understanding of twistronic phenomena is vital in the search for new material platforms for localized moiré electrons. △ Less

Submitted 13 July, 2020; v1 submitted 3 March, 2018; originally announced March 2018.

Comments: 13 pages, 9 figures

Journal ref: Phys. Rev. Research 2, 033162 (2020)

Universal behavior of dispersive Dirac cone in gradient-index plasmonic metamaterials

Authors: Matthias Maier, Marios Mattheakis, Efthimios Kaxiras, Mitchell Luskin, Dionisios Margetis

Abstract: We demonstrate analytically and numerically that the dispersive Dirac cone emulating an epsilon-near-zero (ENZ) behavior is a universal property within a family of plasmonic crystals consisting of two-dimensional (2D) metals. Our starting point is a periodic array of 2D metallic sheets embedded in an inhomogeneous and anisotropic dielectric host that allows for propagation of transverse-magnetic (… ▽ More We demonstrate analytically and numerically that the dispersive Dirac cone emulating an epsilon-near-zero (ENZ) behavior is a universal property within a family of plasmonic crystals consisting of two-dimensional (2D) metals. Our starting point is a periodic array of 2D metallic sheets embedded in an inhomogeneous and anisotropic dielectric host that allows for propagation of transverse-magnetic (TM) polarized waves. By invoking a systematic bifurcation argument for arbitrary dielectric profiles in one spatial dimension, we show how TM Bloch waves experience an effective dielectric function that averages out microscopic details of the host medium. The corresponding effective dispersion relation reduces to a Dirac cone when the conductivity of the metallic sheet and the period of the array satisfy a critical condition for ENZ behavior. Our analytical findings are in excellent agreement with numerical simulations. △ Less

Submitted 12 January, 2018; v1 submitted 6 November, 2017; originally announced November 2017.

Journal ref: Phys. Rev. B 97, 035307 (2018)

Ultracompact amplitude modulator by coupling hyperbolic polaritons over a graphene-covered gap

Authors: Matthias Maier, Andrei Nemilentsau, Tony Low, Mitchell Luskin

Abstract: The hyperbolic phonon-polaritons within the Reststrahlen band of hBN are of great interest for applications in nanophotonics as they are capable of propagating light signals with low losses over large distances. However, due to the phononic nature of the polaritons in hBN, amplitude modulation of its signal proves to be difficult and has been underexplored. In this paper, we propose theoretically… ▽ More The hyperbolic phonon-polaritons within the Reststrahlen band of hBN are of great interest for applications in nanophotonics as they are capable of propagating light signals with low losses over large distances. However, due to the phononic nature of the polaritons in hBN, amplitude modulation of its signal proves to be difficult and has been underexplored. In this paper, we propose theoretically a broadband efficient amplitude modulator for hyperbolic rays in hBN operating in the frequency range between 1450 cm$^{-1}$ and 1550 cm$^{-1}$. The modulating region comprises a few tens of nanometers wide gap carved within the hBN slab and covered by a graphene layer, where electrostatically gated graphene serves as a mediator that facilitates the coupling between phonon-polaritons on each side of the gap through plasmonic modes within graphene. We demonstrate that such an ultra compact modulator has insertion losses as low as 3 dB and provides modulation depth varying between 14 and 20 dB within the type-II hyperbolicity region of hBN. △ Less

Submitted 13 November, 2017; v1 submitted 19 September, 2017; originally announced September 2017.

Incommensurate Heterostructures in Momentum Space

Authors: Daniel Massatt, Stephen Carr, Mitchell Luskin, Christoph Ortner

Abstract: To make the investigation of electronic structure of incommensurate heterostructures computationally tractable, effective alternatives to Bloch theory must be developed. In Massatt2017, we developed and analyzed a real space scheme that exploits spatial ergodicity and near-sightedness. In the present work, we present an analogous scheme formulated in momentum space, which we prove have significant… ▽ More To make the investigation of electronic structure of incommensurate heterostructures computationally tractable, effective alternatives to Bloch theory must be developed. In Massatt2017, we developed and analyzed a real space scheme that exploits spatial ergodicity and near-sightedness. In the present work, we present an analogous scheme formulated in momentum space, which we prove have significant computational advantages in specific incommensurate systems of physical interest, e.g., bilayers of a specified class of materials with small rotation angles. We use our theoretical analysis to obtain estimates for improved rates of convergence with respect to total CPU time for our momentum space method that are confirmed in computational experiments. △ Less

Submitted 7 August, 2017; originally announced August 2017.

Spin-Diffusions and Diffusive Molecular Dynamics

Authors: Brittan A Farmer, Mitchell Luskin, Petr Plecháč, Gideon Simpson

Abstract: Metastable condensed matter typically fluctuates about local energy minima at the femtosecond time scale before transitioning between local minima after nanoseconds or microseconds. This vast scale separation limits the applicability of classical molecular dynamics methods and has spurned the development of a host of approximate algorithms. One recently proposed method is diffusive molecular dynam… ▽ More Metastable condensed matter typically fluctuates about local energy minima at the femtosecond time scale before transitioning between local minima after nanoseconds or microseconds. This vast scale separation limits the applicability of classical molecular dynamics methods and has spurned the development of a host of approximate algorithms. One recently proposed method is diffusive molecular dynamics which aims to integrate a system of ordinary differential equations describing the likelihood of occupancy by one of two species, in the case of a binary alloy, while quasistatically evolving the locations of the atoms. While diffusive molecular dynamics has shown to be efficient and provide agreement with observations, it is fundamentally a model, with unclear connections to classical molecular dynamics. In this work, we formulate a spin-diffusion stochastic process and show how it can be connected to diffusive molecular dynamics. The spin-diffusion model couples a classical overdamped Langevin equation to a kinetic Monte Carlo model for exchange amongst the species of a binary alloy. Under suitable assumptions and approximations, spin-diffusion can be shown to lead to diffusive molecular molecular dynamics type models. The key assumptions and approximations include a well defined time scale separation, a choice of spin exchange rates, a low temperature approximation, and a mean field type approximation. We derive several models from different assumptions and show their relationship to diffusive molecular dynamics. Differences and similarities amongst the models are explored in a simple test problem. △ Less

Submitted 23 September, 2017; v1 submitted 5 February, 2017; originally announced February 2017.

Comments: 26 Pages, corrected a typo from the previous version

MSC Class: 60G20 82C80 37N15 74H15

Generation of surface plasmon-polaritons by edge effects

Authors: Matthias Maier, Dionisios Margetis, Mitchell Luskin

Abstract: By using numerical and analytical methods, we describe the generation of fine-scale lateral electromagnetic waves, called surface plasmon-polaritons (SPPs), on atomically thick, metamaterial conducting sheets in two spatial dimensions (2D). Our computations capture the two-scale character of the total field and reveal how each edge of the sheet acts as a source of an SPP that may dominate the diff… ▽ More By using numerical and analytical methods, we describe the generation of fine-scale lateral electromagnetic waves, called surface plasmon-polaritons (SPPs), on atomically thick, metamaterial conducting sheets in two spatial dimensions (2D). Our computations capture the two-scale character of the total field and reveal how each edge of the sheet acts as a source of an SPP that may dominate the diffracted field. We use the finite element method to numerically implement a variational formulation for a weak discontinuity of the tangential magnetic field across a hypersurface. An adaptive, local mesh refinement strategy based on a posteriori error estimators is applied to resolve the pronounced two-scale character of wave propagation and radiation over the metamaterial sheet. We demonstrate by numerical examples how a singular geometry, e.g., sheets with sharp edges, and sharp spatial changes in the associated surface conductivity may significantly influence surface plasmons in nanophotonics. △ Less

Submitted 2 February, 2017; originally announced February 2017.

MSC Class: 65N30; 78M10; 78M30; 78A45

On the Wiener-Hopf method for surface plasmons: Diffraction from semi-infinite metamaterial sheet

Authors: Dionisios Margetis, Matthias Maier, Mitchell Luskin

Abstract: By formally invoking the Wiener-Hopf method, we explicitly solve a one-dimensional, singular integral equation for the excitation of a slowly decaying electromagnetic wave, called surface plasmon-polariton (SPP), of small wavelength on a semi-infinite, flat conducting sheet irradiated by a plane wave in two spatial dimensions. This setting is germane to wave diffraction by edges of large sheets of… ▽ More By formally invoking the Wiener-Hopf method, we explicitly solve a one-dimensional, singular integral equation for the excitation of a slowly decaying electromagnetic wave, called surface plasmon-polariton (SPP), of small wavelength on a semi-infinite, flat conducting sheet irradiated by a plane wave in two spatial dimensions. This setting is germane to wave diffraction by edges of large sheets of single-layer graphene. Our analytical approach includes: (i) formulation of a functional equation in the Fourier domain; (ii) evaluation of a split function, which is expressed by a contour integral and is a key ingredient of the Wiener-Hopf factorization; and (iii) extraction of the SPP as a simple-pole residue of a Fourier integral. Our analytical solution is in good agreement with a finite-element numerical computation. △ Less

Submitted 10 January, 2017; originally announced January 2017.

MSC Class: 35Q61; 45E10; 45K05

Generalized Kubo Formulas for the Transport Properties of Incommensurate 2D Atomic Heterostructures

Authors: Eric Cancès, Paul Cazeaux, Mitchell Luskin

Abstract: We give an exact formulation for the transport coefficients of incommensurate two-dimensional atomic multilayer systems in the tight-binding approximation. This formulation is based upon the C* algebra framework introduced by Bellissard and collaborators to study aperiodic solids (disordered crystals, quasicrystals, and amorphous materials), notably in the presence of magnetic fields (quantum Hall… ▽ More We give an exact formulation for the transport coefficients of incommensurate two-dimensional atomic multilayer systems in the tight-binding approximation. This formulation is based upon the C* algebra framework introduced by Bellissard and collaborators to study aperiodic solids (disordered crystals, quasicrystals, and amorphous materials), notably in the presence of magnetic fields (quantum Hall effect). We also present numerical approximations and test our methods on a one-dimensional incommensurate bilayer system. △ Less

Submitted 23 November, 2016; originally announced November 2016.

Twistronics: Manipulating the Electronic Properties of Two-dimensional Layered Structures through their Twist Angle

Authors: Stephen Carr, Daniel Massatt, Shiang Fang, Paul Cazeaux, Mitchell Luskin, Efthimios Kaxiras

Abstract: The ability in experiments to control the relative twist angle between successive layers in two-dimensional (2D) materials offers a new approach to manipulating their electronic properties; we refer to this approach as "twistronics". A major challenge to theory is that, for arbitrary twist angles, the resulting structure involves incommensurate (aperiodic) 2D lattices. Here, we present a general m… ▽ More The ability in experiments to control the relative twist angle between successive layers in two-dimensional (2D) materials offers a new approach to manipulating their electronic properties; we refer to this approach as "twistronics". A major challenge to theory is that, for arbitrary twist angles, the resulting structure involves incommensurate (aperiodic) 2D lattices. Here, we present a general method for the calculation of the electronic density of states of aperiodic 2D layered materials, using parameter-free hamiltonians derived from ab initio density-functional theory. We use graphene, a semimetal, and MoS$_2$, a representative of the transition metal dichalcogenide (TMDC) family of 2D semiconductors, to illustrate the application of our method, which enables fast and efficient simulation of multi-layered stacks in the presence of local disorder and external fields. We comment on the interesting features of their Density of States (DoS) as a function of twist-angle and local configuration and on how these features can be experimentally observed. △ Less

Submitted 3 November, 2016; v1 submitted 2 November, 2016; originally announced November 2016.

Comments: 7 pages, 4 figures, to be submitted to Phys. Rev. B

Journal ref: Phys. Rev. B 95, 075420 (2017)

Electronic Density of States for Incommensurate Layers

Authors: Daniel Massatt, Mitchell Luskin, Christoph Ortner

Abstract: We prove that the electronic density of states (DOS) for 2D incommensurate layered structures, where Bloch theory does not apply, is well-defined as the thermodynamic limit of finite clusters. In addition, we obtain an explicit representation formula for the DOS as an integral over local configurations. Next, based on this representation formula, we propose a novel algorithm for computing electron… ▽ More We prove that the electronic density of states (DOS) for 2D incommensurate layered structures, where Bloch theory does not apply, is well-defined as the thermodynamic limit of finite clusters. In addition, we obtain an explicit representation formula for the DOS as an integral over local configurations. Next, based on this representation formula, we propose a novel algorithm for computing electronic structure properties in incommensurate heterostructures, which overcomes limitations of the common approach to artificially strain a large supercell and then apply Bloch theory. △ Less

Submitted 5 August, 2016; originally announced August 2016.

Comments: 20 pages, 9 figures

MSC Class: 91

Analysis of rippling in incommensurate one-dimensional coupled chains

Authors: Paul Cazeaux, Mitchell Luskin, Ellad B. Tadmor

Abstract: Graphene and other recently developed 2D materials exhibit exceptionally strong in-plane stiffness. Relaxation of few-layer structures, either free-standing or on slightly mismatched substrates occurs mostly through out-of-plane bending and the creation of large-scale ripples. In this work, we present a novel double chain model, where we allow relaxation to occur by bending of the incommensurate c… ▽ More Graphene and other recently developed 2D materials exhibit exceptionally strong in-plane stiffness. Relaxation of few-layer structures, either free-standing or on slightly mismatched substrates occurs mostly through out-of-plane bending and the creation of large-scale ripples. In this work, we present a novel double chain model, where we allow relaxation to occur by bending of the incommensurate coupled system of chains. As we will see, this model can be seen as a new application of the well-known Frenkel-Kontorova model for a one-dimensional atomic chain lying in a periodic potential. We focus in particular on modeling and analyzing ripples occurring in ground state configurations, as well as their numerical simulation. △ Less

Submitted 8 June, 2016; originally announced June 2016.

MSC Class: 65Z05; 70C20; 74E15; 70G75

Dipole excitation of surface plasmon on a conducting sheet: finite element approximation and validation

Authors: Matthias Maier, Dionisios Margetis, Mitchell Luskin

Abstract: We formulate and validate a finite element approach to the propagation of a slowly decaying electromagnetic wave, called surface plasmon-polariton, excited along a conducting sheet, e.g., a single-layer graphene sheet, by an electric Hertzian dipole. By using a suitably rescaled form of time-harmonic Maxwell's equations, we derive a variational formulation that enables a direct numerical treatment… ▽ More We formulate and validate a finite element approach to the propagation of a slowly decaying electromagnetic wave, called surface plasmon-polariton, excited along a conducting sheet, e.g., a single-layer graphene sheet, by an electric Hertzian dipole. By using a suitably rescaled form of time-harmonic Maxwell's equations, we derive a variational formulation that enables a direct numerical treatment of the associated class of boundary value problems by appropriate curl-conforming finite elements. The conducting sheet is modeled as an idealized hypersurface with an effective electric conductivity. The requisite weak discontinuity for the tangential magnetic field across the hypersurface can be incorporated naturally into the variational formulation. We carry out numerical simulations for an infinite sheet with constant isotropic conductivity embedded in two spatial dimensions; and validate our numerics against the closed-form exact solution obtained by the Fourier transform in the tangential coordinate. Numerical aspects of our treatment such as an absorbing perfectly matched layer, as well as local refinement and a-posteriori error control are discussed. △ Less

Submitted 27 April, 2017; v1 submitted 26 May, 2016; originally announced May 2016.

MSC Class: 65N30; 78M10; 78M30; 78A45

Journal ref: J. Comput. Phys. 339 (2017) 126-145

Analysis of a Predictor-Corrector Method for Computationally Efficient Modeling of Surface Effects in 1D

Authors: Andrew J. Binder, Mitchell Luskin, Christoph Ortner

Abstract: The regular Cauchy--Born method is a useful and efficient tool for analyzing bulk properties of materials in the absence of defects. However, the method normally fails to capture surface effects, which are essential to determining material properties at small length scales. In this paper, we present a corrector method that improves upon the prediction for material behavior from the Cauchy--Born me… ▽ More The regular Cauchy--Born method is a useful and efficient tool for analyzing bulk properties of materials in the absence of defects. However, the method normally fails to capture surface effects, which are essential to determining material properties at small length scales. In this paper, we present a corrector method that improves upon the prediction for material behavior from the Cauchy--Born method over a small boundary layer at the surface of a 1D material by capturing the missed surface effects. We justify the separation of the problem into a bulk response and a localized surface correction by establishing an error estimate, which vanishes in the long wavelength limit. △ Less

Submitted 18 May, 2016; originally announced May 2016.

Comments: 28 pages, 5 figures

MSC Class: 65Z05 (Primary); 70C20 (Secondary)

On solutions of Maxwell's equations with dipole sources over a thin conducting film

Authors: Dionisios Margetis, Mitchell Luskin

Abstract: We derive and interpret solutions of time-harmonic Maxwell's equations with a vertical and a horizontal electric dipole near a planar, thin conducting film, e.g. graphene sheet, lying between two unbounded isotropic and non-magnetic media. Exact expressions for all field components are extracted in terms of rapidly convergent series of known transcendental functions when the ambient media have equ… ▽ More We derive and interpret solutions of time-harmonic Maxwell's equations with a vertical and a horizontal electric dipole near a planar, thin conducting film, e.g. graphene sheet, lying between two unbounded isotropic and non-magnetic media. Exact expressions for all field components are extracted in terms of rapidly convergent series of known transcendental functions when the ambient media have equal permittivities and both the dipole and observation point lie on the plane of the film. These solutions are simplified for all distances from the source when the film surface resistivity is large in magnitude compared to the intrinsic impedance of the ambient space. The formulas reveal the analytical structure of two types of waves that can possibly be excited by the dipoles and propagate on the film. One of these waves is intimately related to the surface plasmon-polariton of transverse-magnetic (TM) polarization of plane waves. △ Less

Submitted 14 September, 2015; originally announced September 2015.

Comments: 48 pages, 4 figures

A Theoretical Examination of Diffusive Molecular Dynamics

Authors: Gideon Simpson, Mitchell Luskin, David J. Srolovitz

Abstract: Diffusive molecular dynamics is a novel model for materials with atomistic resolution that can reach diffusive time scales. The main ideas of diffusive molecular dynamics are to first minimize an approximate variational Gaussian free energy of the system with respect to the mean atomic coordinates (averaging over many vibrational periods), and to then to perform a diffusive step where atoms and va… ▽ More Diffusive molecular dynamics is a novel model for materials with atomistic resolution that can reach diffusive time scales. The main ideas of diffusive molecular dynamics are to first minimize an approximate variational Gaussian free energy of the system with respect to the mean atomic coordinates (averaging over many vibrational periods), and to then to perform a diffusive step where atoms and vacancies (or two species in a binary alloy) flow on a diffusive time scale via a master equation. We present a mathematical framework for studying this algorithm based upon relative entropy, or Kullback-Leibler divergence. This adds flexibility in how the algorithm is implemented and interpreted. We then compare our formulation, relying on relative entropy and absolute continuity of measures, to existing formulations. The main difference amongst the equations appears in a model for vacancy diffusion, where additional entropic terms appear in our development. △ Less

Submitted 8 June, 2015; originally announced June 2015.

Comments: 19 pages

MSC Class: 65Z05; 65C20; 82C80

Analysis of an optimization-based atomistic-to-continuum coupling method for point defects

Authors: Derek Olson, Alexander V. Shapeev, Pavel Bochev, Mitchell Luskin

Abstract: We formulate and analyze an optimization-based Atomistic-to-Continuum (AtC) coupling method for problems with point defects. Near the defect core the method employs a potential-based atomistic model, which enables accurate simulation of the defect. Away from the core, where site energies become nearly independent of the lattice position, the method switches to a more efficient continuum model. The… ▽ More We formulate and analyze an optimization-based Atomistic-to-Continuum (AtC) coupling method for problems with point defects. Near the defect core the method employs a potential-based atomistic model, which enables accurate simulation of the defect. Away from the core, where site energies become nearly independent of the lattice position, the method switches to a more efficient continuum model. The two models are merged by minimizing the mismatch of their states on an overlap region, subject to the atomistic and continuum force balance equations acting independently in their domains. We prove that the optimization problem is well-posed and establish error estimates. △ Less

Submitted 14 November, 2014; originally announced November 2014.

Analysis of Transition State Theory Rates upon Spatial Coarse-Graining

Authors: Andrew Binder, Mitchell Luskin, Danny Perez, Arthur F. Voter

Abstract: Spatial multiscale methods have established themselves as useful tools for extending the length scales accessible by conventional statics (i.e., zero temperature molecular dynamics). Recently, extensions of these methods, such as the finite-temperature quasicontinuum (hot-QC) or Coarse-Grained Molecular Dynamics (CGMD) methods, have allowed for multiscale molecular dynamics simulations at finite t… ▽ More Spatial multiscale methods have established themselves as useful tools for extending the length scales accessible by conventional statics (i.e., zero temperature molecular dynamics). Recently, extensions of these methods, such as the finite-temperature quasicontinuum (hot-QC) or Coarse-Grained Molecular Dynamics (CGMD) methods, have allowed for multiscale molecular dynamics simulations at finite temperature. Here, we assess the quality of the long-time dynamics these methods generate by considering canonical transition rates. Specifically, we analyze the transition state theory (TST) rates in CGMD and compare them to the corresponding TST rate of the fully atomistic system. The ability of such an approach to reliably reproduce the TST rate is verified through a relative error analysis, which is then used to highlight the major contributions to the error and guide the choice of degrees of freedom. Finally, our analytical results are compared with numerical simulations for the case of a 1-D chain. △ Less

Submitted 22 September, 2014; originally announced September 2014.

Comments: 25 pages, 7 figures