[step:Pass to the limit $\varepsilon \to 0$ to obtain plurisubharmonicity of $\phi^M$]Fix $M > 0$. We claim $\phi^M_\varepsilon \searrow \phi^M$ pointwise on $\Omega$ as $\varepsilon \searrow 0$ (where for each fixed $z \in \Omega$ the statement is interpreted for $\varepsilon < \operatorname{dist}(z, \partial \Omega)$).
For a radial mollifier $\eta_\varepsilon$, write $\eta(w) = h(|w|)$ for some $h \in C_c^\infty([0,1))$ with $h \ge 0$. Switching to polar coordinates on $\mathbb{R}^{2n}$ with $w = r u$, $r \ge 0$, $u \in S^{2n-1}$, the polar formula $d\mathcal{L}^{2n}(w) = r^{2n-1}\, d\mathcal{H}^{2n-1}(u)\, d\mathcal{L}^1(r)$ gives
\begin{align*}
\phi^M_\varepsilon(z) = \int_0^\varepsilon \varepsilon^{-2n} h(r/\varepsilon)\, r^{2n-1} \biggl( \int_{S^{2n-1}} \phi^M(z - r u)\, d\mathcal{H}^{2n-1}(u) \biggr) d\mathcal{L}^1(r).
\end{align*}
Set $\sigma_{2n-1} := \mathcal{H}^{2n-1}(S^{2n-1})$ and define the spherical mean
\begin{align*}
\mathcal{M}(\phi^M, \cdot, \cdot): \{(z, r) : z \in \Omega,\ 0 < r < \operatorname{dist}(z, \partial\Omega)\} &\to \mathbb{R} \\
(z, r) &\mapsto \frac{1}{\sigma_{2n-1}} \int_{S^{2n-1}} \phi^M(z + r u)\, d\mathcal{H}^{2n-1}(u).
\end{align*}
By radial symmetry of $\eta$ (substituting $u \mapsto -u$),
\begin{align*}
\phi^M_\varepsilon(z) = \int_0^\varepsilon \varepsilon^{-2n} h(r/\varepsilon)\, r^{2n-1}\, \sigma_{2n-1}\, \mathcal{M}(\phi^M, z, r)\, d\mathcal{L}^1(r).
\end{align*}
The normalisation $\int_{\mathbb{C}^n} \eta_\varepsilon\, d\mathcal{L}^{2n} = 1$ rewrites as $\int_0^\varepsilon \varepsilon^{-2n} h(r/\varepsilon)\, r^{2n-1}\, \sigma_{2n-1}\, d\mathcal{L}^1(r) = 1$, so $\phi^M_\varepsilon(z)$ is a weighted average of $\mathcal{M}(\phi^M, z, \cdot)$ over $r \in (0, \varepsilon)$.
[claim:Monotonicity of spherical means]
For each $z \in \Omega$, the map $r \mapsto \mathcal{M}(\phi^M, z, r)$ is non-decreasing on $(0, \operatorname{dist}(z, \partial \Omega))$.
[/claim]
[proof]
Identify $S^{2n-1} \subset \mathbb{C}^n$ as the disjoint union of circles $\{e^{i\theta} u : \theta \in [0, 2\pi)\}$ as $u$ ranges over the base of the Hopf fibration $S^{2n-1} \to \mathbb{CP}^{n-1}$ given by the $S^1$-action $\theta \cdot u = e^{i\theta} u$. Let $\nu$ denote the unique $S^1$-invariant Borel probability measure on $\mathbb{CP}^{n-1}$. Disintegration of $\mathcal{H}^{2n-1}$ along this fibration (each fibre carries arc-length measure on a circle of radius $1$, total mass $2\pi$) gives the equality
\begin{align*}
\mathcal{M}(\phi^M, z, r) = \int_{\mathbb{CP}^{n-1}} \biggl( \frac{1}{2\pi} \int_0^{2\pi} \phi^M(z + r e^{i\theta} u)\, d\mathcal{L}^1(\theta) \biggr) d\nu(u),
\end{align*}
where the normalising constants on the two sides match because $\nu$ is a probability measure and the inner integral is normalised by $1/(2\pi)$.
For each fixed $u \in S^{2n-1}$, the inner integral
\begin{align*}
A_u(r) := \frac{1}{2\pi} \int_0^{2\pi} \phi^M(z + r e^{i\theta} u)\, d\mathcal{L}^1(\theta)
\end{align*}
is the circular mean at radius $r$ of the slice $\zeta \mapsto \phi^M(z + \zeta u)$ along the complex line through $z$ in direction $u$. By Step 3 (slice subharmonicity of $\phi^M$), this slice is subharmonic on the open disc $\{\zeta \in \mathbb{C} : |\zeta| < \operatorname{dist}(z, \partial \Omega)\}$ (which is contained in $S_{z, u}$); circular means of a subharmonic function are non-decreasing in the radius. Thus $r \mapsto A_u(r)$ is non-decreasing for $r < \operatorname{dist}(z, \partial \Omega)$. Integrating over $u$ against the non-negative probability measure $\nu$ preserves the monotonicity, so $r \mapsto \mathcal{M}(\phi^M, z, r)$ is non-decreasing.
[/proof]
For fixed $z$, the weighting in $\phi^M_\varepsilon(z)$ concentrates near $r = 0$ as $\varepsilon \searrow 0$, and the integrand $\mathcal{M}(\phi^M, z, r)$ is non-decreasing in $r$; hence $\varepsilon \mapsto \phi^M_\varepsilon(z)$ is non-decreasing in $\varepsilon$ (equivalently, non-increasing as $\varepsilon \searrow 0$). Indeed, for $0 < \varepsilon_1 < \varepsilon_2 < \operatorname{dist}(z, \partial \Omega)$, applying the substitution $r = \varepsilon_i s$ (so $d\mathcal{L}^1(r) = \varepsilon_i\, d\mathcal{L}^1(s)$, and the domain $r \in (0, \varepsilon_i)$ maps to $s \in (0, 1)$) in each integral gives:
\begin{align*}
\phi^M_{\varepsilon_i}(z) = \int_0^1 h(s)\, s^{2n-1}\, \sigma_{2n-1}\, \mathcal{M}(\phi^M, z, \varepsilon_i s)\, d\mathcal{L}^1(s),
\end{align*}
and since $\mathcal{M}(\phi^M, z, \varepsilon_1 s) \le \mathcal{M}(\phi^M, z, \varepsilon_2 s)$ pointwise by monotonicity, $\phi^M_{\varepsilon_1}(z) \le \phi^M_{\varepsilon_2}(z)$.
To identify the limit, upper semicontinuity of $\phi^M$ gives $\limsup_{r \to 0^+} \mathcal{M}(\phi^M, z, r) \le \phi^M(z)$. We apply the reverse [Fatou lemma](/theorems/510) to the family $u \mapsto \phi^M(z + r u)$ on $(S^{2n-1}, \mathcal{H}^{2n-1})$ as $r \to 0^+$. The hypothesis of reverse Fatou is a uniform integrable upper bound: by local boundedness of $\phi^M$, there exist $\rho > 0$ and $C \in \mathbb{R}$ such that $\phi^M(y) \le C$ for every $y \in \overline{B}(z, \rho)$; the constant function $C$ is $\mathcal{H}^{2n-1}$-integrable on $S^{2n-1}$ (which has finite measure $\sigma_{2n-1}$); and for $r < \rho$ we have $\phi^M(z + r u) \le C$ for every $u \in S^{2n-1}$. Reverse Fatou then gives
\begin{align*}
\limsup_{r \to 0^+} \frac{1}{\sigma_{2n-1}} \int_{S^{2n-1}} \phi^M(z + r u)\, d\mathcal{H}^{2n-1}(u)
\le \frac{1}{\sigma_{2n-1}} \int_{S^{2n-1}} \limsup_{r \to 0^+} \phi^M(z + r u)\, d\mathcal{H}^{2n-1}(u)
\le \phi^M(z),
\end{align*}
where the last inequality uses $\limsup_{r \to 0^+} \phi^M(z + r u) \le \phi^M(z)$ for every $u$ by upper semicontinuity of $\phi^M$ at $z$. Combined with the sub-mean-value inequality $\phi^M(z) \le \mathcal{M}(\phi^M, z, r)$ (which holds for subharmonic functions of one complex variable applied to slices, then averaged as above), we get $\lim_{r \to 0^+} \mathcal{M}(\phi^M, z, r) = \phi^M(z)$. Hence $\phi^M_\varepsilon(z) \to \phi^M(z)$ as $\varepsilon \to 0^+$. Combined with monotonicity, $\phi^M_\varepsilon \searrow \phi^M$.
By [Stability Properties of PSH Functions](/theorems/3404), the decreasing limit of plurisubharmonic functions on a fixed [open set](/page/Open%20Set) is either plurisubharmonic on that set or identically $-\infty$ on each connected component of that set. We apply this fact one $\varepsilon$-shell at a time. Fix $\varepsilon_0 > 0$. For every $\varepsilon \in (0, \varepsilon_0]$, $\Omega_{\varepsilon_0} \subseteq \Omega_\varepsilon$, and $\phi^M_\varepsilon$ is plurisubharmonic on $\Omega_\varepsilon$ by Step 4, hence its restriction is plurisubharmonic on $\Omega_{\varepsilon_0}$. The sequence $(\phi^M_\varepsilon|_{\Omega_{\varepsilon_0}})_{\varepsilon \le \varepsilon_0}$ is decreasing with pointwise limit $\phi^M|_{\Omega_{\varepsilon_0}}$. Since $\phi^M \ge -M > -\infty$, the identically $-\infty$ alternative is excluded on every component, so $\phi^M$ is plurisubharmonic on $\Omega_{\varepsilon_0}$. As $\varepsilon_0 > 0$ was arbitrary and $\Omega = \bigcup_{\varepsilon_0 > 0} \Omega_{\varepsilon_0}$, plurisubharmonicity at every point of $\Omega$ follows from the local nature of plurisubharmonicity (each $z \in \Omega$ has a neighbourhood contained in some $\Omega_{\varepsilon_0}$). Thus $\phi^M$ is plurisubharmonic on $\Omega$.[/step]