A comprehensive model atom was developed for Si i-ii-iii using the most up-todate atomic data available so far. Based on the non-local thermodynamic equilibrium (NLTE) line formation for Si i, Si ii, and Si iii and high-resolution observed spectra, we determined the NLTE abundances for a sample of nine unevolved A9 to B3-type stars with well determined atmospheric parameters. For each star, NLTE reduces substantially the line-to-line scatter for Si ii compared with the LTE case and leads to consistent mean abundances from lines of different ionisation stages. In the hottest star of our sample, ι Her, Si ii is subject to overionisation that drives emission in the lines arising from the high-excitation doublet levels. Our NLTE calculations reproduced 10 emission lines of Si ii observed in ι Her. The same overionisation effect leads to greatly weakened Si ii lines, which are observed in absorption in ι Her. Large positive NLTE abundance corrections (up to 0.98 dex for 5055Å) were useful for achieving consistent mean abundances from lines of the two ionisation stages, Si ii and Si iii. It was found that the NLTE effects are overestimated for the Si ii 6347, 6371Å doublet in ι Her, while the new model atom works well for the cooler stars. At this stage, we failed to understand this problem. We computed a grid of the NLTE abundance corrections for lines of Si i, Si ii, and Si iii in the model atmospheres with effective temperatures and surface gravities characteristic of unevolved A-B type stars.