Heavy neutral gauge bosons at the LHC in an extended MSSM

Searching for heavy neutral gauge bosons Z^′$Z^{\prime }$, predicted in extensions of the Standard Model based on a U(1)^′$\mathrm{U}{(1)}^{\prime }$ gauge symmetry, is one of the challenging objectives of the experiments carried out at the Large Hadron Collider. In this paper, we study Z^′$Z^{\prime }$ phenomenology at hadron colliders according to several U(1)^′$\mathrm{U}{(1)}^{\prime }$-based models and in the Sequential Standard Model. In particular, possible Z^′$Z^{\prime }$ decays into supersymmetric particles are included, in addition to the Standard Model modes so far investigated. We point out the impact of the U(1)^′$\mathrm{U}{(1)}^{\prime }$ group on the MSSM spectrum and, for a better understanding, we consider a few benchmarks points in the parameter space. We account for the D-term contribution, due to the breaking of U(1)^′$\mathrm{U}{(1)}^{\prime }$, to slepton and squark masses and investigate its effect on Z^′$Z^{\prime }$ decays into sfermions. Results on branching ratios and cross sections are presented, as a function of the MSSM and U(1)^′$\mathrm{U}{(1)}^{\prime }$ parameters, which are varied within suitable ranges. We pay special attention to final states with leptons and missing energy and make predictions on the number of events with sparticle production in Z^′$Z^{\prime }$ decays, for a few values of integrated luminosity and centre-of-mass energy of the LHC.     

Can density functional theory make use of experimentally determined ground-state electron densities via the Dirac density matrix for molecules of biological interest?

More than four decades ago, March and Murray gave a perturbation theory of the single-particle(s) Dirac density matrix γs(r,r^′)${\gamma }_s(r,r^{\prime })$ to all orders in a given one-body potential energy V(r)$V(r)$. However, for density functional theory in orbital-free form, one requires the functional γs[ρ]${\gamma }_s[\rho ]$ where ρ(r)$\rho (r)$ is the ground-state electron density. Therefore, in the present study, a first-order non-linear differential equation is proposed for γs${\gamma }_s$ in terms of ρ(r)$\rho (r)$ and ∇ρ(r)$\mathrm{\nabla }\rho (r)$, plus the single-particle kinetic energy. Since this latter quantity is itself known to be a functional of ρ  , the existence of such an equation for γs${\gamma }_s$ would be a significant step along the road to determining the desired functional γs[ρ]${\gamma }_s[\rho ]$. As yet, we have succeeded in giving a rigorous proof of the proposed differential equation for γs(r,r^′)${\gamma }_s(r,r^{\prime })$ only for one- and two-level molecules. If it is subsequently proved for an arbitrary number of levels, which we believe should be possible, it would then allow γs${\gamma }_s$ to be calculated for molecules of biological interest, from experimentally measured ground-state densities ρ(r)$\rho (r)$, as the approach is entirely orbital-free.     

A note on embedding non-Abelian finite flavor groups in continuous groups

A non-Abelian finite flavor group G⊂SO(3)$G\subset \mathit{SO}(3)$ can have double covering G^′⊂SU(2)$G^{\prime }\subset \mathit{SU}(2)$ such that G⊄G^′$G\not\subset G^{\prime }$. This situation is not contradictory, but quite natural, and we give explicit examples such as G=Dn$G=D_n$, G^′=Q_2n$G^{\prime }=Q_{2n}$ and G=T$G=T$, G^′=T^′$G^{\prime }=T^{\prime }$. This observation can be crucial in particle theory model building.     

Hadro-charmonium

We argue that relatively compact charmonium states, J/ψ$J/\psi$, ψ(2S)$\psi (2S)$, χc${\chi }_c$, can very likely be bound inside light hadronic matter, in particular inside higher resonances made from light quarks and/or gluons. The charmonium state in such binding essentially retains its properties, so that the bound system decays into light mesons and the particular charmonium resonance. Thus such bound states of a new type, which we call hadro-charmonium, may explain the properties of some of the recently observed resonant peaks, in particular of Y(4.26)$Y(4.26)$, Y(4.32–4.36)$Y(4.32\text{–}4.36)$, Y(4.66)$Y(4.66)$, and Z(4.43)$Z(4.43)$. We discuss further possible implications of the suggested picture for the observed states and existence of other states of hadro-charmonium and hadro-bottomonium.     

Single eta production in heavy quarkonia: Breakdown of multipole expansion

The η   production in the (n,n^′)$(n,n^{\prime })$ bottomonium transitions ?(n)→?(n^′)η$?(n)\to ?(n^{\prime })\eta$, is studied in the method used before for dipion heavy quarkonia transitions. The widths Γη(n,n^′)${\Gamma }_{\eta }(n,n^{\prime })$ are calculated without fitting parameters for n=2,3,4,5$n=2,3,4,5$, n^′=1$n^{\prime }=1$. Resulting Γη(4,1)${\Gamma }_{\eta }(4,1)$ is found to be large in agreement with recent data. Multipole expansion method is shown to be inadequate for large size systems considered.     

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.     

Pion radiative weak decays in nonlocal chiral quark models

We analyze the radiative pion decay π⁺→e⁺νeγ${\pi }^{+}\to e^{+}{\nu }_e\gamma$ within nonlocal chiral quark models that include wave function renormalization. In this framework we calculate the vector and axial-vector form factors FV$F_V$ and FA$F_A$ at q²=0$q^2=0$ — where q²$q^2$ is the e⁺νe$e^{+}{\nu }_e$ squared invariant mass — and the slope a   of FV(q²)$F_V(q^2)$ at q²→0$q^2\to 0$. The calculations are carried out considering different nonlocal form factors, in particular those taken from lattice QCD evaluations, showing a reasonable agreement with the corresponding experimental data. The comparison of our results with those obtained in the (local) NJL model and the relation of FV$F_V$ and a   with the form factor in π⁰→γ^?γ${\pi }^0\to {\gamma }^{?}\gamma$ decays 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.     

Atomic parity violation in the economical 3–3–1 model

The deviation δQW$\delta Q_{\mathrm{W}}$ of the weak charge from its standard model prediction due to the mixing of the W boson with the charged bilepton Y as well as of the Z   boson with the neutral Z^′$Z^{\prime }$ and the real part of the non-Hermitian neutral bilepton X   in the economical 3–3–1 model is established. Additional contributions to the usual δQW$\delta Q_{\mathrm{W}}$ expression in the extra U(1)$\mathrm{U}(1)$ models and the left–right models are obtained. Our calculations are quite different from previous analyzes in this kind of the 3–3–1 models and give the limit on mass of the Z^′$Z^{\prime }$ boson, the Z–Z^′$Z\text{–}{Z}^{\prime }$ and W–Y$W\text{–}Y$ mixing angles with the more appropriate values: MZ^′>564 GeV$M_{Z^{\prime }}>564\text{GeV}$, −0.018<sinφ<0$-0.018<\sin \phi <0$ and |sinθ|<0.043$|\sin \theta |<0.043$.     

Two-dimensional gauge theory and matrix model

We study a matrix model obtained by dimensionally reducing Chern–Simons theory on S³$S^3$. We find that the matrix integration is decomposed into sectors classified by the representation of SU(2)$\mathit{SU}(2)$. We show that the N  -block sectors reproduce SU(N)$\mathit{SU}(N)$ Yang–Mills theory on S²$S^2$ as the matrix size goes to infinity.     

Three ways to solve the Poisson equation on a sphere with Gaussian forcing

Motivated by the needs of vortex methods, we describe three different exact or approximate solutions to the Poisson equation on the surface of a sphere when the forcing is a Gaussian of the three-dimensional distance, ∇²ψ=exp(-2?²(1-cos(θ))-C^Gauss(?)${\nabla }^2\psi =\exp (-2{?}^2(1-\cos (\theta ))-C^{\mathit{Gauss}}(?)$. (More precisely, the forcing is a Gaussian minus the “Gauss constraint constant”, C^Gauss$C^{\mathit{Gauss}}$; this subtraction is necessary because ψ$\psi$ is bounded, for any type of forcing, only if the integral of the forcing over the sphere is zero [Y. Kimura, H. Okamoto, Vortex on a sphere, J. Phys. Soc. Jpn. 56 (1987) 4203–4206; D.G. Dritschel, Contour dynamics/surgery on the sphere, J. Comput. Phys. 79 (1988) 477–483]. The Legendre polynomial series is simple and yields the exact value of the Gauss constraint constant, but converges slowly for large ?$?$. The analytic solution involves nothing more exotic than the exponential integral, but all four terms are singular at one or the other pole, cancelling in pairs so that ψ$\psi$ is everywhere nice. The method of matched asymptotic expansions yields simpler, uniformly valid approximations as series of inverse even powers of ?$?$ that converge very rapidly for the large values of ?  (?>40)$(?>40)$ appropriate for geophysical vortex computations. The series converges to a nonzero O(exp(-4?²))$O(\exp (-4{?}^2))$ error everywhere except at the south pole where it diverges linearly with order instead of the usual factorial order.     

Magnetic properties of the bond and crystal field dilution spin-3/2 Blume–Capel model in an external magnetic field

Magnetic properties of the bond and crystal field dilution spin-3/2 Blume–Capel model in an external magnetic field (h)$(h)$ on simple cubic lattice are studied by using the effective field theory. In the m−T$m-T$ plane, the degeneracy of the magnetization (m)$(m)$ is affected by the concentration of bond or crystal field dilution at low temperature (T)$(T)$. The magnetization curves can appear to fluctuate in certain regions of negative crystal field. In the m−h$m-h$ plane, the initial magnetization curve has an irregular behavior due to the introduction of bond dilution. The crystal field dilution has the influence on the process of magnetic domain displacement. In the χ−h$\chi -h$ plane, there exists one susceptibility (χ)$(\chi )$ shoulder and one step for different negative crystal field. The susceptibility curve takes on the feature of multi-peaks distribution under bond and crystal field dilution conditions.     

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.     

First-order density matrices in one dimension for independent fermions and impenetrable bosons in harmonic traps

To complement existing knowledge of the density matrix γF(x,y)${\gamma }_{\mathrm{F}}(x,y)$ of independent fermions for N   particles in one dimension under harmonic confinement, the corresponding matrix γIB(x,y)${\gamma }_{\mathrm{IB}}(x,y)$ for impenetrable bosons is given for N=2$N=2$ and 3 (with the N=4$N=4$ form available also). For fermions the momentum density is then obtained and illustrated numerically for N=10$N=10$. The boson momentum density is studied analytically at high momentum p  , the coefficients of the p⁻⁴$p^{\mathrm{-4}}$ and p⁻⁶$p^{\mathrm{-6}}$ terms being tabulated for N=2–5$N=2\text{–}5$ inclusive. Their dependence on powers of N   is exhibited numerically. Finally, the functional relationship between γIB(x,y)${\gamma }_{\mathrm{IB}}(x,y)$ and γF(x,y)${\gamma }_{\mathrm{F}}(x,y)$ is formally set out and illustrated.     

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.     

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.     

Exact height distributions for the KPZ equation with narrow wedge initial condition

We consider the KPZ equation in one space dimension with narrow wedge initial condition, h(x,t=0)=−|x|/δ$h(x,t=0)=-|x|/\delta$, δ?1$\delta ?1$, evolving into a parabolic profile with superimposed fluctuations. Based on previous results for the weakly asymmetric simple exclusion process with step initial conditions, we obtain a determinantal formula for the one-point distribution of the solution h(x,t)$h(x,t)$ valid for any x   and t>0$t>0$. The corresponding distribution function converges in the long time limit, t→∞$t\to \infty$, to the Tracy–Widom distribution. The first order correction is a shift of order t−1/3$t^{-1/3}$. We provide numerical computations based on the exact formula.