Fermion masses and mixings in the supersymmetric Pati-Salam landscape from Intersecting D6-Branes

Recently, the complete landscape of three-family supersymmetric Pati-Salam models from intersecting D6-branes on a type IIA $\mathbb{T}^6/(\mathbb{Z}_2\times \mathbb{Z}_2)$ orientifold has been enumerated consisting of 33 independent models with distinct gauge coupling relations at the string scale. Here, we study the phenomenology of all such models by providing the detailed particle spectra and the analysis of the possible 3-point and the 4-point Yukawa interactions in order to accommodate all standard-model fermion masses and mixings. We find that only 17 models contain viable Yukawa textures to explain quarks masses, charged-leptons' masses, neutrino-masses, quarks' mixings and leptons' mixings. These viable models split into four classes, viz. a single model with 3 Higgs fields from the bulk and sixteen models with either 6, 9 or 12 Higgs from the $\mathcal{N}=2$ sector. The models perform successively better with the increasing number of Higgs pairs. Remarkably, the class of models with 12 Higgs naturally predicts the Dirac-type neutrino masses in normal ordering consistent with both the experimental constraints as well as the bounds from the swampland program.