This research work is a part of the project No. 2022/47/P/ST5/01816 within the POLONEZ BIS programme co-funded by the National Science Centre and the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No. 945339. For the purpose of Open Access, the author has applied a CC-BY public copyright licence to any Author Accepted Manuscript (AAM) version arising from this submission. The costs associated with open access publication were supported by the NAWA STARS EU BOOSTER project, funded through the FERS.01.05-IP.08-0219/23 program. Calculations have been partially carried out at the Wroclaw Centre for Networking and Supercomputing http://www.wcss.wroc.pl (Grant no. 171). We would like to express our gratitude to Dr Yurii Kazarinov for help in PL spectra interpretation; Dr Wojciech Piotrowski ( https://orcid.org/0000- 0002-8155-2688 ) for help in PL spectra in NIR region measurement.