# Laser-induced topological superconductivity in cuprate thin films

###### Abstract

We propose a possible way to realize topological superconductivity with application of laser light to superconducting cuprate thin films. Applying Floquet theory to a model of -wave superconductors with Rashba spin-orbit coupling, we derive an effective model and discuss its topological nature. Interplay of the Rashba spin-orbit coupling and the laser light effect induces the synthetic magnetic fields, thus making the system gapped. Then the system acquires the topologically non-trivial nature which is characterized by Chern numbers. The effective magnetic fields do not create the vortices in superconductors, and thus the proposed scheme provides a promising way to dynamically realize a topological superconductor in cuprates. We also discuss an experimental way to detect the signature.

^{†}

^{†}preprint: APS/123-QED

## I Introduction

Topological superconductors (TSCs) have attracted a great deal of interest from the viewpoint of realization of Majorana fermions in solid states and a possible application to quantum computation Qi and Zhang (2011); Sato and Fujimoto (2016). However, it is the current situation that experimental realization of TSCs is still limited. There are two main ways to search for topological superconductors. First one is to engineer a TSC in artificial systems by proximity effect (artificial TSCs). In recent studies, there are substantial developments in some artificial systems, such as a ferromagnetic atomic chain on a superconductor Nadj-Perge et al. (2014) or a nanowire on a superconductor Mourik et al. (2012). In these artificial systems, the realization of TSCs is confirmed to a certain extent. Second way is to find a material which is intrinsically a TSC (intrinsic TSCs). Some candidate materials of intrinsic TSCs have been proposed, for instance, spin-triplet superconductor Kashiwaya et al. (2011) and doped topological insulators Sasaki et al. (2011); Matano et al. (2016); Yonezawa et al. (2016). However, nodal excitation in nearly gapless NishiZaki et al. (2000); Ishida et al. (2000) is harmful for experimental detection of topological response, and topological nontriviality in is still under debate Levy et al. (2013). Therefore, further search of intrinsic TSCs is one of the important issues in the research field.

On the other hand, in recent years, tremendous developments have been achieved in the realization of topological phases by laser light applicationsOka and Aoki (2009); Lindner et al. (2011); Grushin et al. (2014); Titum et al. (2016). A typical example of the laser-induced topological state is a quantum Hall state in graphene Oka and Aoki (2009). This state is induced not by static magnetic field but by circularly polarized dynamical laser light. In this case, laser light effectively induces the next-nearest hopping with complex phase, which makes the system gapped, and thus the topologically non-trivial states similar to Haldane model Haldane (1988) are realized. Though this phenomenon has not been confirmed experimentally in graphene, a similar phenomenon is observed on the surface of the laser-irradiated topological insulators by time-resolved ARPES (Tr-ARPES) experiments Wang et al. (2013). The ARPES image obtained in the experiments shows that the surface Dirac cone becomes gapped when the laser light is applied to the system. Thus, we can regard the laser light as a new tool for controlling the topology of the states of matter.

Motivated by these situations, in this paper, we propose a possible way to realize TSCs with application of the laser light to well-known materials. We discover that the TSCs can be realized in -wave superconductors, such as cuprate, fabricated on a substrate irradiated by circularly polarized laser light. We apply Floquet theory to a model of -wave superconductor and derive an effective model under the irradiation of the laser light. Based on this effective model, we reveal that the system acquires the topologically nontrivial nature, which is characterized by Chern numbers, and show that the laser-induced magnetic field in the effective model play a crucial role in realizing TSCs.

This paper is organized as follows. In Sec. II, we introduce our model and methods. Next we show the derivation of an effective model which describes laser-irradiated cuprate thin films in Sec. III. In Sec. IV, we discuss topological properties of the effective model. We show the topological phase diagram and clarify the nature of each phase. In Sec.V, we discuss the experimental conditions to realize TSCs. Finally, summary and outlook are presented in Sec. VI.

## Ii Model and Methods

A setup of the system is schematically shown in Fig. 1. We consider a cuprate thin film fabricated on a substrate. Because of the asymmetric potential due to the substrate, Rashba spin-orbit coupling appears. In order to describe such a situation, we introduce a Rashba-Hubbard model as

(1) |

where

(2) | ||||

(3) |

with being the annihilation operator of electrons with momentum and spin . We choose the form of , which includes the next-nearest neighbor hopping , in order to reproduce the Fermi surface of typical high- cuprates.

Next we consider the effect of laser light. We treat the laser light as time-dependent classical electromagnetic fields and introduce them as Peierls phases. This treatment is equivalent to substituting with . With this substitution, we obtain the time-dependent model, which describes laser-illuminated cuprate thin films, as

(4) |

with , which corresponds to circularly () or elliptically () polarized laser light. Generally speaking, it is difficult to solve time-dependent quantum many-body problems. However, as for time-periodic problems, we can solve them by Floquet theory, which is known as a useful tool for analyzing time-periodic systems. The model Hamiltonian (4) is periodic in time and therefore we can apply Floquet theory to it.

Floquet theory is based on Floquet theorem, which is, so to speak, “Bloch’s theorem for time direction”: If the Hamiltonian is time-periodic , the eigenfunction can be written by a product of an exponential function and a time periodic function . is called “pesudo energy” that is defined in the range of . Then, we can define the effective Hamiltonian, of which eigenvalues are pesudo energy, as

(5) |

where the time-evolution operator . By definition, the effective Hamiltonian has only the information of and gives “stroboscopic description” of this system. In other words, the information of the time range of is neglected. Intuitively, if the time period is short enough, this Hamiltonian tells us the asymptotic behavior of periodically-driven systems. Although it is not obvious whether the effective Hamiltonian fully describes the quantum states while the system is driven D’Alessio and Rigol (2014), it is recently shown for many-body quantum systems that nonequilibrium steady states can appear with finite life time when the frequency is sufficiently highKuwahara et al. (2016). It is also shown that the nonequilibrium steady states are described by an effective Hamiltonian in Floquet theoryKuwahara et al. (2016). Thus we can use the effective Hamiltonian for understanding the nature of the non-equilibrium steady states of laser-irradiated materials including strongly correlated electron systems.

Though it is difficult to directly calculate the effective Hamiltonian from the definition (5), there are some useful methods to derive the effective Hamiltonian. We adopt a perturbative expansion in Mikami et al. (2016). Using this approach, we can write down the effective Hamiltonian as

(6) |

where . The second term in the case of the laser-irradiated systems represents the second order perturbation process of the -photon absorption and -photon emission in off-resonant light. If the laser intensity is sufficiently small, it is reduced to virtual processes related to one-photon absorption and emission .

## Iii Derivation of the effective model

In this section, we calculate the effective Hamiltonian (6) up to the first order in . The effective model describes the nonequilibrium steady states of the irradiated cuprate thin film. The model is obtained as

(7) |

where

(8) | ||||

(9) | ||||

(10) | ||||

(11) |

and represents the -th Bessel function. We see two effects induced by the laser light. The first one is so-called “dynamical localization”Dunlap and Kenkre (1986). With this effect, the hopping amplitude and the coupling constant are renormalized by the 0-th Bessel function. This effect induces the deformation of Fermi surface, resulting in topological phase transitions as we mention below. The second one is “laser-induced magnetic field”. It causes the Zeeman splitting of which splitting-width varies in momentum space. It plays a crucial role in realizing TSC in this system.

Next we consider the Bogoliubov-de Gennes (BdG) Hamiltonian which describes superconducting states. The order parameter of cuprate superconductors is known to be -wave Yanase et al. (2003). However, in our system, the Rahba spin-orbit coupling breaks the inversion symmetry and thus -wave pairing may be admixed with -wave pairing Tada et al. (2009). Therefore we investigate the +-wave superconducting state by adopting a simple form with and . We assume that is much larger than . With the + wave superconducting order parameter, we write down the BdG Hamiltonian as , where

(12) | ||||

(13) |

and . Without laser light (), this model represents the original wave superconductor and thus it has point nodes (shown in Fig. 2 (a)). With finite intensity of laser light (), the point nodes are gapped out (shown in Fig. 2 (b)). Later we show that the TSC is realized and the chiral Majorana edge modes appear in the laser-induced gap (shown in Fig. 3). This gap opening is caused by the laser-induced magnetic field (10). This term breaks the time-reversal symmetry and changes the symmetry class of the BdG Hamiltonian to class D.

Before closing this section, we remark on differences from the case of usual magnetic field applied to cuprate superconductors. When the orbital depairing effect is neglected, this case is described by the model similar to ours, which has already been studied by Yoshida and Yanase Yoshida and Yanase (2016); Daido and Yanase (2016). However, there are two important differences from their studies. First, the laser-induced magnetic fields do not induce vortices in superconductors. Usual magnetic fields induce vortices and easily suppress the superconducting states. On the other hand, the synthetic magnetic field induced by the laser light leads to only the Zeeman-type energy splitting (shown in Eq. (7) ) and oscillating gauge fields do not induce the vortices in superconductors. This is an advantage in realizing the TSC. Second one is that the laser light induces hopping renormalization (dynamical localization). As we mention below, the structure of superconducting gap is modified by the dynamical localization, and thus the topological states, which cannot be stabilized by usual magnetic fields, are realized.

## Iv Topological properties

### iv.1 Chern number and Phase diagram

As we mentioned above, our BdG Hamiltonian belongs to the symmetry class D in two dimensions. Therefore, the gapped state of this model is specified by the Chern number Thouless et al. (1982) , which is defined by

(14) |

We calculate the Chern number by two methods. The first one is an analytical method, which is proposed by Daido and Yanase Daido and Yanase (2016). Their derivation can be straightforwardly applied to our model. The analytic form is obtained as

(15) |

where is a unit vector in the -direction, and . The summation is taken over all the gapped nodes on Fermi surfaces of the bands. The energy spectrum is written as Daido and Yanase (2016)

(16) | ||||

(17) |

Therefore, the gap nodes appear in the absence of the layer light at satisfying

(18) |

The gapped nodes at are intersections of a Fermi surface and zeros of order parameter .

Based on the analytic formula (15), we can evaluate the Chern number with counting the contribution from gapped nodes. Each gapped node gives a contribution or . Sign of each contribution can be estimated as follows. First we define the direction parallel to the Fermi surface of bands as . Next we check the change of the sign of , which is in the argument of the function of Eq. (15). When it changes from negative to positive along direction at gapped nodes, the contribution is , and vice versa. Summing up all the contributions, we obtain the Chern number of the total bands.

This analytic formula is very useful for understanding the origin of the Chern number. However, it is not convenient for systematic calculations in a broad range of parameters. Therefore, for systematic calculation in a broad range of parameters, we use another method to calculate the Chern number. This is called Fukui-Hatsugai-Suzuki method Fukui et al. (2005), which is an efficient numerical method to calculate the Chern number of the model defined on discretized momentum space. With this approach, we calculate the Chern number for each point and obtain the topological phase diagram which is shown in Fig. 4. In the phase diagram, the number of electrons is fixed by tuning the chemical potential . In some regions, the Chern number is finite, implying the TSCs. In the following subsections, Secs. IVB and IVC, we clarify the nature of the superconducting phases in the low intensity region () and in the high intensity region (), respectively.

### iv.2 Weak intensity region

In the weak intensity region, we find that a TSC specified by is realized in a broad range of parameters. Even with infinitesimally weak intensity of laser light, the TSC is realized, and thus it is possible to experimentally realize TSCs with relatively weak laser light. This is one of the main results of this study. The energy spectrum of the ribbon-shaped system is shown in Fig.(3). As expected from the bulk-edge correspondence, four chiral Majorana modes appear near the edge of the system. The number of chiral Majorana modes corresponds to the Chern number .

The reason why the phase is realized is understood with the analytic formula (15). As we mentioned above, the Chern number is determined by the gapped nodes, which are defined as intersections of a Fermi surface and zeros of order parameters (they are shown in Fig. 5). Dividing the contributions from two Fermi surfaces of the and bands, we write as .

We here evaluate . In the case of circularly polarized laser (), all of the four nodes on each Fermi surface are crystallographically equivalent since they are transformed by the four-fold rotation. Therefore, they give the same contributions Daido and Yanase (2016) and thus must be either 2 or -2. Moreover, the superconducting gaps of each Fermi surface can be adiabatically deformed to each other without closing the gap, and thus the contributions to the Chern number is equivalent. Therefore and we conclude that the Chern number is either 4 or -4. With the procedure mentioned in Sec. IVA, we find that the phase is realized. In the case of elliptically polarized laser(), the superconducting gap can be adiabatically deformed to that in the case of circularly polarized laser, and thus the Chern number is not changed. Indeed, the phase is realized in a broad parameter range of the laser light as shown in Fig. 4.

At the end of this subsection, we discuss the effect of laser light on the superconducting order parameter. Though the superconducting order parameter is assumed to be wave in this study, the nature of the superconducting order can be changed by laser light through two effects, the paramagnetic effect and the deformation of the Fermi surface. However, the laser-induced magnetic field is very small and its effect on the superconducting order parameter is negligible. In the weak intensity region, deformation of Fermi surface is also small. Therefore the assumption of the -wave superconducting order is valid in this region.

### iv.3 Strong intensity region

Even with strong intensity of laser light, the system shows the topologically non-trivial phase with . Furthermore, under elliptic light (), different topological phases, such as and are realized (The phases of and appear out of the region of Fig. 4). For example, the phase with appear in a finite region around the point , which is shown in Fig. 6(b). The appearance of and reflects the fact that the rotation symmetry of the superconducting gap is reduced from four-fold to two-fold under the elliptic light.

The reason why these non-trivial Chern numbers are realized can also be understood with the analytic formula (15). Among them, we explain the phase of and . In Figs. 6 (a) and 6 (c), we show the case of . We can see that the symmetry of superconducting gap is reduced due to the elliptic laser light, but the system still has two-fold rotational symmetry. By this symmetry, is restricted to or . As we mentioned above, the sign of each contribution can be estimated from the sign of , which is shown in Fig. 6 (a) by shading. The change of the sign is opposite between the band and band. Therefore and have opposite signs and thus . It is a topologically trivial state, which can be realized by strong laser irradiation.

Next we discuss the phase in Fig. 6 (b). The figure shows twelve nodes (four nodes in the left panel and eight nodes in the right panel). Due to two-fold rotational symmetry, is limited to or and must be or . Estimating the contribution from each node, we obtain and and thus the Chern number is equal to . Owing to numerical difficulties, the global phase diagram including the phases of and is not shown. However, we have confirmed that the phases of and appear in a finite region of the phase diagram as we mentioned above.

Finally, we remark on the effects of strong laser light on the superconducting order. As we mentioned at the end of Sec. IVB, there are two laser-induced effects, paramagnetic effect and deformation of Fermi surfaces. Even in the strong intensity region, laser-induced magnetic field is still small for cuprates. However, the Fermi surfaces are drastically deformed and thus it is possible that the superconducting order is modified. However, the present system is very likely to remain topologically-nontrivial as long as the nodal spin-singlet component is dominantDaido and Yanase (2016). Therefore, we expect TSCs even when the superconducting gap is more or less deformed from the original -wave one.

## V Experimental conditions

Our proposal has two advantages in experimentally realizing TSCs. First one is that the topologically-non trivial states can be realized by infinitesimal intensity of laser light as we mentioned above. Indeed, any fine tuning of parameters is not required for TSCs. Second one is that the laser-induced magnetic field gives rise to only the Zeeman effect (paramagnetic effect) and thus does not induce vortices. Because of these advanteges, a laser-irradiated cuprate thin film is a promising candidate for TSCs. In the following, we discuss the experimental conditions about materials, laser light and experimental methods.

Concerning candidate materials, a cuprate superconductor has advantages since its superconducting state is expected to be robust to the perturbations of laser light due to its high critical temperature. With a slight modification (e.g. the form of ), our calculation can be applied to any -wave superconductors with Rashba spin-orbit coupling. Thus some of the heavy fermion superconductors, such as Shimozawa et al. (2014), are also a candidate material for laser-induced topological superconductors.

In our model, the laser light is characterized by its frequency and intensity. Strictly speaking, the frequency must be sufficiently high and off-resonant since our calculation is based on the high-frequency expansion. Within the scope of our calculation, frequency must be higher than the band width 8. In solid state systems, there are so many unoccupied bands above Fermi energy, then we have to choose an appropriate frequency so as to make it off-resonant. On the other hand, even if the frequency is small, the signatures of the effective Hamiltonian may remain. For example, even by the low frequency laser light, a nodal point of Dirac cone is gapped out Wang et al. (2013), which is the behavior analogous to the case of off-resonant driving.

Regarding the intensity of the laser light, the topological superconducting states () can be realized by infinitesimal intensity, but weak intensity laser opens only a small gap, which is difficult to detect. The largest energy gap is realized at , where the 1-st Bessel function takes a maximum value(). The energy gap is estimated with eqs. (16) and (17) as,

(19) |

where the renormalized coupling constant is . The amplitude of admixed -wave component is estimated as Bauer and (eds.) (2012); Fujimoto (2007) and thus the superconducting gap is evaluated as when , and are adopted. In this case, the laser intensity corresponds to the electric field amplitude . In order to realize a large gap, we should prepare a system which has large and apply laser with intermediate frequency as long as the high-frequency expansion is appropriate and the off-resonant condition is satisfied.

To experimentally probe the phenomena proposed here, we need to use a method suitable for observing transient phenomena because we may use pulse laser in order to obtain strong intensity of laser light. Then, the TSC is realized only when laser light is applied. The most promising experimental tool is Tr-ARPES. The nodal structure of cuprate superconductors has been already observed in ARPES experimentsHashimoto et al. (2014). Therefore, we believe that the gap-opening at the nodes will be detected with Tr-ARPES.

## Vi Summary and Outlook

In this paper, we have suggested a possible way to realize TSCs with application of laser light to superconducting cuprate thin films Bollinger et al. (2011); Leng et al. (2011). Using Floquet theory, we have analyzed the model of -wave superconductors with Rashba spin-orbit coupling. We have derived an effective model and discussed its topological nature. The effective model includes the laser-induced magnetic fields, which make the system fully gapped and lead the system to the topologically non-trivial states characterized by Chern numbers. The laser-induced magnetic fields do not create vortices in the superconductors, and thus our proposal has an advantage in experimentally realizing TSCs in cuprates.

We have chosen the laser frequency as being off-resonant in this study, but using different frequencies enables us to realize a different effective model. For example, it is proposed that a driving resonant to the interaction induces the correlated hopping and changes the superconducting natureBukov et al. (2016). Engineering the quantum states beyond off-resonant regime is an interesting direction for future work.

###### Acknowledgements.

This work is supported by JSPS KAKENHI (KAKENHI Grants No. JP16J05078, No. JP15H05884, No. JP15K05164, No. JP15H05745, No. JP16H00991, No. JP15H05855, and No. JP16K05501 ). KT thanks JSPS for the support from a Research Fellowship for Young Scientists.## References

- Qi and Zhang (2011) X.-L. Qi and S.-C. Zhang, Rev. Mod. Phys. 83, 1057 (2011).
- Sato and Fujimoto (2016) M. Sato and S. Fujimoto, J. Phys. Soc. Jpn. 85, 072001 (2016).
- Nadj-Perge et al. (2014) S. Nadj-Perge, I. K. Drozdov, J. Li, H. Chen, S. Jeon, J. Seo, A. H. MacDonald, B. A. Bernevig, and A. Yazdani, Science 346, 602 (2014).
- Mourik et al. (2012) V. Mourik, K. Zuo, S. M. Frolov, S. R. Plissard, E. P. a. M. Bakkers, and L. P. Kouwenhoven, Science 336, 1003 (2012).
- Kashiwaya et al. (2011) S. Kashiwaya, H. Kashiwaya, H. Kambara, T. Furuta, H. Yaguchi, Y. Tanaka, and Y. Maeno, Phys. Rev. Lett. 107, 077003 (2011).
- Sasaki et al. (2011) S. Sasaki, M. Kriener, K. Segawa, K. Yada, Y. Tanaka, M. Sato, and Y. Ando, Phys. Rev. Lett. 107, 217001 (2011).
- Matano et al. (2016) K. Matano, M. Kriener, K. Segawa, Y. Ando, and G.-q. Zheng, Nat. Phys. 12, 852â (2016).
- Yonezawa et al. (2016) S. Yonezawa, K. Tajiri, S. Nakata, Y. Nagai, Z. Wang, K. Segawa, Y. Ando, and Y. Maeno, Nat. Phys. advance online publication (2016).
- NishiZaki et al. (2000) S. NishiZaki, Y. Maeno, and Z. Mao, J. Phys. Soc. Jpn. 69, 572 (2000).
- Ishida et al. (2000) K. Ishida, H. Mukuda, Y. Kitaoka, Z. Q. Mao, Y. Mori, and Y. Maeno, Phys. Rev. Lett. 84, 5387 (2000).
- Levy et al. (2013) N. Levy, T. Zhang, J. Ha, F. Sharifi, A. A. Talin, Y. Kuk, and J. A. Stroscio, Phys. Rev. Lett. 110, 117001 (2013).
- Oka and Aoki (2009) T. Oka and H. Aoki, Phys. Rev. B 79, 081406 (2009).
- Lindner et al. (2011) N. H. Lindner, G. Refael, and V. Galitski, Nat. Phys. 7, 490 (2011).
- Grushin et al. (2014) A. G. Grushin, A. Gómez-León, and T. Neupert, Phys. Rev. Lett. 112, 156801 (2014).
- Titum et al. (2016) P. Titum, E. Berg, M. S. Rudner, G. Refael, and N. H. Lindner, Phys. Rev. X 6, 021013 (2016).
- Haldane (1988) F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988).
- Wang et al. (2013) Y. H. Wang, H. Steinberg, P. Jarillo-Herrero, and N. Gedik, Science 342, 453 (2013).
- D’Alessio and Rigol (2014) L. D’Alessio and M. Rigol, Phys. Rev. X 4, 041048 (2014).
- Kuwahara et al. (2016) T. Kuwahara, T. Mori, and K. Saito, Ann. Phys. 367, 96 (2016).
- Mikami et al. (2016) T. Mikami, S. Kitamura, K. Yasuda, N. Tsuji, T. Oka, and H. Aoki, Phys. Rev. B 93, 144307 (2016).
- Dunlap and Kenkre (1986) D. H. Dunlap and V. M. Kenkre, Phys. Rev. B 34, 3625 (1986).
- Yanase et al. (2003) Y. Yanase, T. Jujo, T. Nomura, H. Ikeda, T. Hotta, and K. Yamada, Phys. Rep. 387, 1 (2003).
- Tada et al. (2009) Y. Tada, N. Kawakami, and S. Fujimoto, New J. Phys. 11 (2009).
- Yoshida and Yanase (2016) T. Yoshida and Y. Yanase, Phys. Rev. B 93, 054504 (2016).
- Daido and Yanase (2016) A. Daido and Y. Yanase, Phys. Rev. B 94, 054519 (2016).
- Thouless et al. (1982) D. J. Thouless, M. Kohmoto, M. P. Nightingale, and M. Den Nijs, Phys. Rev. Lett. 49, 405 (1982).
- Fukui et al. (2005) T. Fukui, Y. Hatsugai, and H. Suzuki, J. Phys. Soc. Jpn. 74, 1674 (2005).
- Shimozawa et al. (2014) M. Shimozawa, S. K. Goh, R. Endo, R. Kobayashi, T. Watashige, Y. Mizukami, H. Ikeda, H. Shishido, Y. Yanase, T. Terashima, T. Shibauchi, and Y. Matsuda, Phys. Rev. Lett. 112, 156404 (2014).
- Bauer and (eds.) (2012) E. Bauer and M. S. (eds.), Non-Centrosymmetric Superconductors: Introduction and Overview, Vol. 847 (Springer, Berlin, 2012).
- Fujimoto (2007) S. Fujimoto, J. Phys. Soc. Jpn. 76, 051008 (2007).
- Hashimoto et al. (2014) M. Hashimoto, I. M. Vishik, R.-H. He, T. P. Devereaux, and Z.-X. Shen, Nat. Phys. 10, 483 (2014).
- Bollinger et al. (2011) a. T. Bollinger, G. Dubuis, J. Yoon, D. Pavuna, J. Misewich, and I. Božović, Nature 472, 458 (2011).
- Leng et al. (2011) X. Leng, J. Garcia-Barriocanal, S. Bose, Y. Lee, and A. M. Goldman, Phys. Rev. Lett. 107, 027001 (2011).
- Bukov et al. (2016) M. Bukov, M. Kolodrubetz, and A. Polkovnikov, Phys. Rev. Lett. 116, 125301 (2016).