Second quantization and bosons in Tomonaga's "The story of spin"

TL; DR: According to a book by Tomonaga in a chapter introducing second quantization, the formula
$$e^{pm i Theta/hbar} ,psi(N) = psi(Npm 1)$$
is supposed to prove that the second quantization procedure using commutation relations is suitable for bosons. I am able to follow Tomonaga’s mathematics up to this formula, but I do not understand how he can conclude that it proves the correctness of his theory for bosons. Could you please explain it to me or point me towards an easy/intuitive explanation? Unfortunately I am not an expert in these topics…

---

Full question:
I am self-studying "The story of spin" by Sin-itiro Tomonaga and I am struggling with lecture 6, where he introduces the second quantization formalism. In this lecture, the author starts from the single particle Schrodinger equation and recast it in terms of canonical equation of motions for Hamiltonian mechanics. Then, he moves to the generalization of these ideas for many particles.

Considering a single particle observable $G$ with eigenvalues $g_n$ each one with a probability $P_n = |a_n|^2$, if $N$ is the number of particles considered then we can write the expectation value $N_n$ for the number of particles of the ensemble with a value $g_n$ for $G$ as
$$N_n = NP_n = N |a_n|^2 = A_n A_n^*$$
from which
$$ A_n = sqrt{N} a_n,quad A_n^* = sqrt{N} a_n^*. $$
From this he then defines the conjugated momentum as
$$ Pi_n = i hbar A_n^* $$
from which follows the expectation value of the hamiltonian of the system
$$ bar{H} = -frac{i}{hbar} sum_{n,m} Pi_n H_{n,m} A_m $$
where $H_{n,m}$ is the matrix element of the single particle hamiltonian. We are now able to derive from these equations the canonical equations for the ensemble
$$ frac{dA_n}{dt} = frac{partialbar{H}}{partial Pi_n},quadfrac{dPi_n}{dt} = -frac{partialbar{H}}{partial A_n} $$
where the following normalization condition holds
$$ sum_n |A_n|^2 = N. $$
At this point, Tomonaga recognizes the important contribution from Dirac to the second quantization formalism, since he redefined $A_n$ and $Pi_n$ as quantum numbers and he introduced the following commutation relations
$$ [A_n,Pi_n] = ihbar delta_{n,m},quad [A_n, A_m] = [Pi_n, Pi_m] = 0. $$
He also states that if we are worried by the fact that $A_n$ and $Pi_n$ are complex numbers, we can use the following definitions
$$ A_n = sqrt{N_n} e^{iTheta_n /hbar},quad A_n^* = sqrt{N_n} e^{-iTheta_n/hbar} $$
from which we obtain the following hamiltonian
$$bar{H} = sum_{n,m} sqrt{N_n} e^{-iTheta_n /hbar} H_{n,m} sqrt{N_m} e^{iTheta_m/hbar}. $$
After this mathematics, Tomonaga introduces the concept of "virtual ensembles" and says that in the second quantization formalism we are assuming that we are able to describe such a system of $N$ (non-interacting) particles in terms of the quantum numbers $N_n$ and $Theta_n$ or $A_n$ and $Pi_n$. To prove that this is correct, we need to demonstrate that the "usual" description in terms of a wavefunction in the coordinate space satisfying the Schrodinger equation
$$ left[ H(mathbf{x}_1, mathbf{p}_1) + H(mathbf{x}_2, mathbf{p}_2) + dots + H(mathbf{x}_N, mathbf{p}_N) -ihbar frac{partial}{partial t} right] psi(mathbf{x}_1, mathbf{x}_2,dots,mathbf{x}_N) = 0 $$
is equivalent to the description in the second quantization formalism
$$ left[sum_{n,m} sqrt{N_n} e^{-iTheta_n/hbar} H_{n,m} e^{iTheta_m/hbar } sqrt{N_m} -ihbar frac{partial}{partial t}right] psi(N_1, N_2, dots, N_n,dots) = 0. $$
To do this, Tomonaga uses an argument by Dirac. He starts considering that $Theta$ as we have defined it is the momentum conjugated to $N$, so it must hold that
$$ Theta = -ihbarfrac{partial}{partial N} .$$
If this holds, we can make the following series expansion
$$ e^{pm iTheta/hbar} = e^{pm partial/partial N} = 1pmfrac{partial}{partial N} pm frac{1}{2!}frac{partial^2}{partial N^2} + dots $$
and applying its right-hand side to the wavefunction and recognizing Taylor’s theorem on the series expansion of a function
$$ psi(N) pm psi'(N) + frac{1}{2!} psi”(N) + dots = psi(Npm 1) $$
we obtain the final equation
begin{equation} e^{pm i Theta/hbar} ,psi(N) = psi(Npm 1).quadquad(star)end{equation}
At this point, Tomonaga’s demonstration stops. Quoting his own words:

Since the newly discovered theory has now been proved to be correct for bosons, the use of heuristic theory notwithstanding, there is no reason for us to hesitate to use it. We can proceed with it with confidence.
and he starts applying this new formalism.

My problem is that to me is not clear at all why we can conclude from Eq. $(star)$ that this theory is suitable for bosons and that this approach makes the formulation in the coordinate space using the Schrodinger equation equivalent to the second quantization approach. I am (more or less) able to follow the mathematics, what I miss is probably the physical meaning of Eq. $(star)$. Can you please help me understanding how is it possible to draw such a conclusion? What is the meaning of Eq. $(star)$?

It also seems that the very same problem was first discussed in Dirac’s paper The Quantum Theory of the Emission and Absorption of Radiation in Eqs. 10 to 13 at pages 251 and 252. In the case of this new demonstration, I am stuck in the demonstration of Eq. 10 starting from the definition in Eq. 8 and the commutation relations reported in the unnumbered equations preceding Eq. 10. Do you have any hint on this?