Quantum spherical model with competing interactions

We analyse the phase diagram of a quantum mean spherical model in terms of the temperature T$T$, a quantum parameter g$g$, and the ratio p=−_J2/_J1$p=-J_2/J_1$, where _J1>0$J_1>0$ refers to ferromagnetic interactions between first-neighbour sites along the d$d$ directions of a hypercubic lattice, and _J2<0$J_2<0$ is associated with competing antiferromagnetic interactions between second neighbours along m≤d$m\le d$ directions. We regain a number of known results for the classical version of this model, including the topology of the critical line in the g=0$g=0$ space, with a Lifshitz point at p=1/4$p=1/4$, for d>2$d>2$, and closed-form expressions for the decay of the pair correlations in one dimension. In the T=0$T=0$ phase diagram, there is a critical border, _gc=_gc(p)$g_c=g_c\left(p\right)$ for d≥2$d\ge 2$, with a singularity at the Lifshitz point if d<(m+4)/2$d<(m+4)/2$. We also establish upper and lower critical dimensions, and analyse the quantum critical behavior in the neighborhood of p=1/4$p=1/4$.     

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

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.     

Three PT-symmetric Hamiltonians with completely different spectra

We discuss three Hamiltonians, each with a central-field part H₀$H_0$ and a PT-symmetric perturbation igz$igz$. When H₀$H_0$ is the isotropic Harmonic oscillator the spectrum is real for all g$g$ because H$H$ is isospectral to H₀+g²/2$H_0+g^2/2$. When H₀$H_0$ is the Hydrogen atom then infinitely many eigenvalues are complex for all g$g$. If the potential in H₀$H_0$ is linear in the radial variable r$r$ then the spectrum of H$H$ exhibits real eigenvalues for 0<g<g_c$0<g<g_c$ and a PT phase transition at g_c$g_c$.     

How much information is needed to be the majority during the binary-state opinion formation?

In this paper, we try to propose a toy model, which follows the majority rule with the Fermi function, to uncover the role of the heterogeneous interaction between individuals in opinion formation. In order to do this, we define the impact factor I_Fi$I{F}_i$, says individual i$i$, as the exponential function of its connectivity _ki$k_i$ with the tunable parameter β$\beta$. β$\beta$ also shows the public information that can be collected by individuals in the system. We realize our model in scale-free networks with mean connectivity 〈k〉$〈k〉$. We find that much more public information (β>_β2$\beta >{\beta }_2$) and less public information (β<_β1$\beta <{\beta }_1$) cannot let either of the two opinions be the majority during the opinion formation. Furthermore, _β1${\beta }_1$ is a constant and equal to −0.76(±0.04)$-0.76(\pm 0.04)$, and _β2${\beta }_2$ decreases as a power-law function of the mean connectivity 〈k〉$〈k〉$ of the network. Our work can provide some perspectives and tools to understand the diversity of opinion in social networks.     

Surface and bulk properties of ballistic deposition models with bond breaking

We introduce a new class of growth models, with a surface restructuring mechanism in which impinging particles may dislodge suspended particles, previously aggregated on the same column in the deposit. The flux of these particles is controlled through a probability p$p$. These systems present a crossover, for small values of p$p$, from random to correlated (KPZ) growth of surface roughness, which is studied through scaling arguments and Monte Carlo simulations on one- and two-dimensional substrates. We show that the crossover characteristic time _t×$t_{\times }$ scales with p$p$ according to _t×∼p^−y$t_{\times }\sim p^{-y}$ with y=(n+1)$y=(n+1)$ and that the interface width at saturation _Wsat$W_{sat}$ scales as _Wsat∼p^−δ$W_{sat}\sim p^{-\delta }$ with δ=(n+1)/2$\delta =(n+1)/2$, where n$n$ is either the maximal number of broken bonds or of dislodged suspended particles. This result shows that the sets of exponents y=1$y=1$ and δ=1/2$\delta =1/2$ or y=2$y=2$ and δ=1$\delta =1$ found in all previous works focusing on systems with this same type of crossover are not universal. Using scaling arguments, we show that the bulk porosity P$P$ of the deposits scales as P∼p^y−δ$P\sim p^{y-\delta }$ for small values of p$p$. This general scaling relation is confirmed by our numerical simulations and explains previous results present in literature.     

Is space-time symmetry a suitable generalization of parity-time symmetry?

We discuss space-time symmetric Hamiltonian operators of the form H=H₀+igH^′$H=H_0+ig{H}^{\prime }$, where H₀$H_0$ is Hermitian and g$g$ real. H₀$H_0$ is invariant under the unitary operations of a point group G$G$ while H^′$H^{\prime }$ is invariant under transformation by elements of a subgroup G^′$G^{\prime }$ of G$G$. If G$G$ exhibits irreducible representations of dimension greater than unity, then it is possible that H$H$ has complex eigenvalues for sufficiently small nonzero values of g$g$. In the particular case that H$H$ is parity-time symmetric then it appears to exhibit real eigenvalues for all 0<g<g_c$0<g<g_c$, where g_c$g_c$ is the exceptional point closest to the origin. Point-group symmetry and perturbation theory enable one to predict whether H$H$ may exhibit real or complex eigenvalues for g>0$g>0$. We illustrate the main theoretical results and conclusions of this paper by means of two- and three-dimensional Hamiltonians exhibiting a variety of different point-group symmetries.     

Superconducting cosmic strings and one dimensional extended supersymmetric algebras

In this article we study in detail the supersymmetric structures that underlie the system of fermionic zero modes around a superconducting cosmic string. Particularly, we extend the analysis existing in the literature on the one dimensional N=2$N=2$ supersymmetry and we find multiple N=2$N=2$, d=1$d=1$ supersymmetries. In addition, compact perturbations of the Witten index of the system are performed and we find to which physical situations these perturbations correspond. More importantly, we demonstrate that there exists a much more rich supersymmetric structure underlying the system of fermions with N_f$N_f$ flavors and these are N$N$-extended supersymmetric structures with non-trivial topological charges, with “N$N$” depending on the fermion flavors.     

Singular unitarity in “quantization commutes with reduction”

Let M$M$ be a connected compact quantizable Kähler manifold equipped with a Hamiltonian action of a connected compact Lie group G$G$. Let M//G=?⁻¹(0)/G=_M0$M//G={?}^{-1}\left(0\right)/G=M_0$ be the symplectic quotient at value 0 of the moment map ?$?$. The space _M0$M_0$ may in general not be smooth. It is known that, as vector spaces, there is a natural isomorphism between the quantum Hilbert space over _M0$M_0$ and the G$G$-invariant subspace of the quantum Hilbert space over M$M$. In this paper, without any regularity assumption on the quotient _M0$M_0$, we discuss the relation between the inner products of these two quantum Hilbert spaces under the above natural isomorphism; we establish asymptotic unitarity to leading order in Planck’s constant of a modified map of the above isomorphism under a “metaplectic correction” of the two quantum Hilbert spaces.     

Post-breakthrough scaling in reservoir field simulation

We study the oil displacement and production behavior in an isothermal thin layered reservoir model subjected to water flooding. We use the CMG’s (Computer Modelling Group  ) numerical simulators to solve mass balance equations. The influences of the viscosity ratio (m≡_μoil/_μwater$m\equiv {\mu }_{oil}/{\mu }_{water}$) and the inter-well (injector-producer) distance r$r$ on the oil production rate C(t)$C\left(t\right)$ and the breakthrough time _tbr$t_{br}$ are investigated. Two types of reservoir configuration are used, namely one with random porosities and another with a percolation cluster structure. We observe that the breakthrough time follows a power-law of m$m$ and r$r$, _tbr∝r^αm^β$t_{br}\propto r^{\alpha }{m}^{\beta }$, with α=1.8$\alpha =1.8$ and β=−0.25$\beta =-0.25$ for the random porosity type, and α=1.0$\alpha =1.0$ and β=−0.2$\beta =-0.2$ for the percolation cluster type. Moreover, our results indicate that the oil production rate is a power law of time. In the percolation cluster type of reservoir, we observe that P(t)∝t^γ$P\left(t\right)\propto t^{\gamma }$, with γ=−1.81$\gamma =-1.81$, where P(t)$P\left(t\right)$ is the time derivative of C(t)$C\left(t\right)$. The curves related to different values of m$m$ and r$r$ may be collapsed suggesting a universal behavior for the oil production rate.     

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.     

Fluxmetric and magnetometric demagnetizing factors for cylinders

Fluxmetric and magnetometric demagnetizing factors, _Nf$N_{\mathrm{f}}$ and _Nm$N_{\mathrm{m}}$, for cylinders along the axial direction are numerically calculated as functions of material susceptibility χ$\chi$ and the ratio γ$\gamma$ of length to diameter. The results have an accuracy better than 0.1% with respect to min(_Nf,m,1-_Nf,m)$\min (N_{\mathrm{f},\mathrm{m}},1-N_{\mathrm{f},\mathrm{m}})$ and are tabulated in the range of 0.01?γ?500$0.01?\gamma ?500$ and -1?χ<∞$-1?\chi <\infty$. _Nm$N_{\mathrm{m}}$ along the radial direction is evaluated with a lower accuracy from _Nm$N_{\mathrm{m}}$ along the axis and tabulated in the range of 0.01?γ?1$0.01?\gamma ?1$ and -1?χ<∞$-1?\chi <\infty$. Some previous results are discussed and several applications are explained based on the new results.     

The Bethe lattice treatment of sound attenuation for a spin- 3/2 Ising model

The sound attenuation phenomena is investigated for a spin- 3/2 Ising model on the Bethe lattice in terms of the recursion relations by using the Onsager theory of irreversible thermodynamics. The dependencies of sound attenuation on the temperature (T$T$), frequency (w$w$), Onsager coefficient (γ$\gamma$) and external magnetic field (H$H$) near the second-order (_Tc)$\left(T_c\right)$ and first-order (_Tt)$\left(T_t\right)$ phase transition temperatures are examined for given coordination numbers q$q$ on the Bethe lattice. It is assumed that the sound wave couples to the order-parameter fluctuations which decay mainly via the order-parameter relaxation process, thus two relaxation times are obtained and which are used to obtain an expression for the sound attenuation coefficient (α)$\left(\alpha \right)$. Our investigations revealed that only one peak is obtained near _Tt$T_t$ and three peaks are found near _Tc$T_c$ when the Onsager coefficient is varied at a given constant frequency for q=3$q=3$. Fixing the Onsager coefficient and varying the frequency always leads to two peaks for q=3,4$q=3,4$ and 6 near _Tc$T_c$. The sound attenuation peaks are observed near _Tt$T_t$ at lower values of external magnetic field, but as it increases the sound attenuation peaks decrease and eventually disappear.     

A hidden analytic structure of the Rabi model

The Rabi model describes the simplest interaction between a cavity mode with a frequency ω_c${\omega }_c$ and a two-level system with a resonance frequency ω₀${\omega }_0$. It is shown here that the spectrum of the Rabi model coincides with the support of the discrete Stieltjes integral measure in the orthogonality relations of recently introduced orthogonal polynomials. The exactly solvable limit of the Rabi model corresponding to Δ=ω₀/(2ω_c)=0$\Delta ={\omega }_0/\left(2{\omega }_c\right)=0$, which describes a displaced harmonic oscillator, is characterized by the discrete Charlier polynomials in normalized energy ?$?$, which are orthogonal on an equidistant lattice. A non-zero value of Δ$\Delta$ leads to non-classical discrete orthogonal polynomials ?_k(?)${?}_k\left(?\right)$ and induces a deformation of the underlying equidistant lattice. The results provide a basis for a novel analytic method of solving the Rabi model. The number of ca. 1350 calculable energy levels per parity subspace obtained in double precision (cca 16 digits) by an elementary stepping algorithm is up to two orders of magnitude higher than is possible to obtain by Braak’s solution. Any first n$n$ eigenvalues of the Rabi model arranged in increasing order can be determined as zeros of ?_N(?)${?}_N\left(?\right)$ of at least the degree N=n+n_t$N=n+n_t$. The value of n_t>0$n_t>0$, which is slowly increasing with n$n$, depends on the required precision. For instance, n_t?26$n_t?26$ for n=1000$n=1000$ and dimensionless interaction constant κ=0.2$\kappa =0.2$, if double precision is required. Given that the sequence of the l$l$th zeros x_nl$x_{nl}$’s of ?_n(?)${?}_n\left(?\right)$’s defines a monotonically decreasing discrete flow with increasing n$n$, the Rabi model is indistinguishable from an algebraically solvable model in any finite precision. Although we can rigorously prove our results only for dimensionless interaction constant κ<1$\kappa <1$, numerics and exactly solvable example suggest that the main conclusions remain to be valid also for κ≥1$\kappa \ge 1$.     

On the integrability of symplectic Monge–Ampère equations

Let u$u$ be a function of n$n$ independent variables x¹,…,xⁿ$x^1,\dots ,x^n$, and let U=(_uij)$U=\left(u_{ij}\right)$ be the Hessian matrix of u$u$. The symplectic Monge–Ampère equation is defined as a linear relation among all possible minors of U$U$. Particular examples include the equation detU=1$detU=1$ governing improper affine spheres and the so-called heavenly equation, _u13u24−_u23u14=1$u_{13}{u}_{24}-u_{23}{u}_{14}=1$, describing self-dual Ricci-flat 4$4$-manifolds. In this paper we classify integrable symplectic Monge–Ampère equations in four dimensions (for n=3$n=3$ the integrability of such equations is known to be equivalent to their linearisability). This problem can be reformulated geometrically as the classification of ‘maximally singular’ hyperplane sections of the Plücker embedding of the Lagrangian Grassmannian. We formulate a conjecture that any integrable equation of the form F(_uij)=0$F\left(u_{ij}\right)=0$ in more than three dimensions is necessarily of the symplectic Monge–Ampère type.     

Tenth Peregrine breather solution to the NLS equation

In this paper we construct a particularly important solution to the focusing NLS equation, namely a Peregrine breather of the rank 10 which we call, P₁₀$P_{10}$ breather. The related explicit formula is given by the ratio of two polynomials of degree 110 with integer coefficients times trivial exponential factor. This formula drastically simplifies for the “initial values” namely for t=0$t=0$ or x=0$x=0$. This formula confirms a general conjecture saying that between all quasi-rational solutions of the rank N$N$ fixed by the condition that its absolute value tends to 1 at infinity and its highest maximum is located at the point (x=0,t=0)$(x=0,t=0)$, the P_N$P_N$ breather is distinguished by the fact that P_N(0,0)=2N+1$P_N(0,0)=2N+1$ and, in the aforementioned class of quasi-rational solutions, it is an absolute maximum. At the end we also make a few remarks concerning the rational deformations of P₁₀$P_{10}$ breather involving 2N−2$2N-2$ free real parameters chosen in a way that P_N$P_N$ breather itself corresponds to the zero values of these parameters although we have no intention to discuss the properties of these deformations here.