\section{Computational Analysis}
\label{sec:analysis}
A look at equations (\ref{FOC_d}), (\ref{FOC_a}) and (\ref{NeutralOpt}) and the corresponding expressions for elasticities in the Appendix suggests that a meaningful theoretical analysis does not look feasible. We, therefore, do a computational analysis of $[OPT1]$, $[OPT2]$ and $[OPT3]$. Recall that our objective is to determine whether equilibrium welfare in one regime is higher than in the other.

\begin{figure*}[t!]
\begin{center}
\subfloat[]{\includegraphics[width=0.33\linewidth]{Figures/Welfare_bar_v_eq_r}\label{fig:Wel_v_eq_r}}
\subfloat[]{\includegraphics[width=0.33\linewidth]{Figures/Consumer_density_v_eq_r}\label{fig:Con_v_eq_r}}
\subfloat[]{\includegraphics[width=0.325\linewidth]{Figures/CP_density_v_eq_r}\label{fig:CP_v_eq_r}}\\
\subfloat[]{\includegraphics[width=0.33\linewidth]{Figures/Welfare_bar_v_gr_r}\label{fig:Wel_v_gr_r}}
\subfloat[]{\includegraphics[width=0.33\linewidth]{Figures/Consumer_density_v_gr_r}\label{fig:Con_v_gr_r}}
\subfloat[]{\includegraphics[width=0.325\linewidth]{Figures/CP_density_v_gr_r}\label{fig:CP_v_gr_r}}\\
\caption{\small{(a) Equilibrium Welfare (b) Consumer Market Size and (c) Provider Market Size for $v=r=10$, and (d) Equilibrium Welfare (e) Consumer Market Size and (f) Provider Market Size for $v=10$ and $r=2$, for the cases of Non-neutral Equilibrium, Neutral Equilibrium and Welfare Optimum. Other parameters used for these results are $t = 0.5, f = 0.25,\text{ and }\mu = 3$.}}
\label{fig:VandR}
\end{center}
\end{figure*}

Extensive computational analysis reveals that we get an unambiguous welfare ranking of the non-neutral and neutral equilibrium for a large range of parameter values. Specifically, \emph{the non-neutral equilibrium welfare-dominates the neutral equilibrium and is close to the welfare optimum}. Figure~\ref{fig:VandR} demonstrates this using example results for when the consumer valuation $v$ is equal to the advertising surplus $r$ (Figures~\ref{fig:Wel_v_eq_r}-\ref{fig:CP_v_eq_r}) and when $v > r$ (Figures~\ref{fig:Wel_v_gr_r}-\ref{fig:CP_v_gr_r}). Qualitatively, the results (not shown) are similar for when $v < r$.

Figures~\ref{fig:Con_v_eq_r} and~\ref{fig:Con_v_gr_r} compare the consumer market size $\hat x$ under the neutral equilibrium, non-neutral equilibrium, and the welfare optimum. For a fixed $\lambda$, we see that the consumer market size under the non-neutral equilibrium is larger than in the neutral regime. Also, it is similar to the size under the welfare optimum. Further note from Figures~\ref{fig:CP_v_eq_r} and~\ref{fig:CP_v_gr_r} that the provider market size $n$ under the neutral regime is significantly larger than under the welfare optimum and non-neutral regimes. Since both the consumer and the provider market size in the non-neutral equilibrium are closer to the welfare optimum, the non-neutral equilibrium dominates the neutral equilibrium in terms of welfare, which can be seen in Figures~\ref{fig:Wel_v_eq_r} and~\ref{fig:Wel_v_gr_r}.

Lastly, note that, in line with expectation, as $\lambda$ increases, the resulting increase in congestion at the ISP causes the size of the consumer market to steadily decrease under all three regimes. A similar trend with respect to $\lambda$ is seen for the provider market size.

We summarize the takeaway from computational analysis in the following proposition.
\begin{proposition}  \thlabel{WelfareCompare}
	[Welfare Comparison of Equilibria]. For a large selection of parameter values, the non-neutral regime is welfare superior to neutral regime.
\end{proposition}

The intuitive reason why the result holds is because a switch from a non-neutral regime to a neutral regime leads to welfare losses on both counts -- a decrease in the consumer market size that it entails decreases the economic value generated by consumers (net of any gains achieved from slightly reduced congestion costs); and an explosive increase in the provider market size (way above the welfare optimum) imposes costs on society that far outweigh any benefits to consumers from increased content diversity (number of providers).

\section{Appendix}
Let us define
\begin{align*}
A &= \mu + \lambda \big(v-d - \frac{t (a+f)}{2r} \big) \\
B &= \sqrt{\big[\mu - \lambda \big(v-d - \frac{t (a+f)}{2r} \big) \big]^2 + 4 \lambda}
\end{align*}
Then, we can compute
\begin{align*}
e(\hat{x}(d, a), d) &= \frac{\lambda d \Big[ \mu - \lambda \big(v-d - \frac{t (a+f)}{2r} \big) + B \Big]}{B(A-B)} \\
e(\hat{x}(d, a), a) &= \frac{\lambda ta \Big[ \mu - \lambda \big(v-d - \frac{t (a+f)}{2r} \big) + B \Big]}{2r B(A-B)}
\end{align*}
Equations (\ref{FOC_d}) and (\ref{FOC_a}) can now be reduced to the following pair of equations in $(d, a)$.
\begin{align}
e(\hat{x}(d, a), d) &= \frac{da+df}{da+df+ra}   \label{eq1} \\
e(\hat{x}(d, a), a) &= \frac{f}{a+f}   \label{eq2}
\end{align}
For the neutral regime, we have~(\ref{eq3}).
Equation (\ref{eq3}) can now be used in Equation (\ref{NeutralOpt}) to solve for $d$.
\section{Conclusions and Future Directions}
\label{sec:conclusion}
This paper melds a competition model from the industrial organization literature with a network congestion model, in which equilibrium analysis of non-neutral and neutral regimes as well as a standard economic welfare analysis of equilibrium is carried out. We restrict ourselves to the case of a monopolistic ISP and a competitive provider market. Our results indicate a surprisingly clean welfare analysis. While there are both costs and benefits to increasing consumer and provider market size, the neutral regime is welfare dominated by the non-neutral regime. It remains to be seen whether this conclusion changes when the ISP's investment incentives are taken into account.

We focused on the market for a particular content type  with symmetric providers with the same fixed cost. However, the provider landscape is very rich, populated both by independent monopolists (big players like Google, Amazon, Facebook etc.) streaming different content types and competitors in each content type (for instance Amazon, Flipkart, Alibaba etc. in the e-commerce landscape). It is not immediately clear what is a good way to model the richness of the provider landscape. Another important extension of our work would be to the case of duopoly ISPs as is done, for instance, in \cite{economides2012network} and \cite{armstrong2006competition}.

One difficulty with modeling is that while the neutral regime is relatively clearly defined, the non-neutral regime throws up several possibilities for the ISP. Apart from the feasibility of pricing access to the CPs, the ISP can also introduce lane prioritization and differentially price the lanes. Analyzing all the richness of possibilities in a single model seems elusive and analysis of a battery of models may be our best bet.
\section{Introduction}
\label{sec:introduction}
Network neutrality prohibits an ISP from discriminating between packets of similar type of content from different providers by either pricing their transit differently or by providing them different qualities of service. As has been well articulated by~\cite{coucheney2013impact}, ISP(s) that are not directly connected to CPs and provide a transit to the content generated by the CPs have argued that CPs are cornering good amounts of dividends on an investment that was made by the ISPs. While ad-revenues are rising for the CPs, the transit prices that are a source of revenue for such ISPs are decreasing. The ISPs would thus like the CPs to pay for transit of their content or suffer loss in quality-of-service.

On the other hand, the proponents of net neutrality point to the detrimental effect that a non-neutral regime can have on the CPs and by implication on the consumers. For example, as pointed out in~\cite{economides2012network}, a non-neutral regime could impact innovation at the edge of the network, as it could now be influenced by the ISPs that constitute the core of the network. Allowing discriminatory pricing or provisioning of quality-of-service may enable an ISP to adversely impact access of content by consumers of some among competing CPs.

In this work, we will analyze net neutrality from the point of view of a regulator whose mandate is to maximize the sum welfare of all consumers, CPs, and the ISPs. Note that the larger the number of CPs providing similar content, the higher is the intensity of competition amongst providers. More competition must increase consumer welfare. However, too much competition amongst providers may squeeze their revenues so much that they end up making negative returns on the cost of their investment to set up business. Our regulator can't merely maximize consumer welfare but must contend with the impact of too much competition on the welfare of CPs. A natural line of inquiry is then to include such aspects in any welfare analysis of net neutrality regulation. For example, models used for analysis must capture competition in the provider market. However, we are not aware of any previous work that looks into this direction.

By nature, the ISP operates in a two-sided market. It sells network access to both consumers and CPs. Consumers can then choose to consume content from providers. Our model is therefore a two-sided market model in the style of~\cite{economides2012network}. In this paper, we restrict ourselves to analyzing the case of a monopolist ISP.  We model the interaction of the ISP, the CPs and the consumers as a multi-stage game. The ISP moves first to set access prices for both consumers and CPs. The CPs move next making a binary choice of whether to enter the market. Finally, consumers choose whether to connect to the ISP. This structured interaction via a multi stage game permits an analysis of the influence of the pricing decisions of the ISP on the size of the CP market and the size of the consumer market. Another feature of our modeling is that CPs provide all their content free to consumers\footnote{This is motivated by the providers' business model and agrees with the literature.}; so their only source of revenue is advertising.

In our model, all CPs provide similar content. However, differences in consumer tastes imply that even at the same price a consumer may prefer a certain CP over another. That is the content providers are not perfect substitutes and their content offerings are horizontally differentiated~\cite{tirole1988theory}. This guides how the providers arrange themselves among the consumers of the ISP. Similarly, consumers have non-price related preferences with regards to the ISP. Additionally, we make the reasonable assumption that it is more difficult for a consumer to find an alternate ISP than it is to find an alternate CP.

To capture the aforementioned we use two workhorse models -- namely the linear city model due to \emph{Hotelling}~\cite{hotelling1929} and the circular city model due to \emph{Salop}~\cite{salop1979monopolistic} -- from the industrial organization literature~\cite{tirole1988theory}. Note that a consumer consumes the two products of network access sold by the ISP and content sold by a CP. Technological constraints imply that a consumer must first have network access to subsequently consume content. We, therefore, face two product dimensions -- one for the ISPs and one for the CPs. To capture these, we use a novel way to model our market, which we will call the \emph{Hotelling-Salop} model, where in we employ the \emph{Hotelling} line for the ISP dimension that determines how many consumers connect to the ISP and then let the resulting segment of that line to be the circumference of the \emph{Salop} circle on which CPs align to provide content to the consumers.

In our model, the ISP has a fixed amount of bandwidth that it must use to service requests of content made by its consumers. Thus one would expect the quality-of-service (average delay in this paper) obtainable by any given consumer to be impacted by the demands for content made by the other consumers. To capture this fact, as in~\cite{pil2010net}, we use the M/M/1 queueing model, which assumes that the aggregate consumer demand at the ISP is modeled as a Poisson process and that the time it takes for the ISP to service a consumer's request for content is well modeled by an exponential distribution. The resulting average delay is treated as a cost borne by every consumer. While the M/M/1 model is simplistic, it allows us to capture the important fact that average delay (cost) suffered by a consumer increases -- slowly, when the aggregate demand is small and doesn't congest the available bandwidth, and rapidly as the aggregate demand makes the network congested -- as the number of consumers of the ISP increase.

Finally, in our model, as in~\cite{economides2012network}, network neutrality is defined as the constraint that the ISP cannot charge the CPs.

\thref{WelfareCompare} in Section~\ref{sec:Welfare} summarizes our main result, which is that the non-neutral equilibrium dominates the net neutral equilibrium. As is borne out by our analysis, this is because welfare depends on the two variables of consumer and provider market size and not on the prices charged of the consumers and providers. The prices do not matter for welfare as they end up being mere transfers from one player to another. A switch from a non-neutral to a neutral regime leads to welfare losses on both counts -- a decrease in the consumer market size that it entails decreases the economic value generated by consumers (net of any gains achieved from reduced congestion costs); and an explosive increase in the provider market size (way above the welfare optimum) imposes costs that far outweigh any benefits to consumers from increased diversity of content providers.

The rest of this paper is organized as follows. Section \ref{sec:related} provides a limited overview of the growing literature on net neutrality. In Section~\ref{sec:spatialModels}, we provide a brief on the use of spatial models to model competition and consumer preferences. In Section \ref{sec:model}, we set up the model and pursue the equilibrium analysis under non-neutral and neutral regimes. We then go on to show how to do a welfare analysis of equilibria and also how to compute a welfare optimum. Computational analysis is carried out in Section \ref{sec:analysis}. Lastly, we make a few concluding remarks in Section \ref{sec:conclusion}.
\section{Basic Model}
\label{sec:model}
There is a monopoly ISP and many potential CPs for whom this ISP is not the originating ISP. Introduce the interval $[0,1]$ as the \emph{Hotelling} line of the standard horizontal product differentiation model for the ISPs. Let the ISP be located at the left extremity, location $0$, of this line. There is a continuum of consumers uniformly distributed on the Hotelling line. Let the product differentiation parameter for the ISP be normalized to $1$. Non-price factors therefore impose a cost of $x$ on a customer at location $x$ on the Hotelling line.

Let $t \in (0,1)$ be the product differentiation parameter for the CPs. Note that we restrict $t<1$. This is to capture the fact that it is easier to substitute between CPs than it is to substitute between ISPs. Let $y$ be the distance travelled by a consumer to connect to a CP. This imposes a cost of $ty$ on the consumer. Consumers cannot connect to the CPs without first connecting to the ISP. The CPs space themselves equidistantly on a Salop circle whose perimeter, as we detail soon, is determined by the customers that choose to connect to the ISP.

All CPs provide free content to the consumers; their only source of revenues is from advertising. Let $r$ be the advertising revenue per consumer, same for all CPs. Let $d$ be the access price charged by the ISP to the consumer. Let $f$ be the fixed cost of entering the provider market for any CP.

Consumers pay a congestion cost, which is the average time the ISP takes to fulfill a consumer's request for content. To capture this, we model the ISP as a M/M/1~\cite{gross1998fundamentals} service facility. Specifically, following~\cite{pil2010net}, we assume that the aggregate demand for content from consumers created at the ISP is well modeled by a Poisson process of rate $\lambda$ per unit length of the single dimension containing consumers connected to the ISP. Further, the time taken by the ISP to service a consumer's request for content is exponentially distributed with a mean service time of $1/\mu$. We assume $\lambda < \mu$, which is required for the average time to fulfill a request to stay bounded. Finally, consumer requests are serviced by the ISP in a first-come-first-served manner. This means that a request that arrives at the ISP must wait for requests that arrived earlier to finish obtaining service from the ISP.

Consumer monetary valuation of the total requested content is $v$. In the non-neutral regime, the ISP can charge the CPs for servicing requests for their content at an access price of $a$. In contrast, the neutrality regime prohibits the ISP from charging the CPs a price when it is not the originating ISP for them, which we assume is the case here. The two regimes therefore induce different market games through which the entities interact.

\subsection{Non-neutral Regime}

\begin{figure}
	\centering
	\includegraphics[width=0.9\columnwidth]{Figures/Fig_1.pdf}
	\caption{\small{Left Panel shows the \emph{Hotelling} line for the ISP dimension. Right panel shows the folded \emph{Hotelling} line segment of length $\hat{x}$ as the \emph{Salop} circle with the same circumference. It also shows the symmetric alignment of CPs when $n=4$ for when the customer at $\hat{x}$ travels a distance $\hat{x}/(2n) = \hat{x}/8$ to reach a CP.}}
	\label{fig:ModFig}
\end{figure}

In the non-neutral regime, the ISP, the CPs and the consumers interact through the following \emph{multi-stage game}:
\begin{enumerate
\item ISP chooses the access prices: $a$ for CPs and $d$ for consumers.
\item Consumers then choose whether to connect to the ISP.
\item CPs who are considering entering the market then simultaneously decide whether to enter or not. Suppose $n$ is the number of CPs who chose to enter. The CPs locate equidistantly\footnote{We follow \cite{salop1979monopolistic} in making the assumption of symmetric alignment of CPs on the Salop circle to focus on the entry decisions of CPs.} on a \emph{Salop} circle whose perimeter is the length $\hat{x}$ of the initial segment of the \emph{Hotelling} line on which consumers connect to the ISP.
\item Consumers who have connected to the ISP then connect to a CP and consume content.
\end{enumerate}

We specify utility of a consumer who travels a distance $x$ on the \emph{Hotelling} line to connect to the ISP and distance $y$ on the \emph{Salop} circle to connect to the CP as
\begin{equation}\label{eqn:consumerUtil}
u = v - x  - ty - d - \frac{1}{\mu- \lambda \hat{x}},
\end{equation}
where $\lambda \hat{x}$ is the aggregate rate of request for content from the length $\hat x$ of consumers. The last term in the expression for consumer utility $u$ is the congestion cost faced by a consumer. It is the resulting expected time the ISP takes to fulfill a consumer's request (called waiting time in~\cite{gross1998fundamentals}), given the demand $\lambda \hat{x}$.

We specify ISP's profit as $\pi(d, a) = d \hat{x} + na$. We specify profits of $\text{CP}_i$ as $\pi_i = r D_i - a - f$, where $D_i$ is the demand that $\text{CP}_i$ faces.

To determine this demand $D_i$ we must decide how the CPs are aligned on the Salop circle relative to the marginal consumer at $\hat{x}$ who connects to the ISP. This is important for it determines the distance traveled by this consumer on the Salop circle and therefore affects the location $\hat{x}$ on the Hotelling line. Two extreme cases are possible. In one, illustrated in Figure~\ref{fig:ModFig}, the marginal consumer travels a distance $\frac{\hat{x}}{2n}$; in the other, the consumer travels no distance, that is a CP locates itself at $\hat{x}$. %
One could introduce an extra parameter $\beta \in [0,1]$, and specify the distance traveled by the marginal consumer as $\frac{\beta \hat{x}}{2n}$. The extreme cases correspond to the values $\beta = 1$ and $\beta = 0$. We will, however, choose $\beta = 1$ in our analysis for two reasons. First, any $\beta \in (0, 1]$ retains the structural feature that demand on either side (whether $\hat{x}$ or $n$) responds to prices ($d$ and $a$) on both sides. Therefore $\beta = 1$ carries the flavor of analysis of any $\beta \in (0,1)$. Second, as will be clearer from analysis that follows, $\beta = 0$, kills any dependence of $\hat{x}$ on $a$ and in that sense is a special case.

To summarize, the parameters of the model are the consumer valuation $v$, the advertising revenue $r$, the normalized product differentiation parameter $t$, fixed CP costs $f$, the rate of content request $\lambda$, and the service rate $\mu$. The choice variables are the prices $a$, $d$, the entry decisions of the CPs, and the consumer decisions with regards to connecting to the ISP and the CPs. The solution concept used to analyze the model is subgame perfect Nash equilibrium. We will use the backward induction algorithm to find subgame perfect Nash equilibria of the multi-stage game.

Suppose we are in the subgame in which the mass of consumers who chose to connect to the ISP is $\hat{x}$, the number of CPs who chose to enter is $n$ and they place themselves symmetrically on the \emph{Salop} circle whose perimeter is $\hat{x}$. So $CP_i$ faces a demand $D_i = \frac{\hat{x}}{n}$ with consequent profits
\begin{equation}
\pi_i = r \frac{\hat{x}}{n} - f - a.
\end{equation}
Under free entry, the optimal number of CPs who enter is determined by the zero-profit condition yielding
\begin{equation}\label{CPDemand}
n = \frac{r \hat{x}}{a+f}.
\end{equation}
Ceteris paribus, an increase in ISP's access price $a$ to CPs just adds to their fixed costs and reduces the number of CPs which enter the market. Similarly, an increase in total demand $\hat{x}$ encourages entry.

The marginal consumer at $\hat{x}$ travels a distance $\frac{\hat{x}}{2n}$ and is indifferent between connecting to the ISP and not connecting to it, that is it has a utility of $0$. Setting the utility $u$ in~(\ref{eqn:consumerUtil}) to zero, we can write
\begin{equation*}
v - \hat{x} - \frac{t \hat{x}}{2n} - d - \frac{1}{\mu- \lambda \hat{x}} = 0.
\end{equation*}

Using $n$ from Equation~(\ref{CPDemand}) in the above equation, we have
\begin{equation}   \label{ConDemand}
v - \hat{x} - \frac{t (a+f)}{2r} - d - \frac{1}{\mu- \lambda \hat{x}} = 0.
\end{equation}
Let $\hat{x}(d, a)$ solve Equation (\ref{ConDemand}). Then
\begin{align}
\hat{x}(d, a)  &= \frac{1}{2\lambda}\left(\mu + \lambda \big(v-d - \frac{t (a+f)}{2r} \big)\right.\nonumber\\
&\left.- \sqrt{[\mu - \lambda \big(v-d - \frac{t (a+f)}{2r} \big)]^2 + 4 \lambda}\right). \label{Demand1}
\end{align}
Further, using~(\ref{CPDemand}), we can calculate the number of CPs
\begin{align}
n(d, a) &= \frac{r \hat{x}(d, a)}{a+f}.  \label{Demand2}
\end{align}
So the total demand on either side is decreasing in both the access prices.

The ISP chooses the access prices $d$ and $a$ to maximize its profits subject to the constraints that the resulting consumer demand is nonnegative and the demand on the CP side is at least 1. Formally it faces
\begin{align*}
[OPT1] \quad \quad  \max_{d, a} \quad \pi(d, a) &= d \hat{x}(d, a) + a n(d,a) \\
\text{subject to } \quad 0 \le \hat{x}(d, a) &\le 1,   \\
n(d,a) &\ge 1.
\end{align*}

Define $e(\hat{x}(d, a), d) := - \frac{\ud \hat{x}(d, a)/ \hat{x}(d, a)}{\ud d/ d}$ to be the elasticity of consumer demand $\hat{x}(d, a)$ with respect to $d$. Similarly, define $e(\hat{x}(d, a), a)$ as the elasticity with respect to $a$. Let $e(n(d, a), d)$ and $e(n(d, a), a)$ be the corresponding elasticities of the number of CPs. Then an unconstrained optimum is a solution to the following pair of equations in $(d, a)$.

\begin{align}
e(\hat{x}(d, a), d) + \frac{an}{d \hat{x}} e(n(d, a), d) &= 1,   \label{FOC_d} \\
e(n(d, a), a) + \frac{d \hat{x}}{an} e(\hat{x}(d, a), a) &= 1.   \label{FOC_a}
\end{align}

Equations (\ref{FOC_d}) and (\ref{FOC_a}) can be interpreted as follows. If the monopolist ISP wants to increase $d$, he faces the traditional margin vs volume tradeoff. That is, he loses some consumers but makes more on the existing consumers. However, since the ISP is operating in a two-sided market, an additional effect arises via network externalities. A decreases in consumer market size makes entry unattractive for CPs. So the provider market size decreases which reduces the ISP's revenues from that market. Similarly, if the ISP wants to increase $a$, in addition to facing the traditional margin vs volume tradeoff \SKK{The trade-off is unclear and uses terminology that is unfamiliar to me. Also, while the interpretation of the equations is clear, it is not apparent how the interpretation follows from the equations.} in the provider market, it has to consider that increasing $a$ also decreases the consumer market size which has a second order effect on further reducing volumes in the provider market. The optimal access prices $(d, a)$ balance these effects simultaneously.

\subsection{Net Neutrality Regime}
Net neutrality regulation constrains the access price $a$ that the ISP can charge CPs to be zero. In the neutrality regime, the ISP, the CPs and the consumers interact through the same multi-stage game as in free market regime but with the constraint $a = 0$.

The ISP faces the following optimization problem for choosing the access price to consumers
\begin{align*}
[OPT2] \quad \quad \max_{d} \quad \pi(d, a=0) &= d \hat{x}(d, a=0)   \\
\text{subject to } \quad 0 \le \hat{x}(d, a=0) &\le 1,   && (C_1)\\
n(d, a=0) &\ge 1,     && (C_2)
\end{align*}
where $\hat{x}(d, a=0)$ and $n(d, a=0)$ are found by setting $a = 0$ in Equations (\ref{Demand1}) and (\ref{Demand2}).

The unconstrained optimum of $[OPT2]$ is given by the value of $d$ at which the standard margin vs volume tradeoff \SKK{Unclear terminology; related to previous comment} in the consumer market is resolved optimally at
\begin{equation}   \label{NeutralOpt}
e(\hat{x}(d, 0), d) = 1
\end{equation}
and is the solution if it lies in the constraint set. The corresponding demands $\hat{x}(d, a=0)$ and $n(d, a=0)$ on both sides can be computed by setting $a=0$ and $d$ at its optimal value in Equations (\ref{Demand1}) and (\ref{Demand2}).

\subsection{Welfare}
\label{sec:Welfare}
The model we are studying is a transferable utility model. This implies that the cardinal welfare $W$ at any outcome is defined to be the sum total of consumers' surplus ($CS$), profits ($\pi_i$) of all CPs who enter and profit $\pi$ of the ISP. We have these expressions
\begin{align}
\pi &= d \hat{x} + na,    \nonumber  \\
\sum_{i=1}^n \pi_i &= r \hat{x} - nf - na,  \nonumber \\
CS &= \int_{0}^{\hat{x}} \big(v-d - \frac{1}{\mu- \lambda \hat{x}} -x \big) \ud x - \nonumber\\
&\quad n \int_{0}^{\hat{x}/n} ty \ud y  = \big(v-d-\frac{1}{\mu- \lambda \hat{x}} \big) \hat{x} - \frac{\hat{x}^2}{2} - \frac{t \hat{x}^2}{n}, \nonumber \\
W(\hat{x}, n) &= \pi + \sum_{i=1}^n \pi_i + CS = (v+r) \hat{x} -  \frac{\hat{x}^2}{2} - \frac{t \hat{x}^2}{n}\nonumber\\ &\quad - nf - \frac{\hat{x}}{\mu- \lambda \hat{x}}.  \label{Welfare}
\end{align}
The access prices cancel out because these are mere transfers from one party to another. The overall welfare in an equilibrium depends on the consumer market size and the number of CPs that populate the provider market. If the consumer market size $\hat{x}$ increases, consumer surplus and advertising surplus increase on account of more consumers being served but welfare decreases on account of two factors -- first, an increase in consumer market size increases the heterogeneity in consumer preferences which causes the product diversity along the ISP and the CP dimension to decrease relatively; second, an increase in consumer market size also increases the congestion cost for everyone. Similarly, an increase in the number of CPs who enter increases consumer surplus because there is more product diversity on the CP dimension for consumers but imposes greater fixed costs on the society.

As is standard in economic analysis, we compute the welfare optimum as a solution to the following problem.
\begin{align*}
[OPT3] \quad \quad \max_{\hat{x}, n} \quad & W(\hat{x}, n)   \\
\text{subject to } \quad 0 \le \hat{x} &\le 1,   \\
n &\ge 1.
\end{align*}
The solution to $[OPT3]$ gives the size of consumer and provider market that is best from society's point of view when we have taken everyone's interest into account.
\section{Basic Model}
\label{sec:model}
There is a monopoly ISP and many potential CPs for whom this ISP is not the originating ISP. Introduce the interval $[0,1]$ as the \emph{Hotelling} line of the standard horizontal product differentiation model for the ISPs. Let the ISP be located at the left extremity, location $0$, of this line. There is a continuum of consumers uniformly distributed on the Hotelling line. Let the product differentiation parameter for the ISP be normalized to $1$. Non-price factors therefore impose a cost of $x$ on a consumer at location $x$ on the Hotelling line.

Let $t \in (0,1)$ be the product differentiation parameter for the CPs. Note that we restrict $t<1$. This is to capture the fact that it is easier to substitute between CPs than it is to substitute between ISPs. Let $y$ be the distance travelled by a consumer to connect to a CP. This imposes a cost of $ty$ on the consumer. Consumers cannot connect to the CPs without first connecting to the ISP. The CPs space themselves equidistantly on a Salop circle whose perimeter, as we detail soon, is determined by the consumers that choose to connect to the ISP.

All CPs provide free content to the consumers; their only source of revenues is from advertising. Let $r$ be the advertising revenue per consumer, same for all CPs. Let $d$ be the access price charged by the ISP to the consumer. Let $f$ be the fixed cost of entering the provider market for any CP.

Consumers pay a congestion cost, which is the average time the ISP takes to fulfill a consumer's request for content. To capture this, we model the ISP as a M/M/1~\cite{gross1998fundamentals} service facility. Specifically, following~\cite{pil2010net}, we assume that the aggregate demand for content from consumers created at the ISP is well modeled by a Poisson process of rate $\lambda$ per unit length of the single dimension containing consumers connected to the ISP. Further, the time taken by the ISP to service a consumer's request for content is exponentially distributed with a mean service time of $1/\mu$. We assume $\lambda < \mu$, which is required for the average time to fulfill a request to stay bounded. Finally, consumer requests are serviced by the ISP in a first-come-first-served manner. This means that a request that arrives at the ISP must wait for requests that arrived earlier to finish obtaining service from the ISP.

Consumer monetary valuation of the total requested content is $v$. In the non-neutral regime, the ISP can charge the CPs for servicing requests for their content at an access price of $a$. In contrast, the neutrality regime prohibits the ISP from charging the CPs a price when it is not the originating ISP for them, which we assume is the case here. The two regimes therefore induce different market games through which the entities interact.

\subsection{Non-neutral Regime}
\begin{figure}
	\centering
	\includegraphics[width=0.9\columnwidth]{Figures/Fig_1.pdf}
	\caption{\small{Left Panel shows the \emph{Hotelling} line for the ISP dimension. Right panel shows the folded \emph{Hotelling} line segment of length $\hat{x}$ as the \emph{Salop} circle with the same circumference. It also shows the symmetric alignment of CPs when $n=4$ for when the consumer at $\hat{x}$ travels a distance $\hat{x}/(2n) = \hat{x}/8$ to reach a CP.}}
	\label{fig:ModFig}
\end{figure}
In the non-neutral regime, the ISP, the CPs and the consumers interact through the following \emph{multi-stage game}:
\begin{enumerate
\item ISP chooses the access prices: $a$ for CPs and $d$ for consumers.
\item Consumers then choose whether to connect to the ISP.
\item CPs who are considering entering the market then simultaneously decide whether to enter or not. Suppose $n$ is the number of CPs who chose to enter. The CPs locate equidistantly\footnote{We follow \cite{salop1979monopolistic} in making the assumption of symmetric alignment of CPs on the Salop circle to focus on the entry decisions of CPs.} on a \emph{Salop} circle whose perimeter is the length $\hat{x}$ of the initial segment of the \emph{Hotelling} line on which consumers connect to the ISP.
\item Consumers who have connected to the ISP then connect to a CP and consume content.
\end{enumerate}

We specify utility of a consumer who travels a distance $x$ on the \emph{Hotelling} line to connect to the ISP and distance $y$ on the \emph{Salop} circle to connect to the CP as
\begin{equation}\label{eqn:consumerUtil}
u = v - x  - ty - d - \frac{1}{\mu- \lambda \hat{x}},
\end{equation}
where $\lambda \hat{x}$ is the aggregate rate of request for content from the length $\hat x$ of consumers. The last term in the expression for consumer utility $u$ is the congestion cost faced by a consumer. It is the resulting expected time the ISP takes to fulfill a consumer's request (called waiting time in~\cite{gross1998fundamentals}), given the demand $\lambda \hat{x}$.

We specify ISP's profit as $\pi(d, a) = d \hat{x} + na$. We specify profits of $\text{CP}_i$ as $\pi_i = r D_i - a - f$, where $D_i$ is the demand that $\text{CP}_i$ faces.

To determine this demand $D_i$ we must decide how the CPs are aligned on the Salop circle relative to the marginal consumer at $\hat{x}$ who connects to the ISP. This is important for it determines the distance traveled by this consumer on the Salop circle and therefore affects the location $\hat{x}$ on the Hotelling line. Two extreme cases are possible. In one, illustrated in Figure~\ref{fig:ModFig}, the marginal consumer travels a distance $\frac{\hat{x}}{2n}$; in the other, the consumer travels no distance, that is a CP locates itself at $\hat{x}$. %
One could introduce an extra parameter $\beta \in [0,1]$, and specify the distance traveled by the marginal consumer as $\frac{\beta \hat{x}}{2n}$. The extreme cases correspond to the values $\beta = 1$ and $\beta = 0$. We will, however, choose $\beta = 1$ in our analysis for two reasons. First, any $\beta \in (0, 1]$ retains the structural feature that demand on either side (whether $\hat{x}$ or $n$) responds to prices ($d$ and $a$) on both sides. Therefore $\beta = 1$ carries the flavor of analysis of any $\beta \in (0,1)$. Second, as will be clearer from analysis that follows, $\beta = 0$, kills any dependence of $\hat{x}$ on $a$ and in that sense is a special case.

To summarize, the parameters of the model are the consumer valuation $v$, the advertising revenue $r$, the normalized product differentiation parameter $t$, fixed CP costs $f$, the rate of content request $\lambda$, and the service rate $\mu$. The choice variables are the prices $a$, $d$, the entry decisions of the CPs, and the consumer decisions with regards to connecting to the ISP and the CPs. The solution concept used to analyze the model is subgame perfect Nash equilibrium. We will use the backward induction algorithm to find subgame perfect Nash equilibria of the multi-stage game.

Suppose we are in the subgame in which the mass of consumers who chose to connect to the ISP is $\hat{x}$, the number of CPs who chose to enter is $n$ and they place themselves symmetrically on the \emph{Salop} circle whose perimeter is $\hat{x}$. So $CP_i$ faces a demand $D_i = \frac{\hat{x}}{n}$ with consequent profits
\begin{equation}
\pi_i = r \frac{\hat{x}}{n} - f - a.
\end{equation}
Under free entry, the optimal number of CPs who enter is determined by the zero-profit condition yielding
\begin{equation}\label{CPDemand}
n = \frac{r \hat{x}}{a+f}.
\end{equation}
Ceteris paribus, an increase in ISP's access price $a$ to CPs just adds to their fixed costs and reduces the number of CPs that enter the market. Similarly, an increase in total demand $\hat{x}$ encourages entry.

The marginal consumer at $\hat{x}$ travels a distance $\frac{\hat{x}}{2n}$ and is indifferent between connecting to the ISP and not connecting to it, that is it has a utility of $0$. Setting the utility $u$ in~(\ref{eqn:consumerUtil}) to zero, we can write
\begin{equation*}
v - \hat{x} - \frac{t \hat{x}}{2n} - d - \frac{1}{\mu- \lambda \hat{x}} = 0.
\end{equation*}

Using $n$ from Equation~(\ref{CPDemand}) in the above equation, we have
\begin{equation}   \label{ConDemand}
v - \hat{x} - \frac{t (a+f)}{2r} - d - \frac{1}{\mu- \lambda \hat{x}} = 0.
\end{equation}
Let $\hat{x}(d, a)$ solve Equation (\ref{ConDemand}). Then
\begin{align}
\hat{x}(d, a)  &= \frac{1}{2\lambda}\left(\mu + \lambda \big(v-d - \frac{t (a+f)}{2r} \big)\right.\nonumber\\
&\left.- \sqrt{[\mu - \lambda \big(v-d - \frac{t (a+f)}{2r} \big)]^2 + 4 \lambda}\right). \label{Demand1}
\end{align}
Further, using~(\ref{CPDemand}), we can calculate the number of CPs
\begin{align}
n(d, a) &= \frac{r \hat{x}(d, a)}{a+f}.  \label{Demand2}
\end{align}
So the total demand on either side is decreasing in both the access prices.

The ISP chooses the access prices $d$ and $a$ to maximize its profits subject to the constraints that the resulting consumer demand is nonnegative and the demand on the CP side is at least 1. Formally it faces
\begin{align*}
[OPT1] \quad \quad  \max_{d, a} \quad \pi(d, a) &= d \hat{x}(d, a) + a n(d,a)\\
\text{subject to } \quad 0 \le \hat{x}(d, a) &\le 1,\\
n(d,a) &\ge 1.
\end{align*}
Define $e(\hat{x}(d, a), d) := - \frac{\ud \hat{x}(d, a)/ \hat{x}(d, a)}{\ud d/ d}$ to be the elasticity of consumer demand $\hat{x}(d, a)$ with respect to $d$. Similarly, define $e(\hat{x}(d, a), a)$ as the elasticity with respect to $a$. Let $e(n(d, a), d)$ and $e(n(d, a), a)$ be the corresponding elasticities of the number of CPs. Then an unconstrained optimum is a solution to the following pair of equations in $(d, a)$.
\begin{align}
e(\hat{x}(d, a), d) + \frac{an}{d \hat{x}} e(n(d, a), d) &= 1,\label{FOC_d}\\
e(n(d, a), a) + \frac{d \hat{x}}{an} e(\hat{x}(d, a), a) &= 1.\label{FOC_a}
\end{align}
Equations (\ref{FOC_d}) and (\ref{FOC_a}) can be interpreted as follows. If the monopolist ISP wants to increase $d$, he faces the traditional margin vs volume tradeoff\footnote{A profit maximizing monopolist facing a downward sloping demand curve faces a tradeoff when choosing an optimal price. While a higher price implies he makes more on the units that he sells (this is the margin effect), it also implies a lower sales volume (this is the volume effect) because a higher price leads to a lower demand. }. That is, he loses some consumers but makes more on the existing consumers. However, since the ISP is operating in a two-sided market, an additional effect arises via network externalities. A decrease in consumer market size makes entry unattractive for CPs. So the provider market size decreases which reduces the ISP's revenues from that market\footnote{That the aforementioned is indeed an interpretation of Equation (\ref{FOC_d}) is established by looking at the more fundamental first order condition with respect to $d$, which is
\begin{align}
&\frac{\ud \pi}{\ud d}= \underbrace{\hat{x}}_{\text{margin effect}} + \underbrace{d \frac{\ud \hat{x}}{\ud d}}_{\text{volume effect}} + \underbrace{a \frac{\ud n}{\ud d}}_{\text{network externalities effect}} = 0.\nonumber
\end{align
Equation (\ref{FOC_d}) is equivalent to the above equation under the definition of elasticities.}. Similarly, if the ISP wants to increase $a$, in addition to facing the traditional margin vs volume tradeoff in the provider market, it has to consider that increasing $a$ also decreases the consumer market size which has a second order effect on further reducing volumes in the provider market. The optimal access prices $(d, a)$ balance these effects simultaneously.

\subsection{Net Neutrality Regime}
Net neutrality regulation constrains the access price $a$ that the ISP can charge CPs to be zero. In the neutrality regime, the ISP, the CPs and the consumers interact through the same multi-stage game as in the non-neutral regime but with the constraint $a = 0$.

The ISP faces the following optimization problem for choosing the access price to consumers
\begin{align*}
[OPT2] \quad \quad \max_{d} \quad \pi(d, a=0) &= d \hat{x}(d, a=0)   \\
\text{subject to } \quad 0 \le \hat{x}(d, a=0) &\le 1,   && (C_1)\\
n(d, a=0) &\ge 1,     && (C_2)
\end{align*}
where $\hat{x}(d, a=0)$ and $n(d, a=0)$ are found by setting $a = 0$ in Equations (\ref{Demand1}) and (\ref{Demand2}).

The unconstrained optimum of $[OPT2]$ is given by the value of $d$ at which the standard margin vs volume tradeoff in the consumer market is resolved optimally at
\begin{equation}   \label{NeutralOpt}
e(\hat{x}(d, 0), d) = 1
\end{equation}
and is the solution if it lies in the constraint set. The corresponding demands $\hat{x}(d, a=0)$ and $n(d, a=0)$ on both sides can be computed by setting $a=0$ and $d$ at its optimal value in Equations (\ref{Demand1}) and (\ref{Demand2}).

\subsection{Welfare}
\label{sec:Welfare}
The model we are studying is a transferable utility model. This implies that the cardinal welfare $W$ at any outcome is defined to be the sum total of consumers' surplus ($CS$), profits ($\pi_i$) of all CPs who enter and profit $\pi$ of the ISP. We have these expressions
\begin{align}
\pi &= d \hat{x} + na,    \nonumber  \\
\sum_{i=1}^n \pi_i &= r \hat{x} - nf - na,  \nonumber \\
CS &= \int_{0}^{\hat{x}} \big(v-d - \frac{1}{\mu- \lambda \hat{x}} -x \big) \ud x - \nonumber\\
&\quad n \int_{0}^{\hat{x}/n} ty \ud y  = \big(v-d-\frac{1}{\mu- \lambda \hat{x}} \big) \hat{x} - \frac{\hat{x}^2}{2} - \frac{t \hat{x}^2}{n}, \nonumber \\
W(\hat{x}, n) &= \pi + \sum_{i=1}^n \pi_i + CS = (v+r) \hat{x} -  \frac{\hat{x}^2}{2} - \frac{t \hat{x}^2}{n}\nonumber\\ &\quad - nf - \frac{\hat{x}}{\mu- \lambda \hat{x}}.  \label{Welfare}
\end{align}
The access prices cancel out because these are mere transfers from one party to another. The overall welfare in an equilibrium depends on the consumer market size and the number of CPs that populate the provider market. If the consumer market size $\hat{x}$ increases, consumer surplus and advertising surplus increase on account of more consumers being served but welfare decreases on account of two factors -- first, an increase in consumer market size increases the heterogeneity in consumer preferences which causes the product diversity along the ISP and the CP dimension to decrease relatively; second, an increase in consumer market size also increases the congestion cost for everyone. Similarly, an increase in the number of CPs who enter increases consumer surplus because there is more product diversity on the CP dimension for consumers but imposes greater fixed costs on the society.

As is standard in economic analysis, we compute the welfare optimum as a solution to the following problem.
\begin{align*}
[OPT3] \quad \quad \max_{\hat{x}, n} \quad & W(\hat{x}, n)   \\
\text{subject to } \quad 0 \le \hat{x} &\le 1,   \\
n &\ge 1.
\end{align*}
The solution to $[OPT3]$ gives the size of consumer and provider market that is best from society's point of view when we have taken everyone's interest into account.

\section{Related Work}
\label{sec:related}
The contours of the debate on net neutrality are by now well known and well articulated in works like~\cite{kramer2013net} and~\cite{greenstein2016net}. As we state in the introduction, our modeling objective is to compare in welfare terms the equilibriums in neutral and non-neutral regimes. Non-neutrality introduces many feasible strategies for the ISP, including two-sided pricing and the possibility of prioritization of network lanes. The literature has progressed on both these dimensions. In the direction of work on lane prioritization, in~\cite{ma2017paid} it is found that paid prioritization can be welfare-superior to net neutrality. In~\cite{zou2017optimal} the authors investigate optimal service differentiation for an ISP subject to capacity constraints in an optimal control framework. In~\cite{gharakheili2016economic} non-neutrality is studied in an economic model that allows lane pricing by ISP and revenue generation for CP as well as QoE improvement for the consumers. The economic benefits for ISP and CP are analyzed using traffic traces from real networks. However, lane prioritization is not the focus of our work. In doing this, our work is closer in spirit to~\cite{economides2012network} where too authors abstract from lane prioritization issues.

In the context of a two-sided pricing model, wherein ISP charges both consumers and content providers,~\cite{altman2011network} studies non-neutrality in an economic model including an ISP, a CP, consumers and advertisers and shows that side-payments are beneficial for both the ISP and the CP. In~\cite{andrews2013economic} the interaction between an ISP and a CP is modeled as a Stackelberg game with the objective of maximizing the total profit of both players. In~\cite{coucheney2013impact} the authors study the impact of competition between ISPs on network neutrality. Authors in~\cite{kamiyama2014feasibility} analyze a model with competing ISPs, multiple non-competing CPs and a fixed number of consumers and study the feasibility of imposing side-payments on CPs and its effect on the revenue of both ISP and CP. Other works such as~\cite{ma2013public} and~\cite{mitra2017} study the impact of net-neutrality on consumer utility and investments in the network focusing on incremental bandwidth allocation and caching, respectively.

While prior works consider competing ISPs and multiple non-competing content providers, to the best of our knowledge, our work is the first attempt to study non-neutrality with a focus on competition issues in the provider market for a particular content type. We are interested in how the potential of two-sided pricing introduced by non-neutrality affects the competitive landscape of the provider market, as well as the consumer market base, with consequences for welfare.

In terms of modeling choices, our work is closely related to~\cite{economides2012network} and~\cite{pil2010net}. However, while we study competition in the content provider market, authors in~\cite{economides2012network} model content providers as many independent monopolies. Also, in our work consumers face a congestion cost that is absent in~\cite{economides2012network}.

\section{Modeling Consumer Choice and Competition}
\label{sec:spatialModels}
We need to model consumer preference for competing content providers and for the ISP connection. We will think of both content and the ISP connection as economic goods of which the consumers are the buyers and the content providers and the ISP are respectively the sellers.

Consumer preference is often guided by a mix of price and a myriad of non-price factors. To exemplify, while both Flipkart and Amazon are content providers in the e-commerce space in India and they have similar transaction costs, some Indian consumers may prefer one over the other. Such non-price factors are modeled in the industrial organization~\cite{tirole1988theory} literature using \emph{spatial models of product differentiation}. Spatial models use an abstraction of a single dimension -- a \emph{Hotelling} line or the perimeter of a \emph{Salop} circle -- to capture all non-price factors. Buyers and sellers are identified by locations on this dimension. A buyer's distance from a seller on the single dimension captures the non-price aspects of the buyer's preference for the seller. On this dimension, for our example, a consumer that prefers Flipkart will be closer to Flipkart than to Amazon.

In addition to the above, spatial models have a product differentiation parameter that captures the ease with which a buyer can substitute between offerings of competing sellers. The smaller this parameter, the easier it is for the buyer to substitute one seller for another. Specifically, the cost incurred by a buyer to change its position (relative to the sellers) on the single dimension is smaller.

To summarize, we use a richer model of price competition with differentiated products that is closer to reality than the simplistic model of price competition with no differentiation. We remark that in~\cite{coucheney2013impact} an alternative price based discrete choice model is used to address many of the same modeling issues. However, we find the spatial models attractive for their relative simplicity.


\section*{Supplementary Note}
\if0
\subsection{Suppression of the superstructure in the Bi-O layer}
We performed the X-ray diffraction (XRD) measurement on \ce{Bi_{1.7}Pb_{0.5}Sr_{1.9}CaCu2O_{8+$\delta$}} (Pb-Bi2212) which is the very sample measured in the ARPES experiment in order to check whether the superstructural modulation of \ce{Bi_{2}Sr_{2}CaCu2O_{8+$\delta$}} (Bi2212) was suppressed by substituting Pb for Bi \cite{MUSOLINO20031}. We used a High Luminance In-plane Type X-ray Diffractometer SmartLab 9kW (Rigaku) at room temperature.
The absence of the superstructure modulation was confirmed by the absence of the satellite peaks which indicate the superstructure modulation around the main (020)${}_\text{o}$ and (040)${}_\text{o}$ peaks. Here, the orthorhombic notation for unit cell size of $\sqrt{2}a \times \sqrt{2}a$ is employed.

\begin{figure}[h]
\includegraphics[width = 11cm, bb = 0 0 820 537]{20181114_S1SA.png
\caption{\label{fig:S1} In-plane XRD for Pb-Bi2212. Black arrows indicate the (020)${}_\text{o}$ and (040)${}_\text{o}$ Bragg peaks in Pb-Bi2212. There are no satellite peaks which reflect the superstructure modulation around the main peaks. Extrinsic four peaks marked by * are supposed to come from other elements of the mechanical strain device.}
\end{figure}
\fi

\subsection{Strain estimated from the lattice parameters }
After the ARPES measurements, we carried out XRD measurements in order to estimate the magnitude of the tensile strain applied to the sample by the screw shown in Fig. \ref{fig:S2}(c). We estimated the strain from the anisotropy of the $a$ and $b$ lattice parameters because it is difficult to measure the strain directly because of the plasticity of the device. As shown in Fig. \ref{fig:S2}(a), we focused on the (200) and (020) Bragg peaks of Pb-Bi2212 and each peak was fitted using a single Gaussian. Figure \ref{fig:S2}(b) shows the $a$ and $b$ lattice parameters estimated as $a = 3.814$  \AA \, and $b = 3.812$ \AA. The anisotropy of the crystal structure is $(a-b)/(a+b) < 0.03$ \%. The anisotropy of the electronic structure  induced by the distortion $\delta_\text{crystal} = 7 \times (a-b)/(a+b) < 0.2$ \% expected from Eq. (S4.3) is much smaller than the asymmetry parameter $\delta$ of order 1\% which we observed in the ARPES measurements shown in Fig. 3. Therefore, the enhancement of the nematicity in the pseudogap state cannot be explained due merely to the anisotropy of the crystal structure.
A picture and a schematic figures of the strain device are shown in Fig \ref{fig:S2}(c).

\begin{figure}[h]
\includegraphics[width = 15cm, bb = 0 0 722 677]{20181114_S1SA.png
\caption{\label{fig:S2} Strength of the tensile strain estimated from lattice parameter. (a) X-ray diffraction pattern around the (200) and (020) peaks of Pb-Bi2212. Data were fitted using a single Gaussian. (b) Lattice parameters of $a$ and $b$ estimated from the Bragg peaks using Bragg's law.  The difference between $a$ and $b$ was less than 0.03 \%. (c) Sample holder designed to apply strain to the sample by tightening the stainless screw. The top view is shown on the right.}
\end{figure}

\clearpage
\subsection{Matrix-element effect}

In vacuum, the vector potential $\vec{A}(\vec{r},t)$ satisfies the wave equation,
\begin{equation}
\nabla ^2 \vec{A}(\vec{r},t) = \frac{1}{c^2} \frac{\partial ^2 \vec{A}(\vec{r},t)}{\partial t^2}. \label{eq:2.5.1}
\end{equation}
One of the solutions of this equation is
\begin{equation}
\vec{A}(\vec{r},t) = \vec{\varepsilon}e^{i(\vec{k}\cdot\vec{r}-\omega t)}, \label{eq:2.5.2}
\end{equation}
where $\vec{\varepsilon}$ denotes the polarization (the oscillating direction of the electric field), and $\vec{k}$ and $\omega$ are the wave vector and the frequency of light, respectively ($\omega = c k$). The probability of photoemission is given by
\begin{equation}
w_{\text{i} \rightarrow \text{f}} \propto \left|\braket{f|H'|i}\right|^2,\label{eq:2.5.3}
\end{equation}
where
\begin{equation}
H' = \frac{e}{m}\vec{A}(\vec{r})\cdot \vec{p} = \frac{e}{m}\vec{\varepsilon}e^{i\vec{k}\cdot\vec{r}}\cdot \vec{p},\label{eq:2.5.4}
\end{equation}
according to Fermi's golden rule. In ARPES measurement, we often use vacuum-ultra violet, whose photon energy is about 100 eV i.e., the wave length is about 100 \AA { } which is much longer than the distance between neighboring atoms in the material, so that one obtains
\begin{equation}
e^{i\vec{k}\cdot\vec{r}} = 1 + i\vec{k}\cdot\vec{r} + \cdots  \simeq 1  .\label{eq:2.5.5}
\end{equation}
Then the matrix element in eq (\ref{eq:2.5.3}) can be expressed as
\begin{eqnarray}
\braket{f|H'|i} &\simeq& \braket{f|\frac{e}{m}\vec{\varepsilon}\cdot \vec{p}|i}\\ \label{eq:2.5.6}
&=& \frac{e}{i\hbar} \braket{f|\vec{\varepsilon}\cdot [\vec{x},H_0]|i}\\
&=& \frac{e}{i\hbar} \braket{f|\vec{\varepsilon}\cdot \vec{x}H_0^\dagger - H_0\vec{\varepsilon}\cdot \vec{x}|i}\\
&=& \frac{e}{i\hbar}(E_\text{f}-E_\text{i}) \braket{f|\vec{\varepsilon}\cdot \vec{x}|i},\label{eq:2.5.9}
\end{eqnarray}
using the Hermite property of $H_0$ and the commutation relationship
\begin{equation}
[\vec{x},H_0]=\frac{i\hbar \vec{p}}{m}.\label{eq:2.5.10}
\end{equation}
In summary, the transition probability can be approximated by
\begin{equation}
w_{\text{i} \rightarrow \text{f}} \propto \left|\braket{f|\vec{\varepsilon}\cdot \vec{x}|i}\right|^2,\label{eq:2.5.11}
\end{equation}
and the matrix element can be described by that of the dipole transition.\par
The $k$-dependence of the intensity in $k$-space obtained by ARPES can be estimated by calculating the matrix element (eq. (\ref{eq:2.5.11})) in the case that $\ket{i} = \ket{d_{x^2-y^2}}$, $\vec{E} \propto (+1,-1,0)$, which represents our experimental condition, and $\ket{f} = e^{i\vec{k}\cdot\vec{r}}$ was assumed for simplicity. Thus, we obtain
\begin{eqnarray}
|\braket{f|\vec{E}\cdot \vec{r}|i}|^2 &=& \left|\braket{e^{i\vec{k}\cdot\vec{r}}|\vec{E}\cdot \vec{r}|d_{x^2-y^2}}\right|^2\\
&\propto& \left|\int d^2\vec{r} (x^2-y^2)(x-y)\cos(k_xx+k_yy)\right|^2\\
&& +\left|\int d^2\vec{r} (x^2-y^2)(x-y)\sin(k_xx+k_yy)\right|^2.
\end{eqnarray}
Figure \ref{fig:2.5.1} is the result of the $k$-dependence of the transition probability in the first Brillouin zone.
\begin{figure}[h]
\begin{center}
\includegraphics[width=9cm,bb = 0 0 634 524]{ME4}
\caption{Intensity plot in the first Brillouin zone for $\ket{i} = \ket{d_{x^2-y^2}}$, $\vec{E} \propto (+1,-1,0)$, and $\ket{f} = e^{i\vec{k}\cdot\vec{r}}$ is assumed.}
\label{fig:2.5.1}
\end{center}
\end{figure}

\clearpage

\subsection{Energy dependence of the FWHM of the momentum distribution curves in the normal, pseudogap, and superconducting states}

Figures S3 - S5 show the MDC for the off-nodal cut and  the FWHM of the two MDC peaks as functions of binding energy for the superconducting, pseudogap, and normal states, respectively. The figures show that the FWHM at $E-E_\text{F}$ = -30 meV is small compared to other binding energies, because around $E-E_\text{F}$ = -10 meV the superconducting gap opens, and in the higher binding energy region, where the quasi-particle scattering rate increases.

\begin{figure}[h]
\includegraphics[width = 15cm, bb = 0 0 610 579]{20181005_S3SA.png
\caption{\label{fig:S3} Energy dependence of the FWHM of MDC from ARPES spectra taken in the superconducting state (5 K). (a) FWHM of the MDCs shown in (b) at off-nodal k points where the cut shown in the inset crosses the Fermi surface. The solid curves are fitted ones using a quadratic function. The MDC is broad around -10 meV due to the superconducting gap and around -50 meV because of the higher scattering rate of quasi-particles. (b) MDCs at $E-E_\text{F} =$ -30 meV fitted to  two Lorentzians. (c) Band dispersion along the cut shown in the inset.}
\end{figure}

\begin{figure}[h]
\includegraphics[width = 15cm, bb = 0 0 608 574]{20181005_S4SA.png
\caption{\label{fig:S4} The same as Fig. S4 except for  taken in the pseudogap state (100 K). }
\end{figure}

\begin{figure}[h]
\includegraphics[width = 15cm, bb = 0 0 501 459]{20181005_S5SA.png
\caption{\label{fig:S5} The same as Fig. S4 except for taken in the normal state (200 K).}
\end{figure}

\clearpage
\subsection{Effect of strain on hopping parameters}

One may suspect that the anisotropy arises only from the effect of the uniaxial strain and that the temperature dependence of the anisotropy is caused by thermal expansion. Such a possibility can be excluded by the following consideration:
Here, we consider a transfer integral between atomic orbitals of different atoms {\it i} and {\it j}. When the axis connecting both atoms is taken as the $z$-axis, the transfer integral between the i and j atoms $-t_{ij}(0,0,1) = \braket{i|h|j}$ is given by the Slater-Koster parameter $(l_il_j\mu)$, where $l_A$ is the azimuthal quantum number of atom $A$ $(A=a, b)$ and $\mu$ is the magnetic quantum number of atoms $i$ and $j$.  In the following, $l_A = 0, 1, 2,...$ and $\mu = 0, \pm1, \pm2,...$ are noted as  {\it s, p, d},... and $\sigma, \pi, \delta,...$, respectively.\par
As for the transfer integral between the $2p$ orbital of O atom and the $3d$ orbital of Cu atom, there is approximate relationship between the Slater-Koster parameter and the distance between two atoms $d$ \cite{Harrison}
\begin{equation}
(pd \mu ) \propto d^{-3.5}. \label{eq:3.2.3}
\end{equation}
Hopping between nearest-neighbor Cu atoms can be regarded as consisting of two hopping processes from Cu to O and from O to another Cu and, therefore, the nearest-neighbor hopping parameter $t$ in eq (1) as a function of $d$ is given by
\begin{equation}
t { \propto (pd \mu )^2} \propto d^{-7}. \label{eq:3.2.4}
\end{equation}
If we assume that $t$ depends only on $d$ for simplicity, one can obtain
\begin{equation}
\frac{\Delta t}{t} \simeq -7\frac{\Delta d}{d}. \label{eq:3.2.5}
\end{equation}
If the anisotropy of $t$ which we observed is only due to the uniaxial strain, the anisotropy parameter $\delta \left(= \frac{\Delta t}{t} \right)$ should change monotonically because the atomic distance changes monotonically with temperature by thermal expansion. However, the anisotropy parameter $\delta$ observed in our measurement was enhanced only in the pseudogap state. This means that the stress works to align the nematic domain properly and ensures that the anisotropy does not come only from the distortion of the sample.

\clearpage

\input{Bi2212_nematicity_supplemenal.bbl}

\end{document}

\section{\label{sec:level1}Introduction}

\section*{Introduction}
The pseudogap in cuprate superconductors is characterized by the suppression of the density of states around the Fermi level above the superconducting transition temperature $T_\text{c}$.
Broadly speaking, two possible scenarios for the pseudogap have been discussed, that is, a precursor to the superconducting state and distinct order which competes with superconductivity \cite{Keimer2015From-quantum-ma}.
The latter scenario has been put forward by several measurements which detect distinct orders with spontaneous symmetry breaking such as charge order and time-reversal symmetry breaking at or below the pseudogap temperature $T^*$ \cite{Ghiringhelli821,Chang2012Direct-obs,Comin1335,Le-Tacon2013Inelastic-,Comin390,PhysRevLett.96.197001,PhysRevB.89.094523,PhysRevLett.118.097003,0953-8984-26-50-505701,PhysRevLett.88.137005,Hinkov597,Lawler2010Intra-unit,Sato2017Thermodyna,Daou2010Broken-rot}.

Recently, nematicity that the electronic structure preserves the translational symmetry but breaks the rotational symmetry of the underlying crystal lattice has been found to exist in the pseudogap state of the cuprate superconductors \cite{PhysRevLett.88.137005,Hinkov597,Lawler2010Intra-unit,Sato2017Thermodyna,Daou2010Broken-rot,PhysRevB.92.224502,Wu2017Spontaneou}.
As a possible mechanism, it has been proposed that the nematic order arises from fluctuations of stripe order \cite{Kivelson1998Electronic, RevModPhys.75.1201} or from the instability of the Fermi surface (so-called Pomeranchuk instability) \cite{PhysRevB.73.214517, PhysRevB.64.195109, PhysRevB.72.024502,PhysRevLett.85.5162}.
From experimental perspectives, the nematicity in the cuprate superconductors was first pointed out by the transport measurements of lightly-doped \ce{La_{2-$x$}Sr_$x$CuO4 } and \ce{YBa2Cu3O_$y$} (YBCO) \cite{PhysRevLett.88.137005}.
Anisotropic signals in the spin excitation measured by neutron scattering was observed for untwinned underdoped YBCO \cite{Hinkov597}.
Inequivalence between the oxygen atoms in the $a$ and $b$ directions was detected by scanning tunneling spectroscopy measurements for underdoped \ce{Bi2Sr2CaCu2O_{8+$\delta$}} (Bi2212) \cite{Lawler2010Intra-unit}.
Nernst effect and magnetic torque measurements on underdoped and optimally doped YBCO showed a systematic temperature dependence of the nematicity \cite{Daou2010Broken-rot,Sato2017Thermodyna}.
The onset temperature of the in-plane anisotropic signals in the Nernst coefficient and that in the magnetic susceptibility coincide with $T^*$.
Thus the pseudogap state can be considered as a distinct thermodynamic phase characterized by nematicity.
Although the orthorhombic distortion of the \ce{CuO2} plane caused by the Cu-O chains along the $b$-axis of the untwinned YBCO samples already breaks the four-fold rotational symmetry even above $T^*$, the orthorhombicity is considered to help to form a microscopic nematic domain in one particular direction and enables us to detect nematicity in macroscopic measurements.
That is, this weak orthorhombicity of the \ce{CuO2} plane for nematicity induced by the Cu-O chain has the same role as a weak external magnetic field for ferromagnets \cite{Vojta2009Lattice-sy}.

Motivated by those previous studies on nematicity in cuprate superconductors, we have performed angle-resolved photoemission spectroscopy (ARPES) measurements on \ce{Bi_{1.7}Pb_{0.5}Sr_{1.9}CaCu2O_{8+$\delta$}} (Pb-Bi2212) by applying a uniaxial strain in the Cu-O bond direction in order to detect nematicity below $T^*$.
While most of the studies on nematicity in the cuprate superconductors have been done on YBCO, Pb-Bi2212, whose \ce{CuO2} plane has the tetragonal ($C_4$) symmetry, is a more appropriate material to study nematicity than YBCO.
Owing to the presence of a natural cleavage plane between the BiO layers and rich information accumulated from the previous ARPES studies, Bi2212 is an ideal material for ARPES experiment to investigate novel phenomena in cuprates \cite{RevModPhys.75.473}.
In the present work, we analyze the electronic structure without the assumption that the symmetry of the electronic structure is the same as the symmetry of the \ce{CuO2} lattice in order to shed light on the possibly lowered symmetry of the electronic structure.
Furthermore, using ARPES, one can investigate the behavior of the nematicity not only in the pseudogap state but also in the superconducting state, where transport and magnetic measurements cannot be performed due to the zero resistivity and the Meissner effect, respectively \cite{Sato2017Thermodyna,Daou2010Broken-rot,PhysRevB.92.224502}.

\section*{Results}

\begin{figure}[t]
\includegraphics[width = 8.5cm, bb = 0 0 606 835]{20181123_1SA.png}

  \caption{Energy distribution curves (EDCs) of Pb-Bi2212 at various temperatures. The measurements were performed with decreasing temperature for (A)-(D) and increasing temperature for (E)-(G). Each EDC around the X (Y) point has been integrated from $k_x = -1.1 $ to $-0.8$ (integrated from $k_y = 0.8$ to $1.1$). The difference curves between the intensities of EDCs around the X and Y points are also shown. In the pseudogap state ($T =$ 100 K and 125 K), the spectral intensities between the X and Y points are clearly different while they are identical in the normal ($T =$ 200 K and 210 K) and superconducting states ($T =$ 5 K, 25 K and 50 K). Inset shows the experimental geometry that ensures the spectroscopic equivalence of the $a$- and $b$-axis directions.}
    \label{fig:EDC}
\end{figure}

\begin{figure*}[t]
\includegraphics[width=17cm, bb = 0 0 1920 1007]{20181005_2SA.png}

\caption{\label{fig:Mapping} Constant energy-surface mapping of the ARPES spectra of Pb-Bi2212  taken at various temperatures with decreasing temperature (A)-(D) and increasing temperature (E)-(G). The constant energy is $E$ = - 30 meV and the intensity is integrated within $\pm$ 10 meV of $E$. Red circles show the constant energy surface determined from the peaks of the momentum distribution curves (MDCs). Gray curves are the fitted curves using the tight-binding model of Eq. (2).}

\end{figure*}

In Fig. \ref{fig:EDC}, energy distribution curves (EDCs) around $\vec{k} = (-\pi,0)$ (X point) and $(0,\pi)$ (Y point) at various temperatures are shown. The data were obtained from one sample but from two cleavages: We performed ARPES measurements with decreasing temperature at $T =$ 5 K, 50 K, 100 K, and 200 K on one cleaved surface and with increasing temperature at $T =$ 25 K, 125 K, and 210 K on the other cleaved surface. In both measurements, the tensile strain was applied in the $x$-direction ($||a$-axis). In the normal state (200 K), the line shapes of the EDCs were almost identical between the X and Y points, as  expected from the four-fold rotational symmetry of the \ce{CuO2} plane. With decreasing temperature, the line shapes of the EDCs around the X and Y points became different below $T^*$ in the pseudogap state (100 K) and then below $T_\text{c}$  identical again in the superconducting state (5 K and 50 K).  In the heating series, $T$ = 25 K $\rightarrow$ 125 K $\rightarrow$ 210 K, the spectral changes of the cooling series was reproduced as shown in Fig. 1.
Let us focus on the pseudogap state, where the difference of the line shapes of the EDCs at the X and Y points can be seen.  In Bi2212, it is well known that the Cu-O band is split into the anti-bonding and bonding bands due to the bilayer structure \cite{RevModPhys.75.473}. From a systematic temperature-dependent study of the dispersions near $\vec{k} = (\pi,0)$ in the overdoped sample ($T_\text{c}=$ 87 K and 91 K) whose doping level is close to and the same as the present samples, we consider that the binding energies of the bottoms of the anti-bonding band and bonding band are located around 25 meV  and 110 meV below the Fermi level, respectively, in the pseudogap state. Thus, we conclude that the intensity of the bonding (anti-bonding) band is higher around the X (Y) point than around the Y (X) point, reflecting the in-plane anisotropy of the electronic structure in the pseudogap state. Here, we would like to emphasize that this inequivalence between the X and Y points is not caused by matrix-element effect because the $x$- and $y$-directions are equivalent in our measurement geometry.

In Fig. \ref{fig:Mapping}, the intensity map of the constant energy surface at $E-E_\text{F} =$ -30 meV at various temperatures are displayed. The intensity around $k_y=-k_x$ almost vanishes due to matrix-element effect (See Supplementary Information for the matrix-element effect). From systematic photon energy-dependent-studies of the intensities from the anti-bonding and bonding bands around $\vec{k} = (\pi,0)$, the intensity of the anti-bonding band is stronger than that of the bonding band for $h\nu =$ 60 eV \cite{PhysRevLett.89.077003,PhysRevLett.90.207001} and, therefore, we believe that the constant energy surfaces shown in Fig. \ref{fig:Mapping} largely arise from the anti-bonding band.

The band structure near the Fermi level of the cuprate superconductor can be fitted using the tight-binding model
\begin{eqnarray}
 \varepsilon (k_x,k_y)- \mu = \varepsilon_0 &-&2t(\cos k_x + \cos k_y) \nonumber\\
 -4 t' \cos k_x \cos k_y &-&2t'' (\cos 2k_x + \cos 2k_y)  , \label{eq:3.2.1}
\end{eqnarray}
where $t$, $t'$, and $t''$ are the nearest-neighbor, second-nearest-neighbor, and third-nearest-neighbor hopping parameters, respectively, and $\mu$ is the chemical potential \cite{PhysRevB.95.075109}. The model has four-fold rotational ($C_4$) symmetry.
In order to examine the possibility of the rotational symmetry breaking, we introduce an anisotropy parameter $\delta$ which represents the orthrombicity of the hopping parameters ($t$ and $t''$) \cite{PhysRevB.73.214517} as
\begin{eqnarray}
&&\varepsilon (k_x,k_y) - \mu=\varepsilon_0 -2t[(1-\delta)\cos k_x +(1+\delta) \cos k_y] \nonumber\\
&& -4 t' \cos k_x \cos k_y -2t'' [(1-\delta)\cos 2k_x +(1+\delta) \cos 2k_y] . \nonumber\\ \label{eq:3.2.2}
\end{eqnarray}
By fitting the constant energy surface using this formula, we estimated its deviation from the $C_4$ symmetry through finite $\delta$ values.
As shown in Fig. \ref{fig:Mapping}, constant energy surfaces have been determined from the peak positions of the momentum distribution curves (MDCs) and fitted to the tight-binding model Eq. (\ref{eq:3.2.2}). In the same way, we also evaluated the anisotropy of the constant energy surfaces at $E-E_\text{F} =$ -10 meV and -50 meV although the MDCs are relatively broad due to the superconducting gap opening for $E - E_\text{F} =$ -10 meV and the higher scattering rate for $E - E_\text{F} =$ -50 meV. (See Supplementary Information for the energy dependence of the line width of the MDCs).

The anisotropy parameter values deduced from the series of measurement with decreasing temperature were generally smaller than those deduced from the series of measurement with increasing temperature by a factor of 0.5. Therefore, the former values have been multiplied by 0.5 and are plotted in Fig. \ref{fig:delta}. This procedure would be justified because the uniaxial strain to the sample may change after cleaving the sample depending on the degree of alignment of nematic domains under different magnitudes of the uniaxial strain. Figure 3 indicates that, in the normal state ($T > T_\text{c}$), the anisotropy parameter $\delta$ was almost zero, which was expected from the four-fold rotational symmetry of the \ce{CuO2} plane. In the pseudogap state, $\delta$ became finite in both measurements. The sign of $\delta$ indicates that the hopping parameters along the tensile strain direction became small in the pseudogap state. The sign of $\delta$ is also the same as the anisotropic Fermi surface of \ce{YBa2Cu4O8} under the strain from the Cu-O chains \cite{PhysRevB.80.100505}. Thus, the uniaxial strain served as an external perturbation to align a majority of the nematic domains in one direction (See also Supplementary Information for the effect of strain on hopping parameters).

Figure 3 also shows that $\delta$ was suppressed in the superconducting state. One may suspect that the temperature dependence of $\delta$ is inconsistent with the elongation of the spin excitation peak along the perpendicular direction to the Cu-O chain observed by the neutron scattering measurement which was monotonically enhanced with decreasing temperature even into the superconducting state \cite{Hinkov597}. However, the doping level in that work was much lower than that in the present study and a quantum phase transition was proposed at the hole concentration of  $ \simeq 0.75$ between the neutron study and the present study. Therefore, it is reasonable to consider that the nematicity observed in the present work is a different type from that detected by the neutron scattering measurement, whose onset temperature was much lower than $T^*$.
In the normal and pseudogap states, our result is consistent with the previous magnetic and transport measurements on YBCO in that the nematicity becomes finite below $T^*$\cite{Sato2017Thermodyna,Daou2010Broken-rot}. As for superconducting state, the magnetic and transport measurement do not give information, while our result may indicate competition between nematicity and superconductivity.

\section*{Discussion}

From the theoretical side, dynamical mean-field theory (DMFT) combined with the fluctuation exchange (FLEX) approximation has shown that a Pomeranchuk instability where the  anisotropy of the electronic structure becomes finite, appears in the overdoped region \cite{PhysRevB.95.075109}. The calculation showed nematicity in the overdoped regime and indicated the coexistence of nematicity and superconductivity \cite{PhysRevB.95.075109}. In contrast, according to a cellular dynamical mean-field theory (CDMFT) study, four-fold rotational symmetry breaking was shown in the underdoped pseudogap regime rather than in the overdoped regime \cite{PhysRevB.86.094522,PhysRevB.82.180511}.
According to mean-field model calculations, the Pomerachuk instability competes with superconductivity \cite{doi:10.1143/JPSJ.69.2151,PhysRevB.75.155117}. Thus, it is not clear at present whether the Pomeranchuk instability consistently explains our result or not.

\begin{figure}[t]
\includegraphics[width = 9cm, bb = 0 0 737 530]{Delta_180705_2.png
\caption{\label{fig:delta} Temperature dependence of nematicity represented by $t$. The deduced anisotropy parameter $\delta$ is defined by Eq. (2). The error bars are standard errors of the fitting parameter at each temperature.}
\end{figure}

Nematicity was derived also from fluctuations of charge-density wave (CDW) or stripes \cite{Nie03062014,PhysRevB.92.174505}. In YBCO, CDW resides inside the pseudogap phase and competes with superconductivity \cite{Ghiringhelli821,Chang2012Direct-obs,Comin1335,Le-Tacon2013Inelastic-}. Recently, CDW in Bi2212 was also observed by resonant x-ray scattering experiments and found to compete with superconductivity \cite{PhysRevB.89.220511,Silva-Neto2014Ubiquitous} as in the case of YBCO. More recent work on Bi2212 showed that the CDW signal survived up to the pseudogap temperature $T^*$ although the hole concentration of that work was in the underdoped region ($T_\text{c} =$ 40 K)  \cite{Chaix2017Dispersive}.  From the fact that CDW in Bi2212 competes with superconductivity and possibly survives up to $T^*$ including our case, it is also likely that the nematicity observed in our measurements arises from CDW fluctuations.
Although at present it is difficult to identify the microscopic origin of the nematicity, our result provides evidence that the pseudogap state shows nematicity and competes with superconductivity. The apparent competition between the nematicity and the superconductivity may possibly be caused by the $d$-wave symmetry of the superconducting gap, whose magnitude is identical between the $k_x$ and $k_y$ directions.

However, there are still unsolved issues regarding the pseudogap, that is, how the nematicity is related to the other proposed orders inside the pseudogap phase, e.g., CDW \cite{Ghiringhelli821,Chang2012Direct-obs,Comin1335,Le-Tacon2013Inelastic-,Comin390} and loop current order \cite{PhysRevLett.96.197001,PhysRevB.89.094523,PhysRevLett.118.097003,0953-8984-26-50-505701} which correspond to translational and time-reversal symmetry breaking, respectively and how the $\vec{Q}=0$ nematic order opens pseudogap. It is difficult to discuss the origin of nematicity  and further work is necessary to identify the nature of nematicity.

\section*{Materials and Methods}
Pb-Bi2212 single crystals were grown by the floating-zone method.
The hole concentration was optimimal to slightly overdoped one ($T_\text{c} =$ 91 K) after annealing the samples in a \ce{N2} flow, which allowed us to measure samples above and below $T^*$ ($\sim$160 K) rather easily compared to the underdoped samples whose $T^*$'s are too high \cite{Vishik06112012}.
We measured Pb-doped samples in order to suppress the superstructure modulation present in the BiO layers of Bi2212 and shifted by  $\vec{k} = \pm(0.21\pi ,0.21\pi)$ caused by the diffraction of photoelectrons \cite{RevModPhys.75.473}. The Umklapp replicas make the spectral intensity of the Fermi surface complicated, especially around $\vec{k} = (\pi,0)$ \cite{PhysRevB.64.180505}. \par
YBCO has Cu-O chains in addition to the \ce{CuO2} plane, which helps nematic domains to form in one particular direction, while Pb-Bi2212 has no Cu-O chains and the in-plane crystal structure is tetragonal.
Therefore, we applied a tensile strain to the sample in the Cu-O bond direction in attempt to align nematic domains in one direction using a device similar to that used for Fe-based superconductors, a picture and schematic figure of which is shown in Supplementary Fig. S2(c) \cite{PhysRevLett.115.027005}.
The strain was applied at room temperature and the stressed sample was cooled down to measurement temperatures before cleaving {\it{in situ}} (See Supplementary Information for the strain estimated from the lattice parameters by x-ray diffraction measurement). Such an operation was necessary to obtain anisotropic signals from spectroscopic data because the size of the nematic domains is not necessarily larger than the beam size; otherwise, one can measure only ``averaged isotropic'' signals \cite{Yi26042011,PhysRevB.90.121111}.\par

ARPES measurements were carried out at the undulator beamline BL5U of UVSOR using an MBS A-1 analyzer. The photon energy $h\nu$ was fixed at 60 eV. The energy resolution was set at 30 meV. The linear polarization of the incident light was chosen perpendicular to the analyzer slit and the tilt axis parallel to the analyzer slit, which realizes the unique experimental configuration that preserves the equivalence of the $a$- and $b$-axis directions with respect to the light polarization $\vec{E}$ and the strain direction, as shown in the inset of Fig. \ref{fig:EDC}. This setting ensures that the difference of the spectroscopic data between $a$- and $b$-axis directions does not originate from  matrix-element effects but from intrinsic inequivalence in the electronic structure. The measurements were performed in the normal, pseudogap, and superconducting states in two processes, that is, with increasing temperature and decreasing temperature in order to check the reproductivity. The samples were cleaved {\it in situ} under the pressure better than $2\times 10^{-8}$ Pa.
\section*{Acknowledgements}
Informative discussion with T. Shibauchi and K. Ishida is gratefully acknowledged.  A part of this work was conducted at Advanced Characterization Nanotechnology Platform of the University of Tokyo, supported by "Nanotechnology Platform" of the Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan. ARPES experiments were performed at UVSOR (Proposal Nos. 28-813 and 29-821). This work was supported by KAKENHI Grant No. 15H02109.
\section*{Author contributions} S.N., M.H., K.K., K.H., C.L., and M.S. performed ARPES measurements with the assistance of S.I. and K.T., S.N. analyzed the data. D.S., Y.Y., and E.H. synthesized and characterized single crystals. S.N. and A.F. wrote the manuscript with suggestions by M.H., D.S., and all other coauthors. A.F. was responsible for overall project direction and planning.
\section*{Competing interests} The authors declare that they have no competing interests.
\section*{Data and materials availability} All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

\input{Bi2212_nematicity.bbl}

\end{document}