Boundary CFT and tensor network approach to surface critical phenomena of the tricritical 3-state Potts model

It is widely known that there is no spontaneous symmetry breaking in 1d classical systems at finite temperature, which suggests that 1d classical many body problems are almost trivial and boring.
Then, let us consider 1d classical `edges’ in 2d systems, or 1d classical systems attached to a 2d bulk.
Does something interesting happen in these 1d edges?
In fact, the situation is a bit different from the simple 1d cases.

When the 2d bulk is at criticality, in particular, the strong correlation at the bulk helps the edge ordering, which results in a variety of phase transitions on the 1d edges.
For example, the 2d classical tricritical Ising model exhibits a finite-surface-temperature symmetry breaking on its 1d edge, when the bulk is fine-tuned at the tricritical point [1].
It is also known that another surface transition occurs driven by the surface external fields.

As a theoretical way of studying surface criticality, boundary conformal field theory (BCFT) is very powerful, particularly for 2d classical / 1+1d quantum systems at criticality.
It could be essential to understand surface critical behavior precisely, since the restriction of the conformal invariance may enable us not only to classify the conformal boundary conditions but also to investigate stability of boundary fixed points and exact values of scaling dimensions.

Now, let us consider an extension of the above classical tricritical Ising model in 2d, called the tricritical 3-state Potts model, whose symmetry of the spin is \(S_3\) rather than \(Z_2\).
The Monte Carlo simulation of this model suggests, similarly to the tricritical Ising, rich surface phase transitions when the bulk is at tricriticality [2].
However, the precise study from the view point of BCFT has been missing for a long time.

In order to investigate the more detailed surface phase diagram of this model, we utilize the \(ADE\) classification of the minimal-series BCFTs and perform numerical simulation with tensor network renormaization (TNR) [3,4,5].
The tricritical 3-state Potts model can be described by the \(D\)-type minimal CFT with \(c=6/7\), whose modular invariant partition function is labeled as the pair of Lie algebra \((D_4, A_6)\) in the \(ADE\) classification.
Using the complete classification of the conformal boundary states in Ref. [3], one can derive the twelve conformal boundary conditions labeled by the nodes of the Dynkin diagrams \(D_4\) and \(A_6\).
The triality of \(D_4\) implies that those can be classified into three \(Z_3\)-symmetric boundary states and nine \(Z_3\)-broken ones.

To understand the physical picture of the above twelve boundary fixed points, we simulate the 3-state dilute Potts model on lattices with the TNR technique.
The extended TNR method suited for open-boundary systems enables us to extract accurate conformal spectrum from lattice models, by which one can study the correspondence between the boundary fixed points from BCFT and the surface phases in the phase diagram of the tricritical 3-state Potts model.
Our numerical analysis reveals that the eleven boundary fixed points among the twelve can be realized on the lattice by controlling the external field and coupling strength at the boundary.
The last unfound fixed point would be out of the physically sound region in the parameter space, and we leave it an open question to consider the realization of that boundary condition on the lattice model.
(by Shumpei Iino)

Fig. The surface phase diagram of the tricritical 3-state Potts model obtained by the TNR simulation.

References:
[1] I. Affleck, J. Phys. A 33(37), 6473 (2000).
[2] Y. Deng and H. W. J. Blöte, Phys. Rev. E 70, 035107 (2004); Phys. Rev. E 71, 026109 (2005).
[3] R. E. Behrend, P. A. Pearce, V. B. Petkova, and J.-B. Zuber, Nucl. Phys. B 579(3), 707 (2000).
[4] G. Evenbly and G. Vidal, Phys. Rev. Lett. 115, 180405 (2015); S. Iino, S. Morita, and N. Kawashima, Phys. Rev. B 101, 155418 (2020).
[5] S. Iino, arXiv:2007.03182; J. Stat. Phys. 182, 56 (2021).

SU(N) Heisenberg model with multicolumn representations

The SU(N) symmetric antiferromagnetic Heisenberg model with multicolumn representations on the two-dimensional square lattice is investigated by quantum Monte Carlo simulations. For the representation of a Young diagram with two columns, we confirm that a valence-bond solid (VBS) order appears as soon as the Néel order disappears at N=10, indicating no intermediate phase. In the case of the representation with three columns, there is no evidence for either the Néel or the VBS ordering for N15. This is actually consistent with the large-N theory, which predicts that the VBS state immediately follows the Néel state, because the expected spontaneous order is too weak to be detected.

(by Tsuyoshi Okubo)

Reference:

T. Okubo, K. Harada, J. Lou and N. Kawashima: Phys. Rev. B 92, 134404 (2015).

First-Order Transition to Incommensurate Phase in Frustrated Heisenberg Model

We study a finite-temperature phase transition in the two-dimensional classical Heisenberg model on a triangular lattice with a ferromagnetic nearest-neighbor interaction J1 and an antiferromagnetic third-nearest-neighbor interaction J3 using a Monte Carlo method. Apart from a trivial degeneracy corresponding to O(3) spin rotations, the ground state for J3 ≠0 has a threefold degeneracy corresponding to 120° lattice rotations. We find that this model exhibits a first-order phase transition with the breaking of the threefold symmetry when the interaction ratio is J3 / J1 = -3.

(by Ryo TAMURA)

Reference:

R. Tamura and N. Kawashima,  J. Phys. Soc. Jpn. 77, 103002 (2008).

Asymptotic Dimensional Reduction near a Quantum Critical Point of BCT XY Model

The principal effect of frustration is cancellation among competing couplings. In quasi-two-dimensional systems, when the frustration exists in inter-layer couplings, it tends to enhance murual independence of layers. At finite temperature, however, this cancellation can never be perfect so as to make the system purely two-dimensional because thermal fluctuation generates effective net interlayer couplings. Even at zero temperature, the geometric frustration fails to decouple layers completely since the zero-point fluctuations restore the inter-layer coupling as discussed by Maltseva and Coleman. As a result, the critical phenomena in any real quasi-two-dimensional system belongs to the universality class of some three dimensional theory no matter how perfectly the frustration may seem to cancel the couplings. However, we find that this is not necessarily true at quantum critical points. [1]

The first experimental evidence of this phenomenon was found by measuring critical exponents of a field induced QCP in BaCuSi2O6. [2]  This is a spin dimer system whose highly symmetric crystal structure gives it unique advantages for tackling the fundamental role of dimensionality in quantum criticality. The material consists of layers stacked on top of each other forming a body-centered-cubic (BCT) lattice of dimers. By neglecting two high-energy triplet states of each dimer, the system can be represented by an S=1/2 XY model on BCT lattice. We examined the magnon excitations by using the spin-wave treatment..[1] The effective coupling between layers is proportional to the number of excited magnons, which means that the fluctuation responsible for the effective inter-layer coupling dies out as we approach zero temperature and the quantum critical point can possess purely two dimensional characteristics. Thermal or quantum fluctuations at finite temperature induce a crossover to d = 3 away from the QCP. The two-dimensional value of the exponent characterizes the phase boundary near the zero temperature, in agreement with the experiment [2].

Figure 1: A BCT lattice

The agreement between the pure two-dimensional behavior and the real experiment is not only qualitative, but it is quantitatively good as can be seen in Fig. 2. The upper panel shows the critical field measured in the experiment on BaCuSi2O6 plotted with the phase boundary obtained by quantum Monte Carlo simulation assuming that the inter-layer coupling is zero (filled circles) . Both the results show the two-dimensional quantum criticality that is characterized by the asymptotic linear dependence. The correspondence is further confirmed by other quantities: the field-dependence of the magnetization at zero (or very low) temperature  (lower left panel), and the temperature-dependence of the magnetization at critical magnetic field (lower right panel).

Figure 2: Static properties of BCT XY model. The phase boundary (top), the magnetization vs the field (bottom left) and the temperature (bottom right).

(by Naoki KAWASHIMA)

Reference:

[1] C. D. Batista,  J. Schmalian,  N. Kawashima, P. Sengupta, S. E. Sebastian, N. Harrison, M. Jaime and I. R. Fisher, Phys. Rev. Lett. 98, 257201 (2007).
[2] S. E. Sebastian, N. Harrison, C. D. Batista, L. Balicas, M. Jaime, P. A. Sharma, N. Kawashima and I. R. Fisher, Nature 441, 617 (2006).