Átomos alcalinos são boas sondas para compreender as propriedades do He líquido. Considerável atenção experimental tem sido empregada para analisar as mudanças da posição e da largura da linha do espectro de absorção de átomos alcalinos imersos em um ambiente de He líquido. No lado teórico, vários estudos têm usado modelos simplificados como o modelo de bolhas e o modelo de agregado. Considerações de modelos mais realista agora são oportunas e relevantes. Neste trabalho, nós usamos a combinação da simulação de Monte Carlo (MC) e cálculos ab-initio de mecânica quântica (MQ). As configurações do líquido foram geradas para cálculos posteriores de MQ.Umimportante aspecto é a complexa interação interatômica do par He-He. Usando potenciais parametrizados, as simulações clássicas de MC foram efetuadas para sistemas alcalinos (Na, Rb, Cs e Na2) em He líquido e as condições foram T = 3 K e p = 1 atm. Estruturas estatisticamente descorrelacionadas formadas por um elemento alcalino central, envolvido pela primeira camada de solvatação completa, são amostradas e submetidas em um cálculo do espectro com DFT dependente do tempo usando diferentes funcionais híbridos e conjuntos de bases. Usando os funcionais PBE1PBE e O3LYP com conjuntos de bases extensos obtemos o deslocamento espectral em excelente concordância com os resultados experimentais para os sistemas de um único átomo alcalino. Para comparação, também usamos um modelo de agregado com 14 átomos de He em volta do átomo alcalino obtendo excelentes resultados também. O raio do modelo de agregado convergiu para perto do máximo da primeira camada de solvatação da função de distribuição radial. Um ponto adicional a ser considerado é o cálculo da largura da linha obtido com a simulação em He líquido que é discutida neste trabalho. Para o átomo de Rb, a energia de excitação em He líquido é em torno de -18,9 nm. Com a simulação em ambiente de He líquido obtivemos os melhores resultados entre -16,3 nm e -23,3 nm. O valor do deslocamento espectral usando o modelo de agregado ficou entre os -17,3 nm e - 22,3 nm. Os dois modelos apresentam o mesmo raio da bolha, por volta de 6-7Å. Para outros sistemas, como Na e Cs, encontramos a mesma convergência entre o modelo de agregado, a simulação do He líquido e os resultados experimentais. Para o sistema contendo Na2, os valores obtidos ficaram em boa concordância com os valores experimentais.