Calculating work in adiabatic two-level quantum Markovian master equations: A characteristic function method

We present a characteristic function method to calculate the probability density functions of the inclusive work in adiabatic two-level quantum Markovian master equations. These systems are steered by some slowly varying parameters and the dissipations may depend on time. Our theory is based on the interpretation of the quantum jump for the master equations. In addition to the calculation, we also find that the fluctuation properties of the work can be described by the symmetry of the characteristic functions, which is exactly the same as in the case of isolated systems. A periodically driven two-level model is used to demonstrate the method.