Majorana zero modes (MZMs) obey non-Abelian statistics and are considered building blocks for constructing topological qubits 1,2 . Iron-based superconductors with topological band structures have emerged as promising hosting materials, since isolated candidate MZMs in the quantum limit have been observed inside the topological vortex cores 3-9 . However, these materials suffer from issues related to alloying-induced disorder, uncontrolled vortex lattices 10-13 and a low yield of topological vortices 5-8 . Here, we report the formation of an ordered and tunable MZM lattice in naturally-strained stoichiometric LiFeAs by scanning tunneling microscopy/spectroscopy (STM/S). We observe biaxial charge density wave (CDW) stripes along the Fe-Fe and As-As directions in the strained regions. The vortices are pinned on the CDW stripes in the As-As direction and form an ordered lattice. We detect more than 90 percent of the vortices to be topological and possess the characteristics of isolated MZMs at the vortex center, forming an ordered MZM lattice with the density and the geometry tunable by an external magnetic field. Remarkably, with decreasing the spacing of neighboring vortices, the MZMs start to couple with each other. Our findings provide a new pathway towards tunable and ordered MZM lattices as a platform for future topological quantum computation.The long quest for Majorana excitations [14][15][16][17][18][19] has been fueled by recent experimental progress in material platforms where the localized MZM has been observed 16,[20][21][22][23] . Among these platforms, ironbased superconductors are considered promising for the observation of clean and robust MZMs 3-5 . FeTe0.55Se0.45 was the first iron-based superconductor found to possess topological surface states (TSS) 24,25 . The largeratio, where Δ is the superconducting (SC) energy gap and 𝐸 𝐹 is the Fermi energy, allows the observation of an isolated MZM in the center of a topological vortex 4 , which exhibits nearly quantized plateaus in differential conductance 9 . Later on, stoichiometric iron pnictides were investigated and the existence of MZMs has been verified in CaKFe4As4 (ref. 8 ) and LiFeAs (ref. 6 ). Very recently, 1D dispersing Majorana modes and MZMs have been observed in the layered 2D superconductors with topological band structures [26][27][28] . However, these platforms suffer from alloyinginduced disorder, uncontrolled vortex lattice and a low yield of topological vortices, hindering their potential applications.We start with LiFeAs crystal due to its homogeneous bulk electronic structure, well-defined cleaving surface and rich topological band structures 29 . LiFeAs belongs to the tetragonal crystal system, with lattice constants a=b=3.8 Å and c=6.3 Å (Fig. 1a). The As-As atomic direction defines the high-