How to create interaction variable between only 1 variable with all other variables
You could just make them yourself no? PolynomialFeatures doesn't do anything particularly innovative.
When using the SciKit learn PolynomialFeatures package it is possible to do the following:

Take features: x1, x2, x3
Create interaction variables: [x1, x2, x3, x1x2, x2x3, x1x3] using PolynomialFeatures(interaction_only=True)

My problem is that I only want the interaction terms between x1 and all other terms meaning:
[x1, x2, x3, x1x2, x1x3], I do not want x2x3.
Is it possible to do this?