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.     

Spin-1 Blume–Capel model with random crystal field effects

The random-crystal field spin-1 Blume–Capel model is investigated by the lowest approximation of the cluster-variation method which is identical to the mean-field approximation. The crystal field is either turned on randomly with probability p$p$ or turned off with q=1−p$q=1-p$ in a bimodal distribution. Then the phase diagrams are constructed on the crystal field (Δ$\Delta$)–temperature (kT/J)$(kT/J)$ planes for given values of p$p$ and on the (kT/J,p$kT/J,p$) planes for given Δ$\Delta$ by studying the thermal variations of the order parameters. In the latter, we only present the second-order phase transition lines, because of the existence of irregular wiggly phase transitions which are not good enough to construct lines. In addition to these phase transitions, the model also yields tricritical points for all values of p$p$ and the reentrant behavior at lower p$p$ values.     

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.     

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.     

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.     

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.     

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

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.     

On the nuclear modification factor at RHIC and LHC

We show that pQCD factorization incorporated with pre-hadronization energy-loss effect naturally leads to flatness of the nuclear modification factor RAA$R_{\mathit{AA}}$ for produced hadrons at high transverse momentum pT$p_T$. We consider two possible scenarios for the pre-hadronization: In scenario 1, the produced gluon propagates through dense QCD medium and loses energy. In scenario 2, all gluons first decay to quark–antiquark pairs and then each pair loses energy as propagating through the medium. We show that the estimates of the energy-loss in these two different models lead to very close values and is able to explain the suppression of high-pT$p_T$ hadrons in nucleus–nucleus collisions at RHIC. We show that the onset of the flatness of RAA$R_{\mathit{AA}}$ for the produced hadron in central collisions at midrapidity is about pT≈15$p_T\approx 15$ and 25 GeV at RHIC and the LHC energies, respectively. We show that the smallness (RAA<0.5$R_{\mathit{AA}}<0.5$ ) and the high-pT$p_T$ flatness of RAA$R_{\mathit{AA}}$ obtained from the kT$k_T$ factorization supplemented with the Balitsky–Kovchegov (BK) equation is rather generic and it does not strongly depend on the details of the BK solutions. We show that energy-loss effect reduces the nuclear modification factor obtained from the kT$k_T$ factorization about 30–50% at moderate pT$p_T$.     

Ranking the spreading influence in complex networks

Identifying the node spreading influence in networks is an important task to optimally use the network structure and ensure the more efficient spreading in information. In this paper, by taking into account the shortest distance between a target node and the node set with the highest k$k$-core value, we present an improved method to generate the ranking list to evaluate the node spreading influence. Comparing with the epidemic process results for four real networks and the Barabási–Albert network, the parameterless method could identify the node spreading influence more accurately than the ones generated by the degree k$k$, closeness centrality, k$k$-shell and mixed degree decomposition methods. This work would be helpful for deeply understanding the node importance of a network.     

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

Target mass corrections for spin-dependent structure functions in collinear factorization

We derive target mass corrections (TMC) for the spin-dependent nucleon structure function g₁$g_1$ and polarization asymmetry A₁$A_1$ in collinear factorization at leading twist. The TMCs are found to be significant for g₁$g_1$ at large xB$x_B$, even at relatively high Q²$Q^2$ values, but largely cancel in A₁$A_1$. A comparison of TMCs obtained from collinear factorization and from the operator product expansion shows that at low Q²$Q^2$ the corrections drive the proton A₁$A_1$ in opposite directions.     

On the structural properties of small-world networks with range-limited shortcut links

We explore a new variant of Small-World Networks (SWNs), in which an additional parameter (r$r$) sets the length scale over which shortcuts are uniformly distributed. When r=0$r=0$ we have an ordered network, whereas r=1$r=1$ corresponds to the original Watts–Strogatz SWN model. These limited range SWNs have a similar degree distribution and scaling properties as the original SWN model. We observe the small-world phenomenon for r?1$r?1$, indicating that global shortcuts are not necessary for the small-world effect. For limited range SWNs, the average path length changes nonmonotonically with system size, whereas for the original SWN model it increases monotonically. We propose an expression for the average path length for limited range SWNs based on numerical simulations and analytical approximations.     

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

Intertwining operators for Sklyanin algebra and elliptic hypergeometric series

Intertwining operators for infinite-dimensional representations of the Sklyanin algebra with spins ?$?$ and −?−1$-?-1$ are constructed using the technique of intertwining vectors for elliptic L$L$-operator. They are expressed in terms of elliptic hypergeometric series with operator argument. The intertwining operators obtained (W$W$-operators) serve as building blocks for the elliptic R$R$-matrix which intertwines tensor product of two L$L$-operators taken in infinite-dimensional representations of the Sklyanin algebra with arbitrary spin. The Yang–Baxter equation for this R$R$-matrix follows from simpler equations of the star–triangle type for the W$W$-operators. A natural graphic representation of the objects and equations involved in the construction is used.     

Particle diagrams and statistics of many-body random potentials

We present a method using Feynman-like diagrams to calculate the statistical properties of random many-body potentials. This method provides a promising alternative to existing techniques typically applied to this class of problems, such as the method of supersymmetry and the eigenvector expansion technique pioneered in Benet et al. (2001). We use it here to calculate the fourth, sixth and eighth moments of the average level density for systems with m$m$ bosons or fermions that interact through a random k$k$-body Hermitian potential (k≤m$k\le m$); the ensemble of such potentials with a Gaussian weight is known as the embedded Gaussian Unitary Ensemble   (eGUE) (Mon and French, 1975). Our results apply in the limit where the number l$l$ of available single-particle states is taken to infinity. A key advantage of the method is that it provides an efficient way to identify only those expressions which will stay relevant in this limit. It also provides a general argument for why these terms have to be the same for bosons and fermions. The moments are obtained as sums over ratios of binomial expressions, with a transition from moments associated to a semi-circular level density for m<2k$m<2k$ to Gaussian moments in the dilute limit k?m?l$k?m?l$. Regarding the form of this transition, we see that as m$m$ is increased, more and more diagrams become relevant, with new contributions starting from each of the points m=2k,3k,…,nk$m=2k,3k,\dots ,nk$ for the 2n$2n$th moment.