Unusual magnetic field-dependence of a possible hidden order phase

URu2Si2 exhibits a second-order phase transition at 17.5 K. Initially it was thought that the transition was to a spin density wave phase, however, subsequent measurements do not support this assignment. Despite the unknown nature of the order parameter, many experimental results can be described in terms of the formation of a generic density wave. Here, we report calculations on an unusual phase of the underscreened Anderson lattice model, the so called spin-dependent inter-orbital density wave that has been proposed as describing the “hidden order” phase of URu2Si2. We determine the effects of an applied magnetic field. Since the order parameter describes an ordering in the x–y plane, the electronic properties of the system are anisotropic below the critical temperature THO. We show that the magnetic susceptibility becomes anisotropic below THO. Furthermore, for fields applied along a spontaneously chosen hard axis, THO decreases towards zero and that the HO transition changes from second order to first order at a large value of the magnetic field. Also, we find that the system undergoes a cascade of field-induced Lifshitz transitions and also show how these properties originate from the dependence of the quasi-particle bands on the orientation of the applied field. The good qualitative agreement with experimental findings provides strong support for the proposed description of the HO phase as a spin-dependent inter-orbital density wave phase. Support for an explanation of a quantum phase that has remained mysterious for decades is provided by researchers in Brazil and the USA. Peter Riseborough from Temple University and colleagues determine the magnetic properties of a model proposed to describe the so-called ‘hidden order’. Macroscopic signatures of quantum phases are identifiable in experiments but the underlying microscopic origin is not always clear. One prominent example is uranium ruthenium silicide, which exhibits ordering below a temperature of 17.5 Kelvin. It is known that an applied pressure induces a phase change and yet, contradictorily, experiments indicate that there is no structural change. The team show that one model proposed to solve this problem, the underscreened Anderson lattice model, gives rise to an anisotropic magnetic susceptibility. This agrees with experimental results and provides strong backing for the model.


INTRODUCTION
Despite the 30 years that have elapsed since its discovery, the nature of the phase found in URu 2 Si 2 below the critical temperature T HO = 17.5 K 1, 2 remains unknown. This enigmatic low-temperature phase has come to be known as the "hidden order" (HO) phase. The HO transition is marked by the formation of a gap, which extends over 60% of the Fermi-surface. 2,3 Initially, the transition was thought to lead to a conventional spin density wave phase. However, inelastic neutron scattering measurements ruled out this possibility. 4 Nevertheless, the application of pressure causes a large moment antiferromagnetic phase to appear. 5,6 To complete the puzzling scenario, the preponderance of evidence indicates that no structural transition occurs. [7][8][9] Over the years that have elapsed since its discovery, many theories have been advanced to describe the HO phase. The proposals include quadrupolar, 10, 11 octapolar 12, 13 and higherorder multipolar ordering. [14][15][16] The predictions based on quadrupolar orderings are inconsistent with resonant X-ray diffraction experiments. 17,18 Observation of high-order multipolar ordering via either neutron or X-ray diffraction experiments may prove to be very difficult since, presently, this requires high-precision measurement of the form factors at high momentum transfers. Other recent proposed description of the ordered phase include, unconventional or d spin-density waves, 19 Pomeranchuk instabilities, 20 modulated spin liquids, 21 hybridization waves, 22 spinnematic ordering 23 and hastatic ordering. 24 Since an l = 1 Pomeranchuk instability is subject to severe limitations 25 and since the Fermi-surface measured by ARPES and de Haas-van Alphen measurements show evidence of a doubling of the unit cell, [26][27][28] the proposal that the HO transition in URu 2 Si 2 is a Pomeranchuk instability 20 may require revision to include either Q ≠ 0 or l > 1. If the ordering involves a hybridization wave state, 22 this should be accompanied by a charge density wave, 29 which could have been expected to be detected in one of the many high-resolution scattering experiments that have been performed. [7][8][9] Hastatic ordering 24,30 is accompanied by the formation of a static, in-plane, magnetization of the conduction electrons. The hastatic picture of URu 2 Si 2 and the phenomenological picture of Rau and Kee, 31 which also predicts a static transverse moment, are not supported by the results of neutron scattering experiments 32-34 that have placed an upper limit on the value of a transverse moment which is an order of magnitude smaller than the predicted hastatic value of 0.015 μ B . Due to the absence of a definitive identification of the symmetries that are broken at the HO transition, the nature of the HO phase remains an outstanding issue in heavy-fermion physics. 35,36 Here, we report an unusual magnetic field dependence of the properties found for a novel phase of the underscreened Anderson lattice (UAL) model 37 that was recently proposed to describe the HO phase of URu 2 Si 2 . The choice of an itinerant model is motivated by the absence of evidence of well-defined crystal field levels in inelastic neutron scattering data 38 and highresolution inelastic X-ray scattering experiments, which indicate that the widths of any crystal field excitations are simply too large to be considered as excitations of a localized system. 39 The itinerant description is supported by the determination 40-42 that uranium 5f electrons are delocalized and have an occupation of about 2.7 electrons per uranium and the participation of heavy itinerant 5f electrons in the Fermi-surface. 28 Although there is experimental evidence of hybridization between the 5f electrons and the conduction electrons, 43,44 electronic structure calculations indicate that the bands near the Fermi-energy are primarily of 5f character. 45 Nevertheless, small admixtures of conduction band character are found. Since there appears to be a scarcity of conduction electron weight in the vicinity of the Fermi-energy which is required for the hybridization to either completely screen 46 or delocalize the 5f moments, 47 the UAL model seems apt. 48 However, the UAL model needs to be modified by the inclusion of hopping processes whereby electrons are directly transferred between the 5f orbitals on neighboring atoms. 37 It is the delocalization of all the 5f states caused by the inter-site hopping process that counteracts the tendency for residual 5f local moments to order.
The new feature of our approach is that the candidate order parameter (OP) for the HO phase 37 reflects a very peculiar type of broken symmetry in the UAL model, i.e., a broken spin-rotational symmetry without the formation of a magnetic moment. The broken symmetry is driven by the spin-flip part of the Hund's rule interaction of the UAL model. The candidate OP is given by the expectation value of the non-hermitian operatorẑ χ 0 χ Q;σ (Q is the fixed commensurate nesting vector) defined aŝ which leads to a low-temperature phase that can be described as a spin-dependent inter-orbital density wave. More specifically, the ordering produces a state characterized by an inhomogeneous, spin-dependent mixing of f states with different orbital characters given by χ. It should be remarked the non-trivial spin dependence of the OP, 37 The OP is the expectation value of a non-Hermitean operator and, therefore, does not correspond to a measurable quantity. The proposed HO state does globally break a gauge-invariance of the Hamiltonian. The transition spontaneously generates a phase difference between the different 5f components of the wave function and, as such, can be considered as breaking chiral (axial) gauge invariance. Due to the relation displayed in Eq. (2) and the quenching of the orbital angular momentum, the phase with a purely imaginary OP, z χ 0 ;χ Q;σ ¼ exp ± i π 2 Â Ã r results in the HO state being time-reversal invariant. This is easily seen since the Wigner time-reversal operator T defined by where κ denotes complex conjugation, transforms the up-spin states into the down-spin states. However, general complex values of z χ 0 ;χ Q;σ produce identical magnetic properties since the orbital angular momentum is quenched by the hopping. In this respect, the itinerant UAL model is different from the model considered by Haule and Kotliar. 14 The two time-reversal invariant chiral states have properties similar to those found in Raman scattering experiments 49 since the two degenerate states both preserve time-reversal invariance, have no net spin, and have the same charge character. However, unlike the chiral states of Kung et al., 49 the relation between UAL chiral states to reflection symmetry is, as yet, undetermined. The time-inversion invariance of the HO order phase is not only supported by the results of Kung et al. 49 but also by recent Polar Kerr effect (PKE) measurements. 50 The PKE experiments identified a large field-trainable effect in the superconducting state but only showed a small signal at the HO transition which was attributed to the existence of parasitic phases. [51][52][53] The present work is aimed at the investigation of the consequences that the application of a magnetic field have on the ordered phase. The results of this work support the suggestion that the OP z χ;χ 0 Q , which appears as a result of the spontaneous breaking of spin-rotational invariance, is at work in URu 2 Si 2 . All our results can be traced directly to the behaviour of the 5f quasiparticle dispersion relations. The unusual nature of these dispersion relations follows from the peculiar spin-rotational symmetry breaking that is driven by the x and y components of the Hund's rule interaction in the UAL model. Since the ordering involves the x and y components of the Hund's rule interaction, the ordering produces a hard-axis and an easy-plane but does not produce an ordered magnetic moment. As a consequence, the quasi-particle dispersion relations of the HO phase are found to be affected differently by magnetic fields applied in the easy magnetic direction (x-axis) as opposed to the hard magnetic direction (z-axis).
The temperature versus magnetic field phase diagram and the anisotropy found in magnetic torque measurements are believed to be crucial ingredients for the discover of the identity of the HO phase. Therefore, they are important constraints for any candidate theory for the elucidation of the HO phase in URu 2 Si 2 . In the first case, the experimental field dependence of T HO displays different behaviour for fields applied along the tetragonal c and a axes. For the c-axis, T HO decreases when the field is increased. However, for the b-axis T HO is only weakly affected. 54 Moreover, the application of large magnetic fields along the tetragonal c-axis of URu 2 Si 2 shows that the HO phase is destroyed at around 35 T. The secondorder boundary line between the non-magnetic HO and the normal-paramagnetic phase is replaced by a first-order one around 25 T. 55,56 For fields between 35 and 40 T, the T-H phase diagram contains a dome of an unknown phase. In the case of magnetic torque measurements, 57 it has been reported that, below the ordering temperature, the linear susceptibility becomes anisotropic in the tetragonal a-b plane suggesting that the fourfold rotational symmetry is broken in the low-temperature phase. Since the underlying crystal apparently remains tetragonal, 9 the nematic phase is expected to have a purely electronic origin. Moreover, the experimental findings 58 for the non-linear susceptibility χ 3 as function of temperature shows a sharp anomaly at T = T HO for the fields parallel to the tetragonal c-axis whereas it has a very weak temperature dependence for fields perpendicular to the c-axis.
It is to be noted that since this work uses cubic symmetry, the results only depend on the direction of the magnetic field relative to the hard magnetic direction and not to the crystallographic axes (denoted by a, b and c). This hinders direct comparision with experimental results on URu 2 Si 2 . Nevertheless, the particular form of rotational symmetry breaking in our theory produces a spontaneous selection of an easy-plane and hard axis which can lead to anisotropic responses for applied magnetic field resembling those found in the URu 2 Si 2 .
The paper is structured as follows: in section II we introduce the model and calculate the free energy, and present a detailed discussion of the results. The conclusions are presented in the last section.

RESULTS
We use the UAL model which describes two degenerate localized 5f bands (labeled by χ = α, β), which acquire itinerant character by direct hopping between neighboring 5f shells and by their mixing with one itinerant conduction band with dispersion relation ϵ d k ð Þ. The two degenerate narrow f bands exhibit intraband and where N is the number of sites. At the mean-field level, where the fluctuations inẑ χχ 0 q;σ andn χ f ;σ (the f-electron occupation number operator) are considered as small, the free energy is given as: and κ χ 0 χ Q;σ is a complex number defined by where the non-zero z χ 0 χ ÀQ;σ gives the novel spin-dependent interorbital density wave reflecting the breaking of the spin rotationally invariance within the Hund's rule part of the interaction.
The quasi-particle dispersion relations E ðjÞ k;σ for electrons of spin σ satisfy where I is the identity matrix and the matrix Λ σ ω; k ð Þ is given below in which s =f or d, and a is the lattice parameter. The gyromagnetic factors are assumed to be equal, i.e., g f = g d = 1. The OP z αβ Q;σ can be obtained from the free energy on minimizing F with respect to z βα ÀQ;Àσ . We define the susceptibility g σ for the OP via z αβ Q;σ ¼ Àg σ κ βα Q;σ . The tight-binding parameter values used in the calculations in the next section correspond to t d = −1/6, t f = 1/20 and V = 1/10 eV. These values were chosen as generic values but give rise to a width 2W f of the 5f band due to direct 5f hopping of 0.6 eV which is a little smaller than the 5f band widths of ≈1 and 1.5 eV determined, respectively, by DFT + DMFT 14 and DFT 45 calculations. Apart from raising the degeneracy of the α and β bands, the value of V has little effect on the results. The value of the effective f level energy E f was chosen to be 0.3 eV and the value of μ(T,H) was determined to be around 0.32, so that the 5f bands will be approximately half-filled. The values of U and J were chosen to be equal, in order to avoid an unnecessary increase the number of parameters of the theory. However, estimates of the ratio J/U found in the literature range from 0.49 45 to 0.14. 14 We find that smaller values of the ratio have the effect of increasing the stability of the HO phase wrt a neighbouring AF phase that is stable at smaller W f values. The value of U = 0.074 eV was chosen in order to give T HO = 17.5 K and corresponds to a ratio U/W f ≈ 0.25 which is close to the range 0.26-0.4 used in (J = 0) DFT + DMFT calculations. 45 Throughout the calculations, a commensurate value of Q was used. This fixed value was used since it was observed that at H = 0, (i) the HO instability criterion is dominated by electron-hole processes with energies greater than k B T and, hence, is insensitive to the precise form of the dispersion relations at the Fermi-energy, (ii) commensurate ordering produces greater condensation energies than the neighboring incommensurate phases due to the effect of multiple scattering processes and (iii) at finite fields, the HO instability involves equal spin electron-hole pairs and, thus, μ B H is expected to act as a low-energy cut-off in the gap equation much in the same way as k B T does. Figure 1 shows the specific heat C v obtained from the free energy given in Eq. (5). The specific heat is shown for different magnitudes of the magnetic field H z . The jump in the specific heat marks the temperature at which the transition to the HO phase occurs. For finite applied magnetic fields H z > 0.0 T, the jump in C v is shifted to lower temperatures and for H z ≳ 28.0 T the anomaly becomes sharper indicating that the system is approaching a regime of first-order transitions.
The first-order nature of the transition inferred from the specific heat shown in Fig. 1 is confirmed by the analysis of the free energy as function of the HO OP z χ 0 χ Q;σ . For the particular case in which U = J, a non-trivial solution for z χ 0 χ Q;σ must satisfy the condition J 2 F σ = 1 in which F σ = g σ g − σ . Figure 2a shows the quantity J 2 F σ versus z χ 0 χ Q;σ for different values of H z . For H z = 20.0 T; there is only one non-trivial solution (z Q,σ ≠ 0) s 1 indicating that the system is in a regime of second-order transitions. However, for H z = 33.0 T there are two non-trivial solutions, s 1 and s 2 , characterizing a regime of first-order transitions. Figure 2b shows the behavior of the free energy per unit cell (f = F/N) for H z = 32.0 T and for the temperature T = 2.0 K (<T HO = 3.48 K). The solution s 2 always corresponds to a maximum of the free energy, therefore, it  Figure 2c shows the free energy f for H z = 32.0 T and T = 5.0 K. Since T > T HO , the solution s 0 is the more stable one.
From the minima of the free energy, one can construct the T-H z phase diagram. The phase diagram is displayed in Fig. 3. The solid line marks second-order transitions between the HO (HO) and the paramagnetic (PM) phases. The dashed line at low temperature and high magnetic field marks first-order transitions. The dotted lines are the spinodal lines which mark the boundaries of the regions where multiple solutions for the OP z Q;σ are found (see Fig. 2a). The values of H z and T which describe the first-order and spinodal lines were determined from the free energy, which has forms similar to those shown in Fig. 2. The field scale that determines the switch from second-order to first-order is determined by twice the energy difference between the energy at which the commensurate nesting would be optimal and the actual value of the Fermi-energy. At this field, the condensation energy due to the commensurate scattering no longer sufficiently balances the entropy change. Since the mean-field approximation is not expected to be valid within the Ginzberg critical regime, the effects of commensurate and incommensurate fluctuations should be investigated further.   z-axis, T HO decreases when H z increases. By contrast, T HO is only weakly affected when H x is applied. The linear susceptibility is obtained as the second-derivative of the free energy. Figure 5 shows χ 1 for both the H x and H z magnetic fields. The derivatives were performed numerically. Our results shows that, at low temperatures, χ 1,z has a weaker temperature dependence than χ 1,x . However, as T → T HO , χ 1,z increases rapidly in such way that, at T = T HO , χ 1,z = χ 1,x . This temperature behaviour is in agreement with experimental results 57 for URu 2 Si 2 . The green region in Fig. 5 shows the difference χ 1,x − χ 1,z stressing the Fig. 6 The spin-up quasi-particle bands for α electrons at T = 0 and H z = 20 T for the direction A À Γ À A with A ¼ ðÀπ; Àπ; ÀπÞ, Γ = (0, 0, 0) and A = (π, π, π). The inset shows in detail the gap between the folded β-band and the hybridized α-band   A(ω, k). The panels c, f show in details the gapped region of the quasi-particle bands along the direction M-Γ Unusual magnetic field-dependence of a possible hidden EJ Calegari et al.
dependence of the linear susceptibility on the field orientation and temperature.

DISCUSSION
The results for T HO and χ 1 discussed above can also be explained by the analysis of the quasi-particle dispersion relations. Figure 6 shows the spin-up quasi-particle bands E ðjÞ " k ð Þ for the α electrons at T = 0.0 K and H z = 20 T over the entire diagonal of the Brillouin zone (−π, −π, −π) -(π, π, π). The circle highlights the gapped region that is shown in the inset. Figure 7 shows the polarized quasi-particle bands E ðjÞ σ k ð Þ for the α electrons, in the vicinity of the HO gap for temperatures close to T HO . The up arrows indicate the spin-up sub-bands represented by the red lines while the down arrows indicate the spin-down sub-bands denoted by the blue lines. The Fig. 7a, c, e shows the effect of fields oriented along the z-direction on the HO gap for T = 16.5 K. Inspection of the quasiparticle bands reveals that the behavior of T HO observed in Fig. 4 comes from the combined effects of the magnetic field and temperature. By increasing the field H z , the Zeeman splitting between the spin-up and spin-down sub-bands increases, while the closing of the gap is mainly caused by increasing temperature and, to a lesser extent, the raising of the magnetic field. We note that as the field is increased beyond H z = 10 T the up spin band falls below the Fermi-energy, indicating a change in Fermi-surface topology known as a Lifshitz transition. For H z = 15 T the gap between the spin-up and spin-down band has closed and the HO state has been destabilized. The initial slow field-dependence of the gap reflects the fact that, close to the (gapped) nesting point, the contributions of the Zeeman energy to the free energy cancel exactly. This cancellation occurs as long as the upper spin-up band and the lower spin-down band do not cross the Fermi-energy and produce electron and hole pockets. Hence, the dominant initial H 2 z -dependence of the gap originates from the states with other k values at which the reconstructed bands cross the Fermi-energy. However, since these other states are only weakly coupled to the states at the nesting points k ¼ ± Q/2, the initial field-dependence of the gap is small but accelerates when electron and hole pockets are formed. Likewise, since the resulting H 2 z term in the free energy originating from the quasi-particle bands only comes from a restricted range of k, the linear susceptibility χ 1,z is correspondingly small. The Fig. 7b, d, f shows the change in the quasi-particle bands caused by a field oriented along the x-direction when T = 16.5 K. By contrast to the H z case, the main effect of the H x field is to not to produce a Zeeman splitting but instead produces a spindependent momentum-shift of the quasi-particle bands. The spindependent momentum shifts are produced by H x which couples the spin-up and spin-down states that are comprised of spindependent admixtures of states with k and k ± Q. Therefore, to second-order in H x , the field produces a lowering of the quasiparticle energies that is most prominent where the admixture associated with the HO OP is large. Thus, as seen in Fig. 7, the magnitudes of the gaps Δ ↑ and Δ ↓ and T HO are only slightly dependent on H x . As a result, the bottom of the up-spin band reaches the Fermi-energy at practically the same temperature as at H x = 0.0 T. Therefore, the formation of electron and hole pockets are dominantly caused by temperature variation and hardly depend of H x . The above observations explain the distinct behaviors of T HO seen in Fig. 4. Likewise, since the energy lowering from the H x field is proportional to H 2 x and since it comes from the both the gapped and non-gapped regions of k space, the linear susceptibility χ 1,x is greater than χ 1,z .
The Fermi surfaces for H z = H x = 0 T are shown for two different temperatures (above and below T HO = 17.5 K) in Fig. 8. In 8a, T = 18.0 K, the black lines show the Fermi surface for the α electrons, while the red line shows the Fermi surface for the β electrons. In Fig. 8b, the quasi-particle bands for the α (black) and β (red) electrons are shown for the same temperature as in Fig. 8a. The color map shows the spectral function A(ω,k) for each one of the quasi-particle bands. Figure 8c shows in detail the region where the quasi-particle bands intercept the chemical potential at the points 2, 3 and 4 along the direction M-Γ. Point 2 is associated to the Fermi surface of the β electrons while points 3 and 4 correspond to the Fermi surfaces of the α electrons (see Fig. 8a). Figure 8d show the Fermi surfaces for HO state with T = 16.5 K, while Fig. 8e exhibits the quasi-particle bands and the spectral function for the same temperature as in Fig. 8d. The comparison of the Fig. 8a, d show that when the system enters in the HO state the Fermi surface is reduced significantly. Above T HO (see Fig. 8a) there are three branches of Fermi surface, two from the hybridized α − f and d band and one from the unhybridized β − f band. Below T HO (see Fig. 8d), the hybridized band α − f with wave vector k mixes with the unhybridized band β − f at k þ Q as shown in Fig. 8e. Two bands become gapped as in Fig. 8f and, therefore, there are only two branches of Fermi surface in the HO state. Moreover, the spectral weight of one of the branches (the red one) is very small as shown at the point 1 of Fig. 8f. Therefore, it can be concluded that about 60% of the Fermi-surface has been lost when the system enters in the HO state.
The reconstruction of the Fermi surface in applied field is shown in the field-induced changes of the Fermi-surface topology are Fig. 9 In a, c, e, a cross-section of the spin-up Fermi surfaces in the k z = 0 plane for the same magnetic fields (H z ) as in Fig. 7. In b, d, f, the spin-down Fermi surfaces. The temperature and the model parameters are the same as in Fig. 7 Unusual magnetic field-dependence of a possible hidden EJ Calegari et al. seen in Fig. 9. This process can be better understood through the analysis of the quasi-particle bands. Figure 10 show the quasiparticle bands for a small range of the wave vector k along the direction M-Γ. In Fig. 10a, b, the situation is the same as in the case of H z = 0 and temperature below T HO , i.e., two bands become gaped, therefore there are only two branches of Fermi surface (see Fig. 9a, b). Figure 10c show the spin-up quasi-particle bands for H z = 10.0 T. Notice that the magnetic field H z shifted the bands to lower energies and as a consequence the bands are intercepted by the chemical potential also at points 2 and 3, in addition to points 1 and 4. Therefore, there are four branches for the spin-up Fermi surface as shown in Fig. 9c, but, taking into account that the branch associated to the point 1 has very small spectral weight, it can be considered that only three branches remain as in the normal state. On the other hand, in Fig. 10d the bands are shifted up by the magnetic field, but not sufficiently to intercept the chemical potential at the points 2 and 3. Thus the spin-down Fermi surface for H z = 10.0 T has only two branches as shown in Fig. 9d. The comparison of Fig. 9c, d shows that under an applied magnetic field, due to the band polarization, the reconstruction of the Fermi surface is asymmetric, that is, the spin-up Fermi surface reconstructs before the spin-down Fermi surface. For H z = 15.0 T the magnetic field suppresses the HO state and the structure of the quasi-particle bands of the normal state is recomposed as shown in Fig. 10e, f, and therefore the Fermi surface is completely reconstructed (see Fig. 9e, f). However, due to the polarization of the bands by the magnetic field, the spin-up Fermi surface is more populated than the spin-down one, consequently, the spin-up Fermi surface has a volume greater than the spin-down Fermi surface. To summarize, we have found the application of a magnetic field to the HO state may result in a cascade of changes in the topology of the Fermi surface. The sequence of field-induced Lifshitz transitions has been reported in the literature. 59,60 In addition to describing the above effects, the theory predicts that the density of states (DOS) measured in scanning tunneling microscope (STM) tunneling experiments in field should show marked differences when the field is applied either parallel or perpendicular to the director. The change in the DOS is intimately related to the dependence of the quasi-particle bands on the orientation of the field. Our zero-field calculated DOS shows asymmetric peaks around the partial gap caused by the HO, in accordance with measurements. 61 For fields along the hard-axis, our results show that the asymmetric peaks will spin-split, thereby lowering the Kramer's degeneracy. On the other hand, fields in the easy-plane should not produce spin-splitting and the gap structure should closely resembles the H = 0 structure. This type of STM experiment requires a highly-symmetrical tip and the capability of rotating either the field in the a-b plane or the sample's surface and, therefore, is quite demanding.
To complete our discussion, in Fig. 11a, b are displayed, respectively, the non-linear susceptibility χ 3 and the ratio χ 3 /χ 1 for Unusual magnetic field-dependence of a possible hidden EJ Calegari et al.
both z and x axes. Our results show that the χ 3,z and χ 3,x have, respectively, a large sharp discontinuity and a very small step at T = T HO . This behavior can be traced directly from χ 1 (see Fig. 5). From a broad view, such anisotropies are ultimately related to the way that the orientation of the magnetic field affects the OP introduced in Eq. (1). That can be seen quite clearly in the behaviour of T HO shown in Fig. 4. From that result, one can conclude that there is a sharp discontinuity in χ 3 wherever T HO is dependent on the field. On the other hand, such discontinuity tends to disappear in the direction of the field wherever T HO has weak dependence on the field.
We have performed calculations on the UAL model that may describe the HO phase of URu 2 Si 2 . The magnetic properties, such as the susceptibility, become anisotropic below T HO . We have shown that, for fields applied along a spontaneously chosen hard axis, T HO decreases towards zero and that the HO transition changes from second order to first order at a large value of the magnetic field. This result is consistent with experimentally determined specific heat of URu 2 Si 2 . 55,56 These features are consequences that the transition, driven by the spin-flip part of the Hund's rule interaction, spontaneously breaks the spinrotational invariance. Furthermore, since the ordering involves the x and y components of the Hund's rule exchange, the ordering produces a preferential direction. Therefore, below T HO , the quasiparticle bands depend on the orientation of the magnetic field. This feature produces the anisotropy in the system's properties, such as the magnetic susceptibility, which qualitatively agrees with experimental results 57 and gives support to the OP proposed 37 to describe the HO phase of URu 2 Si 2 . Results for the Fermi surface show that, under an applied field H z , the reconstruction of the Fermi surface is anisotropic due to the band polarization. Finally, the non-linear susceptibility χ 3 exhibits a discontinuity at T = T HO , similar to that observed in the specificheat. The discontinuity is strongly dependent on the orientation of the magnetic field.
To conclude, the missing ingredient in our theory is the Ising anisotropy, which is not a unique feature of the HO phase since it already appears in the PM phase and large moment AFM phases. 62 Our simplified model cannot describe the anisotropy, without the addition of spin-orbit coupling and other ingredients such as the introduction of the set of Slater-Koster parameters t σ , t π , t δ and t ϕ that describe the anisotropic hopping of the 5f electrons. Nevertheless, even while not providing a detailed description of all the properties of URu 2 Si 2 , our results suggest that a variety of experimental results on this compound can be qualitatively understood in terms of one mechanism. Specifically, the spontaneous breaking of spin-rotational invariance driven by the Hund's rule interaction can qualitatively account for the unusual features found in experiments.

METHODS
All the numerical results were computed within a mean-field treatment of Coulomb and Hund's rule exchange interactions in the UAL Model. The conduction band and the two 5f bands were treated in the tight-binding approximation and the hybridization was treated as being independent of k. The HO gap was evaluated from a self-consistency equation by an iterative method. The first-order phase boundary was determined by a direct comparison of the free energies of the two phases.

Data availability
The authors declare that all data supporting the findings of this study are available from the authors. Fig. 11 a χ 3z and χ 3x as function of the temperature. b (χ 3,z )/(χ 1,z ) as a function of the temperature Unusual magnetic field-dependence of a possible hidden EJ Calegari et al.