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.     

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$.     

Non-Abelian fusion rules from an Abelian system

We demonstrate the emergence of non-Abelian fusion rules for excitations of a two dimensional lattice model built out of Abelian degrees of freedom. It can be considered as an extension of the usual toric code model on a two dimensional lattice augmented with matter fields. It consists of the usual C(Z_p)$\mathbb{C}\left({\mathbb{Z}}_p\right)$ gauge degrees of freedom living on the links together with matter degrees of freedom living on the vertices. The matter part is described by a n$n$ dimensional vector space which we call H_n$H_n$. The Z_p${\mathbb{Z}}_p$ gauge particles act on the vertex particles and thus H_n$H_n$ can be thought of as a C(Z_p)$\mathbb{C}\left({\mathbb{Z}}_p\right)$ module. An exactly solvable model is built with operators acting in this Hilbert space. The vertex excitations for this model are studied and shown to obey non-Abelian fusion rules. We will show this for specific values of n$n$ and p$p$, though we believe this feature holds for all n>p$n>p$. We will see that non-Abelian anyons of the quantum double of C(S₃)$\mathbb{C}\left(S_3\right)$ are obtained as part of the vertex excitations of the model with n=6$n=6$ and p=3$p=3$. Ising anyons are obtained in the model with n=4$n=4$ and p=2$p=2$. The n=3$n=3$ and p=2$p=2$ case is also worked out as this is the simplest model exhibiting non-Abelian fusion rules. Another common feature shared by these models is that the ground states have a higher symmetry than Z_p${\mathbb{Z}}_p$. This makes them possible candidates for realizing quantum computation.     

Constraints on kinematic models from the latest observational data

Kinematical models are constrained by the latest observational data from geometry-distance measurements, which include 557 type Ia supernovae (SNIa) Union2 data and 15 observational Hubble data. Considering two parameterized deceleration parameter, the values of current deceleration parameter q₀$q_0$, jerk parameter j₀$j_0$ and transition redshift zT$z_T$, are obtained. Furthermore, we show the departures for two parameterized kinematical models from ΛCDM model according to the evolutions of jerk parameter j(z)$j(z)$. Also, it is shown that the constraint on jerk parameter j(z)$j(z)$ is weak by the current geometrical observed data.     

Thermodynamics properties of an anisotropic Ising antiferromagnet in both external longitudinal and transverse fields

We have studied the anisotropic two-dimensional nearest-neighbor Ising model with competitive interactions in both uniform longitudinal field H$H$ and transverse magnetic field Ω$\Omega$. Using the effective-field theory (EFT) with correlation in cluster with N=1$N=1$ spin we calculate the thermodynamic properties as a function of temperature with values H$H$ and Ω$\Omega$ fixed. The model consists of ferromagnetic interaction _Jx$J_x$ in the x$x$ direction and antiferromagnetic interaction _Jy$J_y$ in the y$y$ direction, and it is found that for H/_Jy∈[0,2]$H/J_y\in [0,2]$ the system exhibits a second-order phase transition. The thermodynamic properties are obtained for the particular case of λ=_Jx/_Jy=1$\lambda =J_x/J_y=1$ (isotropic square lattice).     

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.     

A note on embedding non-Abelian finite flavor groups in continuous groups

A non-Abelian finite flavor group G⊂SO(3)$G\subset \mathit{SO}(3)$ can have double covering G^′⊂SU(2)$G^{\prime }\subset \mathit{SU}(2)$ such that G⊄G^′$G\not\subset G^{\prime }$. This situation is not contradictory, but quite natural, and we give explicit examples such as G=Dn$G=D_n$, G^′=Q_2n$G^{\prime }=Q_{2n}$ and G=T$G=T$, G^′=T^′$G^{\prime }=T^{\prime }$. This observation can be crucial in particle theory model building.     

Pion radiative weak decays in nonlocal chiral quark models

We analyze the radiative pion decay π⁺→e⁺νeγ${\pi }^{+}\to e^{+}{\nu }_e\gamma$ within nonlocal chiral quark models that include wave function renormalization. In this framework we calculate the vector and axial-vector form factors FV$F_V$ and FA$F_A$ at q²=0$q^2=0$ — where q²$q^2$ is the e⁺νe$e^{+}{\nu }_e$ squared invariant mass — and the slope a   of FV(q²)$F_V(q^2)$ at q²→0$q^2\to 0$. The calculations are carried out considering different nonlocal form factors, in particular those taken from lattice QCD evaluations, showing a reasonable agreement with the corresponding experimental data. The comparison of our results with those obtained in the (local) NJL model and the relation of FV$F_V$ and a   with the form factor in π⁰→γ^?γ${\pi }^0\to {\gamma }^{?}\gamma$ decays are discussed.     

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.     

A new entropy based on a group-theoretical structure

A multi-parametric version of the nonadditive entropy S_q$S_q$ is introduced. This new entropic form, denoted by S_a,b,r$S_{a,b,r}$, possesses many interesting statistical properties, and it reduces to the entropy S_q$S_q$ for b=0$b=0$, a=r:=1−q$a=r:=1-q$ (hence Boltzmann–Gibbs entropy S_BG$S_{BG}$ for b=0$b=0$, a=r→0$a=r\to 0$). The construction of the entropy S_a,b,r$S_{a,b,r}$ is based on a general group-theoretical approach recently proposed by one of us, Tempesta (2016). Indeed, essentially all the properties of this new entropy are obtained as a consequence of the existence of a rational group law, which expresses the structure of S_a,b,r$S_{a,b,r}$ with respect to the composition of statistically independent subsystems. Depending on the choice of the parameters, the entropy S_a,b,r$S_{a,b,r}$ can be used to cover a wide range of physical situations, in which the measure of the accessible phase space increases say exponentially with the number of particles N$N$ of the system, or even stabilizes, by increasing N$N$, to a limiting value.     

Generalized Chaplygin gas model: Constraints from Hubble parameter versus redshift data

We examine observational constraints on the generalized Chaplygin gas (GCG) model for dark energy from the 9 Hubble parameter data points, the 115 SNLS Sne Ia data and the size of baryonic acoustic oscillation peak at redshift, z=0.35$z=0.35$. At a 95.4% confidence level, a combination of three data sets gives 0.67?As?0.83$0.67?A_s?0.83$ and −0.21?α?0.42$-0.21?\alpha ?0.42$, which is within the allowed parameters ranges of the GCG as a candidate of the unified dark matter and dark energy. It is found that the standard Chaplygin gas model (α=1$\alpha =1$) is ruled out by these data at the 99.7% confidence level.     

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.     

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.     

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}|$.     

Complex scale-free networks with tunable power-law exponent and clustering

We introduce a network evolution process motivated by the network of citations in the scientific literature. In each iteration of the process a node is born and directed links are created from the new node to a set of target nodes already in the network. This set includes m$m$ “ambassador” nodes and l$l$ of each ambassador’s descendants where m$m$ and l$l$ are random variables selected from any choice of distributions _pl$p_l$ and _qm$q_m$. The process mimics the tendency of authors to cite varying numbers of papers included in the bibliographies of the other papers they cite. We show that the degree distributions of the networks generated after a large number of iterations are scale-free and derive an expression for the power-law exponent. In a particular case of the model where the number of ambassadors is always the constant m$m$ and the number of selected descendants from each ambassador is the constant l$l$, the power-law exponent is (2l+1)/l$(2l+1)/l$. For this example we derive expressions for the degree distribution and clustering coefficient in terms of l$l$ and m$m$. We conclude that the proposed model can be tuned to have the same power law exponent and clustering coefficient of a broad range of the scale-free distributions that have been studied empirically.     

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.