Nematic order in iron superconductors – who is in the driver’s seat?
Although the existence of nematic order in iron-based superconductors is now a well-established experimental fact, its origin remains controversial. Nematic order breaks the discrete lattice rotational symmetry by making the and directions in the Fe plane non-equivalent. This can happen because of (i) a tetragonal to orthorhombic structural transition, (ii) a spontaneous breaking of an orbital symmetry, or (iii) a spontaneous development of an Ising-type spin-nematic order – a magnetic state that breaks rotational symmetry but preserves time-reversal symmetry. The Landau theory of phase transitions dictates that the development of one of these orders should immediately induce the other two, making the origin of nematicity a physics realization of a “chicken and egg problem”. The three scenarios are, however, quite different from a microscopic perspective. While in the structural scenario lattice vibrations (phonons) play the dominant role, in the other two scenarios electronic correlations are responsible for the nematic order. In this review, we argue that experimental and theoretical evidence strongly points to the electronic rather than phononic mechanism, placing the nematic order in the class of correlation-driven electronic instabilities, like superconductivity and density-wave transitions. We discuss different microscopic models for nematicity in the iron pnictides, and link nematicity to other ordered states of the global phase diagram of these materials – magnetism and superconductivity. In the magnetic model nematic order pre-empts stripe-type magnetic order, and the same interaction which favors nematicity also gives rise to an unconventional superconductivity. In the charge/orbital model magnetism appears as a secondary effect of ferro-orbital order, and the interaction which favors nematicity gives rise to a conventional superconductivity. We explain the existing data in terms of the magnetic scenario, for which quantitative results have been obtained theoretically, including the phase diagram, transport properties of the nematic phase, scaling of nematic fluctuations, and the feedback of the nematic order on magnetic and electronic spectra.
The discovery of iron-based superconductors (FeSCs) with transition temperatures as high as K has signaled the beginning of a new era in the investigation of unconventional superconductivity (for a review, see review ()). The key first step to unveil the nature of the superconducting phase is to understand the normal state from which superconductivity arises. In most FeSCs, superconductivity is found in proximity to a magnetically ordered state (transition temperature ), which led early on to the proposal that magnetic fluctuations play the key role in promoting the superconducting pairing review_Mazin (); review_Andrey (). A more careful examination of the phase diagram, however, revealed that there is another non-superconducting ordered state besides magnetism. Namely, at a certain temperature , the system spontaneously breaks the symmetry between the and directions in the Fe plane, reducing the rotational point group symmetry of the lattice from tetragonal to orthorhombic, while time-reversal symmetry is preserved. In some materials, such as hole-doped (BaK)FeAs, the tetragonal-to-orthorhombic and magnetic transitions are simultaneous and first-order (), whereas in electron-doped Ba(FeCo)As and isovalent-doped BaFe(AsP), they are split () and second order Kim11 (); Birgeneau11 (); matsuda_t () (see Fig. 1). As doping increases, the line tracks the line across the phase diagram, approaching the superconducting dome. It is therefore essential to understand the origin of this new order as it may support or act detrimentally to superconductivity.
The order parameter for a transition in which a rotational symmetry is broken but time-reversal symmetry is preserved is a director (i.e. a vector without an arrow), similar to the order parameter in the nematic phase of liquid crystals Fradkin_review (). By analogy, the orthorhombic state in FeSCs has been called a “nematic state”. Unlike isotropic liquid crystals, however, in FeSCs the lattice symmetry forces the director to point only either along or directions, what makes the nematic order parameter Ising-type (Ising-nematic).
At first sight, one might view this tetragonal-to-orthorhombic transition as a regular structural transition driven by lattice vibrations (phonons). However, experiments find anisotropies in several electronic properties, such as the dc resistivity Chu10 (); Tanatar10 (), to be much larger than the anisotropy of the lattice parameters. This led to the idea that the tetragonal-to-orthorhombic transition may be driven by electronic rather than lattice degrees of freedom. If this is the case, then the transition into the nematic phase is driven by the same fluctuations that give rise to superconductivity and magnetic order, and therefore is an integral part of a global phase diagram of FeSCs. Electronic nematic phases have been recently proposed in other unconventional superconductors, such as high- cuprates and heavy-fermion materials Fradkin_review (). An electronically driven nematic state in FeSCs would be in line with a generic reasoning that the pairing in all these correlated electron systems has the same origin.
The discussion on the “nematicity” in FeSCs has been largely focused on two key issues: (i) Can the experiments distinguish “beyond reasonable doubt” between phonon-driven and electron-driven tetragonal symmetry breaking? (ii) If this transition is driven by electrons, which of their collective degrees of freedom are driving it - charge/orbital fluctuations or spin fluctuations? Answering the last question is crucial for the understanding of superconductivity in FeSCs because we argue below that charge/orbital fluctuations favor a sign-preserving -wave state () whereas spin fluctuations favor a sign-changing -wave () or a -wave state. Here we give our perspective on these issues, discuss the phenomenology of the nematic state, its experimental manifestations, and the underlying microscopic models.
Ii Phenomenology of the nematic phase
To describe the nematic state, the first task is to identify the appropriate order parameter. The experimental manifestations of nematic order can be clustered into three classes. Taken alone, each class points to a different origin of the nematic phase (see Fig. 2 for schematic representation):
Structural distortion – the lattice parameters and along the and directions become different Kim11 (). Such an order is normally associated with a phonon-driven structural transition;
Charge orbital order – the occupations and (and on-site energies) of the and Fe-orbitals become different Yi2011 (). The appearance of such an order is normally associated with divergent charge fluctuations;
Spin-nematic order – the static spin susceptibility becomes different along the and directions of the Brillouin zone before a conventional SDW state develops matsuda_t (). The appearance of such an order is normally associated with divergent quadrupole magnetic fluctuations.
The fact that these three order parameters are non-zero in the nematic phase leads to a dilemma, which can be best characterized as the physics realization of a “chicken and egg problem”: all three types of order (structural, orbital, and spin-nematic) must be present no matter who drives the nematic instability. This follows from the fact that bi-linear combinations of the order parameters which break the same symmetry (in our case, the tetragonal symmetry of the system) are invariant under symmetry transformations and must therefore appear in the Landau free energy. Suppose that one of the three order parameters is the primary one, i.e. its fluctuations drive the nematic instability. Let’s call it and the other two and . The free energy has the generic form
Because the nematic transition is driven by , the coefficient , which corresponds to the order parameter susceptibility in the disordered state, diverges at and becomes negative for , while and remain finite and positive (although fluctuations of and may shift slightly ). For , orders on its own: . If in (1) were zero, the other two fields and would not order, but once are finite, a non-zero instantly induces finite values of the secondary order parameters , . As a consequence, there is only one nematic transition temperature at which all three become non-zero (e.g., lattice symmetry is broken at the same temperature where electronic nematic order emerges), and it is not possible to determine who causes the instability by looking solely at equilibrium order parameters. An additional experimental complication is the presence of nematic twin domains below , what effectively averages to zero. This problem can be circumvented by applying a small detwinning uniaxial stress Chu10 (); Tanatar10 (), which acts as a conjugate field to and breaks the tetragonal symmetry at all temperatures, making an ill-defined quantity.
One way to select the primary order is to carefully study fluctuations in the symmetry-unbroken phase at . Because the primary order parameter acts as an external field for the secondary order parameters, and , fluctuations of the former renormalize the susceptibilities of the latter to
where is the susceptibility of the primary field. The renormalized susceptibilities of the secondary fields do diverge at the nematic transition, however for small enough and , and begin to grow only in the immediate vicinity of , where is already large. If one can measure the three susceptibilities independently, Eq. (2) in principle provides a criterion to decide which order parameter drives the instability. The implementation of this procedure is possible (see next section), but is complicated by two factors. First, it only works if and are relatively weak, what normally implies that the systems falls into the weak/moderate coupling category. If the coupling is large, all three order parameters become so inter-connected that the question “who is in the driver’s seat?” becomes meaningless. Second, in some FeSCs the nematic transition is first order, in which case all three susceptibilities jump from one finite value to another, even before the susceptibility of the primary field gets enhanced.
Iii Experimental evidence for electronic nematicity
iii.1 Measurements in the nematic phase
The first evidence for the electronic character of the tetragonal-to-orthorhombic transition came from resistivity measurements in detwinned samples Chu10 (); Tanatar10 (), which revealed that resistivity anisotropies are significantly larger than relative lattice distortions and also display a nontrivial dependence on doping and disorder Nakajima2011 (). Other non-equilibrium quantities, such as thermopower Jiang2013 () and optical conductivityDusza2011 (); Nakajima2011 (), were also found to display large anisotropies, which in optical measurements were observed to extend to energies of several hundreds of . Anisotropies in observables related to charge and spin were also seen: angle-resolved photoemission spectroscopy (ARPES) found a splitting between the on-site energies of the and orbitals, indicative of ferro-orbital order Yi2011 () and torque magnetometry revealed different uniform magnetic susceptibilities along the and directions matsuda_t (). The onset of magnetic anisotropy coincides with the observation of a non-zero orthorhombic distortion, in agreement with the discussions of the previous section. Strong signatures of emerging magnetic anisotropy were also found in the behavior of the nuclear magnetic resonance (NMR) lines across Imai12 ().
The direct observation of electronic anisotropy in the nematic state was made possible by scanning tunneling microscopy (STM). The first measurements, performed deep inside the magnetic phase, found that the local density of states around an impurity is characterized by a dimer-like structure extended along the magnetic ordering vector direction Davis10 (). Subsequent measurements showed that these dimers persist above , in the temperature regime of the nematic state Rosenthal13 (). An additional piece of evidence in favor of the electronic character of the nematic transition came, ironically, from x-ray measurements of the orthorhombic distortion inside the SC phase. These measurements found a strong suppression of the distortion below Nandi10 (), what is a characteristic signature of the competition for the same electronic states between two electronically-driven orders.
iii.2 Measurements in the tetragonal phase
A few recent measurements focused on fluctuations in the tetragonal state, in particular, on the shear modulus , which is the inverse susceptibility of the structural order parameter Fernandes2010 (); Yoshizawa2012 (); Boehmer2013 (). If the structural transition is driven not by the lattice but by some other electronic degree of freedom, Eq. (2) provides a natural way to connect the shear modulus to the electronic nematic susceptibility . An experimentally observed softening of the shear modulus above was successfully fitted by Eq. (2) using both magnetic Fernandes2010 () and charge/orbital Yoshizawa2012 () phenomenological models for , indicating that structural distortion is very likely not the primary order.
Perhaps the strongest evidence that the nematic transition is electronically-driven came from the recent measurements of the anisotropy of the resistivity Chu2012 (). Using a piezoelectric, the measurements were performed by using strain (the structural distortion) as the control parameter, rather than stress, as in previous setups. The strain is one of the order parameter fields in the free energy Eq. (1). Using the resistivity anisotropy as a proxy of the nematic order parameter, it was experimentally shown that the susceptibility diverges near the nematic transition. This is only possible if the structural distortion is a conjugate field to the primary order parameter, rather than the primary order parameter itself – otherwise would be simply proportional to the order parameter , with a constant prefactor.
Iv Microscopic models for electronic nematicity
A successful microscopic theory for electronic nematicity must describe the global phase diagram of FeSCs, i.e. not only the nematic order but also magnetism and superconductivity. A popular starting point is the multi-orbital Hubbard model, which describes hopping between all Fe-As orbitals and local interactions, such as intra-band and inter-band Hubbard repulsions and Hund’s exchange review_Mazin (). There is a general agreement among researches that this model does contain all information about the phase diagram. The model has been analyzed at both weak/intermediate coupling, when the system is a metal, and at strong coupling, when electrons on at least some orbitals were assumed to be localized or “almost localized”. The nematic order has been obtained in both limits, what is yet another indication that it is a generic property of FeSCs. We adopt the itinerant approach, since most FeSCs are metals. In this itinerant scenario, the low-energy electronic states lie around hole-like Fermi-surface pockets at the center of the Fe-square lattice Brillouin zone and electron-like Fermi-surface pockets at the borders of the Brillouin zone, see Fig. 5a. The microscopic reasoning for either magnetic or orbital scenarios of electronic nematicity follows from two different assumptions about the sign of the effective inter-pocket interaction Eremin10 (), which is a combination of the Hubbard and Hund interactions dressed up by coherence factors associated with the transformation from the orbital to the band basis. As we will see, each scenario leads to a prediction of a particular superconducting pairing state.
iv.1 Magnetic scenario
The magnetic mechanism for the nematic order follows from the observation that in most FeSCs the observed magnetic order on the Fe atoms is of stripe type, with ordering vectors or i.e. spins are parallel to each other along one direction and anti-parallel along the other Dai_review () (see Fig. 3a). This order breaks not only the spin-rotational symmetry (and time-reversal symmetry), but it also breaks the lattice rotational symmetry down to by choosing the ordering vectors to be either or . This additional tetragonal symmetry breaking enhances the order parameter manifold to Sachdev (); Kivelson (). In terms of the two magnetic order parameters and , associated with the ordering vectors and , the breaking of the symmetry implies while the breaking of the symmetry implies Fernandes12 (). In a mean-field approach both and symmetries are broken simultaneously at . However, fluctuations split the two transitions and give rise to an intermediate phase at where tetragonal symmetry is broken but the spin-rotational symmetry is not, i.e. while . This is by definition a nematic order, which, viewed this way, is an unconventional magnetic order which preserves time-reversal symmetry (a spin nematic). In real space, the stripe magnetic state can be viewed as two inter-penetrating Neel sublattices with staggered magnetizations and . In terms of these quantities, the nematic state is characterized by while (see Fig. 3).
Within a microscopic description, the instability towards a stripe magnetic order is associated with the divergence of the static spin susceptibility . Without any interactions, the bare particle-hole susceptibility is by itself sizable at and because these wave-vectors connect electronic states at the hole and electron pockets. When the inter-pocket interaction is positive (repulsive), there is an additional RPA-type enhancement of the spin susceptibility, roughly as , and at some , diverges. This however does not guarantee that the magnetic order is of stripe type as the latter emerges only if below and , or vise versa. To determine which magnetic state develops, one needs to calculate higher order terms in the magnetic free energy Fernandes12 (); Eremin10 (); Brydon11 (). The result is that at least for small dopings, the system selects the stripe order. The static nematic susceptibility (the correlator of ) can be obtained by including fluctuations of the nematic order parameter , yielding:
where is the temperature, and is the composite coupling which, when positive, sets the magnetic order to be of stripe type. In dimensions , diverges at (assuming that the magnetic transition is second order). Eq. (3) then shows that the nematic susceptibility diverges at a higher , when , i.e. at sufficiently large but still finite magnetic correlation length. This mechanism naturally ties the nematic and magnetic ordering temperatures to each other over the entire phase diagram. In between and , the symmetry is broken but is not, i.e., but . The difference between and is stronger in quasi-2D systems where is further decreased by thermal fluctuations, while remains unaffected Kivelson (); batista ().
More detailed microscopic calculations show that for some system parameters the nematic transition is second order, but for other input parameters it becomes first-order Fernandes12 (). In the latter case, a jump in the nematic order parameter induces a jump in the magnetic correlation length, which may instantaneously trigger a first-order magnetic transition. In any case, when the Fermi pockets are decomposed into their orbital characters, one finds within the same microscopic model that the emergence of spin-nematic order gives rise to orbital order , since the electron pocket at has mostly character, whereas the electron pocket at has mostly character. Similarly, a spin-nematic order induces a structural distortion Qi09 (); Cano10 ().
We see therefore that the repulsive inter-pocket interaction enhances spin fluctuations, which gives rise to both magnetism and nematicity. To describe the global phase diagram of FeSCs, one needs also to investigate superconductivity. Spin fluctuations peaked at and strongly enhance inter-pocket repulsion, which becomes larger than intra-pocket repulsion. In this situation, the system is known to develop either an unconventional superconductivity, in which the gap functions have different signs in the hole and in the electron pockets, or a superconductivity review_Mazin (); review_Andrey (). We emphasize that spin-nematic order and superconductivity are both intrinsic consequences of the same magnetic scenario.
Other microscopic models also find nematic order in proximity to a magnetic instability. For instance, explicit evaluation of the ferro-orbital susceptibility using the multi-orbital Hubbard model finds that it is enhanced only in the presence of spin fluctuations Kontani12 (), similarly to what is described by Eq. (3). Studies of models with both localized and itinerant orbitals also found Phillips11 (); Dagotto13 (); w_ku10 () that the proximity to magnetism is an important ingredient for orbital order. In purely localized-spin models the interplay between magnetism and ferro-orbital order is blurred by the complicated form of the effective Hamiltonian, which deviates from a simpler Kugel-Khomskii type kruger09 (); Applegate11 ().
iv.2 Charge/orbital scenario
In its simplest form, the charge/orbital scenario for the nematic order parallels the magnetic scenario, the only difference being the sign of the interaction between electron and hole pockets. If this interaction turns out to be negative, it is the charge/orbital susceptibility rather than the spin susceptibility that is enhanced as , diverging at and at a certain . This divergence would signal the onset of a charge density-wave state with ordering vectors or (or both) and order parameters and . This order breaks translational symmetry and, like in the magnetic scenario, breaks also an additional symmetry if only one order parameter becomes non-zero. It is natural to expect, although no explicit calculations have been done to the best of our knowledge, that fluctuations split the temperatures at which the translational and the symmetries are broken, in a manner similar to Eq. (3). Then, in the intermediate temperature range , the system spontaneously develops ferro-orbital order in which while . A structural distortion and the difference between and appear instantly once ferro-orbital order sets in. However, magnetic order only appears at a smaller temperature, presumably via changes in the magnetic correlation length induced by the ferro-orbital order.
For the Cooper pairing, the orbital scenario implies that the inter-pocket interaction is attractive and enhanced. Once this interaction exceeds the intra-pocket repulsion, the system develops a superconducting instability towards an state – a conventional pairing state where the gap functions have the same sign in all pockets review_Andrey ().
What we described above is the simplest scenario for orbital order. More complex models have been also proposed to account for the nematic transition without involving magnetic degrees of freedom. In Ref. Littlewood (), it was proposed that nematicity could arise as an unequal hybridization between localized orbitals and itinerant orbitals. In Ref. Tesanovic11 () it was suggested that both spin and charge interactions are present and that the larger interaction in the spin channel gives rise to magnetic order at, say, . However, before this happens, a weaker charge interaction gives rise to charge order at the other momentum (a pocket density-wave state), which would break the tetragonal symmetry of the system. Whether such a pocket density-wave is experimentally realized in FeSCs remains to be seen.
V Comparison between theory and experiment
Although the experimental evidence presented in Section III favors an electronic nematic instability, disentangling the orbital and magnetic scenarios is difficult on a qualitative level, what begs for a more direct comparison between microscopic models and experimental results. In this regard, the doping evolution of the magnetic and structural transitions is an important benchmark. BaFeAs, one of the compounds most extensively investigated, displays a second-order nematic transition at followed by a a first-order “meta-nematic transition” at a lower , where the system simultaneously undergoes a first-order magnetic transition. The meta-nematic transition has been observed by x-ray Kim11 (); Birgeneau11 () and torque magnetometry matsuda_t (), although the data disagree on the precise value of . As charge carriers are introduced in the system via Co substitution in the Fe sites, the splitting between the two transitions increases, and eventually the meta-nematic transition disappers and the magnetic transition becomes second-order.
How does this compare to theory? For the magnetic scenario, a detailed theoretical analysis Fernandes12 () shows that three types of system behavior are possible in systems that are moderately anisotropic, depending on the value of the nematic coupling (see Fig. 4). At large , nematic and magnetic transitions are simultaneous and first order. At intermediate , nematic order develops via a second-order transition, and there is a meta-nematic transition at a lower , where magnetic order also develops discontinuously. At smaller , nematic and magnetic transitions are separate and second-order, with an intermediate spin-nematic phase between and . The microscopic calculations found that decreases with electron doping, and the theoretical phase diagram in Fig. 4 is fully consistent with the one for the electron-doped Ba(FeCo)As if we place the point in the region II in Fig. 4. No calculations of how the nematic and magnetic transitions evolve with carrier concentration have been done within the charge/orbital scenario.
One can take the comparison with the data even further and compare the two versions of the magnetic scenario – for itinerant and for localized spins. Both predict stripe magnetic order and pre-emptive symmetry-breaking but differ in the details. In particular, in localized models is generally small and is unaffected by carrier concentration Chandra90 (); Kivelson (). This makes the description of the doping dependence in the localized spin approach somewhat problematic, although not impossible batista (). A more essential difference is that in localized models the coupling is always positive, while in itinerant models may become negative at large enough hole doping Lorenzana08 (); Brydon11 (). For negative , there is no tetragonal symmetry breaking either above or below the magnetic transition as the system selects a tetragonally-symmetric combination of both and orders. A symmetry-preserving magnetic state with orders at and has been recently observed in BaNaFeAs Osborn13 () and Ba(FeMn)As Kim10 () at large enough doping – a strong argument in favor of the itinerant magnetic scenario.
Another key quantity to compare experiment and theory is the resistivity anisotropy. Deep in the magnetically ordered state, the anisotropic folding of the Fermi surface plays the major role in determining the resistivity anisotropy Tanatar10 (); bascones (). In the nematic state, , orbital order and spin-nematic order have different effects on the dc resistivity: while the former causes an anisotropy in the Drude weight Devereaux10 (); Phillips11 (), the latter gives rise to anisotropy in the scattering rate Fernandes11 (). The calculated anisotropy in the Drude weight has the opposite sign to the one observed experimentally Devereaux10 (); Phillips11 (), whereas the calculated anisotropy in the magnetic scattering rate was shown to agree with experiments, including a sign-change of the anisotropy between electron-doped and hole-doped materials Blomberg13 ().
One can also compare theoretical and experimental results for the feedback effects from the nematic order on the electronic and the magnetic spectrum Phillips11 (); Fernandes12 (); Goswami11 (). In the magnetic scenario, nematic order enhances the magnetic correlation length, what gives rise to strong magnetic fluctuations and a possible pseudogap in the electronic spectrum. A significant increase of magnetic fluctuations below has been observed via NMR in compounds where and are well separated Ma_NaFeAs (). Also, recent ARPES experiments found the pseudogap behavior (a suppression in the density of states at low energies) whose onset coincides with the nematic transition pseudogap_Matsuda (). Within the orbital scenario, the key feedback from the orbital order is a Pomeranchuk distortion of the Fermi surface induced by orbital order Phillips11 ().
Nematic fluctuations above have also been used to compare experiment and theory. Orbital fluctuations have been argued to affect the density of states at the Fermi levelPhillips12 () and leave signatures in point-contact spectroscopy consistent with the data Greene12 (). Alternatively, one can employ Eq. (2) to compare the renormalized lattice susceptibility (the shear modulus ), assumed to be non-critical, with the susceptibility associated with either the orbital or the spin-nematic order parameter. Eq. (2) must be satisfied if the corresponding electronic order drives the nematic instability. In Ref. Gallais13 (), a quasi-elastic peak in the Raman response was attributed to charge/orbital fluctuations and used to extract the corresponding orbital susceptibility. On the other hand, the spin-nematic susceptibility, being proportional to (see Eq. 3), can be measured via the NMR spin-lattice relaxation rate . Comparison with shear modulus data for a family of electron-doped FeSCs found that there is a robust scaling between and data Fernandes13_shear (). This provides strong support to the idea that the nematic transition is magnetically-driven.
The bulk of experimental and theoretical results which we presented in this mini-review supports the idea that nematic order in FeSCs is of electronic origin, what places it at par with other known electronic instabilities such as superconductivity or density-wave orders. It is likely that magnetic fluctuations drive the nematic instability. In any case, all three orders – spin-nematic, orbital, and structural, appear simultaneously below . The important question not addressed until very recently is the role of nematicity for high-temperature superconductivity. It is unlikely that nematic fluctuations can mediate superconductivity as spin or charge fluctuations do, but nematic fluctuations may nevertheless enhance by reducing the bare intra-pocket repulsion. Below , however, nematic order has been found to compete with superconductivity Nandi10 (); Moon12 (), like density-wave orders do. A special case in which nematicity strongly affects is when -wave and -wave superconducting instabilities are nearly degenerate, what was suggested to be the case for strongly hole-doped and strongly electron-doped FeSCs Taillefer12 (). In this situation nematic order leads to a sizeable enhancement of by lifting the frustration associated with the competing pairing states Livanas12 (); Fernandes_Millis13 (); DungHaiLee13 (). These results clearly point to the need of additional investigations of the interplay between nematicity and superconductivity.
We acknowledge useful discussions with E. Abrahams, J. Analytis, E. Bascones, A. Böhmer, J. van den Brink, P. Brydon, S. Bud’ko, P. Canfield, P. Chandra, P. Dai, M. Daghofer, L. Degiorgi, I. Eremin, I. Fisher, Y. Gallais, A. Goldman, A. Kaminski, V. Keppens, D. Khalyavin, M. Khodas, S. Kivelson, J. Knolle, H. Kontani, A. Kreyssig, F. Krüger, W. Ku, W.C. Lee, J. Lorenzana, W. Lv, S. Maiti, D. Mandrus, R. McQueeney, Y. Matsuda, I. Mazin, C. Meingast, A. Millis, R. Osborn, A. Pasupathy, I. Paul, P. Phillips, R. Prozorov, S. Sachdev, Q. Si, T. Shibauchi, L. Taillefer, M. Takigawa, M. Tanatar, M. Vavilov, P. Wölfle, and M. Yoshizawa. The authors benefited a lot from the discussions with our great colleague Z. Tesanovic who unexpectedly passed away last year. A.V.C. is supported by the Office of Basic Energy Sciences U.S. Department of Energy under the grant #DE-FG02-ER46900.
- (1) J. Paglione and R.L. Greene, “High-temperature superconductivity in iron-based materials”, Nature Phys. 6, 645 (2010).
- (2) P. J. Hirschfeld, M. M. Korshunov, and I. I. Mazin, “Gap symmetry and structure of Fe-based superconductors”, Rep. Prog. Phys. 74, 124508 (2011).
- (3) A. V. Chubukov, “Pairing Mechanism in Fe-Based Superconductors”, Annu. Rev. Cond. Mat. Phys. 3, 57 (2012).
- (4) S. Avci, O. Chmaissem, D.-Y. Chung, S. Rosenkranz, E. A. Goremychkin, J.-P. Castellan, I. S. Todorov, J. A. Schlueter, H. Claus, A. Daoud-Aladine, D. D. Khalyavin, M. G. Kanatzidis, and R. Osborn, “Phase Diagram of (BaK)FeAs”, Phys Rev B 85, 184507 (2012).
- (5) M. G. Kim, R. M. Fernandes, A. Kreyssig, J. W. Kim, A. Thaler, S. L. Bud’ko, P. C. Canfield, R. J. McQueeney, J. Schmalian, and A. I. Goldman, “Character of the structural and magnetic phase transitions in the parent and electron-doped BaFeAs compounds”, Phys. Rev. B 83, 134522 (2011).
- (6) C. R. Rotundu and R. J. Birgeneau, “First- and second-order magnetic and structural transitions in Ba(FeCo)As”, Phys. Rev. B 84, 092501 (2011).
- (7) S. Kasahara, H.J. Shi, K. Hashimoto, S. Tonegawa, Y. Mizukami, T. Shibauchi, K. Sugimoto, T. Fukuda, T. Terashima, A. H. Nevidomskyy, and Y. Matsuda, “Electronic nematicity above the structural and superconducting transition in BaFe(AsP)” Nature 486, 382 (2012).
- (8) E. Fradkin, S. A. Kivelson, M. J. Lawler, J. P. Eisenstein, and A. P. Mackenzie, “Nematic Fermi Fluids in Condensed Matter Physics”,Annu. Rev. Condens. Matter Phys. 1, 153 (2010).
- (9) J.-H. Chu, J. G. Analytis, K. De Greve, P. L. McMahon, Z. Islam, Y. Yamamoto, and I. R. Fisher, “In-Plane Resistivity Anisotropy in an Underdoped Iron Arsenide Superconductor”, Science 329, 824 (2010).
- (10) M. A. Tanatar, E. C. Blomberg, A. Kreyssig, M. G. Kim, N. Ni, A. Thaler, S. L. Bud’ko, P. C. Canfield, A. I. Goldman, I. I. Mazin, and R. Prozorov, “Uniaxial-strain mechanical detwinning of CaFeAs and BaFeAs crystals: Optical and transport study”, Phys. Rev. B 81, 184508 (2010).
- (11) M. Yi et al., “Symmetry-breaking orbital anisotropy observed for detwinned Ba(FeCo)As above the spin density wave transition”, PNAS 108, 6878 (2011).
- (12) M. Nakajima et al., “Effect of Co doping on the in-plane anisotropy in the optical spectrum of underdoped Ba(FeCo)As”, Phys. Rev. Lett. 109, 217003 (2012).
- (13) S. Jiang et al., “Thermopower as a sensitive probe of electronic nematicity in iron pnictides”, Phys. Rev. Lett. 110, 067001 (2013).
- (14) A. Dusza et al., “Anisotropic charge dynamics in detwinned Ba(FeCo)As”, Europhys. Lett 93, 37002 (2011).
- (15) M. Nakajima et al., “Unprecedented anisotropic metallic state in undoped iron arsenide BaFeAs revealed by optical spectroscopy”, PNAS 108, 12238 (2011).
- (16) M. Fu, D. A. Torchetti, T. Imai, F. L. Ning, J.-Q. Yan, and A. S. Sefat, “NMR Search for the Spin Nematic State in LaFeAsO Single Crystal”, Phys. Rev. Lett. 109, 247001 (2012).
- (17) C. Dhital, Z. Yamani, W. Tian, J. Zeretsky, A. S. Sefat, Z. Wang, R. J. Birgeneau, S. D. Wilson, “Effect of Uniaxial Strain on the Structural and Magnetic Phase Transitions in BaFeAs”, Phys. Rev. Lett. 108, 087001 (2012).
- (18) T.-M. Chuang, M. P. Allan, J. Lee, Y. Xie, N. Ni, S. L. Bud’ko, G. S. Boebinger, P. C. Canfield, and J. C. Davis, “Nematic Electronic Structure in the “Parent” State of the Iron-Based Superconductor Ca(FeCo)As”, Science 327, 181 (2010).
- (19) E. P. Rosenthal, E. F. Andrade, C. J. Arguello, R. M. Fernandes, L. Y. Xing, X. C. Wang, C. Q. Jin, A. J. Millis, and A. N. Pasupathy, “Visualization of Electron Nematicity and Unidirectional Antiferroic Fluctuations at High Temperatures in NaFeAs”, arxiv:1307.3526 (2013).
- (20) S. Nandi, M. G. Kim, A. Kreyssig, R. M. Fernandes, D. K. Pratt, A. Thaler, N. Ni, S. L. Bud’ko, P. C. Canfield, J. Schmalian, R. J. McQueeney, and A. I. Goldman, “Anomalous Suppression of the Orthorhombic Lattice Distortion in Superconducting Ba(FeCo)As Single Crystals”, Phys. Rev. Lett. 104, 057006 (2010).
- (21) R. M. Fernandes, L. H. VanBebber, S. Bhattacharya, P. Chandra, V. Keppens, D. Mandrus, M. A. McGuire, B. C. Sales, A. S. Sefat, and J. Schmalian, “Effects of nematic fluctuations on the elastic properties of iron arsenide superconductors”, Phys. Rev. Lett. 105, 157003 (2010).
- (22) M. Yoshizawa et al., “Structural quantum criticality and superconductivity in iron-based superconductor Ba(FeCo)As”, J. Phys. Soc. Jpn. 81 (2012) 024604
- (23) A.E. Böhmer et al., “Nematic susceptibility of hole- and electron-doped BaFe2As2 iron-based superconductors”, arxiv:1305.3515 (2013).
- (24) J.-H. Chu et al., “Divergent nematic susceptibility in an iron arsenide superconductor”, Science 337, 710 (2012).
- (25) H.-H. Kuo, S.C. Riggs, and I.R. Fisher, “Measurement of the elastoresistivity coefficients of the underdoped iron-arsenide Ba(FeCo)As”, Phys. Rev. B 88, 085113 (2013).
- (26) I. Eremin and A. V. Chubukov, “Magnetic degeneracy and hidden metallicity of the spin-density-wave state in ferrophictides”, Phys. Rev. B 81, 024511 (2010).
- (27) P. Dai, J. Hu, and E. Dagotto, “Magnetism and its microscopic origin in iron-based high-temperature superconductors”, Nature Phys. 8, 709-718 (2012).
- (28) C. Fang, H. Yao, W.-F. Tsai, J. Hu, and S. A. Kivelson, “Theory of electron nematic order in LaFeAsO”, Phys. Rev. B 77, 224509 (2008).
- (29) C. Xu, M. Muller, and S. Sachdev, “Ising and spin orders in the iron-based superconductors”, Phys. Rev. B 78, 020501(R) (2008).
- (30) R. M. Fernandes, A. V. Chubukov, J. Knolle, I. Eremin, and J. Schmalian, “Preemptive Nematic order pseudogap, and orbital order in the iron pnictides”, Phys. Rev. B 85, 024534 (2012).
- (31) Y. Kamiya, N. Kawashima, and C. D. Batista, “Dimensional crossover in the quasi-two-dimensional Ising-O(3) model”, Phys. Rev. B 84, 214429 (2011).
- (32) Y. Qi and C. Xu, “Global phase diagram for magnetism and lattice distortion of iron-pnictide materials”, Phys. Rev. B 80, 094402 (2009).
- (33) A. Cano, M. Civelli, I. Eremin, and I. Paul, “Interplay of magnetic and structural transitions in Fe-based pnictide superconductors”, Phys. Rev. B 82, 020408(R) (2010).
- (34) S. Onari H. and H. Kontani, “Self-consistent Vertex Correction Analysis for Iron-based Superconductors: Mechanism of Coulomb Interaction-Driven Orbital Fluctuations”, Phys. Rev. Lett. 109, 137001 (2012).
- (35) F. Krüger, S. Kumar, J. Zaanen, and J. van den Brink, “Spin-orbital frustrations and anomalous metallic state in iron-pnictide superconductors”, Phys. Rev. B 79, 054504 (2009).
- (36) R. Applegate, R. R. P. Singh, C.-C. Chen, and T. P. Devereaux, “Phase transitions in spin-orbital models with spin-space anisotropies for iron pnictides: Monte Carlo simulations”, Phys. Rev. B 85, 054411 (2012).
- (37) W. Lv and P. Phillips, “Orbitally and magnetically induced anisotropy in iron-based superconductors”, Phys. Rev. B 84, 174512 (2011).
- (38) S. Liang, A. Moreo, and E. Dagotto, “Nematic State of Pnictides Stabilized by Interplay between Spin, Orbital, and Lattice Degrees of Freedom”, Phys. Rev. Lett. 111, 047004 (2013).
- (39) C. C. Lee, W. G. Yin, and W. Ku, “Ferro-Orbital Order and Strong Magnetic Anisotropy in the Parent Compounds of Iron-Pnictide Superconductors”, Phys. Rev. Lett. 103, 267001 (2009).
- (40) V. Stanev and P. B. Littlewood, “Nematicity driven by hybridization in iron-based superconductors”, Phys. Rev. B 87, 161122(R) (2013).
- (41) J. Kang and Z. Tesanovic, “Theory of the valley-density wave and hidden order in iron pnictides”, Phys. Rev. B 83, 020505 (2011).
- (42) P. Chandra, P. Coleman, and A. I. Larkin, “Ising transition in frustrated Heisenberg models”, Phys. Rev. Lett. 64, 88 (1990).
- (43) J. Lorenzana, G. Seibold, C. Ortix, and M. Grilli, “Competing Orders in FeAs Layers”, Phys. Rev. Lett. 101, 186402 (2008).
- (44) P. M. R. Brydon, J. Schmiedt, and C. Timm, “Microscopically derived Ginzburg-Landau theory for magnetic order in the iron pnictides”, Phys. Rev. B 84, 214510 (2011).
- (45) S. Avci, O. Chmaissem, S. Rosenkranz, J. M. Allred, I. Eremin, A. V. Chubukov, D.-Y. Chung, M. G. Kanatzidis, J.-P. Castellan, J. A. Schlueter, H. Claus, D. D. Khalyavin, P. Manuel, A. Daoud-Aladine, and R. Osborn, “The Origin of Nematic Order in the Iron-Based Superconductors”, arXiv:1303.2647 (2013).
- (46) M. G. Kim, A. Kreyssig, A. Thaler, D. K. Pratt, W. Tian, J. L. Zarestky, M. A. Green, S. L. Bud’ko, P. C. Canfield, R. J. McQueeney, and A. I. Goldman,”Antiferromagnetic ordering in the absence of a structural distortion in Ba(FeMn)As”, Phys. Rev. B 82, 220503(R) (2010).
- (47) B. Valenzuela, E. Bascones, and M. J. Calderon, “Conductivity Anisotropy in the Antiferromagnetic State of Iron Pnictides”, Phys. Rev. Lett. 105, 207202 (2010).
- (48) C.-C. Chen, J. Maciejko, A. P. Sorini, B. Moritz, R. R. P. Singh, and T. P. Devereaux, “Orbital order and spontaneous orthorhombicity in iron pnictides”, Phys. Rev. B 82, 100504 (2010).
- (49) R. M. Fernandes, E. Abrahams, and J. Schmalian, “Anisotropic In-Plane Resistivity in the Nematic Phase of the Iron Pnictides”, Phys. Rev. Lett. 107, 217002 (2011).
- (50) E. C. Blomberg, M. A. Tanatar, R. M. Fernandes, I. I. Mazin, B. Shen, H.-H. Wen, M. D. Johannes, J. Schmalian, and R. Prozorov, “Sign-reversal of the in-plane resistivity anisotropy in hole-doped iron pnictides”, Nature Comm. 4, 1914 (2013).
- (51) P. Goswami, R. Yu, Q. Si, and E. Abrahams, “Spin dynamics of a antiferromagnet and its implications for iron pnictides”, Phys. Rev. B 84, 155108 (2011).
- (52) L. Ma, G. F. Chen, D.-X. Yao, J. Zhang, S. Zhang, T.-L. Xia, and W. Yu, “Na and As NMR study of antiferromagnetism and spin fluctuations in NaFeAs single crystals”, Phys. Rev. B 83, 132501 (2011).
- (53) T. Shimojima, T. Sonobe, W. Malaeb, K. Shinada, A. Chainani, S. Shin, T. Yoshida, S. Ideta, A. Fujimori, H. Kumigashira, K Ono, Y. Nakashima, H. Anzai, M. Arita, A. Ino, H. Namatame, M. Taniguchi, M. Nakajima, S. Uchida, Y. Tomioka, T.Ito, K. Kihou, C. H. Lee, A. Iyo, H. Eisaki, K. Ohgushi, S. Kasahara, T. Terashima, H. Ikeda, T. Shibauchi, Y. Matsuda, and K. Ishizaka, “Pseudogap formation above the superconducting dome in iron-pnictides”, arXiv:1305.3875 (2013).
- (54) W.-C. Lee and P. W. Phillips, “Non-Fermi liquid due to orbital fluctuations in iron pnictide superconductors”, Phys. Rev. B 86, 245113 (2012).
- (55) H. Z. Arham, C. R. Hunt, W. K. Park, J. Gillett, S. D. Das, S. E. Sebastian, Z. J. Xu, J. S. Wen, Z. W. Lin, Q. Li, G. Gu, A. Thaler, S. Ran, S. L. Bud’ko, P. C. Canfield, D. Y. Chung, M. G. Kanatzidis, and L. H. Greene, “Detection of orbital fluctuations above the structural transition temperature in the iron pnictides and chalcogenides”, Phys. Rev. B 85, 214515 (2012).
- (56) Y. Gallais, R. M. Fernandes, I. Paul, L. Chauviere, Y.-X. Yang, M.-A. Measson, M. Cazayous, A. Sacuto, D. Colson, and A. Forget, “Observation of incipient charge nematicity in Ba(FeCo)As”, arXiv:1302.6255
- (57) R. M. Fernandes, A. E. Böhmer, C. Meingast, and J. Schmalian, “Scaling between magnetic and lattice fluctuations in iron-pnictide superconductors”, Phys. Rev. Lett. 111, 137001 (2013).
- (58) E. G. Moon and S. Sachdev, “Competition between superconductivity and nematic order: Anisotropy of superconducting coherence length”, Phys. Rev. B. 85, 184511 (2012).
- (59) J.-Ph. Reid, M. A. Tanatar, A. Juneau-Fecteau, R. T. Gordon, S. Rene de Cotret, N. Doiron-Leyraud, T. Saito, H. Fukazawa, Y. Kohori, K. Kihou, C. H. Lee, A. Iyo, H. Eisaki, R. Prozorov, and L. Taillefer, “Universal heat conduction in the iron-arsenide superconductor KFeAs : Evidence of a d-wave state”, Phys. Rev. Lett. 109, 087001 (2012).
- (60) R. M. Fernandes and A. J. Millis, “Nematicity as a probe of superconducting pairing in iron-based superconductors”, Phys. Rev. Lett. 111, 127001 (2013).
- (61) G. Livanas, A. Aperis, P. Kotetes, and G. Varelogiannis, “Nematicity from mixed states in iron-based superconductors”, arXiv:1208.2881 (2012).
- (62) F. Yang, F. Wang, and D.-H. Lee, “Fermiology, Orbital Order, Orbital Fluctuation and Cooper Pairing in Iron-based Superconductors”, arXiv:1305.0605 (2013).