Conventional superconductors are strong diamagnets that through the Meissner effect expel magnetic fields. It would therefore be surprising if a superconducting ground state would support spontaneous magnetics fields. Such time-reversal symmetry broken states have been proposed for the high-temperature superconductors, but their identification remains experimentally controversial. Here we show a route to a low-temperature superconducting state with broken time-reversal symmetry that may accommodate currently conflicting experiments. This state is characterised by an unusual vortex pattern in the form of a necklace of fractional vortices around the perimeter of the material, where neighbouring vortices have opposite current circulation. This vortex pattern is a result of a spectral rearrangement of current carrying states near the surfaces.The phase-sensitive experiments [1,2] carried out in the early 1990's showed that the high-temperature superconductors have predominantly d-wave pairing symmetry. Ever since, there has been an ongoing debate whether there exists a low-temperature phase that breaks time-reversal (T ) symmetry, in addition to the reflection symmetry of the crystal broken by the d-wave order parameter itself. Several experiments on tunnelling and charge transport support a phase transition into a state with broken T -symmetry at low temperatures [3][4][5][6][7][8]. On the other hand efforts to detect the concomitant spontaneous magnetic field have failed, or at best have put severe restrictions on how strong the subdominant order can be [9][10][11][12][13]. Below we show how to reconcile these two sets of experiments and how the route to a low-temperature superconducting state with broken Tsymmetry may occur. We emphasize that the superconducting state with broken T -symmetry we describe here does not necessarily imply a multicomponent superconducting order-parameter [14,15].The superconducting state in unconventional superconductors, such as the cuprates, is fragile to scattering off impurities, defects and surfaces [16]. Scattering leads to pair-breaking and formation of so-called Andreev states [16][17][18]. These states are formed by constructive normal and Andreev reflection processes, have energies ε A within the superconducting energy gap ∆, and are spatially bound to the scattering centers. For superconducting grains, when the size of the superconducting material is comparable with the superconducting coherence length, the whole superconducting state of the grain is affected by boundary scattering and the formation of Andreev states leads to properties that are not fully understood yet. For the realization of real devices, a deeper understanding of the ground state of grains is therefore called for.We consider meso-scaled grains of a d-wave superconductor and relax the assumption of translational invariance along the surface. The typical grain sizes we consider correspond to side lengths D ∼ 10ξ 0 to 100ξ 0 , where ξ 0 = v F /(2πk B T c ) is the superconducting coherence length (v F is the...