Path integral Monte Carlo simulation of the absorption spectra of an Al atom embedded in helium

We use a multilevel path integral Monte-Carlo (PIMC) method to simulate the arrangement of He atoms around a single Al atom doped in a He cluster. High-level ab initio Al-He pair potentials and a Balling and Wright pairwise Hamiltonian model are used to describe the full potential and the electronic asymmetry arising from the open-shell character of the Al atom in its ground and excited electronic states. Our calculations show that the doping of the Al 3p electron strongly influences the He packing. The results of the PIMC simulation are used to predict the electronic excitation spectrum of an Al atom embedded in He clusters. With inclusion of tail corrections for the ground and excited states potentials, the calculated 3d<--3p spectrum agrees reasonably well with the experimental spectrum. The blueshift of the calculated spectrum associated with the 4s<--3p transition of solvated Al is about 25 nm (2000 cm-1) larger than seen in experiments on Al embedded in bulk liquid He. We predict that the spectrum associated with the 4p<--3p transition will be blueshifted by approximately 7000 cm-1 (nearly 1 eV).