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

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.     

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

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.     

Phonon contribution to the shear viscosity of a superfluid Fermi gas in the unitarity limit

We present a detailed analysis of the contribution of small-angle Nambu–Goldstone boson (phonon) collisions to the shear viscosity, η$\eta$, in a superfluid atomic Fermi gas close to the unitarity limit. We show that the experimental values of the shear viscosity coefficient to entropy ratio, η/s$\eta /s$, obtained at the lowest reached temperature can be reproduced assuming that phonons give the leading contribution to η$\eta$. The phonon contribution is evaluated considering 1↔2$1\leftrightarrow 2$ processes and taking into account the finite size of the experimental system. In particular, for very low temperatures, T?0.1T_F$T?0.1T_F$, we find that phonons are ballistic and the contribution of phonons to the shear viscosity is determined by the processes that take place at the interface between the superfluid and the normal phase. This result is independent of the detailed form of the phonon dispersion law and leads to two testable predictions: the shear viscosity should correlate with the size of the optical trap and it should decrease with decreasing temperature. For higher temperatures the detailed form of the phonon dispersion law becomes relevant and, within our model, we find that the experimental data for η/s$\eta /s$ can be reproduced assuming that phonons have an anomalous dispersion law.     

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.     

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

A modified Vicsek model for self-propelled agents with exponential neighbor weight and restricted visual field

Large flocks of wild beings can have coordinated motion with neither leading center nor global information. The Vicsek model and its new versions explained such collective behavior of self-propelled agents with the mechanism of un-weighted influence of neighbors and global visual field (GVF), which may not always fit the reality. In the present work, we take both exponential (λ$\lambda$) neighbor weight (ENW) and restricted visual field (RVF) into account, and investigate the combinatory effect of them. Based on this mechanism, in most cases of our simulation, stationary direction consensus _Va$V_a$ exhibits a maximum or decreases with VF$VF$, and a minimum convergence time _Tmin$T_{min}$ appears in RVF, which demonstrates that GVF is not necessary for coordinated motion, while RVF matched ENW would be reasonable. Moreover, we simulated various recipes for optimizing _Va$V_a$ and T$T$ in unfavorable cases, and found critical points separating inverse λ$\lambda$-dependent variations of T$T$. The present work may be useful in improving the efficiency of direction consensus of large flocks of artificial individuals besides the understanding of biological collective motion in nature.     

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.     

Damage spreading in a driven lattice gas model

We studied damage spreading in a Driven Lattice Gas (DLG) model as a function of the temperature T$T$, the magnitude of the external driving field E$E$, and the lattice size. The DLG model undergoes an order–disorder second-order phase transition at the critical temperature _Tc(E)$T_c\left(E\right)$, such that the ordered phase is characterized by high-density strips running along the direction of the applied field; while in the disordered phase one has a lattice-gas-like behavior. It is found that the damage always spreads for all the investigated temperatures and reaches a saturation value _Dsat$D_{sat}$ that depends only on T$T$. _Dsat$D_{sat}$ increases for T<_Tc(E=∞)$T<T_c(E=\infty )$, decreases for T>_Tc(E=∞)$T>T_c(E=\infty )$ and is free of finite-size effects. This behavior can be explained as due to the existence of interfaces between the high-density strips and the lattice-gas-like phase whose roughness depends on T$T$. Also, we investigated damage spreading for a range of finite fields as a function of T$T$, finding a behavior similar to that of the case with E=∞$E=\infty$.     

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.     

Entropy and stability of phase synchronisation of oscillators on networks

I examine the role of entropy in the transition from incoherence to phase synchronisation in the Kuramoto model of N$N$ coupled phase oscillators on a general undirected network. In a Hamiltonian ‘action-angle’ formulation, auxiliary variables J_i$J_i$ combine with the phases θ_i${\theta }_i$ to determine a conserved system with a 2N$2N$ dimensional phase space. In the vicinity of the fixed point for phase synchronisation, θ_i≈θ_j${\theta }_i\approx {\theta }_j$, which is known to be stable, the auxiliary variables J_i$J_i$ exhibit instability  . This manifests Liouville’s Theorem in the phase synchronised regime in that contraction in the θ_i${\theta }_i$ parts of phase space are compensated for by expansion in the auxiliary dimensions. I formulate an entropy rate based on the projection of the J_i$J_i$ onto eigenvectors of the graph Laplacian that satisfies Pesin’s Theorem. This leads to the insight that the evolution to phase synchronisation of the Kuramoto model is equivalent to the approach to a state of monotonically increasing entropy. Indeed, for unequal intrinsic frequencies on the nodes, the networks that achieve the closest to exact phase synchronisation are those which enjoy the highest entropy production. I compare numerical results for a range of networks.     

Rigidity of noncompact complete Bach-flat manifolds

Let (M,g)$(M,g)$ be a noncompact complete Bach-flat manifold with positive Yamabe constant. We prove that (M,g)$(M,g)$ is flat if (M,g)$(M,g)$ has zero scalar curvature and sufficiently small _L2$L_2$ bound of curvature tensor. When (M,g)$(M,g)$ has nonconstant scalar curvature, we prove that (M,g)$(M,g)$ is conformal to the flat space if (M,g)$(M,g)$ has sufficiently small _L2$L_2$ bound of curvature tensor and _L4/3$L_{4/3}$ bound of scalar curvature.     

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.     

Entanglement thresholds for displaying the quantum nature of teleportation

A protocol for transferring an unknown single qubit state evidences quantum features when the average fidelity of the outcomes is, in principle, greater than 2/3$2/3$. We propose to use the probabilistic and unambiguous state extraction scheme   as a mechanism to redistribute the fidelity in the outcome of the standard teleportation when the process is performed with an X$X$-state as a noisy quantum channel. We show that the entanglement of the channel is necessary but not sufficient in order for the average fidelity f_X$f_X$ to display quantum features, i.e., we find a threshold C_X$C_X$ for the concurrence of the channel. On the other hand, if the mechanism for redistributing fidelity is successful then we find a filterable outcome with average fidelity f_X,0$f_{X,0}$ that can be greater than f_X$f_X$. In addition, we find the threshold concurrence of the channel C_X,0$C_{X,0}$ in order for the average fidelity f_X,0$f_{X,0}$ to display quantum features and surprisingly, the threshold concurrence C_X,0$C_{X,0}$ can be less than C_X$C_X$. Even more, we find some special cases for which the threshold values become zero.     

Coherent superpositions of states in coupled Hilbert-space using step by step Morris–Shore transformation

Creation of coherent superpositions in quantum systems with N_a$N_a$ states in the lower set and N_b$N_b$ states in the upper set is presented. The solution is drived by using the Morris–Shore transformation, which step by step reduces the fully coupled system to a three-state Λ$\Lambda$-like system and a set of decoupled states. It is shown that, for properly timed pulse, robust population transfer from an initial ground state (or superposition of M$M$ ground states) to an arbitrary coherent superposition of the ground states can be achieved by coincident pulses and/or STIRAP techniques.