Kinetically modified non-minimal inflation with exponential frame function

We consider supersymmetric (SUSY) and non-SUSY models of chaotic inflation based on the \(\phi ^n\) potential with \(n=2\) or 4. We show that the coexistence of an exponential non-minimal coupling to gravity \(f_\mathcal{R}=\mathrm{e}^{c_\mathcal{R}\phi ^{p}}\) with a kinetic mixing of the form \(f_{\mathrm{K}}=c_{\mathrm{K}}f_\mathcal{R}^m\) can accommodate inflationary observables favored by the Planck and Bicep2/Keck Array results for \(p=1\) and 2, \(1\le m\le 15\) and \(2.6\times 10^{-3}\le r_{\mathcal {R}\mathrm{K}}=c_\mathcal{R}/c_{\mathrm{K}}^{p/2}\le 1,\) where the upper limit is not imposed for \(p=1\). Inflation is of hilltop type and it can be attained for subplanckian inflaton values with the corresponding effective theories retaining the perturbative unitarity up to the Planck scale. The supergravity embedding of these models is achieved employing two chiral gauge singlet supefields, a monomial superpotential and several (semi)logarithmic or semi-polynomial Kähler potentials.     

$$\varvec{B^0\rightarrow K^{*0}\mu ^+\mu ^-}$$ decay in the aligned two-Higgs-doublet model

In the aligned two-Higgs-doublet model, we perform a complete one-loop computation of the short-distance Wilson coefficients \(C_{7,9,10}^{(\prime )}\), which are the most relevant ones for \(b\rightarrow s\ell ^+\ell ^-\) transitions. It is found that, when the model parameter \(\left| \varsigma _{u}\right| \) is much smaller than \(\left| \varsigma _{d}\right| \), the charged scalar contributes mainly to chirality-flipped \(C_{9,10}^\prime \), with the corresponding effects being proportional to \(\left| \varsigma _{d}\right| ^2\). Numerically, the charged-scalar effects fit into two categories: (A) \(C_{7,9,10}^\mathrm {H^\pm }\) are sizable, but \(C_{9,10}^{\prime \mathrm {H^\pm }}\simeq 0\), corresponding to the (large \(\left| \varsigma _{u}\right| \), small \(\left| \varsigma _{d}\right| \)) region; (B) \(C_7^\mathrm {H^\pm }\) and \(C_{9,10}^{\prime \mathrm {H^\pm }}\) are sizable, but \(C_{9,10}^\mathrm {H^\pm }\simeq 0\), corresponding to the (small \(\left| \varsigma _{u}\right| \), large \(\left| \varsigma _{d}\right| \)) region. Taking into account phenomenological constraints from the inclusive radiative decay \(B\rightarrow X_{s}{\gamma }\), as well as the latest model-independent global analysis of \(b\rightarrow s\ell ^+\ell ^-\) data, we obtain the much restricted parameter space of the model. We then study the impact of the allowed model parameters on the angular observables \(P_2\) and \(P_5'\) of \(B^0\rightarrow K^{*0}\mu ^+\mu ^-\) decay, and we find that \(P_5'\) could be increased significantly to be consistent with the experimental data in case B.     

Light deflection and Gauss–Bonnet theorem: definition of total deflection angle and its applications

In this paper, we re-examine the light deflection in the Schwarzschild and the Schwarzschild–de Sitter spacetime. First, supposing a static and spherically symmetric spacetime, we propose the definition of the total deflection angle \(\alpha \) of the light ray by constructing a quadrilateral \(\varSigma ^4\) on the optical reference geometry \({\mathscr {M}}^\mathrm{opt}\) determined by the optical metric \(\bar{g}_{ij}\). On the basis of the definition of the total deflection angle \(\alpha \) and the Gauss–Bonnet theorem, we derive two formulas to calculate the total deflection angle \(\alpha \); (1) the angular formula that uses four angles determined on the optical reference geometry \({\mathscr {M}}^\mathrm{opt}\) or the curved \((r, \phi )\) subspace \({\mathscr {M}}^\mathrm{sub}\) being a slice of constant time t and (2) the integral formula on the optical reference geometry \({\mathscr {M}}^\mathrm{opt}\) which is the areal integral of the Gaussian curvature K in the area of a quadrilateral \(\varSigma ^4\) and the line integral of the geodesic curvature \(\kappa _g\) along the curve \(C_{\varGamma }\). As the curve \(C_{\varGamma }\), we introduce the unperturbed reference line that is the null geodesic \(\varGamma \) on the background spacetime such as the Minkowski or the de Sitter spacetime, and is obtained by projecting \(\varGamma \) vertically onto the curved \((r, \phi )\) subspace \({\mathscr {M}}^\mathrm{sub}\). We demonstrate that the two formulas give the same total deflection angle \(\alpha \) for the Schwarzschild and the Schwarzschild–de Sitter spacetime. In particular, in the Schwarzschild case, the result coincides with Epstein–Shapiro’s formula when the source S and the receiver R of the light ray are located at infinity. In addition, in the Schwarzschild–de Sitter case, there appear order \({\mathscr {O}}(\varLambda m)\) terms in addition to the Schwarzschild-like part, while order \({\mathscr {O}}(\varLambda )\) terms disappear.     

Charmed hadrons from coalescence plus fragmentation in relativistic nucleus-nucleus collisions at RHIC and LHC

In a coalescence plus fragmentation approach we calculate the heavy baryon/meson ratio and the \(p_T\) spectra of charmed hadrons \(D^{0}\), \(D_{s}\) and \(\varLambda _{c}^{+}\) in a wide range of transverse momentum from low \(p_T\) up to about 10 GeV and discuss their ratios from RHIC to LHC energies without any change of the coalescence parameters. We have included the contribution from decays of heavy hadron resonances and also the one due to fragmentation of heavy quarks which do not undergo the coalescence process. The coalescence process is tuned to have all charm quarks hadronizing in the \(p_T\rightarrow 0\) limit and at finite \(p_T\) charm quarks not undergoing coalescence are hadronized by independent fragmentation. The \(p_T\) dependence of the baryon/meson ratios are found to be sensitive to the masses of coalescing quarks, in particular the \(\varLambda _{c}/D^{0}\) can reach values of about \(\mathrm 1\div 1.5 \) at \(p_T \approx \, 3\) GeV, or larger, similarly to the light baryon/meson ratio like \(p/\pi \) and \(\varLambda /K\), however a marked difference is a quite weak \(p_T\) dependence with respect to the light case, such that a larger value at intermediate \(p_T\) implies a relatively large value also for the integrated yields. A comparison with other coalescence model and with the prediction of thermal model is discussed.     

Likelihood analysis of the minimal AMSB model

We perform a likelihood analysis of the minimal anomaly-mediated supersymmetry-breaking (mAMSB) model using constraints from cosmology and accelerator experiments. We find that either a wino-like or a Higgsino-like neutralino LSP, \(\tilde{\chi }^0_{1}\), may provide the cold dark matter (DM), both with similar likelihoods. The upper limit on the DM density from Planck and other experiments enforces \(m_{\tilde{\chi }^0_{1}} \lesssim 3 \,\, \mathrm {TeV}\) after the inclusion of Sommerfeld enhancement in its annihilations. If most of the cold DM density is provided by the \(\tilde{\chi }^0_{1}\), the measured value of the Higgs mass favours a limited range of \(\tan \beta \sim 5\) (and also for \(\tan \beta \sim 45\) if \(\mu > 0\)) but the scalar mass \(m_0\) is poorly constrained. In the wino-LSP case, \(m_{3/2}\) is constrained to about \(900\,\, \mathrm {TeV}\) and \(m_{\tilde{\chi }^0_{1}}\) to \(2.9\pm 0.1\,\, \mathrm {TeV}\), whereas in the Higgsino-LSP case \(m_{3/2}\) has just a lower limit \(\gtrsim 650\,\, \mathrm {TeV}\) (\(\gtrsim 480\,\, \mathrm {TeV}\)) and \(m_{\tilde{\chi }^0_{1}}\) is constrained to \(1.12 ~(1.13) \pm 0.02\,\, \mathrm {TeV}\) in the \(\mu >0\) (\(\mu <0\)) scenario. In neither case can the anomalous magnetic moment of the muon, \((g-2)_\mu \), be improved significantly relative to its Standard Model (SM) value, nor do flavour measurements constrain the model significantly, and there are poor prospects for discovering supersymmetric particles at the LHC, though there are some prospects for direct DM detection. On the other hand, if the \(\tilde{\chi }^0_{1}\) contributes only a fraction of the cold DM density, future LHC Open image in new window -based searches for gluinos, squarks and heavier chargino and neutralino states as well as disappearing track searches in the wino-like LSP region will be relevant, and interference effects enable \(\mathrm{BR}(B_{s, d} \rightarrow \mu ^+\mu ^-)\) to agree with the data better than in the SM in the case of wino-like DM with \(\mu > 0\).     

On the Small Mass Limit of Quantum Brownian Motion with Inhomogeneous Damping and Diffusion

We study the small mass limit (or: the Smoluchowski–Kramers limit) of a class of quantum Brownian motions with inhomogeneous damping and diffusion. For Ohmic bath spectral density with a Lorentz–Drude cutoff, we derive the Heisenberg–Langevin equations for the particle’s observables using a quantum stochastic calculus approach. We set the mass of the particle to equal \(m = m_{0} \epsilon \), the reduced Planck constant to equal \(\hbar = \epsilon \) and the cutoff frequency to equal \(\varLambda = E_{\varLambda }/\epsilon \), where \(m_0\) and \(E_{\varLambda }\) are positive constants, so that the particle’s de Broglie wavelength and the largest energy scale of the bath are fixed as \(\epsilon \rightarrow 0\). We study the limit as \(\epsilon \rightarrow 0\) of the rescaled model and derive a limiting equation for the (slow) particle’s position variable. We find that the limiting equation contains several drift correction terms, the quantum noise-induced drifts, including terms of purely quantum nature, with no classical counterparts.     

Likelihood analysis of the sub-GUT MSSM in light of LHC 13-TeV data

We describe a likelihood analysis using MasterCode of variants of the MSSM in which the soft supersymmetry-breaking parameters are assumed to have universal values at some scale \(M_\mathrm{in}\) below the supersymmetric grand unification scale \(M_\mathrm{GUT}\), as can occur in mirage mediation and other models. In addition to \(M_\mathrm{in}\), such ‘sub-GUT’ models have the 4 parameters of the CMSSM, namely a common gaugino mass \(m_{1/2}\), a common soft supersymmetry-breaking scalar mass \(m_0\), a common trilinear mixing parameter A and the ratio of MSSM Higgs vevs \(\tan \beta \), assuming that the Higgs mixing parameter \(\mu > 0\). We take into account constraints on strongly- and electroweakly-interacting sparticles from \(\sim 36\)/fb of LHC data at 13 TeV and the LUX and 2017 PICO, XENON1T and PandaX-II searches for dark matter scattering, in addition to the previous LHC and dark matter constraints as well as full sets of flavour and electroweak constraints. We find a preference for \(M_\mathrm{in}\sim 10^5\) to \(10^9 \,\, \mathrm {GeV}\), with \(M_\mathrm{in}\sim M_\mathrm{GUT}\) disfavoured by \(\Delta \chi ^2 \sim 3\) due to the \(\mathrm{BR}(B_{s, d} \rightarrow \mu ^+\mu ^-)\) constraint. The lower limits on strongly-interacting sparticles are largely determined by LHC searches, and similar to those in the CMSSM. We find a preference for the LSP to be a Bino or Higgsino with \(m_{\tilde{\chi }^0_{1}} \sim 1 \,\, \mathrm {TeV}\), with annihilation via heavy Higgs bosons H / A and stop coannihilation, or chargino coannihilation, bringing the cold dark matter density into the cosmological range. We find that spin-independent dark matter scattering is likely to be within reach of the planned LUX-Zeplin and XENONnT experiments. We probe the impact of the \((g-2)_\mu \) constraint, finding similar results whether or not it is included.     

CP violations in predictive neutrino mass structures

We study the CP-violation effects from two types of neutrino mass matrices with (i) \((M_\nu )_{ee}=0\), and (ii) \((M_\nu )_{ee}=(M_\nu )_{e\mu }=0\), which can be realized by the high-dimensional lepton number violating operators \(\bar{\ell }_R^c\gamma ^\mu L_L (D_\mu \Phi )\Phi ^2\) and \(\bar{\ell }_R^c l_R (D_\mu {\Phi })^2\Phi ^2\), respectively. In (i), the neutrino mass spectrum is in the normal ordering with the lightest neutrino mass within the range \(0.002\,\mathrm{eV}\lesssim m_0\lesssim 0.007\,\mathrm{eV}\). Furthermore, for a given value of \(m_0\), there are two solutions for the two Majorana phases \(\alpha _{21}\) and \(\alpha _{31}\), whereas the Dirac phase \(\delta \) is arbitrary. For (ii), the parameters of \(m_0\), \(\delta \), \(\alpha _{21}\), and \(\alpha _{31}\) can be completely determined. We calculate the CP-violating asymmetries in neutrino–antineutrino oscillations for both mass textures of (i) and (ii), which are closely related to the CP-violating Majorana phases.     

Holographic screening length in a hot plasma of two sphere

We study the screening length \(L_{\mathrm{max}}\) of a moving quark–antiquark pair in a hot plasma, which lives in a two sphere, \(S^2\), using the AdS/CFT correspondence in which the corresponding background metric is the four-dimensional Schwarzschild–AdS black hole. The geodesic of both ends of the string at the boundary, interpreted as the quark–antiquark pair, is given by a stationary motion in the equatorial plane by which the separation length L of both ends of the string is parallel to the angular velocity \(\omega \). The screening length and total energy H of the quark–antiquark pair are computed numerically and show that the plots are bounded from below by some functions related to the momentum transfer \(P_c\) of the drag force configuration. We compare the result by computing the screening length in the reference frame of the moving quark–antiquark pair, in which the background metrics are “Boost-AdS” and Kerr–AdS black holes. Comparing both black holes, we argue that the mass parameters \(M_{\mathrm{Sch}}\) of the Schwarzschild–AdS black hole and \(M_{\mathrm{Kerr}}\) of the Kerr–AdS black hole are related at high temperature by \(M_{\mathrm{Kerr}}=M_{\mathrm{Sch}}(1-a^2l^2)^{3/2}\), where a is the angular momentum parameter and l is the AdS curvature.     

Nuclear suppression of the $$\phi $$ meson yields with large $$p_T$$pT at the RHIC and the LHC

We calculate \(\phi \) meson transverse momentum spectra in \(\mathrm{p}+\mathrm{p}\) collisions as well as their nuclear suppressions in central \(\mathrm{A}+\mathrm{A}\) collisions both at the RHIC and the LHC in LO and NLO with the QCD-improved parton model. We have included the parton energy loss effect in a hot/dense QCD medium with the effectively medium-modified \(\phi \) fragmentation functions in the higher-twist approach of jet quenching. The nuclear modification factors of the \(\phi \) meson in central \(\mathrm{Au}+\mathrm{Au}\) collisions at the RHIC and central \(\mathrm{Pb}+\mathrm{Pb}\) collisions at the LHC are provided, and nice agreement of our numerical results at NLO with the ALICE measurement is observed. Predictions of the yield ratios of neutral mesons such as \(\phi /\pi ^0\), \(\phi /\eta \) and \(\phi /\rho ^0\) at large \(p_T\) in relativistic heavy-ion collisions are also presented for the first time.     

On quantum analogue of dynamical stabilisation of inverted harmonic oscillator by time periodical uniform field

Quantum analogue of stabilised forced oscillations around an unstable equilibrium position is explored by solving the non-stationary Schrödinger equation (NSE) of the inverted harmonic oscillator (IHO) driven periodically by spatial uniform field of frequency \(\Omega \), amplitude \(F_{0}\) and phase \(\phi \), i.e. the system with the Hamiltonian of \(\hat{{H}}=(\hat{{p}}^{2}/2m)-(m\omega ^{2}x^{2}/2)-F_0 x\sin \) \(\left( {\Omega t+\phi } \right) \). The NSE has been solved both analytically and numerically by Maple 15 in dimensionless variables \(\xi = x\sqrt{m\omega /\hbar }\hbox {, }f_0 =F_0 /\omega \sqrt{\hbar m\omega }\) and \(\tau =\omega t\). The initial condition (IC) has been specified by the wave function (w.f.) of a generalised Gaussian type which suits well the corresponding quantum IC operator. The solution obtained demonstrates the non-monotonous behaviour of the coordinate spreading \(\sigma \left( \tau \right) \hbox { =}\sqrt{\big ( {\overline{\Delta \xi ^{2}\big ( \tau \big )} } \big )}\) which decreases first from quite macroscopic values of \(\sigma _{0} =2^{12,\ldots ,25}\) to minimal one of \(\sim \!(1/\sqrt{2})\) at times \(\tau <\tau _0 =0.125\ln \!\left( {16\sigma _0^4 +1} \right) \) and then grows back unlimitedly. For certain phases \(\phi \) depending on the \(\Omega /\omega \) ratio and \(n=\log _2\!\sigma _0 \), the mass centre of the packet \(\xi _{\mathrm {av}}( \tau )= \overline{\hat{{x}}(\tau )} \cdot \sqrt{m\omega /\hbar }\) delays approximately two natural ‘periods’ \(\sim \!(4\pi /\omega )\) in the area of the stationary point and then escapes to ‘\(+\)’ or ‘?’ infinity in a bifurcating way.  For ‘resonant’ \(\Omega =\omega \), the bifurcation phases \(\phi \) fit well with the regression formula of Fermi–Dirac type of argument n with their asymptotic \(\phi ( {\Omega ,n\rightarrow \infty } )\) obeying the classical formula \(\phi _{\mathrm {cl}} ( \Omega )=-\hbox {arctg} \, \Omega \) for initial energy \(E = 0\) in the wide range of \(\Omega =2^{-4},...,2^{7}\).     

CGC/saturation approach: re-visiting the problem of odd harmonics in angular correlations

In this paper we demonstrate that the selection of events with different multiplicities of produced particles, leads to the violation of the azimuthal angular symmetry, \(\phi \rightarrow \pi - \phi \). We find for LHC and lower energies, that this violation can be so large for the events with multiplicities \(n \ge 2 \bar{n}\), where \(\bar{n}\) is the mean multiplicity, that it leads to almost no suppression of \(v_n\), with odd n. However, this can only occur if the typical size of the dipole in DIS with a nuclear target is small, or \(Q^2 \,>\,Q^2_s\left( A; Y_{\mathrm{min}},b\right) \), where \(Q_s\) is the saturation momentum of the nucleus at \(Y = Y_{\mathrm{min}}\). In the case of large sizes of dipoles, when \(Q^2 \,<\,Q^2_s\left( A; Y_{\mathrm{min}},b\right) \), we show that \(v_n =0\) for odd n. Hadron-nucleus scattering is discussed.     

Systematic study of rigid triaxiality in Ba–Pt nuclei and role of $$$$ subshell effect

A systematic variation of shape variables \(\beta \) and \(\gamma \) with N and \(N_\mathrm {p}N_\mathrm {n}\) is studied in the framework of an asymmetric rotor model of Davydov and Filippov for the \(Z=\) 50–82, \(N=\) 82–126 major shell space. The role of the \(Z=64\) subshell in producing smooth systematics has been discussed. The quadrupole moments are extracted after considering both axially symmetric and axially asymmetric nuclei. The correlation of \(\beta \) with \(\gamma \) together with the measured quadrupole moments indicates that \(\gamma \)-rigidity is better observed in nuclei with modest deformation.     

Scattering of kinks of the sinh-deformed $$\varphi ^4$$ model

We consider the scattering of kinks of the sinh-deformed \(\varphi ^4\) model, which is obtained from the well-known \(\varphi ^4\) model by means of the deformation procedure. Depending on the initial velocity \(v_\mathrm {in}\) of the colliding kinks, different collision scenarios are realized. There is a critical value \(v_\mathrm {cr}\) of the initial velocity, which separates the regime of reflection (at \(v_\mathrm {in}>v_\mathrm {cr}\)) and that of a complicated interaction (at \(v_\mathrm {in}<v_\mathrm {cr}\)) with kinks’ capture and escape windows. Besides that, at \(v_\mathrm {in}\) below \(v_\mathrm {cr}\) we observe the formation of a bound state of two oscillons, as well as their escape at some values of \(v_\mathrm {in}\).     

Stochastic Averaging Principle for Spatial Birth-and-Death Evolutions in the Continuum

We study a spatial birth-and-death process on the phase space of locally finite configurations \({\varGamma }^+ \times {\varGamma }^-\) over \({\mathbb {R}}^d\). Dynamics is described by an non-equilibrium evolution of states obtained from the Fokker-Planck equation and associated with the Markov operator \(L^+(\gamma ^-) + \frac{1}{\varepsilon }L^-\), \(\varepsilon > 0\). Here \(L^-\) describes the environment process on \({\varGamma }^-\) and \(L^+(\gamma ^-)\) describes the system process on \({\varGamma }^+\), where \(\gamma ^-\) indicates that the corresponding birth-and-death rates depend on another locally finite configuration \(\gamma ^- \in {\varGamma }^-\). We prove that, for a certain class of birth-and-death rates, the corresponding Fokker-Planck equation is well-posed, i.e. there exists a unique evolution of states \(\mu _t^{\varepsilon }\) on \({\varGamma }^+ \times {\varGamma }^-\). Moreover, we give a sufficient condition such that the environment is ergodic with exponential rate. Let \(\mu _{\mathrm {inv}}\) be the invariant measure for the environment process on \({\varGamma }^-\). In the main part of this work we establish the stochastic averaging principle, i.e. we prove that the marginal of \(\mu _t^{\varepsilon }\) onto \({\varGamma }^+\) converges weakly to an evolution of states on \({\varGamma }^+\) associated with the averaged Markov birth-and-death operator \({\overline{L}} = \int _{{\varGamma }^-}L^+(\gamma ^-)d \mu _{\mathrm {inv}}(\gamma ^-)\).     

On Blow-up Profile of Ground States of Boson Stars with External Potential

We study minimizers of the pseudo-relativistic Hartree functional \({\mathcal {E}}_{a}(u):=\Vert (-\varDelta +m^{2})^{1/4}u\Vert _{L^{2}}^{2}+\int _{{\mathbb {R}}^{3}}V(x)|u(x)|^{2}\mathrm{d}x-\frac{a}{2}\int _{{\mathbb {R}}^{3}}(\left| \cdot \right| ^{-1}\star |u|^{2})(x)|u(x)|^{2}\mathrm{d}x\) under the mass constraint \(\int _{{\mathbb {R}}^3}|u(x)|^2\mathrm{d}x=1\). Here \(m>0\) is the mass of particles and \(V\ge 0\) is an external potential. We prove that minimizers exist if and only if a satisfies \(0\le a<a^{*}\), and there is no minimizer if \(a\ge a^*\), where \(a^*\) is called the Chandrasekhar limit. When a approaches \(a^*\) from below, the blow-up behavior of minimizers is derived under some general external potentials V. Here we consider three cases of V: trapping potential, i.e. \(V\in L_{\mathrm{loc}}^{\infty }({\mathbb {R}}^3)\) satisfies \(\lim _{|x|\rightarrow \infty }V(x)=\infty \); periodic potential, i.e. \(V\in C({\mathbb {R}}^3)\) satisfies \(V(x+z)=V(x)\) for all \(z\in \mathbb {Z}^3\); and ring-shaped potential, e.g. \( V(x)=||x|-1|^p\) for some \(p>0\).     

Equilibration in the Kac Model Using the GTW Metric $$d_2$$

We use the Fourier based Gabetta–Toscani–Wennberg metric \(d_2\) to study the rate of convergence to equilibrium for the Kac model in 1 dimension. We take the initial velocity distribution of the particles to be a Borel probability measure \(\mu \) on \(\mathbb {R}^n\) that is symmetric in all its variables, has mean \(\vec {0}\) and finite second moment. Let \(\mu _t(dv)\) denote the Kac-evolved distribution at time t, and let \(R_\mu \) be the angular average of \(\mu \). We give an upper bound to \(d_2(\mu _t, R_\mu )\) of the form \(\min \left\{ B e^{-\frac{4 \lambda _1}{n+3}t}, d_2(\mu ,R_\mu )\right\} ,\) where \(\lambda _1 = \frac{n+2}{2(n-1)}\) is the gap of the Kac model in \(L^2\) and B depends only on the second moment of \(\mu \). We also construct a family of Schwartz probability densities \(\{f_0^{(n)}: \mathbb {R}^n\rightarrow \mathbb {R}\}\) with finite second moments that shows practically no decrease in \(d_2(f_0(t), R_{f_0})\) for time at least \(\frac{1}{2\lambda }\) with \(\lambda \) the rate of the Kac operator. We also present a propagation of chaos result for the partially thermostated Kac model in Tossounian and Vaidyanathan (J Math Phys 56(8):083301, 2015).     

Likelihood analysis of the pMSSM11 in light of LHC 13-TeV data

We use MasterCode to perform a frequentist analysis of the constraints on a phenomenological MSSM model with 11 parameters, the pMSSM11, including constraints from \(\sim 36\)/fb of LHC data at 13 TeV and PICO, XENON1T and PandaX-II searches for dark matter scattering, as well as previous accelerator and astrophysical measurements, presenting fits both with and without the \((g-2)_\mu \) constraint. The pMSSM11 is specified by the following parameters: 3 gaugino masses \(M_{1,2,3}\), a common mass for the first-and second-generation squarks \(m_{\tilde{q}}\) and a distinct third-generation squark mass \(m_{\tilde{q}_3}\), a common mass for the first-and second-generation sleptons \(m_{\tilde{\ell }}\) and a distinct third-generation slepton mass \(m_{\tilde{\tau }}\), a common trilinear mixing parameter A, the Higgs mixing parameter \(\mu \), the pseudoscalar Higgs mass \(M_A\) and \(\tan \beta \). In the fit including \((g-2)_\mu \), a Bino-like \(\tilde{\chi }^0_{1}\) is preferred, whereas a Higgsino-like \(\tilde{\chi }^0_{1}\) is mildly favoured when the \((g-2)_\mu \) constraint is dropped. We identify the mechanisms that operate in different regions of the pMSSM11 parameter space to bring the relic density of the lightest neutralino, \(\tilde{\chi }^0_{1}\), into the range indicated by cosmological data. In the fit including \((g-2)_\mu \), coannihilations with \(\tilde{\chi }^0_{2}\) and the Wino-like \(\tilde{\chi }^\pm _{1}\) or with nearly-degenerate first- and second-generation sleptons are active, whereas coannihilations with the \(\tilde{\chi }^0_{2}\) and the Higgsino-like \(\tilde{\chi }^\pm _{1}\) or with first- and second-generation squarks may be important when the \((g-2)_\mu \) constraint is dropped. In the two cases, we present \(\chi ^2\) functions in two-dimensional mass planes as well as their one-dimensional profile projections and best-fit spectra. Prospects remain for discovering strongly-interacting sparticles at the LHC, in both the scenarios with and without the \((g-2)_\mu \) constraint, as well as for discovering electroweakly-interacting sparticles at a future linear \(e^+ e^-\) collider such as the ILC or CLIC.