Wiki defines an almost-prime zeta function as a sum of inverse powers of the k-primes (the integers which are a products of $k$ not necessarily distinct primes): $$P_k(s)=\sum_{n: \Omega(n)=k} \frac{1}{n^s}$$ There is also an obvious identity that decomposes the Riemann zeta function into an infinite sum of the $P_{k}$:
$$\zeta(s)=1+\sum_{k=1}^{\infty} P_k(s)$$ My question is about partial (finite) sums of the $P_k(s)$: $$\zeta(s)_N=1+\sum_{k=1}^{N} P_k(s)$$ These $\zeta(s)_N$ are a sort of intermediate function beetwin prime zeta function ($N=1$) and Riemann zeta function ($N= \infty$). If there is a lot of literature on the prime zeta function, then I could not find anything about the sums of $N>1$. What can be said about the zeros and poles of $\zeta(s)_N$? How does their location change with $N$?
Here goes (my solution using symmetric polynomials and residues): First, observe that by an arithmetic argument, $$P_k(s) = [u^k] \prod_p \left(1 + \sum_{r \geq 1} p^{-rs} u^r\right) = \prod_p \left(1-\frac{u}{p^s}\right)^{-1},$$ where setting $s \mapsto 1$ gives the usual (non-split) Euler product for $\zeta(s)$. We can interpret the RHS product in terms of the complete homogeneous symmetric polynomials for $x_j := p_j^{-s}$: $$P_k(s) = h_k(x_1,x_2,\ldots).$$ Now by generating function expansions for these polynomials, we have that $$h_k = [z^k] \exp\left(-\sum_{k \geq 1} \frac{p_k}{k} z^k\right),$$ where $$p_k = \sum_p p^{-ks} = P(sk).$$ Here, $P(s)$ is the prime zeta function. So by generating function arithmetic, we have that $$\zeta(s)_N = [z^N] \frac{1}{1-z} \exp\left(-\sum_{k \geq 1} \frac{p_k}{k} z^k\right).$$ The RHS has poles at $s := 1/k$ for $k \geq 1$ by the nature of singularities of $P(s)$. From an inverse Z-transform, we can recover these coefficients as $$\begin{align} \zeta(s)_N & = \frac{1}{2\pi\imath} \oint_C \frac{s^{N-1}}{1-s} \exp\left(-\sum_{k \geq 1} \frac{P(sk)}{k} s^k\right) ds \\ & = \sum_{k \geq 1} \operatorname{Res}_{s=\frac{1}{k}}\left[\frac{s^{N-1}}{1-s} \exp\left(-\sum_{k \geq 1} \frac{P(sk)}{k} s^k\right)\right]. \end{align}$$ Since we can expand $P(s)$ near $s \rightarrow 1^{+}$ by $$P(1+\varepsilon) \sim -\log(\varepsilon) + B + O(\varepsilon),$$ my back-of-the-envelope calculations suggest that these residues correspond to (NOTE: I haven't checked these out numerically) $$R_k = \begin{cases} -e^{-B} \exp\left(-\sum\limits_{j \geq 2} \frac{P(j)}{j}\right), & k = 1 \\ \frac{1}{k^N (k-1)} \exp(-B / k^{k+1}) \exp\left(-\sum\limits_{\substack{j \geq 1 \\ j \neq k}} \frac{P(kj)}{j \cdot k^j}\right), & k \geq 2. \end{cases}$$ The constant $B$ in the previous equations is defined by $$B = \sum_{n \geq 2} \frac{\mu(n)}{n} \log \zeta(n) \approx -0.315718452.$$ Now you should be able to proceed with evaluating zeros for the function.
N.b.: With respect to there not being much in the way of literature on $P(s)$, my friend has pointed out to me offline that the prime zeta function is essentially meaningless due to the essential and messy, complicated structure of its poles. It cannot be analytically continued past zero, and hence, is not nicely modular, nor does it have a nice functional equation like most zeta functions do. It's easier just not to work with the function unless you are forced.