To study the nuclear ( 1 kpc) dust of nearby (z < 0.1) type 1 Quasi Stellar Objects (QSOs) we obtained new near-infrared (NIR) high angular resolution (∼ 0.3 arcsec) photometry in the H and Ks bands, for 13 QSOs with available mid-infrared (MIR) high angular resolution spectroscopy (∼ 7.5 − 13.5 µm). We find that in most QSOs the NIR emission is unresolved. We subtract the contribution from the accretion disk, which decreases from NIR (∼ 35%) to MIR (∼ 2.4%). We also estimate these percentages assuming a bluer accretion disk and find that the contibution in the MIR is nearly seven time larger. We find that the majority of objects (64%, 9/13) are better fitted by the Disk+Wind H17 model (Hönig & Kishimoto 2017), while others can be fitted by the Smooth F06 (14%, 2/13, Fritz et al. 2006), Clumpy N08 (7%, 1/13, Nenkova et al. 2008a,b), Clumpy H10 (7%, 1/13, Hönig & Kishimoto 2010b), and Two-Phase media S16 (7%, 1/13, Stalevski et al. 2016) models. However, if we assume the bluer accretion disk, the models fit only 2/13 objects. We measured two NIR to MIR spectral indexes, α N IR−M IR(1.6,8.7 µm) and α N IR−M IR(2.2,8.7 µm) , and two MIR spectral indexes, α M IR(7.8,9.8 µm) and α M IR(9.8,11.7 µm) , from models and observations. From observations, we find that the NIR to MIR spectral indexes are ∼ −1.1 and the MIR spectral indexes are ∼ −0.3. Comparing the synthetic and observed values, we find that none of the models simultaneously match the measured NIR to MIR and 7.8 − 9.8 µm slopes. However, we note that measuring the α M IR(7.8,9.8 µm) on the starburst-subtracted Spitzer/IRS spectrum, gives values of the slopes (∼ −2) that are similar to the synthetic values obtained from the models.