To understand the nematicity in Fe-based superconductors, nontrivial k-dependence of the orbital polarization (∆Exz(k), ∆Eyz(k)) in the nematic phase, such as the sign reversal of the orbital splitting between Γ-and X,Y-points in FeSe, provides significant information. To solve this problem, we study the spontaneous symmetry breaking with respect to the orbital polarization and spin susceptibility self-consistently. In FeSe, due to the sign-reversing orbital order, the hole-and electron-pockets are elongated along the ky-and kx-axes respectively, consistently with experiments. In addition, an electron-pocket splits into two Dirac cone Fermi pockets with increasing the orbital polarization. The orbital-order in Fe-based superconductors originates from the strong positive feedback between the nematic orbital order and spin susceptibility. The spontaneous symmetry breaking from C 4 -to C 2 -symmetry, so called the electronic nematic transition, is one of the fundamental unsolved electronic properties in Fe-based superconductors. To explain this nematicity, both the spin-nematic scenario [1][2][3][4][5][6] and the orbital order scenario [7][8][9][10][11][12][13] have been studied intensively. Above the structural transition temperatures T str , large enhancement of the electronic nematic susceptibility predicted by both scenarios [1,10] is actually observed by the measurements of the softening of the shear modulus C 66 [1,10,14,15], Raman spectroscopy, [16][17][18], and in-plane resistivity anisotropy ∆ρ [19].To investigate the origin of the nematicity, FeSe (T c = 9 K) is a favorable system since the electronic nematic state without magnetization is realized below T str = 90 K down to 0 K. Above T str , the antiferromagnetic fluctuations is weak and T -independent according to the NMR [20,21] and neutron scattering [22][23][24] The nontrivial electronic state below T str gives a crucial test for the theories proposed so far. In the orthorhombic phase with (a − b)/(a + b) ∼ 0.3%, large orbital-splitting |E xz − E yz | of order 50 meV is observed at X,Y-points by ARPES studies in BaFe 2 As 2 [30], NaFeAs [31], and FeSe [32][33][34][35][36][37][38][39][40]. Especially, noticeable deformation of the Fermi surfaces (FSs) with C 2 -symmetry is realized in FeSe, because of the smallness of the Fermi momenta. In FeSe, Ref. [38] reports that the orbital splitting E xz − E yz is positive at Γ-point, whereas it is negative at X,Y-points. This sign-reversing orbital splitting is not realized in the non-magnetic orthorhombic phase in NaFeAs [31]. In addition, the e-FS1 at X-point is deformed to two Dirac cone Fermi pockets in thin-film FeSe [37,40]. The aim of this study is to explain these nontrivial electronic states in the orbital-ordered states based on the realistic multiorbital Hubbard model. Microscopically, the orbital order is expressed by the symmetry breaking in the self-energy. In the mean-field level approximations, however, the self-energy is constant in k-space unless large inter-site Coulomb interactions are in...