I gently introduce the diagrammatic birdtrack notation, first for vector algebra and then for permutations. After moving on to general tensors I review some recent results on Hermitian Young operators, gluon projectors, and multiplet bases for SU(N) color space.

Search for new physics events at the LHC mostly rely on the assumption that the events are characterized in terms of standard-reconstructed objects such as isolated photons, leptons, and jets initiated by QCD-partons. While such strategy works for a vast majority of physics beyond the standard model scenarios, there are examples aplenty where new physics give rise to anomalous objects (such as collimated and equally energetic particles, decays due to long lived particles etc.) in the detectors, which can not be classified as any of the standard-objects. Varied methods and search strategies have been proposed, each of which is trained and optimized for specific models, topologies, and model parameters. Further, as LHC keeps excluding all expected candidates for new physics, the need for a generic method/tool that is capable of finding the unexpected can not be understated. In this paper, we propose one such method that relies on the philosophy that all anomalous objects are $\it{not}$ standard-objects. The anomaly finder, we suggest, simply is a collection of vetoes that eliminate all standard-objects upto a pre-determined acceptance rate. Any event containing at least one anomalous object (that passes all these vetoes), can be identified as a candidate for new physics. Subsequent offline analyses can determine the nature of the anomalous object as well as of the event, paving a robust way to search for these new physics scenarios in a model-independent fashion. Further, since the method relies on learning only the standard-objects, for which control samples are readily available from data, one can build the analysis in an entirely data-driven way.

Jul 25 2017

hep-ph arXiv:1707.07666v1

Recently LHCb discovered the first doubly-charmed baryon $\Xi_{cc}^{++} = ccu$ at $3621.40 \pm 0.78$ MeV, very close to our theoretical prediction. We use the same methods to predict a doubly-bottom tetraquark $Tq(bb\bar u\bar d)$ with $J^P{=}1^+$ at $10,389\pm 12$ MeV, 215 MeV below the $BB^*$ threshold. $Tq(bb\bar u\bar d)$ is therefore stable under strong and EM interactions and can only decay weakly, the first exotic hadron with such a property. On the other hand, the mass of $Tq(cc\bar u\bar d)$ with $J^P{=}1^+$ is predicted to be $3882\pm12$ MeV, 7 MeV above the $D^0 D^{*+}$ threshold. $Tq(bc\bar u\bar d)$ with $J^P{=}0^+$ is predicted at $7134\pm13$ MeV, 11 MeV below the $\bar B^0 D^0$ threshold. Our precision is not sufficient to determine whether $cc\bar u\bar d$ and $bb\bar u\bar d$ are actually above or below the threshold. They could manifest themselves as narrow resonances just at threshold.

The determination of the neutrino mass ordering is currently pursued by several experiments and proposals. A very challenging one is its evaluation from reactor experiments based on the tiny interference effect between the $\Delta m^2_{31}$ and $\Delta m^2_{32}$ oscillations. Current analyses require several years of data taking and an extreme energy resolution to achieve anyhow less than 5 $\sigma$. Referring to the JUNO experimental conditions we developed a completely new technique that will allow to definitively perform a robust 5 $\sigma$ measurement in less than six years of running. Evaluation and inclusion of systematic errors and backgrounds have been performed, the most relevant among them being the addition of the two remote reactor plants 250 km away. Baselines of each contributing reactor core and its spatial resolution have been taken into account. Possible results after two years of running and the foreseen initially-reduced available reactor power have been studied, too. These results confirm the very positive perspectives for JUNO to determine the mass ordering in a vacuum-oscillation dominated region. As a side result the new technique provides a $\Delta m^2_{32(13)}$ evaluation with an unprecedented accuracy at reactors.

Jul 25 2017

hep-ph arXiv:1707.07606v1

We discuss the structure of rapidity divergences that are presented in the soft factors of transverse momentum dependent (TMD) factorization theorems. To provide the discussion on the most general level we consider soft factors for multi-parton scattering. We show that the rapidity divergences are result of the gluon exchanges with the distant transverse plane, and are structurally equivalent to the ultraviolet divergences. It allows to formulate and to prove the renormalization theorem for rapidity divergences. The proof is made with the help the conformal transformation which maps rapidity divergences to ultraviolet divergences. The theorem is the systematic form of the factorization of rapidity divergences, which is required for the definition of TMD parton distributions. In particular, the definition of multi parton distributions is presented. The equivalence of ultraviolet and rapidity divergences leads to the exact relation between soft and rapidity anomalous dimensions. Using this relation we derive the rapidity anomalous dimension at the three-loop order.

We review how Quark-Hadron Duality (QHD) for (u,d) flavors at high energies and in the scaling regime suggests a radial and angular behaviour of mesonic and baryonic resonance masses of the Regge form $M^2_{nJ} = \mu^2 n + \beta^2 J + M_0^2 $. The radial mass dependence is asymptotically consistent with a common two-body dynamics for mesons and baryons in terms of the quark-antiquark ($q \bar q)$ and quark-diquark ($qD$) degrees of freedom, respectively. This formula is validated phenomenologically within an uncertainty determined by half the width of the resonances, $\Delta M_{nJ}^2 \sim \Gamma_{nJ} M_{nJ}$. With this error prescription we find from the non-strange PDG hadrons different radial slopes $\mu^2_{\bar qq}=1.34(4) {\rm GeV}^2$ and $\mu^2_{qD}=0.75(3) {\rm GeV}^2$, but similar angular slopes $\beta_{\bar q q}^2 \sim \beta_{ q D}^2 \sim 1.15 {\rm GeV}^2$.

Solving the Dirac equation for a positron in the Coulomb field of a nucleus with $Z>Z_{\rm cr}$ we observe (confirm) the resonance dependence of the scattering phase $\delta_\varkappa$ on the energy $\varepsilon_{\rm p}$ of the positron. A physical interpretation of these resonances is suggested.

Jul 25 2017

hep-ph arXiv:1707.07364v1

We propose new mechanisms for small neutrino masses based on clockwork mechanism. The Standard Model neutrinos and lepton number violating operators communicate through the zero mode of clockwork gears, one of the two couplings of the zero mode is exponentially suppressed by clockwork mechanism. Different types of models are realized depending on the profile and chirality of the zero mode fermion. Each type of realization would have phenomenologically distinctive features with the accompanying heavy neutrinos.

We compute modifications to the jet spectrum in the presence of a dense medium. We show that in the large-$N_c$ approximation and at leading logarithmic accuracy the jet nuclear modification factor factorizes into a quenching factor associated to the total jet color charge and a Sudakov suppression factor which accounts for the energy loss of jet substructure fluctuations. This factor, called the jet collimator, implements the fact that subjets, that are not resolved by the medium, lose energy coherently as a single color charge, whereas resolved large angle fluctuations suffer more quenching. For comparison, we show that neglecting color coherence results in a stronger suppression of the jet nuclear modification factor.

We discuss a novel approach to systematically determine the long-distance contribution to $B\to K^*\ell\ell$ decays in the kinematic region where the dilepton invariant mass is below the open charm threshold. This approach provides the most consistent and reliable determination to date and can be used to compute Standard Model predictions for all observables of interest, including the kinematic region where the dilepton invariant mass lies between the $J/\psi$ and the $\psi(2S)$ resonances. We illustrate the power of our results by performing a New Physics fit to the Wilson coefficient $C_9$. This approach is systematically improvable from theoretical and experimental sides, and applies to other decay modes of the type $B\to V\ell\ell$, $B\to P\ell\ell$ and $B\to V\gamma$.

Jul 25 2017

hep-ph arXiv:1707.07300v1

We extend the standard model with three right-handed singlet neutrino and a real singlet scalar. We imposed two $Z_2$ and $Z_2^\prime$-symmetries. We explain the tiny neutrino mass-squared differences with two $Z_2$- and $Z_2^\prime$-even right-handed neutrinos using the type-I seesaw mechanism. The $Z_2$-odd fermion and the $Z_2^\prime$-odd scalar can both serve as viable dark matter candidates. We identify new region of parameter spaces, which are consistent with relic density of the dark matter, recent direct search experiment LUX-2016 and LHC data.

We argue that the large suppression of the $\psi(2S)$ inclusive cross-section relative to the $J/\psi$ inclusive cross-section in proton-nucleus (p+A) collisions can be attributed to factorization breaking effects in the formation of quarkonium. These factorization breaking effects arise from soft color exchanges between charm-anticharm pairs undergoing hadronization and comoving partons that are long-lived on time scales of quarkonium formation. We compute the short distance pair production of heavy quarks in the Color Glass Condensate (CGC) effective field theory and employ an improved Color Evaporation model (ICEM) to describe their hadronization into quarkonium at large distances. The combined CGC+ICEM model provides a quantitative description of $J/\psi$ and $\psi(2S)$ data in proton-proton (p+p) collisions from both RHIC and the LHC. Factorization breaking effects in hadronization, due to additional parton comovers in the nucleus, are introduced heuristically by imposing a cutoff $\Lambda$, representing the momentum kick from soft color exchanges, in the ICEM model. Such soft exchanges have no perceptible effect on $J/\psi$ suppression in p+A collisions. In contrast, the interplay of the physics of these soft exchanges at large distances, with the physics of semi-hard rescattering at short distances, causes a significant additional suppression of $\psi(2S)$ yields relative to that of the $J/\psi$. A good fit of all RHIC and LHC $J/\psi$ and $\psi(2S)$ data, for transverse momenta $P_\perp\leq 5$ GeV in p+p and p+A collisions, is obtained for $\Lambda\sim 10$ MeV.

The elastic scattering of an atomic nucleus plays a central role in the dark matter direct detection experiments. The recoil signatures are detected through ionization, scintillation, and the production of heat in the detectors. In those experiments, it is usually assumed that the atomic electrons around the nucleus immediately follow the motion of the recoiled nucleus. However, it takes some time for the electrons to catch up, which results in ionizations and excitations of the atom. In previous studies, the final state ionizations/excitations are treated separately from the nuclear recoil and are estimated by using the so-called Migdal's approach. In this paper, we reformulate the Migdal's approach so that the "atomic recoil" cross section is obtained coherently, where we make transparent the energy-momentum conservation and the probability conservation. We show that the final state ionizations/excitations can enhance the detectability of a rather light dark matter in the GeV mass range even through the \it nuclear scattering. We also discuss the neutrino-nucleus coherent scattering, where the same effects are expected.

Jul 25 2017

hep-ph arXiv:1707.07228v1

Using the chiral unitary approach in coupled channels and $SU(3)$ symmetry we describe the production of $f_0(500)$, $f_0(980)$ and $a_0(980)$ in the $\chi_{c1} \to \eta \pi^+ \pi^-$ reaction, recently performed by the BESIII collaboration. A very strong peak for the $a_0(980)$ can be seen in the $\eta\pi$ invariant mass, while clear signals for the $f_0(500)$ and $f_0(980)$ appear in the one of $\pi^+\pi^-$. Next, we make predictions for the analogous decay $\eta_c \to \eta \pi^+ \pi^-$, which could also be measured experimentally. We discuss the differences of these reactions which are interesting to test the picture where these scalar mesons are dynamically generated from the interaction of pairs of pseudoscalars.

We use Light-Front Holography to estimate the valence quark and the meson cloud contributions to the nucleon axial form factor. The free couplings of the holographic model are determined by the empirical data and by the information extracted from lattice QCD. The holographic model provides a very good description of the empirical data when we consider a meson cloud mixture of 28% in the physical nucleon state, and also a excellent description of the lattice QCD data for pion masses in the range 200--300 MeV. Our estimate of the meson cloud contribution to the nucleon axial form factor has a very slow falloff with the square momentum transfer when compared to typical estimates from quark models with meson cloud dressing.

Jul 25 2017

hep-ph arXiv:1707.07155v1

An analysis is performed to constrain the $W'$ boson production using the measurement of the s-channel single top quark production cross section. Both phenomenological and statistical approaches are examined and results are presented. In the best case, $W'$ bosons that interact only to the right-handed fermions are excluded below 1390 GeV. To our knowledge, it is the first time that the measured cross section of the s-channel single top quark production from the colliders is used to rule out part of the phase space of the $W'$ boson.

Quasi-PDFs provide a path toward an ab initio calculation of parton distribution functions (PDFs) using lattice QCD. One of the problems faced in calculations of quasi-PDFs is the renormalization of a nonlocal operator. By introducing an auxiliary field, we can replace the nonlocal operator with a pair of local operators in an extended theory. On the lattice, this is closely related to the static quark theory. In this approach, we show how to understand the pattern of mixing that is allowed by chiral symmetry breaking, and obtain a master formula for renormalizing the nonlocal operator that depends on three parameters. We present an approach for nonperturbatively determining these parameters, and show results for two lattice spacings using Wilson twisted mass fermions and for different discretizations of the Wilson line in the nonlocal operator. Using these parameters we show the effect of renormalization on nucleon matrix elements with pion mass approximately 370 MeV, and compare renormalized results for the two lattice spacings. The renormalized matrix elements are consistent among the different Wilson line discretizations and lattice spacings.

Jul 25 2017

hep-ph arXiv:1707.07144v1

The mass spectra and decay properties of heavy quarkonia are computed in nonrelativistic quark-antiquark Cornell potential model. We have employed the numerical solution of Schrödinger equation to obtain their mass spectra using only four parameters namely quark mass ($m_c$, $m_b$) and confinement strength ($A_{c\bar c}$, $A_{b\bar b}$). The spin hyperfine, spin-orbit and tensor components of the one gluon exchange interaction are computed perturbatively to determine the mass spectra of excited $S$, $P$, $D$ and $F$ states. Digamma, digluon and dilepton decays of these mesons are computed using the model parameters and numerical wave functions. The predicted spectroscopy and decay properties for quarkonia are found to be consistent with available experimental observations and results from other theoretical models. We also compute mass spectra and life time of the $B_c$ meson without additional parameters. The computed electromagnetic transition widths of heavy quarkonia and $B_c$ mesons are in tune with available experimental data and other theoretical approaches.

It was established that distribution of the near-zero modes of the Dirac operator is consistent with the Chiral Random Matrix Theory (CRMT) and can be considered as a consequence of spontaneous breaking of chiral symmetry (SBCS) in QCD. The higher-lying modes of the Dirac operator carry information about confinement physics and are not affected by SBCS. We study distributions of the near-zero and higher-lying modes of the overlap Dirac operator within $N_F = 2$ dynamical simulations. We find that distributions of both near-zero and higher-lying modes are the same and follow the Gaussian Unitary Ensemble of Random Matrix Theory. This means that randomness, while consistent with SBCS, is not a consequence of SBCS and is related to some more general property of QCD in confinement regime.

The dissipation of small-scale perturbations in the early universe produces a distortion in the blackbody spectrum of cosmic microwave background photons. In this work, we propose to use these distortions as a probe of the microphysics of dark matter on scales $1\,\textrm{Mpc}^{-1}\lesssim k \lesssim 10^{4}\,\textrm{Mpc}^{-1}$. We consider in particular models in which the dark matter is kinetically coupled to either neutrinos or photons until shortly before recombination, and compute the photon heating rate and the resultant $\mu$-distortion in both cases. We show that the $\mu$-parameter is generally enhanced relative to $\Lambda$CDM for interactions with neutrinos, and may be either enhanced or suppressed in the case of interactions with photons. The deviations from the $\Lambda$CDM signal are potentially within the sensitivity reach of a PRISM-like experiment if $\sigma_{\textrm{DM}-\gamma} \gtrsim 1.1\times10^{-30} \left(m_{\textrm{DM}}/\textrm{GeV}\right) \textrm{cm}^{2}$ and $\sigma_{\textrm{DM}-\nu} \gtrsim 4.8\times 10^{-32} \left(m_{\textrm{DM}}/\textrm{GeV}\right) \textrm{cm}^{2}$ for time-independent cross sections, and $\sigma^{0}_{\textrm{DM}-\gamma} \gtrsim 1.8 \times 10^{-40} \left(m_{\textrm{DM}}/\textrm{GeV}\right) \textrm{cm}^{2}$ and $\sigma^{0}_{\textrm{DM}-\nu} \gtrsim 2.5 \times 10^{-47} \left(m_{\textrm{DM}}/\textrm{GeV}\right) \textrm{cm}^{2}$ for cross sections scaling as temperature squared, coinciding with the parameter regions in which late kinetic decoupling may serve as a solution to the small-scale crisis. Furthermore, these $\mu$-distortion signals differ from those of warm dark matter (no deviation from $\Lambda$CDM) and a suppressed primordial power spectrum (strongly suppressed or a negative $\mu$-parameter), demonstrating that CMB spectral distortion can potentially be used to distinguish between solutions to the small-scale crisis.

Approaches based on viscous hydrodynamics for the hot and dense stage and hadronic transport for the final dilute rescattering stage are successfully applied to the dynamic description of heavy ion reactions at high beam energies. One crucial step in such hybrid approaches is the so called particlization, the transition between the hydrodynamic description to microscopic degrees of freedom. For this purpose, individual particles are sampled on the Cooper-Frye hypersurface. In this work, 4 different realizations of sampling algorithms are compared, where three of them incorporate global conservation laws of quantum numbers in each event. The algorithms are compared within two types of scenarios: simple "box" hypersurface consisting of only one static cell and a typical particlization hypersurface for Au+Au collisions at $\sqrt{s_\mathrm{NN}} = 200$ GeV. For all algorithms the mean multiplicities (or particle spectra) remain unaffected by global conservation laws in the case of large volumes. In contrast, the fluctuations of the particle numbers are affected considerably. The fluctuations of the newly developed SPREW algorithm based on exponential weights and the recently suggested SER algorithm based on ensemble rejection are smaller than without conservation laws and agree with the expectation from the canonical ensemble. The previously applied mode sampling algorithm produces dramatically larger fluctuations, than it is expected in the corresponding microcanonical ensemble, and therefore should be avoided in fluctuation studies. This study might be of interest for investigations of particle fluctuations and correlations, e.g. the suggested signatures for a phase transition or a critical endpoint, in hybrid approaches that are affected by global conservation laws.

Jul 25 2017

hep-ph arXiv:1707.07019v1

Using the Fenchel-Eggleston theorem for convex hulls (an extension of the Caratheodory theorem), we prove that any likelihood can be maximized by either a dark matter 1- speed distribution $F(v)$ in Earth's frame or 2- Galactic velocity distribution $f^{\rm gal}(\vec{u})$, consisting of a sum of delta functions. The former case applies only to time-averaged rate measurements and the maximum number of delta functions is $({\mathcal N}-1)$, where ${\mathcal N}$ is the total number of data entries. The second case applies to any harmonic expansion coefficient of the time-dependent rate and the maximum number of terms is ${\mathcal N}$. Using time-averaged rates, the aforementioned form of $F(v)$ results in a piecewise constant unmodulated halo function $\tilde\eta^0_{BF}(v_{\rm min})$ (which is an integral of the speed distribution) with at most $({\mathcal N}-1)$ downward steps. The authors had previously proven this result for likelihoods comprised of at least one extended likelihood, and found the best-fit halo function to be unique. This uniqueness, however, cannot be guaranteed in the more general analysis applied to arbitrary likelihoods. Thus we introduce a method for determining whether there exists a unique best-fit halo function, and provide a procedure for constructing either a pointwise confidence band, if the best-fit halo function is unique, or a degeneracy band, if it is not. Using measurements of modulation amplitudes, the aforementioned form of $f^{\rm gal}(\vec{u})$, which is a sum of Galactic streams, yields a periodic time-dependent halo function $\tilde\eta_{BF}(v_{\rm min}, t)$ which at any fixed time is a piecewise constant function of $v_{\rm min}$ with at most ${\mathcal N}$ downward steps. In this case, we explain how to construct pointwise confidence and degeneracy bands from the time-averaged halo function. Finally, we show that requiring an isotropic ...

Jul 25 2017

hep-ph arXiv:1707.07016v1

We investigate models of a heavy neutral gauge boson Z' which could explain anomalies in B meson decays reported by the LHCb experiment. In these models, the Z' boson couples mostly to third generation fermions. We show that bottom quarks arising from gluon splitting can fuse into Z' as an essential production mechanism at the LHC, thereby allowing to probe these models. The study is performed within a generic framework for explaining the B anomalies that can be accommodated in well motivated models. The flavor violating b s coupling associated with Z' in such models produces lower bound on the production cross-section which gives rise to a cross-section range for such scenarios for the LHC to probe. Results are presented in Z' -> $\mu \mu$ decays with at least one bottom-tagged jet in its final state. Some parts of the model parameter space become constrained by the existing dimuon-resonance searches by the ATLAS and CMS collaborations. However, the requirement of one or two additional bottom-tagged jets in the final state would allow for probing a larger region of the parameter space of the models at the ongoing LHC program.

Jul 25 2017

hep-ph arXiv:1707.07000v1

Extensions of the Standard Model which address the hierarchy problem and dark matter (DM) often contain top partners and additional resonances at the TeV scale. We explore the phenomenology of a simplified effective model with a vector resonance $Z'$, a fermionic vector-like coloured partner of the top quark $T'$ as well as a scalar DM candidate $\phi$ and provide publicly available implementations in CalcHEP and MadGraph. We study the $pp \to Z' \to T'\overline{T'} \to t\bar{t}\,\phi\phi$ process at the LHC and find that it plays an important role in addition to the $T'\overline{T'}$ production via strong interactions. It turns out that the presence of the $Z'$ can provide a dominant contribution to the $t\bar{t}+E_T^{\text{miss}}$ signature without conflicting with existing bounds from $Z'$ searches in di-jet and di-lepton final states. We find that through this process, the LHC is already probing DM masses up to about 900 GeV and top partner masses up to about 1.5 TeV, thus exceeding the current bounds from QCD production alone almost by a factor of two for both particles.

We provide expressions for the nonperturbative matching of the effective field theory describing dark matter interactions with quarks and gluons to the effective theory of nonrelativistic dark matter interacting with nonrelativistic nucleons. We give the leading and subleading order expressions in chiral counting. In general, a single partonic operator already matches onto several nonrelativistic operators at leading order in chiral counting. Thus, keeping only one operator at the time in the nonrelativistic effective theory does not properly describe the scattering in direct detection. Moreover, the matching of the axial--axial partonic level operator, as well as the matching of the operators coupling DM to the QCD anomaly term, naively include momentum suppressed terms. However, these are still of leading chiral order due to pion poles and can be numerically important. We illustrate the impact of these effects with several examples.

A rare decay $K_L \to \mu^+ \mu^- $ has been measured precisely, while a rare decay $K_S \to \mu^+ \mu^- $ will be observed by an upgrade of the LHCb experiment. Although both processes are almost CP-conserving decays, we point out that an interference contribution between $K_L$ and $K_S$ in the kaon beam emerges from a genuine direct CP violation. It is found that the interference contribution can change $K_S \to \mu^+ \mu^-$ standard-model predictions at $\mathcal{O}(60\%)$. We also stress that an unknown sign of $\mathcal{A}(K_L \to \gamma \gamma)$ can be determined by a measurement of the interference, which can much reduce a theoretical uncertainty of $\mathcal{B}(K_L \to \mu^+ \mu^-)$. We also investigate the interference in a new physics model, where the $\epsilon'_K / \epsilon_K$ tension is explained by an additional $Z$-penguin contribution.

Recent experimental observations of the charged hadron properties in $U+U$ collisions at $193$ GeV contradict many of the theoretical models of particle production including two-component montecarlo Glauber model. The experimental results show a small correlation between the charged hadron multiplicity and the initial geometrical configurations (e.g. body-body, tip-tip etc.) of uranium nuclei. In this article, we have used montecarlo HYDJET++ model to study the charged hadron production in $U+U$ collisions at $193$ GeV center-of-mass energy. This model incorporates glauber type initial conditions along with hard particle production from jet fragmentation and soft particle production from thermal processes. We have calculated the pseudorapidty distribution, transverse momentum distribution and elliptic flow distribution of charged hadrons with different control parameters in various geometrical configurations possible for $U+U$ collision. We find that HYDJET++ model along with Glauber like initial conditions provides the results which qualitatively satisfy the experimental data and thus supports a small correlation between the various properties of charged hadrons and the initial geometrical configuration of $U+U$ collision.