Upper bounds for critical probabilities in Bernoulli Percolation models
We consider bond and site Bernoulli Percolation in both the oriented and the non-oriented cases on $\mathbb{Z}^d$ and obtain rigorous upper bounds for the critical points in those models for every dimension $d \geq 3$.
\section{Introduction}





The study of Bernoulli  percolation on $\mathbb{Z}^d$ is more than 60 years old and the existence of a non-trivial phase transition for $d\geq 2$ is well established for the model and several of its variants, but the exact value of the critical parameter $p_c$ is seldom known. A celebrated result of Kesten (see \cite{K3}) proved that the critical probability in Bernoulli bond percolation on $\mathbb{Z}^2$ is $1/2$. Beyond that, a handful of planar lattices had their critical probability established and planarity was always the key factor. On the other hand, for dimensions $d\geq 3$, there is not much hope of finding exact values for the critical probability and the best we can expect are numerical results via Monte Carlo Methods (see the very efficient algorithm in \cite{NZ}), statistical estimates based on a comparison with dependent percolation (see Section 6.2 of \cite{BR} for an overview) or rigorous bounds (see for instance \cite {W,WP}). In this paper, we focus on finding rigorous upper bounds for site and bond Bernoulli Percolation on $\mathbb{Z}^d$ for every $d\geq 3$, in both the oriented and the non-oriented cases.


The primary tools we use are couplings between the models we seek to understand and models where bounds or precise values for the critical probabilities are known. Although those bounds are still far from the values obtained from Monte Carlo simulations, we believe that most of them are the best rigorous upper bounds in the literature. 


The remainder of the text is organized as follows: in Section \ref{sec:model} we define more precisely the models and state the main results, in Section \ref{sec:couplings} we establish some dynamical couplings, in Section \ref{sec:proofs} we prove the theorems and in Section \ref{sec:comments} we give a numerical table with upper bounds for the critical probabilities in homogeneous Bernoulli percolation models for dimensions up to $d=9$.



\section{The models and main results}\label{sec:model}

For $d \geq 1$, the underlying graph for all the models will be the $d$-dimensional hypercubic lattice in which the set of vertices is $\mathbb{Z}^d$ and the set of edges is the set of non ordered pairs $E(\mathbb{Z}^d) : = \{ \langle v, u \rangle : v, u \in \mathbb{Z}^d ~ \text{and} ~ |v-u| =1\}$. We abuse notation and denote this graph simply by $\mathbb{Z}^d$. 

\subsection{Bond percolation}

Given  $p \in [0,1]$, consider a family of independent random variables $\{X_e\}_{e \in E(\mathbb{Z}^d)}$, where, for each $e \in E(\mathbb{Z}^d)$, $X_e$ has Ber($p$) distribution. Let $\mu_e$ be the law of $X_e$, and let $\mathbb{P}_p := \prod_{e \in E(\mathbb{Z}^d)} \mu_e$ be the resulting product measure. We declare an edge $e$ to be {\it open} if $X_e=1$ and {\it closed} otherwise. 

We first consider the non-oriented case and denote by $\{x \leftrightarrow y\}$ the event where  $x, y \in \mathbb{Z}^d$ are connected by an open path, i.e., there exist $x_0, \dots, x_n$ such that $x_0 = x$, $x_n =y$ and each $\langle x_{j-1}, x_j \rangle$ belongs to $ E(\mathbb{Z}^d)$ and is open for $j=1, \dots, n$. Let  $\mathcal{C}^b_0 := \{ x \in \mathbb{Z}^d:  0 \leftrightarrow x\}$ be the open cluster of the origin, and $|\mathcal{C}^b_0|$ its size.

We define the percolation probability by
$\theta^{b}_d(p) := \mathbb{P}_p(|\mathcal{C}^b_0|= \infty)$. The critical point for the {\it non-oriented bond} Bernoulli percolation model will be denoted by \[p_c^{b}(d)=\sup\{p \geq 0 : \theta^b_d(p)=0\}.\]

We now consider the oriented case. Let $\{e_1, \dots, e_d\}$ be the set of positive unit vectors of $\mathbb{Z}^d$. We denote by $\{x \rightarrow y\}$ the event where  $x, y \in \mathbb{Z}^d$ are connected by an oriented open path, i.e., there exist $x_0, \dots, x_n$ such that $x_0 = x$, $x_n =y$ and for each $j=1, \dots, n$, we have $x_j = x_{j-1} + e$, for some $e \in \{e_1, \dots, e_d\}$ and $\langle x_{j-1}, x_j\rangle$ is open. Let $\vv{\mathcal{C}}^b_0 := \{ x \in \mathbb{Z}^d : 0 \rightarrow x\}$ be the oriented open cluster of the origin, and $|\vv{\mathcal{C}}^b_0|$ its size.

Analogously, we define 
$\vv{\theta}^{\, b}_d(p) :=  \mathbb{P}_p(|\vv{\mathcal{C}}^b_0|= \infty)$ the corresponding oriented percolation probability, and we denote the critical point for the {\it oriented bond} Bernoulli percolation model by \[\vv{p}_c^{\, b}(d)=\sup\{p \geq 0 : \vv{\theta}^{\,b}_d(p)=0\}.\]

Our results are the following:

\begin{theorem}\label{theo:bounds} Consider non-oriented bond Bernoulli percolation and 
let $p^{\ast}(d)$ be the unique solution in $(0,1)$ of 
\[\prod_{i=0}^2 \left(1-(1-p)^{\lfloor\frac{d+i}{3}\rfloor}\right)=2-\sum_{i=0}^2(1-p)^{\lfloor\frac{d+i}{3}\rfloor}.\]

Then, for every $d\geq 3$, we have $p_c(d) \leq p^{\ast}(d)$.
\end{theorem}

\begin{theorem}\label{theo:bondO}
Consider oriented bond Bernoulli percolation on $\mathbb{Z}^d$. 

1) If $d$ is even, then $ \vv{p}^{\,b}_c(d)\leq 1-(1/3)^{2/d}$;


2) For $d \geq 4$, we have that \[\vv{p}^{\,b}_c(d)\leq \frac{1}{d}+\frac{C_d}{d^2},\]

where \[C_d=1+\frac{8}{d}+\frac{d^{5/2}}{(\sqrt{2\pi})^{d-1}}\left[\frac{d-1}{d-3}\right]e^{\frac{1}{12d}}.\]

3) For any dimension $d\geq 2$, we have that  $\vv{p}^{\,b}_c(d+1)\leq f(d)$, where $f(d)$ is the unique solution  in $(0,1)$ of 
\[ p = \vv{p}^{\,b}_c(d) \left[p + (1-p)^{(d+1)/d}\right]. \]

\end{theorem}

{\bf Remark:} It is known (see \cite{CD}) that $\vv{p}^{\,b}_c(d)\sim 1/d$, hence the third upper bound above is asymptotically sharp.






\subsection{Site percolation}

Given a parameter $p \in [0,1]$, we consider a family $\{X_v\}_{v \in \mathbb{Z}^d}$ of independent Bernoulli random variables with parameter $p$. As before, $\mathbb{P}_p$ will denote the resulting product measure. A vertex $v \in \mathbb{Z}^d$ is declared to be {\it open} if  $X_v = 1$ and {\it closed} otherwise. 

Now, the sequence of vertices $(x_0, \dots, x_n)$ is said to be an open path if all the vertices are open and for each $j = 1, \dots, n$, $\langle x_{j-1}, x_j\rangle \in E(\mathbb{Z}^d)$. If in addition to these conditions, for each $j = 1, \dots, n$, we have $x_j = x_{j-1}+ e$ for some $e \in \{e_1, \dots, e_d\}$, the sequence is said to be an oriented open path. Having these definitions we  define $\mathcal{C}^s_0$, $\vv{\mathcal{C}}^s_0$, $\theta_d^{\, s}(p)$ and $\vv{\theta}_d^{\, s}(p)$ accordingly.

Finally, the critical points for {\it non-oriented} and {\it oriented} site Bernoulli percolation are respectively given by
\[ p_c^{s}(d)=\sup\{p \geq 0 : \theta^s_d(p)=0\} \quad \quad \text{and} \quad \quad \vv{p}_c^{\, s}(d)=\sup\{p \geq 0 : \vv{\theta}^{\,s}_d(p)=0\} .\]

For site percolation models, our results are the following

\begin{theorem}\label{theo:siteNO}
Consider non-oriented site Bernoulli percolation on $\mathbb{Z}^d$. 

1) If $d$ is even, then $ p_c^{s}(d)\leq 1-(0,32)^{2/d}$; 

2) If $d$ is divisible by 3, then $p_c^{s}(d)\leq 1-(0,5)^{3/d}$;

3) For any dimension $d\geq 2$, we have that $p_c^{s}(d+1)\leq g(d)$, where $g(d)$ is the unique solution in $(0,1)$ of 
\[ p = p^{\,s}_c(d) \left[p + (1-p)^{2d/(2d-1)}\right]. \]

\end{theorem}

\begin{theorem}\label{theo:siteO}
Consider oriented site Bernoulli percolation on $\mathbb{Z}^d$. 

1) If $d$ is even, then $ \vv{p}_c^{\,s}(d) \leq 1-(0,25)^{2/d}$;

2) For any dimension $d \geq 2$, we have that $\vv{p}_c^{\,s}(d+1) \leq h(d)$, where $h(d)$ is the unique solution in $(0,1)$ of
\[p = \vv{p}^{\,s}_c(d) \left[p + (1-p)^{(d+1)/d}\right]. \]


\end{theorem}






\section{The Dynamical Couplings}\label{sec:couplings}



Although we are only interested in homogeneous percolation, some of the tools we will use are related to couplings between {\it anisotropic} bond percolation models.

We now define the anisotropic (or inhomogeneous) non-oriented and oriented bond percolation models. For each $i=1, \dots, d$, let $E_i = \{ \langle x, x + e_i \rangle: x \in \mathbb{Z}^d\}$ be the set of edges parallel to $e_i$. Given $p_1, \dots, p_d \in [0,1]$,  we consider a family of independent random variables $\{X_e\}_{e \in E(\mathbb{Z}^d)}$, but now, for each $e \in E_i$, $X_e$ has Ber($p_i$) distribution, $i=1, \dots, n$. The open cluster of the origin $\mathcal{C}^b_0$ and the oriented open cluster of the origin $\vv{\mathcal{C}}^b_0$ are  defined analogously. The probabilities that $|\mathcal{C}^b_0|$, $|\vv{\mathcal{C}}^b_0|$ are infinite, will be denoted by $\theta_d(p_1, \dots, p_d)$ and $\vv{\theta}_d (p_1, \dots, p_d)$ respectively.



The first coupling is the content of Proposition 1 in \cite{GPS2}:
\begin{proposition} \label{prop:coupling} 
Consider inhomogeneous non-oriented Bernoulli bond percolation in $\mathbb{Z}^d$. Let $p_1, \dots, p_{d+1} \in [0,1]$ and let $\Tilde{p}_d \in [0,1]$ be such that \[(1-\Tilde{p}_d) = (1-p_d)(1-p_{d+1}).\] Then,
    $\theta_{d+1}(p_1, \dots,p_d,  p_{d+1})  \geq \theta_{d}(p_1, \dots, p_{d-1}, \Tilde{p}_{d}) .$
\end{proposition}

Now we define bond and site percolation models in the triangular lattice $\mathbf{T}$. This lattice is simply $\mathbb{Z}^2$ with  extra edges of the form $\langle v,v+(1,1)\rangle$. That is, $\mathbf{T} = (V_{\mathbf{T}}, E_{\mathbf{T}})$ where $V_{\mathbf{T}} = \mathbb{Z}^2$ and $E_{\mathbf{T}}$ is the set of non ordered pairs $\{\langle v, u\rangle : v-u = (1,0), (0,1) \text{ or } (1,1) \}$. We will denote by $\theta^b_{\mathbf{T}}$ and $\theta^s_{\mathbf{T}}$ the corresponding percolation probability for bond and site models, respectively.

We will consider inhomogeneous bond percolation on $\mathbf{T}$, where the corresponding parameters $(p_1,p_2,p_3)$ will refer to edges of the form $\langle v,v+(1,0)\rangle$, $\langle v,v+(0,1)\rangle$ and   $\langle v,v+(1,1)\rangle$, respectively. 







In the next proposition, we construct a monotonic coupling between inhomogeneous bond percolation on the  triangular lattice $\mathbf{T}$ with parameters $(p_1,p_2,p_3)$ and on $\mathbb{Z}^3$ with the same parameters.

\begin{proposition}\label{prop:triangular}
Let  $(p_1,p_2,p_3)\in [0,1]$ and consider two inhomogeneous bond Bernoulli percolation processes on the triangular lattice  and on $\mathbb{Z}^3$, both with parameters $(p_1,p_2,p_3)$. Then
\[\theta^b_{\mathbf{T}}(p_1,p_2,p_3)\leq\theta^b_3(p_1,p_2,p_3).\]
\end{proposition}


\begin{proof}

We will construct a dynamic coupling between the percolation process on $\mathbb{Z}^{3}$ with parameters $(p_1,p_2, p_3)$  and an infection process over $\mathbf{T}$.  We will do it in such a way that the law of infected sites in $\mathbf{T}$ is the same as the law of the open cluster of the origin for anisotropic percolation on $\mathbb{Z}^3$ and also that, if the infection process survives, the open cluster of the origin of the process in $\mathbb{Z}^3$ must be infinite. 

The coupling will be built based on a susceptible-infected strategy described as follows. First, we declare  the origin of $\mathbf{T}$ as the \textit{initial} infected component. Next, at each time-step, we possibly grow the infected component and associate each new vertex $v$ of the infected component in $\mathbf{T}$ to a vertex $x(v)$  in the open cluster of the origin in $\mathbb{Z}^3$.  More precisely, consider a vertex $v$ in the infected component of $\mathbf{T}$ and a neighbor $v+u$ out of the infected component. The vertex $v$, in the infected component in $\mathbf{T}$, will be associated with some vertex $x(v)$ in the  open cluster of the origin in $\mathbb{Z}^3$. According to $u = \pm(1,0) , \pm(0,1)$ or $\pm(1,1)$ we denote $\tau(u) = \pm(1,0,0)$, $\tau(u) = \pm(0,1,0)$ or $\tau(u) = \pm(0,0,1)$ respectively. If $\langle x(v), x(v)+\tau(u) \rangle$,  is open, we infect $v+u$. (and write $x(v+u) = x(v) + \tau(u)$).


Let's define the sequence of sets $(I_n, x(I_n), R_n, S_n)_{n \geq 0}$. Here, $I_n$ represents the \textit{infected vertices} in $\mathbf{T}$ and $x(I_n)$ represents the vertices in $\mathbb{Z}^{3}$ associated with the infected vertices. $R_n$ represents the \textit{removed edges} of $\mathbf{T}$. Finally, given $I_n$, $x(I_n)$ and $R_n$, the \textit{susceptible edges set} is given by
\begin{equation*}
    S_{n} := \{ \langle v, u \rangle : v \in I_{n} \: \mbox{and} \: u \notin I_{n} \} \cap R_{n}^C.
\end{equation*}
At time $n=0$, we set
\begin{itemize}
    \item $I_0 = \{0\} \subset V_{\mathbf{T}}$; 
    
    \item $R_0 = \emptyset \subset E_{\mathbf{T}}$;
    
    \item $x(0) = 0 \in \mathbb{Z}^{3}$;
    
    \item $S_0 := \{ \langle v, u \rangle : v \in I_{0} \: \mbox{and} \: u \notin I_{0} \} \cap R_{0}^C \:\: = \{ \langle 0, \pm (1,0)\rangle, \langle 0,\pm (1,0)\rangle ,\langle 0,\pm (1,1)\rangle \}$. 
\end{itemize}

This means that, at time $n=0$, only the vertex $0$ is infected, and it can potentially infect any of its neighbours, so all edges containing the origin are susceptible. After that, in each step, an infected vertex tries to infect a non-infected vertex through a susceptible edge (if the latter exists). From now on, we choose an arbitrary, but fixed, ordering of the edges in $\mathbf{T}$. Suppose that $I_n, x(I_n), R_n$ and $S_n$ are already defined. If there is no susceptible edge then the process stops. More specifically, if $S_n = \emptyset$, then for all $k \geq 1$,
    \begin{itemize}
        \item $I_{n+k} = I_n$;
        
        \item $R_{n+k} = R_n$;
        
        \item $S_{n+k} = S_n$.
    \end{itemize}
Otherwise, if there exists some susceptible edge, then the infected vertex in the smallest (in the previously fixed ordering) such edge tries to infect its non-infected neighbour. Let us write it in symbols. Suppose $S_n \neq \emptyset$ and let $g_n$ be the smallest edge in $S_n$.
Since $g_n \in S_n$, it has to be equal to some $\langle v, v+u_n \rangle$, where $v \in I_n$, $v+u_n \notin I_n$, and $u_n \in \{ \pm (0,1), \pm (1,0), \pm (1,1)\}$. 
We set $\tau(\pm(1,0))=\pm (1,0,0)$, $\tau(\pm(0,1))=\pm (0,1,0)$ and $\tau(\pm(1,1))=\pm (0,0,1)$ .
Then, $v$ \textit{infects} $v+u_n$ if $\langle x(v), x(v)+\tau (u_n) \rangle$ is open in $\mathbb{Z}^{3}$. More precisely, if $\langle x(v), x(v)+\tau (u_n) \rangle$ is open in $\mathbb{Z}^{3}$ then we write

\[ I_{n+1} := I_n \cup \{v+u_n\},\]
and define
\[ x(v+u_n) := x(v) + \tau(u_n). \] 
Otherwise, if $\langle x(v), x(v)+ \tau(u_n) \rangle$ is closed in $\mathbb{Z}^{3}$,  we set $I_{n+1}:= I_n$.



Now that we have explored $g_n$, we \textit{remove} it and write
\[ R_{n+1} := R_n \cup \{ g_n \}.\]
Next, to conclude our induction step, we set
\begin{equation*}
  S_{n+1} := \{ \langle v, u \rangle : v \in I_{n+1} \: \mbox{and} \: u \notin I_{n+1} \} \cap R_{n+1}^C.
\end{equation*}
Observe that the function $x: \cup_n I_n \longrightarrow \mathbb{Z}^{3}$ is injective.  
In fact, if $v= (v_1, v_2) \in I_n$, then by construction, we have that $x(v) = (x_1,x_2,x_3)$ satisfies
\[ v_1=x_1+x_3 \mbox{ and }
    v_2=x_2+x_3.\]
Now, observe that the image of $x$ is contained in the  open cluster of the origin $\mathcal{C}^b_0$ of $\mathbb{Z}^{3}$. Since $x$ is injective,  $|\cup_n I_n| \leq |\mathcal{C}^b_0|$.  Also, note that $\cup_{n} I_n$ has the same  law as  $\mathcal{C}^{\mathbf{T}}_0$, where $\mathcal{C}^{\mathbf{T}}_0$ is the  open cluster of the origin in $\mathbf{T}$ with parameters $p_1, p_2,p_3$.  Therefore,
\begin{equation*}
 \theta^b_{\mathbf{T}}(p_1,p_2,p_3) \leq \theta_{3}^b (p_1,p_2,p_3),
\end{equation*}
and the proof of Proposition \ref{prop:triangular} follows.
\end{proof}

To conclude this section, we construct two couplings that will be used to prove Item 3 of Theorem \ref{theo:bondO}, along with its site version, which will be used in the proofs of Item 3  of Theorem \ref{theo:siteNO} and Item 2 of Theorem \ref{theo:siteO}. These couplings are reminiscent of the coupling in \cite{GSS} and give an upper bound for the critical probability for any percolation model in $\mathbb{Z}^{d+1}$ as a function of the corresponding critical probability in $\mathbb{Z}^d$.  Since our goal is to obtain better upper bounds, we need to improve the coupling.

\begin{proposition}\label{prop:crossover}
Consider oriented bond Bernoulli percolation on $\mathbb{Z}^{d+1}$ and suppose that 
\begin{equation} \label{eq:crossbond}
    p > \vv{p}^{\,b}_c(d) \left[p + (1-p)^{(d+1)/d}\right].
\end{equation}
Then, $\vv{\theta}^{\,b}_{d+1}(p)>0$.
\end{proposition}


\begin{proof} 
Let ${\mathcal{E}} = \{e_1, \dots, e_d\}$ denote the set of positive unit vectors of $\mathbb{Z}^d$. To avoid ambiguities, we denote the set of positive unit vectors of $\mathbb{Z}^{d+1}$ by $\{ u_1, \dots, u_{d+1}\}$. Consider the multigraph $\mathbb{Z}^{d+1}_{\mathcal{E}}$ defined as follows: the set of vertices is $\mathbb{Z}^{d+1}$ and the set of edges is given by $E(\mathbb{Z}^{d+1}_{\mathcal{E}}):= \left( \cup_{i=1}^d E_i \right) \cup E_{\mathcal{E}}$ where  $E_i = \{ \langle x, x + u_i \rangle: x \in \mathbb{Z}^{d+1}\}$, $i=1, \dots, d$, and we define $E_{\mathcal{E}} := \{\langle v, v+u_{d+1}\rangle_e : v \in \mathbb{Z}^{d+1} , e \in {\mathcal{E}}\}$. In words, each edge of $\mathbb{Z}^{d+1}$ parallel to $u_{d+1}$ is partitioned into another $|{\mathcal{E}}|$ edges indexed by ${\mathcal{E}}$ in $\mathbb{Z}^{d+1}_{\mathcal{E}}$, while edges parallel to all other directions remain unmodified.

We prove this proposition in two steps. First, we will define a multigraph $\mathbb{Z}^{d+1}_{\mathcal{E}}$ and show that, with certain parameters, the model on $\mathbb{Z}^{d+1}_{\mathcal{E}}$ is equivalent to the homogeneous model on $\mathbb{Z}^{d+1}$ with parameter $p$. To conclude, we will show that if $p$ satisfies Inequality \eqref{eq:crossbond}, then this new model on $\mathbb{Z}^{d+1}_{\mathcal{E}}$ dominates a supercritical model on $\mathbb{Z}^{d}$.

Consider now inhomogeneous oriented bond Bernoulli percolation on $\mathbb{Z}^{d+1}_{\mathcal{E}}$, where edges in $\cup_{i=1}^d E_i$ are open with probability $p$ and edges in $E_{\mathcal{E}}$ are open with probability $q$, where $q$ is such that $(1-p) = (1-q)^{|{\mathcal{E}}|}$. Clearly, the distribution of the  open cluster in this model is the same as in the homogeneous model on $\mathbb{Z}^{d+1}$ with parameter $p$.  We will construct a coupling showing that the model on $\mathbb{Z}_{\mathcal{E}}^{d+1}$ with parameters $p$ and $q$ as above, dominates the homogeneous model on $\mathbb{Z}^d$ with parameter $p/[1 - (1-p)q]$. 

First, for each $e_i \in \mathcal{E}$, let $ \sigma(e_i):= u_i \in \mathbb{Z}^{d+1}$. Then,  for each $v \in \mathbb{Z}^{d+1}$ and each $e \in {\mathcal{E}}$, let $A(v,e)$ be the event where either the edge $\langle v, v+\sigma(e) \rangle$ is open or, for some $k \geq 1$,
\begin{itemize}
    \item the edges $\langle v+ iu_{d+1}, v + (i+1)u_{d+1}  \rangle_e$, $i = 0, \dots, k-1$, are open and
    \item the edges $\langle v + iu_{d+1}, v + iu_{d+1} + \sigma(e) \rangle$, $i = 0, \dots, k-1$ are closed, and 
    \item $\langle v + ku_{d+1}, v + ku_{d+1} + \sigma(e) \rangle$ is open.
\end{itemize}
In the event  $A(v,e)$, we define $u(v,e) = v+\sigma(e)$ if $\langle v, v+\sigma(e) \rangle$ is open, or $u(v,e) = v + k u_{d+1} + \sigma(e)$ if $k \geq 1$ is such that the above three conditions are met. Observe  that
\begin{eqnarray*}
    \mathbb{P}(A(v,e)) &=&
    p\sum_{i = 0}^{\infty} \left[q(1-p)\right]^i
    = \frac{p}{1 - (1-p)q}\\
    &=& \frac{p}{1-(1-p)(1-(1-p)^{1/d})} \\
    &=& \frac{p}{ p + (1-p)^{(d+1)/d} }.
\end{eqnarray*}

Similarly to what was done in Proposition \ref{prop:triangular}, we now build the sequence of sets $(I_n, x(I_n), R_n, S_n)_{n \geq 0}$. For $n =0$, we set
\begin{itemize}
    \item $I_0 = \{0\} \subset \mathbb{Z}^d$;
    \item $R_0 = \emptyset \subset E(\mathbb{Z}^d)$;
    \item $x(0) = 0 \in \mathbb{Z}_{\mathcal{E}}^{d+1}.$
\end{itemize}
Suppose that, for some $n \geq 0$, $I_n, x(I_n)$ and $R_n$ are already defined. Then $S_n \subset E(\mathbb{Z}^d)$ is given by
\begin{equation*
    S_{n} := \{ \langle v, u \rangle : v \in I_{n} \: \mbox{and} \: u \notin I_{n} \} \cap R_{n}^C.
\end{equation*}
If $S_n = \emptyset$, our sequence becomes constant, that is, 
\begin{equation}\label{eq:stop}
    (I_{n+k}, x(I_{n+k}), R_{n+k}, S_{n+k}) = (I_n, x(I_n), R_n, S_n), \quad \forall k \geq 1.
\end{equation} 
Otherwise, let $g_n = \langle v , v + e \rangle$ be the smallest (according to a prefixed ordering, as in Proposition~\ref{prop:triangular}) edge of $S_n$, where $v \in I_n$, $e \in {\mathcal{E}}$, and $v+e \notin I_n$. We set $R_{n+1} = R_n \cup \{g_n\}$. We also set
\begin{equation*}
    I_{n+1} = 
    \begin{cases} I_n \cup \{v+e\}, &\text{if } A(x(v),e) \text{ occurs};\\
                  I_n, & \text{otherwise}.
    \end{cases}
\end{equation*}
In case  $A(x(v),e)$ occurs, we set $x(v+e) = u(x(v),e)$. 

Once our sequence $(I_n, x(I_n), R_n, S_n)_ {n \geq 0}$ is built, note that by construction, the function $x: \cup_n I_n \longrightarrow \mathbb{Z}^{d+1}_{\mathcal{E}}$ is injective. In fact, for each $v \in \cup_n I_n$, the projection of $x(v)$ into $\mathbb{Z}^d$ is equal to $v$. The conclusion of the proof follows in a similar way to that of Proposition \ref{prop:triangular}.
\end{proof}

\begin{proposition}\label{prop:crossoversite}
Consider non-oriented site Bernoulli percolation on $\mathbb{Z}^{d+1}$ and suppose that 
\begin{equation}\label{eq:cross1}
p > p^s_c(d)  \left[p + (1-p)^{2d/(2d-1)}\right].
\end{equation}
Then $\theta^s_{d+1}(p)>0$.
\end{proposition}


\begin{proof}
The strategy of the proof is similar to the previous proposition. As before, let $\{e_1, \dots, e_d\}$ denote the set of positive unit vectors of $\mathbb{Z}^d$ and let $\{ u_1, \dots, u_{d+1}\}$ denote the set of positive unit vectors of $\mathbb{Z}^{d+1}$. We will consider a graph where each vertex of $\mathbb{Z}^{d+1}$ will be partitioned into another $2d-1$ vertices.

For each $v \in \mathbb{Z}^{d+1}$, we define the set of split vertices $V_v := \{v^{(1)}, \dots, v^{(2d-1)}\}$. Let $\mathbb{Z}^{d+1}_V$ be the graph with vertex set  $\cup_{v \in \mathbb{Z}^{d+1}} V_v$ and  edge set
\begin{equation*}
    E(\mathbb{Z}^{d+1}_V) := \left\{ \langle x,y \rangle ~:~ x \in V_v \text{ and } y \in V_u, \text{ for some } \langle v,u \rangle \in E(\mathbb{Z}^{d+1}) \right\}.
\end{equation*}
Given $p$ satisfying \eqref{eq:cross1}, let $q$ be such that $(1-p) = (1-q)^{2d-1}$. We will consider the non-oriented site Bernoulli percolation model on $\mathbb{Z}^{d+1}_V$ with parameter $q$. Note that the distribution of the open cluster of the origin in this new model is the same as in the model on $\mathbb{Z}^{d+1}$ with parameter $p$.

To build the coupling between $\mathbb{Z}^{d+1}_V$ with parameter $q$ and $\mathbb{Z}^d$ with parameter $p/[1  - (1-p)q]$, we again construct a sequence of sets $(I_n, x(I_n), R_n, S_n)_{n \geq 0}$. Since we are  considering a site model, for each $n \geq 0$, the sets $R_n$ and $S_n$ will be, respectively, the sets of {\it removed} and {\it susceptible} vertices (instead of edges) at step $n$. We start with two infected vertices (otherwise the origin should have been split into $2d$ vertices, instead of $2d-1$). For $n = 0$, we set
\begin{itemize}
    \item $I_0 = \{0, e_1\} \subset \mathbb{Z}^d$;
    \item $R_0 = \emptyset \subset \mathbb{Z}^d$;
    \item $x(0) = 0^{(1)} \in \mathbb{Z}_V^{d+1} \text{ and } x(e_1) = u_1^{(1)} \in \mathbb{Z}_V^{d+1}. $
\end{itemize}
Inductively, if for some $n \geq 0$, the sets $I_n, x(I_n)$ and $R_n$ are already defined, we set 
\begin{equation}\label{eq:SnSITE}
    S_n = \{u \in \mathbb{Z}^d : \langle v,u \rangle \in E(\mathbb{Z}^d) \text{ for some } v \in I_n\} \cap (I_n \cup R_n)^C.
\end{equation}
If $S_n = \emptyset$, the sequence becomes constant, as in \eqref{eq:stop}. Otherwise, let $a_n$ be the smallest (in a preset ordering) vertex of $S_n$. We can write $a_n =  v+ e$, such that $v \in I_n$ and $e \in \{\pm e_1, \dots, \pm e_d\}$. In this case, let $j_n$ denote the number of susceptible neighbors of $v$ (including $a_n$), that is
\begin{equation*}
    j_n := \left\vert \{ u \in S_n : \langle v, u \rangle \in  E(\mathbb{Z}^d)  \} \right\vert .
\end{equation*}
Since we start with two infected vertices, necessarily $|S_n| \neq 2d$, and then $1 \leq j_n \leq 2d-1$. 


For each $i \in \{1, \dots, d\}$, we set the notation $\sigma(\pm e_i) = \pm u_i$.

Consider the event $A_n$ where, either, for some $\ell \in \{1, \dots, 2d-1 \}$, the vertex $(x(v)+\sigma(e))^{(\ell)}$ is open or, for some $k \geq 1$ and $\ell \in \{1, \dots, 2d-1\}$,
\begin{itemize}
    \item the vertices $\left( x(v) + iu_{d+1} \right)^{(j_n)}$, $i = 1, \dots, k$, are open, and
    \item the vertices $\left( x(v) + i u_{d+1} + \sigma(e) \right)^{(j)}$,  are closed $\forall j \in \{1, \dots, 2d-1 \}$, $i = 0, \dots. k-1$, and
    \item the vertex $\left( x(v) + ku_{d+1} + \sigma(e) \right)^{(\ell)}$ is open.
\end{itemize}
Note that  $A_n$ has probability   
\begin{eqnarray*}
    \mathbb{P}(A_n) &=& \left(1 - (1-q)^{2d-1}\right) \sum_{i=0}^{\infty} \left[q(1-q)^{2d-1}\right]^i \\
    &=& \frac{p}{1 - (1-p)q}\\
    &=& \frac{p}{1-(1-p)\left(1-(1-p)^{1/(2d-1)}\right)} \\
    &=& \frac{p}{ p + (1-p)^{2d/2d-1} }.
\end{eqnarray*}

To conclude our induction step, we set 
\begin{itemize}
    \item $R_{n+1} = 
    \begin{cases} R_n \cup \{v+e\}, &\text{if } A_n \text{ does not occur};\\
                  R_n, & \text{otherwise};
    \end{cases}$
    \item $I_{n+1} = 
    \begin{cases} I_n \cup \{v+e\}, &\text{if } A_n \text{ occurs};\\
                  I_n, & \text{otherwise}.
    \end{cases}$
\end{itemize}
In the event $A_n$, if for some $\ell \in \{1, \dots, 2d-1 \}$ the vertex $(x(v)+\sigma(e))^{(\ell)}$ is open, we set $x(v+e) =( x(v)+\sigma(e) )^{(\ell)}$. Otherwise, let $k \geq 1$ and $\ell \in \{1, \dots, 2d-1\}$ be such that the above three conditions are satisfied, in this case, we set $x(v+e) =( x(v)+ k u_{d+1} + \sigma(e) )^{(\ell)}$.

By construction, for each $v \in \cup_n I_n$, the projection of $x(v)$ onto $\mathbb{Z}^d$ is equal to $v$. Therefore $x: \mathbb{Z}^d \longrightarrow \mathbb{Z}^{d+1}_V$ is also injective. It follows that the site percolation model in $\mathbb{Z}^{d+1}_V$ (with $e_1^{(1)}$ declared open), dominates the infection process $(I_n)_{n \geq 0}$. 
Since  $|\cup_n I_n|$ has the same distribution as the size of the open cluster of $\{0, e_1\}$ in a supercritical percolation model on $\mathbb{Z}^d$, the proof follows.

\end{proof}





\section{Proof of Theorems}\label{sec:proofs}
In this section, we prove the theorems using the couplings of the last section

\subsection{Proof of Theorem \ref{theo:bounds}}



\begin{proof}
Given Propositions \ref{prop:coupling} and \ref{prop:triangular}, the proof of Theorem \ref{theo:bounds} is quite straightforward. Recall that for bond percolation on the triangular lattice with parameters $(p_1,p_2,p_3)$, assuming $p_1,p_2,p_3<1$, we have (see Theorem 11.116 in \cite{Grim})
\begin{equation*}
\theta_{\mathbf{T}}(p_1,p_2,p_3)>0 \Leftrightarrow p_1+p_2+p_3-p_1p_2p_3>1. 
\end{equation*}

Consider now the partition
\begin{equation*}
    d=\left\lfloor\frac{d}{3}\right\rfloor+\left\lfloor\frac{d+1}{3}\right\rfloor+\left\lfloor\frac{d+2}{3}\right\rfloor.
\end{equation*}

Let $p>p^{\ast}(d)$ and set, for $i=1,2,3$,
\begin{equation*}
    p_i=1-(1-p)^{\lfloor \frac{d+i-1}{3}\rfloor}.
\end{equation*} 
Then, iteratively applying Proposition \ref{prop:coupling} and finally Proposition \ref{prop:triangular}, we have
\begin{equation*}
    \theta_d(p)\geq\theta_3(p_1,p_2,p_3)\geq \theta_{\mathbf{T}}(p_1,p_2,p_3)>0,
\end{equation*}
which concludes the proof.
\end{proof}



\subsection{Proof of Theorem \ref{theo:bondO}}
\begin{proof}
First, we remark that  Proposition \ref{prop:coupling} holds {\it mutatis mutandis} for inhomogeneous oriented bond Bernoulli percolation. As a consequence, we have the following corollary which we state without proof.
\begin{lemma}\label{lem:Cor}
Let $k \in \mathbb{N}$ be such that $d$ is divisible by $k$. Given $p \in [0,1]$, let $\tilde{p}$ be such that 
\[(1-\tilde{p}) = (1-p)^{d/k}.\]
Then,  $\vv{\theta}^{\, b}_d(p) \geq \vv{\theta}_k^{\, b}(\tilde{p})$.
\end{lemma}

Taking $k=2$ in Lemma \ref{lem:Cor} and using  Liggett's upper bound $\vv{p}_c^{\, b}(2) \leq 2/3$ (see \cite{L}), the first item of Theorem \ref{theo:bondO} follows. The third item is equivalent to Proposition \ref{prop:crossover}. 

 Finally, the second item is implicitly proved in \cite{GPS}. 
 There, the authors defined a quantity $\lambda(1/d, \dots, 1/d)$ and showed that (see Equation (3.2)), if $p \in [0,1]$ satisfies
\[\phi(d) :=  \lambda(1/d, \dots, 1/d) \leq dp -1,\]
then $\vv{\theta}^{\, b}_d(p) > 0$. In particular, 
\[ \vv{p}_c^{\, b}(d) \leq \frac{1+\phi(d)}{d}.\]
The conclusion follows from the estimate given in Subsection 3.1 of \cite{GPS} (see the equation above (3.5)), that is
\[\phi(d) \leq \frac{1}{d}+\frac{8}{d^2}+\frac{d^{3/2}}{(\sqrt{2\pi})^{d-1}}\left[\frac{d-1}{d-3}\right]e^{\frac{1}{12d}}. \]
\begin{comment}
 Let $\{S_n^1\}_{n \geq 0}$ and $\{S_n^2\}_{n \geq 0}$ be two independent, oriented, simple and symmetrical random walks on $\mathbb{Z}^d$, starting on the origin. In the paper \cite{GPS} was shown that, if $p \in [0,1]$ is such that
\[\phi(d) := \sum_{n=1}^{\infty} \mathbb{P}\left( S_n^1 = S_n^2 \right)  \leq dp -1,\]
then $\vv{\theta}^{\, b}_d(p) > 0$. In particular, 
\[ \vv{p}_c^{\, b}(d) \leq \frac{1+\phi(d)}{d}.\]
The conclusion follows from the estimate given in Subsection 3.1 of \cite{GPS} (see the equation above 3.5),
\[\phi(d) \leq \frac{1}{d}+\frac{8}{d^2}+\frac{d^{3/2}}{(\sqrt{2\pi})^{d-1}}\left[\frac{d-1}{d-3}\right]e^{\frac{1}{12d}}. \]
\end{comment}
\end{proof}






\subsection{Proof of Theorem \ref{theo:siteNO}}
Inhomogeneous percolation is not well defined for site models, and therefore, we are not able to  formulate a version of Proposition \ref{prop:coupling} for site percolation. But note that Lemma \ref{lem:Cor} only involves homogeneous models. In the following, we give its site version.

\begin{lemma}\label{lem:ksite}
Consider non-oriented site Bernoulli percolation on $\mathbb{Z}^{d}$. Let $k \in \mathbb{N}$ be such that $d$ is divisible by $k$. Given $p \in [0,1]$, let $\tilde{p}$ be such that 
\[(1-\tilde{p}) = (1-p)^{d/k}.\]
Then, we have ${\theta}^s_d(p) \geq {\theta}^s_k(\tilde{p})$.
\end{lemma} 

\begin{proof}
The lemma follows by a coupling between non-oriented site Bernoulli percolation on $\mathbb{Z}^d$, with parameter $p$, and on $\mathbb{Z}^{k}$, with parameter $\tilde{p}$. To establish this coupling, we again construct a sequence of vertex sets $(I_n, x(I_n), R_n, S_n)_{n \geq 0}$.  

First, we recall that $\mathcal{E} = \{e_1, \dots, e_d\}$ is the set of positive unit vectors of $\mathbb{Z}^d$ and let $(D_{u_1}, \dots, D_{u_k})$ be a uniform partition of $\mathcal{E}$ into $k$ subsets indexed by the set $\{u_1, \dots, u_k \}\subset\mathbb{Z}^k$ of positive unit vectors of $\mathbb{Z}^k$. For each $u \in \{u_1, \dots, u_k \} $, let $D_{-u} = -D_u$.

Consider the non-oriented site Bernoulli percolation model on $\mathbb{Z}^d$, with parameter $p$. For each $v \in \mathbb{Z}^d$ and each $u \in \{\pm u_1, \dots, \pm u_k\}$, we define the following event
\begin{equation*}
    B(v,u) := \{ v+e \text{ is open, for some } e \in D_u \}.
\end{equation*}
If the event $B(v,u)$ occurs, we set $e(v,u) = v+e$, where $e$ is an open vertex that guarantee the occurrence of $B(v,u)$. Note that  $\mathbb{P}(B(v,u)) = \tilde{p}.$

For $n = 0$, we define 
\begin{itemize}
    \item $I_0 = \{0\} \subset \mathbb{Z}^k$;
    \item $R_0 = \emptyset \subset \mathbb{Z}^k$;
    \item $x(0) = 0 \in \mathbb{Z}^d$.
\end{itemize}
If for some $n \geq 0$, the sets $I_n, x(I_n)$ and $R_n$ are already defined, let $S_n$ be as given in \eqref{eq:SnSITE}. If $S_n = \emptyset$, the sequence becomes constant as in \eqref{eq:stop}. Otherwise, let $v_n$ be the smallest vertex of $S_n$ (according to a preset ordering of $\mathbb{Z}^k$). In this case, we write $v_n = v+u \notin I_n$, where $v \in I_n$ and $u \in \{\pm u_1, \dots, \pm u_k\}$. To conclude the induction step, we define
\begin{itemize}
    \item $R_{n+1} = 
    \begin{cases} R_n \cup \{v+u\}, &\text{if } B(x(v),u) \text{ does not occurs};\\
                  R_n, & \text{otherwise};
    \end{cases}$
    \item $I_{n+1} = 
    \begin{cases} I_n \cup \{v+u\}, &\text{if } B(x(v),u) \text{ occurs};\\
                  I_n, & \text{otherwise}.
    \end{cases}$
\end{itemize}
If $B(x(v),u)$ occurs, we define $x(v+u) = e(x(v),u)$.

Note  that, by construction, the function $x: \cup_n I_n \longrightarrow \mathbb{Z}^d$ is such that, for each $v = v_{u_1}u_1+ \cdots + v_{u_k}u_k \in \cup_n I_n$, writing $x(v) = x_{e_1} e_1 + \cdots + x_{e_d} e_d$, we have
\begin{equation*}
    \sum_{e \in D_u} x_e = v_u, \quad \forall u \in \{u_1, \dots, u_k\}.
\end{equation*}
Therefore, $x: \cup_n I_n \longrightarrow \mathbb{Z}^d$ is injective. The conclusion of the lemma follows as in previous couplings.
\end{proof}

With Lemma \ref{lem:ksite} we are now able to conclude the goal of this section.

\begin{proof}[Proof of Theorem \ref{theo:siteNO}]
The first item of Theorem \ref{theo:siteNO} follows by taking $k=2$ in Lemma \ref{lem:ksite} toghether with Wierman's upper bound $p_c^s(2)\leq 0.68$ (see \cite{W2}). Taking $k=3$ in Lemma \ref{lem:ksite} and using the upper bound $p_c^s(3) < 1/2$ of Campanino-Russo (see \cite{CR}), the second item follows. The third item follows directly by Proposition \ref{prop:crossoversite}.
\end{proof}



\subsection{Proof of Theorem \ref{theo:siteO}}
\begin{proof}
First, we remark that Lemma \ref{lem:ksite} and Proposition \ref{prop:crossoversite} hold {\it mutatis mutandis} for oriented site Bernoulli percolation. Therefore to conclude the proof of Theorem \ref{theo:siteO}, the last ingredient is  Liggett's upper bound $\vv{p}^{\,s}_c(2) \leq 3/4$ (see \cite{L}).
\end{proof}





\section{Explicit bounds}\label{sec:comments}

In this section, we present a table with upper bounds for bond and site  Bernoulli Percolation in $\mathbb{Z}^ d$, up to $d=9$, in both the oriented and the non-oriented cases. 

In each column we give a numerical upper bound rounded to four decimals for the critical probability in the head of the column.

\begin{center}

\begin{tabular}{|c|c|c|c|c|}  \hline
\textbf{Dimension} & $p_c^b(d)$ & $\vv{p_c}^b(d)$ & $p_c^s(d)$ &   $\vv{p_c}^s(d)$\\ \hline

3 & 0,3473 &0,5680 &0,5000 &0,6422 \\   \hline

4 & 0,2788  &0,4227 &0,4344 &0,5000 \\   \hline

5 & 0,2284 & 0,3926 &0,4156& 0,4615 \\   \hline

6 & 0,1922 & 0,2734 & 0,2929 & 0,3701 \\   \hline

7 & 0,1682  & 0,2028 & 0,2866 & 0,3533\\   \hline

8 & 0,1486 &0,1627  & 0,2479  &0,2929\\   \hline

9 & 0,1326  &0,1371 &0,2063 &0,2844\\   \hline

\end{tabular}

\end{center}

\vspace{0.5cm}

To obtain each bound, we used the following:

\vspace{0.5cm}

\textbf{Non-oriented bond percolation.} All bounds were obtained from Theorem \ref{theo:bounds}. 

\textbf{Oriented bond percolation} All bounds came from Theorem \ref{theo:bondO}. For  $d=4$ we used Item 1), for  $d=5$ we used the upper bound for $d= 4$, and  Item 3). For  $d\geq 6$ we used Item 2). 


\textbf{Non-oriented site percolation.} The bound for  $d=3$ follows from Campanino-Russo (see \cite{CR}). The others  follow from Theorem \ref{theo:siteNO}. For  $d=4$ and $d=8$ we used Item 1). For   $d=6$ and $d=9$ we used  Item 2). For  $d=5$ and $d=7$ we used Item 3) along with the upper bounds obtained for $d=4$ and $d=6$.


\textbf{Oriented site percolation.}   All the bounds came from Theorem \ref{theo:siteO}.
For even dimensions the bounds follow from  Item 1), and for odd dimensions the bounds follow from Item 2) using the upper bounds for even dimensions. 


    


\begin{comment}\begin{tabular}{|c|c|c|c|c|}  \hline
\textbf{Dimension} & $p_c^b(d)$ & $\vv{p_c}^b(d)$ & $p_c^s(d)$ &   $\vv{p_c}^s(d)$\\ \hline

3 & 0,3473 &0,5680 &0,5000 &0,6422 \\   \hline

4 & 0,2787  &0,4227 &0,4344 &0,5528\\   \hline

5 & 0,2283 & 0,3926 &0,4156& 0,5086 \\   \hline

6 & 0,1921 & 0,2733 & 0,2929 & 0,3701 \\   \hline

7 & 0,1681  & 0,2028 & 0,2796 & 0,3533\\   \hline

8 & 0,1485 &0,1627  & 0,2479  &0,2929\\   \hline

9 & 0,1326  &0,1370 &0,2063 &0,2843\\   \hline

\end{tabular}
\end{comment}




\vspace{1.2cm}
{\bf \Large Acknowledgements} 

The authors thank Roger Silva for valuable comments on the manuscript. P.A. Gomes has been supported by São Paulo Research Foundation (FAPESP), grant 2020/02636-3 and grant 2017/10555-0.
R. Sanchis has been partially supported by Conselho Nacional de Desenvolvimento
Científico e Tecnológico (CNPq), CAPES and by FAPEMIG (PPM 00600/16).









\thebibliography{}


\bibitem{BR} Bollob\'as, B., and Riordan, O. {\it Percolation}, Cambridge University Press (2006).

\bibitem{CR} Campanino, M., and Russo, L. {\it An upper bound on the critical percolation probability for the three-dimensional cubic lattice}, Ann. Prob. 13 (2), 478-491, (1985).

\bibitem{CD} Cox, J.T., and  Durrett, R. {\it Oriented percolation in dimensions $d \geq 4$: bounds and asymptotic formulas.} Math Proc. Camb. Phil. Soc. Vol. 93(1), pp.151-162, (1983).

\bibitem{GPS} Gomes, P.A., Pereira, A., and Sanchis, R. {\it Anisotropic oriented percolation in high dimensions}. ALEA Lat. Am. J. Probab. Math. Stat. 17 (1), 531-543, (2020).

\bibitem{GPS2}  Gomes, P.A., Pereira, A., and Sanchis, R. {\it Anisotropic percolation in high dimensions: the non-oriented case}. ArXiv:2106.09083.





\bibitem{GSS} Gomes P.A., Sanchis, R. and Silva, R.W.C. {\it A note on the dimensional crossover critical exponent.} Lett. Math. Phys. 110 (12), 3427–3434, (2020).

\bibitem{Grim} Grimmett, G.R. {\it Percolation.} Second edition. Springer-Verlag, Berlin, (1999).



\bibitem{K3} Kesten, H. {\it The Critical Probability of Bond Percolation
on the Square Lattice Equals 1/2}, Commun. Math. Phys. 74 (1), 41-59, (1980). 

\bibitem{L} Liggett, T.M. {\it Survival of discrete time growth models with applications to oriented percolation.} Ann. Appl. Probab. 5 (3), 613-636, (1995).

\bibitem{NZ} Newman, R.E.J., and Ziff, R.M. {\it Fast Monte Carlo algorithm for site or bond percolation}, Phys. Rev. E, 64 (1), 016706, (2001).

\bibitem{W} Wierman, J.C. {\it Bond percolation critical probability bounds for the Kagomé lattice by a substitution method}. In: Disorder in Physical Systems, Grimmett, G. and Welsh, D. J. A. (eds.),
pp.349-360, (1990).

\bibitem{W2} Wierman, J.C. {\it Substitution method critical probability bounds for the square lattice site percolation model}.
Combin. Probab. Comput. 4 (2), 181-188 (1995).

\bibitem{WP} Wierman, J.C.,  and Parviainen, R.  {\it Ordering Bond Percolation Critical Probabilities}, UUDM Report, (2002)
\url{http://www.ams.jhu.edu/~wierman/Papers/Bond-Pc-ordering.pdf}


\vspace{2cm}
\begin{minipage}{0.45\textwidth}
Pablo Almeida Gomes\\
Universidade de S\~ao Paulo, Brasil\\
E-mail: \url{pagomes@usp.br}
\end{minipage}
\vspace{0.5cm}

\begin{minipage}{0.45\textwidth}
Alan Pereira\\
Universidade Federal de Alagoas, Brasil\\
E-mail: \url{alan.pereira@im.ufal.br}
\end{minipage}

\vspace{0.5cm}
\begin{minipage}{0.5\textwidth}
R\'emy Sanchis\\
Universidade Federal de Minas Gerais, Brasil\\
E-mail: \url{rsanchis@mat.ufmg.br}
\end{minipage}


\end{document}














2021-06-22T02:04:44.000
2106
2106.10388
en
https://arxiv.org/abs/2106.10388
We consider bond and site Bernoulli Percolation in both the oriented and the non-oriented cases on $\mathbb{Z}^d$ and obtain rigorous upper bounds for the critical points in those models for every dimension $d \geq 3$.
Probability (math.PR)
Upper bounds for critical probabilities in Bernoulli Percolation models
Qwen/Qwen-72B
1. YES
2. YES