Again, I’m not really familiar with modular theory, so I can’t really do this fully rigorously. However, since this question still has no answer at this point and this is a bit too long for a comment, I’m writing an answer here to present my understanding of how this probably should work. I advise the OP to take my answer with a grain of salt, and mostly treat this as a proof idea to be filled in with details by themself.
For $x, y$ satisfying suitable assumptions - presumably $x \in \mathfrak{n}_\psi$, $y \in \mathfrak{n}_\psi^\ast$, and both elements are entire (i.e., $\sigma_\alpha(x)$ and $\sigma_\alpha(y)$ are well-defined for all $\alpha \in \mathbb{C}$), we have,
$$\begin{split}
\eta_\psi(x)y &= Jy^\ast J\eta_\psi(x)\\
&= Jy^\ast S\Delta^{-\frac{1}{2}}\eta_\psi(x)\\
&= Jy^\ast S\eta_\psi(\sigma_{\frac{1}{2}i}(x))\\
&= Jy^\ast \eta_\psi(\sigma_{\frac{1}{2}i}(x)^\ast)\\
&= S\Delta^{-\frac{1}{2}}\eta_\psi(y^\ast \sigma_{-\frac{1}{2}i}(x^\ast))\\
&= S\eta_\psi(\sigma_{\frac{1}{2}i}(y^\ast \sigma_{-\frac{1}{2}i}(x^\ast)))\\
&= S\eta_\psi(\sigma_{\frac{1}{2}i}(y^\ast)x^\ast)\\
&= \eta_\psi(x\sigma_{\frac{1}{2}i}(y^\ast)^\ast)\\
&= xS\eta_\psi(\sigma_{\frac{1}{2}i}(y^\ast))\\
&= xJ\Delta^{\frac{1}{2}}\eta_\psi(\sigma_{\frac{1}{2}i}(y^\ast))\\
&= xJ\eta_\psi(y^\ast)\\
&= x\eta_\psi’(y)
\end{split}$$
Using some approximation argument (presumably weak$^\ast$ density of entire elements and the Kaplansky density theorem), one should be able to deduce then that $\eta_\psi(x)y = x\eta_\psi’(y)$ for all $x \in \mathfrak{n}_\psi$, $y \in \mathfrak{n}_\psi^\ast$. The desired result immediately follows.
(Heuristically, $\eta_\psi(x)y = x\eta_\psi’(y)$ is just saying that $(x\psi^{\frac{1}{2}})y = x(\psi^{\frac{1}{2}}y)$.)