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.     

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.     

Transition between SU(4) and SU(2) Kondo effect

Motivated by experiments in nanoscopic systems, we study a generalized Anderson, which consist of two spin degenerate doublets hybridized to a singlet by the promotion of an electron to two conduction bands, as a function of the energy separation δ$\delta$ between both doublets. For δ=0$\delta =0$ or very large, the model is equivalent to a one-level SU(N$N$) Anderson model, with N=4$N=4$ and 2 respectively. We study the evolution of the spectral density for both doublets (ρ_1σ(ω)${\rho }_{1\sigma }(\omega )$ and ρ_2σ(ω)${\rho }_{2\sigma }(\omega )$) and their width in the Kondo limit as δ$\delta$ is varied, using the non-crossing approximation (NCA). As δ$\delta$ increases, the peak at the Fermi energy in the spectral density (Kondo peak) splits and the density of the doublet of higher energy ρ_2σ(ω)${\rho }_{2\sigma }(\omega )$ shifts above the Ferrmi energy. The Kondo temperature T_K (determined by the half-width at half maximum of the Kondo peak in density of the doublet of lower energy ρ_1σ(ω)${\rho }_{1\sigma }(\omega )$) decreases dramatically. The variation of T_K with δ$\delta$ is reproduced by a simple variational calculation.     

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.     

A Monte Carlo model for networks between professionals and society

We propose a network model with a fixed number of nodes and links and with a dynamic which favors links between nodes differing in connectivity. We observe a phase transition and parameter regimes with degree distributions following power laws, P(k)∼k^-γ$P(k)\sim k^{-\gamma }$, with γ$\gamma$ ranging from 0.2$0.2$ to 0.5$0.5$, small-world properties, with a network diameter following D(N)∼logN$D(N)\sim \log N$ and relative high clustering, following C(N)∼1/N$C(N)\sim 1/N$ and C(k)∼k^-α$C(k)\sim k^{-\alpha }$, with α$\alpha$ close to 3. We compare our results with data from real-world protein interaction networks.     

On complete spacelike hypersurfaces with two distinct principal curvatures in Lorentz–Minkowski space

We investigate complete spacelike hypersurfaces in Lorentz–Minkowski space with two distinct principal curvatures and constant m$m$th mean curvature. By using Otsuki’s idea, we obtain the global classification result. As their applications, we obtain some characterizations for hyperbolic cylinders. We prove that the only complete spacelike hypersurfaces in Lorentz–Minkowski (n+1)$(n+1)$-spaces (n≥3$n\ge 3$) of nonzero constant m$m$th mean curvature (m≤n−1$m\le n-1$) with two distinct principal curvatures λ$\lambda$ and μ$\mu$ satisfying inf(λ−μ)²>0$inf{(\lambda -\mu )}^2>0$ are the hyperbolic cylinders. We also obtain a global characterization for hyperbolic cylinder Hⁿ⁻¹(c)×R${\mathbb{H}}^{n-1}\left(c\right)\times \mathbb{R}$ in terms of square length of the second fundamental form.     

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.     

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

Green functions and twist correlators for N branes at angles

We compute the Green functions and correlator functions for N   twist fields for branes at angles on T²$T^2$ and we show that there are N−2$N-2$ different configurations labeled by an integer M   which is roughly associated with the number of reflex angles of the configuration. In order to perform this computation we use an SL(2,R)$\mathit{SL}(2,\mathbb{R})$ invariant formulation and geometric constraints instead of Pochammer contours. In particular the M=1$M=1$ or M=N−1$M=N-1$ amplitude can be expressed without using transcendental functions. We determine the amplitudes normalization from N→N−1$N\to N-1$ reduction without using the factorization into the untwisted sector. Both the amplitudes normalization and the OPE of two twist fields are unique (up to one constant) when the ?↔1−?$?\leftrightarrow 1-?$ symmetry is imposed. For consistency we find also an infinite number of relations among Lauricella hypergeometric functions.     

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

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

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

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.     

Exact height distributions for the KPZ equation with narrow wedge initial condition

We consider the KPZ equation in one space dimension with narrow wedge initial condition, h(x,t=0)=−|x|/δ$h(x,t=0)=-|x|/\delta$, δ?1$\delta ?1$, evolving into a parabolic profile with superimposed fluctuations. Based on previous results for the weakly asymmetric simple exclusion process with step initial conditions, we obtain a determinantal formula for the one-point distribution of the solution h(x,t)$h(x,t)$ valid for any x   and t>0$t>0$. The corresponding distribution function converges in the long time limit, t→∞$t\to \infty$, to the Tracy–Widom distribution. The first order correction is a shift of order t−1/3$t^{-1/3}$. We provide numerical computations based on the exact formula.     

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.     

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.     

Local derivative post-processing for the discontinuous Galerkin method

Obtaining accurate approximations for derivatives is important for many scientific applications in such areas as fluid mechanics and chemistry as well as in visualization applications. In this paper we discuss techniques for computing accurate approximations of high-order derivatives for discontinuous Galerkin solutions to hyperbolic equations related to these areas. In previous work, improvement in the accuracy of the numerical solution using discontinuous Galerkin methods was obtained through post-processing by convolution with a suitably defined kernel. This post-processing technique was able to improve the order of accuracy of the approximation to the solution of time-dependent symmetric linear hyperbolic partial differential equations from order k+1$k+1$ to order 2k+1$2k+1$ over a uniform mesh; this was extended to include one-sided post-processing as well as post-processing over non-uniform meshes. In this paper, we address the issue of improving the accuracy of approximations to derivatives of the solution by using the method introduced by Thomée [19]. It consists in simply taking the α$\alpha$th-derivative of the convolution of the solution with a sufficiently smooth kernel. The order of convergence of the approximation is then independent   of the order of the derivative, |α|$|\alpha |$. We also discuss an efficient way of computing the approximation which does not involve differentiation but the application of simple finite differencing. Our results show that the above-mentioned approximations to the α$\alpha$th-derivative of the exact solution of linear, multidimensional symmetric hyperbolic systems obtained by the discontinuous Galerkin method with polynomials of degree k$k$ converge with order 2k+1$2k+1$ regardless of the order |α|$|\alpha |$ of the derivative.