Arxiv:1707.02457V1 [Cond-Mat.Str-El] 8 Jul 2017

Quantum Order-by-Disorder in Strongly Correlated Metals

Andrew G. Green,1 Gareth Conduit,2 and Frank Krüger1, 3 1London Centre for Nanotechnology, University College London, Gordon St., London, WC1H 0AH, United Kingdom 2TCM, Cavendish Laboratory, JJ Thomson Ave, Cambridge, CB3 0HE 3ISIS Facility, Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire, OX11 0QX, United Kingdom Entropic forces in classical many-body systems, e.g. colloidal suspensions, can lead to the for- mation of new phases. Quantum fluctuations can have similar effects: spin fluctuations drive the superfluidity of Helium-3 and a similar mechanism operating in metals can give rise to supercon- ductivity. It is conventional to discuss the latter in terms of the forces induced by the quantum fluctuations. However, focusing directly upon the free energy provides a useful alternative per- spective in the classical case and can also be applied to study quantum fluctuations. Villain first developed this approach for insulating magnets and coined the term order-by-disorder to describe the observed effect. We discuss the application of this idea to metallic systems, recent progress made in doing so, and the broader prospects for the future.

CONTENTS rubber and the apocryphal attraction between ships in a swell. Such forces arise due a state-dependent restriction I. INTRODUCTION 1 of the spectrum of fluctuations and thus their entropic contribution to the free energy. In many body systems, II. THE FERROMAGNETIC METAL 2 such forces can lead to phases that are favoured for en- A. Simple Model and Historical Development 2 tirely entropic reasons. In colloids, they are responsible B. Fermionic Quantum Order by Disorder 4 for a variety of different phases with characteristic de- 1. Heuristic Description 4 pendence upon the shape of the colloidal particles[1, 2]. 2. Path Integral Formulation 5 Similar effects are responsible for the folding of DNA due C. Phase Reconstruction 6 to the conformational entropy of the surrounding wa- 1. 1st Order transition and Helimagnetism 6 ter, and the celebrated Berezinskii-Kosterlitz-Thouless 2. Superconductivity and Nematic Order 6 transition[3–5] is driven entirely by the entropy of un- 3. Commentary 7 bound vortices. Often, it is revealing to discuss these transitions di- III. EXPERIMENTALLY DRIVEN rectly at the level of the free energy rather than through DEVELOPMENTS 7 the resulting forces. Viewed in this way, the idea of en- A. Ultra-cold Atomic Gases 7 tropically driven order is a powerful unifying concept. B. Metallic Ferromagnets 7 In engineering it comes under the banner of the term C. Experimental Prospects — Beyond the state-dependent noise as exemplified by the noisy in- Ferromagnet 9 verted pendulum[6]. Further afield, in business or evo- lution, fluctuations may make adaptability a favourable IV. RELATION TO OTHER APPROACHES 9 strategy. Essentially, fluctuations can stabilise financial A. QOBD vs DFT, DMFT and Kadanoff-Baym 9 or ecological niches that would not be stable in their ab- B. Functional Renormalisation Group 10 sence. Entropic priors in Bayesian inference bias data C. Monte Carlo 10 analyses in the same way[7]. Quantum fluctuations can generate forces in a very V. BROADER PERSPECTIVE AND similar manner to classical entropic effects, the Van der CONCLUSIONS 11 Waals force being the most famous example[8, 9]. Forces generated by the fluctuations of quantum spins are re- arXiv:1707.02457v1 [cond-mat.str-el] 8 Jul 2017 DISCLOSURE STATEMENT 12 sponsible for the superfluidity of Helium-3[10, 11] and a similar mechanism operating in metals can give rise to ACKNOWLEDGMENTS 12 spin-fluctuation induced superconductivity[12]. In all of these cases, it is conventional to discuss the effect of quan- References 12 tum fluctuations in terms of the forces that they induce. However, focussing directly upon the free energy provides a useful alternative perspective. In the quantum case, it I. INTRODUCTION is the zero-point energy of fluctuations rather than an entropic contribution to the free energy that is at play, The entropic generation of forces is familiar in classical however, the effects are very similar. Indeed, the terms systems. Examples include the elastic forces in stretched Casimir force and Van der Waals force are sometimes 2 generalised beyond their original context of forces due to a heuristic presentation and an overview of a more formal zero-point fluctuations of the electromagnetic vacuum, to field theoretical derivation, emphasising the close parallel refer to entropic forces in colloids. with Villain’s order-by-disorder. As in the classical case, the order-by-disorder approach Villain provided a concrete example of the utility of fo- reveals commonalities between effects that are difficult to cussing directly upon the free-energy[13] and coined the appreciate from other perspectives. This can be partic- term order-by-disorder to describe the observed effect. ularly advantageous when trying to understand experi- He considered a magnetic model whose classical ground- ments. Measurement of the fluctuation spectra and ap- states form a degenerate manifold. Allowing for fluctu- preciation of how it is altered by different types of order ations (either quantum zero-point or thermal[14]) about can lead directly to predictions of what types of instabil- these classical configurations breaks the degeneracy and ity a system is prone to, even in the absence of a detailed picks out a particular ordered state — hence the term microscopic model. We discuss several developments of order-by-disorder. These ideas have been applied in a fermionic quantum order-by-disorder that were directly number of insulating magnets, especially in cases where influenced by experiment. frustration leads to a degenerate manifold of classical The difference between spin-fluctuation theory fo- groundstate configurations that is broken by fluctuations, cussing on forces, and fermionic quantum order-by- though this degeneracy is not necessary for the fluctua- disorder focussing upon the free energy is one of perspec- tion contribution to the free energy to be appreciable. tive. Both encompass the same physics and, in circum- Examples include spin ice pyrochlores[15, 16] and frus- stances where direct comparisons can be made, lead ulti- trated antiferromagnets on the honeycomb lattice[17]. mately to the same equations. However, there are other The main focus in this perspective is the application of methods — especially numerical — that can be used to order-by-disorder to metallic systems. The central phi- analyse strongly correlated quantum systems. We illus- losophy is simple to state: modifications of the Fermi trate the relationship that fermionic quantum order-by- surface, for example by the introduction of some order disorder bears to them; its comparison to various ab initio parameter, modify the electron dispersion and as a result techniques, the potential for future inclusion in density reconfigure the spectrum of low-lying excitations. This functional code, and agreement with Monte Carlo calcu- in turn shifts the zero-point energy (and entropy at fi- lations for the itinerant ferromganet. Finally, we discuss nite temperature) of fluctuations. When fluctuations are the prospects for near- and long-term development of the large, this may self-consistently determine the Fermi sur- technique and how the broader application of the idea has face. The similarity of this to the Casmir effect is clear, resonance with ideas of fluctuation-induced geometry of as is the importance of fermionic statistics. The Fermi entanglement structure and entropic gravity. surface self-consistently provides boundary conditions for the electrons via Pauli exclusion. One could say that II. THE FERROMAGNETIC METAL particle-hole excitations induce generalised Casimir or Van der Waals forces between electrons. We first present these ideas in the context of the critical itinerant ferro- A. Simple Model and Historical Development magnets, where it was initially developed. We give both

The simplest non-trivial fermionic model to which one is approximately quadratic and where screening renders can apply quantum order-by-disorder is that of electrons the effective Coulomb interaction between electrons short 2 with a quadratic dispersion, k = k , and contact inter- ranged — we shall discuss the effects of longer range in- action, g; teractions later. A mean-field analysis allowing for the possibility of a 3 finite magnetization, M = σ nˆx,σ , gives rise to = kcˆk† ,σcˆk,σ + g d x cˆx† , cˆx† , cˆx, cˆx, , (1) σ= H ↑ ↓ ↓ ↑ ± h i k the free energy X,σ Z P

1 β(k σgM µ) 3 2 with cˆx† (k),σ a fermionic creation operator in position MF = ln 1 + e− − − +g d xM (x). (momentum) space. Below, we take the spin label F −β σ = k,σ   Z corresponding to up- and down-spin relative to the lo-± X (2) cal magnetisation. Despite its apparent simplicity, this This is simply the Stoner model of ferromagnetism. It model displays a remarkable range of different phenom- shows a second order transition between paramagnetic ena. Moreover, it finds direct realisation in cold atomic and ferromagnetic phases, the temperature of which gases where atomic interactions are local. It is a good varies with the interaction strength g and which occurs approximation to electrons in solids with a chemical po- at ρFg = 1 at zero temperature (ρF being the density of tential near the bottom of a band so that the dispersion states at the Fermi level). 3

0 Tri-critical Point

-1 Helical FM ) F -2 Ferromagnet Wave Density Pair T/T (

10 -3 log

-4 Ferromagnet & SC

SC -5 0.0985 0.0990 0.0995 0.1000 0.1005 0.1010

1/gvF

FIG. 1. Phase Diagram of the Critical Ferromagnet: Quantum fluctuations in the vicinity of the putative quantum critical point drive a reconstruction of the phase diagram. In the vicinity of the tricritical point, where the transition becomes first order, an expansion of the Ginzburg-Landau function to quartic order is sufficient. Approaching zero temperature, the phase boundary is determined by singularities in the higher order terms of the expansion. Fluctuations drive a superconducting pairing in the p-wave channel by the same mechanism as the pairing in 3He. A novel pair density wave order is found in the region of the phase diagram where superconductivity and helimagnetic order overlap.

In his seminal paper on this model, Hertz[18] realised (a) (b) that fluctuations in certain regions of the phase diagram are profoundly affected by quantum mechanics. This leads to behaviour in a different universality class to that of classical phase transitions; a class that Hertz termed quantum critical. Hertz extended the free en- ergy given in Eq.(2) to an action for(a) the dynamical mag- (b) netization field, the dynamics of which is given by the decay of magnons into particle-hole pairs, i.e. Landau damping[19]. This modification — and the assump- tion that the non-analyticity of the Landau damping does not extend to other terms under renormalization — yields power-law dependence of physical quantities on temperature that are characteristic of the quantum critical point[18, 20, 21]. The presence of strong ferro- magnetic fluctuations near to the zero-temperature limit FIG. 2. Diagrams contributing to non-analytic extensions of the phase transition can also lead to superconducting to Moriya-Hertz-Millis theory: Re-summation of diagrams of pairing in the p-wave channel, via a mechanism trans- these types propagate non-analyticity from Landau damping lated from superfluid[10, 11] 3He to metallic systems by to all diagrams in the expansion. These ultimately lead to Fay and Appel[12]. reconstruction of the phase diagram in the vicinity of the quantum critical point. (a) Contributions to non-analyticity However, this story is not complete. It was first in the order parameter. (b) Contributions to non-analyticity realised by Belitz, Kirkpatrick and Voijta[22] that in gradients of the order parameter. the Moriya-Hertz-Millis theory suffers from an internal inconsistency[23] — the non-analyticity present in Lan- dau damping propagates under renormalization to all 3 terms in an expansion of the action[22, 24–27]. This dimensions and M log M in two dimensions) and its 2 presages the fact that fluctuations in the vicinity of the gradients (of the form q M in three dimensions and q 3/2M 2 in two dimensions).−| | These drive the tran- quantum critical point are in fact so strong that they −| | favour a reconstruction of the phase diagram in the vicin- sition into the ferromagnet first order and favour heli- ity of the quantum critical point. The Ginzburg-Landau magnetic order at low temperatures in the vicinity of the function for the magnetization incurs non-analyticities in quantum critical point. both the magnetization (of the form M 4 log M in three It is natural to ask whether the first order ferromag- 4

(a) nz = Q/Q (b) (c)

nx ny

q

M(r)/M

FIG. 3. Illustration of Helimagnetic Order: In the vicinity of the itinerant ferromagnetic quantum critical point, fluctuations

Monday, 9 May 2011drive the formation of helimagnetic order in which — as illustrated in (a) — the magnetic quantization axis rotates when moving along the direction of the pitch vector q. A realisation of this in the square lattice is shown in (b).

net, helimagnet, and superconducting instabilities have a correction to the free energy is given by common cause. It is not at all obvious from the diagram- matics that they do [See Fig.2]. Moreover, it would be + + 0 f f − (f + f − ) desirable to have a simple heuristic picture that allows 2 k1 k2 k3 k4 (3) fl = 2geff + + , F −  + −  − one to anticipate the helimagnetic instability in advance k1 k2 k3 k4 k1...k4 − − of detailed calculation. As we shall see next, fermionic X quantum order-by-disorder provides this unified descrip- where the summation is taken over momenta such that tion and has a simple heuristic interpretation. k1 + k2 = k3 + k4. The electron dispersion for spin B. Fermionic Quantum Order by Disorder σ is given, for example in the presence of helimag- σ netic order [See Fig.3] with pitch vector q, by k = 1. Heuristic Description 2 2 2 k σ (q k) + g M , corresponding to the disper- − · eff The essence of fermionic quantum order-by-disorder is sion ofq an electron with wave-vector k with spin up or down (σ = + or ) relative to the helimagnetic back- to combine mean-field and fluctuation (zero-point or en- − tropic/thermal) contributions to the free energy whilst ground. A similar replacement is made in the mean field contribution, Eq.(2) and σ is the Fermi distri- self-consistently allowing for the possibility of additional fk order, such as helimagnetic or superconducting order. bution function for occupation of this mode. We have also made a one-loop renormalization of the interaction When the fluctuation contributions are large, they may 2 1 g = g 2g 0 + − + − allowing for the provide the dominant contribution to the free energy and eff k3,k4  +   − k1 k2 − k3 − k4 determine the state adopted by the system. In the case of leading order correctionP to the electron pair wavefunc- the simple model given in Eq.(1), the leading fluctuation tion.

The presence of the background helimagnetic order, with total momentum zero. It turns out that the inte- Mq, modifies the mean-field dispersion of the electrons gral in Eq.(3) is dominated by particle-hole pairs with and hence the fluctuation corrections. These fluctuations momentum near 2kF . The density of such excitations turn the ferromagnetic transition first order[28] and ul- at low energy is enhanced by the presence of ferromag- timately drive helimagnetic order[29] at sufficiently low netic or helimagnetic order (see Fig. 4). Since the sec- temperatures. This can be understood as follows: the ond order perturbative correction is negative, this lowers second order correction to the energy involves the virtual the free energy and drives ferromagnetic or helimagnetic excitation of pairs of particle-hole pairs — one spin-up order in the region of the phase diagram where fluctua- particle-hole pair and one spin-down particle-hole pair — tions contributions to the free energy dominate. Indeed, 5

(a)(a) (b)(b) spin and charge channels using a Hubbard-Stratonovich transformation. At this stage, the Moriya-Hertz-Millis approach (which decouples just the spin part of this inter- action) would be to integrate out the electrons in favour of an effective theory for the Hubbard-Stratonovich spin field. Quantum order-by-disorder, makes a simple, but important modification to this. Anticipating the possibil- (a) (b) (a) (b) ity of static magnetic or charge order, the zero- and finite- frequency parts of the Hubbard-Stratonovich fields are separated. The electrons (whose action is now quadratic) are then integrated out and the resulting effective action for the finite-frequency spin and charge fluctuations in the static background is truncated to quadratic order. Finally, the finite-frequency charge and spin fluctuations FIG. 4. Enhanced low-energy phase space of particle-hole ex- are integrated out to obtain a Ginsburg-Landau function citations: (a) Ferromagnetic or (b) helimagnetic distortions for the spin field[34]. of the Fermi surface enhance the low-energy phase space of By organising the calculation in this way, the com- particle-hole excitations lowering the zero-point energy and parison with Villain’s order-by-disorder[13] is transpar- driving the transition to ferromagnetic order first order or ent. Villain considered magnetic fluctuations about driving helimagnetic order in the vicinity of the quantum crit- some mean-field classical background in an approxima- ical point. tion that treated them as non-interacting. The resulting bosonic Hamiltonian could be diagonalised by a Bogoli- ubov transformation from which the zero-point energy if Eq.(3) is expanded in powers of M and q one finds[30] (and entropic contribution to the free energy) was de- that the leading singular terms M 4 and q 2M 2 both have duced. Here, by separating zero- and finite-frequency coefficients proportional to log T [29].| This| reflects the parts of the Hubbard-Stratonovich fields, the propaga- singularities found in diagrammatic analyses[22, 24–27]. tion of electrons is explicitly calculated in the presence At sufficiently low temperatures near the critical point, of some static, background order, thus accounting (af- these effects win out and drive first order transitions and ter integrating out the electrons) for the modification of helimagnetic order. As temperature is lowered below this spin and charge fluctuations that results. Expansion of point, singularities at higher and higher order in M and the effective action of the finite-frequency fluctuations q dominate[31]. At first glance, Eq.(3) appears to be to quadratic order, similarly mirrors the non-interacting a| term| in an expansion in the interaction strength. Ob- approximation. viously this is not the case as the free energy is a func- Moreover, the following argument demonstrates that tional of the mean-field electron dispersion, which con- there are appreciable zero-point fluctuations: Within the tains the interaction through its dependence upon . M ˆ kF ξ non-interacting approximation, the operator b† = In fact, the expansion is in e− where ξ is a typical q,(σ,p) lengthscale of the interaction potential. As shown by cˆp† +qσcˆpσ with p + q > kF and p < kF can Conduit and Keyserlink in Ref.[32], the fluctuation cor- be considered a| bosonic| creation operator| | for spin- rection Eq.(3) contains a factor in its integrand given by σ particle-hole excitations at momentum q and addi- the square of the Fourier transform of the interaction po- tional internal index p. The mean-field Fermi sur- tential, V (k k ) 2. Since the integrals in Eq.(3) are face provides a vacuum for these excitations. The | 1 − 3 | dominated by k1 k3 2kF , these factors are propor- propagation of these composite bosonic particles is not kF |ξ − | ≈ tional to e− where ξ is a typical lengthscale of the trivial — it is accounted for by precisely the elec- potential. The mean field dispersion does not contain tronic polarization loops that generate Landau damp- this exponential suppression of the interaction since the ing. However, the interaction part of the Hamiltonian potential enters there with near zero momentum. can be expressed simply and contains anomalous terms of

the form g q p0 ,k b† q p bq† k + bq,( ,p)b q,( ,k) , ,( , ) ,( , ) ↑ − ↓ where the prime indicates− ↑ the appropriately↓ restricted 2. Path Integral Formulation P P h i summations over the momenta p and k. It is evident that a Bogoliubov transformation is required to diag- A formal derivation of Eqs.(2,3) from a path integral onalize the bosonic Hamiltonian, suggesting significant was presented in Ref.[33]. Whilst the full details are not zero-point fluctuations and a groundstate of the form[35] important here, aspects are of note because of the con- nection that they reveal to the Moriya-Hertz-Millis the- ory, the approximations used and their limitations. The 0 ψ = exp Uq,k,pb† q,( ,p)bq† ,( ,k) MF | i  − ↑ ↓  | i derivation of proceeds via the following steps: First, a q p,k fermionic path integral is constructed for the partition X X   function, and the contact interaction decoupled in both i.e. dressing the Fermi surface with pairs of spin-up and 6 spin-down particle-hole pairs. the Fermi surface, the pitch of the helimagnetism q en- ters as a directionally dependent magnetization on the Fermi surface. This implies that coefficients of M 2α are 2β 2(α β) C. Phase Reconstruction proportional to those of q M − . When the coeffi- cient of M 4 becomes negative| | at the tricritical point, the coefficient of q 2M 2 also becomes negative indicating a The singular contribution of zero-point fluctuations to minimum of the| | free energy at finite wave-vector and the the free energy of the ferromagnet drives a rich variety formation of magnetisation at that wave-vector. of phases, even in the simple model of Eq.(1). As in- dicated in Fig. 1, these include a first order transition out of the paramagnetic phase at low temperatures, he- 2. Superconductivity and Nematic Order limagnetic order, the coexistence with p-wave supercon- ductivity, and also the possibility of spin-antisymmetric nematic order. Next, we will discuss how all of these can The analysis of how fluctuations drive the ferromag- be accommodated within a fermionic order-by-disorder netic transition first order was helped by the fact that treatment. Including additional, experimentally relevant the interaction had finite weight in the ferromagnetic terms in the Hamiltonian can drive further effects, which channel. Instabilities such as superconductivity or spin we will turn to later. anti-symmetric Pomaranchuk/nematic order do not have finite weight in the bare contact interaction considered in Eq.(1). It is well-understood diagrammatically how in- teractions may be generated in these channels[10]. An 1. 1st Order transition and Helimagnetism appealing picture of how spin fluctuations might drive p-wave superconductivity in a metallic ferromagnet was 1st Order Transition: Fluctuations generate a diverg- given in Ref.[36]: an electron induces a polarization cloud ing negative contribution to the M 4-coefficient, which in the electronic fluid. The net energy of the polar- overcomes the positive mean-field contribution and drives ization around two parallel spins is reduced when the the transition first order at a finite-temperature tri- clouds overlap. In essence this describes a reduction in critical point (See Fig.1), as found diagrammatically[22]. zero-point energy upon the formation of Cooper pairs. As the phase boundary is followed to lower and lower One must think a little to recover the same effects in temperatures, its position is determined by singularities the fermionic quantum order-by-disorder approach[37]. at higher and higher powers of M. In order to track The central idea is to calculate a generating function the transition faithfully down to zero-temperature, these for the new type of order. For example, in the case increasingly singular contributions must be re-summed of p-wave superconductivity in the ferromagnet, a term as in Ref.[31]. Note that the form M 4 log M at zero- k j∆θkck† , c† k, + c.c is added to the Hamiltonian, temperature is non-analytical in M which invalidates a ↑ − ↑ where j is a source field for the superconducting or- conventional Ginzburg-Landau expansion. It does not P  ∆  however invalidate the self-consistent calculation of fluc- der parameter, ∆ = k θkck† , c† k, . The analysis h ↑ − ↑i tuations at finite magnetisation employed in fermionic then proceeds as before, diagonalising the quadratic parts P quantum order-by-disorder. of Hamiltonian and calculating the second order per- Helimagnetism: The fluctuation contribution to the turbation to the free energy. Diagonalisation of the free energy Eq.(3) also contains the seeds of helimag- Hamiltonian in the presence of these sources requires netic order. Heuristically, the ferromagnetic distortion a Bogoliubov transformation, which modifies the inter- increases the low-energy phase space for pairs of particle- action vertex. It is this modification that accounts for the generation of new interaction channels found in the hole excitations with momentum near 2kF and opposite spin. In fact, any spin anti-symmetric distortion will also diagrammatic analysis. Once the fluctuation correc- enhance this low-energy phase space and be favoured by tions to the free energy have been calculated, a Leg- the resulting lowering of zero-point energy (see Fig.4). endre transformation gives the Ginzburg Landau func- tion. A similar approach has been used to study spin- These arguments can be made quantitative by noting antisymmetric Pomeranchuk or nematic distortions of that the free energy is a functional of the mean-field the Fermi surface[33]. electron dispersion, F F [ ]. Aside from a classical k A trick[37] may be used to calculate the Ginzburg- contribution d3xM 2,≡ the dependence upon the order parameter comes entirely from this dependence. The Landau function for continuous transitions: Adding R mean-field electron dispersion in the presence of a heli- and subtracting a term k ∆θkck† , c† k, + c.c to the ↑ − ↑ magnetic order Mq = M(cos q x, sin q x, 0) is given Hamiltonian. Diagonalising the quadratic parts of the · · P   — for a quadratic bare dispersion of electrons — by Hamiltonian alongside the added ∆c†c†-term and treat- σ 2 2  = k σ (k q) + M , describing the propaga- ing the remaining terms alongside the subtracted ∆c c - k,Mq − · † † tion of electrons with momentum k and spin parallel or term perturbatively, recovers the result of adding a source p antiparallel (σ = + or ) to the background helimag- followed by Legendre transformation up to quadratic or- netic order. When the− dispersion can be linearised at der. The saddle-point equations that result from this 7 analysis up to quadratic order are precisely the Eliash- ory. Amongst the first applications of fermionic quantum berg equations found in diagrammatic, spin-fluctuations order-by-disorder was to understand the ferromagnetic calculations. instability of cold fermionic atoms. Early predictions of the scattering length at which this would occur were cru- cially important[38, 39]. Moreover, since atomic gases 3. Commentary are naturally in the canonical ensemble, without a mech- anism to change the total spin of the gas, it is natural to Fermionic quantum order-by-disorder provides a com- consider systems with a fixed, spin imbalance[28, 40]. plementary formulation of the physics embodied by di- Of course, the study of atomic gases is not without its agrammatic analyses. Its alternative perspective uni- experimental difficulties. When tuned to positive scatter- fies the particle-instabilities (such as first order transi- ing lengths, there is an unavoidable loss process[41, 42], tions and the formation of helimagnetism) of the itiner- which affected the first attempt to observe ferromagnetic ant ferromagnet, and particle-hole instabilities such as behavior [43]. Subsequent experiments have focussed on superconductivity. In conventional, diagrammatic anal- the polaron limit [44] for which variational wavefunc- yses, these effects are treated very differently, the former tions may be presented[45]. Going forward the cold atom through non-analyticities of the Moriya-Hertz-Millis free gas offers a chance to explore not only ferromagnetism energy, and the latter through spin-fluctuation generated with the contact interaction, but also systems with in- pairing in the appropriate channel. The discussion above teractions with an effective range [32], a two-dimensional shows how these effects are borne out in a simple model. gas [46], and a mass imbalance between up and down-spin The approach also provides more direct access to new particles [47]. phenomena. Just as in the case of classical, entropically- driven phase transitions, focussing upon the free energy of fluctuations directly can be very useful in deducing B. Metallic Ferromagnets the possibility of other phases and instabilities before a detailed calculation has been performed. Knowing the Unlike ultra-cold atomic gases, for which precise quan- spectrum of fluctuations and how they are affected by titative predictions from realistic microscopic models are the introduction of certain order can be sufficient to pre- often desirable, in the solid-state we often require min- dict the formation of new phases — even in the absence imal models that capture the essence of new emergent of a microscopic model. The formulation has the advan- phenomena. Fermionic quantum order-by-disorder has tage that calculations follow directly from it. In the fol- provided a good guide in a variety of cases. Ferro- lowing section, we outline a number of developments to magnetic transitions are thought to be generically first- fermionic quantum order-by-disorder driven by unusual order at low temperature for precisely the reasons illumi- experimental observation. nated above and this behaviour has been seen in a wide range of materials including Sr1 xCaxRuO3 [48], CoO2 − [49], UGe2[50] and URhGe [51]. Other cases in which III. EXPERIMENTALLY DRIVEN fermionic quantum order-by-disorder has proved illumi- DEVELOPMENTS nating include: Hard Axis Magnetism: In a few itinerant metals, such The experimental systems to which fermionic quan- as YbRh2Si2 and YbNi4P2, hard axis magnetic suscepti- tum order-by-disorder can be applied fall into two main bility in the paramagnetic phase gives way to ferromag- classes, ultra-cold atomic gasses and solid-state itinerant netism along the hard direction in the ordered phase[52– magnets. Experimental capabilities for these two types of 55][56]. Whilst one can fit such behaviour with an insul- system are very different, placing consequently different tating model of frustrated magnetism[57] it is more nat- demands upon theoretical models. ural to understand it from the point of view of quantum order by disorder[58]. Establishing ferromagnetic order in the direction anti-favoured at the mean field level costs A. Ultra-cold Atomic Gases mean field energy, but leads to a flattened dispersion of low energy excitations and so lower fluctuation contribu- Ultra-cold atomic gases offer a pristine arena in which tion to the free energy. When fluctuations dominate, this to study collective quantum phenomena. Interactions determines the direction of magnetic order — i.e. in the can be controlled with great precision by tuning Fesh- direction anti-favoured by mean field considerations (See bach resonances, requiring accurate predictions from the- Fig.5).

Partial order in MnSi: A similar argument can be MnSi[59, 60]. In this material, spin orbit coupling com- used to explain the unusual partially ordered phase of bined with the absence of inversion symmetry conspire 8 2

(a) 2nd order (c) Elastic neutron scattering: (0,0,2+δ) shown samples we have studied (at least 4), with neutron scat- is an intra-cell antiferromagnetic component to the order Q =(0, 0, 2+) tering7 (at D23, ILL and SPINS, NIST) and with resonant (equivalent to that in the FM state), which for clarity is X-ray scatteringparamagnet (BM28, ESRF). Data obtained with res- not shown in the figure. The neutron cross-sections were 1st order onant[K] Temperature X-ray scattering have the highest q-resolution and taken to be averages over opposite tilts and helices. A Easy Axis 6 are shown inspiral FIG 2 1. The incommensurate di↵raction sig- slightly larger b-axis/a-axis moment ratio was found for

T/µ FM nal is visible for X-ray energy at the Pr L2 resonance edge SDW1 (both for q1 and q10 ). Further measurements at (6.4445 keV),spiral which 1 fluorescence and absorbanceq 3 measure-q di↵erent Bragg positions were however not made in this ments suggest is a simple dipole transition. The observed state. intensity Hard Axis intensityhomogeneous in neutron scattering at the same wavevector- FM 4 transfer proveFM that the satellites are of magnetic origin. In both SDW states the modulation vectors are precisely 1/⇢U 3 along the c-axis. (b) 1.8 1.9 2 2.1 2.2

6 250 25 Repeated(a) for 3 samples on 2 diffractometers,(b) 200 20 D235.5 (ILL) and SPINS (NIST) shown. 150 15 Temperature Partial Order 5 T(K) q 100 3q 10 4.5 T ⇤ 50 5

[111] Helical 4 0 0 2.06 2.07 2.08 2.09 2.1 2.2 2.22 2.24 2.26 FIG. 2. Two views of the moment structure for the SDW2 Ferromagnet Non Fermi Liquid L(rlu) L(rlu) state over a length of approximately 1/2 the modulation pe-

6 (c) (d) 3.5 (e) ± q ’ riod, consistent with the measured spin dependence of neu- p⇤ pc 2 Pressure 5.5 q ’/q 3 2 2 tron scattering. Only the Pr Atoms are shown. These can ± q 5 1 ± q ’

1 Ratio be viewed as two bc-plane zig-zag chains of atoms running FIG. 5. Three experimental situations where the fermionic quantumT(K) ± order-by-disorderq approach has proven2.5 fruitful,q ’/q (a) Hard 4.5 2 1 1 along the c-axis with di↵erent a-coordinate (one comprises A axis magnetic order: In materials such as YbRh2Si2 and Yb Ni4P2, hard axis magnetic susceptibility in the paramagnetic phase 4 2 sites and the other B sites). There is an additional intracell gives way to ferromagnetic order along this direction in the ordered0.06 0.07 phase.0.08 0.09 Fermionic0.1 0.2 quantum0.21 0.22 0.23 order-by-disorder0.24 4.5 5 captures5.5 6 this q(rlu) q(rlu) T(K) canting of the moments along the c-axis that is not shown for by an inversion of the magnetic anisotropic in the zero-point fluctuation| contributions| to the| free| energy. The figure indicates the phase diagram calculated for this effect in a simple itinerant model. (b) Partial Order in MnSi: A variety of explanations clarity, but which tilts the moments out of the shaded disks. have been posited for the partial order in MnSi, including the incipient formation of Skyrmion lattices. Fermionic quantum order-by-disorder provides a relatively simple explanation; orienting the helimagnetic order in directions disfavoured in mean field (away from [111]) lowers the energy of fluctuations andFIG. is thus1. (a) favoured and (b) when are fluctuations colour-scale are image strong of near di↵ theracted quantum X- critical point. The phase diagram indicates that this partial ordered phase occurs in the same region of the phase diagram as The above establishes unambiguously that b-axis mo- ray intensity as a function of reciprocal coordinate (0, 0,L) the helimagnetic phase of the simple model for the ferromanget. (c) Helimagnetic Order in PrPtPl: This material is one of the ments are present in the modulated states.The presence best examples so far for a fluctuation-induced helimagneticand phase. temperature The figure for summarises scattering neutron at the Pr(top) L2 and(L x-ray1) resonance (bottom) edge (6.444 keV). (c) and (d) are the modulation q-vectors of the third order harmonic for SDW2 provides indirect scattering data obtained from PrPtAl. This supports the notion that helimagnetic order appears in region of the phase diagram evidence that the SDW2 state is in fact a spiral state. where quantum fluctuations are strong. The theory mustas be a modified function to of incorporate temperature the effects deduced of local from spins. (a) & After (b). doingq1 and so, a range of thermodynamic and scattering measurements areq10 accountedare the wavectors for by the associatedmodel. with SDW1. q2 and q20 are For a spiral along z with a constant magnitude of mo- associated with SDW2. denote peak positions relative to ment m , the magnetisation in the presence of crystalline ± 0 (0,0,2). (e) The ratio of modulation vectors. q20 /q2 =3shows anisotropy can be written to favour helimagnetic order, with a very long pitch [61–that q20foris the it in third ZrZn harmonic2[64] and of NbFeq2. There2[65], isthis no has simple not ratio been 63]. Crystalline anisotropy further favours the pitch vec-betweenconclusive.q1 and q10 . Indeed, it was noted in Ref.[66] that the cos[q z (z)] m = m m . (1) tors lining up along the [111] directions. In the partially, fluctuation-induced helimagnet might be rather suscep- 0 sin[q z (z)] ordered phase near to the quantum critical point, neu- tible to the effects of disorder, possibly explaining the ✓ m ◆ tron scattering shows a smooth distribution of orienta- apparent spin glass behaviour in the region of the phase tions of the pitch vector over the sphere that if any-To furtherdiagram ofexplore CeFePO(Ref.[67]) the nature where of magnetic helimagnetism order might we The inclusion of the phase (z) allows the pitch of the thing anti-favours the [111] directions with a slightmeasured in- otherwise the spin-dependencebe predicted. More of recently, the neutron two materials cross- spiral to adjust to reduce the anisotropy energy. This crease in the pitch of the helimagnetic order. This be-sectionhave at wave-vector emerged as compelling transfer (0 candidates, 0, 2 q) with for the the emergence SPINS energy is minimised for (z)= sin(2q z)withq the ± m m haviour arises naturally with a minimal modificationinstrument of of fluctuation-induced at NIST. This revealed helimagnetic the presence order, PrPtAl[68] of spin- principal modulation vector. A significant third order the ferromagnetic Hamiltonian, Eq.(1), that allows forflip (SF)and scattering LaCrGe3[69, of 70]. intensity To capture 1/10 the the properties non spin-flip of PrP- amplitude (/2 for small ) is predicted. The intensity of the Dzyalosinksii-Moriya spin orbit coupling and fluctua-(NSF)tAl, scattering the simple for modelincoming of Eq.(1) neutron must moments be supplemented aligned the 3rd order harmonic can be calculated from the mea- tions (See Fig.5). Directing the helimagnetism away from with local spins. A range of thermodynamic and scatter- along the a-axis for SDW2. Significant SF scattering in sured ratio of SF/NSF (SF/NSF 0.1 for (0,0,2.07)) the mean-field favoured directions flattens the dispersion ing measurements are accounted for by the predictions of ⇡ of excitations, lowering their zero-point energy. this geometryquantum shows order-by-disorder that the modulated (Fig.5). state has a mag- scattering from EQN 1. This gives a third harmonic in- Fluctuation-induced helimagnetism: The experimentalnetic moments along the b-axis as well as along a.This tensity 7% of the fundamental, which agrees well with search for fluctuation-induced helimagnetism has provenwould be expected for a spiral. The examination of that that observed. more difficult. Whilst there is circumstantial evidenceSDW2 satellite close to other Bragg positions revealed The changes of magnetic structure are also seen in ther- there is also a moment component along the c-axis (par- modynamic measurements. The heat capacity has cusp- allel q) [see suppl info]. like features at all three transitions (FIG 2). The scatter A magnetic structure consistent with these measure- in values at the FM-SDW2 and SDW1-SDW2 transitions ments for the SDW2 state is shown in FIG 2. In the (depending on the temperature history) indicate that figure the moments rotate in a plane whose normal is in- these transitions are almost certainly first order. The dc clined from c-axis (about the a-axis). In addition there susceptibility, namely magnetisation divided by the field HANNAPPEL, PEDDER, KRUGER,¨ AND GREEN PHYSICAL REVIEW B 93,235105(2016) HANNAPPEL, PEDDER, KRUGER,¨ AND GREEN PHYSICAL REVIEW B 93,235105(2016)

(a) q =0 (a) q =0 A+ A+

q/η =0 q/η =0.q/3 η =0 q/η =1q/η.0=0.3 q/η =1.0

λx = λ↑ λ↓ A x x λx = λx↑ λx↓ A − − − − 9

λy = λy↑ λy↓ λy = λ↑ λ↓ − y − y (a) (b)(b)q = q(1, 1)/√2 (b) q = q(1, 1)/√2 ω , ω , ↓ ± ↑ ↓ ↑ ± ↑ ↑

FIG. 5. Fermi surfaces A of the electronic bands ϵ (k) in FIG. 5. Fermi surfaces A of the electronic bands± ϵ (k) in ± the presence of helical spin-tripletthe± presenced-density of helical wave spin-triplet order with± d-densityq wave order with q q(1,1,0)/√2. Red and blue colors denote the spin-up and spin-down= q(1,1,0)/√2. Red and blue colors denote the spin-up and spin-down= character of the lobes. As we increase the value of q,movingfrom character of the lobes. As we increase the value of q,movingfrom left to right, we see that this spin character gets mixed, and the Fermi q left to right, we see that this spinsurfaces character deform gets along mixed, the (1 and,1,0) the direction. Fermi surfaces deform along the (1,1,0) direction. spin x-direction spin y-direction B. Visualization in momentum space and real spacespin x-direction spin y-direction B. Visualization in momentum space and real space FIG. 6. Visualization of the spin-triplet nematic order parameter The helical spin-triplet d-density wave is not easy to FIG. 6. Visualizationon of a squarethe spin-triplet lattice. (a) nematic The homogeneous order parameter state corresponds to bond The helical spin-tripletvisualize.d-density For wavesmall q isvectors, not easy corresponding to to a long period on a square lattice. (a) Theorder homogeneous which breaks the state rotation corresponds symmetry to bondof the square lattice. The FIG. 6. visualize.Twisted For Spin-Antisymmetric small q vectors,of the corresponding modulation Nematic in Order: to real a long space,Translating period a Wigner representation the mechanism as of fermionic quantum order-by-disorder to order which breaks theorder rotation parameter symmetry changes of the sign square under lattice. 90◦ rotation, The as well as under the spin-antisymmetric nematicused reveals in Fig. an unusual1 is the most intertwining convenient of depiction. helimagnetic This is and a spin-antisymmetric nematic order whereby of the modulation in real space, a Wigner representation as order parameter changesspin sign inversion, under 90 and isrotation, invariant as under well the as twounder combined operations. the spin quantization axis of themixed spin-antisymmetric real/momentum space nematic representation. order Over rotates a subsystem in space. a) illustrates◦ this order in the continuum used in Fig. 1 is the most convenient depiction. This is a spin inversion, and is invariant(b) Bond-density under the wave two order combined corresponding operations. to the helical spin-triplet showing howmixed electrons real/momentum with opposite spacewhose representation. spinsize is have less then Over d-wave the a subsystem wavelength distortions of the of modulation, their fermi-surfaces the d-density with wave opposite with q along sign. the The (1,1) quantization direction. The two panels of this spin axis rotates upon movingorder is inapproximately the direction uniform of the and onepitch may(b) vector Bond-density defineq a. quasi- b) wave illustrates order corresponding the projection to the helical of this spin-triplet order onto a whose size is less then the wavelength of the modulation, the d-density wave with qshowalong the thex and (1,y1)spin direction. components The of two the modulated panels order parameter, square lattice, demonstrating itsFermi-surface correspondence equivalent to a to bond that of order the related modulation. homogeneous respectively. order is approximately uniform and one may define a quasi- show the x and y spin components of the modulated order parameter, Fermi-surface equivalent tospin-triplet that of the nematic. related The homogeneous helical modulation in spin space implies that the spin direction rotates fromrespectively. subsystem to spin-triplet nematic. The helical modulation in spin space ν subsystem with a period 2π/q. in terms of expectation values of bond operators, λx(y) C.implies Experimental that the spin Prospects direction rotates — Beyond from subsystem the to tum order-by-disorder, that in fact fluctuations promote = Apurelymomentumspacepictureisalsousefulasit ψr†,ν ψr xˆ ( ˆy),ν .Theorderparameter˜νη is shown in Fig. 6(a). subsystem with a period 2π/q. in terms of expectation⟨ values+ of⟩ bond operators, λx(y) Ferromagnethelps reveal how spatial modulation mightthe beformation favored by of aIt new changes phase sign near under to spin the inversion, quantum as= wellcritical as under 90◦ Apurelymomentumspacepictureisalsousefulasit ψr†,ν ψr xˆ ( ˆy),ν .Theorderparameter˜η isˆ s shown in Fig. 6(a). the softening of fluctuations. In Fig. 5,weshowtheFermipoint that intertwinesrotation. magnetic Because modulationR1(r) 0, the and straind-wave components of ⟨ + ⟩ ⟨ ⟩ = helps reveal how spatialsurfaces modulationA and mightA befor favored the twoby electronicIt changes bands ϵ sign(k) underspin-up spin inversion, and spin-down as well electrons as under exactly 90◦ cancel each other, The physics underlying fermionic quantum+ − order-by- orbital order± in= a continuumˆ s version of bond density wave the softening of fluctuations.k2 In√( Fig.k q)52,weshowtheFermiη2$2(k), with the waverotation. vector q Becausein the R(λ1↑(r) λ↑)0,(λ the↓ strainλ↓) 0. components of 1 order(see Fig.6).⟨ Thesex ⟩ =y unusualx y phases have some in- disordersurfaces is veryA generaland A andfor isfavored the expected∓ two (1 electronic,·1,0) to+direction. find bandsA applica-ϵand(kA) are thespin-up Fermi surfaces and spin-down for In electrons− the helical+ exactly spin-triplet− cancel= d-density each other, wave, the spin direction tion far beyond2 √ the+ 2 itinerant2− 2 ferromagnet to which+ it± has−= triguing observable consequences. For example, whereas k (k q) η $1(k),electrons with the with wave spin parallel vector andq in antiparallel the ( toλx↑ the backgroundλy↑) (λx↓ λy↓rotates) 0. in a plane in spin space, e.g., between the x and y initiallyfavored been∓ applied. (1,·1,0)+ direction. For example,Ahelimagneticand A theare ordering, thespin-triplet Fermi respectively. surfaces ne- for In thethe limitIn− static theq helical+0, spin-triplet we spin-triplet− directions,= d nematic-density as specified wave, responds by the the spin order todirection parameter a uniformη(r)(9). This = matic stateelectronsfirst with proposed spin parallel in [71]recover+ and and antiparallel the− studied elliptical to the in Fermi abackground simple surfaces ofmagneticrotates the homogeneous in a field plane by in generatingcan spin again space, be expressedan e.g., anisotropic between in terms the of strain[81],x expectationand y values the of bond model byhelimagnetic a number ordering, of authors[72–82] respectively.spin-triplet nematic.was In the first limit These shownq deformations0, to we oftripletdirections, the Fermi d-density surfaces as specified waveoperators, by the generates order parameter a spatiallyη(r)(9). modulated This recover the elliptical Fermichange surfaces the spectrum of the of homogeneous electronic= particle-holestrain.can excitations again This be expressed and offers in new terms possibilities of expectation for values experimentally of bond have its transition driven firstenhance order the by phase fluctuations space for fluctuations. in spin-triplet nematic. These deformations of the Fermi surfaces operators, η˜x (r) %r†σx %r xˆ %r†σx %r ˆy η˜ cos(qr), Ref.[83]. It was later shown [84], using fermionic quan- isolating multipolar order.= ⟨ + ⟩−⟨ + ⟩ = change the spectrum of electronicWe conclude particle-hole this section excitations by providing and a real-space picture of the homogeneous spin-triplet nematic and the modulated η˜y (r) %r†σy %r xˆ %r†σy %r ˆy η˜ sin(qr). enhance the phase space for fluctuations. = ⟨ + ⟩−⟨ + ⟩ = triplet d-density wave states when projected ontoη˜x a(r lattice.) % Thisr†σx %r xˆ %r†σx %r ˆy η˜ cos(qr), = ⟨ + ⟩−⟨The order-parameter+ ⟩ = componentsη ˜x (r)and˜ηy (r)areshown We conclude this sectionillustrates by providing the connection a real-space of our picture continuum model to lattice- η˜y (r) %r†σy %r inxˆ Fig. 6(b)%r†σfory % arq vectorˆy η˜ alongsin(qr (1,)1). that is commensurate with of the homogeneous spin-tripletbased models nematic of bondand the density modulated wave order. For simplicity,= ⟨ we + ⟩−⟨ + ⟩ = triplet d-density wave states when projected onto a lattice. This the underlying square lattice. Figure 6 is in essence a lattice consider a two-dimensional square lattice. WeThe discretize order-parameter the projection components of theη ˜ Wignerx (r)and˜ representationηy (r)areshown shown in Fig. 1. The interplayillustrates the of connection fluctuations of our with continuum band model effectsˆ t tois lattice- of1 ment.2 2 order parameter η(r) R1(r) 2 %†(rin)σ Fig.(∂x 6(b)∂y )%for(r a) q, vector along (1,1) that is commensurate with particularbased importance models of bond in real-life densitywhich wavesystems. (for order.fixed α The For=1) simplicity, band⟨ is a three-dimensional⟩ dis-= we⟨ vector− in spin⟩ the underlying square lattice. FigureV. EXPERIMENTAL6 is in essence SIGNATURES a lattice OF persion itselfconsider may a two-dimensional lead to phasespace. square reconstruction For lattice. the homogeneous We= discretize in some spin-triplet the nematic state along t projection of the Wigner representationTHE SPIN-TRIPLETshown in Fig. 1 NEMATIC. order parameter η(r) Rˆthe(rz)spin1 direction,% (r)σ we(∂2 obtain∂2) the%(r lattice) , order parameter ways similar to that caused by1 fluctuations2 † drivingx y tran- IV. RELATION TO OTHER APPROACHES which (for fixed α =1)⟨ is a three-dimensional⟩ = ⟨ vector− in spin⟩ Spin-triplet nematic order simultaneously breaks spatial sitions first order when= the Wolfharth-Rhodes criterion V. EXPERIMENTAL SIGNATURES OF space. For the homogeneous spin-triplet nematicη˜ (λx↑ stateλy↑) along(λx↓ λy↓) rotation symmetry and spin-rotation symmetry. This entan- is satisfied[85] or leading to spin-antisymmetric= nematic− − − THE SPIN-TRIPLET NEMATIC the z spin direction, we obtain the lattice order parameter As we have emphasised, fermionic quantum order-by- or spin density wave due to weak nesting as found in disorderSpin-triplet provides235105-8 nematic a ordercomplementary simultaneously view breaks on spatial the role of Sr RuO [86, 87]. A proper account of this interplay has 3 3 η˜ (λx↑ λy↑) (λx↓ λy↓) fluctuationsrotation symmetry in forming and spin-rotation exotic quantum symmetry. phases. This entan- We have not yet been given, though= − a preliminary− − analysis was emphasised the relation to spin-fluctuation theories and presented in [88]. The extreme of such interplay occurs in235105-8 — in the application to the itinerant ferromagnet — non- the critical antiferromagnet in which spin fluctuations are analyitic extensions of Moriya-Hertz-Millis theory. There thought to drive d-wave superconductivity[89–91], and are a variety of other approaches to which it is related in in which the subtle interplay of band effects and spin spirit, and it is worth mentioning some of them here. fluctuations still leaves open questions. Analysis using fermionic quantum order-by-disorder may be useful[92]. A range of other systems of experimental interest may also be analysed using fermionic quantum order-by- A. QOBD vs DFT, DMFT and Kadanoff-Baym disorder. For example, the lattice Dirac fermions found in graphene harbour a number of instabilities that might be Whilst analytically tractable models play an important captured in this way[93–95]. Similarly, the gradual addi- role in the qualitative understanding of new phenomena, tion of new features, such as spin-orbit coupling, topolog- sometimes numerically precise predictions are needed, ical bandstructures etc., to the toy Hamiltonian, Eq.(1) particularly when guiding experiment through a delicate might reveal new and unusual features. The possibilities balance of competing possibilities. Though progress has are very rich and best explored in concert with experi- been made towards first principle calculations[96, 97], 10 presently, the calculation of fluctuation-induced effects, tronic excitations and order-parameter fluctuations. One such as superconductivity, is carried out in parallel with of the main advantages of the fRG-based one-loop com- ab initio band structure calculations; the latter provid- putation of the two-particle vertex, compared to other ing spin susceptibilities that are used in field theoretical weak-coupling approaches, is that particle-particle and calculations. It would be appealing if these calculations particle-hole channels are treated on an equal footing, could be performed in concert. QOBD suggests a nat- without artificial bias towards a particular channel of in- ural way to do so. It sits naturally in the framework stability [102]. described in Ref.[98], which places density functional the- Practical implementations require approximations. ory, dynamical mean field and Baym-Kadanoff theory in Since the flow of a given vertex depends upon higher- the unified context of Legendre transformation of an ap- order vertices, the effective action must be truncated at propriate generating functional. The potential for sys- some order. Moreover, an infrared cut-off Λ that serves tematically including these effects into ab intio codes has as RG flow parameter needs to be introduced in the great promise. free-fermion action. Typical choices are a momentum Fermionic quantum order-by-disorder is essentially a cut-off that defines a small region around the Fermi sur- restriction of Baym-Kadanoff theory where the varia- face, similar to the Wilsonian RG approach to interact- tional parameter — the Green’s function — is itself char- ing fermions [104–107], or a frequency cut-off. Once self- acterised by a restricted set of parameters given by the energy effects are taken into account, the Fermi surface order parameters of the fields that we wish to study. is usually deformed in the course of the RG flow. Using a Indeed, the additional fluctuation corrections resulting frequency cut-off has the advantage that it does not inter- from the inclusion of superconducting order (see Ref.[37] fere with Fermi-surface deformations and that particle- where these results are obtained as outlined in Section III processes with small momentum transfer are captured C), are essentially the modfications to the density func- smoothly by the flow [108]. Finally, numerical integra- tional used for ab initio calculations of electronically me- tion requires a discretisation of the fRG equations — diated superconductivity in Refs.[96, 97]. This presents usually achieved by introducing a mesh in momentum broader possibilities for the inclusion of novel order, such space. Because of this, subtle changes associated with as spin-antisymmetric nematic order, in density func- the onset of an exotic order parameter could be missed. tional code. If successful, this would complement the For example, the fluctuation-driven helical magnetic or- analytical description of minimal models with numerical der close to ferromagnetic quantum-critical points has a accuracy when applied to real materials, and would be a very large periodicity, corresponding to very small Fermi- potentially fruitful avenue of future study. surface deformations. Because of the computational complexity, fRG has so far only been successfully applied to one- and two- B. Functional Renormalisation Group dimensional electron systems. In this context, it is inter- esting to mention that early fRG studies of graphene-type Fluctuations-driven ordering phenomena are naturally systems predicted a topological quantum spin-Hall state captured in renormalisation-group (RG) approaches. As for strong next-nearest neighbor interactions [109]. More discussed previously, certain instabilities are not directly recent fRG studies with a much higher momentum reso- supported by short range Coulomb repulsion, but fluc- lution, however, show that this transition is pre-empted tuations renormalise the low-energy theory and generate by the formation of a three-sublattice, charge-modulated (at one-loop order) new interaction vertices that support state [110]. It would be interesting to determine the pre- dictions that QOBD makes for this system. While the the instability. The p-wave superconductivity in metal- lic ferromagnets[10] is a good example. A more recent fRG equations in their full form are exact, the commonly example are the spin-fluctuation driven nematic and su- used truncation schemes restrict their applications to the perconducting instabilities of iron-based superconductors study of weak-coupling instabilities. At present it re- [99, 100]. That the same physics can be described by the mains an open question which truncation schemes of the QOBD mechanism [37] shows the close relation between fRG equations are required to study the full range of the two approaches. fluctuation-driven ordering phenomena captured by the QOBD mechanism. The functional renormalization group (fRG) is, in prin- ciple, an exact method that can be applied to both bosonic [101] and fermionic [102] systems. A comprehen- sive summary can be found in the review by Salmhofer C. Monte Carlo [103]. Instead of the scale dependence of a finite number of coupling constants, fRG keeps track of the full mo- The exquisite precision of experiments on ultra-cold mentum and frequency dependence of the single-particle gases alluded to in Section IV. A, places exacting de- Green function and interaction vertices. It is a flexi- mands upon theoretical predictions. Precise determina- ble and unbiased tool to study scale dependent behav- tion of scattering lengths (the preferred measure of in- ior in electron systems, e.g. competing magnetic, charge, teraction strength) at which phase transitions and insta- and pairing instabilities and the interplay between elec- bilities occur must be met with similarly precise calcu- 11 lations. Various quantum Monte-Carlo calculations have restricted only by the fixed node approximation. The been performed to achieve this — in broad agreement Quantum Monte Carlo approach also offers access to two with QOBD. important extensions: the consequences of effective range Quantum Monte Carlo [111] provides an accurate nu- interactions [32], two-dimensional geometry [46], and dif- merical method to probe the effects of quantum fluc- ferent masses for the spin-up and spin-down particles [47]. tuations in strongly correlated systems. The method However, the numerical calculations are performed in a not only offers a complementary analysis underpinned finite sized system and so cannot capture all long wave- by different approximations to the QOBD formalism, length fluctuations. but also a precise level of control over the Hamilto- nian and in particular the interaction potential. The starting point for fermionic systems is a Slater-Jastrow V. BROADER PERSPECTIVE AND J trial wavefunction [111], ψ = e D D . Here D , CONCLUSIONS is a Slater determinant of the up-/down-spin↑ ↓ particle{↑ ↓} plane-wave states corresponding to the non-interacting Fermionic QOBD has had a number of successes — ground state, and ensures the correct fermion antisym- particularly in its application to ferromagnetic metals. metry. Different polarizations are addressed by changing Development of the method to improve its range and ac- the number of up-/down-spin particles in their respec- curacy of application suggests several near- and medium- tive determinants. J is a Jastrow factor that contains term goals. additional variational terms in particle-particle separa- Near-term goals might include considering the forma- tion [112], J = α (r r )n. In Variational n,i,j n,σi,σj i − j tion of phases near the quantum critical points of differ- Monte Carlo the parameters αn,σiσj are optimized to ent background order — we discussed above extensions P { } minimize the ground state energy. Further refinement of the original application to ferromagnetic order to spin- of this wavefunction is can be obtained using Diffusion anti-symmetric nematic order and anti-ferromagnetic or- Monte Carlo; an accurate Green’s function method that der. Charge density wave order is another interesting projects out the ground state [111]. possibility. However, the possibilities for QOBD effects Most ultra-cold gases have contact interactions with to operate in bands with non-trivial topological order is zero effective range. These are characterised by the scat- perhaps more intriguing. Implicit in the latter is the ef- tering phase shift δ, which can be related to the scatter- fect of more intricate band structure and spin-orbit cou- ing length a 0 according to cot δ(k) = 1/(ka) (for pling. Conventional band structure can itself generate a zero-range≥ interaction). Several different− approaches instabilities towards different types of order through nest- have been used to capture such interactions in Monte- ing and other features in the density of states. A system- Carlo: i. An explicitly repulsive top hat potential [29]. atic way to treat these alongside fluctuations does not yet The strength and radius of the top hat can be adjusted to exist. Addition of more ingredients in this way makes as- give zero effective range at the cost of introducing higher sessment of which effects win out more and more delicate order terms in the scattering phase shift [46]. ii. The — a systematic, numerical way to accommodate this (for contact interaction limit can be reached by using an at- example by inclusion in density functional theory) along- tractive interaction potential with short range [113, 114]. side analytical calculations for minimal toy models is de- However, this potential harbors a bound state. To avoid sirable. the system entering the bound state the trial wave func- Beyond these considerations, the inclusion of new tion must be restricted to have no variational parameters, types of instabilities and exotic order presents a num- meaning that the ground state cannot be reached. iii. ber of medium-term goals. The description of local crit- The short range behavior of the wave function induced ical fluctuations[117] from this perspective and indeed by the contact interaction can be directly imprinted upon whether Mott physics can be accessed by QOBD raise the Jastrow factor [113] However, this gives incorrect important questions. There is some hope for the latter scattering properties in the ` 1 scattering channels. given the similarities between QOBD and the description Notwithstanding these subtleties,≥ it has proven possible of the Mott physics using dynamical mean field theory. to extract a general, transferable pseudopotential that Indeed a rather direct translation of the simple Landau gives the correct scattering properties, and which does theory developed in Ref.[118] may indeed be able to re- not harbor a bound state, allowing it to be used in a express the Mott transition in terms of QOBD. variationally [115]. Unconventional order such as that found in gapless The Quantum Monte Carlo approaches have been used spin-liquid phases presents even greater challenges[119]. to confirm the presence of a first order transition and he- The fermionic version of QOBD might possibly be ap- limagnetic phase [29] Moreover, these results have also plied to understand instabilities of the spinon-fermi sur- been used to define a density functional[116], thus ful- face. In its wider application there have been recent hints filling part of the aim of bringing QOBD effects into the that QOBD might be useful to understand deconfined realm of density functional theory. These results provide quantum criticality in terms of fluctuation-induced tran- independent verification of the QOBD approach: they in- sitions in entanglement structure[120–122]. The 1/N- clude all possible correlations, and deliver a ground state effects of fluctuations in anti-ferromagnets[123] can be 12 reframed as a zero-point energy breaking the degener- way to proceed. In some cases, QOBD allows a partic- acy amongst a set of states with different entanglement ularly simple, heuristic picture of the physical processes structure[122]. This brings us full circle to the work of promoting certain collective behaviour, one that can very Villain[13] and shows that his perspective can be used to directly be turned to concrete calculation. It is an idea understand a large range of phenomena. This broader whose ramifications have yet to be fully realised. perspective is closely aligned to one of the most fasci- nating developments in fundamental physics; the notion that gravity itself might be an entropic force[124–126] emerging from zero-point fluctuations. The idea being DISCLOSURE STATEMENT that the space-time metric can be associated with entan- glement structure, and that this entanglement structure The authors are not aware of any affiliations, member- is itself determined by quantum fluctuations. In the case ships, funding, or financial holdings that might be per- of fundamental physics, these ideas arise from the du- ceived as affecting the objectivity of this review. ality between gauge theory and string theory. Villain’s simple notion of order-by-disorder applied to condensed matter systems, may yet provide an alternative route to emergent geometry of entanglement. ACKNOWLEDGMENTS Quantum order-by-disorder is one amongst many per- spectives from which to view strongly correlated quan- We acknowledge funding from the EPSRC and Royal tum systems. Each perspective privileges a particular Society

[1] Onsager L. 1949. Annals of the New York Academy of [20] Moriya T. 2012. Spin fluctuations in itinerant electron Sciences 51:627–659 magnetism. vol. 56. Springer Science & Business Media [2] Adams M, Dogic Z, Keller SL, Fraden S. 1998. Nature [21] Millis A. 1993. Physical Review B 48:7183 393:349–352 [22] Belitz D, Kirkpatrick T, Vojta T. 1997. Physical Review [3] Berezinskiˇı VL. 1971. Soviet Journal of Experimental B 55:9452 and Theoretical Physics 32:493 [23] In fact, the existence of these non-analyticities had been [4] Berezinskiˇı VL. 1972. Soviet Journal of Experimental noted earlier[127] but their physical significance was not and Theoretical Physics 34:610 understood until the work of Belitz-Kirkpatrick-Vojta. [5] Kosterlitz JM, Thouless DJ. 1973. Journal of Physics [24] Betouras J, Efremov D, Chubukov A. 2005. Physical C: Solid State Physics 6:1181 Review B 72:115112 [6] Simons YB, Meerson B. 2009. Phys. Rev. E 80:042102 [25] Rech J, Pepin C, Chubukov AV. 2006. Physical Review [7] MacKay DJ. 2003. Information theory, inference and B 74:195126 learning algorithms. Cambridge university press [26] Efremov DV, Betouras JJ, Chubukov A. 2008. Physical [8] London F. 1930. Zeitschrift fur Physik 63:245–279 Review B 77:220401 [9] Casimir HB. 1948. On the attraction between two per- [27] Maslov DL, Chubukov AV. 2009. Physical Review B fectly conducting plates. In Proceedings of the KNAW, 79:075112 vol. 51 [28] Conduit G, Simons B. 2009. Physical Review A [10] Anderson PW, Brinkman W. 1973. Physical Review Let- 79:053606 ters 30:1108 [29] Conduit G, Green A, Simons B. 2009. Physical review [11] Brinkman W, Serene J, Anderson P. 1974. Physical Re- letters 103:207201 view A 10:2386 [30] The first calculation of the integrals in Eq.(3) was car- [12] Fay D, Appel J. 1980. Physical Review B 22:3173 ried out in Refs.[128–130] at zero magnetic field. To [13] Villain J, Bidaux R, Carton JP, Conte R. 1980. Journal date, no closed form for the integral at finite M has de Physique 41:1263–1272 been found. Initial studies[28, 38] calculated the integral [14] The very similar effects of zero-point fluctuations and numerically and demonstrated how fluctuations could thermal fluctuations is revealed particularly clearly; a drive the transition in to the ferromagnet first order. factor of nB + 1/2 for each spinwave mode accounts for Later, leading singularities where calculated analytically both zero-point (as T 0) and thermal (as T ) and re-summed[31]. occupation. → → ∞ [31] Pedder C, Krüger F, Green A. 2013. Physical Review B [15] Tchernyshyov O, Moessner R, Sondhi SL. 2002. Phys. 88:165109 Rev. Lett. 88:067203 [32] von Keyserlingk C, Conduit G. 2013. Physical Review [16] Zhitomirsky ME, Gvozdikova MV, Holdsworth PCW, B 87:184424 Moessner R. 2012. Phys. Rev. Lett. 109:077204 [33] Karahasanovic U, Krüger F, Green AG. 2012. Physical [17] Mulder A, Ganesh R, Capriotti L, Paramekanti A. 2010. Review B 85:165111 Phys. Rev. B 81:214419 [34] And potentially the charge field, although in applica- [18] Hertz JA. 1976. Physical Review B 14:1165 tions to date, it has been assumed that there is no static [19] Landau LD. 1946. Zh. Eksp. Teor. Fiz. 10:25 spatial modulation of the charge field, which can then be absorbed into the chemical potential 13

[35] Compare, for example to the dressing of the Neél state [65] Crook M, Cywinski R. 1995. Journal of magnetism and by spinwave excitations: magnetic materials 140:71–72 " # [66] Thomson S, Krüger F, Green A. 2013. Physical Review X ˆ ˆ B 87:224203 ψ = exp (uk/vk)bk† b† k ...... , | i − | ↑↓↑↓ i [67] Lausberg" S. 2012. Physical Review Letters 109 k [68] Abdul-Jabbar G, Sokolov DA, O’neill CD, Stock C, where uk and vk are the coefficients used in the Bogoli- Wermeille D, et al. 2015. Nature Physics 11:321–327 ubov diagonalisation of the spinwave Hamiltonian. [69] Lin X, Taufour V, Bud’ko SL, Canfield PC. 2013. Phys- [36] Julian S. 2012. Physics 5:17 ical Review B 88:094405 [37] Conduit G, Pedder C, Green A. 2013. Physical Review [70] Taufour V, Kaluarachchi US, Khasanov R, Nguyen B 87:121112 MC, Guguchia Z, et al. 2016. Physical Review Letters [38] Duine R, MacDonald A. 2005. Physical review letters 117:037207 95:230403 [71] Kivelson SA, Fradkin E, Emery VJ. 1998. Nature [39] Conduit G, Simons B. 2009. Physical review letters 393:550–553 103:200403 [72] Pomeranchuk IY, et al. 1958. Sov. Phys. JETP 8:361 [40] Parish MM, Marchetti FM, Lamacraft A, Simons BD. [73] Oganesyan V, Kivelson SA, Fradkin E. 2001. Physical 2007. Nat Phys 3:124–128 Review B 64:195109 [41] Conduit G, Altman E. 2011. Physical Review A [74] Quintanilla J, Schofield A. 2006. Physical Review B 83:043618 74:115126 [42] Pekker D, Babadi M, Sensarma R, Zinner N, Pollet L, [75] Kee HY, Kim EH, Chung CH. 2003. Physical Review B et al. 2011. Physical review letters 106:050402 68:245109 [43] Jo GB, Lee YR, Choi JH, Christensen CA, Kim TH, [76] Khavkine I, Chung CH, Oganesyan V, Kee HY. 2004. et al. 2009. Science 325:1521–1524 Physical Review B 70:155110 [44] Kohstall C, Zaccanti M, Jag M, Trenkwalder A, Massig- [77] Kee HY, Kim YB. 2005. Physical Review B 71:184402 nan P, et al. 2012. Nature 485:615–618 [78] Yamase H, Oganesyan V, Metzner W. 2005. Physical [45] Bugnion P, Conduit G. 2013. Physical Review A Review B 72:035114 87:060502 [79] Doh H, Friedman N, Kee HY. 2006. Physical Review B [46] Conduit G. 2013. Physical Review B 87:184414 73:125117 [47] von Keyserlingk C, Conduit G. 2011. Physical Review [80] Lawler MJ, Barci DG, Fernández V, Fradkin E, Oxman A 83:053625 L. 2006. Physical Review B 73:085101 [48] Uemura Y, Goko T, Gat-Malureanu I, Carlo J, Russo [81] Wu C, Sun K, Fradkin E, Zhang SC. 2007. Physical P, et al. 2007. Nature Physics 3:29–35 Review B 75:115103 [49] Otero-Leal" M. 2008. Physical Review B 78 [82] Lawler MJ, Fradkin E. 2007. Physical Review B [50] Taufour" V. 2010. Physical Review Letters 105 75:033304 [51] Yelland EA, Barraclough JM, Wang W, Kamenev KV, [83] Kirkpatrick T, Belitz D. 2011. Physical review letters Huxley AD. 2011. Nat Phys 7:890–894 106:105701 [52] Lausberg S, Hannaske A, Steppke A, Steinke L, Gruner [84] Hannappel G, Pedder C, Krüger F, Green A. 2016. T, et al. 2013. Phys. Rev. Lett. 110:256402 Physical Review B 93:235105 [53] Steppke A, Küchler R, Lausberg S, Lengyel E, Steinke [85] Wohlfarth E, Rhodes P. 1962. Philosophical Magazine L, et al. 2013. Science 339:933–936 7:1817–1824 [54] Ishida K, Okamoto K, Kawasaki Y, Kitaoka Y, [86] Grigera S, Gegenwart P, Borzi R, Weickert F, Schofield Trovarelli O, et al. 2002. Physical review letters A, et al. 2004. Science 306:1154–1157 89:107202 [87] Berridge A, Green A, Grigera S, Simons B. 2009. Phys- [55] Klingner C, Krellner C, Brando M, Geibel C, Steglich ical review letters 102:136404 F, et al. 2011. Physical Review B 83:144405 [88] Conduit G. 2009. Collective phenomena in correlated [56] The susceptibilities along the different primary axes semiconductors, degenerate fermi gases, and ferro- cross just above the magnetic transition temperature, electrics. Ph.D. thesis, University of Cambridge, De- the hard axis becoming an easy axis. partment of Physics [57] Andrade EC, Brando M, Geibel C, Vojta M. 2014. Phys- [89] Fradkin E. 2013. Field theories of condensed matter ical Review B 90:075138 physics. Cambridge University Press [58] Krüger F, Pedder C, Green A. 2014. Physical review [90] Zaanen J, Chakravarty S, Senthil T, Anderson P, Lee letters 113:147001 P, et al. 2006. Nature Physics 2:138–143 [59] Pfleiderer C, Reznik D, Pintschovius L, Löhneysen Hv, [91] Nagaosa N. 1999. Quantum field theory in strongly Garst M, Rosch A. 2004. Nature 427:227–231 correlated electronic systems. Springer Science & [60] Krüger F, Karahasanovic U, Green AG. 2012. Physical Business Media review letters 108:067003 [92] Hannappel G. 2016. The role of fluctuations near an- [61] Ishikawa eY, Tajima K, Bloch D, Roth M. 1976. Solid tiferromagnetic and spin-triplet nematic quantum crit- State Communications 19:525–528 ical points. Ph.D. thesis, University College London, [62] Bak P, Jensen MH. 1980. Journal of Physics C: Solid http://discovery.ucl.ac.uk/1489651/ State Physics 13:L881 [93] Lemonik Y, Aleiner I, Fal’ko V. 2012. Physical review b [63] Pfleiderer C, Julian SR, Lonzarich GG. 2001. Nature 85:245451 414:427–430 [94] Dahal HP, Wehling TO, Bedell KS, Zhu JX, Balatsky [64] Uhlarz M, Pfleiderer C, Hayden S. 2004. Physical review A. 2010. Physica B: Condensed Matter 405:2241–2244 letters 93:256404 14

[95] Scherer MM, Uebelacker S, Honerkamp C. 2012. Phys- [129] Lee T, Yang C. 1957. Physical Review 105:1119 ical Review B 85:235408 [130] Abrikosov A, Khalatnikov I. 1958. Sov. Phys. JETP [96] Essenberger F, Sanna A, Linscheid A, Tandetzky F, 6:888 Profeta G, et al. 2014. Physical Review B 90:214504 [97] Essenberger F, Sanna A, Buczek P, Ernst A, Sandratskii L, Gross E. 2016. Physical Review B 94:014503 [98] Kotliar G, Savrasov SY, Haule K, Oudovenko VS, Par- collet O, Marianetti C. 2006. Reviews of Modern Physics 78:865 [99] Fang C, Yao H, Tsai WF, Hu J, Kivelson SA. 2008. Phys. Rev. B 77:224509 [100] Fernandes RM, Chubukov AV, J. Knolle IE, Schmalian J. 2012. Phys. Rev. B 85:024534 [101] Wetterich C. 1993. Phys. Lett. B 301:90 [102] Metzner W, Salmhofer M, Honerkamp C, Meden V, Schönhammer K. 2012. Rev. Mod. Phys. 84:299 [103] Salmhofer M. 1999. Renormalization, an introduction. Berlin, Heidelberg: Springer-Verlag [104] Benfatto G, Gallavotti G. 1990. Phys. Rev. B 42:9967 [105] Shankar R. 1991. Physica A (Amsterdam) 177:530 [106] Shankar R. 1994. Rev. Mod. Phys. 66:129 [107] Polchinski J. 1992. Nucl. Phys., B 422:617 [108] Husemann C, Salmhofer M. 2009. Phys. Rev. B 79:195125 [109] Raghu S, Qi XL, Honerkamp C, Zhang SC. 2008. Phys. Rev. Lett. 100:156401 [110] de la Peña DS, Lichtenstein J, Honerkamp C. 2016. Competing electronic instabilities of extended hubbard models on the honeycomb lattice: A functional renor- malization group calculation with high wavevector res- olution. ArXiv:1606.01124 [111] Needs R, Towler M, Drummond N, Ríos PL. 2009. Jour- nal of Physics: Condensed Matter 22:023201 [112] Whitehead T, Michael M, Conduit G. 2016. Physical Review B 94:035157 [113] Pilati S, Bertaina G, Giorgini S, Troyer M. 2010. Phys- ical review letters 105:030405 [114] Chang SY, Randeria M, Trivedi N. 2011. Proceedings of the National Academy of Sciences 108:51–54 [115] Bugnion P, Ríos PL, Needs R, Conduit G. 2014. Physical Review A 90:033626 [116] Ma PN, Pilati S, Troyer M, Dai X. 2012. Nature Physics 8:601–605 [117] Si Q, Rabello S, Ingersent K, Smith JL. 2001. Nature 413:804–808 [118] Kotliar G, Lange E, Rozenberg M. 2000. Physical review letters 84:5180 [119] Lee PA, Nagaosa N, Wen XG. 2006. Reviews of modern physics 78:17 [120] Senthil T, Vishwanath A, Balents L, Sachdev S, Fisher M. 2004. Science 303:1490–1494 [121] Kuklov A, Matsumoto M, Prokof’Ev N, Svistunov B, Troyer M. 2008. Physical Review Letters 101:050405 [122] Green A, Hooley C, Keeling J, Simon S. 2016. arXiv preprint arXiv:1607.01778 [123] Read N, Sachdev S. 1989. Physical Review Letters 62:1694 [124] Maldacena J. 1998. Int. J. Theor. Phys 38:1113 [125] Maldacena J. 2005. Scientific American 293:56–63 [126] Verlinde E. 2011. Journal of High Energy Physics 2011:1–27 [127] Geldart DJW, Rasolt M. 1977. Phys. Rev. B 15:1523– 1532 [128] Huang K, Yang CN. 1957. Physical review 105:767