Quoting from Green-Tao, "Linear equations in primes" (especially Cor. 1.9 in https://arxiv.org/pdf/math/0606088.pdf), any system of linear forms of finite complexity and without any local obstructions will assume simultaneous prime values infinitely often.

(E.g., $(X,Y,X+Y)$ is of finite complexity, but with a local obstruction at 2, since one of these is always even. On the other hand $(X,Y,X+2Y)$ has no local obstruction, so there are infinitely many primes $x,y$ such that $x+2y$ is also prime.)

My question is: Does this result remain true when, instead of asking for prime values, one asks for values in some positive density *subset* of the primes? This kind of strengthening is true e.g. for the famous special case $X, X+Y, X+2Y,...$ of arithmetic progressions in primes, but I couldn't find it in the general case.

Again, to illustrate what I mean: Let's take the set $S$ of all primes $\equiv 1$ mod $4$. Then the system $(X,Y,3X+4Y)$ will be obstructed, since for $x,y\equiv 1$ mod $4$, one has $3x+4y\equiv 3$ mod $4$. On the other hand, $(X,Y,5X+4Y)$ would not be obstructed, so should be expected to have infinitely many values in $S^3$. Is this known in generality?

[I should add that for the last example, one may of course just replace X and Y by 4X+1 and 4Y+1; but of course there are positive density subsets of primes which are not just residue classes. I'm interested in those as well.]

[EDIT: To be more explicit: Let S be a Chebotarev set of primes, i.e. the set of all primes with some given Frobenius in some given Galois extension $K/\mathbb{Q}$. Then there will usually be a congruence condition (but not only!) on $S$ coming from the largest cyclotomic subfield of $K$. So there is some $N$ and some smallest possible finite union $U$ of mod $N$ residue classes containing all of $S$. A local ostruction for a system of linear forms now definitely occurs if there are no integer values simultaneously lying in $U$. I'm asking if that is the only obstruction for values to simultaneously lie in $S$.

Example: if S is the set of primes which decompose completely in the splitting field of $x^3-x^2+1$, then the largest abelian subfield is $\mathbb{Q}(\sqrt{-23})$. To be decomposed in this field means to be a square mod 23 (leaving 11 of 23 possible residue classes). I'm tempted to conjecture that a system of linear forms of finite complexity which has integer values simultaneously lying in one of these 11 classes (i.e. not obstructed at 23), and which for any other prime $p$ has integer solutions simultaneously coprime to $p$ (i.e. not obstructed at $p$) will have prime values simultaneously lying in $S$.]