We report infrared studies of AFe 2 As 2 (A¼Ba, Sr), two representative parent compounds of iron-arsenide superconductors, at magnetic fields (B) up to 17.5 T. Optical transitions between Landau levels (LLs) were observed in the antiferromagnetic states of these two parent compounds. Our observation of a ffiffiffi ffi B p dependence of the LL transition energies, the zero-energy intercepts at B ¼ 0 T under the linear extrapolations of the transition energies and the energy ratio (∼2.4) between the observed LL transitions, combined with the linear band dispersions in two-dimensional (2D) momentum space obtained by theoretical calculations, demonstrates the existence of massless Dirac fermions in the antiferromagnet BaFe 2 As 2 . More importantly, the observed dominance of the zeroth-LL-related absorption features and the calculated bands with extremely weak dispersions along the momentum direction k z indicate that massless Dirac fermions in BaFe 2 As 2 are 2D. Furthermore, we find that the total substitution of the barium atoms in BaFe 2 As 2 by strontium atoms not only maintains 2D massless Dirac fermions in this system, but also enhances their Fermi velocity, which supports that the Dirac points in iron-arsenide parent compounds are topologically protected. DOI: 10.1103/PhysRevLett.119.096401 In condensed matter, massless Dirac fermions (MDF), which represent a type of quasiparticles with linear energymomentum dispersions, provide a cornerstone for various quantum phenomena, such as the novel quantum Hall effect [1-3], Klein tunneling [4], and giant linear magnetoresistance [5]. Because of their crucial role in diverse quantum phenomena, searching for MDF in new systems is one of the most active research areas in condensed matter science. Generally, massless Dirac fermions with their valence and conduction band touching at a degenerate point (i.e., Dirac point) are protected by symmetries in solids [6][7][8][9][10]. However, the formation of magnetic order is always accompanied with time-reversal symmetry breaking and sometimes followed by crystalline symmetry lowering. Therefore, it is a challenge to achieve MDF in magnetic ground states of solid-state electronic systems [11,12]. Moreover, two-dimensional (2D) MDF, which were realized in 2D materials and on the surfaces of threedimensional (3D) topological insulators [1][2][3][13][14][15][16][17], have rarely been observed in the bulk of 3D crystals [11]. The discovery of iron-arsenide superconductors not only brings people a new class of unconventional superconductivity [18], but also sheds light on seeking 2D MDF in magnetic compounds [19,20].The parent compounds of iron-arsenide superconductors (PCIS) exhibit collinear antiferromagnetic (AFM) order at low temperature (T) [21,22]. After the AFM phase transition, the paramagnetic Brillouin zone of PCIS would be folded along the AFM wave vector connecting the hole-type Fermi surfaces (FSs) surrounding the Γ point and the electron-type FSs at the M point, which leads to the intersection between the electron...