Partition function zeros of the antiferromagnetic Ising model on triangular lattice in the complex temperature plane for nonzero magnetic field

The grand partition functions Z(T,B)$Z(T,B)$ of the Ising model on L×L$L\times L$ triangular lattices with fully periodic boundary conditions, as a function of temperature T and magnetic field B  , are evaluated exactly for L<12$L<12$ (using microcanonical transfer matrix) and approximately for L?12$L?12$ (using Wang–Landau Monte Carlo algorithm). From Z(T,B)$Z(T,B)$, the distributions of the partition function zeros of the triangular-lattice Ising model in the complex temperature plane for real B≠0$B\ne 0$ are obtained and discussed for the first time. The critical points aN(x)$a_N(x)$ and the thermal scaling exponents yt(x)$y_t(x)$ of the triangular-lattice Ising antiferromagnet, for various values of x=e−2βB$x=e^{-2\beta B}$, are estimated using the partition function zeros.     

Worm Monte Carlo study of the honeycomb-lattice loop model

We present a Markov-chain Monte Carlo algorithm of worm   type that correctly simulates the O(n)$O(n)$ loop model on any (finite and connected) bipartite cubic graph, for any real n>0$n>0$, and any edge weight, including the fully-packed limit of infinite edge weight. Furthermore, we prove rigorously that the algorithm is ergodic and has the correct stationary distribution. We emphasize that by using known exact mappings when n=2$n=2$, this algorithm can be used to simulate a number of zero-temperature Potts antiferromagnets for which the Wang–Swendsen–Kotecký cluster algorithm is non-ergodic, including the 3-state model on the kagome lattice and the 4-state model on the triangular lattice. We then use this worm algorithm to perform a systematic study of the honeycomb-lattice loop model as a function of n?2$n?2$, on the critical line and in the densely-packed and fully-packed phases. By comparing our numerical results with Coulomb gas theory, we identify a set of exact expressions for scaling exponents governing some fundamental geometric and dynamic observables. In particular, we show that for all n?2$n?2$, the scaling of a certain return time in the worm dynamics is governed by the magnetic dimension of the loop model, thus providing a concrete dynamical interpretation of this exponent. The case n>2$n>2$ is also considered, and we confirm the existence of a phase transition in the 3-state Potts universality class that was recently observed via numerical transfer matrix calculations.     

Hadro-charmonium

We argue that relatively compact charmonium states, J/ψ$J/\psi$, ψ(2S)$\psi (2S)$, χc${\chi }_c$, can very likely be bound inside light hadronic matter, in particular inside higher resonances made from light quarks and/or gluons. The charmonium state in such binding essentially retains its properties, so that the bound system decays into light mesons and the particular charmonium resonance. Thus such bound states of a new type, which we call hadro-charmonium, may explain the properties of some of the recently observed resonant peaks, in particular of Y(4.26)$Y(4.26)$, Y(4.32–4.36)$Y(4.32\text{–}4.36)$, Y(4.66)$Y(4.66)$, and Z(4.43)$Z(4.43)$. We discuss further possible implications of the suggested picture for the observed states and existence of other states of hadro-charmonium and hadro-bottomonium.     

Weak mirror symmetry of complex symplectic Lie algebras

A complex symplectic structure on a Lie algebra h$\mathfrak{h}$ is an integrable complex structure J$J$ with a closed non-degenerate (2,0)$(2,0)$-form. It is determined by J$J$ and the real part Ω$\Omega$ of the (2,0)$(2,0)$-form. Suppose that h$\mathfrak{h}$ is a semi-direct product g?V$\mathfrak{g}?V$, and both g$\mathfrak{g}$ and V$V$ are Lagrangian with respect to Ω$\Omega$ and totally real with respect to J$J$. This note shows that g?V$\mathfrak{g}?V$ is its own weak mirror image in the sense that the associated differential Gerstenhaber algebras controlling the extended deformations of Ω$\Omega$ and J$J$ are isomorphic.     

Zero finite-temperature charge stiffness within the half-filled 1D Hubbard model

Even though the one-dimensional (1D) Hubbard model is solvable by the Bethe ansatz, at half-filling its finite-temperature T>0$T>0$ transport properties remain poorly understood. In this paper we combine that solution with symmetry to show that within that prominent T=0$T=0$ 1D insulator the charge stiffness D(T)$D\left(T\right)$ vanishes for T>0$T>0$ and finite values of the on-site repulsion U$U$ in the thermodynamic limit. This result is exact and clarifies a long-standing open problem. It rules out that at half-filling the model is an ideal conductor in the thermodynamic limit. Whether at finite T$T$ and U>0$U>0$ it is an ideal insulator or a normal resistor remains an open question. That at half-filling the charge stiffness is finite at U=0$U=0$ and vanishes for U>0$U>0$ is found to result from a general transition from a conductor to an insulator or resistor occurring at U=U_c=0$U=U_c=0$ for all finite temperatures T>0$T>0$. (At T=0$T=0$ such a transition is the quantum metal to Mott-Hubbard-insulator transition.) The interplay of the η$\eta$-spin SU(2)$SU\left(2\right)$ symmetry with the hidden U(1)$U\left(1\right)$ symmetry beyond SO(4)$SO\left(4\right)$ is found to play a central role in the unusual finite-temperature charge transport properties of the 1D half-filled Hubbard model.     

If Gauss–Bonnet interaction plays the role of dark energy

A cosmological model has been constructed with Gauss–Bonnet-scalar interaction, where the Universe starts with exponential expansion but encounters infinite deceleration, q→∞$q\to \infty$ and infinite equation of state parameter, w→∞$w\to \infty$. During evolution it subsequently passes through the stiff fluid era, q=2$q=2$, w=1$w=1$, the radiation dominated era, q=1$q=1$, w=1/3$w=1/3$ and the matter dominated era, q=1/2$q=1/2$, w=0$w=0$. Finally, deceleration halts, q=0$q=0$, w=−1/3$w=-1/3$, and it then encounters a transition to the accelerating phase. Asymptotically the Universe reaches yet another inflationary phase q→−1$q\to -1$, w→−1$w\to -1$. Such evolution is independent of the form of the potential and the sign of the kinetic energy term, i.e., even a non-canonical kinetic energy is unable to phantomize (w<−1)$(w<-1)$ the model.     

Right unitarity triangles,stable CP-violating phases and approximate quark–lepton complementarity

Current experimental data indicate that two unitarity triangles of the CKM quark mixing matrix V   are almost the right triangles with α≈90°$\alpha \approx 90°$. We highlight a very suggestive parametrization of V and show that its CP-violating phase ? is nearly equal to α   (i.e., ?−α≈1.1°$?-\alpha \approx 1.1°$). Both ? and α   are stable against the renormalizaton-group evolution from the electroweak scale MZ$M_Z$ to a superhigh energy scale MX$M_X$ or vice versa, and thus it is impossible to obtain α=90°$\alpha =90°$ at MZ$M_Z$ from ?=90°$?=90°$ at MX$M_X$. We conjecture that there might also exist a maximal CP-violating phase φ≈90°$\phi \approx 90°$ in the MNS lepton mixing matrix U. The approximate quark–lepton complementarity relations, which hold in the standard parametrizations of V and U, can also hold in our particular parametrizations of V and U   simply due to the smallness of |V_ub|$|V_{ub}|$ and |V_e3|$|V_{e3}|$.     

Renormalization-group flow for the field strength in scalar self-interacting theories

We consider the renormalization-group coupled equations for the effective potential V(?)$V(?)$ and the field strength Z(?)$Z(?)$ in the spontaneously broken phase as a function of the infrared cutoff momentum k  . In the k→0$k\to 0$ limit, the numerical solution of the coupled equations, while consistent with the expected convexity property of V(?)$V(?)$, indicates a sharp peaking of Z(?)$Z(?)$ close to the end points of the flatness region that define the physical realization of the broken phase. This might represent further evidence in favor of the non-trivial vacuum field renormalization effect already discovered with variational methods.     

Bogolyubov approximation for diagonal model of an interacting Bose gas

We study, using the Bogolyubov approximation, the thermodynamic behavior of a superstable Bose system whose energy operator in the second-quantized form contains a nonlinear expression in the occupation numbers operators. We prove that for all values of the chemical potential satisfying μ>λ(0)$\mu >\lambda (0)$, where λ(0)?0$\lambda (0)?0$ is the lowest energy value, the system undergoes Bose–Einstein condensation.     

Vortices and Jacobian varieties

We investigate the geometry of the moduli space of N$N$ vortices on line bundles over a closed Riemann surface Σ$\Sigma$ of genus g>1$g>1$, in the little explored situation where 1≤N<g$1\le N<g$. In the regime where the area of the surface is just large enough to accommodate N$N$ vortices (which we call the dissolving limit), we describe the relation between the geometry of the moduli space and the complex geometry of the Jacobian variety of Σ$\Sigma$. For N=1$N=1$, we show that the metric on the moduli space converges to a natural Bergman metric on Σ$\Sigma$. When N>1$N>1$, the vortex metric typically degenerates as the dissolving limit is approached, the degeneration occurring precisely on the critical locus of the Abel–Jacobi map of Σ$\Sigma$ at degree N$N$. We describe consequences of this phenomenon from the point of view of multivortex dynamics.     

First-order density matrices in one dimension for independent fermions and impenetrable bosons in harmonic traps

To complement existing knowledge of the density matrix γF(x,y)${\gamma }_{\mathrm{F}}(x,y)$ of independent fermions for N   particles in one dimension under harmonic confinement, the corresponding matrix γIB(x,y)${\gamma }_{\mathrm{IB}}(x,y)$ for impenetrable bosons is given for N=2$N=2$ and 3 (with the N=4$N=4$ form available also). For fermions the momentum density is then obtained and illustrated numerically for N=10$N=10$. The boson momentum density is studied analytically at high momentum p  , the coefficients of the p⁻⁴$p^{\mathrm{-4}}$ and p⁻⁶$p^{\mathrm{-6}}$ terms being tabulated for N=2–5$N=2\text{–}5$ inclusive. Their dependence on powers of N   is exhibited numerically. Finally, the functional relationship between γIB(x,y)${\gamma }_{\mathrm{IB}}(x,y)$ and γF(x,y)${\gamma }_{\mathrm{F}}(x,y)$ is formally set out and illustrated.     

Comparative study between a two-dimensional anisotropic Heisenberg antiferromagnet with easy-axis single-ion anisotropy and one with easy-axis exchange anisotropy

Generally, in literature, easy-axis single ion anisotropy and easy-axis exchange anisotropy was treated in indistinct way. In this work we propose to perform a comparative study of the effects of these two easy-axis anisotropies on the behavior of the magnetization and the critical temperature (_Tc)$(T_{\mathrm{c}})$ in the 2D classical Heisenberg antiferromagnetic model. In order to study the low-temperature thermodynamics of this model, we should consider the contribution of anisotropic spin waves, using a self-consistent harmonic approximation (SCHA) theory. We compare the predictions of SCHA with numerical simulations on L×L$L\times L$ square lattices using Monte Carlo (MC) simulations, which include effects due to all thermodynamically allowed excitations. Our SCHA results are in good agreement with our MC simulations results and have shown that the strong K$K$ limit gives two different Ising-like behavior. In the exchange anisotropic case, the dependence of _Tc$T_{\mathrm{c}}$ on anisotropic parameter K$K$ becomes linear and in the single-ion anisotropic case, _Tc$T_{\mathrm{c}}$ becomes independent of K$K$. Also, using MC simulations and finite size scaling, we show that the critical exponents in the two anisotropic case are compatible with the 2D Ising values α=0.125$\alpha =0.125$ and γ=1.75$\gamma =1.75$.     

Lifshitz-point correlation length exponents from the large-n expansion

The large-n expansion is applied to the calculation of thermal critical exponents describing the critical behavior of spatially anisotropic d-dimensional systems at m  -axial Lifshitz points. We derive the leading non-trivial 1/n$1/n$ correction for the perpendicular correlation-length exponent _νL2${\nu }_{L2}$ and hence several related thermal exponents to order O(1/n)$O(1/n)$. The results are consistent with known large-n expansions for d  -dimensional critical points and isotropic Lifshitz points, as well as with the second-order epsilon expansion about the upper critical dimension d^?=4+m/2$d^{?}=4+m/2$ for generic m∈[0,d]$m\in [0,d]$. Analytical results are given for the special case d=4$d=4$, m=1$m=1$. For uniaxial Lifshitz points in three dimensions, 1/n$1/n$ coefficients are calculated numerically. The estimates of critical exponents at d=3$d=3$, m=1$m=1$ and n=3$n=3$ are discussed.     

The nature of ZZ branes

In minimal non-critical string theory we show that the principal (r,s)$(r,s)$ ZZ brane can be viewed as the basic (1,1)$(1,1)$ ZZ boundary state tensored with the (r,s)$(r,s)$ Cardy boundary state. In this sense there exists only one ZZ boundary state, the basic (1,1)$(1,1)$ boundary state.     

The puzzle of apparent linear lattice artifacts in the 2d non-linear σ-model and Symanzik's solution

Lattice artifacts in the 2d O(n) non-linear σ  -model are expected to be of the form O(a²)$\mathrm{O}(a^2)$, and hence it was (when first observed) disturbing that some quantities in the O(3)$\mathrm{O}(3)$ model with various actions show parametrically stronger cutoff dependence, apparently O(a)$\mathrm{O}(a)$, up to very large correlation lengths. In a previous letter Balog et al. (2009) [1] we described the solution to this puzzle. Based on the conventional framework of Symanzik's effective action, we showed that there are logarithmic corrections to the O(a²)$\mathrm{O}(a^2)$ artifacts which are especially large (ln³a${\ln }^{3}a$) for n=3$n=3$ and that such artifacts are consistent with the data. In this paper we supply the technical details of this computation. Results of Monte Carlo simulations using various lattice actions for O(3)$\mathrm{O}(3)$ and O(4)$\mathrm{O}(4)$ are also presented.