The µeV-mass axion is one of the most promising candidates for cold dark matter, and remains to be a well-motivated solution to the CP problem of Quantum Chromodynamics (QCD) via the Peccei-Quinn mechanism. In this paper, we propose a novel method to detect the dark-matter axions in our galaxy via the resonant emission |e → |g + γ + γ + a (or absorption a + |e → |g + γ + γ ) in an atomic system with superradiance, where |e and |g stand for the excited and ground energy levels of atoms, respectively. A similar process via |e → |g + γ + a (or a + |e → |g + γ) is also put forward to probe the axion-electron coupling. For the nominal experimental setup assuming a background-free environment, most of the parameter space for typical QCD axion models can be covered with parahydrogen molecules or ytterbium atoms. However, the background in a realistic experimental setup remains to be a major issue that needs to be solved in future studies. Searching for better atomic or molecular candidates may be required for a bigger signal-to-noise ratio. PACS numbers: 93.35.+d, 98.35.Gi, 21.60.CsIntroduction.-More than forty years ago, Peccei and Quinn (PQ) proposed an appealing solution to the CP problem of Quantum Chromodynamics (QCD) by introducing a dynamical scalar field and imposing a global U(1) PQ symmetry on the whole Lagrangian [1,2]. It was Weinberg [3] and Wilczek [4] who shortly discovered that a Nambu-Goldstone boson, i.e., the axion, arose from the spontaneous breaking of the PQ symmetry at some high-energy scale. Although the original model with the PQ symmetry spontaneously broken at the electroweak scale Λ EW ≡ 10 2 GeV has been ruled out, the "invisible" axion models, such as the KSVZ model [5,6] and the DFSZ model [7,8], are still attracting a lot of attention. Apart from providing a solution to the strong CP problem, these models also indicate that axions can be a good candidate for cold dark matter in our Universe [9][10][11][12][13], and can be detected in realistic experiments [14][15][16][17]]. An excellent overview of possible experimental methods to probe axions has recently been presented in Ref. [18].Due to the instanton effects, the axion acquires a small mass m a from the explicit PQ symmetry breaking at low energies. In a generic axion model, the axion mass m a and decay constant f a (i.e., the energy scale of spontaneous PQ symmetry breaking) are related to each other via m a · f a 6.0 µeV · 10 12 GeV. The experimental searches for dark matter axions mainly rely on their couplings to photons and fermions [19], i.e.,