Define
\begin{equation}
h(z) := \prod_{n=1}^\infty (1 - z/\rho_n) e^{z/\rho_n};\tag{9.2.3}
\end{equation}
the product converges absolutely uniformly on any closed disc because the multiplicand is
\begin{equation*}
1 + \frac{1}{2} \left(\frac{z}{\rho_n} \right)^2 + O\left(\left(\frac{z}{\rho_n}
\right)^3\right)
\end{equation*}
and the fact that the number of roots of norm
\(\leq R\) is
\(O(R^{1+\epsilon})\) (by
Remark 9.6) implies that
\(\sum 1/|\rho_n|^2\) converges (by partial summation).
We now argue (following
[3], section 11) that
\(f/h\) is also of order
\(\leq 1\text{.}\) It will suffice to establish an upper bound on
\(|1/h(z)|\) of the form
\(\exp(R^{1+\epsilon})\) when
\(|z| = R\) for some infinite divergent sequence of values of
\(R\text{,}\) as then this will imply a corresponding bound for the entire function
\(f/h\) and we can use the maximum modulus principle to fill in other values of
\(R\text{.}\) (Note that we cannot take this last step directly for
\(1/h\) because it is not entire.) Put
\(r_n := |\rho_n|\text{.}\) Since
\(\sum 1/r_n^2\) converges, the total lengths of the intervals
\([r_n - r_n^{-2}, r_n + r_n^{-2}]\) is finite, so we can choose an infinite divergent sequence of values of
\(R\) for each of which we have
\begin{equation*}
|R-r_n| > r_n^{-2} \qquad (n=1,2,\dots).
\end{equation*}
To prove the claimed bound for a particular such
\(R\text{,}\) we factor
\(h = h_1 h_2 h_3\) by partitioning the product
(9.2.3) according to which range
\(r_n\) lies in:
\begin{equation*}
h_1\colon [0, R/2), \qquad h_1\colon [R/2, 2R], \qquad h_3: (2R, \infty).
\end{equation*}
It will now suffice to prove that
\begin{equation}
|z| = R \Longrightarrow |h_i(z)| > \exp(-R^{1+\epsilon}) \qquad (i=1,2,3);\tag{9.2.4}
\end{equation}
At this point, we know that \(f/h\) is entire of order \(\leq 1\) with no zeroes at all. Consequently, \(g(z) = \log (f(z)/h(z))\) is entire and satisfies \(|g(z)| = O(|z|^{1+\epsilon})\text{.}\) This in turn means that
\begin{equation*}
g_2(z) = \frac{g(z) - g(0) - g'(0)z}{z^2}
\end{equation*}
is entire and bounded, hence constant by Liouville's theorem. This yields the desired result.