We present a detection scheme to search for QCD axion dark matter, that is based on a direct interaction between axions and electrons explicitly predicted by DFSZ axion models. The local axion dark matter field shall drive transitions between Zeeman-split atomic levels separated by the axion rest mass energy $m_a c^2$. Axion-related excitations are then detected with an upconversion scheme involving a pump laser that converts the absorbed axion energy ($\sim $ hundreds of $\mu$eV) to visible or infrared photons, where single photon detection is an established technique. The proposed scheme involves rare-earth ions doped into solid-state crystalline materials, and the optical transitions take place between energy levels of $4f^N$ electron configuration. Beyond discussing theoretical aspects and requirements to achieve a cosmologically relevant sensitivity, especially in terms of spectroscopic material properties, we experimentally investigate backgrounds due to the pump laser at temperatures in the range $1.9-4.2$ K. Our results rule out excitation of the upper Zeeman component of the ground state by laser-related heating effects, and are of some help in optimizing activated material parameters to suppress the multiphonon-assisted Stokes fluorescence.

Buffer gas cooling with a $^4$He gas is used to perform laser-absorption spectroscopy of the $^{12}$C$_2$H$_2$ ($\nu_1+\nu_3$) band at cryogenic temperatures. Doppler thermometry is first carried out to extract translational temperatures from the recorded spectra. Then, rotational temperatures down to 20 K are retrieved by fitting the Boltzmann distribution to the relative intensities of several ro-vibrational lines. The underlying helium-acetylene collisional physics, relevant for modeling planetary atmospheres, is also addressed. In particular, the diffusion time of $^{12}$C$_2$H$_2$ in the buffer cell is measured against the $^4$He flux at two separate translational temperatures; the observed behavior is then compared with that predicted by a Monte Carlo simulation, thus providing an estimate for the respective total elastic cross sections: $\sigma_{el}(100\ {\rm K})=(4\pm1)\cdot 10^{-20}$ m$^{2}$ and $\sigma_{el}(25\ {\rm K})=(7\pm2)\cdot 10^{-20}$ m$^{2}$.