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.     

Transition between SU(4) and SU(2) Kondo effect

Motivated by experiments in nanoscopic systems, we study a generalized Anderson, which consist of two spin degenerate doublets hybridized to a singlet by the promotion of an electron to two conduction bands, as a function of the energy separation δ$\delta$ between both doublets. For δ=0$\delta =0$ or very large, the model is equivalent to a one-level SU(N$N$) Anderson model, with N=4$N=4$ and 2 respectively. We study the evolution of the spectral density for both doublets (ρ_1σ(ω)${\rho }_{1\sigma }(\omega )$ and ρ_2σ(ω)${\rho }_{2\sigma }(\omega )$) and their width in the Kondo limit as δ$\delta$ is varied, using the non-crossing approximation (NCA). As δ$\delta$ increases, the peak at the Fermi energy in the spectral density (Kondo peak) splits and the density of the doublet of higher energy ρ_2σ(ω)${\rho }_{2\sigma }(\omega )$ shifts above the Ferrmi energy. The Kondo temperature T_K (determined by the half-width at half maximum of the Kondo peak in density of the doublet of lower energy ρ_1σ(ω)${\rho }_{1\sigma }(\omega )$) decreases dramatically. The variation of T_K with δ$\delta$ is reproduced by a simple variational calculation.     

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.     

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.     

Dynamical properties of a dissipative discontinuous map: A scaling investigation

The effects of dissipation on the scaling properties of nonlinear discontinuous maps are investigated by analyzing the behavior of the average squared action 〈I²〉$〈I^2〉$ as a function of the n-th iteration of the map as well as the parameters K and γ  , controlling nonlinearity and dissipation, respectively. We concentrate our efforts to study the case where the nonlinearity is large; i.e., K?1$K?1$. In this regime and for large initial action _I0?K$I_0?K$, we prove that dissipation produces an exponential decay for the average action 〈I〉$〈I〉$. Also, for _I0≅0$I_0\cong 0$, we describe the behavior of 〈I²〉$〈I^2〉$ using a scaling function and analytically obtain critical exponents which are used to overlap different curves of 〈I²〉$〈I^2〉$ onto a universal plot. We complete our study with the analysis of the scaling properties of the deviation around the average action ω.     

Rational r-matrices,higher rank Lie algebras and integrable proton–neutron BCS models

We study integrable cases of pairing BCS hamiltonians containing several types of fermions. We prove that there exist three classes of such integrable models associated with classical rational r  -matrices and Lie algebras gl(2m)$\mathit{gl}(2m)$, sp(2m)$\mathit{sp}(2m)$ and so(2m)$\mathit{so}(2m)$ correspondingly. We diagonalize the constructed hamiltonians by means of the algebraic Bethe ansatz. In the partial case of two types of fermions (m=2$m=2$) the obtained models may be interpreted as N=Z$N=Z$ proton–neutron integrable models. In particular, in the case of sp(4)$\mathit{sp}(4)$ we recover the famous integrable proton–neutron model of Richardson.     

Local derivative post-processing for the discontinuous Galerkin method

Obtaining accurate approximations for derivatives is important for many scientific applications in such areas as fluid mechanics and chemistry as well as in visualization applications. In this paper we discuss techniques for computing accurate approximations of high-order derivatives for discontinuous Galerkin solutions to hyperbolic equations related to these areas. In previous work, improvement in the accuracy of the numerical solution using discontinuous Galerkin methods was obtained through post-processing by convolution with a suitably defined kernel. This post-processing technique was able to improve the order of accuracy of the approximation to the solution of time-dependent symmetric linear hyperbolic partial differential equations from order k+1$k+1$ to order 2k+1$2k+1$ over a uniform mesh; this was extended to include one-sided post-processing as well as post-processing over non-uniform meshes. In this paper, we address the issue of improving the accuracy of approximations to derivatives of the solution by using the method introduced by Thomée [19]. It consists in simply taking the α$\alpha$th-derivative of the convolution of the solution with a sufficiently smooth kernel. The order of convergence of the approximation is then independent   of the order of the derivative, |α|$|\alpha |$. We also discuss an efficient way of computing the approximation which does not involve differentiation but the application of simple finite differencing. Our results show that the above-mentioned approximations to the α$\alpha$th-derivative of the exact solution of linear, multidimensional symmetric hyperbolic systems obtained by the discontinuous Galerkin method with polynomials of degree k$k$ converge with order 2k+1$2k+1$ regardless of the order |α|$|\alpha |$ of the derivative.     

Semiempirical fine-tuning for Hartree–Fock ionization potentials of atomic ions with non-integral atomic number

Amovilli and March (2006) [8] used diffusion quantum Monte Carlo techniques to calculate the non-relativistic ionization potential I(Z)$I(Z)$ in He-like atomic ions for the range of (fractional) nuclear charges Z   lying between the known critical value _Zc=0.911$Z_{\mathrm{c}}=0.911$ at which I(Z)$I(Z)$ tends to zero and Z=2$Z=2$. They showed that it is possible to fit I(Z)$I(Z)$ to a simple quadratic expression. Following that idea, we present here a semiempirical fine-tuning of Hartree–Fock ionization potentials for the isoelectronic series of He, Be, Ne, Mg and Ar-like atomic ions that leads to excellent estimations of _Zc$Z_{\mathrm{c}}$ for these series. The empirical information involved is experimental ionization and electron affinity data. It is clearly demonstrated that Hartree–Fock theory provides an excellent starting point for determining I(Z)$I(Z)$ for these series.     

On DBI Textures with generalized Hopf fibration

In this Letter we show numerical existence of O(4)$O(4)$ Dirac–Born–Infeld (DBI) Textures living in (N+1)$(N+1)$ dimensional spacetime. These defects are characterized by S^N→S³$S^N\to S^3$ mapping, generalizing the well-known Hopf fibration into _πN(S³)${\pi }_N(S^3)$, for all N>3$N>3$. The nonlinear nature of DBI kinetic term provides stability against size perturbation and thus renders the defects having natural scale.     

The puzzle of apparent linear lattice artifacts in the 2d non-linear σ-model and Symanzik's solution

Lattice artifacts in the 2d O(n) non-linear σ  -model are expected to be of the form O(a²)$\mathrm{O}(a^2)$, and hence it was (when first observed) disturbing that some quantities in the O(3)$\mathrm{O}(3)$ model with various actions show parametrically stronger cutoff dependence, apparently O(a)$\mathrm{O}(a)$, up to very large correlation lengths. In a previous letter Balog et al. (2009) [1] we described the solution to this puzzle. Based on the conventional framework of Symanzik's effective action, we showed that there are logarithmic corrections to the O(a²)$\mathrm{O}(a^2)$ artifacts which are especially large (ln³a${\ln }^{3}a$) for n=3$n=3$ and that such artifacts are consistent with the data. In this paper we supply the technical details of this computation. Results of Monte Carlo simulations using various lattice actions for O(3)$\mathrm{O}(3)$ and O(4)$\mathrm{O}(4)$ are also presented.     

Strange attractor in the Potts spin glass on hierarchical lattices

The spin-glass q-state Potts model on d  -dimensional diamond hierarchical lattices is investigated by an exact real space renormalization group scheme. Above a critical dimension _dl(q)$d_l(q)$ for q>2$q>2$, the coupling constants probability distribution flows to a low-temperature strange attractor   or to the high-temperature paramagnetic fixed point, according to the temperature is below or above the critical temperature _Tc(q,d)$T_c(q,d)$. The strange attractor was investigated considering four initial different distributions for q=3$q=3$ and d=5$d=5$ presenting strong robustness in shape and temperature interval suggesting a condensed phase with algebraic decay.     

The fate of conformal symmetry in the non-linear Schrödinger theory

The free Schrödinger theory in d   space dimensions is a non-relativistic conformal field theory. The interacting non-linear theory preserves this symmetry in specific numbers of dimensions at the classical (tree) level. This holds in particular for the |Φ|⁴${|\Phi |}^4$-theory in d=2$d=2$. We compute the full quantum corrections to the 1PI 4-point function in d=2−?$d=2-\mathrm{?}$ dimensions and find a non-trivial β  -function completely given by the 1-loop result. We exhibit an explicit Ward-identity showing that scale-invariance is broken in the limit d=2$d=2$ by an anomalous contribution proportional to the β-function.     

Long-range memory elementary 1D cellular automata: Dynamics and nonextensivity

We numerically study the dynamics of elementary 1D cellular automata (CA), where the binary state _σi(t)∈{0,1}${\sigma }_i(t)\in \{0,1\}$ of a cell i   does not only depend on the states in its local neighborhood at time t-1$t-1$, but also on the memory of its own past states _σi(t-2),_σi(t-3),…,_σi(t-τ),…${\sigma }_i(t-2),{\sigma }_i(t-3),\dots ,{\sigma }_i(t-\tau ),\dots$ . We assume that the weight of this memory decays proportionally to τ^-α${\tau }^{-\alpha }$, with α?0$\alpha ?0$ (the limit α→∞$\alpha \to \infty$ corresponds to the usual CA). Since the memory function is summable for α>1$\alpha >1$ and nonsummable for 0?α?1$0?\alpha ?1$, we expect pronounced changes of the dynamical behavior near α=1$\alpha =1$. This is precisely what our simulations exhibit, particularly for the time evolution of the Hamming distance H   of initially close trajectories. We typically expect the asymptotic behavior H(t)∝t^1/(1-q)$H(t)\propto t^{1/(1-q)}$, where q   is the entropic index associated with nonextensive statistical mechanics. In all cases, the function q(α)$q(\alpha )$ exhibits a sensible change at α?1$\alpha ?1$. We focus on the class II rules 61, 99 and 111. For rule 61, q=0$q=0$ for 0?α?_αc?1.3$0?\alpha ?{\alpha }_c?1.3$, and q<0$q<0$ for α>_αc$\alpha >{\alpha }_c$, whereas the opposite behavior is found for rule 111. For rule 99, the effect of the long-range memory on the spread of damage is quite dramatic. These facts point at a rich dynamics intimately linked to the interplay of local lookup rules and the range of the memory. Finite size scaling studies varying system size N   indicate that the range of the power-law regime for H(t)$H(t)$ typically diverges ∝N^z$\propto N^z$ with 0?z?1$0?z?1$.