## 1. Introduction

Superconductivity is a macroscopic quantum phenomenon [1]. Therefore it shows quite interesting properties because of its quantum nature. Such properties are described by a macroscopic complex wave function of the superconductivity. Especially, a phase of the macroscopic wave function play an important role in these properties. For example, superconducting devices, such as, superconducting charge, flux and phase qubits, superconducting single flux quantum device, and intrinsic Josephson junction Terahertz emitter of high Tccuprate superconductors, use such quantum nature of superconductivity. They have attracted much attention recently.

In conventional superconductors, there is only single phase

The high-Tccuprate superconductors are typical example of d-wave superconductors. Since the discovery of the high-Tc superconductors, the pairing symmetry, as well as its origin, has been controversial, but phase-sensitive experiment is crucial for determining the symmetry of the Cooper pairs. One such phase-sensitive experiment is a corner junction experiment [4,5]. In this experiment, a square-shaped high-Tc superconductor is connected with a conventional s-wave superconductor with two Josephson junctions A and B, around a corner and each junction is perpendicular to either of the x or the y axes of the high-Tc superconductor, as shown in Fig. 2.

In such geometry if the Cooper pair tunnels through junction A has phase 0, then the Cooper pair tunnels through junction B has phase π. The critical current between the s- and the d-wave superconductors through these two junctions is zero under zero external magnetic field because each supercurrent cancels with each other. This is in apparent contrast with a corner junction between both s-wave superconductors for which the critical current is largest under zero external field. Therefore, from this experiment, the symmetry of the Cooper pairs in high-Tc superconductors was determined as d-wave, especially

In addition, this experiment showed that the critical current becomesmaxwhen the total magnetic flux through the corner junction was

where C is a contour around the corner junction and *S* is the area surrounded by the C and

This result shows, the stable superconducting state under zero external current and zero external field becomes nontrivial. For such stable state, the free energy of whole system, especially the superconducting condensation energies of both d- and s-wave superconductors should be low. Therefore the order parameter should become continuous across both of junctions. And then the phase of order parameter cannot be uniform because of the phase differenceπbetween two junctions. This phase difference π should be compensated by changing the phase of the order parameter spatially in both s- and d-wave superconducting regions because of singlevaluedness of the order parameter. This spatial variation of the order parameter causes the supercurrent around the corner junction and then spontaneous magnetic flux is created at the corner without an external field. Because the associated phase change of this supercurrent is not 2π but only π, the spontaneous magnetic flux is not singly quantized magnetic flux

An experiment showing this property was done by Higenkamp et al., whomade a zigzag junction between conventional s-wave superconductor Nb and high-Tc d-wave superconductor YBCO, which consists of successive corner junctions [6]. And they observed spontaneous magnetic fluxes at every corner under a zero external field, using scanning SQUID microscope. The spontaneous magnetic fluxes aligned antiferromagnetically, because of the attractive vortex-anti-vortex interaction. They also made small ring with two junctions between Nb and YBCO, which is called π-rings, and controlled the spontaneous half-quantum magnetic fluxes [7].

When spontaneous magnetic flux appears in the zero external magnetic field, this state does not have the time reversal symmetry, which the original system has. Therefore there are always two degenerate stable states. In these states, supercurrent directions are opposite and henceforth directions of spontaneous magnetic fluxes are opposite. This spontaneously appeared magnetic flux is useful and can be used as a spin or a bit by it self. Ioffe et al. [8] proposed a quantum bit using this half-quantum flux. Using the spontaneous magnetic flux as an Ising spin system, Kirtley et al. made a frustrated triangular lattice of π-rings [7]. In their systems, the π-rings were isolated and interacted with each other purely by the electromagnetic force.

In contrast to the previous approach for using the spontaneous half-flux quantum, we considered nano-sized d-wave superconductor embedded in an s-wave matrix, as shown in Fig. 3 [9-11]. We want to consider the whole d-dot system as a single element, not the individual half-quantum fluxes. As in the single half-quantum system, our d-dot has two degenerate states if the spontaneous magnetic fluxes appear, because the state with spontaneous magnetic fluxes under zero external field also breaks time-reversal symmetry. This property is independent from the shape, and the d-dot in any shape always has two degenerate stable states. Therefore, the d-dot as a whole can be considered as a single element with two level states and it might be used as a spin or a bit also as a qubit. It has better properties than those of single flux quantum element, which will be shown infollowing sections.

In the following, we first show a phenomenological superconducting theory, which describes the spontaneous magnetic fluxes in these composite structures, especially in d-dot systems and then we discuss the basic properties of this d-dot, based on this phenomenological theory. Also, we discuss the difference between a single half-quantum flux system and our d-dot system.

## 2. Model: two components Ginzburg-Landau free energy

In order to discuss the basic properties of the d-dot, especially to describe the appearance of spontaneous magnetic fluxes, we use the phenomenological Ginzburg-Landau (GL) theory. However, for anisotropic superconductors, such as the

Multi-components GL equations were used for exotic superconductors, e.g. heavy fermion superconductors[16-18] and Sr_{2}RuO_{4}[19,20]. Also, recently, two components GL equationsare studied for two-band or two-gap superconductors, such as MgB_{2} [21,22].

The model byRen.et al. especially emphasize the anisotropy of d-wave superconductivity. Therefore, we use the following two-component Ginzburg-Landau (GL) free energy for d-wave superconductors;

where

Also, for s-wave superconductors, we use the following two-component GL equationwith attractive and repulsive interactions for the s-wave and the d-wave channels, respectively:

Here

In these free energies, the anisotropy of the d-wave superconductivity appears in the coupling terms of the gradient of both components of order parameters, *x* and the *y* directions have different signs.

By minimizingthe sum of these free energies, we obtain the Ginzburg-Landau (GL) equations. For this purpose, we use the finite-element method [23-25], because we want to investigate variously shaped d-dots. Hereafter we consider two-dimensional system and ignore the variation of physical quantities along the direction perpendicular to the two dimensional system. In the finite-element method for two-dimensional system, we divide the superconductor region into small triangular elements (see Fig. 4 (a)) and then expand the order parameters and the vector potential in each element by using the area coordinate, which is defined as,

where_{j} and P_{k} in Fig. 5 (b), where

Using these properties, the order parameters are expanded as

and also the magnetic vector potentialis expanded as,

where

And for the s-wave superconducting order parameter,

Coefficients for d-wave region and s-wave region are different and are given in Appendix.

Also for the vector potential following equations are obtained as follows,

Here coefficients are also given in Appendix. At the boundary of d- and s-wave superconductors, because the wave function is continuous, following boundary conditions are applied:

In this model, d- and s-wavecomponents of the order parameter interfere anisotropicallywith each other in the both of d-wave and s-wave superconductors. And the boundaries between both superconductors are assumed clean. But we cannot take into account the roughness of the boundary.

For treating roughness of the junctions, we consider second model. In contrast to the first model, in the second model, the s-wave superconductor is connected to the d-wave superconductor through a thin metal or an insulator layers. In the second model, we consider only the d-wave (s-wave) component of the order parameter in the d-wave (s-wave) superconducting region, but in the thin metal layer, we take both components. So, free energies are given as,

## 3. Appearance of the half-quantum magnetic flux in d-dot

Using the two models in Section 2, the spontaneous magnetic flux can be described.In Fig.7, the s-wave and d-wave order parameters and magnetic field distribution are shown for a corner junction between s-wave and d-wave superconductors under zero external magnetic field.

In this calculation, the physical parameters, such as the GL parameters and the transition temperatures, for Nb and YBCO in the s- and the d-wave superconducting regions are used, and system size is set as

How such spontaneous magnetic flux appears, when a square d-wave superconductor is embedded in an s-wave superconductor? In Fig.8, the typical magnetic flux structure of such square d-dot in which edge of square is parallel to the x- or y-directions, under zero magnetic field is shown. In this figure, red color means *H*_{ z}< 0. Magnetic fluxes appear at four corners of the square d-dot, and they order antiferromagnetically. Such antiferromagnetic order is already obtained by the experiment for zigzag junctions by Hilgencamp et al. [6]. The reason of this antiferromagnetic order can be understood by considering the interaction between vortices. Usual theory of vortex interaction tells us that parallel vortices repel each other and antiparallel vortices attract each other because of the current around the vortices and Lorentz force by this current. For the d-dot in Fig.8, the spontaneous current flows form top-center of the d-wave superconductors, turns right and left and flows into left and right sides of d-wave superconductor. So this current naturally creates upward flux at the left-top corner and downward flux at the right-top corner. Fluxes at the lower two corners can be explained similarly.

The second model of d-dots (Eqs.21-23) also describes these spontaneous half-quantum magnetic fluxes. Using this model, we argue the size-dependence of these magnetic fluxes. Question is: The d-dot in Fig.8 have four half-quantum magnetic fluxes, but if the size of d-dot becomes small, then what happens? In Fig.9, the typical size dependences of the order parameter structure and the magnetic fluxes are shown. In these figures, for smaller d-dots the amplitudes of d-wave order parameter and magnetic field distributions are shown. When the size of the d-dot becomes small, then the amplitude of d-wave component of the order parameter becomes small, and the spontaneous magnetic field becomes small. This is because the spontaneous current, which flows around each corner, does not decrease much at the center of the edge of the square when the size of the d-dot is comparable to the coherence length. Then, the total magnetic flux becomes less than the half-flux quantum, although the fluxoid is still a half-flux quantum.

The temperature dependences of amplitude of order parameters and the spontaneous magnetic fluxes also show similar tendency, because in this case, the coherence length increases with increasing temperature.

These spontaneous magnetic fluxes do not appears for the geometry for a square d-wave superconductor embedded in an s-wave superconducting matrix, but they appear for a square s-wave superconductor embedded in d-wave superconducting matrix and various shaped d-wave superconductors embedded in the s-wave superconducting matrix. The s-wave superconducting dot case is called “s-dot”. For an s-dot, order parameter structures and field distribution are shown in Fig. 10. There appear spontaneous magnetic fluxes around the corners antiferromagnetically. These magnetic fluxes also are explained similarly by the phase anisotropy of d-wave superconductivity. But the distribution of magnetic field is opposite to the d-dot case, where magnetic field appears mainly outside of inner superconducting region, but in this case, magnetic field appears mainly inside of inner superconducting region. The s-dot is also useful but in the following we focus on the d-dot case, which is more easily fabricated, we think.

Next we discuss the shape dependence of the d-dot. Even if the shape of the d-wave superconducting region is different from the square that is parallel to the crystal axis or x- and y-axis, the spontaneous magnetic field is also expected. In Fig. 11, distributions of the order parameters and magnetic field for an equilateral triangle plate are shown. For this case the spontaneous magnetic fluxes appear along the upper edges connected to the top corner. The spontaneous current flows the top corner and return to the intermediate points of upper edges. Also spontaneous magnetic fluxes appear around the lower corners. These spontaneous magnetic field along the edge also appear for rotated or diagonal squares, as shown in Fig. 12. In these figures, the spontaneous currents across the junction mainly flow along x- or y- directions. And direction of junction between d- and s-wave superconductors is important for the appearance of magnetic flux.

Also we note that not only square d-dots that is parallel to the x- and y-axis but also arbitrary shaped d-dots that show spontaneous magnetic field, such as in Figs. 11 and 12, have doubly degenerate stable states. And the shape of d-dot controls 15 the magnetic field distribution.

These doubly degenerate states have good properties for applications. First they remain under weak external magnetic field. In Fig. 13, external field dependence of spontaneous magnetic field distributions for a square d-dot is shown. The spontaneous magnetic flux parallel (anti-parallel) to the external field becomes large (small), respectively. Although the degeneracy from broken time reversal symmetry is lifted under the external magnetic field, the state with π/2 rotated magnetic field distributions are equally stable. These doubly degenerate states come from the broken four-fold symmetry of the square shape. This property depends on the shape of d-dots. For asymmetric shaped d-dots, one of the doubly degenerate states becomes more stable than another state. This means that we can control these degenerate states using the magnetic field for asymmetric d-dots.

Second property is the stability of the degenerate stable states. For the square d-dots, most stable states show antiferromagnetic order of spontaneous magnetic fluxes, and we call these state udud (up down up down) (Fig. 4 (a)). There are other states, which have higher free energy. In Fig. 14 (b), one of such states is shown. In this state, spontaneous magnetic fluxes do not show antiferromagnetic order, but parallel magnetic fluxes align at the upper or lower edges. We call this state uudd (up up down down). The free energies of the udud and uudd states are shown in Fig. 15. Well below the critical temperature of the s-wave superconductor*T*

_{cs}, free energy difference between udud (a) and uudd (b) states becomes comparable to the condensation energy of the superconductor. Therefore we can treat them as two-level systems.

## 4. Interaction between d-dots

As shown in previous section, the d-dots have double degenerate stable states. So we can use them as bits or 1/2 spins. In order to use them as artificial spins, the d-dots will be placed periodically or randomly. Then the interaction between them is important for these spin systems. For using the d-dots as computational bits, they are also placed to transform the information. Therefore interaction between d-dots also important for these applications.

How the d-dots interact with each other? Interaction between d-dots basically comes from the interaction between spontaneous magnetic fluxes or vortices. If the spontaneous vortices are independent from each other, that is, if there is no current flow between the vortices, then via a purely electromagnetic interaction, they interact. This is the case of the π-ring system of Kirtley et al. [4]. If the vortices are interacting in the same superconductors, there is a supercurrent flow around the vortices and ordinary vortices interact with each other through this current. The current distribution around a singly quantized vortex is given by the first order Bessel function. And therefore the interaction force to vortex 1 from vortex 2 is given as,

Where

For d-dots, there is an s-wave region between d-wave islands, and the spontaneous currents around the corners affect each other as usual supercurrent around singly quantized vortices, mentioned above. Then, an interaction between d-dots arises through the spontaneous currents in the s-wave region. If two square d-dots are in a line (Fig. 16 (a)), nearest vortex pair should be antiparallel (e.g. right-upper flux in the left d-dot and left upper vortex in the right d-dot should be antiparallel). And then we can expect that two d-dots will have a same spontaneous magnetic flux distribution, as shown in Fig. 16 (a). Therefore when we regard a d-dot as a spin, they interact ferromagnetically. In contrast to this configuration, if two square d-dots are placed diagonally (Fig. 16 (b)), then we expect an antiferromagnetic interaction by the same argument. In Fig. 17, stable states for parallel two d-dots are shown. For short distance ((e)), the nearest magnetic fluxes disappear and the flux distributions of two d-dots are same. Increasing the distance, the nearest magnetic fluxes appear and becomes gradually large ((f)-(h)) and the states of the d-dots are still same.Therefore ferromagnetic states are always stable, as we expected, and this is independent from the distance.

However, for diagonally placed two d-dots, the magnetic field distribution is not so simple. In Fig. 18, the stable state is shown for short distance between two diagonally placed d-dots. Unlike to the expectation from the usual vortex-vortex interaction, ferromagnetic state becomes stable. This is because two adjacent half-quantum fluxes are connected and form broad single quantum flux and this can occur when two fluxes are parallel. Therefore the ferromagnetic state becomes stable. It seems that the energy of this configuration becomes lower when two half-quantum fluxes are closer to each other. This property does not appear for ordinary vortices and we call this a fusion of half-quantum vortices. In Fig. 19, distance dependence of the free energies for the ferromagnetic and antiferromagnetic states are plotted. When distance between two d-dots is short, the ferromagnetic state has much lower free energy than that of the antiferromagnetic state. For longer distance, free energies of these states are almost same, because the spontaneous magnetic fluxes become almost independent as shown in the insets of Fig. 19.

These interactions between d-dots can be used for changing the state of the d-dots [26]. Using these interactions, the d-dots have a potential applicability to the superconducting devices[27].

## 5. Conclusions

We showed that the anisotropic pairing of superconductivity causesinteresting phenomena when two different superconductors are combined. Spontaneous half-quantized magnetic flux around the junction between d-wave and s-wave superconductors is one of such phenomena. It is simulated by the two-componentsGinzburg-Landau equations, which can treat the anisotropy of d-wave superconductivity. This half-quantum magnetic state for the d-wave superconductors embedded in the s-wave superconductor has good properties for applications to superconducting devices, especially for classical bits or qubits.

AppendixIn this appendix, the coefficients in Eqs.13-18 aregiven. For d-wave superconducting regions, they are defined as,

For s-wave superconducting regions, they are defined similarly.

The coefficients in Eqs. 17 and 18 in the d-wave superconducting region are given as,

One of the authors (M. K.) thanks his former and present graduate students, M. Ako, M. Hirayama, S. Nakajima, H. Suematsu, T.Minamino,S. Tomita, Y. Niwa, and D. Fujibayashi for useful discussions. This work is partly supported by “The Faculty Innovation Research Project” of Osaka Prefecture University and it was supported by the CREST-JST project.