Again, apologies for the initial incorrect answer. I hope I got it right this time. :-)
Here's a way to determine rigorously whether a given linear combination of spherical harmonics is chiral. The functions with different values of $l$ don't get mixed under rotations. Thus, a combination is achiral if and only if there's a single rotation that rotates the mirror images of all $l$ components into their originals. So you can look at the $l$ values individually and see whether there's any intersection among the sets of rotations.
For your example, $x+y+z+xy$, the $l=1$ component $x+y+z$ is a rotated version of $Y_{10}$ oriented towards $(1,1,1)$. The rotations that rotate the mirror image of this into the original are rotations through $\pi$ about an axis perpendicular to $(1,1,1)$. (Subsequently rotating about $(1,1,1)$ doesn't add anything; it just corresponds to a different choice of axis for the first rotation.)
The $l=2$ component $xy$ is invariant under inversion, so we need to find the rotations that leave it invariant. Its angular dependence is $\sin^2\theta\sin\phi\cos\phi$, and the only rotations that leave the $\theta$ dependence invariant are ones with axes in the $x$-$y$ plane or orthogonal to it (since the poles have to end up at the poles). There's one axis in the $x$-$y$ plane that's orthogonal to $(1,1,1)$, namely $(1,-1,0)$. Coincidentally, rotating through $\pi$ about that axis corresponds to $x\to-y$, $y\to-x$, $z\to-z$, and thus works for $x+y+z+xy$. However, this is only one of two axes in the $x$-$y$ plane for which rotating through $\pi$ leaves $xy$ invariant, since rotating the rotation axis about the $z$ axis by $\alpha$ amounts to rotating the result about the $z$ axis by $2\alpha$, and that only leaves $xy$ invariant if $2\alpha$ is a multiple of $\pi$. So you can pick almost any other combination of the $l=1$ functions, as long as the associated vector isn't orthogonal to either the $z$ axis or the axes $(1,-1,0)$ and $(1,1,0)$, and add $xy$ to get a chiral combination.