We discuss the use of Wilson fermions with twisted mass for simulations of QCD thermodynamics. As a prerequisite for a future analysis of the finite-temperature transition making use of automatic O(a) improvement, we investigate the phase structure in the space spanned by the hopping parameter κ, the coupling β , and the twisted mass parameter µ. We present results for N f = 2 degenerate quarks on a 16 3 × 8 lattice, for which we investigate the possibility of an Aoki phase existing at strong coupling and vanishing µ, as well as of a thermal phase transition at moderate gauge couplings and non-vanishing µ.