SEARCH · Engineering Papers
Results for “bilayer films”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Nematic superconductivity in twisted bilayer graphene
Twisted bilayer graphene displays insulating and superconducting phases caused by exceptional flattening of its lowest energy bands. Superconductivity with highest T c appears near half filling of the valence band (n ≈ –2). The data show that in the large part of the superconducting dome near n = –2, threefold lattice rotation symmetry is broken in the superconducting phase, i.e., a superconductor is also a nematic. We perform a comprehensive analysis of superconductivity in twisted-bilayer graphene within an itinerant approach and present a mechanism for nematic superconductivity. We take as an input the fact that at dopings, where superconductivity has been observed, the Fermi energy lies in the vicinity of twist-induced Van Hove singularities in the density of states. We argue that the low-energy physics can be properly described by patch models with either six or twelve Van Hove points. We obtain pairing interactions for the patch models in terms of parameters of the microscopic model for the flat bands, which contains both local and twist-induced nonlocal interactions. Here, we show that the latter give rise to attraction in different superconducting channels. For six Van Hove points, there is just one attractive d-wave channel, and we find chiral d ± id superconducting order, which breaks time-reversal symmetry but leaves the lattice rotation symmetry intact. For twelve Van Hove points, we find two attractive channels, g and i waves, with almost equal coupling constants. We show that both order parameters are nonzero in the ground state and explicitly demonstrate that in this co-existence state the threefold lattice rotation symmetry is broken, i.e., the superconducting state is also a nematic. We find two possible nematic states, one is time-reversal symmetric, the other additionally breaks time-reversal symmetry. We apply our scenario to twisted bilayer graphene near n = –2 and argue that it is applicable also to other systems with two (or more) attractive channels with similar couplings as our reasoning for a nematic superconductivity is based on generic symmetry considerations.
Ultrafast formation of a transient two-dimensional diamondlike structure in twisted bilayer graphene
Due to the absence of matching carbon atoms at honeycomb centers with carbon atoms in adjacent graphene sheets, theorists predicted that a sliding process is needed to form AA, AB', or ABC stacking when directly converting graphite into sp 3 bonded diamond. Here, using twisted bilayer graphene, which naturally provides AA and AB' stacking configurations, we report the ultrafast formation of a transient two-dimensional diamondlike structure (which is not observed in aligned graphene) under femtosecond laser irradiation. This photoinduced phase transition is evidenced by the appearance of bond lengths of 1.94 and 3.14 Å in the time-dependent differential pair distribution function using MeV ultrafast electron diffraction. Molecular dynamics and first-principles calculation indicate that sp 3 bonds nucleate at AA and AB' stacked areas in a moiré pattern. This work sheds light on the direct graphite-to-diamond transformation mechanism, which has not been fully understood for more than 60 years.
Efficient simulation of moiré materials using the density matrix renormalization group
Here we present an infinite density-matrix renormalization group (DMRG) study of an interacting continuum model of twisted bilayer graphene (tBLG) near the magic angle. Because of the long-range Coulomb interaction and the large number of orbital degrees of freedom, tBLG is difficult to study with standard DMRG techniques - even constructing and storing the Hamiltonian already poses a major challenge. To overcome these difficulties, we use a recently developed compression procedure to obtain a matrix product operator representation of the interacting tBLG Hamiltonian which we show is both efficient and accurate even when including the spin, valley, and orbital degrees of freedom. To benchmark our approach, we focus mainly on the spinless, single-valley version of the problem where, at half filling, we find that the ground state is a nematic semimetal. Remarkably, we find that the ground state is essentially a k-space Slater determinant, so that Hartree-Fock and DMRG give virtually identical results for this problem. Our results show that the effects of long-range interactions in magic angle graphene can be efficiently simulated with DMRG and open up a new route for numerically studying strong correlation physics in spinful, two-valley tBLG, and other moiré materials in future work.
Twisted bilayer graphene. I. Matrix elements, approximations, perturbation theory, and a k · p two-band model
We investigate the twisted bilayer graphene (TBG) model of Bistritzer and MacDonald (BM) [Bistritzer and MacDonald, Proc. Natl. Acad. Sci. 108, 12233 (2011)] to obtain an analytic understanding of its energetics and wave functions needed for many-body calculations. We provide an approximation scheme for the wave functions of the BM model, which first elucidates why the BM K M -point centered original calculation containing only four plane waves provides a good analytical value for the first magic angle (θ M ≈ 1°). The approximation scheme also elucidates why most of the many-body matrix elements in the Coulomb Hamiltonian projected to the active bands can be neglected. By applying our approximation scheme at the first magic angle to a Γ M -point centered model of six plane waves, we analytically understand the reason for the small Γ M -point gap between the active and passive bands in the isotropic limit w 0 = w 1 . Furthermore, we analytically calculate the group velocities of the passive bands in the isotropic limit, and show that they are almost doubly degenerate, even away from the Γ M point, where no symmetry forces them to be. Furthermore, moving away from the Γ M and K M points, we provide an explicit analytical perturbative understanding as to why the TBG bands are flat at the first magic angle, despite the first magic angle is defined by only requiring a vanishing K M -point Dirac velocity. We derive analytically a connected “magic manifold” w 1 = $2\sqrt{1 + w^{2}_{0}}$ $-\sqrt{2 + 3w^2_0}$, on which the bands remain extremely flat as w 0 is tuned between the isotropic (w 0 = w 1 ) and chiral (w 0 = 0) limits. We analytically show why going away from the isotropic limit by making w 0 less (but not larger) than w 1 increases the Γ M -point gap between the active and the passive bands. Finally, by perturbation theory, we provide an analytic Γ M point k ∙ p two-band model that reproduces the TBG band structure and eigenstates within a certain w 0 , w 1 parameter range. Further refinement of this model are discussed, which suggest a possible faithful representation of the TBG bands by a two-band Γ M point k ∙ p model in the full w 0 , w 1 parameter range.
Strange metal behavior of the Hall angle in twisted bilayer graphene
Twisted bilayer graphene (TBG) with interlayer twist angles near the magic angle ≈ 1.08° hosts flat bands and exhibits correlated states including Mott-like insulators, superconductivity, and magnetism. Additionally, a linear-in-temperature normal state resistivity in TBG has been attributed to an exotic Planckian dissipation mechanism but can be equally well explained in terms of conventional electron-phonon scattering. To address this issue, we perform combined temperature-dependent transport measurements of both the longitudinal and Hall resistivities in near-magic-angle TBG. While the observed longitudinal resistivity follows linear temperature T dependence consistent with previous reports, the Hall resistance shows an anomalous T dependence with the cotangent of the Hall angle cot Θ H ∝T 2 . Boltzmann theory for quasiparticle transport predicts that both the resistivity and cot Θ H should have the same T dependence, contradicting the observed behavior. This failure of quasiparticle-based theories is reminiscent of other correlated strange metals such as cuprates.
Realization of semiconducting layered multiferroic heterojunctions via asymmetrical magnetoelectric coupling
Two-dimensional (2D) semiconducting multiferroics that can effectively couple magnetic and polarization (P) orders have great interest for both fundamental research and technological applications in nanoscale, which are, however, rare in nature. In this work, we propose a general mechanism to realize semiconducting 2D multiferroics via van der Waals (vdW) heterojunction engineering, as demonstrated in a typical heterostructure consisting of magnetic bilayer CrI 3 (bi-CrI 3 ) and ferroelectric monolayer In 2 Se 3 . Interestingly, the novel indirect orbital coupling between Se 4p and Cr 3d orbitals, intermediated by the interfacial I 5p orbitals, is switchable in the opposite P configurations, resulting in an unexpected mechanism of strong asymmetrical magnetoelectric coupling. Therefore, along with the noticeable ferroelectric energy barrier induced by In 2 Se 3 , the realization of opposite magnetic orders in opposite P configurations can eventually result in the novel multiferroicity in bi– CrI 3 /In 2 Se 3 . Finally, we demonstrate that our mechanism can generally be applied to design other vdW multiferroics even with tunable layer thickness.
Collective excitations and quantum incompressibility in electron-hole bilayers
In this study, we apply quantum continuum mechanics to the calculation of the excitation spectrum of a coupled electron-hole bilayer. The theory expresses excitation energies in terms of ground-state intra- and interlayer pair correlation functions, which are available from quantum Monte Carlo calculations. The final formulas for the collective modes deduced from this approach coincide with the formulas obtained in the “quasilocalized particle approximation” by Kalman et al. [G. Kalman, V. Valtchinov, and K. I. Golden, Phys. Rev. Lett. 82 , 3124 (1999) ], and likewise the theory predicts the existence of gapped excitations in the charged channel, with the gap arising from electron-hole correlation. An immediate consequence of the gap is that the static density-density response function of the charged channel vanishes as q 2 for wave vector q → 0 , rather than linearly in q , as commonly expected. In this sense, the system is incompressible . This feature, which has no analog in the classical electron-hole plasma, is consistent with the existence of an excitonic ground state and implies the existence of a discontinuity in the chemical potential of electrons and holes when the numbers of electrons and holes are equal. It should be experimentally observable by monitoring the densities of electrons and holes in response to potentials that attempt to change these densities in opposite directions.
Heterobilayer moiré magnets: Moiré skyrmions and commensurate-incommensurate transitions
In this work, we study untwisted heterobilayers of ferromagnetic and antiferromagnetic van der Waals materials, with in particular a Dzyaloshinskii-Moriya interaction in the ferromagnetic layer. A continuum low-energy field theory is utilized to study such systems. We develop a phase diagram as a function of the strength of interlayer exchange and Dzyaloshinskii-Moriya interactions, combining perturbative and strong coupling analyses with numerical simulations using Landau-Lifshitz-Gilbert equations. Various moiré-periodic commensurate phases are found, and the commensurate-incommensurate transition is discussed. Among the commensurate phases, we observe an interesting skyrmion lattice phase wherein each moiré unit-cell hosts one skyrmion.
$O(N)$ ab initio calculation scheme for large-scale moiré structures
Here we present a two-step method specifically tailored for band structure calculation of the small-angle moiré-pattern materials which contain tens of thousands of atoms in a unit cell. In the first step, the self-consistent field calculation for the ground state is performed with the O(N) Krylov subspace method implemented in openmx. Second, the crystal momentum-dependent Bloch Hamiltonian and overlap matrix are constructed from the results obtained in the first step and only a small number of eigenvalues near the Fermi energy are solved with shift-invert and Lanczos techniques. By systematically tuning two key parameters, the cutoff radius for electron hopping interaction and the dimension of the Krylov subspace, we obtained the band structures for both rigid and corrugated twisted bilayer graphene structures down to the first magic angle (θ = 1.08°) with high enough accuracy at affordable costs. The band structures are in good agreement with those from tight-binding models, continuum models, plane-wave pseudopotential based ab initio calculations, and experimental observations. This method is also shown to be efficient in twisted double-bilayer graphene and bilayer WSe 2 . We think this two-step method can play a crucial role in other twisted two-dimensional materials, especially those with much more complex band structure and where the effective model is hard to construct.
Trions in twisted bilayer graphene
The strong-coupling phase diagram of magic-angle twisted bilayer graphene (TBG) predicts a series of exact one-particle charge ±1 gapped excitations on top of the integer-filled ferromagnetic ground states. Finite-size exact diagonalization studies showed that these are the lowest charge ±1 excitations in the system (for 10 nm screening length), with the exception of charge +1 at filling -1 in the chiral limit. In the current paper we show that this “trion bound state,” a 3-particle, charge 1 excitation of the insulating ferromagnetic ground state of the projected Hamiltonian of TBG, is the lowest charge +1 overall excitation at ν =-1, and also for some large (≈20 nm) screening lengths at ν =-2 in the chiral limit and with very small binding energy. At other fillings, we show that trion bound states do exist, but only for momentum ranges that do not cover the entire moiré Brillouin zone. The trion bound states (at different momenta) exist for finite parameter range w 0 /w 1 but they all disappear in the continuum far below the realistic values of w 0 /w 1 = 0.8. Thus we find the conditions for the existence of the trion bound state, a good variational wave function for it, and investigate its behavior for different screening lengths, at all integer fillings, on both the electron and hole sides.
First-principles study of bilayer polymeric manganese phthalocyanine
Here we study bilayer manganese phthalocyanine (MnPc) molecules and MnPc polymeric sheets using first-principles simulations with a focus on the magnetic interactions between Mn atoms. We find that the most stable position of the upper layer with respect to the lower layer is shifted about 1/8 of a lattice vector from the center of the bottom layer along the direction toward a nearest-neighbor N atom. The magnetic ground state is the Néel antiferromagnetic (AFM) configuration within a layer and ferromagnetic (FM) between Mn atoms in adjacent layers. In this state, the system becomes a semiconductor with an indirect band gap of 11meV. The strongest interaction is the interlayer coupling between the closest Mn atoms. A maximally localized Wannier analysis suggests that the dominant coupling pathway is Mn-(N,C)-Mn rather than a direct Mn-Mn coupling. The maximum calculated magnetic anisotropy energy is found to be 1.0meV per Mn atom. We also find that the bilayer molecule shows a significant stacking angle change from FM to AFM configurations accompanied by a change of orbital filling ordering.
Magnetic states of graphene proximitized Kitaev materials
Single-layer α-ruthenium trichloride (α–RuCl 3 ) has been proposed as a potential quantum spin liquid. Graphene/RuCl 3 heterobilayers have been extensively studied with a focus on the large interlayer electron transfer that dopes both materials. Here we examine the interplay between the competing magnetic state of RuCl 3 layer and graphene electronic properties. We perform self-consistent Hartree-Fock calculations on a Hubbard-Kanamori model of the 4d 5 t 2g electrons of α–RuCl 3 and confirm that out-of-plane ferromagnetic and zigzag antiferromagnetic states are energetically competitive. We show that the influence of hybridization between graphene and RuCl 3 bands is strongly sensitive to the magnetic configuration of RuCl 3 and the relative orientations of the two layers. We argue that strong hybridization leads to graphene magnetoresistance and that it may tilt the balance between closely competing magnetic states. Furthermore, Our analysis can be applied to any van der Waals heterobilayer system with weak interlayer hybridization and allows for arbitrary lattice constant mismatch and relative orientation.
Impurity effects on local electronic structure of twisted bilayer cuprates
It has been recently proposed that a fully gapped d-wave superconducting phase emerges in twisted bilayer cuprates when twist angle is in the vicinity of 45°. By considering a twisted bilayer model with a moiré cell, we study the local electronic structure around a single-site impurity and impurity dimer in both topologically trivial and nontrivial superconducting states. Here we find that the topological nontrivial superconducting phase is robust against single scalar-potential impurity scattering but is very sensitive to the interlayer impurity dimer scattering, which avoids and enforces the creation of the subgap quasiparticle state. This behavior is intrinsic to the bilayer structure and topological nature, in contrast to the other types of impurities that inevitably bring subgap quasiparticle states. We expect that our results can contribute to the understanding of electronic properties of twisted bilayer cuprates systems.
Stacking Domain Wall Magnons in Twisted van der Waals Magnets
In this report using bilayer CrI 3 as an example, we demonstrate that stacking domain walls in van der Waals magnets can host one-dimensional (1D) magnon channels, which have lower energies than bulk magnons. Interestingly, some magnon channels are hidden in magnetically homogeneous background and can only be inferred with the knowledge of stacking domain walls. Compared to 1D magnons confined in magnetic domain walls, 1D magnons in stacking domain walls are more stable against external perturbations. We show that the relaxed moiré superlattices of small-angle twisted bilayer CrI 3 is a natural realization of stacking domain walls and host interconnected moiré magnon network. Our Letter reveals the importance of stacking domain walls in understanding magnetic properties of van der Waals magnets and extends the scope of stacking engineering to magnetic dynamics.
Fragile Topology and Flat-Band Superconductivity in the Strong-Coupling Regime
In flat bands, superconductivity can lead to surprising transport effects. The superfluid “mobility”, in the form of the superfluid weight D s , does not draw from the curvature of the band but has a purely band-geometric origin. In a mean-field description, a nonzero Chern number or fragile topology sets a lower bound for D s , which, via the Berezinskii-Kosterlitz-Thouless mechanism, might explain the relatively high superconducting transition temperature measured in magic-angle twisted bilayer graphene (MATBG). For fragile topology, relevant for the bilayer system, the fate of this bound for finite temperature and beyond the mean-field approximation remained, however, unclear. Here, we numerically use exact Monte Carlo simulations to study an attractive Hubbard model in flat bands with topological properties akin to those of MATBG. We find a superconducting phase transition with a critical temperature that scales linearly with the interaction strength. Then, we investigate the robustness of the superconducting state to the addition of trivial bands that may or may not trivialize the fragile topology. Our results substantiate the validity of the topological bound beyond the mean-field regime and further stress the importance of fragile topology for flat-band superconductivity.
Topological Charge Pumping in Excitonic Insulators
Here, we show that in excitonic insulators with s-wave electron-hole pairing, an applied electric field (either pulsed or static) can induce a p-wave component to the order parameter, and further drive it to rotate in the s + ip plane, realizing a Thouless charge pump. In one dimension, each cycle of rotation pumps exactly two electrons across the sample. Higher dimensional systems can be viewed as a stack of one-dimensional chains in momentum space in which each chain crossing the Fermi surface contributes a channel of charge pumping. Physics beyond the adiabatic limit, including in particular dissipative effects is discussed.
Second-Order Josephson Effect in Excitonic Insulators
We show that in electron-hole bilayers with excitonic orders arising from conduction and valence bands formed by atomic orbitals that have different parities, nonzero interlayer tunneling leads to a second-order Josephson effect. This means the interlayer electrical current is related to the phase of the excitonic order parameter as J = J c sin 2θ instead of J = J c sin θ and that the system has two degenerate ground states at θ = 0, π that can be switched by an interlayer voltage pulse. When generalized to a three dimensional stack of alternating electron-hole planes or a two dimensional stack of chains, the ac Josephson effect implies that electric field pulses perpendicular to the layers and chains can steer the order parameter phase between the two degenerate ground states, making these devices ultrafast memories. Here, the order parameter steering also applies to the excitonic insulator candidate Ta 2 NiSe 5 .