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.     

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.     

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.     

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.     

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.     

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.     

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

Domain wall propagation in Permalloy nanowires with a thickness gradient

The magnetization reversal behavior of Permalloy nanowires has been investigated using a magneto-optic Kerr effect setup. Nanowires with various widths, w=250$w=250$ nm to 3 μm and a thickness of t=10$t=10$ nm were fabricated by electron-beam lithography and subsequent lift-off. Furthermore, similar nanowires but with a thickness gradient along the nanowire axis have been prepared to investigate the influence of the gradient on the magnetic domain wall propagation. Magnetization hysteresis loops recorded on individual nanowires without a gradient are compared to corresponding wires with a thickness gradient. The dependence of the coercive field, _Hc$H_c$ vs. t/w$t/w$ shows a linear behavior for wires without a gradient. However, wires with a gradient display a more complex crossover behavior. We find a plateau in the _Hc$H_c$ vs. t/w$t/w$ curve at values of w$w$, where a transformation from transverse to vortex domain wall type is expected.     

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.     

Generalized entanglement constraints in multi-qubit systems in terms of Tsallis entropy

We provide generalized entanglement constraints in multi-qubit systems in terms of Tsallis entropy. Using quantum Tsallis entropy of order q$q$, we first provide a generalized monogamy inequality of multi-qubit entanglement for q=2$q=2$ or 3$3$. This generalization encapsulates the multi-qubit CKW-type inequality as a special case. We further provide a generalized polygamy inequality of multi-qubit entanglement in terms of Tsallis-q$q$ entropy for 1≤q≤2$1\le q\le 2$ or 3≤q≤4$3\le q\le 4$, which also contains the multi-qubit polygamy inequality as a special case.     

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.     

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.     

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.     

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.     

Universal behavior in a model of city traffic with unequal green/red times

The effect of green/red asymmetry is studied for the single-car traffic model proposed in [B.A. Toledo, V. Muñoz, J. Rogan, C. Tenreiro, J.A. Valdivia, Modeling traffic through a sequence of traffic lights, Phys. Rev. E 70 (1) (2004) 016107], on two different signal synchronization strategies, namely, all signals in phase, and a green wave. The asymmetry is characterized by the parameter g=_tgr/T$g=t_{\text{gr}}/T$, where _tgr$t_{\text{gr}}$ is the green time and T$T$ the signal period. Although the car dynamics turns simpler or more complex, as compared with the equivalent situation for the symmetric case g=0.5$g=0.5$, critical behavior around resonance is shown to be preserved. However, unlike the case g=0.5$g=0.5$, critical exponents at both sides of the resonance are not equal and depend on g$g$. Analytical expressions for them are found, and shown to be both consistent with simulation results and independent of the distribution of distances between signals for the green wave case. Also, it is found that the green wave strategy is more robust to changes in g$g$, with respect to the synchronized lights strategy, in the sense that larger departures from g=0.5$g=0.5$ are needed to have noticeable effects on the car dynamics.     

Random matrix theory for transition strength densities in finite quantum systems: Results from embedded unitary ensembles

Embedded random matrix ensembles are generic models for describing statistical properties of finite isolated interacting quantum many-particle systems. For the simplest spinless fermion (or boson) systems, with say m$m$ fermions (or bosons) in N$N$ single particle states and interacting via k$k$-body interactions, we have EGUE(k$k$) [embedded GUE of k$k$-body interactions] with GUE embedding and the embedding algebra is U(N)$U\left(N\right)$. A finite quantum system, induced by a transition operator, makes transitions from its states to the states of the same system or to those of another system. Examples are electromagnetic transitions (then the initial and final systems are same), nuclear beta and double beta decay (then the initial and final systems are different), particle addition to or removal from a given system and so on. Towards developing a complete statistical theory for transition strength densities (transition strengths multiplied by the density of states at the initial and final energies), we have derived formulas for the lower order bivariate moments of the strength densities generated by a variety of transition operators. Firstly, for a spinless fermion system, using EGUE(k$k$) representation for a Hamiltonian that is k$k$-body and an independent EGUE(t$t$) representation for a transition operator that is t$t$-body and employing the embedding U(N)$U\left(N\right)$ algebra, finite-N$N$ formulas for moments up to order four are derived, for the first time, for the transition strength densities. Secondly, formulas for the moments up to order four are also derived for systems with two types of spinless fermions and a transition operator similar to beta decay and neutrinoless beta decay operators. In addition, moments formulas are also derived for a transition operator that removes k₀$k_0$ number of particles from a system of m$m$ spinless fermions. In the dilute limit, these formulas are shown to reduce to those for the EGOE version derived using the asymptotic limit theory of Mon and French (1975). Numerical results obtained using the exact formulas for two-body (k=2$k=2$) Hamiltonians (in some examples for k=3$k=3$ and 4$4$) and the asymptotic formulas clearly establish that in general the smoothed (with respect to energy) form of the bivariate transition strength densities take bivariate Gaussian form for isolated finite quantum systems. Extensions of these results to bosonic systems and EGUE ensembles with further symmetries are discussed.     

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

MAP,MAC, and vortex-rings configurations in the Weinberg–Salam model

We report on the presence of new axially symmetric monopoles, antimonopoles and vortex-rings solutions of the SU(2)×U(1) Weinberg–Salam model of electromagnetic and weak interactions. When the ?$?$-winding number n=1$n=1$, and 2, the configurations are monopole–antimonopole pair (MAP) and monopole–antimonopole chain (MAC) with poles of alternating sign magnetic charge arranged along the z$z$-axis. Vortex-rings start to appear from the MAP and MAC configurations when the winding number n=3$n=3$. The MAP configurations possess zero net magnetic charge whereas the MAC configurations possess net magnetic charge of 4πn/e$4\pi n/e$.     

Transition between monopole and antimonopole for electrically charged monopole–antimonopole chain system of solutions

The dependence of physical properties of the electrically charged monopole–antimonopole pair (MAP) solutions in the Higgs self-coupling constant is previously investigated. In this paper we study the three-poles monopole–antimonopole chain (MAC) solutions. The study includes ?$?$-winding number n=2,3$n=2,3$, and 4. For the case of n=2$n=2$, no bifurcation and geometrical transition is detected for the interval of the study. For the case of n=3$n=3$, two geometrical transitions happen along the fundamental solution. Also two bifurcations and one joining point is detected for the interval of study. The case of n=4$n=4$ includes one bifurcation. There is also a geometrical transition along the fundamental solution and two transitions along the Higher energy bifurcating branch. This study implies that during some specific kind of geometrical transition, a magnetic and electric charge transition occurs for the pole which is located at the origin.