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

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.     

Reduction of generalized complex structures

We study reduction of generalized complex structures. More precisely, we investigate the following question. Let J$J$ be a generalized complex structure on a manifold M$M$, which admits an action of a Lie group G$G$ preserving J$J$. Assume that _M0$M_0$ is a G$G$-invariant smooth submanifold and the G$G$-action on _M0$M_0$ is proper and free so that _MG?_M0/G$M_G?M_0/G$ is a smooth manifold. Under what condition does J$J$ descend to a generalized complex structure on _MG$M_G$? We describe a sufficient condition for the reduction to hold, which includes the Marsden–Weinstein reduction of symplectic manifolds and the reduction of the complex structures in Kähler manifolds as special cases. As an application, we study reduction of generalized Kähler manifolds.     

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.     

Effects of colored noise and noise delay on a calcium oscillation system

As a calcium oscillations system is in steady state, the effects of colored noise and noise delay on the system is investigated using stochastic simulation methods. The results indicate that: (1) the colored noise can induce coherence bi-resonance phenomenon. (2) there exist three peaks in the R–_τ0$R\text{–}{\tau }_0$ (R$R$ is the reciprocal coefficient of variance, and _τ0${\tau }_0$ is the self-correlation time of the colored noise) curves. For the same noise intensity Q=1$Q=1$, the Gaussian colored noise can induce calcium spikes but the white noise cannot do this. (3) the delay time can improve noise induced spikes regularity as _τ0${\tau }_0$ is small, and R$R$ has a significant minimum with increasing τ$\tau$ as _τ0${\tau }_0$ is large. (4) large values of ζ$\zeta$ reduce noise induced spikes regularity.     

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.     

Properties of reachable sets in the sub-Lorentzian geometry

The aim of this paper is to develop local theory of future timelike, nonspacelike and null reachable sets from a given point _q0$q_0$ in the sub-Lorentzian geometry. In particular, we prove that if U$U$ is a normal neighbourhood of _q0$q_0$ then the three reachable sets, computed relative to U$U$, have identical interiors and boundaries with respect to U$U$. Further, among other things, we show that for Lorentzian metrics on contact distributions on R²ⁿ⁺¹${\mathbb{R}}^{2n+1}$, n≥1$n\ge 1$, the boundary of reachable sets from _q0$q_0$ is, in a neighbourhood of _q0$q_0$, made up of null future directed curves starting from _q0$q_0$. Every such curve has only a finite number of non-smooth points; smooth pieces of every such curve are Hamiltonian geodesics. For general sub-Lorentzian structures, contrary to the Lorentzian case, timelike curves may appear on the boundary. It turns out that such curves are always Goh curves. We also generalize a classical result on null Lorentzian geodesics: every null future directed Hamiltonian sub-Lorentzian geodesic initiating at _q0$q_0$ is contained, at least to a certain moment of time, in the boundary of the reachable set from _q0$q_0$.     

The first Chevalley–Eilenberg Cohomology group of the Lie algebra on the transverse bundle of a decreasing family of foliations

In [L. Lebtahi, Lie algebra on the transverse bundle of a decreasing family of foliations, J. Geom. Phys. 60 (2010), 122–133], we defined the transverse bundle V^k$V^k$ to a decreasing family of k$k$ foliations _Fi$F_i$ on a manifold M$M$. We have shown that there exists a (1,1)$(1,1)$ tensor J$J$ of V^k$V^k$ such that J^k≠0$J^k\ne 0$, J^k+1=0$J^{k+1}=0$ and we defined by _LJ(V^k)$L_J\left(V^k\right)$ the Lie Algebra of vector fields X$X$ on V^k$V^k$ such that, for each vector field Y$Y$ on V^k$V^k$, [X,JY]=J[X,Y]$[X,JY]=J[X,Y]$.     

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.     

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

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.     

On leafwise conformal diffeomorphisms

For every diffeomorphism φ:M→N$\phi :M\to N$ between 3-dimensional Riemannian manifolds M$M$ and N$N$, there are locally two 2-dimensional distributions _D±$D_{\pm }$ such that φ$\phi$ is conformal on both of them. We state necessary and sufficient conditions for a distribution to be one of _D±$D_{\pm }$. These are algebraic conditions expressed in terms of the self-adjoint and positive definite operator induced from _φ∗${\phi }_{\ast }$. We investigate the integrability condition of _D+$D_{+}$ and _D−$D_{-}$. We also show that it is possible to choose coordinate systems in which leafwise conformal diffeomorphism is holomorphic on leaves.     

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.     

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.