Criticality and Phase Diagram of Quantum Long-Range O(N) models

# Criticality and Phase Diagram of Quantum Long-Range O(n) models

Nicolò Defenu Institut für Theoretische Physik, Universität Heidelberg, D-69120 Heidelberg, Germany    Andrea Trombettoni CNR-IOM DEMOCRITOS Simulation Center, Via Bonomea 265, I-34136 Trieste, Italy SISSA, Via Bonomea 265, I-34136 Trieste, Italy INFN, Sezione di Trieste, I-34151 Trieste, Italy    Stefano Ruffo SISSA, Via Bonomea 265, I-34136 Trieste, Italy INFN, Sezione di Trieste, I-34151 Trieste, Italy
July 3, 2019
###### Abstract

Several recent experiments in atomic, molecular and optical systems motivated a huge interest in the study of quantum long-range systems. Our goal in this paper is to present a general description of their critical behavior and phases, devising a treatment valid in dimensions, with an exponent for the power-law decay of the couplings in the presence of an symmetry. By introducing a convenient ansatz for the effective action, we determine the phase diagram for the -component quantum rotor model with long-range interactions, with corresponding to the Ising model. The phase diagram in the plane shows a non trivial dependence on . As a consequence of the fact that the model is quantum, the correlation functions are anisotropic in the spatial and time coordinates for smaller than a critical value and in this region the isotropy is not restored even at criticality. Results for the correlation length exponent , the dynamical critical exponent and a comparison with numerical findings for them are presented.

###### pacs:
Valid PACS appear here

## I Introduction

A series of remarkable advancements in the last decade lead to a striking development of the experimental techniques for the control and manipulation of atomic, molecular and optical systems such as Rydberg atoms Saffman et al. (2009), dipolar quantum gases Lahaye et al. (2009), polar molecules Ritsch et al. (2013), multimode cavities Carr et al. (2009) and trapped ions Blatt and Roos (2012). Such progress paved the way for the simulation of a variety of quantum models Bloch et al. (2008), and in this direction one of the most active field of research is the implementation and study of equilibrium and dynamical properties of quantum long-range (LR) models Britton et al. (2012); Schauß et al. (2012); Aikawa et al. (2012); Lu et al. (2012); Yan et al. (2013); Islam et al. (2013); Richerme et al. (2014); Jurcevic et al. (2014); Douglas et al. (2015); Schempp et al. (2015); Landig et al. (2015, 2016). Ising and/or quantum spin chains with tunable LR interactions can currently be realized using ions stored in a Penning trap Britton et al. (2012), neutral atoms coupled to photonic modes of a cavity Douglas et al. (2015); Landig et al. (2015, 2016) or with trapped ions coupled to motional degrees of freedom Islam et al. (2013); Richerme et al. (2014); Jurcevic et al. (2014). A key property in these studies is that the resulting interactions decay algebraically with the distance and that the decay exponent can be experimentally tuned Islam et al. (2013); Richerme et al. (2014); Jurcevic et al. (2014). As an example, Rydberg gases have been used to observe and study spatially ordered structures Schauß et al. (2012) and correlated transport Schempp et al. (2015). Dipolar spin-exchange interactions with lattice-confined polar molecules were observed as well Yan et al. (2013). The possibility of control LR interactions in spin chains simulated with trapped ions was also at the basis of the recent experimental simulation of the Schwinger model Martinez et al. (2016).

An important reason for the interest in the equilibrium and non-equilibrium properties of quantum LR systems is the connection with typical themes of classical LR physics. The traditional interest in the behavior of LR interacting statistical mechanics models has been largely due by the outnumbering possible applications in condensed matter, plasma physics, astrophysics and cosmology Codello et al. (2015a); Campa et al. (2014). Therefore, the general question is how quantum fluctuations modify the traditional picture of LR interactions into complex systems. We focus on the study of the criticality in quantum LR systems and the development of a general renormalization group (RG) approach for their study.

From the point of view of critical behavior, the main difference between classical and quantum systems is due to the presence of unitary quantum dynamics in the latter case. The critical behavior in the time domain is usually described in terms of the dynamical critical exponent , which is related to the critical scaling of the gap in the thermodynamic limit , where is the linear size of the system. This is different from classical systems, where the value of depends on the chosen dynamics.

In several short-range (SR) systems the value of this exponent is strictly one, leading to the well known equivalence between the universality class of a dimensional quantum system and its classical equivalent in dimensions, where is the dimension of the corresponding classical system. For non-unitary values of the dynamical critical exponent , which is the case of LR interacting quantum systems, the relation between classical and quantum critical behavior is more subtle. The general expectation is that classical and quantum universalities should be connected in general for Sachdev (2011).

The aim of our investigation is focused on the derivation of the critical exponent and of its spatial counterpart, i.e. the correlation length critical exponent . Various relations, exact or not, which connect both classical and quantum, SR and LR systems, are also investigated. In order to accomplish this task we employ the functional renormalization group (FRG) formalism already developed for classical isotropic and anisotropic LR models Defenu et al. (2015, 2016).

The interplay between LR interactions and quantum effects is a traditional topic in condensed matter physics. A major example is provided by the Ising model Thouless (1969), displaying a behavior related to the spin- Kondo problem with the occurrence of a topological phase transition of the Berezinskii-Kosterlitz-Thouless (BKT) type Yuval and Anderson (1970); Cardy (1981). Moreover, the experimental realization of dipolar Ising spin glasses in the ’s led to investigations of the behavior of quantum LR Ising and rotor models Dutta and Bhattacharjee (2001).

The application of functional RG techniques to classical LR spin systems produced a nice all-in-one picture of the phase transitions occurring in these models, contributing to clarify the behavior of the anomalous dimension of these systems Defenu et al. (2015); Sak (1973); Picco (2012); Blanchard et al. (2013); Brezin et al. (2014). Moreover functional RG formalism was successfully applied to the case of anisotropic LR interactions Defenu et al. (2016). Good agreement with Monte Carlo (MC) simulations was found, not only in the two dimensional Ising model Luijten and Blöte (2002); Angelini et al. (2014); Mori (2010); Horita et al. (2016), but also in -dimensional LR bond percolation by means of effective dimension relations Gori et al. (2017).

Comparing with traditional perturbative RG analysis, functional RG was able to produce, at least for SR models, numerical values for the critical exponents, whose accuracy remains stable as a function of all relevant system parameters, i. e. the dimension and number of components . In the classical SR models the accuracy never falls below 20%, while remaining well below 10% for all continuos symmetries Codello and D’Odorico (2013); Codello et al. (2015b). Such accuracy is expected to hold also in the LR case, where correlation effects should never be stronger than in the two dimensional classical SR Ising model.

The issue we address is the general description of the criticality in quantum LR spin models and the development of a RG framework for this investigation. As a first, one would naively think that the problem is reduced to the study of a LR classical problem, but in this paper it is shown explicitly that this is the case only for large enough. A main aim of this paper is therefore to widen and unify the theoretical picture for the phase transitions quantum LR models and to construct the full landscape of their phase diagram within a single formalism. We also determine the dynamical critical exponent , which characterizes the critical quantum dynamics.

The structure of the paper is the following. In Section II we introduce the models used in the rest of the paper. In Section III we discuss the field theoretical description and the RG methods, and we present our ansatz for the effective action. We also derive the mean field results and study the corrections to the phase diagram due to the effects of the presence of anomalous dimension. Estimates for the critical exponents in and are presented in Section IV, together with a discussion of the comparison with some numerical results available in literature. Our conclusions are presented in Section VI, together with a summary of the relevance and implications of our results.

## Ii Long-Range models

In this paper we consider two models, namely the Ising model in a transverse field and the quantum rotor model. The Ising Hamiltonian in a transverse field reads

 HI=−∑ijJij2σziσzj−H∑iσxi, (1)

where the and are the Pauli matrices defined in the sites of a -dimensional hyper-cubic lattice. The quantum rotor model can be written as

 HR=−∑ijJij2^ni⋅^nj+λ2∑iL2i, (2)

where the are components unit length vector operators (), is a real constant and is the invariant operator formed from the asymmetric rotor space angular momentum tensor Sachdev (2011). The interaction matrix in both cases is power-law decaying with the distance

 Jij=Jrd+σij (3)

with a positive constant, the spatial dimension of the system and any real number, while is the distance between the sites and .

The interaction potential in Eq.(3) induces in the model radically different behaviors depending on the value of the exponent . In particular, for the internal energy of the system may diverge in the thermodynamic limit, leading to an ill defined model. However, a suitable redefinition of the interaction strength Uhlenbeck et al. (1963) produces finite interaction energies, still preserving many interesting results typical of non additive systems Campa et al. (2014).

For fast enough decaying interactions the system is always additive and thermodynamics is well defined. For lying into a certain range, the finite temperature system may present a phase transition with spontaneous symmetry breaking at a finite critical temperature .

For a classical spin system the phase diagram in the plane is therefore divided into three regions Sak (1973). For the universal behavior is the one obtained at mean field level, while for larger than a critical value the system has the same critical behavior of its SR analogue. Finally for LR interactions are relevant and the system has a phase transition with peculiar LR critical behavior.

The determination of the boundary where LR interactions become irrelevant respect to the usual SR ones has been at the center of intense investigations in last few years Brezin et al. (2014); Angelini et al. (2014); Defenu et al. (2015), see as well the recent works Behan et al. (2017a, b). The functional RG framework Tetradis and Wetterich (1992); Morris (1994); Berges et al. (2002) proved to be useful to explore the effect of LR interactions Defenu et al. (2015, 2016) in a compact way. The accuracy of such results remains reliable for all values of the system parameters, i.e. the spatial dimension , the number of components and the criticality index .

Rotor models are the straightforward generalization to the quantum case of the classical models Sachdev (2011). Their low energy behavior describes the physics of many relevant physical models, in particular the case is in the same universality of antiferromagnetic quantum Heisenberg spin systems, while the case is related to the quantum critical behavior of the Bose-Hubbard model.

In this paper we only consider the case in order to have a well defined thermodynamics Codello et al. (2015a); Campa et al. (2014). The system undergoes a paramagnetic to ferromagnetic quantum phase transition (QPT) at zero temperatures for all and in the Ising case, equation (1), i.e. the one component case of the rotor model in equation (2). On the other hand, for , due to the Mermin-Wagner theorem (MWT) the system manifests spontaneous symmetry breaking (SSB) only for in , where is the threshold value above which short-range (SR) behavior is recovered.

At finite temperature the system has a finite temperature classical phase transition (CPT), which lies in the same universality class of standard LR classical spin systems. The results for the finite temperature transition behavior were subject of intense, long-lasting investigations in the last decades Sak (1973); Defenu et al. (2015) and they will be used, but not further discussed, in this paper.

## Iii Field theoretical representation and phase diagram

The universal behavior of condensed matter and statistical mechanics models can be investigated by means of field theoretical techniques Mussardo (2010). In order to employ the functional RG approach, one may consider the Wetterich equation Wetterich (1993)

 ∂tΓk=12Tr[∂tRkΓ(2)+Rk], (4)

where is an infrared regulator function, a finite scale proportional to the inverse system size, , and the logarithmic scale. is the effective action of the system and plays the role of an exact Ginzburg-Landau free energy, while is its second derivative with respect to the system order parameters.

In order to solve equation (4) it is useful to restrict ourselves to a functional space spanned by a finite number of functions and/or couplings. Using the field theory representation and the Trotter decomposition, it can be shown that the models defined in equations (1) and (2) display the same universal behavior of a classical systems in dimension where the interaction is LR in directions and SR in the remaining Trotter dimension Dutta and Bhattacharjee (2001).

Thus a convenient ansatz for the effective action of an quantum rotor model is

 Γk=∫dτ∫ddx{ Kk∂τφi∂τφi−ZkφiΔσ2φi − Z2,kφiΔφi+Uk(ρ)} (5)

where is the spatial Laplacian in dimensions, is the ”Trotter”/imaginary time direction, is the -th component () of the system magnetization density, is the system order parameter, is the scale dependent effective potential and (,,) are three scale dependent wave-function renormalization terms. In equation (III) the summation over repeated indexes is intended.

The effective action ansatz (III) is sufficient to investigate, at least at approximate level, the low energy behavior of the models described in Hamiltonians (2) and (1). The presence of two kinetic terms in the spatial directions are necessary to take into account for the competition between the LR non analytic momentum term in the propagator and the usual term in the region.

The frequency and momentum dependence of the propagator at criticality are connected by the dynamical critical exponent , defined by . We also expect the momentum dependence of the propagator for large wavelength to obey the scaling form

 limq→0G(q)∝q2−η. (6)

Equation (6) defines the anomalous dimension as the deviation of the long wavelength behavior of the propagator from the standard SR mean field behavior . Finally, as for CPT, the correlation length of the quantum system diverges close to the critical point. Such divergence is power-law, at least for standard SSB, and the scaling behavior reads

 ξ∝(λ−λc)−ν (7)

where is the correlation length, is the coupling appearing in Hamiltonian (2) and is the critical value of the coupling at which the system undergoes SSB. Then equation (7) can be considered the definition of the critical exponent . The universal behavior of any second order QCP can be described only in terms of the three exponents with all the other exponents given by scaling relations Sachdev (2011).

### iii.1 Local potential approximation

The lowest order approximation for the ansatz (III) is obtained considering only the actions parametrized by the effective potential . Imposing

 Zk=Kk=1, Z2,k=0, (8)

it is possible to derive the flow equation in the local potential approximation (LPA). The result is equivalent to the one obtained in the case of classical anisotropic spin system Defenu et al. (2016). One gets

 ∂tUk(ρ) =(d+z)~Uk(~ρ)−(d+z−σ)~ρ~U(1)(~ρ) −11+~U(1)(~ρ)+2~ρ~U(2)(~ρ)−N−11+~U(1)(~ρ). (9)

The superscript indicates rescaled quantities, which are defined in order to ensure scale invariance of the effective action at the critical point

 U(ρ) =kDU~U(~ρ) (10) φ =kDφ~φ, (11)

with

 DU =d+z (12) Dφ =d+z−σ. (13)

At LPA level the kinetic term is fixed and does not renormalize, thus anomalous dimension effects vanish in both the frequency and momentum sectors. Then both the anomalous dimension and the dynamical critical exponent attain their mean field values

 η =2−σ, (14) z =σ2. (15)

Also equation (III.1) is equivalent to the flow for the effective potential in the LPA approximation for a classical SR spin system Codello and D’Odorico (2013); Codello et al. (2015b) in an effective fractional dimension

 d′′SR=2(d+z)σ, (16)

with the latter equivalence does not account for anomalous dimension effects. Substituting the mean field result for the exponent one obtains

 d′SR=2dσ+1, (17)

which is exact in the vanishing anomalous dimension cases Joyce (1966). Similar effective dimension relations have already been introduced in literature for the diluted model Dutta (2003), where the random disorder is introduced to simulate LR effects. The critical exponents of the LR model can be obtained from the ones of the SR model in dimension via some simple linear transformations and they fulfill all the usual scaling relations included hyper-scaling, as long as . We refer to Defenu et al. (2015) for a discussion on the reliability of the effective dimension, see as well Behan et al. (2017b).

We are then in position to predict the upper critical dimension, which is obtained from the condition , in such a way that

 duc=32σ, (18)

in agreement with the one obtained by relevance arguments in Dutta and Bhattacharjee (2001). For the system undergoes SSB with mean field critical exponents given in equations (14) and (15), with

 ν =σ−1. (19)

For continuous symmetries in the SR interacting case the anomalous dimension vanishes also at the lower critical dimension as ensured by the MWT. Thus one can again employ the effective dimension approach to obtain

 dlc=σ2. (20)

It is worth noting that, while relation (18) is valid for all , included the Ising case, the expression (20) is limited to continuos symmetries .

Summarizing, the final picture emerging from LPA approximation is rather simple. Equation (III.1) has, in general, two fixed point solutions, one is the Gaussian fixed point, which represents the mean field universality, characterized by the critical exponents in equations (14), (15) and (19). The other solution has finite renormalized mass and represents the interacting Wilson-Fisher (WF) universality.

For the Gaussian universality is the attractive one and the quantum LR system will display SSB with mean field exponents. On the other hand for the attractive fixed point is the WF one and no analytic expression for the critical exponents is known. Finally for no SSB is possible for continuous symmetries and the system will have a single phase in the case.

At the border region the LR and SR momentum terms are competing in the propagator and anomalous dimension effects become increasingly important Codello et al. (2015b). In this section we discarded anomalous dimension effects and the threshold values at which LR interactions become irrelevant with respect to SR ones is as can be detected by mean field arguments Dutta and Bhattacharjee (2001). However, it seems by now established that the actual boundary between LR and SR critical behavior is given by Sak’s expression Sak (1973), where is the anomalous dimension of the system in the case. In order to investigate this effect we shall rely on the complete ansatz (III).

### iii.2 Anomalous dimension effects

In order to introduce anomalous dimension effects it is necessary to consider the flow of the kinetic sectors. In the simplified approach considered in this paper it is sufficient to compute the running of the coefficients , and . The flow for the LR wave-function is given by

 ∂tZk=(2−σ−η)Zk, (21)

where the anomalous dimension is defined with respect to the SR term

 η=∂tZ2,kZ2,k. (22)

At the fixed point the running of the couplings should vanish and thus we have only two possibilities, either

 η=2−σ (23)

or

 limk→0Zk=0. (24)

The latter result will imply that LR interactions are irrelevant and do not influence the critical behavior. Latter analysis is consistent with Sak’s result. Indeed for LR interactions are relevant and , while for the LR term vanishes at criticality and the system recovers SR universal behavior .

Therefore the boundary value is obtained when the coefficient in the r.h.s of equation (21) vanishes and attains its SR value :

 σ∗=2−ηSR. (25)

For the effective action of the system at the fixed point becomes isotropic and the dynamical critical exponent is . Then, according to quantum to classical correspondance Mussardo (2010); Sachdev (2011) the QCP of the SR system is in the same universality of the CPT of its classical analogous in .

The phase diagram of a quantum LR rotor model in the plane is reported in figure 1. In the discrete symmetry case, figure 1(a), the system has a mean field validity region, indicated by the cyan shaded area, where the exponents are the ones obtained by mean field approximation. Otherwise for the system has peculiar non mean field exponents. The boundary after which the system recovers SR behavior is indicated by the red solid line, the line has been computed using expression (25) with the expression for obtained by ansatz (III) without any non analytic terms. For discrete symmetries SSB appears also in .

For continuous symmetries, figure 1(b), the scenario is the same as for discrete symmetries, except that, due to the MWT, the boundary recovers its mean field value not only at but also for at . For , gray shaded area, the system displays only a single phase.

## Iv The critical exponents

As shown in the previous section the system has a phase diagram similar to the one already depicted for a classical LR spin systems Angelini et al. (2014); Defenu et al. (2015) with non trivial LR behavior occurring only for and . In this region no analytical expressions are possible for the critical exponents and one shall resort to numerical approximated techniques.

In order to compute the critical exponents in the non trivial region one can apply the same procedure as done in Defenu et al. (2016) for anisotropic LR spin systems. Again we consider the ansatz (III) in the case of vanishing analytic corrections in the momentum sector. Indeed the latter term is relevant only close to the boundary and it has been shown in Defenu et al. (2015, 2016) to produce only minor corrections even for

The flow of the effective potential for vanishing analytic correction reads

 ∂t¯Uk=(d+z)¯Uk(¯ρ)−(d+z−σ)¯ρ¯U′k(¯ρ)−σ2(N−1)1−ητz3σ+2d1+¯U′k(¯ρ)−σ21−ητz3σ+2d1+¯U′k(¯ρ)+2¯ρ¯U′′k(¯ρ), (26)

where the anomalous dimension has been introduced. This quantity introduces a correction in the frequency dependence of the system propagator

 limω→0G(ω,1)−1∝ω2−ητ. (27)

The expression for the frequency anomalous dimension is readily obtained by the flow of

 ητ=f(~ρ0,~U(2)(~ρ0))(3σ+2d)d+(3σ+d)(1+f(~ρ0,~U(2)(~ρ0))), (28)

where the function is the expression for the spatial anomalous dimension of the correspondent SR range model

 f(~ρ0,~U(2)(~ρ0))=4~ρ0~U(2)(~ρ0)2(1+2~ρ0~U(2)(~ρ0))2 (29)

as is found in Codello and D’Odorico (2013) after rescaling an unessential geometric coefficient.

Merging the definition of the dynamical critical exponent with the definition of in equation (27) it is readily obtained that

 z=σ2−ητ, (30)

which has to be compared with the mean field expression (15).

Eq. (26) allows for an higher order effective dimension relations which connects the universal behavior of a dimensional quantum LR system to its classical SR counterpart in dimensions with

 dSR=(2−ηSR)(d+z)σ, (31)

where is the anomalous dimension of the SR system in dimension . Such relation can be easily demonstrated using the procedure depicted in Defenu et al. (2015) to relate Eq. (26) to the potential flow shown in Codello et al. (2015b). However the exact demonstration only holds if we neglect the coefficients proportional to appearing in the second line of equation (26). These coefficients are small regulator dependent quantities which produce contributions to the universal quantities. It has been verified that neglecting such coefficients produces numerical errors below 1% in the solution of flow equations like (26) Berges et al. (2002); Defenu et al. (2015).

In figure 2 the dynamical critical exponent as a function of the decay exponent is shown in the , panel (a) and (b) respectively. The data have been obtained solving the expression for the fixed point effective potential, equation (26) with the l.h.s posed to zero, and equation (28) in a self consistent cycle.

The dynamical critical exponent attains its mean field value for , while it becomes a non trivial curve in the non Gaussian region. It should be noted that correlation effects always increase the value of the dynamical critical exponent with respect to the mean field prediction. However in the continuos symmetry cases , for , due to the MWT the curves bend down at some finite value and the dynamical critical exponent recovers its SR value only at . On the other hand in the case the Ising model curve (red solid line) is monotonically increasing and met the SR line (gray solid line in figure 2) at the value .

The dynamical critical exponent values obtained by the present approach for the Ising model, red solid line in 1(a), are in very good agreement with MC simulation on an Ising spin system coupled to an anomalous Bosonic bath, which lies in the same universality, red dots in 1(a). The agreement is poorer in the case, blue solid line theoretical values and blue dots MC simulation. The MC data are taken from reference Sperstad et al. (2012).

The lower accuracy found for the can be due to the presence of the BKT mechanism for this model. In passing we note that the MC points seems to provide a value with even for the , which is consistent with the BKT scenario.

The presented computation, despite having the merit to be valid for any , and , is anyway not able to fully capture the effects of the topological phase transition and no qualitative difference is found between the and the for in . However, close to for it appears that MC simulations may be as well plagued by finite size effects as it will become evident discussing the results of figure 3, since they do not reproduce the expected limiting behavior for the correlation length exponent, see also the discussion in Sperstad et al. (2012). On the other hand our calculation of the quantity , see figure 3, is fully consistent with the expected exact behavior coming from expansions Br??zin and Zinn-Justin (1976) and the discrepancy found between the theoretical prediction and the MC data for in the case is perhaps due to finite size errors of the MC data rather than to unexpected BKT effects for .

In order to completely characterize the critical behavior of quantum LR models it is necessary to derive the values for the correlation length critical exponent in the non trivial region. In order to accomplish this task it is sufficient to make a linear perturbation of the fixed point potential with the form , where the denotes the any scaling solution of equation (26).

Employing the technique already introduced in Codello et al. (2015b); Morris and Turner (1998) it is possible to obtain the discrete spectrum of the eigenvalues , which describe the scaling behavior of any physical perturbation of the fixed point effective action. As it is well known for standard second order phase transitions only one relevant perturbation of the critical Ginzburg-Landau free energy exist, which either drives the system to its high or low temperature phases. Then the eigenvalues spectrum will contain only one positive solution, namely the inverse correlation exponent .

Another equivalent way to obtain the correlation length exponent is to take advantage of the effective dimension relation (31) in order to map the values of in the classical SR case to the quantum LR case via the relation

 ν=2−ηSRσνSR, (32)

where is the anomalous dimension of the SR system in dimension . This procedure is equivalent to solving the stability spectrum of equation (26) once the regulator dependent coefficients in the non-linear part have been neglected Defenu et al. (2015). In figure 3 the results for the exponent are reported in the case of a quantum LR rotor model in dimension and panel 3(a) and 3(b) respectively. The results for the Ising criticality (red solid line) are in good agreement with the numerical findings (red dots).

On the other hand the case [blue solid line in figure 3(a)] is, once again, in disagreement with numerical findings (blue dots). In this case, it appears that MC simulations are not completely reliable since they do not reproduce the exact behavior of the correlation length exponent in the limit, where it is expected to diverge, see also the discussion in Sperstad et al. (2012). Such divergence is in agreement with the BKT behavior, where exponential divergence of the correlation length is found, as it can be derived exactly generalizing the expansion technique of the SR non linear -model in reference Br??zin and Zinn-Justin (1976) to the present case.

Latter generalization is more readily obtained by means of effective dimension approach. One should consider equation (21) of reference Br??zin and Zinn-Justin (1976) at lowest order in

 ν−1SR=(d−2)+O(~ε2) (33)

employing the relation (31) with the effective dimension given by equation (17), in the limit of vanishing anomalous dimension. The result is

 ν≃1−σ2 (34)

which when multiplied by gives

 limσ→2(zν)−1=2σ−1,∀N≥2. (35)

It should be noted that the derivation above does not consider anomalous dimension effects which we expect to be sub-leading with respect to the main contribution to . The latter result (black dot dashed line in 3(a)) is in perfect agreement with our findings, but it cannot be reproduced by MC simulations possibly due to finite size effects close to the region.

While we expect numerical simulations to be very accurate only for well below , it is also found that our results for should be slightly underestimated in the intermediate region, even if reproducing all the limiting behaviors. This is in agreement with previous investigation on classical LR and SR O(N) models Codello and D’Odorico (2013); Codello et al. (2015b); Defenu et al. (2015), where the FRG results for the correlation length exponent is found to slightly overestimate high precision numerical estimation in the whole dimension range. Here the same effect is present as a function of the exponent . Thus it is fair to conclude that the exact correlation length exponent curve for the case in the case should be located in between the theoretical curve and the reported numerical estimates. Apart from the case, which is rather special exhibiting the BKT transition, we expect our findings to be very accurate for all , since for large components number () the ansatz (III) becomes exact.

## V The Spherical model case

In the infinite components limit () the classical model lies in the same universality of the spherical model, which is exactly solvable. Within the classical analogue of ansatz (III) the exact expressions for the critical exponents can be recovered Codello and D’Odorico (2013); Codello et al. (2015b). Such property also remains valid in the quantum case Nieuwenhuizen (1995); Vojta (1996) leading to

 z =σ2, (36) ν =22d−σ. (37)

These expressions, obtained from the exact solution of the flow equations in the limit, give the exact values of the critical exponents for the quantum spherical model Vojta (1996).

It should be noted that the same expression for the critical exponents could be obtained via the effective dimension relations. Indeed in the classical SR models we have and , then employing the effective dimension relation, equation (17), and multiplying by we immediately get (37). Thus the effective dimension relation is exact in the limit, as it is expected from the classical spherical model case Joyce (1966).

## Vi Conclusions

The results in figure 3 completes the necessary information to derive the scaling behavior of all the thermodynamic quantities of quantum LR rotor models. Indeed we determined both the dynamical critical exponent and the inverse correlation length exponent from which all the critical exponents can be obtained via the scaling relations. It should be noted that, in principle, one should also know the values for the anomalous dimension , which however in the LR system is found to be always given by the mean field result . In table 1 the values of the critical exponents for various values of and are reported.

Our motivations for the study presented in this paper were two-fold: on one side the progress in the control of atomic, molecular and optical systems made it possible to experimentally implement quantum long-range (LR) systems with tunable parameters, including the range of the interactions. These advancements pave the way to the study of the phase diagram and of the criticality quantum LR models, calling for the development of a unified treatment of their properties. On the other side, this motivates the search for the analogies and differences between quantum and classical LR systems, to understand what phenomena in classical systems typically due to the long-rangeness are also present in quantum systems and what are specifically connected to the presence of quantum fluctuations. Finally, a last motivation for the present work was to have a compact formalism to study different quantum LR systems at once, allowing also to clarify the relation between quantum LR models in dimensions and the corresponding anisotropic classical models in dimensions.

We therefore developed a general description of quantum LR models based on renormalization group (RG) techniques and capable to work for different dimensions , power-law decay exponents and groups of symmetry. We focused on the derivation of universal exponents describing the critical behavior of power-law decaying interacting components quantum rotor models. The case corresponds to a quantum LR Ising model, which is described by Hamiltonian (1), while the general cases, described by Hamiltonian (2), are quantum generalizations of the celebrated symmetric models.

LR power decaying interactions deeply affect the critical behavior of quantum models and their critical dynamics. We remind that the usual short-range (SR) quantum model with local propagator has a quantum critical point at zero temperature which lies in the same universality class of the classical phase transition in dimension Sachdev (2011). As discussed in the text, in the LR interacting case the low energy behavior of the propagator is not analytic in the momentum, being , and the quantum field theory obtained using Trotter decomposition is anisotropic in the extended dimensional space. However our analysis clearly shows that the problem is reduced to the study of a LR classical problem only for large enough, and that this is not the case for smaller than a critical value. For such small values of the correlation functions are genuinely strongly anisotropic in the spatial and time coordinates and the isotropy is not restored even at the criticality.

The mapping of a quantum SR universality into its classical equivalent obtained by Trotter decomposition is exact and so it is the relation between the universal behavior of a quantum LR model and its anisotropic classical equivalent. On the other hand, from the first line of equation (26) a correspondance emerges between the quantum LR system in dimension and its classical analogous in dimension .

Such correspondance naturally emerges in our framework, even if it appears not to be exact due to the difference in the kinetics sectors of the two models. This scenario is schematically drawn in figure 4, where the dashed lines stand for a non-exact mapping. Indeed the coefficients proportional to appearing in the second and third lines of equation (26) are substantially different from the ones of the classical case Defenu et al. (2015).

These coefficients are regulator dependent quantities whose precise value can be shown to not affect much the final result for the critical exponents, indeed both coefficients are proportional to which is never larger than and they can be shown to produce corrections to the critical exponents which are of order . We also explicitly solved the equations both in presence and absence of the two coefficients in the classical case Defenu et al. (2015) showing that the error arising from neglecting these two quantities is well below 1% in all the range, in the quantum case such error is expected to be smaller.

The results depicted in figure 3 are thus consistent with the naive expectation of equal universality class for quantum systems in dimension and classical systems in dimension . As discussed before even at our simple approximation level such equivalence is spoiled by the presence of the regulator dependent coefficients multiplying the non linear term in (26). However, our results for the correlation length exponent can be exactly mapped into the ones of classical SR models via effective dimension relation (31), since we disregarded the regulator dependents factors on the second line of equation (26). These results are in good agreement with the MC data taken from Sperstad et al. (2012), see figures 2 and 3, showing, once again, the reliability of the effective dimension approach Defenu et al. (2015).

The phase diagram of a LR system can be then represented in the plane, as in figure 1. For small enough values of the decay exponent the system undergoes spontaneous symmetry breaking with mean field exponents given by relations (14), (15) and (19), cyan shaded region in figure 1. For intermediate values of the system has peculiar LR exponents which continuously merge with the SR values at the threshold.

In analogy with the classical case Angelini et al. (2014); Defenu et al. (2015), such threshold value is given by , where is the anomalous dimension found in the pure SR system.

It is worth noting that for continuous symmetries it is also possible to identify a region, gray shaded area in figure 1(b), where no phase transition is possible. Such identification does not have a counterpart in the discrete symmetry case of the Ising , where no exact lower critical dimension is known even in the traditional SR case.

This theoretical picture advocates for more detailed comparison with ongoing experiments on LR power decaying interactions in quantum systems. Their understanding can pave the way to the comprehension of the role of dimension in quantum phase transitions and in the reliability of effective dimensional mappings.

Finally, we mention that a topic we did not address in detail in this paper is the emergence of the properties of the BKT phase transition in the , limit. We think this is a deserving subject of investigation, also in connection with the properties of the -dimensional XXZ with LR couplings Laflorencie et al. (2005); Bermudez et al. (2016).

Acknowledgements. N. Defenu would like to thank I. Frerot for the useful discussion about the relation between classical and quantum systems during the BEC 2016 conference. We are also pleased to acknowledge discussions with M. Ibáñez Berganza, G. Parisi, L. Tagliacozzo and S. Whitlock.

## References

• Saffman et al. (2009) M. Saffman, T. G. Walker,  and K. Molmer, “Quantum information with Rydberg atoms,” Rev. Mod. Phys. 82, 2313–2363 (2009)arXiv:0909.4777 .
• Lahaye et al. (2009) T Lahaye, C Menotti, L Santos, M Lewenstein,  and T Pfau, ‘‘The physics of dipolar bosonic quantum gases,” Reports Prog. Phys. 71, 126401 (2009)arXiv:0905.0386 .
• Ritsch et al. (2013) Helmut Ritsch, Peter Domokos, Ferdinand Brennecke,  and Tilman Esslinger, “Cold atoms in cavity-generated dynamical optical potentials,” Rev. Mod. Phys. 85, 553–601 (2013).
• Carr et al. (2009) Lincoln D Carr, David DeMille, Roman V Krems,  and Jun Ye, “Cold and ultracold molecules: science, technology and applications,” New J. Phys. 11, 055049 (2009).
• Blatt and Roos (2012) R. Blatt and C. F. Roos, “Quantum simulations with trapped ions,” Nat. Phys. 8, 277–284 (2012).
• Bloch et al. (2008) Immanuel Bloch, Jean Dalibard,  and Wilhelm Zwerger, “Many-body physics with ultracold gases,” Rev. Mod. Phys. 80, 885–964 (2008).
• Britton et al. (2012) Joseph W. Britton, Brian C. Sawyer, Adam C. Keith, C.-C. Joseph Wang, James K. Freericks, Hermann Uys, Michael J. Biercuk,  and John J. Bollinger, “Engineered two-dimensional Ising interactions in a trapped-ion quantum simulator with hundreds of spins,” Nature 484, 489–492 (2012).
• Schauß et al. (2012) Peter Schauß, Marc Cheneau, Manuel Endres, Takeshi Fukuhara, Sebastian Hild, Ahmed Omran, Thomas Pohl, Christian Gross, Stefan Kuhr,  and Immanuel Bloch, “Observation of spatially ordered structures in a two-dimensional Rydberg gas,” Nature 491, 87–91 (2012).
• Aikawa et al. (2012) K. Aikawa, A. Frisch, M. Mark, S. Baier, A. Rietzler, R. Grimm,  and F. Ferlaino, “Bose-Einstein Condensation of Erbium,” Phys. Rev. Lett. 108, 210401 (2012).
• Lu et al. (2012) Mingwu Lu, Nathaniel Q. Burdick,  and Benjamin L. Lev, “Quantum Degenerate Dipolar Fermi Gas,” Phys. Rev. Lett. 108, 215301 (2012).
• Yan et al. (2013) Bo Yan, Steven A. Moses, Bryce Gadway, Jacob P. Covey, Kaden R. A. Hazzard, Ana Maria Rey, Deborah S. Jin,  and Jun Ye, “Observation of dipolar spin-exchange interactions with lattice-confined polar molecules,” Nature 501, 521–525 (2013).
• Islam et al. (2013) R. Islam, C. Senko, W. C. Campbell, S. Korenblit, J. Smith, A. Lee, E. E. Edwards, C.-C. J. C. J. Wang, J. K. Freericks,  and C. Monroe, “Emergence and Frustration of Magnetism with Variable-Range Interactions in a Quantum Simulator,” Science (80-. ). 340, 583–587 (2013).
• Richerme et al. (2014) Philip Richerme, Zhe-Xuan Gong, Aaron Lee, Crystal Senko, Jacob Smith, Michael Foss-Feig, Spyridon Michalakis, Alexey V. Gorshkov,  and Christopher Monroe, “Non-local propagation of correlations in quantum systems with long-range interactions,” Nature 511, 198–201 (2014)arXiv:1401.5088 .
• Jurcevic et al. (2014) P. Jurcevic, B. P. Lanyon, P. Hauke, C. Hempel, P. Zoller, R. Blatt,  and C. F. Roos, “Quasiparticle engineering and entanglement propagation in a quantum many-body system,” Nature 511, 202–205 (2014).
• Douglas et al. (2015) J. S. Douglas, H. Habibian, C.-L. Hung, A. V. Gorshkov, H. J. Kimble,  and D. E. Chang, “Quantum many-body models with cold atoms coupled to photonic crystals,” Nat. Photonics 9, 326–331 (2015).
• Schempp et al. (2015) H. Schempp, G. Günter, S. Wüster, M. Weidemüller,  and S. Whitlock, “Correlated Exciton Transport in Rydberg-Dressed-Atom Spin Chains,” Phys. Rev. Lett. 115, 093002 (2015)arXiv:1504.01892 .
• Landig et al. (2015) Renate Landig, Ferdinand Brennecke, Rafael Mottl, Tobias Donner,  and Tilman Esslinger, “Measuring the dynamic structure factor of a quantum gas undergoing a structural phase transition,” Nat. Commun. 6, 7046 (2015)arXiv:1503.05565 .
• Landig et al. (2016) Renate Landig, Lorenz Hruby, Nishant Dogra, Manuele Landini, Rafael Mottl, Tobias Donner,  and Tilman Esslinger, “Quantum phases from competing short- and long-range interactions in an optical lattice,” Nature 532, 476–479 (2016)arXiv:1511.00007 .
• Martinez et al. (2016) Esteban A. Martinez, Christine A. Muschik, Philipp Schindler, Daniel Nigg, Alexander Erhard, Markus Heyl, Philipp Hauke, Marcello Dalmonte, Thomas Monz, Peter Zoller,  and Rainer Blatt, “Real-time dynamics of lattice gauge theories with a few-qubit quantum computer,” Nature 534, 516–519 (2016).
• Codello et al. (2015a) Alessandro Codello, Nicoló Defenu,  and Giulio D’Odorico, “Critical exponents of O (N) models in fractional dimensions,” Phys. Rev. D - Part. Fields, Gravit. Cosmol. 91, 105003 (2015a)arXiv:1410.3308v2 .
• Campa et al. (2014) A Campa, Thierry Dauxois, D Fanelli,  and S Ruffo, Physics of Long-Range Interacting Systems (2014) p. 410.
• Sachdev (2011) Subir Sachdev, Quantum Phase Transitions (Cambridge University Press, Cambridge, 2011) p. 517, arXiv:9811058 [cond-mat] .
• Defenu et al. (2015) Nicoló Defenu, Andrea Trombettoni,  and Alessandro Codello, “Fixed-point structure and effective fractional dimensionality for O( N ) models with long-range interactions,” Phys. Rev. E 92, 052113 (2015)arXiv:1211.3991 .
• Defenu et al. (2016) Nicolò Defenu, Andrea Trombettoni,  and Stefano Ruffo, “Anisotropic Long-Range Spin Systems,” Phys. Rev. B 94, 224411 (2016)arXiv:1606.07756 .
• Thouless (1969) D. J. Thouless, “Long-Range Order in One-Dimensional Ising Systems,” Phys. Rev. 187, 732–733 (1969).
• Yuval and Anderson (1970) G. Yuval and P. W. Anderson, “Exact results for the kondo problem: One-body theory and extension to finite temperature,” Phys. Rev. B 1, 1522–1528 (1970).
• Cardy (1981) J L Cardy, “One-dimensional models with interactions,” J. Phys. A. Math. Gen. 14, 1407–1415 (1981).
• Dutta and Bhattacharjee (2001) Amit Dutta and J. Bhattacharjee, ‘‘Phase transitions in the quantum Ising and rotor models with a long-range interaction,” Phys. Rev. B 64, 1–7 (2001).
• Sak (1973) J. Sak, “Recursion relations and fixed points for ferromagnets with long-range interactions,” Phys. Rev. B 8, 281–285 (1973).
• Picco (2012) Marco Picco, “Critical behavior of the Ising model with long range interactions,”  (2012)arXiv:1207.1018 .
• Blanchard et al. (2013) Thibault Blanchard, Marco Picco,  and M. A. Rajabpour, “Influence of long-range interactions on the critical behavior of the Ising model,” EPL 101, 56003 (2013)arXiv:1211.6758 .
• Brezin et al. (2014) E. Brezin, G. Parisi,  and F. Ricci-Tersenghi, “The Crossover Region Between Long-Range and Short-Range Interactions for the Critical Exponents,” J. Stat. Phys. 157, 855–868 (2014).
• Luijten and Blöte (2002) Erik Luijten and Henk W. J. Blöte, “Boundary between Long-Range and Short-Range Critical Behavior in Systems with Algebraic Interactions,” Phys. Rev. Lett. 89, 25703 (2002).
• Angelini et al. (2014) Maria Chiara Angelini, Giorgio Parisi,  and Federico Ricci-Tersenghi, “Relations between short-range and long-range Ising models,” Phys. Rev. E 89, 062120 (2014).
• Mori (2010) Takashi Mori, “Analysis of the exactness of mean-field theory in long-range interacting systems,” Phys. Rev. E 82, 060103 (2010)arXiv:1004.3622 .
• Horita et al. (2016) Toshiki Horita, Hidemaro Suwa,  and Synge Todo, “Upper and Lower Critical Decay Exponents of Ising Ferromagnets with Long-range Interaction,” Phys. Rev. E 95, 012143 (2016)arXiv:1605.09496 .
• Gori et al. (2017) G. Gori, M. Michelangeli, N. Defenu,  and A. Trombettoni, “One-dimensional long-range percolation: a numerical study,” Phys. Rev. E 96, 012108 (2017)arXiv:1610.00200 .
• Codello and D’Odorico (2013) Alessandro Codello and Giulio D’Odorico, “O(N) -Universality Classes and the Mermin-Wagner Theorem,” Phys. Rev. Lett. 110, 141601 (2013).
• Codello et al. (2015b) Alessandro Codello, Nicoló Defenu, Giulio D’Odorico,  and Giulio D’Odorico, “Critical exponents of models in fractional dimensions,” Phys. Rev. D 91, 105003 (2015b).
• Uhlenbeck et al. (1963) G. E. Uhlenbeck, P. C. Hemmer,  and M. Kac, “On the van der Waals Theory of the Vapor-Liquid Equilibrium. I. Discussion of a One-Dimensional Model,” J. Math. Phys. 4, 216 (1963).
• Behan et al. (2017a) Connor Behan, Leonardo Rastelli, Slava Rychkov,  and Bernardo Zan, “A scaling theory for the long-range to short-range crossover and an infrared duality,” To Appear  (2017a)arXiv:1703.05325v1 .
• Behan et al. (2017b) Connor Behan, Leonardo Rastelli, Slava Rychkov,  and Bernardo Zan, “Long-range critical exponents near the short-range crossover,” eprint arXiv:1703.03430  (2017b), 10.1103/PhysRevLett.118.241601arXiv:1703.03430 .
• Tetradis and Wetterich (1992) N. Tetradis and C. Wetterich, “Scale dependence of the average potential around the maximum in 4 theories,” Nucl. Phys. B 383, 197–217 (1992).
• Morris (1994) Tim R. Morris, “Derivative expansion of the exact renormalization group,” Phys. Lett. B 329, 241–248 (1994)arXiv:0108031 [hep-th] .
• Berges et al. (2002) Jürgen Berges, Nikolaos Tetradis,  and Christof Wetterich, “Non-perturbative renormalization flow in quantum field theory and statistical physics,” Phys. Rep. 363, 223–386 (2002).
• Mussardo (2010) G. Mussardo, Statistical field theory : an introduction to exactly solved models in statistical physics (Oxford University Press, 2010) p. 755.
• Wetterich (1993) Christof Wetterich, “Exact evolution equation for the effective potential,” Phys. Lett. B 301, 90–94 (1993).
• Joyce (1966) G. S. Joyce, “Spherical model with long-range ferromagnetic interactions,” Phys. Rev. 146, 349–358 (1966).
• Dutta (2003) Amit Dutta, “Effect of long-range interactions on the pure and random quantum Ising transitions,” in Phys. A Stat. Mech. its Appl., Vol. 318 (2003) pp. 63–71.
• Sperstad et al. (2012) Iver Bakken Sperstad, Einar B. Stiansen,  and Asle Sudbø, “Quantum criticality in spin chains with non-Ohmic dissipation,” Phys. Rev. B 85, 214302 (2012)arXiv:1204.2538 .
• Br??zin and Zinn-Justin (1976) E. Brézin and J. Zinn-Justin, ‘‘Renormalization of the nonlinear model in 2+ dimensions-application to the heisenberg ferromagnets,” Phys. Rev. Lett. 36, 691–694 (1976).
• Fey and Schmidt (2016) Sebastian Fey and Kai Phillip Schmidt, “Critical behavior of quantum magnets with long-range interactions in the thermodynamic limit,” Phys. Rev. B 94, 075156 (2016).
• Morris and Turner (1998) Tim R Morris and Michael D Turner, “Derivative Expansion of the Renormalization Group in O(N) Scalar Field Theory,” Nucl. Phys. B 509, 637–661 (1998)arXiv:9704202 [hep-th] .
• Nieuwenhuizen (1995) Th. M. Nieuwenhuizen, “Quantum Description of Spherical Spins,” Phys. Rev. Lett. 74, 4293–4296 (1995).
• Vojta (1996) Thomas Vojta, “Quantum version of a spherical model: Crossover from quantum to classical critical behavior,” Phys. Rev. B 53, 710 (1996).
• Laflorencie et al. (2005) Nicolas Laflorencie, Ian Affleck,  and Mona Berciu, “Critical phenomena and quantum phase transition in long range Heisenberg antiferromagnetic chains,” J. Stat. Mech. Theory Exp. 2005, P12001–P12001 (2005)arXiv:0509390 [cond-mat] .
• Bermudez et al. (2016) A. Bermudez, L. Tagliacozzo, G. Sierra,  and P. Richerme, “Long-range Heisenberg models in quasi-periodically driven crystals of trapped ions,” Phys. Rev. B 95, 024431 (2016)arXiv:1607.03337 .
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters