I didn't test my expand function extensively but it seems it works and can help you. You can also take a look at various links in the topic it was created in:
MMA.SE: Isolating cross terms
expand[
y1[z],
"SmallTerms" -> {Subscript[\[Xi], 0], Subscript[\[Rho], 0]}
] // Normal // Expand
