We provide evidence that spin ferroquadrupolar (FQ) order is the likely ground state in the nonmagnetic nematic phase of stoichiometric FeSe. By studying the variational mean-field phase diagram of a bilinearbiquadratic Heisenberg model up to the 2nd nearest neighbor, we find the FQ phase in close proximity to the columnar antiferromagnet commonly realized in iron-based superconductors; the stability of FQ phase is further verified by the density matrix renormalization group. The dynamical spin structure factor in the FQ state is calculated with flavor-wave theory, which yields a qualitatively consistent result with inelastic neutron scattering experiments on FeSe at both low and high energies. We verify that FQ can coexist with C 4 breaking environments in the mean-field calculation, and further discuss the possibility that quantum fluctuations in FQ act as a source of nematicity. Superconductivity in the iron-based superconductors [1,2] is widely recognized to have spin fluctuations at its origin [3,4], as it develops after the suppression of columnar antiferromagnetism (CAFM) by doping or applied pressure on the parent compounds [5][6][7][8]. The CAFM phase is characterized by the magnetic Bragg peaks at wave vectors Q 1,2 = (π, 0)/(0, π) in the one-iron Brillouin zone, seen ubiquitously in different families of the iron pnictides and chalcogenides [5,9,10]. The discovery of superconductivity in stoichiometric FeSe thus came as a surprise, because the long-range magnetic order is conspicuously absent in this material [11][12][13][14][15][16]. Another important feature, universally observed across different families of iron-based superconductors, is the appearance of an electronic nematic phase [17][18][19][20], which spontaneously breaks the lattice C 4 rotational symmetry. Usually, nematicity appears in close proximity to magnetism above the Néel temperature; however, in FeSe, the nematic phase appears without any accompanying magnetism and coexists with superconductivity [12][13][14][15]. It is thus important to understand the origin of this nonmagnetic nematic phase, in particular, to gain insight into its effect on superconductivity.It turns out that magnetic order can be induced by applying hydrostatic pressure to FeSe [12][13][14]. It has also been suggested based on ab initio calculations that the nonmagnetic phase in FeSe lies in close proximity to the CAFM phase [21][22][23]. Further evidence of proximity to long-range magnetic order comes from inelastic neutron scattering (INS) experiments, which found large spectral weight at wave vectors Q 1,2 [24][25][26][27]. Two natural questions arise: In the theoretical phase diagram, is there a nonmagnetic phase that neighbors on the CAFM? And, furthermore, how does such a nonmagnetic phase give rise to nematicity?In an attempt to answer these questions, several theoretical scenarios have been proposed for nonmagnetic ground states that may appear as a result of frustration: a nematic quantum paramagnet [28], a spin quadrupolar state with wave vectors Q 1,2...