  A procedure to find optimal regimes for quantum thermal engines (QTMs) is
described and demonstrated. The QTMs are modelled as the periodically-driven
non-equilibrium steady states of open quantum systems, whose dynamics is
approximated in this work with Markovian master equations. The action of the
external agent, and the couplings to the heat reservoirs can be modulated with
control functions, and it is the time-dependent shape of these control
functions the object of optimisation. Those functions can be freely
parameterised, which permits to constrain the solutions according to
experimental or physical requirements.
