Bipartite entanglement dynamics of two-level systems in sub-Ohmic reservoirs

The quantum dynamics of pairs of two level systems immersed in dissipative reservoirs with sub-Ohmic spectral distributions is studied by means of numerically exact path integral Monte Carlo methods. It is shown that this class of reservoirs, relevant for generic properties of strongly non-Markovian environments with possible realizations in solid state structures at cryogenic temperatures, supports bipartite entanglement in broad ranges of parameter space and even for broken symmetries and in presence of thermal noise. The sensitivity of creation of entanglement in non-equilibrium on initial preparations of the reservoir is investigated in detail.