The Ising Model on a Quenched Ensemble of c=−5 Gravity Graphs

We study with Monte Carlo methods an ensemble of c=–5 gravity graphs, generated by coupling a conformal field theory with central charge c=–5 to two-dimensional quantum gravity. We measure the fractal properties of the ensemble, such as the string susceptibility exponent _s and the intrinsic fractal dimension d _H. We find _s=–1.5(1) and d _H=3.36(4), in reasonable agreement with theoretical predictions. In addition, we study the critical behavior of an Ising model on a quenched ensemble of the c=–5 graphs and show that it agrees, within numerical accuracy, with theoretical predictions for the critical behavior of an Ising model coupled dynamically to two-dimensional quantum gravity, with a total central charge of the matter sector c=–5.     

Symmetric vertex models on planar random graphs

We solve a 4-(bond)-vertex model on an ensemble of 3-regular (Φ³) planar random graphs, which has the effect of coupling the vertex model to 2D quantum gravity. The method of solution, by mapping onto an Ising model in field, is inspired by the solution by Wu et.al. of the regular lattice equivalent – a symmetric 8-vertex model on the honeycomb lattice, and also applies to higher valency bond vertex models on random graphs when the vertex weights depend only on bond numbers and not cyclic ordering (the so-called symmetric vertex models).The relations between the vertex weights and Ising model parameters in the 4-vertex model on Φ³ graphs turn out to be identical to those of the honeycomb lattice model, as is the form of the equation of the Ising critical locus for the vertex weights. A symmetry of the partition function under transformations of the vertex weights, which is fundamental to the solution in both cases, can be understood in the random graph case as a change of integration variable in the matrix integral used to define the model.Finally, we note that vertex models, such as that discussed in this paper, may have a role to play in the discretisation of Lorentzian metric quantum gravity in two dimensions.     

Fat Fisher zeroes

We show that it is possible to determine the locus of Fisher zeroes in the thermodynamic limit for the Ising model on planar (“fat”) φ⁴ random graphs and their dual quadrangulations by matching up the real part of the high and low temperature branches of the expression for the free energy. The form of this expression for the free energy also means that series expansion results for the zeroes may be obtained with rather less effort than might appear necessary at first sight by simply reverting the series expansion of a function g(z) which appears in the solution and taking a logarithm.Unlike regular 2D lattices where numerous unphysical critical points exist with non-standard exponents, the Ising model on planar φ⁴ graphs displays only the physical transition at c=exp(−2β)=1/4 and a mirror transition at c=−1/4 both with KPZ/DDK exponents (α=−1, β=1/2, γ=2). The relation between the φ⁴ locus and that of the dual quadrangulations is akin to that between the (regular) triangular and honeycomb lattices since there is no self-duality.     

Correlation functions in the multiple Ising model coupled to gravity

The model of p Ising spins coupled to 2d gravity, in the form of a sum over planar φ³ graphs, is studied and in particular the two-point and spin-spin correlation functions are considered. We first solve a toy model in which only a partial summation over spin configurations is performed and, using a modified geodesic distance, various correlation functions are determined. The two-point function has a diverging length scale associated with it. The critical exponents are calculated and it is shown that all the standard scaling relations apply. Next the full model is studied, in which all spin configurations are included. Many of the considerations for the toy model apply for the full model, which also has a diverging geometric correlation length associated with the transition to a branched polymer phase. Using a transfer function we show that the two-point and spin-spin correlation functions decay exponentially with distance. Finally, by assuming various scaling relations, we make a prediction for the critical exponents at the transition between the magnetized and branched polymer phases in the full model.     

Effective spin susceptibility, electron mass and Landé factor in two-dimensional (2D) Si inversion layers

We have studied the silicon (Si) band-structure, electron–electron and electron-ionized donor interaction effects on our accurate and approximate results (AcR and ApR) for renormalized effective spin susceptibitity (RESS), electron mass (EEM), Landé factor and spin polarization in the impure 2D Si (electron system), showing that:(i) our ApR, being strongly deviated from our AcR, reproduces approximately all the data obtained recently by Pudalov et al. (Phys. Rev. Lett. 88 (2002) 196404) [in particular, RESS =4.7 at the critical value of Wigner–Seitz radius r_s: r_s=r_c≈8.5 at which occur the “apparent” metal–insulator transition (MIT)] and can also be compared with other ApRs found in the recent literature,(ii) both the RESS and EEM produce physical singularities at the same critical value: r_s=r_c11.05661 (weakly disordered samples) at which occurs the “true” MIT; the existence of such two “apparent and true” critical values in this impure system agrees with a recent discussion by Abrahams et al. (Rev. Mod. Phys. 73 (2001) 251), and(iii) at r_s=r_c=8.5, at which occurs the “apparent” MIT, our AcR for effective spin polarization and the corresponding result, obtained using a disordered Hubbard model and a determinant quantum Monte Carlo method by Denteneer and Scalettar (Phys. Rev. Lett. 90 (2003) 246401), both give the same result: ξ_eff.=ξ_c0.31 at B0.4 T, which is found to be lower than the critical parallel magnetic field for full spin polarization, B_c=1.29 T, supporting thus the existence of such an “apparent” MIT.     

Disordered loop model and coupled conformal field theories

A family of models for fluctuating loops in a two-dimensional random background is analyzed. The models are formulated as O(n) spin models with quenched inhomogeneous interactions. Using the replica method, the models are mapped to the M→0 limits of M-layered O(n) models coupled each other via _1,3 primary fields. The renormalization group flow is calculated in the vicinity of the decoupled critical point, by an epsilon expansion around the Ising point (n=1), varying n as a continuous parameter. The one-loop beta function suggests the existence of a strongly coupled phase (0<n<n_*) near the self-avoiding walk point (n=0) and a line of infrared fixed points (n_*<n<1) near the Ising point. For the fixed points, the effective central charges are calculated. The scaling dimensions of the energy operator and the spin operator are obtained up to two-loop order. The relation to the random-bond q-state Potts model is briefly discussed.     

Wang-Landau study of the 3D Ising model with bond disorder

We implement a two-stage approach of the Wang-Landau algorithm to investigate the critical properties of the 3D Ising model with quenched bond randomness. In particular, we consider the case where disorder couples to the nearest-neighbor ferromagnetic interaction, in terms of a bimodal distribution of strong versus weak bonds. Our simulations are carried out for large ensembles of disorder realizations and lattices with linear sizes L in the range L=8-64L=8{-}64. We apply well-established finite-size scaling techniques and concepts from the scaling theory of disordered systems to describe the nature of the phase transition of the disordered model, departing gradually from the fixed point of the pure system. Our analysis (based on the determination of the critical exponents) shows that the 3D random-bond Ising model belongs to the same universality class with the site- and bond-dilution models, providing a single universality class for the 3D Ising model with these three types of quenched uncorrelated disorder.     

Diagrammatic decompositions on lattice for correlation functions of Ising spin models by iterative method

In the present paper the iterative method enables us to calculate correlation functions of Ising spin models by two approximate ways. The equations resulting from two models A and B, based on different physical considerations concerning correlations in surroundings of the reference spin, are solved for one, two, and three-dimensional Ising models with nearest neighbour interactions. Except of anomalies occuring in low-dimensional systems the model A leads to the critical pointA _c=0.191 for the three-dimensional cubic Ising model.     

HighT c superconductivity in the itinerant ising model: Spin singlet pairing

We study the spin singlet superconductivity exhibited in an itinerant Ising model Hamiltonian. This Hamiltonian models the Cu–O layers in highT _c oxide superconductors. Electrons are itinerant through nearest neighbor hopping. An Ising term is introduced to describe the antiferromagnetic superexchange interaction between electrons nominally on nearest neighbor Cu sites. We discuss various symmetry states allowed by the model, and give detailed predictions of the superconducting energy gap, specific heat, susceptibility, andT _c variation with carrier concentration. Results are compared to experimental data on highT _c superconductors and reasonable agreement is obtained.     

Spin glasses on thin graphs

In a recent paper [C. Baillie, D.A. Johnston and J.-P. Kownacki, Nucl. Phys. B 432 (1994) 551] we found strong evidence from simulations that the Ising antiferromagnet on “thin” random graphs — Feynman diagrams — displayed a mean-field spin-glass transition. The intrinsic interest of considering such random graphs is that they give mean-field theory results without long-range interactions or the drawbacks, arising from boundary problems, of the Bethe lattice. In this paper we reprise the saddle-point calculations for the Ising and Potts ferromagnet, antiferromagnet and spin glass on Feynman diagrams. We use standard results from bifurcation theory that enable us to treat an arbitrary number of replicas and any quenched bond distribution. We note the agreement between the ferromagnetic and spin-glass transition temperatures thus calculated and those derived by analogy with the Bethe lattice or in previous replica calculations.

We then investigate numerically spin glasses with a ±J bond distribution for the Ising and Q = 3, 4, 10, 50 state Potts models, paying particular attention to the independence of the spin-glass transition from the fraction of positive and negative bonds in the Ising case and the qualitative form of the overlap distribution P(q) for all of the models. The parallels with infinite-range spin-glass models in both the analytical calculations and simulations are pointed out.     

Neutron scattering study of two-dimensional Ising nature of K2CoF4

By neutron elastic and quasi-elastic scattering from K₂CoF₄, sublattice magnetization and spin correlations are examined quantitatively. The obtained critical exponents β, ν, γ, and η coincide precisely with those from Onsager's exact solution for the two-dimensional Ising model.     

The spin-1/2XXZ Heisenberg chain,the quantum algebra Uq[sl(2)], and duality transformations for minimal models

The finite-size scaling spectra of the spin-1/2XXZ Heisenberg chain with toroidal boundary conditions and an even number of sites provide a projection mechanism yielding the spectra of models with a central chargec < 1, including the unitary and nonunitary minimal series. Taking into account the half-integer angular momentum sectors—which correspond to chains with an odd number of sites—in many cases leads to new spinor operators appearing in the projected systems. These new sectors in theXXZ chain correspond to new types of frustration lines in the projected minimal models. The corresponding new boundary conditions in the Hamiltonian limit are investigated for the Ising model and the 3-state Potts model and are shown to be related to duality transformations which are an additional symmetry at their self-dual critical point. By different ways of projecting systems we find models with the same central charge sharing the same operator content and modular invariant partition function which, however, differ in the distribution of operators into sectors and hence in the physical meaning of the operators involved. Related to the projection mechanism in the continuum there are remarkable symmetry properties of the finiteXXZ chain. The observed degeneracies in the energy and momentum spectra are shown to be the consequence of intertwining relations involvingU _q [sl(2)] quantum algebra transformations.     

The O(n) Loop Model on the 3–12 Lattice

The partition function of the O(n) loop model on the honeycomb lattice is mapped to that of the O(n) loop model on the 3–12 lattice. Both models share the same operator content and thus critical exponents. The critical points are related via a simple transformation of variables. When n = 0 this gives the recently found exact value = 1.711041... for the connective constant of self-avoiding walks on the 3–12 lattice. The exact critical points are recovered for the Ising model on the 3–12 lattice and the dual asanoha lattice at n = 1.     

Renormalization group for matrix models with branching interactions

We develop a method to obtain the large-N renormalization group flows for matrix models of two-dimensional gravity plus branched polymers. This method gives precise results for the critical points and exponents for one-matrix models. We show that it can be generalized to two-matrix models and we recover the Ising critical points.     

Computer simulation studies of magnetic phase transitions

The advancements which have been made in the use of computer simulations to study magnetic-phase transitions and critical phenomena are reviewed. We describe how the use of a combination of sophisticated Monte Carlo simulation algorithms and reweighting (histogram) techniques have allowed the determination of the static critical behavior with unprecedented precision. The study of “dynamic” critical behavior in simple spin models by both Monte Carlo and spin dynamics methods is also reviewed. Recent estimates for dynamic critical exponents are given including those for true dynamics.     

A hierarchically constrained kinetic Ising model

A concrete model for hierarchically constrained dynamics in the sense proposed by Palmer et al. (Phys. Rev. Lett.53, 958 (1984)) is presented. The model is a kinetic Ising chain with an asymmetric kinetic constraint, allowing a spin to flip only if its neighbour to the right is in the up spin state. The spin autocorrelation function is obtained by numerically exact calculation for finite chain length up toL=9 and by Monte Carlo simulation for effectively infinite chain length. The Kohlrausch-Williams-Watts formula is found to fit the results only with limited accuracy, and within limited time intervals. We also performed an analytical calculation using an effective-medium approximation. The approximation leads to a spurious blocking transition at a critical up spin concentrationc=0.5.     

Study of the three-dimensional Ising model on film geometry with the cluster Monte Carlo method

Topological properties of clusters are used to extract critical parameters. This method is tested for the bulk properties ofd=2 percolation and thed=2, 3 Ising model. For the latter we obtain an accurate value of the critical temperatureJ/k _B T _c=0.221617(18). In the case of thed=3 Ising model with film geometry the critical value of the surface coupling at the special transitions is determined as J_1c/J=1.5004(20) together with the critical exponents ₁ ^m =0.237(5) and=0.461(15).     

Renormalization group approach to spin glass systems

A renormalization group transformation suitable for spin glass models and, more generally, for disordered models, is presented. The procedure is non-standard in both the nature of the additional interactions and the coarse graining transformation, that is performed on the overlap probability measure. Universality classes are thus naturally defined on a large set of models, going from and Gaussian spin glasses to Ising and fully frustrated models, and others. The proposed analysis is tested numerically on the Edwards-Anderson model in d = 4. Good estimates of the critical index ν and of T _c are obtained, and an RG flow diagram is sketched for the first time. Received 17 November 2000     

Complex Gaussian Multiplicative Chaos

In this article, we study complex Gaussian multiplicative chaos. More precisely, we study the renormalization theory and the limit of the exponential of a complex log-correlated Gaussian field in all dimensions (including Gaussian Free Fields in dimension 2). Our main working assumption is that the real part and the imaginary part are independent. We also discuss applications in 2D string theory; in particular we give a rigorous mathematical definition of the so-called Tachyon fields, the conformally invariant operators in critical Liouville Quantum Gravity with a c = 1 central charge, and derive the original KPZ formula for these fields.     

2D quantum gravity from quantum entanglement

In quantum systems with many degrees of freedom the replica method is a useful tool to study the entanglement of arbitrary spatial regions. We apply it in a way that allows them to backreact. As a consequence, they become dynamical subsystems whose position, form, and extension are determined by their interaction with the whole system. We analyze, in particular, quantum spin chains described at criticality by a conformal field theory. Its coupling to the Gibbs' ensemble of all possible subsystems is relevant and drives the system into a new fixed point which is argued to be that of the 2D quantum gravity coupled to this system. Numerical experiments on the critical Ising model show that the new critical exponents agree with those predicted by the formula of Knizhnik, Polyakov, and Zamolodchikov.