We revisit a supersymmetric flavor model based on the symmetries $SU(2)_L \times A_4 \times Z_3 \times U(1)_R$, which extends the original Altarelli and Feruglio construction by introducing flavon and driving superfields responsible for the spontaneous breaking of the flavor symmetry in order to obtain non-zero reactor angle. The vacuum alignments of flavon fields are achieved through the minimization of the scalar potential derived from the superpotential. This setup leads to specific mass matrices for the charged leptons and neutrinos that are consistent with current experimental data, including the measured values of the lepton mixing angles and neutrino mass squared differences. We investigate whether the model can simultaneously accommodate successful thermal leptogenesis. In particular, we analyze the CP asymmetry generated in the decay of heavy Majorana neutrinos, the resulting lepton asymmetry, and its conversion to the baryon asymmetry through the electroweak sphalerons. However the CP asymmetry is zero, since the Dirac neutrino mass matrix is simple texture in the leading order for our model. Then we consider the next-to-leading order in Yukawa interactions of the Dirac neutrinos. Therefore, we can realize the baryon asymmetry of the universe at the present universe. By numerically scanning the parameter space, we identify the regions consistent with both neutrino oscillation data and the observed baryon asymmetry. In the specific case such that one of the couplings for the right-handed Majorana neutrinos is real parameter, the predicted lightest neutrino mass is at least $5$ meV and $15$ meV for the normal and inverted neutrino mass hierarchies, respectively. In addition, the range of the Majorana phases may be tested in future experiments.

We study baryogenesis in a hybrid inflation model which is embedded to the minimal supersymmetric model with right-handed neutrinos. Inflation is induced by a linear combination of the right-handed sneutrinos and its decay reheats the universe. The decay products are stored in conserved numbers, which are transported under the interactions in equilibrium as the temperature drops down. We find that at least a few percent of the initial lepton asymmetry is left under the strong wash-out due to the lighter right-handed (s)neutrinos. To account for the observed baryon number and the active neutrino masses after a successful inflation, the inflaton mass and the Majorana mass scale should be $10^{13}\,{\rm GeV}$ and ${\cal O}(10^{9}$-$10^{10})\,{\rm GeV}$, respectively.

We study a relic density of the flavino dark matter in modified Altarelli and Feruglio $A_4$ model which is respecting the $SU(2)_L\times A_4\times Z_3\times U(1)_R$ symmetry. We calculate the Lagrangian from the superpotential in the model. In estimating relic density, we consider the relevant interactions from the Lagrangian that realize the vacuum expectation value alignments and charged lepton masses where we assume that the supersymmetry breaking effects are small for flavon sector. As a result, we find the degenerate masses between the lightest ``flavon'' and ``flavino'', and only two parameters in the potential determines the relic density. Then allowed parameter space of these parameters are estimated from relic density calculation taking a constraint from lepton flavor violation into account. We also briefly discuss other dark matter physics such as direct detection, indirect detection and collider search.