Properties of numerous selenium derivatives of maltol (3-hydroxy-2-methyl-4H-pyran-4-one), which is an important ligand in medicinal chemistry, were studied. Selenium derivatives were obtained using consecutive in silico replacements of maltol oxygen atoms, which produced three mono-substituted, three di-substituted and one tri-substituted derivatives. The geometry optimisations for these compounds in neutral, anionic and cationic forms were performed at the B1LYP/6-311++G** level of theory. The aromatic properties of the considered selenium compounds were studied in their heterocyclic rings using several indices of aromaticity: harmonic oscillator model of aromaticity, nucleus independent chemical shift(0), nucleus independent chemical shift(1), aromatic stabilisation energy, H index, para delocalisation index and multicentre index. The methods cover all of the main types of aromaticity indices: geometric, magnetic, energetic and electronic. The calculated data show the identical aromaticity order (cation > neutral molecule > anion) for maltol and its thio and seleno derivatives. The predicted aromaticity levels of selenomaltols were compared with the previously studied maltol and various thiomaltols.