Eccentric-orbit extreme-mass-ratio-inspiral radiation II: 1PN correction to leading-logarithm and subleading-logarithm flux sequences and the entire perturbative 4PN flux

In a recent paper we showed that for eccentric-orbit extreme-mass-ratio inspirals the analytic forms of the leading-logarithm energy and angular momentum post-Newtonian (PN) flux terms (radiated to infinity) can, to arbitrary PN order, be determined by sums over the Fourier spectrum of the Newtonian quadrupole moment. We further showed that an essential part of the eccentricity dependence of the related subleading-logarithm PN sequences, at lowest order in the symmetric mass ratio $\nu$, stems as well from the Newtonian quadrupole moment. Once that part is factored out, the remaining eccentricity dependence is more easily determined by black hole perturbation theory. In this paper we show how the sequences that are the 1PN corrections to the entire leading-logarithm series, namely terms that appear at PN orders $x^{3k+1} \log^k(x)$ and $x^{3k+5/2} \log^k(x)$ (for PN compactness parameter $x$ and integers $k\ge 0$), at lowest order in $\nu$, are determined by the Fourier spectra of the Newtonian mass octupole, Newtonian current quadrupole, and 1PN part of the mass quadrupole moments. We also develop a conjectured (but plausible) form for 1PN correction to the leading logs at second order in $\nu$. Further, in analogy to the first paper, we show that these same source multipole moments also yield nontrivial parts of the 1PN correction to the subleading-logarithm series, and that the remaining eccentricity dependence (at lowest order in $\nu$) can then more easily be determined using black hole perturbation theory. We use this method to determine the entire analytic eccentricity dependence of the perturbative (i.e., lowest order in $\nu$) 4PN non-log terms, $\mathcal{R}_4(e_t)$ and $\mathcal{Z}_4(e_t)$, for energy and angular momentum respectively.