Prominent Josephson tunneling between twisted single copper oxide planes of Bi2Sr2-xLaxCuO6+y

Josephson tunneling in twisted cuprate junctions provides a litmus test for the pairing symmetry, which is fundamental for understanding the microscopic mechanism of high temperature superconductivity. This issue is rekindled by experimental advances in van der Waals stacking and the proposal of an emergent d+id-wave. So far, all experiments have been carried out on Bi2Sr2CaCu2O8+x (Bi-2212) with double CuO2 planes but show controversial results. Here, we investigate junctions made of Bi2Sr2-xLaxCuO6+y (Bi-2201) with single CuO2 planes. Our on-site cold stacking technique ensures uncompromised crystalline quality and stoichiometry at the interface. Junctions with carefully calibrated twist angles around 45° show strong Josephson tunneling and conventional temperature dependence. Furthermore, we observe standard Fraunhofer diffraction patterns and integer Fiske steps in a junction with a twist angle of 45.0±0.2°. Together, these results pose strong constraints on the d or d+id-wave pairing and suggest an indispensable isotropic pairing component.

The recent application of van der Waals (vdW) stacking technique to high temperature (  ) cuprate superconductors [1][2][3][4][5] revives a dormant debate over the pairing symmetry of the order parameter [6][7][8].Phase sensitive experiments using inplane Josephson junctions showed evidence for d-wave pairing nearly three decades ago [9,10].However, experiments [11][12][13] of the c-axis Josephson effect, which is also able to determine the phase component of the pairing, revealed a large supercurrent between two cuprate crystals at a twist angle of 45°.The results were against the dwave pairing scenario [14,15] and suggested isotropic pairing instead.Since Josephson tunneling depends sensitively on the details of the interface, poor sample quality may account for the inconsistency [16].Of late, the vdW stacking of two-dimensional materials allows ever-expanding material pool [17,18] with improved interface quality [19][20][21] and precision control [22,23].These advancements benefit the study of high-  superconductivity too [1-4, 24, 25], because cuprates such as Bi-2212 are layered materials that can be exfoliated, even down to the monolayer level [1,5].Twisted Bi-2212 Josephson junctions are successfully realized by restacking the ultrathin flakes [6][7][8].They show atomically flat interfaces, closely matching the theoretically considered situation.This advancement allows a revisit to the observation of isotropic pairing in twisted cuprates [6].
At the twist angle of 45°, while the complete suppression of Josephson tunneling is demanded by pure d-wave symmetry, it is recently proposed that a d+id or d+is order parameter [27] can account for the finite Josephson tunneling.This emergent order parameter can open up a nodeless gap [27] and gives rise to a non-monotonic temperature dependence of the Josephson critical current   () or critical current density   () [28,29].The Fraunhofer diffraction pattern is expected to show doubling in frequency or halving of the period, owning to the co-tunneling of paired Cooper pairs with a net charge of 4e instead of 2e [26,28,29].Theory also predicts that the 45°-twist bilayer can produce Shapiro steps at half-integer positions under microwave irradiation.Although the proposed full gap was not observed in postannealed junctions [6], an experiment using cryogenically fabricated junctions exhibited non-monotonic   () and half-integer Shapiro steps [7], hinting at exotic pairing.Moreover, there is an apparent inconsistency in the angular dependences of the Josephson coupling strength     (  is the normal state resistance) in these experiments.Further experimental and theoretical investigations [30,31] are therefore necessary to elucidate these discrepancies.
So far, all the existing experiments on c-axis twisted cuprates have been carried out on a single type of cuprate superconductors-Bi-2212 [6][7][8][11][12][13].Whether the anomalous tunneling phenomena can be extended to other high temperature cuprate superconductors remains unaddressed.In this work, we choose Bi-2201 to fabricate twisted cuprates and study the Josephson tunneling (Fig. 1a).Although a series of scanning tunneling and photoemission experiments have been carried out on Bi-2201 [33][34][35], this material has not been employed in either the in-plane or out-of-plane Josephson tunneling experiments for phase-sensitive detection.Bi-2201 has the advantage that it avoids the unnecessary complexity brought by the bonding and antibonding effects of double CuO2 layers in Bi-2212 [36,37].While reduction from bilayer to single CuO2 layer enhances the two-dimensional effect [33], this situation better matches that considered by the theories [27].Technically, the doping level of ultrathin Bi-2201 can be more stable since it is governed by La substitution of Sr and interstitial O.The stoichiometry can be better preserved during the sample fabrication process [1,3,5].In the case of Bi-2212, the junctions often host a relatively short Josephson penetration depth λ  (0.2 μm) in comparison to the typical junction width (>1 μm).It gives rise to nonuniform phase distribution [38,39], which may further induce fractional Shapiro steps [40].To enhance uniformity, one can use the underdoped Bi-2212 with reduced   (<100 A/cm 2 ) such that λ  >1 μm [39,41].
Here, we fabricate multiple twisted Bi-2201 junctions with their twist angles calibrated to be close to 45°.We demonstrate their atomically flat interfaces and uncompromised stoichiometry by atomically resolved structural analysis and low temperature transport studies.All these junctions show prominent Josephson tunneling, indicating a systematic deviation from the expected behavior of d-wave pairing.We further reveal the conventional   () and the standard Fraunhofer pattern in one of the junctions with a calibrated twist angle of 45.0±0.2°.The absence of frequency-doubling as a function of temperature, together with the integer Fiske steps, strongly speaks against the existence of a subdominant Cooper pair cotunneling component-a crucial ingredient of d+id-wave pairing.Our study establishes that there exists pronounced isotropic pairing component in the 45°-twist cuprate bilayers.

Results
Device fabrication and structural analysis.We employ a carefully designed on-site cold stacking method to fabricate Bi-2201 junctions [3,[5][6][7].The whole fabrication process (Fig. 1b) assures that the cleaved flakes never leave the cold stage before protection-layer capping.As such, the crystallinity and stoichiometry of the bulk are well preserved in the artificial junctions.Figure 1b shows that the substrate holder and the cantilever with the polydimethylsiloxane (PDMS) stamp are cooled to a temperature below -50 °C.Two SiO2/Si substrates-one bare and another with prepatterned electrodes-are anchored on the same holder.We first cleave a Bi-2201 flake with a typical thickness of 300 nm on the bare substrate into two pieces.The flake left on the substrate, which is usually thicker, is subsequently rotated against its thinner counterpart on the stamp.We then use the stamp to pick up the rotated flake by making sure that it has a small overlap with the initially cleaved flake.The overlapping region therefore constitutes the twisted artificial junction.An ultraclean interface is guaranteed because there exists no glue residue from the tape or the stamp on the two freshly cleaved surfaces.The whole stack consisting of both top and bottom Bi-2201 is consequently transferred onto the other substrate with prepatterned electrodes.The release procedure is taken at about -30 °C.Finally, a flake of hexagonal boron nitride (h-BN) was deposited onto the junction region for protection.
Figure 1c demonstrates the high quality of our junctions.We obtain atomically resolved structures (after the transport measurements) in a large spatial span of about The supermodulations are usually attributed to the misfit between the CuO2/SrO planes and the BiO plane [46,47] and can locally modulate the energy gap [48].At the twisted interface, only one of the two BiO layers shows the persisting supermodulation.For sample S2, we observe that the Bi atoms in the bottom Bi-2201 bunch into pairs whereas the top Bi-2201 shows no such behavior.It indicates that the bottom part is viewed from [100], which is orthogonal to that of the bottom halves of S1, S3 and S4.At the interface, only one of the two BiO layers shows bunching.For sample S5, both the top and bottom Bi-2201 show bunching but the bunching in the top Bi-2201 is shifted by half a period in comparison to the bottom Bi-2201.Notably, the twist angles of samples S1 and S2 must be very close to 45°.This can be immediately inferred from the simultaneous atomic resolution of both the top and bottom parts of Bi-2201 with distinct crystalline orientations.A more precise determination by using the Kikuchi patterns (details in Supplementary Note 2 and Supplementary Fig. 1, 2) [6] yields a deviation (Δ) from 45° as small as 0.2° for S1 and 1.0°, 1.5°, 3.0° for S2, S3 and S4.For simplicity, we use  = 45°−  as the twist angles for these samples throughout the paper, since  = 45°+  is equivalent.
Figure 1d shows the normalized intensity along the c-axis for samples S1 to S5.One sees that the BiO and SrO planes (red arrows) at the junction interface possess comparable integrated intensities to those of the bulk (black arrows).More importantly, the two rotated CuO2 planes exhibit the overlapping signal intensities of Cu atoms as those CuO2 planes far away from the interface (Supplementary Figure 3), attesting to the uncompromised quality of the two superconducting layers involved in the Josephson tunneling process.We note that the distance between the two nearestneighboring BiO layers ( Bi−Bi ) at around  = 0 is larger in samples S1 to S4 than that of sample S5.The corresponding distance between the two CuO2 planes ( Cu−Cu ) across the interface also increases by about 5% from the value in S5 to that in S1-S4 (Supplementary Figure 4).
Junction resistance and tunneling.We now study the transport properties in the same batch of samples (S1-S4) with their crystalline structures well clarified.All junctions exhibit a single superconducting transition (Fig. 2a1-a4 and Supplementary Fig. 5).For S3, we further demonstrate that the transition temperature of the junction is exactly the same as that measured from the individual flakes.These results demonstrate that the doping concentration is preserved throughout the complete stack of Bi-2201 flakes.Figure 2b1-b4 presents the tunneling I-V characteristics.These samples possess a single branch of supercurrent before transitioning to the normal state.In Supplementary Fig. 6, we show the I-V characteristics for two samples in a wide current range.The single jump in one current direction clearly distinguishes from those hysteretic I-V behaviors caused by phase slip lines [49,50].We argue that the observed Josephson tunneling solely occurs between the two twisted CuO2 planes since they have the smallest tunneling area (numbers listed in Fig. 2b1-b4).The intrinsic junctions in the bottom Bi-2201, connected in series in the electrical circuit, have much larger tunneling area (3 to 10 times larger) such that their critical current is not reached in our measured regime (further support is given in Supplementary Notes 3 and 4 and Supplementary Fig. 8).For samples S3, S4 and OP1, the I-V characteristics show jumps from the zero-bias branch to the branch with finite resistances.There also exists prominent hysteresis in the two opposite sweeping directions.For samples S1 and S2 (Supplementary Fig. 5b), however, the switching seems continuous.These two samples are patterned by focused ion beam (FIB) into long strips (insets of Fig. 2a and Supplementary Fig. 5a) for the investigation of Fraunhofer patterns (to be discussed in the following section).This additional step seems to suppress the hysteresis.Nevertheless, we will show that data from these two samples also stand against the expectation from either d-or d+id-wave pairing.
Figures 2c1-c4 show that our samples exhibit the conventional temperature dependence of   : a rapid rise of   when  just decreases from   followed by saturation of   at low temperatures ( ≤ 0.5  ).This is similar to that prescribed by the Ambegaokar-Baratoff (AB) formula for the Josephson tunneling between two conventional superconductors [51].By contrast, the recent theories for d+id -wave pairing [28,29] predict a gradual increase of   from  =   to about 0.5  followed by a pronounced rise at lower  at twist angles close to 45°.They also predict strong non-monotonic   () around 30°.These features distinguish the higher-order term due to co-tunneling of Cooper pairs from the standard Josephson term.The conventional   () observed in samples S1-S4 and OP1 therefore suggest the dominance of the first-order Josephson tunneling process.
Fraunhofer pattern and Fiske steps.Further insight into the Josephson effect in twisted bilayer cuprates can be garnered from the magnetic field dependence.For sample S1, we measure the Josephson critical current as a function of magnetic field  ∥ applied parallel to the longer direction of the junction (Fig. 3a). Figure 3b displays the critical Josephson currents at zero bias from the  −  characteristics at different  ∥ , showing the well-defined Fraunhofer diffraction pattern.Figure 3c further shows the colored plot of / as a function of current and  ∥ at different temperatures.It captures both the Fraunhofer pattern as well as the AC Josephson effect.Solid curves in the figures are theoretical fits by using the standard formula: , where  ,0 is the zero-field critical current and  =  ∥  0 ⁄ with  0 being the oscillation period.This standard formula nicely reproduces the experimental features.By contrast, if the d+id-wave is present as a subdominant term, the co-tunneling of Cooper pairs (with 4e) will cause suppression of   ( ∥ ) at  0 /2, with  = 1, 3, 5 ⋯ (dashed curves in Fig. 3) [26][27][28][29].This frequency-doubling effect is expected to be most pronounced at low temperatures and gets smeared out when  →   .Such a phase transition is clearly absent in S1 here and in S2 in Supplementary Fig. 5. Instead, we observe only a subtle change, as indicated by the decreasing value of  0 with increasing .We defer further discussions on this evolution to the end of the paper.In addition, noises in the color plot become suddenly prominent upon increasing  ∥ to a certain value (marked by the red arrows).This onset field becomes larger with decreasing temperatures, reminiscent of the temperature dependence of critical magnetic field (Supplementary Fig. 11).We therefore attribute the noises to the entrance of Abrikosov vortices above the lower critical field.They are created by a small out-of-plane magnetic field due to slight misalignment between the magnetic field and the sample plane.
Notably, there exist additional regular branches in the Fraunhofer patterns (marked by horizontal arrows).Between these branches and the major lobe of the Fraunhofer pattern (following the solid curves) are regions hosting local minima of / (dark blue color).These regions are in fact related to the AC Josephson effect.Typically, the AC Josephson effect is demonstrated by current steps at discrete voltages-Shapiro steps-in a junction under microwave irradiation.Here, no external microwave source is necessary because we utilize the rectangular shape of our junction.Such a cavity possesses certain electromagnetic modes which can be excited under a bias voltage.
Interference between the electromagnetic fields resonating in the cavity and the AC Josephson effect of the junction gives rise to another type of current steps at discrete voltages-Fiske steps (inset of Fig. 4a) [52].Conventionally, Fiske steps occur at   =  0   , where  = 1,2,3 ⋯ ,  0 = ℎ/2 is the flux quantum and   is the fundamental frequency of the junction cavity mode.We replot the data in Fig. 3 in terms of dI/dV as a function of the junction voltage and  ∥ in Fig. 4b and c.A Fiske step in () (Fig. 4a) gives rise to a local peak in dI/dV (Fig. 4b).Ten horizontal stripes at non-zero voltages, reflecting ten Fiske steps, can be identified in Fig. 4c (marked alphabetically: a, b, c, d, e).The Fiske steps observed in sample S2 are shown in Supplementary Fig. 7.For each Fiske step, its intensity gets modulated by the magnetic field and breaks into several sections, as indicated by the black, gray and white triangles.The maximal intensity in each section drops monotonically with increasing  ∥ .The dotted curve in Fig. 4d reflects this monotonic behavior of the first step (stripea).Notably, the maximal intensity of the first section from each stripe also shows a monotonic decreasing trend with the voltage position (dashed curve in Fig. 4d).These features are all consistent with the observation in intrinsic junctions of Bi-2212 [52].It indicates that the Fiske steps follow a continuously increasing integer of  from stripe-a to stripe-e.Interestingly, the departure from stripe-a to zero bias ( 0 ) is two times the separation between the stripes (  for example).We therefore assign the stripes (a, b, c, d, e) as Fiske steps with  = 2, 3, 4, 5, 6.In this case, the first step ( = 1) is missing at small  ∥ .However, it is worth noting that the first step seems to reemerge at higher fields.This is illustrated in Fig. 4e where peaks in dI/dV occur at the expected positions of  = ±  0 /2 .This reentrant behavior was recently reported in Al/InAs junctions too, but in the topologically trivial regime [53].The missing of the first step and its reemergence may be caused by the smearing effect of the resistive branch in the I-V characteristics [54].The strong non-linearity of the I-V curve causes the highest dissipation at bias voltages close to zero.This can be seen from the darker color close to zero bias in Fig. 4b and c.With increasing  ∥ , the normal state branch becomes less resistive such that the first step can re-appear.We point out that the monotonic increase of the sequence of Fiske steps alone is compatible with Josephson tunneling of either purely 2e or purely 4e carriers.However, a dominant higher-order tunneling with 4e has a specific temperature dependence, which we have excluded by the conventional   () in Fig. 2. The remaining possibility is that there exists a mixture of a dominant 2e-tunneling with a small fraction of 4etunneling.In this case, the Fiske steps have two sets: one prominent set at integer steps due to 2e-tunneling and one weaker set at half integers due to 4e-tunneling.The intensity in / of the Fiske steps would therefore show non-monotonic evolution with consecutively increasing bias voltages, in sharp contrast to our data in Fig. 4.

Discussion
In Fig. 5, we include data points from 17 samples with similar   (within 3 K) (Supplementary Fig. 9b).All the samples were fabricated by using the same recipe.We outline the upper and lower bounds of the (2θ) dependence demanded by dwave pairing [14,15].They are based on the data points at 0° and 90°.Clearly, a cluster of seven data points around 45° fall outside this demarcation.The slight decrease in     at around 45°, if comparing the maximal values obtained at different angles, can be explained by the orbital effect between two conventional superconductors [55].
The observed increase of interlayer distance for S1-S4 may also play a role.Despite this variation, we highlight that these data points deviate strongly from the meV.This value is two times the gap of Bi-2201 [56] and greatly exceeds the upper bound expected even from the AB formula.It is also one order of magnitude larger than the typical experimental     values we obtained at 0° and 90° (gray band in Fig. 5).In terms of   , a similar calculation by using the data of S2 would yield   (0 ∘ ) =   (44.0 ∘ )/ (2 × 44.0 ∘ ) = 1.99 kA/cm 2 , exceeding the measured   in bulk Bi-2201 (Supplementary Note 4) by one order of magnitude and becoming on par with the record of Bi-2212 intrinsic Josephson junctions (Supplementary Fig. 10).We comment that   in our junctions is actually underestimated since the effective tunneling area may be smaller than the physical area [6].In general, a purely anisotropic d-wave pairing is incompatible with the experimental results, suggesting the existence of a pronounced isotropic pairing component in the twisted junctions.Such a prominent s-wave pairing could be a manifestation of symmetry breaking in this system.The innately present supermodulations may contribute because they break the C4 rotational symmetry required for the d-wave pairing.However, this effect was not revealed in other experiments until now.
We expand on the discussion of the temperature dependence of Josephson current   ().The recent theory on d+id-wave pairing proposed that the critical current hosts a Λ-shaped temperature dependence for junctions with twist angles from 18° to 36°: increases first as  increases from zero to about 0.4  and decreases at higher .Although we obtain predominantly junctions with conventional   () as shown in Fig. 2, we do observe non-monotonic   () occasionally.Supplementary Figure 12a provides such an example at  = 0°.Notably, at this twist angle the theory of d+id-wave pairing predicts standard AB-type behavior.The angular independence suggests a more trivial origin.To shed further insight into this issue, we study the influence of a small perpendicular magnetic field.Supplementary Figure 12b shows that the Josephson tunneling evolves to be a dome-like shape under a small magnetic field of 0.4 mT.Strikingly, for sample S3 with initially AB-type   () at zero field (Fig. 2c2), such a dome-like behavior can be induced by a small magnetic field (Supplementary Fig. 12c).Similar behaviors were reported in intrinsic Bi-2212 Josephson junctions [57].There, the increase in   () with increasing temperature was attributed to the depinning and realignment of the flux pancakes.We speculate that trapped fluxes can play a similar role in junctions even at zero magnetic field.
They can produce the non-monotonic behavior that resembles the predicted feature of a topologically non-trivial origin.
Finally, we speculate on the mechanism for the gradual temperature variation of  0 in the Fraunhofer patterns (Fig. 3c, further summarized in Supplementary Fig. 12d).
It indicates a crossover from the tunneling between two layered superconductors to that between anisotropic ones.The crossover occurs at a temperature  * that   () = ξ 0 /√1 − /  becomes comparable to the interlayer distance between the neighboring CuO2 planes.At  <  * , ξ  is shorter than the interlayer distance such that the twisted junction should be treated as a multilayered system.Such a layered system can be described by the phenomenological theory proposed by Lawrence and Doniach [58].If  < λ  (λ  is the Josephson penetration depth) and there exist an infinite number of layers with equal Josephson currents, the L-D model gives rise to the Fraunhofer pattern with a period [59]:  − =  0 /().Here  is the junction thickness and in Bi-2201 it should be equal to  Cu−Cu .At  * <  <   , ξ  is large enough that the top and bottom Bi-2201 can be individually treated as an anisotropic superconductor following the 3D Ginzburg-Landau (G-L) formula [7].The Fraunhofer period becomes:  − =  0 /[( Sr-Sr + 2λ  )], where λ  = λ 0 /√1 − /  is the caxis penetration depth and  Sr-Sr is the thickness of the tunnel barrier (between two SrO planes that sandwich the double BiO planes).This theoretical analysis can qualitatively reproduce the experimental trend (inset of Supplementary Fig. 12d).
In summary, we successfully realize twisted cuprate bilayers out of a previously unexplored compound-Bi-2201.We demonstrate that our samples host ultraclean interfaces and preserve the stoichiometry from the bulk-critical for investigating

Methods:
The optimally doped Bi2Sr1.63La0.37CuO6+x(Bi-2201) single crystal was grown by the traveling solvent floating zone method [60].The crystal possesses an onset superconducting transition temperature (  ) of 32 K.All junctions included in this work were fabricated from the same piece of crystal (Supplementary Fig. 9).The polydimethylsiloxane (PDMS) used for the van der Waals stacking was prepared by using Dow Corning Sylgard 184, mixing base and curing agent in a ratio of 10:1.It was then transferred onto a clean sapphire slide and quickly baked to 120 C on a hot plate for 10 min for adhesion.The sapphire slide ensured good thermal conduction at cryogenic temperatures.The SiO2/Si wafers (4 inch) with a 285 nm thick SiO2 were used as substrates.For prepatterning the substrates, the entire wafer was processed by employing the ultraviolet photo-lithography followed by metal deposition of Ti/Au (5/40 nm) with E-beam.Both the pristine and prepatterned wafers were then diced into rectangular pieces (3  10 mm 2 ) for fabrication of the Josephson junctions (schematically shown in Fig. 1b).Before loading into the glovebox for junction stacking, the unpatterened/patterned substrates were treated by O2 plasma (power 100 W/60 W, pressure 0.1/0.3mbar, time period 1/3 min), in order to enhance the bonding between Bi-2201 and SiO2 [5].The sample fabrication was carried out in a glovebox filled with Ar gas (H2O < 0.01 ppm, O2 < 0.01 ppm).
For samples S1 and S2, we employed FIB with gallium ions to pattern the narrow strips.The acceleration voltage was set at 30 kV and a small beam current of 2 pA was used to minimize the degradation effect.It was reported in the Bi-2212 intrinsic junctions that there existed a 30 nm thick amorphous layer after the FIB processing [61].This thickness is one order of magnitude smaller than the width of our samples (3.5 μm for S1 and 5 μm for S2) such that the influence of the amorphous layer should be negligible.However, in other strips with their widths smaller than 2 μm, we observed that   of the junction became suppressed to 10 K-15 K.It indicates that when the strip is too narrow FIB patterning may cause oxygen out-diffusion, presumably due to heating.We emphasize that this effect is minimized for samples S1 and S2, because their   values are the same as that of the bulk crystal.Nevertheless, we observe no hysteresis in the I-V characteristics of S1 and weaker hysteresis in that of S2, in contrast to those of samples without FIB processing.To understand this difference, we start from the RCSJ model [62].Whether the hysteresis is pronounced or not depends on the damping factor:   = 2    2 /ℏ, where  is the effective capacitance of the junction.A large   ≫ 1 corresponds to pronounced hysteresis.
The equation can be further written as: where   is the normal state resistivity,  is the effective length of the junction,  is the effective distance between the two superconducting planes that constitute a planar capacitor (essentially  ∼ ), and ϵ  is the relative dielectric constant.From the TEM images (Fig. 1), the atomic structure of the junctions is unchanged after FIB patterning such that  and  should not be affected.A less pronounced hysteresis, i.e., a smaller   , may therefore be attributed to the reduction of     ,   or ϵ  .It indicates that FIB patterning may influence our junctions in a subtle manner [63].
Future improvements can be gained by utilizing helium ions for FIB patterning [64,65].
The samples were wired to the chip carriers within an hour at ambient condition outside the glovebox.Transport measurements were carried out in two different closed-cycle helium-free systems equipped with superconducting magnets (base temperature below 1.6 K).Resistances were measured by the standard low-frequency lock-in technique with a typical ac current of 1 μA (13.3 Hz).For I-V measurements, we used a dc source to supply the current and a nanovoltmeter to measure the voltage in a four-terminal configuration.We measured dV/dI simultaneously by measuring the four-terminal ac voltage with a lock-in amplifier at an ac current of 500 nA (13.3 Hz).
To measure the Fraunhofer pattern, the samples were mounted on a home-built piezo-driven rotator inside the cryogenic system.The rotator has an angular precision of 0.006°.To align the magnetic field to the in-plane direction of the junction, we rotated the sample at 1 T and 14 K until the junction resistance showed a minimum.
We ramped up the magnetic field from zero in one direction and swept the tunneling current at each fixed magnetic field to measure both  and /.To obtain the data in the opposite magnetic field direction, we carried out a thermal cycling up to 50 K and back at zero field in order to expel trapped fluxes.
After the transport measurements, the junctions were investigated by transmission electron microscope (TEM) (S1-S5).The TEM samples were prepared by using FIB 70 nm.Images with equivalently high quality were taken from multiple locations along the direction parallel to the interface.The artificial interface, consisting of the lowest BiO layer of the top Bi-2201 and the uppermost BiO layer of the bottom Bi-2201 (indicated by arrows), patently distinguishes itself from other double BiO layers in either side of Bi-2201.For samples S1, S3 and S4, wavy undulations can be observed in the bottom halves of the images, suggesting that these sections are viewed from the [01 ̅ 0] crystalline direction perpendicular to the uniaxial supermodulations running along the [100] direction.The top parts show no wavy structures because they are viewed from the [11 ̅ 0] direction that is along the diagonal between [01 ̅ 0] and [100].
their intrinsic properties.Transport experiments reveal pronounced     values in these high-quality junctions at the twist angle of 45°, similar to that observed in twisted Bi-2212 junctions.These results show sharp contrast to the pure d-wave scenario and strongly suggest the presence of an isotropic pairing component in the twisted junctions.We carry out a systematic investigation to check the validity of recent theoretical proposals on the emergent d+id wave.We show that a Josephson junction at 45.0±0.2°possesses: 1) conventional temperature dependence of the critical current; 2) standard Fraunhofer pattern with no frequency doubling as a function of temperature; 3) integer Fiske steps.These results strongly disfavor the existence of Cooper pair co-tunneling.Our results shed valuable insight into the twisted cuprate systems and help settle down the debate on their underlying pairing symmetry.
technique with Ga ions.For S1 and S2, HAADF-STEM experiments were carried out in the FEI-Titan Cubed Themis 60-300 operating at 300 keV with double spherical aberration (Cs) correctors and a spatial resolution of 0.6 Å.For S3-S5, TEM experiments were performed in the JEOL-ARM200CF operating at 200 keV with double spherical aberration (Cs) correctors and a spatial resolution of 0.8 Å.To determine the twist angles accurately, we focused the electron beam on either top or bottom Bi-2201 to observe their respective Kikuchi patterns.We then calculated the angle on the separation between the two patterns.