Realization of Photonic Charge-2 Dirac Point by Engineering Super-modes in Topological Superlattices

Quite recently a novel variety of unconventional fourfold linear band degeneracy points has been discovered in certain condensed-matter systems. Contrary to the standard 3-D Dirac monopoles, these quadruple points referred to as the charge-2 Dirac points are characterized by nonzero net topological charges, which can be exploited to delve into hitherto unknown realms of topological physics. Here, we report on the experimental realization of the charge-2 Dirac point by deliberately engineering hybrid topological states called super-modes in a 1-D optical superlattice system with two additional synthetic dimensions. Utilizing direct reflection and transmission measurements, we exhibit the existence of super-modes attributed to the synthetic charge-2 Dirac point, which has been achieved in the visible region for the first time. We also show the experimental approach to manipulating two spawned Weyl points that are identically charged in synthetic space. Moreover, topological end modes uniquely resulting from the charge-2 Dirac point can be delicately controlled within truncated superlattice samples, opening a pathway for us to rationally engineer local fields with intense enhancement.

E ver since the remarkable discovery that fermion-like energy excitations predicted by relativistic quantum field theories can emerge in electronic crystals whose band structures display linear band degeneracy points, a great deal of theoretical and experimental interest has been attracted in exploring such materials known as topological semimetals. The corresponding gapless semimetal phases are regarded as novel topological states, which open a new era in investigating condensed-matter physics. Substantial attention is engaged by Weyl points (WPs) 1-6 and 3-D Dirac points (DPs) 7,8 . WPs that reported actively in electronic systems are identified as synthetic magnetic monopoles in momentum space, carrying topological charges (Chern numbers) of ±1 and featured by "Fermi arc" surface states 6,9 . DPs can be viewed as two overlapping WPs with opposite topological charges, predicted and observed in crystals as well. However, it has been recently demonstrated that unconventional topological points appear in certain crystal structures [10][11][12][13][14][15][16][17][18] , which cannot be described in accordance with an emergent relativistic field theory. On such candidate is the charge-2 Dirac point (CDP), existing as a double-Weyl phonon in transition metal monosilicides 14,15 . It's generated by merging a pair of identically charged WPs, and hence possessing the topological charge of ±2. Consequently, CDPs are radically distinct from traditional DPs and can give rise to novel physical phenomena.
While topological semimetals found in nature exhibit exotic phases of matter, great process in understanding such band topology has also been impelled by the research on engineered systems. The core idea of engineering lattices is to create emergent band structures analogous to those formed in electronic crystals, which can be highly tunable and have fundamentally discriminative properties, providing us unprecedented opportunities of studying topological physics. Recent developments of experimental techniques have propelled ultracold atomic gases 19,20 , photonics 4,[21][22][23][24][25][26] , and acoustics [27][28][29][30] as promising systems to engineer WPs and DPs with novel emergent properties. Constructing complex 3-D structures with certain symmetry broken is perceived as the most common strategy 4,[21][22][23][24][27][28][29] , whereas another route to realize topological points is based on synthetic dimensions [30][31][32][33][34][35] . The initial motivation for employing synthetic dimensions was to explore fundamental physical effects in a space with a higher dimensionality via introducing controllable artificial dimension(s) in addition to the real spatial degree(s), in especial topological effects within systems beyond 3-D space. However, lately the interest of synthetic dimensions is fueled by the capacity to study topological features of 3-D degenerate points in 2-D (1-D) systems, dramatically simplifying experimental designs 19,26,30,34,36,37 . Either of the methods has been extensively exploited for WPs and DPs in the recent years. Nevertheless, as for the CDP, the only engineered system supporting it reported so far is made up of an acoustic metamaterial corresponding to a classical 3-D phononic crystal with a nonsymmorphic structure 38 . To our knowledge, in the visible regime, neither 3-D engineered systems nor 2-D (1-D) structures equipped with synthetic dimensions have been established to realize CDPs.
Here, we propose an experimentally feasible scheme to realize CDPs in a 1-D optical superlattice system with working frequencies lying in the visible region harnessing synthetic dimensions, and manipulate the spawned WPs with the same topological charges. To attain this, we start by designing suitable photonic modes interacting with each other to form a 1-D superlattice. Instructively, interfaces between distinct topological phases of matter host robust and exotic quantum states, the use of which acts as a strong driver of current research in condensed matter [39][40][41][42][43] . Hence, we stack together two kinds of photonic crystals (PCs) belonging to different class of topology to create such topological interface modes (TIMs), and on this basis topological states of photons associated with CDPs can be fully investigated under the introduction of synthetic space, facilitating the experimental realization, which is otherwise elusive at such frequencies. Furthermore, intriguing topological-protected end modes emerge at the termini of the truncated superlattice, guaranteed by the CDP with nonzero topological charge. More precisely, these end modes uniquely result from the bulk-edge correspondence 44 for each of the two WPs producing the CDP in synthetic space, which in turn could be tuned independently. Such topological end modes resemble surface states in Weyl semimetals 2,3,6 , holding great potential for applications in nonlinear optics 45 , quantum optics 46 , and lasers 47 owing to strongly enhanced localized fields.

Results
Design concept of the creation of CDPs. The starting point of our scheme is to construct a 1-D topological superlattice by use of TIMs existing at interface of two PCs with discriminative topological class. Specifically, our lattice consists of these two PCs stacked alternatively, in which each interface supports a TIM that hybridizes with each other to form a novel variety of artificial collective modes, resulting in a 1-D superlattice band structure where a single TIM serves as the photonic orbital. Similar heterostructures have previously been rendered for grapheme nanoribbons 41,42 and topological insulator superlattices 43 .
For CDPs to occur, we require two more dimensions added to the wave vector dimension provided by the existing 1-D superlattice. It's noticeable that the coupling of nearestneighbor TIMs, inclusive of both the magnitude and the sign, can be feasibly tuned by altering the repeated number of the PC's unit cell between adjacent interfaces. Moreover, the on-site resonance frequency of a TIM can be highly controllable if we put a defective unit with adjustable thickness at the interface. Therefore, the modulation of the coupling between adjacent TIMs and the on-site frequency of a single TIM is readily available, which allows us to parameterize these two variables and treats them as two artificial momentum dimensions. Through meticulous design, WPs can thus arise in such 3-D synthetic space owing to the hybrid modes designated as the super-modes, but the realization of CDPs begs for the overlap of two equivalently charged WPs. To this end, we exploit the polarization degree of light. The fact that TIMs response discriminately to transverse-magnetic (TM) and transverse-electric (TE) polarized light appends a so-called 'pseudospin' degree of freedom to the synthetic space, and the appearance of a CDP is finally achieved by merging a pair of WPs with the same topological charge but different pseudospins. Surprisingly, the CDP can conversely be split into two spawned WPs in synthetic space, whose trajectories are tunable via utilizing the pseudospin degree. Such procedure has never been revealed in practice prior to us, offering the evidence that our proposed artificial systems are used to not only explore topological excitations discovered before, but also navigate a way of studying novel phenomena. In particular, we design an applicable and smart strategy to detect the CDP and spawned WPs straightforwardly, which has never been reported before us.
The adjustability of a single TIM. We first provide a detailed introduction to the proposed structure holding a single TIM that is highly adjustable. As shown in Fig. 1a, it consists of two kinds of PCs (PC-p and PC-q), and a defective unit D. The unit cells of PC-p, PC-q, and the defective unit D are represented Here, M1 and M2 denote two dielectric materials and the subscripts stand for the thickness of associated layers. The stacking structure built up of repeated p(q)-type unit cells can thus be described as p 6 (q 6 ), in which the subscripts are employed to show the number of unit cells. Hence, we adopt p 6 Dq 6 to label the structure shown in the bottom of Fig. 1a, which can support a single TIM 48 as discussed in the Supplementary Note 1. Experimentally, such structure is fabricated with e-beam evaporation and we fabricate three samples with identical d 1 = 70 nm and d 2 = 79 nm, but setting d = 0 nm (sample I), d = 5 nm (sample II), and d = 10nm (sample III), respectively. Figure 1b exhibits the scanning electron microscope (SEM) picture of sample II with the highlighted p-type and q-type unit cells. The measured trans-m0ission spectra of the sample II (III) are given in Fig. 1c by black (magenta) circles under normal incident light (k x = 0 μm −1 ), where the common band gaps (bands) of the two PCs are highlighted as the white (gray) regions. It can be seen that sharp peaks inside the gap appear, which are attributed to the excitation of a TIM. To verify this, Fig. 1d exhibits the calculated spatial distribution of the electric field profile for the associated state of Sample II, from which we can see that such state decays rapidly away from the position of D-a distinctive signature of the TIM. We can also see that the TIM peak of sample III lies at the lower frequency than that of sample II in Fig. 1c, showing the ddependent feature of resonance frequency of TIMs. To make it clear, resonance frequencies of TIMs for these three samples are marked by red open circles in Fig. 1e, decreasing significantly with increasing value of d. This confirms the fact that modulating thickness of the defective unit provides us a feasible strategy to adjust on-site frequencies of TIMs.
Furthermore, the TIM exists for both TM and TE polarizations, but they are degenerate in the case of normal incidence (k x = 0 μm −1 ). To lift such degeneracy, we need to use oblique The refractive indices of HfO 2 and SiO 2 are 2 and 1.46, respectively. c Measured transmission spectra of Sample II under normal incidence, oblique incidence of transverse-magnetic (TM) waves, and oblique incidence of transverse-electric (TE) waves are shown by black, red, and blue circles, respectively. The oblique incident angle is 30°. The magenta circles show measured transmission spectrum of Sample III under normal incidence. incident light with k x ≠ 0 μm −1 , and hence we measure the transmission spectra of sample II under the TM (TE) polarized light with oblique incident angles (k x = 6 μm −1 ), as shown by red (TM) and blue (TE) open circles in Fig. 1c, confirming the removal of the degeneracy. Figure 1f shows the measured TIM frequencies as a function of k x for Sample II for both polarizations. It is noticeable that the splitting between TM (red) and TE (blue) polarized TIMs increases monotonically with the increment of k x , matching well with the calculated results (solid lines). Therefore, such splitting between TM and TE polarized TIMs affords us another degree of freedom to manipulate the TIMs.
Coupling signs and magnitudes between two TIMs. Next, we investigate the effects of q-type and p-type PCs as the coupling channel between two TIMs. As shown schematically in Fig. 2a, each structure is made up of stacking PCs (p-type and q-type) separated by two defective units. With the same notation in Fig. 1a, two designs in Fig. 2a can then be denoted as p 6 Dq N Dp 6 and q 6 Dp N Dp 6 , respectively, where the subscript N is the number of unit cells of associated PC. The overlapping of two individual TIMs with the same frequency ω 0 gives rise to two hybridized TIMs, one symmetric mode (S) at ω S and one antisymmetric one (AS) at ω AS . Here, the symmetric types are defined by the symmetry of the electric field which uses the center of p N (q N ) as the reference point. In the following, we demonstrate that the normalized coupling strength J ≡ (ω S − ω AS )/2ω 0 (see Methods), which describes the coupling amplitudes and signs, is directly controlled by N for either q-type or p-type PC in the middle.
We start by considering the q-type PC as the coupling elements, namely p 6 Dq N Dp 6 with normal incidence (k x = 0 μm −1 ). Figure 2b shows the transmission spectra of p 6 Dq 6 Dp 6 and p 6 Dq 7 Dp 6 (black circles). The parameters of both p-type and q-type unit cells are identical to those in Fig. 1b except that d of the defective unit is 0 nm. For both samples, we see two transmission peaks owing to two hybridized TIMs (S and AS). For the sample p 6 Dq 6 Dp 6 , ω S < ω AS such that J < 0. While for the sample p 6 Dq 7 Dp 6 , ω S > ω AS such that J > 0. In Fig. 2c, the N-dependence of J for q-type PC case, which is extracted from experimental data, is shown by blue downward-pointing triangles. We can see that |J| possesses a negative association with N, and the sign of J totally relies on the parity of N. For the samples p 6 Dq N Dp 6 , if Nis odd, J < 0, otherwise J > 0. This is because the accumulated phase for each unit cell is π.
We then explore J for the p-type PC as the coupling channel case, namely the samples q 6 Dp N Dp 6 . The corresponding results are present as red upward-pointing triangles in Fig. 2c, in which J has the same magnitude as that of p 6 Dq N Dp 6 but with opposite signs for a given N. In order to confirm these results, we also plot the results calculated by a commercial finite element simulation software (COMSOL) in Fig. 2c, which shows good agreements with experimental ones. The sign of J is determined by the coupling between two TIMs. More details of two coupled TIMs are provided in the in the Supplementary Note 2.
Moreover, we investigate coupling effects for nonzero in-plane wave vector (k x < 0 μm −1 ). Figure 2b also depicts the measured transmission spectra of the two samples p 6 Dq 6 Dp 6 and p 6 Dq 7 Dp 6 for TM (red open circles) and TE (blue open circles) polarizations with k x = 6 μm −1 , indicating the polarization-dependent characteristic of the hybridized TIMs. To get a further step, we plot J as a function of k x in Fig. 2d, where p 6 (q 6 ) and p 7 (q 7 ) are employed as representatives of the even and odd N cases. The result clearly shows that the sign of J changes in the same way as that of k x = 0 μm −1 . However, given a fixed p N (q N ), the magnitude of J due to the TM mode has discriminated variation tendency compared with J of the TE one as k x increases, which almost remains the same for the former while decreases significantly for the latter. The calculated results achieved by a COMSOL via the optic module are shown in Fig. 2d by lines, which used to confirm the experimental results.

Realization of the CDPs and spawned WPs in synthetic space.
According to the previous analysis, the eigenfrequency of a single TIM is readily controllable, and the coupling (including signs and magnitudes) between two adjacent TIMs is highly tunable. All of these are sufficient for us to construct a 1-D topological superlattice analogous to a dimerized atomic chain, where we regard TIMs as photonic orbitals. The hybridization of them forms hybrid orbitals, which are referred to as super-modes here. As a consequence, we deliberately design an optical superlattice to create a periodic sequence of TIMs, which is built up of alternating structures of p m , D A , q n , and D B , as illustrated in Fig. 3a. The i-th unit dimer with two sublattice sites A i and B i is defined as [p m/2 D A q n D B p m/2 ] marked by the magenta dashed rectangle in Fig. 3a. In this notation, the subscripts A and B denote two different defective units with their respective thickness d A and d B , and the subscripts m(n) labels the number of unit cells of the q (p)-type PC. We then express the on-site resonance frequencies of two adjacent TIMs as ω A,s = ω s + Δ s and ω B,s = ω s − Δ s , where s = ↑↓ denoting two polarizations, ω s ≡ (ω A,s + ω B,s )/2, and Δ s ≡ (ω A,s − ω B,s )/2. Here, Δ s refers to a staggered on-site frequency offset regarding ω s as a reference value. As shown in Fig. 1, the values of Δ s are determined by the difference (d A − d B ). The coupling PCs q n (p m ) directly determine the intra(inter)dimer coupling strength, denoted as J 1,s (J 2,s ). What's more, a remarkable feature of our superlattice system is the adjustability of the coupling sign, since that m and n are simultaneously odd or even numbers leads to J 1,s J 2,s > 0, otherwise J 1,s J 2,s < 0. Taking this into account, we introduce an additional parameter g as g sgn J 2;s J 1;s , and utilize J s ≡ (−gJ 1,s − J 2,s )/2g and δ s ≡ (−gJ 1,s + J 2,s )/2g for further discussion. As shown in Fig. 2, the values of δ s are determined by (m − n). Accordingly, the Hamiltonian for the super-modes formed by multiple TIMs can be written as an effective dimerized model.
Here, a y i;s (b y i;s ) and a i,s (b i,s ) are the creation and annihilation operators of the TIM lying on A i (B i ) site of the chain, respectively. Since D A and D B have negligible effect on the coupling strength, δ s and Δ s can be treated as independent parameters (see details in the Supplementary Note 3). If we merely restrict ourselves to the case of zero in-plane wave vector (k x = 0 μm −1 ), this Hamiltonian represents a novel 1-D Su-Schrieffer-Heeger (Rice-Mele) chain with d A = d B (d A ≠ d B ) and hence Δ s = 0(Δ s ≠ 0), of which band structures and topological properties such as topological end states are analyzed detailedly in the Supplementary Note 4. Based on the fact that the degeneracy of TM and TE polarized TIMs lifts when k x ≠ 0 μm −1 , ω s is a function of k x such that ω s (k x ) = ω 0 + τ s (k x ), where ω 0 ≡ ω s (0) denotes the eigenfrequency of the k x = 0 μm −1 case and τ s (k x ) refers to the frequency shift compared with ω 0 . The Hamiltonian (Eq. 1) can thus be transformed into the Bloch momentum space, and expressed in the pseudospin up (s = ↑) (TM) and down (s = ↓) (TE) representation as Here, we introduce ðd s Þ in which Λ is the length of the unit dimer, ξ serves as the Bloch wave vector in the z direction and σ stands for Pauli matrices. Thereby, the eigenvalue of the Hamiltonian (Eq. 2) could be figured out, denoted asω ω À ω 0 with With respect to the special case at k x = 0 μm −1 , TM and TE polarized super-modes are degenerate sinceτ ¼ 0 andd " ¼d # . Hence we introduce the parameters δ ≡ δ ↑ = δ ↓ (k x = 0 μm −1 ) and Δ ≡ Δ ↑ = Δ ↓ (k x = 0μm −1 ), together with the original Bloch wave vector ξ, to form a 3-D synthetic space (δ,ξ,Δ). The Hamiltonian then can be transformed into H(δ,ξ,Δ) = (H ↑ , 0; 0, H ↓ ), in which H s¼";# d s¼";# Á σ. As a result, the associated four bands cross at the degenerate point (δ,ξ,Δ) = (0, 0, 0). To characterize this degenerate point, we expand the two-by-two Hamiltonian H s around it: where v δx,s = −2, v ξy,s = ΛJ s , and v Δz,s = 1 (see "seehods"). The above Hamiltonian exhibits a standard Weyl Hamiltonian form, and thus the band crossing point for either TM or TE supermodes can be regarded as a WP in the synthetic space. An important characteristic of a WP is the capacity to carry a topological charge, which corresponds to its chirality c s (= ±1). Such a four-fold cone can be detected unambiguously in experiment. We start by making five samples with structural parameters (m,n) = (4,4), (4,6), (4,8), (6,4), and (8,4), respectively, featured by d A = d B such that Δ = 0. We then measure the transmission spectra under normal incidence for these five samples to obtain locations ofω. Figure 3d presents the transmission spectra as a function of Δf (utilizing ω 0 as the reference) for the sample (m,n) = (6,4), where the black dashed lines emphasize the super-modes band edges. In Fig. 3b, c, we employ black squares to mark locations of such band edges, which almost lie on the crossing lines indicating a great agreement with the theory. The locations ofω for other four samples are plotted as well by different dots in Fig. 3b, c, all of which are situated at the crossing lines and thus exhibit the characteristic of linear crossing for the Dirac point, matching with the theory quite well (See details in the Supplementary Note 6). Moreover, we fabricate another two samples with (m,n) = (4,6) and (6,4), characterized by d A ≠ d B and hence Δ ≠ 0. The results gotten from experimental data are also shown in Fig. 3b, well-located at the cone's surface. Consequently, the experimental results support our theory of the CDP, and hence the realization of a CDP in the visible light range is achieved.
When k x ≠ 0 μm −1 , the degeneracy of TM and TE super-modes is removed, lending toτ ≠ 0 andd " ≠d # . Therefore, nonzero k x will split the CDP at k x = 0 μm −1 into two WPs of TM and TE polarized super-modes, respectively. The solid surface in Fig. 3e shows such two WPs in the δ − Δ space at ξ = 0 with k x = 6 μm −1 , and the dispersion in the (ξ,Δ) = (0,0) plane are present in Fig. 3f by red(blue) solid lines for TM(TE) modes. To demonstrate it, we measure the transmission spectra under oblique incident light of the seven samples defined in Fig. 3b, c. We choose to show the transmission spectra of the sample (m,n) = (6,4) in Fig. 3d, and the band edges of the super-modes for TM and TE polarizations are highlighted by red and blue vertical lines. We see that the band edges red-shifted (blue-shifted) for TE (TM) polarization, which agrees with theoretical results in Fig. 3e. We further mark the locations of associatedω of all these samples in Fig. 3e, f with red (blue) color for the TM (TE) super-modes. The consistency between the theory and experiments indubitably corroborates our idea that the CDP is separated into two WPs in the synthetic space with the frequency splitting equals 2τ resulting from the nonzero k x . As a result,τ can be treated as the effective Zeeman term, which increases with the enhancement of the "magnetic field" (that is, the increase of k x ). Note that varying k x has no effect on c s of both WPs, so c ↑ = c ↓ as those of k x = 0 μm −1 (See details in the Supplementary Note 5).
Topological end modes in truncated optical superlattices. In contrast to the conventional 3-D DPs, which carry no net topological charge and thus are lack of topological surface states, the CDPs arising in our system characterized by Chern numbers equal to ±2 imply the existence of intriguing topological end modes. Such end modes originate from each of the two WPs guaranteed by the bulk-edge correspondence, known as one of the most significant experimental properties of WPs. In our synthetic space, these topologically protected modes can be separated into two groups due to TM and TE polarized supermodes, respectively, each of which is supposed to be engineered independently under different polarizations. In what follows we demonstrate their existence in our optical superlattice system when truncated in space. The configurations could be generalized as p s [p m/2 D A q n D B p m/2 ] 5 p s , which are composed of 5 unit dimes with an open condition along z-axis and are extended by s additional unit cells of at both termini to avoid interacting with external environment. Figure 4a sketches a specific structure with (m,n,s) = (4,6,5) and d A ≠ d B such that δ < 0 and Δ ≠ 0, meeting the condition of supporting two nondegenerate end states (See "Methods"). The substrate is made from SiO 2 at the bottom of the structure. We describe the incident light from the front (bottom) as F (B). Topological end modes should come in pairs regardless of the value of Δ, but it is Δ that determines locations of these two end states in synthetic space, as experimentally demonstrated in Fig. 4b, c. To make it clear, in Fig. 4d we depict eigenfrequency surfaces of WPs and corresponding topological ends modes in the δ − Δ space with k x = 6 μm −1 for both TM and TE polarized super-modes. In Fig. 4d, the end modes on purple sheets are located at the front side of truncated chains with ω ¼ ±τ þ Δ and can only be excited by F, whereas those on orange sheets are localized at the end of the chains withω ¼ ±τ À Δ excited only by B, in which the first plus (minus) sign applies to TM (TE) polarized end modes. Notably, the end states belonging to the intersections of the two sheets connect to the WPs, sharing the same mathematical origin as that of the Fermiarc surface states 1,3,6 in Weyl semimetals. They are plotted by magenta dotted lines, with Δ = 0 andω ¼ ±τ in synthetic space, and hence they can be excited by either F or B. It is, however, indispensable to underscore here that the Fermi-arc links two WPs in a periodically arranged system while the Fermi-arc-like end modes in our system connect a WP to the boundary of synthetic space explained by the existence of the net topological charge. Figure 4b presents the measured reflection spectra of the sample illustrated in Fig. 4a for TM and TE polarizations excited by both F and B with an oblique incident angle of 30°, where we observe four dips inside the TIM gaps. These modes are labelled as 1-4, whose distributions of the electric field norm are exhibited in Fig. 4a. It is worth noting that, for both TM and TE polarized super-modes, the modes excited by F (1 and 3) and B (2 and 4) are located at the front and bottom of the sample, respectively. We then investigate another sample with (m,n,s) = (4,6,5) and d A = d B , characterized by δ < 0 and Δ = 0. Its measured reflection spectra achieved from F with an incident angle of 30°for TM and TE polarizations are exhibited in Fig. 4c, where each dip in the gap of super-modes is attributed to two degenerate end modes localized at both termini. In Fig. 4b, c, transparent gray regions correspond to the common bulk band gaps of PC-p and PC-q, and the gray regions with extra inclined downward and upward lines stands for the bands of TM and TE polarized super-modes, respectively. Due to the fact that the introduction of the nonzero k x leads to a removal of the degeneracy between TM and TE polarized end modes, each mode divides into two states with the splitting increasing rapidly as k x increases, as revealed in Fig. 4e for the sample described in Fig. 4a. In addition, Fig. 4f, g provide projected dispersion cones within different Δ planes. The eigenfrequency surfaces of topological end modes are projected as straight dashed line and the locations of end modes mentioned above are plotted in corresponding planes. Therefore, the great conformance between calculations and experiments further verifies our argument of topological end modes based on the established optical superlattice system.

Discussion
As a burgeoning field of topological physics, the study of topological photonics has captured huge attention in recent years 49 . A wide range of photonic systems have been devoted to this field, such as waveguide arrays 50,51 , coupled silicon ring resonators 52 , and polariton superstructures 53 , achieving remarkable accomplishments in various branches. In particular, the importance of 3-D gapless states such as Weyl and multi-Weyl points is a strong driver of current research in topological photonics. The concept of synthetic dimensions, though initially introduced to explore the higher dimensional physics by parametric coupling between internal modes or by dynamically scanning over the parameter space, has been extended and developing rapidly in the realm of Weyl or Weyl-related physics 26,27,30,34,36,37 . Relevant photonic systems include photonic crystals 26 , 1-D circuit-QED lattices 36 , and 2-D ring resonator lattices 37 . The essential phenomena owing to these gapless phases such as bulk-edge correspondence should be viewed in synthetic space accordingly, but still can reflect the topological characters (e.g., Chern number) of these points, just like what we have done with synthetic dimensions in this work.
Thanks to the availability and adjustability of our 1-D superlattice system, we can investigate fundamental topological features of the CDP-a novel kind of multi-Weyl points-in the photonic context by the aid of synthetic dimensions. We demonstrate the highly tunable on-site resonance frequency of each TIM and the controlled periodic coupling of nearestneighbor TIMs within our superlattices. The TIMs play the role of photonic orbitals, and their hybridizations form topological super-modes, whose band structures can be ingeniously engineered to create CDPs in synthetic space with the pseudospin degree originating from the polarized property of light. It is, for the first time, to realize CDPs in the visible region. Without the help of synthetic dimensions, as well as the utilization of pseudospins which fundamentally change the system's behaviors, the creation of CDPs is more demanding, possible only in the infrared range restrained by obstacles to the fabrication of complex structures. Furthermore, the CDP can be artificially split into two spawned WPs with the same Chern number by introducing nonzero horizontal wave vector that removes the degeneracy between TM and TE polarized super-modes. Such amazing process has not been obtained in previous studies, thus opening a new frontier to explore emergent phenomena of topological physics. In addition, the approach of experimental detection we render here is facile and obtainable by measuring transmission and reflection spectra to examine band structures of super-modes unambiguously. It is noteworthy that the bulk-edge correspondence displays itself as topological end modes exclusive to the CDP, residing at boundaries of truncated superlattices, which can be manipulated with ease and hence be applied for local field enhancement in various realms [45][46][47] .
Although our system was engineered to explore CDPs in the optical region at first, it offers a versatile approach to investigating other prevailing topological physics. Notably, the CDPs in our system are, in essence, a sort of secondary topological phases. It arises due to the fact that each TIM itself is of topological origin, and hence the CDPs could be treated as a consequence of coupling among multiple TIMs. Nowadays, the secondary topological effects arouse great interest and curiosity, and various systems have been devoted to this novel field 54-56 . Our work thus presents the core concepts similar to these systems, and exhibits intriguing secondary topological signatures, as shown in Fig. 4 where topological end states lying in the tight-binding gaps of TIMs. Another hot topic on topological photonics is the exploration of non-Hermitian effects 51,57-64 , which highlights itself as the appearance of exceptional points 57,59 , rings 51 , and surfaces 64 . By harnessing absorptive losses, our system has the potential to control imaginary parts of the TIMs' eigenfrequencies (See details in Supplementary Note 7), and then can be exploited to realize exceptional points and associated non-Hermitian effects. In addition, we may research on nonlinear topological phenomena by utilizing materials with intensity-dependent refractive indexes, which are currently a hotspot and can hold great promise for applications 65,66 .

Methods
Tight-binding analysis of coupling effects. In this section we provide the approach we used to obtain normalized coupling strength between nearestneighbor TIMs. Since the TIMs are highly localized modes that decay rapidly into bulks, we can apply the tight-binding method to analyze such coupling effects. For the cases illustrated in Fig. 2a, the Hamiltonian can be written as the following matrix: where ω 0 is the eigenvalue of either TIM and t corresponds to the coupling term.
Diagonalizing the Hamiltonian matrix (Eq. 4), we achieved two eigenstates with eigenfrequencies and wave functions given by: Here, ϕ 1 and ϕ 2 are the wave functions of the TIMs trapped by two identical defective units. φ S j i and φ AS j i are the emergent symmetric (S) and antisymmetric (AS) coupling modes as described in the text. As a result, the normalized coupling strength can be written by definition as: To calculate the value of J as a function of N for both p-type and q-type PCs, as well as the in-plane wave vector k x for both TM and TE polarizations, a series of numerical simulations were performed by means of COMSOL Multiphysics via the optic module. The results are shown in Fig. 2c, d, and the distributions of electric field profiles for some specific configurations are exhibited in Supplementary Fig. 1.
In the experimental measurements, we harnessed the transmission spectra of each sample under normal or oblique incidence, where two transmission peaks attribute to two coupling TIMs (S and AS), and hence found the value of J. The outcomes are present in Fig. 2, matching well with the tight-binding analysis. During this procedure we treated k x as a constant with certain incident angle, since the range of working frequency is so small that k x is little influenced by the frequency variation. Strictly speaking, this treatment could be another reason (additional to the fabrication imperfection and experimental errors) that the measurements are slightly deviated from the simulations on frequencies, but it has minimal impact on the value of J. Therefore, error bars used in Fig. 2c, d only result from the broadening transmission peaks owing to experimental uncertainties.
Eigenfrequency surfaces of topological end modes. The bulk-edge correspondence of spawned WPs formed in our synthetic space is guaranteed by the topologically nontrivial property of the truncated superlattice. Such kind of superlattice must meet the condition that the absolute value of the coupling strength within a unit dimer is less than that between two adjacent dimers, no matter what value Δ takes. Consequently, if both m and n are even numbers such that J s < 0, δ must be negative and the corresponding eigenfrequency surfaces of topological end modes are plotted in Fig. 4d. However, if m and n are simultaneously odd, J s > 0 and δ must be positive, making the eigenfrequency surfaces of topological end modes become the reflection of those shown in Fig. 4d in the plane δ = 0.

Data availability
The data that support the findings of this study are available from the corresponding author on reasonable request.