Bilayer Honeycomb Research Articles

3D printed honeycomb bilayer structures based on polylactic acid as lightweight microwave absorbing materials

Abstract3D printed honeycomb structures constituted by neat polylactic acid (PLA) and conductive PLA (PLAc) with different cell sizes and thicknesses were manufactured through 3D printing technology based on material extrusion (MEX). These structures were arranged into bi‐layer honeycomb structures (BHS) for evaluation of the microwave absorbing properties at the X‐band (8.2–12.4 GHz) and Ku‐band (12.4–18 GHz) frequency ranges. The effects of cell size, sample thickness, layer thickness ratio and the nature of the top layer on the electromagnetic attenuation performance of BHS samples were investigated in terms of reflection loss (RL). The components geometric characteristics and the arrangement of the layers exerted great influence on the RL values and effective absorption bandwidth (EAB). Overall the bilayer structures with 5 mm total thickness presented the best EM response in terms of EAB and minimum RL, where the component with the larger cell size (7.7 mm), PLAc (1 mm) as the top layer and PLA (4 mm) as the bottom layer exhibited the best results with RL values of −41 dB and EAB of 9.5 GHz. The BHS system with 7.7 mm cell size was more efficient in Ku band frequency range, whereas that with 3.5 mm cell size gave better response in the X‐band frequency range. Therefore, RL values and EAB can be adjustable by modifying the building strategy of the printed parts, such as thickness of the sample, cell size and thickness distribution of the conductive and non‐conductive layers. The findings in this work show the importance of the building strategy to fabricate components with controllable BHS and improved properties for microwave absorption application.

Effect of interlayer exchange coupling interaction on topological phase of a bilayer honeycomb Heisenberg ferromagnet

Layered magnetic topological materials are material systems that exhibit both magnetic ordering and topological properties in their smallest two-dimensional units. Studying these systems may lead to the observation of new physical properties and phenomena, which has attracted considerable attention from researchers. The effect of interlayer exchange coupling interactions on bilayer honeycomb Heisenberg ferromagnets with interlayer coupled topological phase is investigated by using linear spin wave theory. The influence of introducing two additional types of interactions, i.e. interlayer exchange coupling interaction and interlayer easy-axis anisotropy interaction, on the topological phase transition are also explored in this work. By calculating the magnon dispersion relations at various interlayer exchange coupling interaction intensities, it is found that the band gaps of high energy band and low energy band both close and reopen at the Dirac points when the system reaches the critical value of interlayer exchange coupling interaction. In magnon systems, such physical phenomena typically relate to topological phase transitions. When calculating the Berry curvature and Chern numbers for the bands in the aforementioned process, it is found that the sign of the Berry curvature reverses and the Chern numbers change when the critical value of interlayer exchange coupling interaction strength is reached, confirming that a topological phase transition occurs indeed. Introducing two other types of interlayer exchange coupling interactions in this process can lead various novel topological phases to occur in the system. The enhancement of interlayer easy-axis anisotropy interactions is likely to impede the topological phase transitions occurring in the system. We find that a major distinction between bilayer honeycomb ferromagnets and their single-layer counterparts lies in the fact that during a topological phase transition, the sign of the magnon thermal Hall coefficient does not change; on the contrary, abrupt shift in the thermal Hall coefficient curve occurs which can be seen as an indicator of topological phase transition of bilayer honeycomb ferromagnets, and is also reflected in the change in magnon Nernst coefficient. The research results of this work can provide theoretical support for developing novel spintronic devices with enhanced information transmission capabilities by using bilayer honeycomb ferromagnetic materials, and can also provide theoretical reference for studing other bilayer ferromagnetic systems.

Magnetization Process in Bilayer Honeycomb Spin Lattice

The magnetization process of the AA-stacked bilayer honeycomb spin lattice in the out-of-plane applied field is examined in the framework of the Ising model and mean field approximation. Competition between different kinds (antiferromagnetic: AF, or ferromagnetic: FM) of intra-layer and inter-layer spin exchange couplings could lead to the first order magnetization process, characterized by sharp jumps in magnetization curves occurring in critical external fields (spin-flop and spin- flip fields). There is only the spin-flip field at very low temperature and its magnitude depends not only on the intra-layer frustration level but also exchange coupling between spin layers. During the magnetization process, the initial AF bilayer honeycomb spin lattice may undergo ferrimagnetic or week-ferromagnetic states before reaching full ferromagnetic saturation state by the spin-flip field. The results of this work can be applied for an explanation of the spin-flip phenomenon registered in the AF bilayer CrI3.

The Ising bilayer honeycomb lattice: A cluster mean-field study

The Ising bilayer honeycomb lattice: A cluster mean-field study

Sliding ferroelectricity in bilayer honeycomb structures: A first-principles study

The sliding ferroelectricity is engineered by artificially stacking the nonpolar two-dimensional (2D) materials, which greatly broadens the 2D ferroelectrics from very few candidate materials to a large family of 2D materials. However, the electric polarizations are generally small due to the weak van der Waals interlayer interaction. The search for 2D sliding ferroelectrics with large polarization presents an ongoing challenge. Here we systematically investigate the sliding ferroelectricity in the bilayer honeycomb structures of $\mathrm{B}X$ $(X=\mathrm{P}, \mathrm{As}, \mathrm{Sb})$, $Y\mathrm{N}$ $(Y=\mathrm{Al}, \mathrm{Ga}, \mathrm{In})$, and $Z\mathrm{C}$ $(Z=\mathrm{Si}, \mathrm{Ge}, \mathrm{Sn})$ based on first-principles calculations. It is shown that the electric polarization decreases with the increase of the interlayer distance, and increases with the difference in electronegativity of the two constituent elements. Such dependence is further corroborated by a simple model. It is interesting to see that GeC can be an ideal sliding ferroelectric material with high polarization and energetically favorable polar stacking. Our results reveal the key factors in determining the electric polarization, which could facilitate the search and design of 2D sliding ferroelectrics with large out-of-plane polarization.

Magnon corner states in twisted bilayer honeycomb magnets

Search for higher-order topological insulators, characterized by topologically protected gapless boundary states of codimension higher than one, in bosonic systems has attracted growing interest. Here, we establish twisted bilayer honeycomb magnets as a new platform for hosting second-order topological magnon insulators (SOTMIs) without fine-tuning. We employ a simple, minimal Heisenberg spin model to describe misaligned bilayer sheets of honeycomb ferromagnetic magnets with a large commensurate twist angle. We found that the higher-order topology in this bilayer system shows a significant dependence on the interlayer exchange coupling. The SOTMI, featuring topologically protected magnon corner states, appears for ferromagnetic interlayer couplings, while the twisted bilayer exhibits a nodal phase in the case of antiferromagnetic interlayer coupling.

Frustrated ferromagnetic transition in AB-stacked honeycomb bilayer

In two-dimensional (2D) ferromagnets, anisotropy is essential for the magnetic ordering as dictated by the Mermin-Wagner theorem. But when competing anisotropies are present, the phase transition becomes nontrivial. Here, utilizing highly sensitive susceptometry of scanning superconducting quantum interference device microscopy, we probe the spin correlations of ABC-stacked CrBr3 under zero magnetic field. We identify a plateau feature in susceptibility above the critical temperature (TC) in thick samples. It signifies a crossover regime induced by the competition between easy-plane intralayer exchange anisotropy versus uniaxial interlayer anisotropy. The evolution of the critical behavior from the bulk to 2D shows that the competition between the anisotropies is magnified in the reduced dimension. It leads to a strongly frustrated ferromagnetic transition in the bilayer with fluctuation on the order of TC, which is distinct from both the monolayer and the bulk. Our observation demonstrates unconventional 2D critical behavior on a honeycomb lattice.

Orbital dependent band degeneracy and edge states in single layer and AA bilayer honeycomb lattice systems with p orbital degeneracy

We have investigated band structure and edge states of single and AA bilayer honeycomb lattice systems with p orbital degeneracy which would be relevant for relatively heavy metals such as Pb and Bi. The p electrons, which are triply degenerate in the atomic limit, form quadruple degeneracy at K point in a single layer honeycomb lattice. The spin-orbit interaction splits the quartet at K point into two singlets and one doublet which are mixtures of p x , p y and p z orbitals. In the AA bilayer honeycomb lattice, the Dirac node at K point is deformed into circular line node around it. With the intermediate interlayer coupling, there is a degenerate point along Γ-K with p z character. A p x /p y type degenerate point exists under the strong interlayer coupling. This degenerate point is robust against the spin-orbit interaction and the antiferromagnetic spin arrangement only along the intralayer bonds. Without SOI, localized edge states with p x /p y or p z character are created at the zigzag edges both in the single and AA bilayer systems. In the AA bilayer system with SOI, only the p x /p y edge states exhibit the Dirac-like dispersion.

Dispersion-tunable photonic topological waveguides

Dispersion-tunable photonic topological waveguides have recently attracted much attention, due to their promising applications on topological devices with tunable operational frequencies. Since dispersions of topological waveguides traverse the whole bandgaps of bulk structures, tuning the dispersions (especially the bandwidths) requires changing the whole bulk of corresponding photonic topological insulators. A previously reported material-modification approach provided a parallel tuning on such numerous lattices; however, the increased material loss deteriorated transmissions of the topological waveguide. Here, a parallel tuning approach on structures is theoretically proposed and demonstrated, which spawns dispersion-tunable photonic topological waveguides without increasing material loss. Based on the bilayer honeycomb model, a topological valley waveguide by utilizing bilayer designer plasmonic structures is constructed, accomplished with dispersion tunings by altering interlayer distance. Experimental results validate the theoretical model and display a 61%-relative-tuning range of frequency, with a tunable relative bandwidth up to 16%. This approach may promise applications in tunable topological lasers, robust delay lines, and intelligent photonic devices.

Frustrated spin-[formula omitted] Heisenberg magnet on an [formula omitted]-stacked honeycomb bilayer: High-order study of the collinear magnetic phases of the [formula omitted]–[formula omitted]–[formula omitted] model

The zero-temperature phase diagram of the frustrated spin-12J1–J2–J1⊥ Heisenberg magnet on an AA-stacked honeycomb bilayer lattice is studied using the coupled cluster method implemented to very high orders. On each monolayer the spins interact via nearest-neighbor (NN) and frustrating next-nearest-neighbor isotropic antiferromagnetic Heisenberg interactions with respective strength parameters J1>0 and J2≡κJ1>0. The two layers are coupled such that NN interlayer pairs of spins also interact via a similar isotropic Heisenberg interaction of strength J1⊥≡δJ1, which may be of either sign. In particular, we locate with high accuracy the complete phase boundaries in the κ–δ half-plane with κ>0 of the two quasiclassical collinear antiferromagnetic phases with Néel or Néel-II magnetic order in each monolayer, and the interlayer NN pairs of spins either aligned (for δ<0) or anti-aligned (for δ>0) to one another. Compared to the two-sublattice Néel order, in which all NN intralayer pairs of spins are antiparallel to one another, the four-sublattice Néel-II order is characterized by NN intralayer pairs of spins on the honeycomb lattice being antiparallel to one another along zigzag (or sawtooth) chains in a specified direction from among the three equivalent honeycomb-lattice directions, and parallel to one another for the corresponding interchain pairs.

Higher-order Topological States in Acoustic Twisted Moiré Superlattices

Twisted moir\'e superlattices (TMSs) are emergent materials with exotic physical properties. Among their properties, higher-order topology is seldom realized and investigated in experiments. Here, we report on the experimental observation of a class of bilayer higher-order topological states in TMSs. We create the physical realization of acoustic TMSs using moir\'e twisting and ultrastrong interlayer couplings in bilayer honeycomb lattices of coupled acoustic cavities. Through such a design we reach an unexplored regime in TMSs that is not available in graphene TMSs and other solid-state TMS systems. We reveal that the ultrastrong interlayer couplings lead to a large acoustic band gap with unique higher-order band topology, which is characterized by unprecedented topological indices and layer-hybridized corner states. The higher-order topological edge and corner states are observed via acoustic pump-probe measurements, which show consistency with theory and simulations. Our study paves the way toward higher-order topological phenomena in TMSs.

Metallic and Deconfined Quantum Criticality in Dirac Systems

Motivated by the physics of spin-orbital liquids, we study a model of interacting Dirac fermions on a bilayer honeycomb lattice at half filling, featuring an explicit global SO(3)×U(1) symmetry. Using large-scale auxiliary-field quantum MonteCarlo (QMC) simulations, we locate two zero-temperature phase transitions as function of increasing interaction strength. First, we observe a continuous transition from the weakly interacting semimetal to a different semimetallic phase in which the SO(3) symmetry is spontaneously broken and where two out of three Dirac cones acquire a mass gap. The associated quantum critical point can be understood in terms of a Gross-Neveu-SO(3) theory. Second, we subsequently observe a transition toward an insulating phase in which the SO(3) symmetry is restored and the U(1) symmetry is spontaneously broken. While strongly first order at the mean-field level, the QMC data are consistent with a direct and continuous transition. It is thus a candidate for a new type of deconfined quantum critical point that features gapless fermionic degrees of freedom.

Theoretical realization of rich magnon topology by symmetry-breaking in honeycomb bilayer ferromagnets

We reveal the rich magnon topology in honeycomb bilayer ferromagnets (HBF) induced by the combined effect of interlayer exchange, Dzyaloshinskii-Moriya interaction (DMI), and electrostatic doping (ED). In particular, we present a systematic study of the Hamiltonian non-adiabatic evolution in the HBF parametric space, spanned by the symmetry-breaking terms (DMI and ED) and interlayer exchange. We determine the band closure manifolds which are found to divide the parametric space into six distinct regions, matched with five distinct topological phases and one topologically trivial phase. The characteristic Chern numbers and thermal Hall conductivities are calculated for the topological phases. Edge spectra, dictated by the bulk-edge correspondence, are also analyzed in the nanoribbon version of the model. Both bulk and edge spectra are found to be nonreciprocal as a consequence of ED and edge magnons are observed to counter propagate on opposite edges. The predicted results offer new insights on the manipulation of magnonic Chern numbers and magnon topological transport via experimentally tunable parameters.

Direct driving of electronic and phononic degrees of freedom in a honeycomb bilayer with infrared light

We study theoretically AB-stacked honeycomb bilayers driven by light in resonance with an infrared phonon within a tight-binding description. We characterize the phonon properties of honeycomb bilayers with group theory and construct an electronic time-dependent tight-binding model for the system following photoexcitation in resonance with an infrared phonon. We adopt an atomically adiabatic approximation [V. Mohanty and E. J. Heller, Proc. Natl. Acad. Sci. 116, 18316 (2019)] to describe classically vibrating nuclei, but obtain the Floquet quasienergy spectrum associated with the time-dependent model exactly. We introduce a general scheme to disentangle the complex low-frequency Floquet spectrum to elucidate the relevant Floquet bands. As a prototypical example, we consider bilayer graphene. We find that light in the low-frequency regime can induce a band gap in the quasienergy spectrum in the vicinity of the $K$ points even if it is linearly polarized, in contrast with the expectations within the Born-Oppenheimer approximation and the high-frequency regime. Finally, we analyze the diabaticity of the driven electron and driven phonon processes and find contrasting effects on the autocorrelation functions at the same driving frequency: Driven phonons preserve the character of the initial state while driven electrons exhibit strong deviations within a few drive cycles. The procedure outlined here can be applied to other materials to describe the combined effects of low-frequency light on phonons and electrons.

Effective field theory approach for the S=32 bilayer honeycomb antiferromagnet

The spin-$\frac{3}{2}$ Heisenberg antiferromagnet on the bilayer honeycomb lattice is a minimal model to describe the magnetic behavior of ${\mathrm{Bi}}_{3}{\mathrm{Mn}}_{4}{\mathrm{O}}_{12}({\mathrm{NO}}_{3})$. We study this model with frustrating interlayer second-neighbor couplings taking into account quantum and thermal fluctuations. We use a path integral formulation in terms of coherent states to describe the low-energy physics of the model. We show that for a particular point in the parameter space, close to the experimental estimated couplings, a continuum classical degeneracy is lifted by both quantum and thermal fluctuations, and a collinear state is then selected by an order by disorder mechanism. Our results provide a global perspective in the understanding of the experimental observations.

Topological Lifshitz transition and novel edge states induced by non-Abelian SU(2) gauge field on bilayer honeycomb lattice

We investigate the SU(2) gauge effects on bilayer honeycomb lattice thoroughly. We discover a topological Lifshitz transition induced by the non-Abelian gauge potential. Topological Lifshitz transitions are determined by topologies of Fermi surfaces in the momentum space. Fermi surface consists of N = 8 Dirac points at π-flux point instead of N = 4 in the trivial Abelian regimes. A local winding number is defined to classify the universality class of the gapless excitations. We also obtain the phase diagram of gauge fluxes by solving the secular equation. Furthermore, the novel edge states of biased bilayer nanoribbon with gauge fluxes are also investigated.

Designing flat-band tight-binding models with tunable multifold band touching points

Being dispersionless, flat bands on periodic lattices are solely characterized by their macroscopically degenerate eigenstates: compact localized states (CLSs) in real space and Bloch states in reciprocal space. Based on this property, this work presents a straightforward method to build flat-band tight-binding models with short-range hoppings \emph{on any periodic lattice}. The method consists in starting from a CLS and engineering families of Bloch Hamiltonians as quadratic (or linear) functions of the associated Bloch state. The resulting tight-binding models not only exhibit a flat band, but also multifold quadratic (or linear) band touching points (BTPs) whose number, location, and degeneracy can be controlled to a large extent. Quadratic flat-band models are ubiquitous: they can be built from any arbitrary CLS, on any lattice, in any dimension and with any number $N\geq2$ of bands. Linear flat-band models are rarer: they require $N\geq3$ and can only be built from CLSs that fulfill certain compatibility relations with the underlying lattice. Most flat-band models from the literature can be classified according to this scheme: Mielke's and Tasaki's models belong to the quadratic class, while the Lieb, dice and breathing Kagome models belong to the linear class. Many novel flat-band models are introduced, among which an $N=4$ bilayer honeycomb model with fourfold quadratic BTPs, an $N=5$ dice model with fivefold linear BTPs, and an $N=3$ Kagome model with BTPs that can be smoothly tuned from linear to quadratic.

Twisted superfluid and supersolid phases of triplons in bilayer honeycomb magnets

We demonstrate that low-lying triplon excitations in a bilayer Heisenberg antiferromagnet provide a promising avenue to realize magnetic analogs of twisted superfluid and supersolid phases that were recently reported for two-component ultracold atomic condensate in an optical lattice. Using a cluster Gutzwiller mean-field theory, we establish that Dzyaloshinskii-Moriya interactions (DMI), that are common in many quantum magnets, stabilize these phases in a magnetic system, in contrast to the pair hopping process that is necessary for ultracold atoms. The critical value of DMI for transition to the twisted superfluid and twisted supersolid phases depends on the strength of the (frustrated) interlayer interactions that can be tuned by applying external pressure on and / or shearing force between the layers. Furthermore, we show that the strength of DMI can be controllably varied by coupling to tailored circularly polarized light. Our results provide crucial guidance for the experimental search of twisted superfluid and supersolid phases of triplons in real quantum magnets.

Gross-Neveu-Heisenberg criticality from competing nematic and antiferromagnetic orders in bilayer graphene

We study the phase diagram of an effective model of competing nematic and antiferromagnetic orders of interacting electrons on the Bernal-stacked honeycomb bilayer, as relevant for bilayer graphene. In the noninteracting limit, the model features a semimetallic ground state with quadratic band touching points at the Fermi level. Taking the effects of short-range interactions into account, we demonstrate the existence of an extended region in the mean-field phase diagram characterized by coexisting nematic and antiferromagnetic orders. By means of a renormalization group approach, we reveal that the quantum phase transition from nematic to coexistent nematic-antiferromagnetic orders is continuous and characterized by emergent Lorentz symmetry. It falls into the $(2+1)$-dimensional relativistic Gross-Neveu-Heisenberg quantum universality class, which has recently been much investigated in the context of interacting Dirac systems in two spatial dimensions. The coexistence-to-antiferromagnetic transition, by contrast, turns out to be weakly first order as a consequence of the absence of continuous spatial rotational symmetry on the honeycomb bilayer. Implications for experiments in bilayer graphene are discussed.

Flat Bands in Magic-Angle Bilayer Photonic Crystals at Small Twists.

The new physics of magic-angle twisted bilayer graphene (TBG) motivated extensive studies of flat bands hosted by moiré superlattices in van der Waals structures, inspiring the investigations into their photonic counterparts with potential applications including Bose-Einstein condensation. However, correlation between photonic flat bands and bilayer photonic moiré systems remains unexplored, impeding further development of moiré photonics. In this work, we formulate a coupled-mode theory for low-angle twisted bilayer honeycomb photonic crystals as a close analogy of TBG, discovering magic-angle photonic flat bands with a non-Anderson-type localization. Moreover, the interlayer separation constitutes a convenient degree of freedom in tuning photonic moiré bands without high pressure. A phase diagram is constructed to correlate the twist angle and separation dependencies to the photonic magic angles. Our findings reveal a salient correspondence between fermionic and bosonic moiré systems and pave the avenue toward novel applications through advanced photonic band or state engineering.