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.     

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.     

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.     

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

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.     

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.     

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.     

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.     

Bertrand curves in the three-dimensional sphere

A curve α$\alpha$ immersed in the three-dimensional sphere S³${\mathbb{S}}^3$ is said to be a Bertrand curve if there exists another curve β$\beta$ and a one-to-one correspondence between α$\alpha$ and β$\beta$ such that both curves have common principal normal geodesics at corresponding points. The curves α$\alpha$ and β$\beta$ are said to be a pair of Bertrand curves in S³${\mathbb{S}}^3$. One of our main results is a sort of theorem for Bertrand curves in S³${\mathbb{S}}^3$ which formally agrees with the classical one: “Bertrand curves in S³${\mathbb{S}}^3$ correspond to curves for which there exist two constants λ≠0$\lambda \ne 0$ and μ$\mu$ such that λκ+μτ=1$\lambda \kappa +\mu \tau =1$”, where κ$\kappa$ and τ$\tau$ stand for the curvature and torsion of the curve; in particular, general helices in the 3-sphere introduced by M. Barros are Bertrand curves. As an easy application of the main theorem, we characterize helices in S³${\mathbb{S}}^3$ as the only twisted curves in S³${\mathbb{S}}^3$ having infinite Bertrand conjugate curves. We also find several relationships between Bertrand curves in S³${\mathbb{S}}^3$ and (1,3)-Bertrand curves in R⁴${\mathbb{R}}^4$.     

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

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

Entanglement entropy and variational methods: Interacting scalar fields

We develop a variational approximation to the entanglement entropy for scalar ?⁴${?}^4$ theory in 1+1$1+1$, 2+1$2+1$, and 3+1$3+1$ dimensions, and then examine the entanglement entropy as a function of the coupling. We find that in 1+1$1+1$ and 2+1$2+1$ dimensions, the entanglement entropy of ?⁴${?}^4$ theory as a function of coupling is monotonically decreasing and convex. While ?⁴${?}^4$ theory with positive bare coupling in 3+1$3+1$ dimensions is thought to lead to a trivial free theory, we analyze a version of ?⁴${?}^4$ with infinitesimal negative bare coupling, an asymptotically free theory known as precarious  ?⁴${?}^4$ theory, and explore the monotonicity and convexity of its entanglement entropy as a function of coupling. Within the variational approximation, the stability of precarious ?⁴${?}^4$ theory is related to the sign of the first and second derivatives of the entanglement entropy with respect to the coupling.     

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 dynamics in the distribution of players’ scoring performance in Rugby Union world cups

The evolution of the scoring performance of Rugby Union players is investigated over the seven rugby world cups (RWC) that took place from 1987 to 2011, and a specific attention is given to how they may have been impacted by the switch from amateurism to professionalism that occurred in 1995. The distribution of the points scored by individual players, _Ps$P_s$, ranked in order of performance were well described by the simplified canonical law _Ps∝(r+?)^−α$P_s\propto {(r+?)}^{-\alpha }$, where r$r$ is the rank, and ?$?$ and α$\alpha$ are the parameters of the distribution. The parameter α$\alpha$ did not significantly change from 1987 to 2007 (α=0.92±0.03)$(\alpha =0.92\pm 0.03)$, indicating a negligible effect of professionalism on players’ scoring performance. In contrast, the parameter ?$?$ significantly increased from ?=1.32$?=1.32$ for 1987 RWC, ?=2.30$?=2.30$ for 1999 to 2003 RWC and ?=5.60$?=5.60$ for 2007 RWC, suggesting a progressive decrease in the relative performance of the best players. Finally, the sharp decreases observed in both α(α=0.38)$\alpha (\alpha =0.38)$ and ?(?=0.70)$?(?=0.70)$ in the 2011 RWC indicate a more even distribution of the performance of individuals among scorers, compared to the more heterogeneous distributions observed from 1987 to 2007, and suggest a sharp increase in the level of competition leading to an increase in the average quality of players and a decrease in the relative skills of the top players. Note that neither α$\alpha$ nor ?$?$ significantly correlate with traditional performance indicators such as the number of points scored by the best players, the number of games played by the best players, the number of points scored by the team of the best players or the total number of points scored over each RWC. This indicates that the dynamics of the scoring performance of Rugby Union players is influenced by hidden processes hitherto inaccessible through standard performance metrics; this suggests that players’ scoring performance is connected to ubiquitous phenomena such as anomalous diffusion.     

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.     

Minimal hypersurfaces in a nearly Kaehler 6-sphere

In this paper we show that for a compact minimal hypersurface M$M$ of constant scalar curvature in the unit sphere S⁶$S^6$ with the shape operator A$A$ satisfying ‖A‖²>5${\Vert A\Vert }^2>5$, there exists an eigenvalue λ>10$\lambda >10$ of the Laplace operator of the hypersurface M$M$ such that ‖A‖²=λ−5${\Vert A\Vert }^2=\lambda -5$. This gives the next discrete value of ‖A‖²${\Vert A\Vert }^2$ greater than 0 and 5.     

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