AbstractThe unification model of active galactic nuclei postulates an accreting supermassive black hole as the central engine, surrounded by a putative dusty torus. This dust absorbs the incoming radiation, re-emits it in the infrared and obscures our view of the central region at certain inclinations. We present a new set of AGN models, in which the torus is modelled as a 3D multiphase medium. These new models can explain the observed spectral energy distribution of AGNs over the entire infrared domain, including the observed silicate feature strength and the level of near-infrared continuum. A new generation of multi-phase models, based on hydrodynamical simulations, is being constructed. We will compute the polarisation structure of these physically motivated 3D torus models, and compare them to simpler smooth torus models and to the available observational data.