First numerical observation of the Berezinskii-Kosterlitz-Thouless transition in language models