abstract: This work proposes a novel framework for the development of new products and services in transportation through an open innovation approach based on automatic content analysis of social media data. The framework is able to extract users comments from Online Social Networks (OSN), to process and analyze text through information extraction and sentiment analysis techniques to obtain relevant information about product reception on the market. A use case was developed using the mobile application Uber, which is today one of the fastest growing technology companies in the world. We measured how a controversial, highly diffused event influences the volume of tweets about Uber and the perception of its users. While there is no change in the image of Uber, a large increase in the number of tweets mentioning the company is observed, which meant a free and important diffusion of its product.
author: Daniela Ulloa, Pedro Saleiro, Rosaldo J. F. Rossetti and Elis Regina Silva [^1] [^2] [^3]
title:  **Mining Social Media for Open Innovation in Transportation Systems**

# INTRODUCTION

The connotation of a smart city and a smart transportation system represents high innovation and disruption in management and policy, as well as technology. In the traditional value creation process, companies and consumers had distinct roles of production and consumption. Now, literature about innovation posits that this process lies in an altogether approach to value creation, based on an individual-centered co-creation of value between consumers and companies \[1\]. Through users' inputs, it is possible to get improvements on the innovation process because their opinions represent substantial value, bringing benefits for companies that adopt them. Analyzing the field of software development, the involvement of the end user during the development process is not new: there is a need for bringing users' inputs by collecting their needs through requirements engineering processes, which can be seen as a way of co-creating with users. This process is mainly done through interviews, questionnaires and surveys. Nowadays, this type of data is no longer enough: companies now need access to digital behavioral data as well \[2\]. With the development of Web 2.0, crowdsourcing emerged as a way of reuniting every possible type of user, allowing a "crowd" of people to collaborate through Web technologies to harness creative collective solutions. While this approach has been used for long time, it is only recently that businesses and other entities have started increasingly turning to web crowdsourcing as a means of obtaining external expertise \[3\]. However, both physical and virtual requirements elicitation techniques require manual collection and/or analysis of data, which is costly and time-consuming.

In this regard it is interesting to make a sort of passive crowdsourcing. On the one hand people today comment almost everything about products, releases or news in OSN and on the other technology for word processing on OSN such as sentiment analysis is already available and accessible. Therefore, picking up the available technology and harnessing the huge quantity of data about people's sentiment, preferences, thinking and opinions available in OSN, we propose a novel framework for the development of new products and services in transportation using an open innovation approach, based on automatic content analysis of social media data. This work follows up previous effort to explore the pontential use of social media in the domain of traffic and transportation engineering \[32, 33, 34\].

Social media applied to co-creation methods has received an increased deal of attention over the past years, despite this update, the literature is mainly theoretical, needing more empirical research, because the role of social media in companies' innovation process is not yet clear and its potential on becoming a source of users' information for firms to improve their NPD process results has not being developed or tested.

To achieve the objective, we apply an architecture pipeline where we extract users' opinions from the Web, according to previously defined heuristics, for being standardized and having a set of meaningful, general information. Next, using filters, we obtain a knowledge database, where the information gathered before is processed by a natural language processing (NLP), which involves natural language understanding, that is, enabling computers to derive meaning from human or natural language input. Once only data related to the field to be researched is held, we perform quantitative and sentiment analysis processes, with the objective of analyzing the opinion conveyed to each social content, assigning a categorical (positive, negative, neutral) or numerical sentiment score. Finally, the results are shown, through easily understandable ways of visualization.

# RELATED WORK

Three main topics complete the structure of this paper: how the transportation paradigm has changed thanks to the emergence of new needs of transportation technologies in the context of smart cities development; the shift of focus on creating value in companies, seeking more involvement of key actors for an open approach to value creation; and previous research in the area of transportation, smart cities and alike who use text processing of social media.

## Development of new Transportation Paradigms

Many studies in the literature have been intensified to introduce new technological solutions in the scope of services in mobility to minimize problems arising from large cities such as traffic congestion, greenhouse gas emissions, fuel consumption, among others \[4\]. Many works propose models of shared services to make transport services more efficient to both business models and reduction of travel costs for passengers \[5\]. Shared transportation services that has been most discussed in the literature are Carpooling and Taxi Sharing, these services are aimed at the sharing of a single vehicle for multiple passengers on a journey, and may or may not have local searches and similar destinations \[6\].

Mobile applications are gaining considerable attention in shared transport services because of their usefulness, not only by the scientific community, but also by many companies. Services such as "Taxi for Two´´ in United kingdom, where the taxi service is shared by two people who have similar travel routes and share the costs of travel \[8\]. Another existing service on shared transport is Caronetas-Carpool Smart, whose goal is to share a taxi or a private car to employees or company employees in São Paulo-Brazil \[9\]. There is also the app "BlaBlaCar´´ a community of sharing private hitchhikers, present in over 22 countries and with over 30 million members, and allow users to travel by car shares both at national and international level, obtaining economic and social benefits \[10\].

One of the most disruptive applications in the last decade appeared in 2009: Uber, today's one of the most valuable technology companies in the world, the one with fastest growing and the most highly rated of all time. It develops, markets and operates the Uber mobile App, released in June 2010, which allows consumers with smartphones to submit a trip request, which is then routed to Uber drivers who use their own cars.

## Open Innovation in the NPD process

Successful NPD requires two essential types of information: (1) information about customer needs and (2) information about how best to solve these needs \[11\].

Users have become a valuable source of innovation for companies. Today's consumers are more connected, informed, empowered and active \[12\], having different experiences and background, which increase competitive advantages for firms and improve efficiency, turning drive down costs.

Because of their characteristics, companies often take advantage of years of accumulated product knowledge and experience about the precise needs and problems that consumers experience to innovate more successfully \[13\]. During the years, the degree of involvement of customers have changed, moving from being integrated in the front end innovation (early stages of the process) to a participation in every stage. Given this evolution it is important for companies to change the way they relate with users, shifting their proposition on value creation from a firm-centric view to a customer-centric or co-creation view \[12\].

Marketing practice and theory increasingly recognize the potential that co-creation has for the firm's performance \[12\]. By successfully implementing and managing co-creation, a firm can create two significant sources of competitive advantages \[14\]: (1) productivity gains through increased efficiency and (2) improved effectiveness.

Co-creation increases productivity and efficiency gains through cost-minimization, since employees' input can be substituted with consumers' input in the product/service development \[15\]. Cost saving arise from various sources: virtually cost-less acquisition of consumer ideas and outsourcing of NPD efforts, which decrease the need for inputs from traditional market research and employees \[16\], reduced risk of product failure and inventory holding costs \[17\], faster speed-to-market \[18\], and post-launch gains through continuous product improvements and exploration of additional usages \[19\]. These outcomes may directly influence organizational performance, increasing the efficiency of operations, product/service turnover, employee satisfaction and ultimately, revenues and profitability \[20\].

## Social Media data processing

In \[26\] the authors presented the domain-agnostic framework for intelligent processing of textual streams from social networks. The framework implemented a pipeline of techniques for semantic representation, sentiment analysis, and automatic content classification. The authors used two scenarios related to smart cities to prove the effectiveness of the framework. In the first scenario they monitored the recovery of the state's capital of L'Aquila's city after the earthquake in April 2009, while in the second scenario was performed a semantic analysis of the content posted on the social network Twitter, building an Italian hate map. With the results of the analysis, the authors confirmed that the methodology used for textual content analysis can provide interesting discoveries to improve understanding of the very complex phenomena.

In the study conducted by Saleiro et al. \[27\] it was presented the POPmine system, an open source system that is able to collect web-based texts from conventional media such as Twitter and blogs and track political views. The authors emphasize that it is the first platform that integrates data collection, information extraction, and opinion mining and public opinion data visualization. The platform processes the text, recognizes the issues and political actors also analyzing relevant language units and generating indicators of dimensions of frequencies and polarities of references to political actors through sources, types of sources and over time.

In \[28\] TweetCred was developed and evaluated, a real-time web-based system to automatically evaluate the credibility of the Twitter content. The system used a lexicon with 45 feature that included the content of tweets, author of the features and information about URLs. The system was trained from opinions obtained by users on a particular area of study. The system computed rankings credibility of 1.1 million tweets and gave feedback to 936 tweets, of this total 43 of tweets, users agreed with the contents of credibility score calculated by TweetCred. TweetCred was developed in the form of a browser extension, a web application and an API available in http://twitdigest.iiitd.edu.in/TweetCred/.

# FRAMEWORK

We design a platform with the objective of extracting, analyzing and aggregating data in order to produce valuable results by performing four major components: data collection, information extraction, sentiment analysis and analytics.

The general architecture of the platform is showed in Figure 1. Next, a description of each component is provided.

-   **Data collection and processing**

For conducting the data collection, the user (company) defines the extraction heuristics, which lead the search and allow connecting to a social network (Facebook, Twitter, Instagram, etc.) The data collection components crawl data from specific data sources, which implement specific web interfaces (e.g. RSS feeds, Twitter API). Each data source must have its own data collection module, which in turn connects to the system using REST services. The system stores data collected in a document oriented NoSQL database (MongoDB). This configuration allows modularity and flexibility, allowing the possibility of developing specific data collection components tailored to specific data sources \[27\].

The heuristics defined by the user can be one of three options: • Search terms: extracts all data that match a specific term • Users pool: extract all data posted by a specific user, given its user name • Geo reference: extract all geolocalized data

Next, it is necessary to process the data extracted. For this, the user defines the processing heuristics, which designate specific product's aspects that are intended to be monitored. Besides, the user also defines the indicators by which the product's aspects are measured and analyzed.

When having the user's input, the data collection components crawl data from specific data sources, in our case, a specific OSN, which implement specific web interfaces.

Data are continuously extracted and stored, for the whole time interval an analysis is running.

-   **Information extraction**

This component addresses two tasks: Named Entity Recognition and Named Entity Disambiguation. This component further processes the content gathered through the data collection process before aggregating, filtering and presenting it in the Analytics UI. This requirement is due to the fact that the extraction process is carried out through a keyword-based matching, so there is a chance that irrelevant content is extracted \[26\]. When monitoring opinions of a given entity it is first necessary to guarantee that all data are relevant to that entity.

We envision an application scenario where we need to track specific products entities. Usually this type of entities are well known therefore we opted to use a knowledge base to provide metadata about the target entities, namely the most common surface forms of their names. Once we had the list of surface forms to search for we applied a sequential classification approach using a prefix tree to detect mentions. This method can result in noisy mentions when applied to Social Media. For instance, an opinion containing the Word "Apple" can be related with more than one entity, such as the technology company or a fruit. Furthermore, posts in OSN are usually short which results in a reduced context for entity disambiguation. When monitoring the opinion of a given entity on OSN, it is first necessary to guarantee that all posts are relevant to that entity.

It is very important to conduct this part of the framework in order to obtain a richer and fine-grained semantic representation, which is very valuable and useful for final users, in virtue of its transparency and readability. Once the filtering process is done, the final data is stored into the Social Media Database \[29\].

Consequently, other processing tasks, such as sentiment analysis will benefit from filtering out noise in the data stream.

-   **Sentiment analysis**

It is the computational study of people's opinions, appraisals, attitudes, and emotions toward entities such as products, services, organizations, individuals, events, and their different aspects \[30\].

Its basic task is to extract and summarize people's opinions expressed on entities and aspects of entities. It consists of three core sub-tasks:

-   Identifying and extracting entities in evaluative texts

-   Identifying and extracting aspects of the entities

-   Determining sentiment polarities on entities and aspects of entities

To this aim, we implement a lexicon-based algorithm for sentiment analysis. Lexicon-based approaches \[31\] infer the sentiment conveyed by a piece of text by relying on (external) lexical resources, which map each term to a categorical (positive, negative, neutral) or numerical sentiment score.

We built a sentiment analysis module to detect and classify opinionated posts, mentioning at least one entity, as expressing a positive, negative or neutral opinion. This module was built with a supervised text classification approach based on the bag-of-words assumption, that is, ignoring word order and representing messages as high-dimensional feature vectors with the size of the vocabulary. Then, a manually labeled dataset was used to train a linear classifier that estimates the probability of a sentiment label, given a message.

One of the main challenges of developing text classification systems for social media, is dealing with the small size of the messages and the large vocabularies needed to account for the significant lexical variation caused by the informal nature of the text. This causes feature vectors to become very sparse, hampering the generalization of the classifiers. Therefore, we first normalized the messages and then enriched the bag-of-words model with additional features based on dense word representations. These word representations, were derived from an unlabeled corpus of 10 Million tweets using two well- known unsupervised feature learning methods: (i) Brown (Brown et al., 1992) hierarchical clustering algorithm, that groups together words that tend to appear in the same contexts. Words are thus represented as the cluster to which they belong; (ii) Mikolov (Mikolov et al., 2013) neural language models that induce dense word vectors (commonly referred to as word embeddings) by training a neural network to maximize the probability that words within a given window size are predicted correctly.

This component enriches the comprehension of the content by analyzing the opinion conveyed by each social content. With this, we aim to determine the attitude of a user with respect to the aspects indicators previously defined.

-   **Analytics**

A set of analysis can be executed with our architecture. This component let the user visualize and interact with the results of the analysis.

-   Mentions to mobility mobile Apps.

-   Sentiment conveyed through opinions regarding mobility mobile Apps aspects.

-   Evaluation of the performance of mentioned mobility mobile Apps aspects.

-   Following the mapping from section 3.1 determine which innovation characteristics the mentioned mobility mobile Apps match.

Through the results the framework provides, companies can improve their NPD process to have clarity about the perception that users have of their products, and what are the ideas for improvement that they propose on the Web, serving as background for companies' decision making.

# USE CASE

As exemplification of the methodology presented, a use case is developed, taking the mobile App Uber. The aim is that the proposed framework is comprehensible, using a popular and known mobile App that has come in a disruptive way, to change the paradigm of transportation systems.

Uber Technologies Inc. was founded in 2009 in San Francisco, USA. It develops, markets and operates the Uber mobile App, released in June 2010, which allows consumers with smartphones to submit a trip request, which is then routed to Uber drivers who use their own cars. As of June 3rd, 2016 the service is available in over 66 countries and 469 cities worldwide.

As in other countries, the arrival of Uber to Portugal has not been without controversy. While users advocate its use indicating that the advantages and ease far outweigh traditional taxis, the latter have not received well the arrival of the mobile App, considering it an unfair competition. On April 29, 2016, about four thousand taxis took to the streets of Lisbon in a slow march to protest Uber. The manifestation began at 9:00 am at Park of Nations, and ended after 13:00 at the Parliament. At the same time there were manifestations also in Porto, where about 1500 taxis participated, and in Faro, where adherence was around 500 taxi drivers.

The first filter applied was the search entity "Uber". To develop the use case we decided to use tweets generated during manifestation's previous and next day, so a date filter is established, taking the tweets of the day before and the day after the protest: April 28th (247 tweets) and 30th (180 tweets), 2016. The purpose of this is to analyze, through sentiment analysis, how the protest influenced Uber perception in the population, by comparing the sentiment attached on the tweets before the protest and after the protest. When there are situations of social revolt, also publicized by the media, people tend to react more in OSN, expressing opinions and thoughts; mainly on Twitter, which has become a tool that actually measures the pulse of humanity: hundreds of global social movements are instantly organized and shared through this social network, influencing the perception about the discussed topic.

The whole atmosphere surrounding the demonstration of taxi drivers against Uber in Portugal had a major impact on OSN, reflected on Twitter in a significant increase in entries and comments.

While the data shows that, during that day, the brand Uber became clearly more visible in OSN (Figure 2), the sentiment of users about the brand practically did not changed. While the positive and negative feelings decreased by 3 and 1 percentage point, respectively, the neutral sentiment increased by 3 percentage points (Figure 3 and Figure 4).

A massive and opinion-conflicting event, which was highly diffused through OSN and media influenced greatly the visibility of the main matter and installed the brand in the minds of users; but it did not influence the perception that users had about it, because to have enough arguments to change that view requires a deeper process of experimentation and a longer period of time. The main topic of discussion was not focused on functional characteristics of Uber, but rather in its dispute with traditional taxi drivers, so the discussion that took place on Twitter had more to do with this than with an assessment about quality service of Uber. In fact, the aim of the taxi drivers was damaging the image of Uber, but they did not achieve it. This event meant a lot of free publicity for Uber.

# CONCLUSIONS AND FUTURE WORK

The objective of this paper was to assess the impact of OSN as a user-centered information provider for companies to improve their NPD processes, proposing a framework for companies, specifically in the transportation area, to harvest information provided by users. The purpose of creating a system with this information is to support the companies' decision-making process when taking into account users' opinions and sentiment about products, for improving their qualities and better meet users' needs.

As exemplification of the methodology proposed, we developed a use case taking as example the mobile App Uber, chosen for being a disruptive technology, which has quickly become popular in over 66 countries around the world, transforming the traditional transportation concept. To collect the data we applied two filters: the search entity "Uber´´ and the geolocation "Portugal".

In order to evaluate the polarity shift of the population in terms of the perception of Uber because of a protest by traditional taxi drivers, we considered the tweets containing the word "Uber" generated during the days before and after the protest. While there was a significant increase in entries and comments on Twitter about Uber, the sentiment of users about the brand practically did not change. The aim of the taxi drivers was damaging the image of Uber, but this did not happen, in fact, this event meant a lot of free publicity for Uber. A massive, highly diffused and opinion-conflicting event influenced greatly the visibility of the matter under discussion; but it did not influence the perception that users had about it, because to have enough arguments to change that view requires a deeper process of experimentation and a longer period of time.

The results of the use case are useful for Uber and taxis, maybe especially up to taxis. In addition to the framework be used to evaluate the products, it can also be used to monitor what is said about competitors. And in this case, people/users in Portugal do not discuss the question of the legality of the Uber but above all the advantages of its service. Perhaps taxi drivers association should use these results to try to improve the services provided by taxis, taking into account what people value about Uber.

We believe this research can be deepen by developing the platform for actually having a system to constantly monitor what user express about a product on the Web. It is also important to validate the framework in a closer to reality environment, using actual cases of product development. To this end, next steps in this research include coupling information from social networks to the MAS-Ter Lab framework , so as to foster and enhance the analysis of new transportation paradigms, in which social simulation and artificial societies are underlying concepts .

[^1]: D. Ulloa is supported by Erasmus Mundus Action 2 Program with an MSc scholarship grant; E. Silva is supported by CNPQ (Conselho Nacional de Desenvolvimento Científico e Tecnológico), Brazil, with doctoral fellowship grant 234933/2014-3.

[^2]: D. Ulloa, R.J.F. Rossetti, and E. Silva are with the Artificial Intelligence and Computer Science Laboratory (LIACC), Department of Informatics Engineering, Faculty of Engineering, University of Porto, Rua Dr. Roberto Frias, s/n, 4200-465 Porto, Portugal. E-mails: `{up201401125, rossetti, firstname.lastname@example.com.`

[^3]: P. Saleiro is with Sapo Labs, Faculty of Engineering, University of Porto, Rua Dr. Roberto Frias, s/n, 4200-465 Porto, Portugal. E-mail: `firstname.lastname@example.com`


abstract: We derive the latest constraints on various simplified models of natural SUSY with light higgsinos, stops and gluinos, using a detailed and comprehensive reinterpretation of the most recent 13 TeV ATLAS and CMS searches with $\sim 15$ fb$^{-1}$ of data. We discuss the implications of these constraints for fine-tuning of the electroweak scale. While the most "vanilla" version of SUSY (the MSSM with $R$-parity and flavor-degenerate sfermions) with 10% fine-tuning is ruled out by the current constraints, models with decoupled valence squarks or reduced missing energy can still be fully natural. However, in all of these models, the mediation scale must be extremely low ($<100$ TeV). We conclude by considering the prospects for the high-luminosity LHC era, where we expect the current limits on particle masses to improve by up to $\sim 1$ TeV, and discuss further model-building directions for natural SUSY that are motivated by this work.

# Introduction and Summary

Recently, at the ICHEP 2016 Conference , the ATLAS and CMS collaborations presented the results of the search for new physics using the first $\sim 15$ fb$^{-1}$ at 13 TeV . This represents a significant milestone for the LHC: with this dataset the sensitivity to new physics at the energy frontier begins to truly surpass that achieved at Run I. Now is therefore the perfect time to assess the implications of these results for well-motivated models of new physics such as supersymmetry (SUSY).

Weak-scale SUSY has long occupied a central place in the theoretical expectations for the LHC, as the addition of superpartners to the Standard Model (SM) particles at or near the scale of electroweak symmetry breaking stabilizes the Higgs mass and solves the hierarchy problem (for a review and original references, see e.g. ). Given that the superpartners must be heavier than their Standard Model counterparts, the supersymmetric cancellation of loops protecting the Higgs mass cannot be perfect. The heavier the superpartners, the more finely tuned the original bare mass must be against the loop contributions. If we require the theory be fully "natural" – that the level of fine-tuning be less than some fixed amount, taken in this work to be the arbitrary threshold of 10% – then an upper limit can be derived on the mass of the Standard Model SUSY partners . For a recent review on naturalness in SUSY with many original references, see e.g. .

While specific realizations of SUSY can have a wide variety of predictions for the spectrum of the superpartner masses, basic requirements of naturalness which hold in a wide class of models include a light higgsino (which directly sets the Higgs mass squared parameter at tree-level), relatively light stop squarks (as the top has an ${\cal O}(1)$ Yukawa, leading to large one-loop-corrections to the Higgs), and a relatively light gluino (which corrects the mass of the stop itself, yielding a two-loop correction to the Higgs) .

As the higgsino is color-neutral, it is difficult to produce directly and detect at the LHC. However, if the gluino and stop are relatively light as required by naturalness, they will be copiously produced at the LHC. As they cascade decay down to the light higgsinos, they will typically yield at least one or more of the following signatures :

-   Significant missing transverse momentum (MET)

-   Top quarks

-   High object multiplicity

The purpose of this paper is to investigate the parameter space of natural SUSY that is still allowed by the latest LHC searches. So far, no search has turned up definitive evidence for new physics, and the null results are phrased in terms of limits on a small set of "simplified models". In order to carry these limits over to more general scenarios, e.g. those motivated by natural SUSY, a detailed reinterpretation ("recasting") of the LHC searches is required. In this work, we have performed a comprehensive recasting of the 13 TeV post-ICHEP analyses which are most relevant for natural SUSY (i.e. target the signatures listed above), see Table  for a complete list. (Previous work that reinterpreted recent 13 TeV LHC results includes .) For more details and the validation of our simulations by comparison with the official ATLAS and CMS limits, see Section  and Appendix .

Using the recasted searches, we will explore the parameter space of natural SUSY, using a carefully chosen set of representative simplified models. Our philosophy here will be similar to that of : we work purely bottom up, motivated to find the most conservative limits on gluinos, stops and higgsinos. In all the models we consider, the lightest MSSM sparticle is the higgsino, with $\mu\le 300$ GeV as suggested by naturalness. However, apart from this assumption, we allow ourselves a great deal freedom in the simplified models.[^1]

First, as an essential baseline model, we will consider "vanilla SUSY" – the minimal supersymmetric Standard Model (MSSM) with $R$-parity conservation and flavor-degenerate sfermions. Next we will examine simplified models of natural SUSY that alleviate the latest LHC constraints by either reducing the signal cross section or by reducing the signal acceptance. We will consider the "effective SUSY" scenario where the 1st/2nd generation squarks are decoupled. Decoupling the valence squarks in particular reduces the total SUSY cross section by factors of ${\cal O}(10)$ or more in the region of the mass plane near the current LHC limits. We also consider two scenarios that trade MET for jets: baryonic $R$-parity violating (RPV) decays of the higgsino (see e.g.  for a review and original references); and a hidden-valley (HV) scenario inspired by Stealth SUSY . By trading MET for jets, the signal becomes more difficult to distinguish from QCD, thereby significantly degrading the acceptance.

By examining the LHC limits on these simplified models, we will attempt to draw more general conclusions on viable directions for SUSY models post-ICHEP. For each simplified model, we will overlay the current LHC limits in the gluino/stop mass plane together with the $\Delta\leq10$ natural regions, for different choices of the higgsino mass $\mu$ and the messenger scale $\Lambda$. Here $\Delta$ is derived from the Barbieri-Giudice fine-tuning measure with respect to a soft SUSY-breaking parameter $M$, reformulated in terms of $m_h$ instead of $m_Z$ in order to better take into account the effect of radiative corrections to the Higgs quartic: $$\Delta_{M^2} = {\partial \log m_h^2\over \partial \log M^2}$$ When multiple sources of tuning are present, we take the maximum tuning as our measure, $\Delta=\max_{\{M_i\}}\Delta_{M_i^2}$.

For the calculation of $\Delta_{M^2}$, it has been conventional in much of the literature to work in the leading-log (LL) approximation (see however for notable exceptions). There the quadratic sensitivity of the Higgs mass-squared parameter to the higgsino, stop and gluino soft masses arises at tree level, one-loop and two-loops respectively:

-   Higgsinos: $$\label{deltamhusqhiggsino}
    \delta m_{H}^2 = \mu^2$$

-   Stops: $$\label{deltamhusqstop}
    \delta m_{H}^2\sim -{3\over8\pi^2}y_t^2 m_{\rm stop}^2 \log {\Lambda\over Q}$$

-    Gluinos: $$\label{deltamhusqgluino}
    \delta m_{H}^2\sim -{g_3^2y_t^2\over 4\pi^4}|M_3|^2\left(\log{\Lambda\over Q}\right)^2$$

Here $\Lambda$ is the messenger scale of SUSY breaking, and $Q$ is the IR scale, conventionally taken to be 1 TeV in many works (see e.g. ). Aside from the naturalness bounds this yields on the higgsino, stop and gluino masses, these LL formulas also demonstrate that the tuning is worsened as the messenger scale is raised. Natural SUSY theories greatly prefer lower values of $\Lambda$.

In this work, we go beyond the leading-log approximation and include a number of important higher-order effects, including the full two-loop RGEs, one and two-loop threshold corrections to stop and gluino masses and threshold corrections to the Higgs potential. A detailed description of these effects (and original references) will appear in a companion paper . Here we will summarize the main idea: we translate the tuning bounds on the UV mass parameters $M^2$ which enter into the Barbieri-Giudice measure into upper limits on pole masses at the IR scale. The physical pole masses are what the LHC sets limits on. Surprisingly, the full set of differences between UV and IR parameters in tuning calculations has largely been neglected in the literature so far, but we find that they have several crucial consequences.

First, they are numerically important and they can raise the tuning bounds on sparticles by $\mathcal{O}(1)$ factors. Second, including these higher-order corrections makes $\Delta_{M_3^2}$ ($\Delta_{m_{Q_3}^2}$) dependent on the stop (gluino) mass. Large gluino masses significantly raise the stop IR mass through the RG, while large stop masses contribute non-negligible threshold corrections to the gluino pole mass. Therefore in a natural spectrum, the stop and gluino mass are actually correlated: a heavy stop implies a heavy gluino, and vice versa. Perhaps it is counterintuitive, but the fact that we have seen neither the stop nor the gluino may be *more* consistent with natural SUSY than the discovery of one and not the other. In any event, this means that the $\Delta\leq10$ natural region is not simply a rectangle in the gluino/stop mass plane, but instead turns out to be wedge-shaped. Finally, the higher-order corrections include effects from the 1st/2nd generations that become very important in effective SUSY scenarios. At one-loop, heavy 1st/2nd generation squarks appreciably lift the gluino pole mass, which helps to relax the tuning bound on the gluino. At two-loops, the RGEs drive the stop mass lower in the IR, which can strengthen the tuning bound on the stops. See for more details.

We will see that vanilla SUSY is strongly constrained by the current searches and cannot be natural at the 10% level for any choice of the messenger scale – this was true already after Run I. With either Effective SUSY or models that trade MET for jets, the LHC limits are greatly reduced, but still eliminate most of the parameter space with $\Delta\leq10$, except at the very lowest messenger scales $\Lambda\lesssim 20$ TeV. Finally, we consider the combination of Effective SUSY with RPV, and show that significant natural parameter space still remains at $\Lambda\lesssim$ 100 TeV. As many of models put forward as alternative solutions to the Hierarchy Problem *start* with at least 10% tuning , the continued survival of natural SUSY serves as a reminder that – despite the lack of discovery – supersymmetry remains one of the least tuned solutions for physics beyond the Standard Model.[^2]

Throughout this work, we will neglect the role of the gravitino, and in particular the possibility that the higgsino LSP decays to it within the detector. This is compatible with the low mediation scale $\Lambda\lesssim 100$ TeV provided that the effective SUSY breaking felt by the messengers was much smaller than the ultimate SUSY-breaking scale $\sqrt{F}$ in the hidden sector. (This is the parameter called $k$ in .) An assumption along these lines is also necessary to make the baryonic RPV scenario compatible with the low messenger scale; otherwise with very light gravitinos, new proton decay channels such as $p\to\psi_{3/2}K^+$ arise, and proton stability bounds ($\lambda''_{ijk}<(10^{-6}-10^{-15})(m_{3/2}/1$ eV), depending on flavor indices ) would preclude the possibility of hiding SUSY by trading MET for RPV jets.

We conclude by showing rough estimates of the LHC reach for sparticles through 300 fb$^{-1}$ and 3 ab$^{-1}$, using an extrapolation based on the method of . As will be seen, we are situated at approximately the middle of the rapid rise in the superpartner reach due to the increased LHC energy; after reaching approximately 50 fb$^{-1}$ of data (tentatively expected for Moriond 2017), the LHC will have spent most of its energy boost and additional coverage will be slower and more incremental. Overall, we expect the asymptotic improvement in the LHC reach to be an across-the-board increase of $\sim 900-1200$ GeV to the current limits, largely independent of the SUSY particle mass. These projections imply that the high-luminosity LHC (HL-LHC) can exclude or discover all models of fully-natural SUSY that we consider in this work.

These extrapolations assume no qualitatively new analysis techniques will be developed, so there could be room for even greater future improvement. In particular, the recasted limits we have derived on some of our simplified models come from searches that were generally not designed with the phenomenology of the natural SUSY in mind. Perhaps by targeting natural SUSY (and specifically the simplified models considered here), ATLAS and CMS could significantly improve their sensitivity. For example, the RPV and multi-jet searches are generally optimized for gluino pair production, but we find them to be also relevant for constraining stop pair production. Since the latter involves a different set of physics objects and object multiplicity, with a greatly reduced cross section, perhaps a re-optimization would better maximize $S/\sqrt{B}$ and further extend the reach in this case. For these reasons, we encourage the LHC collaborations to adopt some or all of our simplified models for natural SUSY, for the purpose of optimizing searches and setting official limits.[^3] The question of naturalness is one of the prime motivations for new physics at the LHC, and as such it is too important to be left solely to amateurs such as ourselves!

Our work suggests many future directions that will be interesting and important to explore. On the collider phenomenology side, given the ability of the LHC searches to powerfully constrain conventional natural SUSY models, as demonstrated here, it will be important to examine models with more exotic final states (see for a recent brief review). Possibilities may include displaced decays (for the signatures we study, these were already well constrained by 8 TeV data , while displaced leptons are less tested ); and collimated particles that fail isolation such as "dirty leptons" and "lepton jets" . Also, we have ignored tau leptons in this work, assuming that they will be at least as stringently constrained by either jets+MET or lepton-based searches. It might be worthwhile to test this assumption more rigorously, e.g. in the case of displaced taus where there are important gaps in coverage .

On the model building side, there are several, well-known extensions of the MSSM (e.g. Dirac gauginos) which provide loopholes to the tuning bounds. These are not considered here, but they are increasingly well-motivated. We will discuss them further in Section . Even within the more conventional context of natural SUSY with light higgsinos, stops and gluinos, there are many interesting model building directions to pursue. As discussed above, we will show that $\Delta\leq10$ requires a very low messenger scale ($\Lambda\lesssim 100$ TeV), and this is an important constraint of future models of natural SUSY. This is especially true for models of Effective SUSY, since these constructions generally tie SUSY breaking to the generation of flavor in the SM. This would mean that flavor must also be generated at an extremely low scale, and it is not at all obvious that this is viable. Some examples of previous attempts include Refs. . Also, we find that adding a HV/Stealth sector to the MSSM can trade MET for jets and greatly reduce the bounds. It is interesting to speculate whether this additional sector could be used for anything else, such as dark matter or raising the Higgs mass. More generally, obtaining a 125 GeV Higgs is a major issue for natural SUSY and requires going beyond the MSSM, and it is interesting to think about whether extensions of the MSSM which succeed in raising the Higgs mass could also help to hide SUSY at the LHC.

This paper is organized as follows. In Section , we describe our methodology for reinterpreting LHC searches (with additional information provided in Appendix ). The models of natural SUSY we consider are described in Section , and the resulting experimental limits on these models in Section . We conclude in Section  with projections for the future reach of the LHC and model-building directions suggested by the existing constraints.

# Recasted searches and methodology

In the following sections, we will consider the status of natural SUSY after the most recent results from the 13 TeV LHC, as mostly reported in the ICHEP 2016 conference . These results, using $12-18$ fb$^{-1}$ of data from CMS and ATLAS, greatly extend the experimental reach of the LHC for gluinos and squarks. We concentrate on the searches listed in Table , each of which has many signal regions (SRs) that target specific mass spectra and supersymmetric production modes.[^4]

| Search | Data (fb$^{-1}$) | Reference |
|:--:|:--:|:--:|
| ATLAS 2-6 jets + MET | 13.3 | [`ATLAS-CONF-2016-078`](https://cds.cern.ch/record/2206252/) |
| ATLAS 8-10 jets + MET | 18.2 | [`ATLAS-CONF-2016-095`](https://cds.cern.ch/record/2212161/) |
| ATLAS $b$-jets+MET | 14.8 | [`ATLAS-CONF-2016-052`](https://cds.cern.ch/record/2206134/) |
| CMS jets + MET ($H_T^{\rm miss}$) | 12.9 | [`CMS-PAS-SUS-16-014`](https://cds.cern.ch/record/2205158/) |
| ATLAS SSL/3L | 13.2 | [`ATLAS-CONF-2016-037`](https://cds.cern.ch/record/2205745/) |
| ATLAS 1L+jets+MET | 14.8 | [`ATLAS-CONF-2016-054`](https://cds.cern.ch/record/2206136/) |
| ATLAS multi-jets (RPV) | 14.8 | [`ATLAS-CONF-2016-057`](https://cds.cern.ch/record/2206149/) |
| ATLAS lepton+many jets | 14.8 | [`ATLAS-CONF-2016-094`](https://cds.cern.ch/record/2211457/) |

Searches most important to our study. All use the 13 TeV LHC data.

As in , we did not recast searches with photons or two or more opposite-sign leptons, under the assumption that any natural SUSY scenario yielding these signatures would be even more constrained than the simplified models we have considered here.

As can be seen from Table , we primarily use ATLAS searches. In most cases, the CMS searches have so many signal regions (100+) that they are difficult to reinterpret. A proper approach would require sophisticated statistical methods combining multiple exclusive bins, using information (the correlation matrix of errors) that is not publicly available. In contrast, the ATLAS searches explicitly provide 95% CL limits on number of events due to new physics for each signal region, and generally have far fewer, coarser bins, allowing us to simply use the most sensitive SR to set a conservative (but reasonably accurate) limit.

One important case where CMS did include aggregate signal region information is the jets+MET search with $H_T^{\rm miss}$ , which we find to be very powerful. The CMS jets+MET search has $b$-tagged categories, low-MET and high-MET categories, few jet and many jet categories. As a result, it is equivalent to the union of several different ATLAS jets+MET searches.[^5]

Our simulation methodology is as follows. We generate hard events using MadGraph 220.127.116.11 , generating pairs of gluinos, squarks, and antisquarks in all possible combinations. Decay and showering is performed via Pythia 8.219 (with which we implement the RPV and stealth decays of the higgsino), and we use Delphes 3.3.2 with an ATLAS-approximating detector geometry for detector simulation. In order to speed up the event generation, we worked with unmatched samples in most cases.[^6] The NLO cross sections for the superparticle production are obtained via Prospino . The cutflows for the ATLAS and CMS analyses are recasted using Root 5.34 , and validated against the published experimental limits on supersymmetric simplified models.

The details of the simulated analyses, along with the results of these validation checks, are shown in Appendix . As seen there, while our simulation technique is largely successful in matching the experimental results, sometimes there are slight differences due to the crudeness of our recasting framework. Sources of error could include: the use of Delphes instead of a complete Geant4  detector simulation; the use of the ATLAS detector geometry for CMS searches; or the use of the 70% $b$-tagging efficiency working point for all searches, while some of the experimental searches use the 77% working point. (This latter choice significantly increases the $c$-quark mistagging rate: thus, the use of lower efficiency working point will not always reduce our limits when compared to the official searches, especially if the search involves $b$-vetoes.) Another possible issue is that each of our validation plots (based on a simplified SUSY model) has limits set by a small number of signal regions, while our general models might be sensitive to different signal regions, which are thus never explicitly validated against the experimental efficiencies.

In any event, as can be seen in Appendix , the addition of a 50% "recasting uncertainty" (i.e. multiplying or dividing the signal efficiencies by a factor of 1.5) around the baseline results is sufficient to bring our recast in line with the official limits, in almost every case. For the models of natural SUSY analyzed below, we will show exclusion limits in which we have taken the lower end of our recasting uncertainty. That is, we consider a point in the mass parameter space excluded if it exceeds the observed limit from any of the SRs for each search in Table , after we have lowered our efficiencies by a factor of 1.5. We use this conservative estimate to make sure that we do not falsely exclude parameter space because of possible inaccuracies in our recasting procedure.

# Overview of simplified models

Here, we will give a brief overview of the various simplified models that we will use in the next section to illustrate the current status of natural SUSY. Although our simplified models are based on the MSSM for the most part (the HV/Stealth scenario is one exception), we expect they are representative of a much broader class of natural SUSY models, and therefore our qualitative conclusions should be much more general.

In all the models we consider, the lightest supersymmetric particle in the MSSM spectrum is the higgsino, which for reasons of naturalness should be lighter than 300 GeV. (We do not consider alternative models where the higgsino mass arises not from $\mu$ but from some SUSY-breaking operator . We will discuss this possibility further in Section .) At present, the best limits on direct higgsino production are set by LEP, with $\mu \geq 95$ GeV for a stable LSP and $\mu\geq103$ GeV for RPV decays of the higgsino into $udd$ quarks . (As discussed in the Introduction, we assume that the gravitino is such that the higgsino does not decay to it within the detector.) The recent CMS search with opposite-sign leptons and MET does not yet set limits on direct higgsino production, but bodes well for the future in finally overtaking LEP limits. We will take two benchmark values of the higgsino mass, either $\mu=300$ GeV (the maximum allowing $\Delta=10$) or $\mu=100$ GeV (at the LEP limit).

One might wonder if a stop or gluino LSP (possibly NLSP with a gravitino LSP) is allowed in any part of the parameter space: if the colored partner is (meta-)stable, searches for $R$-hadrons set limits at 1.6 and 1 TeV for gluinos and stops, respectively. For a colored NLSP with prompt decays to a gravitino LSP, the usual simplified topologies with a massless neutralino apply (e.g. $\tilde{t}\to t\psi_{3/2}$) and both stops and gluinos limits are well above 300 GeV. If the decays violate $R$-parity, stop and gluino LSPs are excluded by pair-produced multi-jet resonances, up to 400 GeV and 800 GeV , respectively. Similarly, the presence of additional light neutralinos and charginos (as might be expected alongside the higgsino LSP) will only add to the decay chains of colored spartners, with more final states more easily picked up by the searches considered here. Therefore, we do not believe that we have introduced significant blind spots by assuming a natural higgsino LSP.

For each choice of $\mu$, we consider the following models:

1.  *"Vanilla SUSY"*: the MSSM with $R$-parity conservation (RPC) and all three generations of squarks degenerate.[^7]

2.  *"Effective SUSY"*: The RPC MSSM with all squarks other than the stops and the left-handed sbottom decoupled from the mass spectrum (see Refs. ). For the light third-generation fields, their (pole) masses are taken to be the same, but our results will still hold if they are of the same order.

    The key advantage this model has over the vanilla MSSM in evading the LHC constraints is the much reduced cross section for colored SUSY sparticle production, as the 1st and 2nd generation squarks are heavy. With a Majorana gluino, the $t$-channel gluino-mediated valence-squark cross sections are enormous, assuming the squarks are kinematically accessible. In Fig. , we show the ratio of the total colored SUSY production cross section in the Effective SUSY scenario over the vanilla SUSY case. In the region of most interest for natural SUSY, when both the gluinos and squarks masses are $\sim 1-2$ TeV, the cross section for colored pair production is reduced by more than a factor of 10. Obviously, this significantly reduces the experimental reach for these models.

    As we briefly discussed in the Introduction (see for details), decoupled 1st/2nd generation squarks tighten (loosen) the naturalness bound on the stops (gluino), because they lower the stop IR masses through the two-loop RGEs (raise the gluino masses through the one-loop pole mass corrections). In the "Effective SUSY" models we consider here, we set 1st/2nd generation squarks to what amounts to a "sweet spot" at 5 TeV: too heavy to be efficiently produced at the LHC, heavy enough to help with the gluino fine tuning, but light enough that they do not reduce the natural stop mass range by much.

3.  *RPV SUSY*: The MSSM with baryonic $R$-parity violation (RPV) in which the lightest higgsino can itself decay promptly to three quarks (see for a review). By turning the MET from the higgsino LSP into more jets via RPV decays, the SUSY signal can be hidden (to some degree) in the larger QCD background. Although there are many possible flavor combinations for $UDD$ RPV, we will focus on the $cds$ operator as the one which is relatively safe from precision constraints , yet will relax experimental limits the most. We note that the limits from the will largely carry over to other flavor combinations as well (except for top quarks).[^8]

4.  *HV/Stealth SUSY*: an R-parity conserving "hidden valley" (HV) extension of the MSSM in which the higgsino can further decay into new gauge-singlet scalar $S$ and its fermionic superpartner $\tilde S$: $$\tilde{H} \to S \tilde{S},~ S \to gg,$$ The scalar $S$ decays with 100% branching ratio into pairs of gluons, while the fermion is stable due to RPC. This model trades MET for jets, provided that $m_S+m_{\tilde S}\approx m_{\tilde H}$ and $m_{\tilde S}\approx 0$. (We take $m_S=m_{\tilde H}-10$ GeV, and $m_{\tilde{S}}=5$ GeV.) This simplified model is also a proxy for a number of different well-motivated scenarios. For instance, it could also arise from GMSB with higgsino NLSP, $\tilde H\to h+\tilde G$, with $m_{\tilde H}\approx m_h$ . It can also be thought of as a "lite" version of Stealth SUSY ; embedding our particle spectrum into an actual Stealth SUSY construction would provide a natural explanation of the required mass degeneracy, while presumably not modifying the limits significantly.

5.  *Effective SUSY with RPV*: a combined Effective-RPV SUSY scenario with first and second generation squarks decoupled and a higgsino LSP decaying to three jets via baryonic RPV. Not surprisingly, this scenario is the least constrained by current searches.

In the models with an unstable higgsino, we assume the conventional cascade decays until reaching the higgsino LSP: for example, the chargino (nearly-mass degenerate with the neutral higgsino) decays to the LSP via an off-shell $W$ and is not allowed to decay directly to jets.

# Results

Finally, we are ready to explore the implications of the latest LHC null results for natural SUSY. We will use our recasting framework to calculate the limits on the simplified models described in the previous section, and then overlay the $\Delta\le 10$ "fully natural" region (as determined using the precision calculations in ) over these limits to see what range of natural gluino and stop masses are still allowed. As discussed in the Introduction and described in detail in , the fully-natural $\Delta\le 10$ region is wedge shaped, because the physical stop mass now depends sensitively on the gluino mass through the RGEs, while the physical gluino mass depends to a lesser extent on the stop masses through finite thresholds. For other choices of $\Delta$, the limits on the masses scale approximately as $\sqrt \Delta$. The extent of these regions depends on the messenger scale $\Lambda$, with lower $\Lambda$ leading to larger allowed masses, see the LL formulas ()-() for the rough, qualitative intuition. In our plots, we show the tuning wedges for $\Lambda = 20$ and 100 TeV.[^9]

Our results from applying the reinterpreted LHC searches to the natural SUSY models, scanning over gluino and stop masses, are summarized in Fig.  for vanilla SUSY, Fig. , for Effective SUSY, Fig.  for both RPV SUSY and HV/Stealth SUSY, and Fig.  for Effective RPV SUSY. Of the eight available experimental analyses in Table , we only show the most constraining for each model, to avoid cluttering the plots. A dashed line indicates the combined nominal limits that we find in absence of any recasting uncertainty; solid lines indicate the limits with the aforementioned conservative reduction in signal efficiency by a factor of 1.5. For reference, we also include the appropriate limit on each model from the 8 TeV data, using the recasting framework developed in Refs. .

## Vanilla SUSY

In Fig.  we see that the parameter space for natural vanilla SUSY (i.e. the MSSM with light higgsinos, gluinos, and flavor-degenerate squarks) with $\Delta \le 10$ is completely excluded by the LHC results, even for $\Lambda = 20$ TeV. In fact, this was basically true even at 8 TeV. Not surprisingly, the most powerful analyses for constraining this scenario are the general-purpose jets+MET searches. The combination of the large production cross section of gluinos and squarks and the large missing momentum carried away by the LSP due to the short decay chain makes vanilla SUSY a rather easy target for these searches.

Though not shown in Fig. , we find the current 13 TeV limits correspond to $\Delta=20$ (i.e. 5% tuning) with a low messenger scale, or $\Delta=200$ with $\Lambda=10^{16}$ GeV (this is discussed further in Section ). Increasing $\Lambda$ only reduces the natural region for the gluino and squark masses, as expected from the LL formulas ()-().

Some of the limits exhibit a sharp discontinuity along $m_{\tilde g}=m_{\tilde q}$, becoming much weaker below the diagonal. Above the diagonal, where $m_{\tilde q}>m_{\tilde g}$, the gluino decays dominantly to Higgsinos via off-shell stops and sbottoms, so the signal is top and bottom rich. Below the diagonal, where $m_{\tilde q}<m_{\tilde g}$, the gluino decays to all flavors of squarks with nearly equal branching ratio. This reduces the jet multiplicity and the number of $b$-jets on average, significantly weakening the limits from the searches that require these signatures.

## Effective SUSY

We show the limits on the Effective SUSY model (with 1st/2nd generation squarks decoupled to 5 TeV) in Fig. : gluinos below 1.8 TeV are excluded while limits on direct stop production are at 900 GeV, significantly raising the previous 8 TeV limits.

Again, we see that the strongest constraints are set by searches targeting large MET, in this case the ATLAS $b$'s+MET and the CMS jets+MET searches. While we did not reinterpret the many dedicated stop searches from ATLAS and CMS for this study, as can be seen, the general SUSY searches are very powerful, excluding stops nearly up to 1 TeV. Indeed, from the CMS official summary plot , one sees that the general purpose CMS jets+MET is nearly as effective as the dedicated stop searches in constraining the basic $\tilde t\to t+\chi_1^0$ simplified model. So we expect that including the dedicated stop searches would not qualitatively change the conclusions here.

Despite these strong limits, there remains a viable (albeit small) range of natural gluino and stop masses in Effective SUSY, but only for extremely low values of $\Lambda$. While $\Lambda =20$ TeV is not yet ruled out, $\Lambda = 100$ TeV is already excluded. Evidently, reducing the SUSY cross section by a factor of $\sim 10$ (see Fig. ) by decoupling the first and second generation squarks is not enough to completely relax the constraints from the latest round of searches.

## RPV and HV/Stealth SUSY

We now turn to SUSY models which trade MET for jets. Obviously, these models are going to be far less constrained by the standard MET-based searches. However, searches which target large multiplicities of high-$p_T$ jets instead of MET (such as the ATLAS 8-10 jets search and the ATLAS RPV search ), are still very powerful. In these scenarios, we have included one additional jet at the generator level (and matched the matrix-element and parton-shower calculations in the MLM scheme ): for squarks, the hard process would have resulted in 8 final partons, and adding an extra parton raises the reach of the ATLAS RPV and ATLAS 8-10 jets searches by approximately 100 GeV.

In Fig , top row, we show the limits on RPV SUSY, allowing the higgsino to decay into a trio of $cds$ quarks (the results would be similar for any $u_{i\neq3}d_jd_k$ operator, in particular, for final states with $b$ quarks the ATLAS RPV limits would increase, while ATLAS 8-10 jets would stay the same as it does not involve b-tagging). As can be seen, while the natural masses are excluded for $\Lambda = 100$ TeV, a small region of the $\Lambda = 20$ TeV gluino and squark mass range remains unexplored, assuming the higgsino mass is 100 GeV. If this mass is raised to 300 GeV, the jets resulting from the RPV decay are more effectively captured by the high-multiplicity searches, and the entire $\Delta < 10$ space is excluded. With a lighter higgsino, the quarks in the final states are more collimated and result in significantly fewer resolved jets, which is more difficult to distinguish from the QCD background. This important characteristic was discussed at length in .

Again, as in the vanilla SUSY case, there is a pronounced shift in the strength of the limits across the $m_{\tilde g}=m_{\tilde q}$ diagonal, because the gluino goes from dominantly decaying to Higgsinos via off-shell stops and sbottoms (above), to decaying to all flavors of squarks equally (below). Above the diagonal, where the gluino decays to Higgsino are top-rich, the ATLAS SS dilepton search sets an equally strong limit as the high-multiplicity searches.[^10]

The same features are also seen in the HV/Stealth SUSY results, bottom row of Fig. . As in the RPV case, $\Lambda = 100$ TeV is already ruled out for HV/Stealth, and only lower values of the messenger scale remain viable. The fact that the limits on the RPV and the HV/Stealth scenarios are so quantitatively similar, despite the scenarios having different kinematics and different number of jets in the final state, is evidence that the LHC limits are fairly robust, and that the simplified models we have chosen are representative of a broader class of scenarios that trade MET for jets.

## RPV Effective SUSY

Finally, having seen the (limited) success of both RPV/HV scenarios and Effective SUSY in evading the LHC bounds on natural SUSY, we consider their combination. In Fig. , we show the experimental reach for models of Effective RPV SUSY, where the higgsino is unstable and only the two stop squarks and the left-handed sbottom are light, while the remaining sbottom and the first and second generation squarks are decoupled at 5 TeV. As expected, the limits are the weakest of all the models considered so far, with limits on gluinos and third-generation squarks respectively at 1.4–1.5 TeV and 600–800 GeV. (As in the previous subsection, the same limits would apply to any $u_{i\neq3}d_jd_k$ final state.) Even here, the $\Delta\leq10$ parameter space for a 100 TeV messenger scale is nearly excluded, but much of the $\Lambda = 20$ TeV parameter space remains viable.

## Summary of results and further implications

In the previous sections we have excluded a wide range of gluino, squark and stop masses for a variety of natural SUSY models, and understood the implications for fine-tuning. In Fig. , we further apply the calculations of in order to show the minimum amount of tuning $\Delta$ compatible with a given messenger scale $\Lambda$, for each of the natural SUSY models we consider in this paper. (Qualitatively, these curves can be understood/extrapolated from the results shown in the previous subsections, using the LL formulas ()-().) As can be seen, even with our most optimistic scenario (Effective SUSY with RPV decay of higgsinos), the scale $\Lambda$ must be less than 100 TeV for $\Delta\leq10$. It should be noted that other choices of "acceptable" levels of fine-tuning allow higher messenger scales. For example every scenario we have considered (except perhaps vanilla SUSY) is only tuned at the percent-level or better, even with messengers at the GUT scale.

Aside from naturalness considerations, the individual recasted limits on each superpartner are noteworthy as they cannot always be obtained from the ATLAS and CMS summary plots (this is particularly true for the RPV/HV/Stealth cases, where the ATLAS RPV and 8-10 jets searches do not consider squark simplified models). For this reason, in Table  we summarize the asymptotic limits on each colored superpartner (gluinos and either mass degenerate squarks or third-generation squarks); these limits are obtained from the plots above by decoupling either the gluinos or the squarks.[^11]

For the models considered so far, we have restricted ourselves to $\mu\leq300$ GeV, due to the tree-level higgsino contribution to fine-tuning. An obvious question is how the limits scale with higher higgsino masses. We find that most of the limits do not change significantly for higgsino masses below $\sim 500$ GeV, but above $\sim 500$ GeV, the MET-based scenarios such as vanilla SUSY and effective SUSY start to see some degradation due to compression, with third generation squarks affected first (near 450 GeV), then degenerate squarks (above 600 GeV), and finally gluinos (above 800 GeV). On the other hand, for models with unstable neutralinos, the limits actually *increase* moderately (by $\sim 100-200$ GeV) with higher higgsino masses, as the jets neutralino decays are less boosted. These features are fully consistent with the neutralino mass dependence of the validation plots shown in the Appendix.

Before we conclude, a comment is required: in both RPV and HV/Stealth models, it can be seen that fully hadronic MET-based searches set limits that are competitive with searches which do not rely on missing energy. A close inspection at individual events reveals that the source of missing energy is neutrinos in the presence of semi-leptonic $W$ decays (for example, from top quarks or $B$ mesons), where the concurring lepton is not isolated and is therefore removed. The resulting events have non-zero missing energy and no isolated leptons, and are accepted by the all-hadronic searches (in particular, CMS jets and ATLAS 8-10 jets, which have moderate MET requirements). While we have closely mirrored the overlap removal procedures in the experimental papers, our cruder simulation framework might be overestimating the reach of these searches. In any case, the more robust limits from other searches (particularly ATLAS RPV and 1L+multijets) in Figs.  and result in similar exclusions, with differences of at most 50 GeV.

```latex
\resizebox{\columnwidth}{!}{
\begin{tabular}{|c|c|c|c|c|c|c|c|c|c|c|c|} \hline
Model & {Vanilla SUSY} & \multicolumn{2}{c|}{Effective SUSY} & \multicolumn{2}{c|}{RPV SUSY} & \multicolumn{2}{c|}{Stealth/HV SUSY} & \multicolumn{2}{c|}{RPV Effective SUSY}
\\\hline
$\mu$ [GeV] & 100 -- 300 &100  & 300 &100 &300 &100 &300 &100 &300
\\\hline
$m_{\tilde{g}}$ [GeV] & 1730 & \ 1690\, & 1690 & 1310 & 1500 & \ \ 1330 \ \,  & 1440 & \ \ 1350\ \, & 1490
\\
$m_{\tilde{q},\tilde{t}}$ [GeV] & 1500 & 975 & 950 & 700 & 810 & 600 & 750 & 550 & 750
\\\hline
\end{tabular}
\caption{Observed lower limits on the mass of the gluino $m_{\tilde{g}}$ and squarks $m_{\tilde{q}},m_{\tilde{t}}$ for each model considered above. The limits are asymptotic in the sense that they refer to the case where all other superpartners decouple. The vanilla SUSY limits are independent on the higgsino mass for $\mu<300$ GeV.
 \label{tab:results}}
}
```

# Conclusions and Future Directions

## Projections to 300 fb$^{-1}$ and 3 ab$^{-1}$

We have found in the previous section that while vanilla SUSY has not been fully natural already since Run I, many flavors of SUSY beyond vanilla are currently still viable with $\Delta\le 10$. Here, we show naive extrapolations of our limits to 300 fb$^{-1}$ and 3 ab$^{-1}$ of integrated luminosity (assuming either 13 and 14 TeV collisions), corresponding approximately to the end of the regular LHC operations and the ultimate end of the HL-LHC runs, respectively.

Our methodology is the same as (we thank A. Weiler for some clarifications in this regard). We assume that the signal cross section is controlled by a single mass scale $m$, and that signal *efficiencies* and background *counts* remain constant as the integrated luminosity ${\mathcal L}$ and/or center-of-mass (CM) energy $\sqrt{s}$ are increased. The latter makes sense, provided that the mass reach grows in a way that one can always cut harder to keep the backgrounds low, while preserving the sensitivity to new physics at higher and higher mass scales. It also assumes that as ${\mathcal L}$ and/or $\sqrt{s}$ are increased, there are no qualitatively new obstacles that cannot be overcome with more clever analysis techniques. For example, our projections ignore the effect of pile-up: while a challenge for the HL-LHC, it should prove less of a barrier to searches in the mass range above 1 TeV, where SUSY decays typically result in many high-$p_T$ jets and/or large MET.

Using this approach, if the previous limit at CM energy $\sqrt{s_1}$ and integrated luminosity ${\mathcal L_1}$ was at $m_1$, the extrapolated limit to CM energy $\sqrt{s_2}$ and integrated luminosity ${\mathcal L_2}$ can be obtained by requiring $$m_1^{-2} f(m_1/\sqrt{s_1}) {\mathcal L_1} = m_2^{-2} f(m_2/\sqrt{s_2}) {\mathcal L_2}$$ where $f$ is the parton luminosity (taken here to be $gg$ for simplicity – the projections for $q\bar q$ initiated production are actually slightly stronger, but qualitatively similar).

Interestingly, under these assumptions, the ultimate improvement in mass reach, compared to current limits, is nearly a constant shift, $m \to m+\Delta m$, across a wide range of masses (1 TeV $\lesssim m \lesssim$ 3 TeV). This can be traced back to the fact that, in this range of $x=m/\sqrt{s}$ values, the parton luminosities at the LHC happen to be dropping nearly exponentially, $f(x)\sim e^{-a x}$. In Fig. , we show the increase in mass reach benchmarked against a limit set with 12 fb$^{-1}$ of 13 TeV, as a function of integrated luminosity (for both 13 and 14 TeV LHC running). After 300 fb$^{-1}$ of data is collected at 13 TeV, the rough rule of thumb is that we expect the limits on both gluinos and stops to improve by $\Delta M\sim 500$ GeV from their current levels, with a further improvement to $\Delta M \sim 900-1200$ GeV after 3 ab$^{-1}$ at 14 TeV. For example, taking the current best limits on gluino (1.8 TeV) and stop (900 GeV) masses, this would correspond to an ultimate reach of $\approx 3$ TeV and $\approx 2$ TeV respectively. Using these projections, we see that the $\Delta\le 10$ regions for all models we consider can be fully explored by the end of the HL-LHC era.

As an aside, while very light particles do not see the full shift in mass reach, the asymptotic improvement would still have a very significant impact, fractionally speaking. This is especially relevant for non-colored particles, where the current limits in many cases are not far beyond the LEP bounds. Here an increase of $\sim 200-300$ GeV to the range of experimental sensitivity to these channels would constitute at least a doubling of the existing bounds, if not more. This provides a strong motivation for the construction of the HL-LHC.

## Future directions for model building

One of the principal values of a recasting work on simplified models is that it motivates further model building to populate the allowed parameter space, which could in turn lead to new correlated signatures to explore (or could suggest that the allowed parameter space cannot be realized in UV complete models). So we will conclude by describing some of the model building avenues that are suggested or motivated by our work. This is not an exhaustive list, but highlights some directions of theoretical and experimental interest.

The first (and least interesting) possibility in light of these results is simply to relax the requirement of naturalness. The weak scale may involve a slight numerical accident; after all there are plenty of other percent-level accidents in Nature. This is a logical possibility (perhaps further motivated by the 125 GeV Higgs mass), but it may not be experimentally testable. Tuning at the level of $\Delta \sim 100$ is not well constrained at the LHC, and the relevant parameter space is unlikely to be excluded in most SUSY models in the foreseeable future.

The remaining alternatives all require extending SUSY beyond the vanilla scenario. For instance, in this paper we have explored two possibilities. The first, effective SUSY, lowers the total SUSY cross section by removing the possibility of valence squark production. As emphasized in our companion paper , without the addition of further particles as in , the squarks cannot be set to arbitrarily high masses without risking tachyonic 3rd generation sparticles. Decoupling some of the squarks also raises questions about the resolution of the supersymmetric flavor problem. All of these difficulties are compounded by the requirement of very low messenger scales ($\Lambda\lesssim 100$ TeV). It remains to be seen whether any viable model of effective SUSY exists with such low messenger scales. Some promising prior work that explored models in this direction includes .

The second avenue we explored is RPV/HV/Stealth SUSY, which reduces the signal to background ratio by trading MET for jets. It would be interesting to investigate the HV/Stealth sector further, to see whether it could be useful for anything else, e.g. dark matter or lifting the Higgs mass. Also, further exploration of RPV scenarios is well-motivated. For instance, a potential option to further reduce the constraints on RPV discussed in this work might be having at least one large coupling, such that squarks decay dominantly to jets instead of higgsinos (searches for paired dijet resonances are then relevant ). Although constraints from flavor physics can be important, this can be accomodated in models with hierarchies dictated by flavor symmetries . At the same time, this also opens the possibility of resonant squark production via the $udd$ operator, with additional signatures and 8 TeV limits recently described in . Nevertheless, a low messenger scale calls for cautious model-building if new flavor-violating interactions are to be present.

As noted previously, in this work we have not considered the possibility of decays to light gravitinos inside the detector. In the high-MET scenarios (vanilla SUSY, effective SUSY), allowing the higgsino to decay to gravitino plus $h/\gamma/Z$ is unlikely to qualitatively change our conclusions on a stable higgsino, with the increase in final state multiplicity possibly reinforcing our limits. One exception would be the compressed scenario $m_{\tilde \chi^0_1}\simeq m_h$, where some of this phenomenology would be covered by our Hidden Valley model, as already noted above. In any event, for what concerns experimental limits, it would be interesting to see if opening a gravitino decay channel would significantly change the limits presented in this work.

Finally, another interesting direction would be to challenge the underlying assumptions going into the tuning calculations: for example the SUSY production cross section can be reduced by introducing "super-safe" Dirac gluinos (which eliminate valence squark diagrams with $t$-channel gluinos). Alternatively, models where the higgsino mass is not set primarily by the $\mu$ term, e.g.  are increasingly well-motivated. A higgsino LSP above 600 GeV can lead to compressed spectra and greatly weakened limits, as discussed in Section ; for instance the gluino could be as light as $\sim 800$ GeV given the current bounds on the $\tilde g\to qq\chi_1^0$ simplified model. Also, if the higgsino was not the LSP, it could lead to weakened bounds if the gluino can decay directly to a HV/Stealth sector .

In this work, we have explored simplified models motivated by natural SUSY where cascade decays of accessible gluinos and stops down to a light higgsino produce collections of ordinary jets, leptons, etc. We have considered scenarios with and without MET. While limits from current searches may still allow fully natural SUSY models, it is a testament to the breadth of the ATLAS and CMS experimental programs that much of the $\Delta\leq10$ parameter space is fully excluded already, while what remains (at $\Lambda\lesssim 100$ TeV) can be completely covered by the end of the HL-LHC run.

#### Note Added

The first preprint version of this paper differed with respect to the present published version in two aspects: first, it did not include matching, which led to an underestimation of some limits, especially for the ATLAS RPV search constrains on unstable neutralino scenarios of Sec. . Second, a bug was found in Delphes 3.3.2 which caused an overestimation of the ATLAS 8-10 jets search limits on the same RPV and Hidden Valley scenarios.[^12] The neutralino, which had a significant boost originating from the large mass splitting with the squarks or gluinos, decays to boosted quarks and in the ensuing parton shower a small (but significant enough) fraction of events had long-lived unstable hadrons, with a nominal lifetime large enough to escape the detector. While at the LHC those would deposit all their energy in the hadronic calorimeter and then decay, Delphes would throw away their daughter particles as it used only the list of final stable particles within the detector volume for each event. Given these missing high-momentum particles in some events, the resulting $E_T^{\rm miss}/\sqrt{H_T}$ spectrum was not as steeply falling as it should be, and the requirement $E_T^{\rm  miss}/\sqrt{H_T}>4\,{\rm GeV}^{1/2}$ of Ref.  was satisfied in a larger fraction of events, resulting in stronger limits from that search. Fixing this bug reduced the ATLAS 8-10 jets limits by $200-400$ GeV, especially for the squarks. The net effect of these two changes is negligible in most of the parameter space of Fig.  with no differences in the combined limits in the natural regions of parameter space and of at most 100 GeV elsewhere.

# Acknowledgements

We are grateful to C. Campagnari, T. Lari, F. Moortgat and I. Vivarelli for their help in clarifying many aspects of the ATLAS and CMS SUSY searches. We also thank N. Arkani Hamed, N. Craig, K. Cranmer, J. A. Evans, Y. Kats, M. Papucci, M. Reece, J. Ruderman, M. Strassler, S. Thomas and A. Weiler for helpful discussions. Finally, we are grateful to J. A. Evans, M. Farina, Y. Kats, M. Reece and M. Strassler for comments on the draft. The work of AM and DS are supported by DOE grant DE-SC0013678. The work of SM is supported by DOE grant DE-SC0010008.

# Recasting Details and Validation

Here, we discuss the details of our reinterpretation of each of the ATLAS and CMS experimental searches used in this paper. Each search has provided details of the cut-flow for the various signal regions, along with the expected and observed number of events, allowing us to apply the results to the supersymmetric models of interest. In order to validate our implementation of each search, we apply our recast to the simplified supersymmetric models used by the experiments themselves. These simplified models typically have only two or three supersymmetric particles kinematically accessible. As such, they are not appropriate for our study of naturalness in SUSY.

As in the study of the full supersymmetric theories, described in Section , we generate hard events for the simplified models using MadGraph5 , with NLO cross sections from Prospino ; we have not noticed appreciable differences in the validation plots between generating matched and unmatched samples. We decay and shower these events with Pythia8 . Detector simulation is via Delphes3 , which makes use of the FastJet package for jet finding. We use the same jet-clustering algorithms as in each experimental search, namely the longitudinally invariant $k_t$ algorithm and the anti-$k_t$ algorithm , as well as jet trimming on reclustered jets . In Delphes, we reconstruct the missing energy vector as the negative sum of all calorimeter deposits and all muons. In Delphes a ParticleFlow algorithm is also present, which combines the tracker and the calorimeter information to define physics objects. The differences in the missing energy are usually small, but we choose the former method as it results in validation plots that are in slightly better agreement with the official results, although the differences are within our "theory error" estimate.

In most cases, recasting the searches requires only a straight-forward implementation of the cut-flow described in the relevant conference note. Where necessary, we note any deviations we were forced to take from the search as reported by the experimental collaboration. In the following, we show our validation plots, comparing the experimental exclusion region with the exclusion region we find on simplified supersymmetric models. As in the original experimental searches, limits are set by taking observed limits on the number of signal events in the signal region which has the best expected limits.

## ATLAS Same-sign Lepton/Three Lepton

The ATLAS note CONF-2016-037 is a search for supersymmetric particles decaying to jets and leptons, requiring either two same sign leptons or three leptons in the final state. A number of signal regions are defined, separated by number of $b$-tagged jets, lepton multiplicity and missing transverse momentum.

We validate our recasting of the search by generating events using one of the supersymmetric simplified models considered by : gluino pair production decaying to top pairs and a neutralino $\tilde{g} \to t\bar{t}\tilde{\chi}_1^0$. The published limits are shown in Fig.  with the results of our recasted search on simulated data, along with a 50% "recasting uncertainty" on the number of events in each signal region.

## ATLAS Lepton Plus Jets

The ATLAS note CONF-2016-054 is a search for gluinos and squarks decaying through $W^\pm$ bosons (via charginos), requiring one lepton in the final state, along with jets and missing transverse momentum. A number of signal regions are defined, separated by number of jets, $b$-tagged jets, and missing transverse momentum.

We validate our recasting of the search by generating events using the supersymmetric simplified model considered by : gluino pair production decaying to light-flavor quarks and a chargino, which itself decays to a neutralino and a $W$ boson: $\tilde{g} \to q\bar{q}'\tilde{\chi}_1^\pm,~\tilde{\chi}^\pm_1 \to W^\pm \tilde{\chi}^0_1$. The chargino mass is set to be the average of the gluino and neutralino mass. The published limits are shown in Fig.  with the results of our recasted search on simulated data, along with a 50% "recasting uncertainty" on the number of events in each signal region.

## ATLAS Multi-$b$

The ATLAS note CONF-2016-052 is a search for gluinos decaying to third generation quarks (tops or bottoms) and missing transverse momentum. At least three $b$-jets must be identified in the final state. Some signal regions further require "fat" jets which have topological similarities to top quarks. We followed the procedure outlined in by re-clustering the $\Delta R =0.4$ jets into jets of radius $0.8$ using the anti-$k_T$ algorithm in Delphes, and then further trimming the resulting jets by removing subjets whose $p_T$ falls below $10\%$ of the $p_T$ of the re-clustered jet.

We validate our recasting of the search by generating events using the supersymmetric simplified model considered by : $(i)$ gluino pair production decaying to top pairs and a neutralino $\tilde{g} \to t\bar{t}\tilde{\chi}_1^0$ and $(ii)$ gluino pair production decaying to bottom pairs and a neutralino $\tilde{g} \to b\bar{b}\tilde{\chi}_1^0$. The published limits are shown in Fig.  with the results of our recasted search on simulated data, along with a 50% "recasting uncertainty" on the number of events in each signal region.

## ATLAS RPV

The ATLAS note CONF-2016-057 is a search for RPV SUSY. A number of signal regions are identified with varying number of jets, $b$-tagged jets, and large radius jets. These "fat" jets are simulated in Delphes by reclustering the calorimeter jets into jets of radius $1.0$ using the anti-$k_T$ algorithm. Then, the resulting large jets are trimmed by re-clustering their components using the $k_T$ algorithm with a sub-jet radius parameter of $0.2$ and discarding sub-jets carrying less than $5\%$ of the original large jet. The surviving sub-jets are used to calculate the "fat" jet energy and momentum, which is then further corrected by the jet energy scale.

We validate our recasting of the search by generating events using the supersymmetric simplified model considered by : $(i)$ gluino pair production decaying to all quark pairs and a neutralino $\tilde{g} \to q\bar{q}\tilde{\chi}_1^0$ followed by neutralino decay via RPV operators into three quarks $\tilde{\chi}_1^0 \to qqq$ (with equal branching ratios to all available flavor combinations), and $(ii)$ gluino pair production decaying directly to three quarks $\tilde{g} \to qqq$. The published limits are shown in Fig.  with the results of our recasted search on simulated data, along with a 50% "recasting uncertainty" on the number of events in each signal region.

## ATLAS 2–6 Jets Plus MET

The ATLAS note CONF-2016-078 is a search for gluinos and squarks decaying to jets and missing energy, requiring between two and six jets, significant missing energy, and vetoing on leptons. Two search strategies are employed in : one using a $m_{\rm eff}$ variable to separate signal and background, and a second using RJR variables . For simplicity, we use the former search signal regions, which sets bounds as competitive as the latter.

We validate our recasting of the search by generating events using the supersymmetric simplified models considered by : $(i)$ gluino pair production decaying to light-flavor quarks and a neutralino $\tilde{g} \to q\bar{q}\tilde{\chi}_1^0$, $(ii)$ light-flavor squark pair production decaying to quarks and a neutralino $\tilde{q} \to q\tilde{\chi}_1^0$, $(iii)$ and gluino pair production decaying to light-flavor quarks and a chargino, which itself decays to a neutralino and a $W$ boson: $\tilde{g} \to q\bar{q}'\tilde{\chi}_1^\pm,~\tilde{\chi}^\pm_1 \to W^\pm \tilde{\chi}_0^1$. In the latter case, the chargino mass is set to be the average of the gluino and neutralino mass. The published limits are shown in Fig.  with the results of our recasted search on simulated data, along with a 50% "recasting uncertainty" on the number of events in each signal region.

## CMS Multi-Jet + MET

There are three searches from the CMS Collaboration which search for gluinos and squarks decaying to jets and missing energy which were presented at ICHEP in 2016 . We chose to work with (CMS-SUS-16-014), which has equivalent reach as the other two searches. This search requires at least three jets, no leptons, and significant missing transverse momentum. The full analysis uses 160 signal regions, separated by minimum jet, $b$-tagged jet, $H_T$, and $H_T^{\rm miss}$. However, this large number of signal regions can be simplified to twelve aggregated regions. In each region, we calculate the maximum number of signal events which can be accommodated at 95% CL from the published background expectation and observation (Appendix A.5 of ) using the CL$_{\rm s}$ method.

We validate our recasting of the search by generating events using the supersymmetric simplified models considered by : $(i)$ gluino pair production decaying to light-flavor quarks and a neutralino $\tilde{g} \to q\bar{q}\tilde{\chi}_1^0$, $(ii)$ light-flavor squark pair production decaying to quarks and a neutralino $\tilde{q} \to q\tilde{\chi}_1^0$, $(iii)$ gluino pair production decaying to top pairs and a neutralino $\tilde{g} \to t\bar{t}\tilde{\chi}_1^0$, and $(iv)$ gluino pair production decaying to bottom pairs and a neutralino $\tilde{g} \to b\bar{b}\tilde{\chi}_1^0$. The published limits are shown in Fig.  with the results of our recasted search on simulated data, along with a 50% "recasting uncertainty" on the number of events in each signal region.

## ATLAS 8–10 Jets Plus MET

The ATLAS note CONF-2016-095 (which is the direct update of the 7–10 jet search in with an increased luminosity of 18.2 fb$^{-1}$) is a search for gluinos decaying to jets and missing energy, requiring between eight and ten jets, some missing energy, and vetoing on leptons. "Fat" jets are used to discriminate over the background, in addition to the ratio $E_T^{\rm miss}/\sqrt{H_T}$. A number of signal regions are defined, with varying jet multiplicity and requirements on the sum of the fat jet masses. These large-$R$ jets are found by reclustering the small-$R$ jets with the anti-$k_t$ algorithm and a radius $R=1.0$ in Delphes. Then, the sum of the masses of the reclustered jets is used to define the signal regions.

We validate our recasting of the search by generating events using the supersymmetric simplified models considered by , $(i)$ gluino pair production decaying to light-flavor quarks and a chargino, which itself decays to a neutralino and a $W$ boson, $\tilde{g} \to q\bar{q}'\tilde{\chi}_1^\pm,~\tilde{\chi}^\pm_1 \to W^\pm \tilde{\chi}_0^1$, with the chargino mass set to be the average of the gluino and neutralino mass, and $(ii)$ gluino pair production decaying to light-flavor quarks and a chargino, which then decays to a $W$ boson and a neutralino $\tilde \chi^0_2$, followed by the neutralino decay to a $Z$ boson and the lightest neutralino $\tilde \chi^0_1$, $\tilde{g} \to q\bar{q}'\tilde{\chi}_1^\pm,~\tilde{\chi}^\pm_1 \to W^\pm \tilde{\chi}_0^2, \tilde{\chi}_0^2\to Z \tilde{\chi}_0^1$. The published limits are shown in Fig.  with the results of our recasted search on simulated data, along with a 50% "recasting uncertainty" on the number of events in each signal region.

## ATLAS Lepton Plus Many Jets

The ATLAS note CONF-2016-094 is a search for gluinos decaying to top-rich final states and little missing transverse momentum, requiring one lepton and multiple jets in the final state. A number of signal regions are defined, separated by number of jets and $b$-tagged jets.

We validate our recasting of the search by generating events using the supersymmetric simplified model considered by : gluino pair production decaying to top quarks and a neutralino, which itself decays via RPV to three light quarks: $\tilde{g} \to t\bar{t}'\tilde{\chi}_1^0,~\tilde{\chi}^0_1 \to uds$. The published limits are shown in Fig.  with the results of our recasted search on simulated data, along with a 50% "recasting uncertainty" on the number of events in each signal region.

[^1]: Since this paper is focused on the implications of the latest LHC direct searches for natural SUSY, we will not require our simplified models to raise the Higgs mass to 125 GeV, instead remaining agnostic as to the source of this mass. As is well understood, if the SUSY Higgs mass corrections arise only from the MSSM stop squark loop, the level of tuning is at the few percent level or worse (see e.g.  for a recent detailed discussion and references).

[^2]: We note that there have been papers in the literature, even during Run I, that have claimed SUSY is at least percent-level tuned in all circumstances. Obviously, given the content of this paper, we believe these claims were premature. A clear point of reference is with the work of : although we consider a rather similar set of simplified models, we come to completely different conclusions about the tuning. A more in-depth comparison reveals the sources of the discrepancy. Aside from using a different measure of fine-tuning (summing in quadrature and multiplying tunings, vs. taking the max of the EW tuning), the main difference is that in , broad conclusions about fine-tuning in SUSY were drawn based on the consideration of a small, limited set of more UV-complete models whose messenger scales never go below $\Lambda=300$ TeV. Whereas in this paper, we consider a broad range of messenger scales down to $\Lambda=20$ TeV, with no attempt at model building. Since we find that only models with $\Lambda<100$ TeV can be better than 10% fine-tuned after ICHEP, there is in fact no contradiction with the work of .

[^3]: In our work we have conservatively adopted a 50% "theory uncertainty" on our recasting efficiencies, based on our validations against official limit plots. If ATLAS and CMS made official limits on our simplified models, that alone might "improve" the limits derived in this work by up to 50% in the effective cross section.

[^4]: We also considered the ATLAS 7-10 jets+MET search with 3.2/fb and the CMS black hole search with 2.2/fb. Due to the strong possibility of control-region contamination, the latter necessitated a conservative reinterpretation along the lines of . Neither search set the strongest limit in any of the simplified models considered in this paper, so they are not included here. However, an update of the CMS BH search to the full $\sim 15$ fb$^{-1}$ dataset (especially if optimized for gluino production as in ) would likely have competitive sensitivity to high-multiplicity, non-MET simplified models such as the RPV-based scenarios considered here.

[^5]: Other CMS general-purpose searches provide an equivalent reach on simplified models studied but provide significantly less information. has no aggregate signal regions and does not even provide observed event counts in the text! does have aggregate signal regions, but adds considerable computational complexity in calculating the clustered jets used in the $M_{T2}$ variable, and was beyond the scope of this work.

[^6]: Because of the choice of models studied in this paper (in particular, given the high number of jets in the unmatched events), matrix/element/parton-shower matching is usually not necessary: the only case where including extra jets at the generator level makes a difference in the experimental acceptance is for the RPV/Stealth models. We refer to Section  for further comments.

[^7]: Note that we set both left- and right-handed squarks at the same scale. A variant of this model would have right-handed down-type squarks decoupled: this does not affect the fine-tuning (which is affected by both chiralities of stops, but not by the right-handed sbottom) but reduces the SUSY cross sections by removing the production of $\tilde d_R, \tilde s_R$. In this case, the limits on squarks are reduced by about 100 GeV with respect to the Vanilla SUSY model (shown in Fig. ).

[^8]: Baryonic RPV scenarios in which the $td_id_j$ operator dominates (e.g. motivated by flavor symmetries ) result in top-rich final states which either give large-radius jets or leptons. We have checked this case is better constrained: the limits are raised by about 200 GeV for both gluinos and squarks/stops.

[^9]: While $\Lambda=10^{16}$ GeV is well-motivated, the $\Delta\le 10$ tuning wedge for this case is so tiny that it would barely show up in the plots, so we do not bother to show it. In fact, for $\Lambda=10^{16}$ GeV, the $\Delta\le 100$ tuning region happens to be qualitatively similar to the $\Lambda=20$ TeV, $\Delta\le 10$ region, see for details.

[^10]: This may come as a slight surprise, as the SS dilepton searches were found to be not as effective for constraining natural SUSY at 8 TeV . The difference is that in the new search, there is a new signal region (SR1b-GG) that does not require any MET. This further highlights the power and importance of doing SUSY searches with low or no MET.

[^11]: We note that, for mass-degenerate squarks, "decoupling" the gluino in practice means taking $m_{\tilde g}\gtrsim10$ TeV, otherwise the gluino $t$-channel contribution to valence squark production remains non-negligible.

[^12]: For more details and the bug-fix appearing in Delphes 3.4.1, see <https://cp3.irmp.ucl.ac.be/projects/delphes/ticket/1084>.