Is there a way I can use JuMP directly to solve this with Gurobi?
Not yet. An upcoming release of Gurobi will offer better support for nonlinear (non-quadratic) functions.
I can introduce another variable and split the trilinear term into two quadratic terms
This is the way.