Search

Dense lattices of photonic crystals can serve as artificial materials, with light propagation in these structures described by effective material parameters that surpass the capabilities of natural materials. In this study, we introduce a metamaterial that supports quadrupole magnetization, a characteristic rarely observed in existing structures. We experimentally demonstrate a magnetic quadrupole metamaterial associated with Mie-resonance-excited stop bands below the Bragg band. Additionally, we develop a theoretical model that addresses both dispersion and boundary conditions within this framework. Using a Fabry-Perot resonator as a case study, we validate our model and reveal that the quadrupole metamaterial can exhibit a markedly different reflection/transmission spectrum, including zero reflection at normal incidence. Our findings underscore the practical potential for both experimental and theoretical investigations of metamaterials that extend beyond the dipole approximation.

We investigate the interplay of generalized global symmetries in 2+1 dimensions by introducing a lattice model that couples a $\mathbb{Z}_N$ clock model to a $\mathbb{Z}_N$ gauge theory via a topological interaction. This coupling binds the charges of one symmetry to the disorder operators of the other, and when these composite objects condense, they give rise to emergent generalized symmetries with mixed 't Hooft anomalies. These anomalies result in phases with ordinary symmetry breaking, topological order, and symmetry-protected topological (SPT) order, where the different types of order are not independent but intimately related. We further explore the gapped boundary states of these exotic phases and develop theories for phase transitions between them. Additionally, we extend our lattice model to incorporate a non-invertible global symmetry, which can be spontaneously broken, leading to domain walls with non-trivial fusion rules., Comment: 84 pages, 13 figures

A charge density wave (CDW) is a phase of matter characterized by a periodic modulation of the valence electron density accompanied by a distortion of the lattice structure. The microscopic details of CDW formation are closely tied to the dynamic charge susceptibility, $\chi(q,\omega)$, which describes the behavior of electronic collective modes. Despite decades of extensive study, the behavior of $\chi(q,\omega)$ in the vicinity of a CDW transition has never been measured with high energy resolution ($\sim$meV). Here, we investigate the canonical CDW transition in ErTe$_3$ using momentum-resolved electron energy loss spectroscopy (M-EELS), a technique uniquely sensitive to valence band charge excitations. Unlike phonons in these materials, which undergo conventional softening due to the Kohn anomaly at the CDW wavevector, the electronic excitations display purely relaxational dynamics that are well described by a diffusive model. The diffusivity peaks around 250 K, just below the critical temperature. Additionally, we report, for the first time, a divergence in the real part of $\chi(q,\omega)$ in the static limit ($\omega \rightarrow 0$), a phenomenon predicted to characterize CDWs since the 1970s. These results highlight the importance of energy- and momentum-resolved measurements of electronic susceptibility and demonstrate the power of M-EELS as a versatile probe of charge dynamics in materials.

The Pauli exclusion principle combined with interactions between fermions is a unifying basic mechanism that can give rise to quantum phases with spin order in diverse physical systems. Transition-metal ferromagnets, with isotropic ordering respecting crystallographic rotation symmetries and with a net magnetization, are a relatively common manifestation of this mechanism, leading to numerous practical applications, e.g., in spintronic information technologies. In contrast, superfluid $^3$He has been a unique and fragile manifestation, in which the spin-ordered phase is anisotropic, breaking the real-space rotation symmetries, and has zero net magnetization. The recently discovered altermagnets share the spin-ordered anisotropic zero-magnetization nature of superfluid $^3$He. Yet, altermagnets appear to be even more abundant than ferromagnets, can be robust, and are projected to offer superior scalability for spintronics compared to ferromagnets. Our Perspective revisits the decades of research of the spin-ordered anisotropic zero-magnetization phases including, besides superfluid $^3$He, also theoretically conceived counterparts in nematic electronic liquid-crystal phases. While all sharing the same extraordinary character of symmetry breaking, we highlight the distinctions in microscopic physics which set altermagnets apart and enable their robust and abundant material realizations., Comment: 12 pages, 4 figures

The interplay between quantum and thermal fluctuations in the presence of quenched random disorder is a long-standing open theoretical problem which has been made more urgent by advances in modern experimental techniques. The fragility of charge density wave order to impurities makes this problem of particular interest in understanding a host of real materials, including the cuprate high-temperature superconductors. To address this question, we consider the quantum version of an exactly solvable classical model of two-dimensional randomly pinned incommensurate charge density waves first introduced by us in a recent work, and use the large-$N$ technique to obtain the phase diagram and order parameter correlations. Our theory considers quantum and thermal fluctuations and disorder on equal footing by accounting for all effects non-perturbatively, which reveals a novel crossover between under-damped and over-damped dynamics of the fluctuations of the charge density wave order parameter., Comment: 49 pages, 7 figures

Thanks to their record high refractive index and giant optical anisotropy, van der Waals (vdW) materials have accelerated the development of nanophotonics. However, traditional high refractive index materials, such as titanium dioxide (TiO2), still dominate in the most important visible range. This is due to the current lack of transparent vdW materials across the entire visible spectrum. In this context, we propose that germanium disulfide (GeS2) could offer a significant breakthrough. With its high refractive index, negligible losses, and biaxial optical anisotropy across the whole visible range, GeS2 has the potential to complement TiO2 and close the application gap of vdW materials in the visible spectrum. The addition of GeS2 could have a profound impact on the design of van der Waals nanophotonic circuits for any operation wavelength from ultraviolet to infrared, emphasizing the significance of the potential impact of GeS2 on the field of nanophotonics., Comment: 15 pages, 4 figures

Predicting molecular impact on cellular function is a core challenge in therapeutic design. Phenomic experiments, designed to capture cellular morphology, utilize microscopy based techniques and demonstrate a high throughput solution for uncovering molecular impact on the cell. In this work, we learn a joint latent space between molecular structures and microscopy phenomic experiments, aligning paired samples with contrastive learning. Specifically, we study the problem ofContrastive PhenoMolecular Retrieval, which consists of zero-shot molecular structure identification conditioned on phenomic experiments. We assess challenges in multi-modal learning of phenomics and molecular modalities such as experimental batch effect, inactive molecule perturbations, and encoding perturbation concentration. We demonstrate improved multi-modal learner retrieval through (1) a uni-modal pre-trained phenomics model, (2) a novel inter sample similarity aware loss, and (3) models conditioned on a representation of molecular concentration. Following this recipe, we propose MolPhenix, a molecular phenomics model. MolPhenix leverages a pre-trained phenomics model to demonstrate significant performance gains across perturbation concentrations, molecular scaffolds, and activity thresholds. In particular, we demonstrate an 8.1x improvement in zero shot molecular retrieval of active molecules over the previous state-of-the-art, reaching 77.33% in top-1% accuracy. These results open the door for machine learning to be applied in virtual phenomics screening, which can significantly benefit drug discovery applications.

The long-sought pair density wave (PDW) is an exotic phase of matter in which charge density wave (CDW) order is intertwined with the amplitude or phase of coexisting, superconducting order \cite{Berg2009,Berg2009b}. Originally predicted to exist in copper-oxides, circumstantial evidence for PDW order now exists in a variety of materials. Recently, scanning tunneling microscopy (STM) studies have reported evidence for a three-component charge density wave (CDW) at the surface of the heavy-fermion superconductor, UTe$_2$, persisting below its superconducting transition temperature. Here, we use hard x-ray diffraction measurements on crystals of UTe$_2$ at $T = 1.9$ K and $12$ K to search for a bulk signature of this CDW. Using STM measurements as a constraint, we calculate the expected locations of CDW superlattice peaks, and sweep a large volume of reciprocal space in search of a signature. We failed to find any evidence for a CDW near any of the expected superlattice positions in many Brillouin zones. We estimate an upper bound on the CDW lattice distortion of $u_{max} \lesssim 4 \times 10^{-3} \mathrm{\AA}$. Our results suggest that the CDW observed in STM is either purely electronic, somehow lacking a signature in the structural lattice, or is restricted to the material surface.

Pair density waves (PDWs) are a inhomogeneous superconducting states whose Cooper pairs possess a finite momentum resulting in a oscillatory gap in space, even in the absence of an external magnetic field. There is growing evidence for the existence of PDW superconducting order in many strongly correlated materials, particularly in the cuprate superconductors and in several other different types of systems. A feature of the PDW state is that inherently it has a CDW as a composite order associated with it. Here we study the structure of the electronic topological defects of the PDW, paying special attention to the half-vortex and its electronic structure that can be detected in STM experiments. We discuss tell-tale signatures of the defects in violations of inversion symmetry, in the excitation spectrum and their spectral functions in the presence of topological defects. We discuss the ``Fermi surface'' topology of Bogoliubov quasiparticle of the PDW phases, and we briefly discuss the role of quasiparticle interference., Comment: A few editorial edits were updated as well as some acknowledgments. A couple minus signs were also missing in equation 3.13 and 3.14. Finally, the orthogonal matrix after equation B.5 needed to be updated. A few references were also updated or removed

Scaling deep learning models has been at the heart of recent revolutions in language modelling and image generation. Practitioners have observed a strong relationship between model size, dataset size, and performance. However, structure-based architectures such as Graph Neural Networks (GNNs) are yet to show the benefits of scale mainly due to the lower efficiency of sparse operations, large data requirements, and lack of clarity about the effectiveness of various architectures. We address this drawback of GNNs by studying their scaling behavior. Specifically, we analyze message-passing networks, graph Transformers, and hybrid architectures on the largest public collection of 2D molecular graphs. For the first time, we observe that GNNs benefit tremendously from the increasing scale of depth, width, number of molecules, number of labels, and the diversity in the pretraining datasets. We further demonstrate strong finetuning scaling behavior on 38 highly competitive downstream tasks, outclassing previous large models. This gives rise to MolGPS, a new graph foundation model that allows to navigate the chemical space, outperforming the previous state-of-the-arts on 26 out the 38 downstream tasks. We hope that our work paves the way for an era where foundational GNNs drive pharmaceutical drug discovery.

Geometric fluctuations of the density mode in a fractional quantum Hall (FQH) state can give rise to a nematic FQH phase, a topological state with a spontaneously broken rotational symmetry. While experiments on FQH states in the second Landau level have reported signatures of putative FQH nematics in anisotropic transport, a realistic model for this state has been lacking. We show that the standard model of particles in the lowest Landau level interacting via the Coulomb potential realizes the FQH nematic transition, which is reached by a progressive reduction of the strength of the shortest-range Haldane pseudopotential. Using exact diagonalization and variational wave functions, we demonstrate that the FQH nematic transition occurs when the system's neutral gap closes in the long-wavelength limit while the charge gap remains open. We confirm the symmetry-breaking nature of the transition by demonstrating the existence of a "circular moat" potential in the manifold of states with broken rotational symmetry, while its geometric character is revealed through the strong fluctuations of the nematic susceptibility and Hall viscosity., Comment: Four figures in main text with supplementary information

In this work we study the interaction between vortices and nematic domain walls within the framework of a Ginzburg Landau approach. The free energy of the system is written in terms of a complex order parameter characteristic of $s$-wave superconductivity and a real (Ising type) order parameter associated to nematicity. The interaction between both order parameters is described by a biquadratic and a trilinear derivative term. To study the effects of these interactions we solve the time-dependent dissipative Ginzburg Landau equations using a highly performant pseudospectral method by which we calculate the trajectories of a vortex that, for different coupling parameters, is either attracted or repelled by a wall, as well as of the wall dynamics. We show that despite its simplicity, this theory displays many phenomena observed experimentally in Fe-based superconductors. In particular we find that the sign of the biquadratic term determines the attractive (pining) or repulsive (antipining) character of the interaction, as observed in FeSe and BaFeCoAs compounds respectively. The trilinear term is responsible for the elliptical shape of vortex cores as well as for the orientation of the axes of the ellipses and vortex trajectories with respect to the axes of the structural lattice. For the case of pining, we show that the vortex core is well described by a heart-shaped structure in agreement with STM experiments performed in FeSe.

Grating couplers are widely used to couple waveguide modes with the far field. Their usefulness is determined not only by energy efficiency but also by additional supported functionality. In this paper, we demonstrate a plasmonic grating on a silicon nitride slab waveguide that couples both TE and TM waveguide modes with circularly polarized light in the far field. Specifically, we experimentally confirmed that circularly polarized light excites TE and TM modes propagating in opposite directions, and the direction is controlled by the handedness. The routing efficiency for normally incident light reaches up to 95%. The same structure operates in the outcoupling regime as well, demonstrating up to 97% degree of circular polarization, where the handedness is determined by the polarization and propagation direction of outcoupled modes. Our results pave the way for the realization of polarization-division multiplexers and demultiplexers, integrated circular polarization emitters, as well as detectors of the polarization state of the incident optical field.

In the face of rapidly accumulating genomic data, our understanding of the RNA regulatory code remains incomplete. Recent self-supervised methods in other domains have demonstrated the ability to learn rules underlying the data-generating process such as sentence structure in language. Inspired by this, we extend contrastive learning techniques to genomic data by utilizing functional similarities between sequences generated through alternative splicing and gene duplication. Our novel dataset and contrastive objective enable the learning of generalized RNA isoform representations. We validate their utility on downstream tasks such as RNA half-life and mean ribosome load prediction. Our pre-training strategy yields competitive results using linear probing on both tasks, along with up to a two-fold increase in Pearson correlation in low-data conditions. Importantly, our exploration of the learned latent space reveals that our contrastive objective yields semantically meaningful representations, underscoring its potential as a valuable initialization technique for RNA property prediction.

Fundamental studies of the interaction of chiral light with chiral matter are important for the development of techniques that allow handedness-selective optical detection of chiral organic molecules. One approach to achieve this goal is the creation of a Fabry-P\'erot cavity that supports eigenmodes with a desired electromagnetic handedness, which interacts differently with left and right molecular enantiomers. In this paper, we theoretically study chiral Fabry-P\'erot cavities with mirrors comprising one-dimensional photonic crystal slabs made of van der Waals As$_2$S$_3$, a material with one of the highest known in-plane anisotropy. By utilizing the anisotropy degree of freedom provided by As$_2$S$_3$, we design Fabry-P\'erot cavities with constitutional and configurational geometrical chiralities. We demonstrate that in cavities with constitutional chirality, electromagnetic modes of left or right handedness exist due to the chirality of both mirrors, often referred to as handedness preserving mirrors in the literature. At the same time, cavities with configurational chirality support modes of both handednesses due to chiral morphology of the entire structure, set by the twist angle between the optical axes of the upper and lower non-chiral anisotropic mirrors. The developed chiral Fabry-P\'erot cavities can be tuned to the technologically available distance between the mirrors by properly twisting them, making such systems a prospective platform for the coupling of chiral light with chiral matter., Comment: 33 pages, 9 figures

Parafermion zero modes can arise in hybrid structures composed of $\nu=1/m$ fractional quantum Hall edges proximitized with an s-wave superconductor. Here we consider parafermion and Cooper pair tunneling, and backscattering in a junction formed in such hybrid structures. We find that the $4\pi m$ periodicity due to parafermion-only tunneling reduces, in the presence of backscattering, to $4\pi$-periodic at zero temperature and $2\pi$-periodic at finite temperature unless the fermion parity is fixed. Nevertheless, a clear signature of parafermion tunneling remains in the shape of the current-phase relation., Comment: v1: 6 pages, 4 figures; v2: 8 pages, 5 figures, appendix added

The nature of the interplay between fluctuations and quenched random disorder is a long-standing open problem, particularly in systems with a continuous order parameter. This lack of a full theoretical treatment has been underscored by recent advances in experiments on charge density wave materials. To address this problem, we formulate an exactly solvable model of a two-dimensional randomly pinned incommensurate charge density wave, and use the large-$N$ technique to map out the phase diagram and order parameter correlations. Our approach captures the physics of the Berezinskii-Kosterlitz-Thouless phase transition in the clean limit at large $N$. We pay particular attention to the roles of thermal fluctuations and quenched random field disorder in destroying long-range order, finding a novel crossover between weakly- and strongly-disordered regimes., Comment: 40 pages, 4 figures

Topological defects are singularities in an ordered phase that can have a profound effect on phase transitions and serve as a window into the order parameter. In this work we use scanning tunneling microscopy to visualize the role of topological defects in the novel magnetic field induced disappearance of an intertwined charge density wave (CDW) in the heavy fermion superconductor, UTe2. By simultaneously imaging the amplitude and phase of the CDW order, we reveal pairs of topological defects with positive and negative phase winding. The pairs are directly correlated with a zero CDW amplitude and increase in number with increasing magnetic field. These observations can be captured by a Ginzburg Landau model of a uniform superconductor coexisting with a pair density wave. A magnetic field generates vortices of the superconducting and pair density wave order which can create topological defects in the CDW and induce the experimentally observed melting of the CDW at the upper critical field. Our work reveals the important role of magnetic field generated topological defects in the melting the CDW order parameter in UTe2 and provides support for the existence of a parent pair density wave order on the surface of UTe2., Comment: 29 Pages, includes manuscript and supplemental information, 4 Main figures and 8 Supplemental Figures

In this chapter I discuss the impact of concepts of Quantum Field Theory in modern Condensed Physics. Although the interplay between these two areas is certainly not new, the impact and mutual cross-fertilization has certainly grown enormously with time, and Quantum Field Theory has become a central conceptual tool in Condensed Matter Physics. In this chapter I cover how these ideas and tools have influenced our understanding of phase transitions, both classical and quantum, as well as topological phases of matter, and dualities., Comment: Revised version to appear as a chapter in the Encyclopedia of Condensed Matter Physics 2e. 126 pages, including 380 references

Multilayer stacks of twisted optical metasurfaces are considered as a prospective platform for chiral nanophotonic devices. Such structures are primarily used for the realization of circularly polarized light sources, artificial optical rotation, and circular dichroism. At the same time, the behavior of their hybrid photonic modes is strongly affected by the moir\'e-pattern of superimposed periodic constituents. In this work, we show that moir\'e-periodicity in bilayer dielectric photonic crystal slabs leads to an arise of unlimitedly narrow optical resonances, which are very sensitive to the relative twist and gap width between the sublayers. We demonstrate the structure providing twist-tuning of the hybrid mode wavelength in the range of 300--600~nm with quality factor varying from~$10^2$~up~to~$10^5$ correspondingly. The obtained results pave the wave for the utilization of moir\'e-assisted effects in multilayer photonic crystal slabs.

The intense interest in triplet superconductivity partly stems from theoretical predictions of exotic excitations such as non-abelian Majorana modes, chiral supercurrents, and half-quantum vortices. However, fundamentally new, and unexpected states may emerge when triplet superconductivity appears in a strongly correlated system. In this work we use scanning tunneling microscopy to reveal an unusual charge density wave (CDW) order in the heavy fermion triplet superconductor, UTe2. Our high-resolution maps reveal a multi-component incommensurate CDW whose intensity get weaker with increasing field, eventually disappearing at the superconducting critical field, Hc2. To explain the origin and phenomenology of this unusual CDW, we construct a Ginzburg-Landau theory for a uniform triplet superconductor coexisting with three triplet pair density wave (PDW) states. This theory gives rise to daughter CDWs which would be sensitive to magnetic field due to their origin in a triplet PDW state, and naturally explains our data. Our discovery of a CDW sensitive to magnetic fields and strongly intertwined with superconductivity, provides important new information for understanding the order parameter of UTe2 and uncovers the possible existence of a new kind of triplet PDW order which has not been previously explored., Comment: 12 pages, 4 figures, Supplementary information

A long-standing problem in the study of topological phases of matter has been to understand the types of fractional topological insulator (FTI) phases possible in 3+1 dimensions. Unlike ordinary topological insulators of free fermions, FTI phases are characterized by fractional $\Theta$-angles, long-range entanglement, and fractionalization. Starting from a simple family of $\mathbb{Z}_N$ lattice gauge theories due to Cardy and Rabinovici, we develop a class of FTI phases based on the physical mechanism of oblique confinement and the modern language of generalized global symmetries. We dub these phases oblique topological insulators. Oblique TIs arise when dyons -- bound states of electric charges and monopoles -- condense, leading to FTI phases characterized by topological order, emergent one-form symmetries, and gapped boundary states not realizable in 2+1-D alone. Based on the lattice gauge theory, we present continuum topological quantum field theories (TQFTs) for oblique TI phases involving fluctuating one-form and two-form gauge fields. We show explicitly that these TQFTs capture both the generalized global symmetries and topological orders seen in the lattice gauge theory. We also demonstrate that these theories exhibit a universal "generalized magnetoelectric effect" in the presence of two-form background gauge fields. Moreover, we characterize the possible boundary topological orders of oblique TIs, finding a new set of boundary states not studied previously for these kinds of TQFTs., Comment: 67 pages, 4 figures, 5 appendices. v2: new appendix, added references, and minor changes. v3: light revisions, additional discussion in Section 6 and Appendix E, and updated acknowledgements and references

In this work we explore the interplay between superconductivity and nematicity in the framework of a Ginzburg Landau theory with a nematic order parameter coupled to the superconductor order parameter, often used in the description of superconductivity of Fe based materials. In particular, we focus on the study of the vortex-vortex interaction in order to determine the way nematicity affects its attractive or repulsive character. To do so, we use a dynamical method based on the solutions of the Time Dependent Ginzburg Landau equations in a bulk superconductor. An important contribution of our work is the implementation of a pseudo-spectral method to solve the dynamics, known to be highly efficient and of very high order in comparison to the usual finite differences/elements methods. The coupling between the superconductor and the (real) nematic order parameters is represented by two terms in the free energy: a biquadratic term and a coupling of the nematic order parameter to the covariant derivatives of the superconductor order parameter. Our results show that there is a competing effect: while the former independently of its competitive or cooperative character generates an attractive vortex-vortex interaction, the latter always generates a repulsive interaction., Comment: 37 pages; 12 figures; corrected typos; added references

Plasmonic metasurfaces form a convenient platform for light manipulation at the nanoscale due to their specific localized surface plasmons. Nevertheless, despite the high degree of light localization in metals, their intrinsic Joule losses are often considered prevention from applications in high-quality dielectric structures. Here, we experimentally demonstrate that in some cases, the capabilities of plasmonic particles for light manipulation prevail over the negative impact of absorption. We show the lattice of plasmonic nanoparticles onto a dielectric waveguide that efficiently couples the light of both circular polarizations to guided modes propagating in opposite directions. We demonstrate 80% degree of circular polarization for the out-coupled emission of GaAs-waveguide-embedded quantum dots. The results allow us to consider the lattice as a circular-polarization-controlled grating coupler operating at normal incidence and make this structure prospective for further implementation as an efficient coupling interface for various integrated devices.

Online reviews are typically written by volunteers and, as a consequence, information about seller quality may be under-provided in digital marketplaces. We study the extent of this under-provision in a large-scale randomized experiment conducted by Airbnb. In this experiment, buyers are offered a coupon to review listings that have no prior reviews. The treatment induces additional reviews and these reviews tend to be more negative than reviews in the control group, consistent with selection bias in reviewing. Reviews induced by the treatment result in a temporary increase in transactions but these transactions are for fewer nights, on average. The effects on transactions and nights per transaction cancel out so that there is no detectable effect on total nights sold and revenue. Measures of transaction quality in the treatment group fall, suggesting that incentivized reviews do not improve matching. We show how market conditions and the design of the reputation system can explain our findings.

Understanding the nature of charge density waves (CDW) in cuprate superconductors has been complicated by material specific differences. A striking example is the opposite doping dependence of the CDW ordering wavevector in La-based and Y-based compounds, the two families where charge ordering is strongest and best characterized. Here we report a combined resonant soft X-ray scattering (RSXS) and neutron scattering study of charge and spin density waves in isotopically enriched La$_{1.8-x}$ Eu$_{0.2}$ Sr$_{x}$ CuO$_{4}$ over a range of doping $0.07 \leq x \leq 0.20$. For all dopings studied by RSXS, we find that the CDW amplitude is approximately temperature-independent and develops well above experimentally accessible temperatures. Surprisingly, the CDW ordering wavevector shows a non-monotonic temperature dependence, with a sudden change occurring at temperatures near the SDW onset temperature. We describe this behavior with a Landau-Ginzburg theory for an incommensurate CDW in a metallic system with a finite charge compressibility and CDW-SDW coupling. Our Landau-Ginzburg analysis suggests that the ordering wavevector at high temperatures decreases with increased doping. This behavior is opposite to the trend at low temperatures and highly reminiscent of the doping dependence seen in YBa$_2$ Cu$_3$ O$_{6+\delta}$ , suggesting a common origin of the CDW in hole-doped cuprate superconductors., Comment: Manuscript submitted to PNAS; ArXiv version formatted in Revtex4-2; Supplementary Material included as appendices; 20 pages, 14 figures (6 in the main text), and 41 references

We investigate the thermalization of Sachdev-Ye-Kitaev (SYK) models coupled via random interactions following quenches from the perspective of entanglement. Previous studies have shown that when a system of two SYK models coupled by random two-body terms is quenched from the thermofield double state with sufficiently low effective temperature, the R\'enyi entropies do not saturate to the expected thermal values in the large-$N$ limit. Using numerical large-$N$ methods, we first show that the R\'enyi entropies in a pair SYK models coupled by two-body terms can thermalize, if quenched from a state with sufficiently high effective temperature, and hence exhibit state-dependent thermalization. In contrast, SYK models coupled by single-body terms appear to always thermalize. We provide evidence that the subthermal behavior in the former system is likely a large-$N$ artifact by repeating the quench for finite $N$ and finding that the saturation value of the R\'enyi entropy extrapolates to the expected thermal value in the $N \to \infty$ limit. Finally, as a finer grained measure of thermalization, we compute the late-time spectral form factor of the reduced density matrix after the quench. While a single SYK dot exhibits perfect agreement with random matrix theory, both the quadratically and quartically coupled SYK models exhibit slight deviations., Comment: v1: 41 pages single column, 14 figures; v2: Additional references; v3: Minor revisions

In recent years twisted bi-layers of 2D materials became very popular in the field due to the possibility to totally change their electronic properties by simple rotation. At the same time, in the wide field of photonic crystals, this idea still remains almost untouched, and only some particular problems were considered. One of the reasons is the computational difficulty of the accurate consideration of Moir\'e superlattices that appear due to the superimposition of misaligned lattices. Indeed, the unit cell of the complex lattice is typically much larger than the original crystals and requires much more computational resources for the computations. Here, we propose a careful adaptation of the Fourier modal method in the form of the scattering matrices for the description of twisted 1D gratings' stacks. Our approach allows us to consider sublattices in close vicinity to each other and account for their interaction via the near-field. In the developed numerical scheme, we utilize the fact that each sublattice is only 1D-periodic and therefore simpler than the resulting 2D superlattice, as well as the fact that even a small gap between the lattices filters out high Fourier harmonics due to their evanescent origin. This accelerates the computations from 1 up to 3 and more orders of magnitude for typical structures depending on the number of harmonics. This paves the way for rigorous study of almost any photonic crystals of the proposed geometry and demonstration of specific Moir\'e-associated effects.