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.