Search

In some models, periodic configurations can be shown to be stable under, both, global $\ell^2$ or local perturbations. This is not the case for aperiodic media. The specific class of aperiodic media we are interested, in arise from taking two 2D periodic crystals and stacking them parallel at a relative twist. In periodic media, phonons are generalized eigenvectors for a stability operator acting on $\ell^2$, coming from a mechanical energy. The goal of our analysis is to provide phonons in the given class of aperiodic media with meaning. As rigorously established for the 1D Frenkel-Kontorova model and previously applied by one of the authors, we assume that we can parametrize minimizing lattice deformations w.r.t. local perturbations via continuous stacking-periodic functions, for which we previously derived a continuous energy density functional. Such (continuous) energy densities are analytically and computationally much better accessible compared to discrete energy functionals. In order to pass to an $\ell^2$-based energy functional, we also study the offset energy w.r.t. given lattice deformations, under $\ell^1$-perturbations. Our findings show that, in the case of an undeformed bilayer heterostructure, while the energy density can be shown to be stable under the assumption of stability of individual layers, the offset energy fails to be stable in the case of twisted bilayer graphene. We then establish conditions for stability and instability of the offset energy w.r.t. the relaxed lattice. Finally, we show that, in the case of incommensurate bilayer homostructures, i.e., two equal layers, if we choose minimizing deformations according to the global energy density above, the offset energy is stable in the limit of zero twist angle. Consequently, in this case, one can then define phonons as generalized eigenvectors w.r.t. the stability operator associated with the offset energy., Comment: 50 pages, 8 figures

We report a twisted triple bilayer graphene platform consisting of three units of Bernal bilayer graphene (BLG) consecutively twisted at 1.49{\deg} and 1.68{\deg}. We observe inter-moir\'e Hofstadter butterflies from two co-existing moir\'e superlattices and a Hofstadter butterfly from reconstructed moir\'e-of-moir\'e 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\'e and moir\'e-of-moir\'e band-insulator state, characterized by measured strong non-local valley Hall effect (VHE) that sensitively depends on the inter-moir\'e 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\'e-of-moir\'e superlattice, including a new type of moir\'e-of-moir\'e band-insulator states and atomically enhanced moir\'e 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.

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., Comment: 28 pages, 4 figures

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\'e barriers in a tBLG device, utilizing the band-insulator gap of the moir\'e superlattice. A single junction can be tuned to host a chosen combination of conventional pseudo barrier and moir\'e tunnel barriers, from which we demonstrate improved collimation efficiency. By measuring transport through two consecutive moir\'e collimators separated by 1 um, we demonstrate evidence of electron collimation in tBLG in the presence of realistic twist-angle inhomogeneity.

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\'e superlattices and the interactions between them. When both moir\'e superlattices contribute equally to electrical transport, we report a new type of inter-moir\'e Hofstadter butterfly. Its Brown-Zak oscillation corresponds to one of the intermediate quasicrystal length scales of the reconstructed moir\'e of moir\'e (MoM) superlattice, shedding new light on emergent physics from competing atomic orders.

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\'e 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., Comment: 27 pages, 8 figures

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., Comment: 67 pages, 9 figures, V4: Streamlined presentation

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., Comment: 35 pages. Small corrections and improvements to exposition

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., Comment: 21 pages, 14 figures

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., Comment: 35 pages, 10 figures

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\'e 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., Comment: 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

Moir\'e 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\'e electrons requires significant technical work specific to each material, impeding large-scale searches for useful moir\'e 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\'e 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\'e electronic structure from the easily computed SD-LDOS of aligned bilayers. This network can satisfactorily predict moir\'e electronic structures, even for materials that are not included in its training data., Comment: 11 pages, 12 figures

We present a general method to analyze the topological nature of the domain boundary connectivity that appeared in relaxed moir\'e superlattice patterns at the interface of 2-dimensional (2D) van der Waals (vdW) materials. At large enough moir\'e lengths, all moir\'e 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\'e 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\'e 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\'e 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\'e structure from experimental TEM data, which provides a quantitative relation between the various components of heterostrain and vortex-antivortex density in moir\'e systems., Comment: 15 pages with 11 figures

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., Comment: 40 pages, 10 figures

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.

We consider the chiral model of twisted bilayer graphene introduced by Tarnopolsky-Kruchkov-Vishwanath (TKV). TKV have proved that for inverse twist angles $\alpha$ such that the effective Fermi velocity at the moir\'e $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 $\alpha \approx .586$. In this work, we give a proof that the Fermi velocity vanishes for at least one $\alpha$ between $.57$ and $.61$ by rigorously justifying TKV's formal expansion of the Fermi velocity over a sufficiently large interval of $\alpha$ 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 $\alpha$. Together with TKV's work our result proves existence of at least one perfectly flat band of the chiral model., Comment: 61 pages (last 7 pages are supplementary material), 11 figures. Improved after review comments: (1) proved that the Fermi velocity perturbation series converges for $|\alpha| < 1/3$ (2) added proof that round-off error cannot affect the conclusions of our numerical computations (3) more extensive discussion of related work

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., Comment: 18 pages, 7 figures

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\'e superlattices.

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.

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., Comment: 27 pages

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.

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\'e superlattice. When the stacking involves more than two layers with independent twist angles between adjacent layers, it generates moir\'e of moir\'e superlattices, with multiple length scales that control the system's behavior. Here we demonstrate these effects of a high-order moir\'e superlattice in twisted trilayer graphene with two consecutive small twist angles. We report correlated insulating states near the half filling of the moir\'e of moir\'e 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\'e flat bands. Our findings shed new insights into the microscopic mechanisms of moir\'e correlated states and provide the impetus for future studies on this material platform, such as the demonstration of phase coherence and Meissner-like effect.

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\'e of moir\'e", as a result of the incommensurate coupling two bilayer moir\'e 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.

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., Comment: 16 pages, 3 figures

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).

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 $\beta$ is sufficiently small compared to the inverse relaxation time $\eta$, namely $\beta \lesssim \eta^{-1/2}$, then the dominant computational cost is $\mathcal{O}\bigl(\eta^{-3/2}\bigr)$ inner products for a suitably truncated Chebyshev series, which significantly improves on the $\mathcal{O}\bigl(\eta^{-2}\bigr)$ inner products required by a naive Chebyshev approximation. Thirdly, we propose a rational approximation scheme for the low temperature regime $\eta^{-1/2} \lesssim \beta$, where the cost of the polynomial method increases up to $\mathcal{O}\bigl(\beta^2\bigr),$ but the rational scheme scales much more mildly with respect to $\beta$., Comment: 40 pages, 10 figures

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.

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.

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., Comment: 18 pages, 8 figures

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., Comment: 20 pages, 13 figures

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., Comment: 26 pages, 4 figures

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., Comment: 30 pages, 4 figures

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., Comment: 6 pages, 4 figures

Control of the interlayer twist angle in two-dimensional (2D) van der Waals (vdW) heterostructures enables one to engineer a quasiperiodic moir\'e superlattice of tunable length scale. In twisted bilayer graphene (TBG), the simple moir\'e 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\'e structure to an array of commensurate domain structures as we decrease the twist angle across the characteristic crossover angle, $\theta_c$ ~1\deg. In the twist regime smaller than $\theta_c$ where the atomic and electronic reconstruction become significant, a simple moir\'e 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.

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\'e 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\'e molecules" or "moir\'e crystals". The molecules generate localized electronic states and moir\'e 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\'e 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\'e electrons., Comment: 13 pages, 9 figures

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.

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.

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.

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., Comment: 26 Pages, corrected a typo from the previous version

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.

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.

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.

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., Comment: 7 pages, 4 figures, to be submitted to Phys. Rev. B

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., Comment: 20 pages, 9 figures

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.

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.