Radiation induced oscillating gap states of nonequilibrium twodimensional superconductors
Abstract
We study effects of infrared radiations on a twodimensional BCS superconductor coupled with a normal metal substrate through a tunneling barrier. The phase transition conditions are analyzed by inspecting stability of the system against perturbations of pairing potentials. We find an oscillating gap phase with a frequency not directly related to the radiation frequency but resulting from the asymmetry of electron density of states of the system as well as the tunneling amplitude. When such a superconductor is in contact with another superconductor, it will give rise to an unusual alternating Josephson current.
 PACS numbers

74.40.Gh
,74.25.N
,74.78.w
,74.50.
r+
pacs:
Valid PACS appear hereThe rapid development of timeresolved spectroscopy technology have drawn growing interests in the study of nonequilibrium phenomena. The conductivity properties of solid materials can be greatly changed in nonequlibrium states induced by radiation. In twodimensional electron gas, radiation induces zero resistance states at high Landau filling factors zudov2003evidence; mani2002zero; shi2003radiation; durst2003radiation. Recent discoveries on high temperature cuprate superconductors reveal that infrared radiations transform nonsuperconducting compounds into transient superconductors or enhance coherent superconducting transport even at temperature above the superconducting transition temperature fausti2011light; hu2014optically. For irradiated BCS superconductors, the enhancement of superconducting gap wyatt1966mic; dayem1967behav; Kom1977meas; beck2013transi or the oscillating amplitude modes in nonadiabatic regime bara2004collec; bara2006synch; yuz2005noneq; yuz2006relax; sch2011reson; matsu2013higgs induced at near gap frequencies have been discussed in many previous works. Recently, experimental studies on twodimensional superconductors have shown superconducting order remains robust in ultra thin crystalline films which are a few atomic layers thick qin2009superconductivity; zhang2010superconductivity; guo2004superconductivity; ozer2006hard; ozer2007tuning. It brings new opportunities to the study of nonequilibrium superconducting orders in reduced dimensions.
In this Letter, we investigate effects of radiations on a twodimensional BCS superconductor coupled with a normal metal through a tunneling barrier and find a new oscillating state induced by radiations which has not been discussed before. Different from the usual (fractional) ac Josephson effect josephson1962possible; shapiro1963josephson, in which oscillation is induced by dc voltages, the frequency of this alternating phase is not directly affected by the radiation frequency , but the radiation intensity, the amplitude of the tunneling interaction and most importantly, the asymmetry of density of states around Fermilevel. Such a state can be a probe of the internal properties of 2D superconductors.
The physical system we concern is a superconducting film coupled to a normal metal substrate by a tunneling barrier. The amplitude of the tunneling matrix element is determined by the thickness of the insulating film. We excite the superconductor with infrared radiations, the wavelength of which is assumed to be resonant with a phonon mode in the superconductor (shown in Fig.1).
As a result, the electrons feel an oscillating crystal field resulting from the excited phonons. The Hamiltonian of this tunneling system can be written as the sum of the Hamiltonians of N, S and the tunneling barrier,
(1) 
where and are annihilation operators of the single particle state in the superconductor and the normal metal substrate respectively. The Hamiltonian of superconductor with a single frequency oscillating field is written in the Nambu matrix form nambu1961dynamical; *gor1959microscopic, where and is the basis of Nambu representation, and is the electronelectron interaction. is the tunneling amplitude between S and N. For simplicity, and are omitted in the discussion.
To determine superconducting phase transition of the irradiated system, we apply the theory of linear response and investigate the system’s stability against a weak external paring potential which perturbs the system with . The response function of paring amplitude to the perturbation is defined by
(2) 
In the presence of external paring potential, the disturbance will reach a steady state in the normal state, while in the superconducting state, however small perturbation can lead to a disturbance increasing with the time. Superconductivity is a result of the instability of the system to the perturbing pairing potential.
To determine the response function of the system, we employ the RPAlike approach. First, we define as the pairing response function of a noninteracting system, i.e. in Eq. (1). When the interaction is included, an electron in the system will feel not only the external pairing potential, but also an induced pairing potential exerted by all other electrons in the system. As a result:
(3) 
The BCS mean field theory bardeen1957theory corresponds to assuming,
(4) 
where is the effective electronelectron interaction potential. Combine Eq. (2  4), we can express of the whole system in the form of operators:
(5) 
where is the operator form of the response function and is the inverse of the operator . Thus we can determine the irradiated superconducting system’s response to the perturbation as long as the noninteracting system’s response function is obtained.
We derive the expression of with the Keldysh Green function technique keldysh1965diagram; haug2008quantum. is proportional to the offdiagonal element of the less Green function defined as , where and . With the perturbation Hamiltonian , the Dyson equation of can be obtained with the Langreth Theorem. For noninteracting system, can be expanded to the first order of as
(6)  
where
(7) 
and () is the retarded (lesser) Green function of the noninteracting irradiated system without the external perturbation. We obtain in the presence of radiation:
(8) 
where is the energy level broadening resulting from the tunneling into the metal substrate. We assume that the selfenergy resulting from the tunneling is nearly a constant around the Fermi level at energy scale of . While the real part of the selfenergy can be absorbed in the electron dispersion, the imaginary part of the selfenergy is approximated to be a constant . can be obtained similarly. By substituting and into Eq.(6), is given as
(9)  
with
(10) 
where is the Bessel function of the first kind and comes from our assumption that the conservation of momentum is kept during the tunneling process.
With the response function derived, we can determine the conditions of the superconducting phase transitions. The stability of a linear system requires all poles of the response function lie on the lower half plane in the frequency space. The transition condition is determined when the first pole comes across the real axis during the change of the system parameters. Driven by periodic radiation field with frequency , and can be expressed in matrices defined by and , where and are integers. Equation (5) in the frequency domain shows that the pole of is determined by the zero point of . We can determine eigenvalues and eigenvectors of . The stability condition proposed by Bode bode1945network requires that for on the real axis, each eigenvalue of satisfies when . Thus the transitional point for superconductivity is that for the largest which satisfies , it should be . By this approach we set up a correspondence between the superconducting transitional temperature and the strength of the effective attractive potential , as well as the phase frequency of the gap phase .
In the approach above, we investigate the most unstable mode of the linear system against perturbation near the phase transition point. To obtain the physical observables at the temperature region far below , such as the steadystate gap value, nonlinear gap equations should be applied. Here we make the assumption that the most unstable solution mode of the linearized equation at the transition point will correspond to the steady state solution of the nonlinear system. The validity of our assumption can be assured at the close region.
The above analysis raises an interesting possibility, i.e., one may find a solution with . There is no a priori reason to believe that must be zero or harmonics of the radiation frequency . Such a solution implies an oscillating state whose frequency is neither zero nor integer (or half integer) multiple of the radiation frequency . Thus it will be a new radiationinduced nonequilibrium effect, different from the usual (fractional) ac Josephson effect. Without radiation, it can be proved that must be zero and Eq. (9) is consistent with the BCS gap equation. However in the irradiated state, there is no such conclusion and the occurrence of nonzero is possible. In the following paragraphs, we will discuss the existing conditions of this oscillating gap state numerically.
First we investigate the effects of radiation on the superconducting system with constant density of state (DOS) around the Fermi energy. In our calculation, to involve the possible subbands motivated by radiations, the effective interaction potential is taken to be , where is the magnitude of the attractive potential and is the step function. We assume that the radiation frequency is the order of the Debye frequency , which is in terahertz regime in most BCS superconductors, and that the radiation energy and the imaginary part of selfenergy due to tunneling are at the order of magnitude of . The electronphonon coupling constant of the superconducting system assumed as , where is the density of state at the Fermi level.
We show the dependence of the transition temperature on the magnitude of radiation energy at different tunneling amplitudes in Fig. 2(a). Here we make the simplification that the radiation energy close to the Fermi energy. We can see the transition temperature is suppressed by the radiation. curves at different values of show that tunneling to the substrate will also suppress the superconductivity. In this case, the alternating phase frequency we obtain is always zero.
To search for cases of nonzero , we investigate the superconducting system with particlehole asymmetry near the Fermi surface. We assume the superconductor in our system has a DOS with a linear inclination at the Fermi surface, which is expressed as , where is the slope and the total range of is limited to to avoid unphysical results. In such a system, we observe the occurrence of the radiationinduced oscillating state. In Fig. 3(a), the dependence of the oscillating frequency on the radiation intensity at different values of slope is shown. We see that increase with the DOS slope as well as the radiation intensity. shows a nearly linear relationship with the asymmetric degree of the DOS [Fig. 3(a) inset]. we obtain is in the regime of GHz, three or more orders of magnitude smaller than the radiation frequency. The significant difference of magnitude between this alternating phase and the radiation frequency could make the experimental observation of the effect easier.
Besides the asymmetry of the DOS, the occurrence of nonzero also depends on the tunneling amplitudes between the superconductor and the substrate. Fig. 3(b) shows with respect to the radiation intensity at different values of . increase rapidly with the tunneling amplitudes. corresponds to the bulk superconducting system, where the interactions with the substrate are screened over the length of penetration. At this case, is always zero, which indicates this radiation induced effect is special to the twodimensional superconductor.
We calculate the relative modulus of the gap by solving the eigenvectors of the operator . The evolution of the gap modulus over one period of radiation at the Fermi surface for systems with constant DOS [Fig. 2(b)] and with nonzero motivated at and [Fig. 2(c)] are plotted respectively. We can see that the radiation induce oscillation of the gap amplitude with the same frequency as the radiation. For system with nonzero , large proportion of the secondary and higher harmonic wave is also excited, as shown in Fig. 2(c).
When the system in Fig. 1 with nonzero is connected to a bulk superconductor or a 2D superconductor with particlehole symmetry, an unusual alternating Josephson current can be observed experimentally. When irradiated with infrared radiation, the phase difference of the two superconductor is . The fast oscillating term due to the infrared radiation is ignored. Thus the tunneling current is , where is affected by the radiation intensity and frequency. Unconventionally, the frequency of this alternating current depends on the radiation intensity and the degree of the electronhole asymmetry of the 2D superconductor.
In conclusion, we find a radiationinduced oscillating state of the gap phase in a twodimensional BCS superconductor with particlehole asymmetry coupled to a normal metal substrate. Its oscillating frequency is determined by the radiation intensity, the tunneling amplitude with the substrates and the electronhole asymmetry of the superconducting system, which is different from the (fractional) ac Josephson effect. When this system is connected to another superconductor in a Josephson junction, alternating current corresponding to this phase occurs.