A robust design method for sizing renewable energy systems in nZEB is studied. Monte Carlo simulation is applied for uncertainty analysis. The relationship between mismatch ratio and probability to be nZEB is identified. Case studies are conducted on the Hong Kong Zero Carbon Building.