Aims: Hydrogen deficient white dwarfs are characterized by very dense, fluid-like atmospheres of complex physics and chemistry that are still poorly understood. The incomplete description of these atmospheres by the models results in serious problems
with the description of spectra of these stars and subsequent difficulties in derivation of their surface parameters. Here, we address the problem of infrared (IR) opacities in the atmospheres of cool white dwarfs by direct $ab$ $initio$ simulations of IR absorption of dense helium. Methods: We applied state-of-the-art density functional theory-based quantum molecular dynamics simulations to obtain the time evolution of the induced dipole moment. The IR absorption coefficients were obtained by the Fourier transform of the dipole moment time autocorrelation function. Results: We found that a dipole moment is induced due to three- and more-body simultaneous collisions between helium atoms in highly compressed helium. This results in a significant IR absorption that is directly proportional to $rm rho_{rm He}^3$, where $rho_{rm He}$ is the density of helium. To our knowledge, this absorption mechanism has never been measured or computed before and is therefore not accounted for in the current atmosphere models. It should dominate the other collisionally induced absorptions (CIA), arising from $rm H-He$ and $rm H_2-He$ pair collisions, and therefore shape the IR spectra of helium-dominated and pure helium atmosphere cool white dwarfs for $rm He/H>10^4$. Conclusions: Our work shows that there exists an unaccounted IR absorption mechanism arising from the multi-collisions between He atoms in the helium-rich atmospheres of cool white dwarfs, including pure helium atmospheres. This absorption may be responsible for a yet unexplained frequency dependence of near- and mid- IR spectra of helium-rich stars.
In Density Functional Theory simulations based on the LAPW method, each self-consistent field cycle comprises dozens of large dense generalized eigenproblems. In contrast to real-space methods, eigenpairs solving for problems at distinct cycles have
either been believed to be independent or at most very loosely connected. In a recent study [7], it was demonstrated that, contrary to belief, successive eigenproblems in a sequence are strongly correlated with one another. In particular, by monitoring the subspace angles between eigenvectors of successive eigenproblems, it was shown that these angles decrease noticeably after the first few iterations and become close to collinear. This last result suggests that we can manipulate the eigenvectors, solving for a specific eigenproblem in a sequence, as an approximate solution for the following eigenproblem. In this work we present results that are in line with this intuition. We provide numerical examples where opportunely selected block iterative eigensolvers benefit from the reuse of eigenvectors by achieving a substantial speed-up. The results presented will eventually open the way to a widespread use of block iterative eigensolvers in ab initio electronic structure codes based on the LAPW approach.
