Quasiparticle excitations in FeSe were studied by means of specific heat (C) measurements on a high-quality single crystal under rotating magnetic fields. The field dependence of C shows threestage behavior with different slopes, indicating the existence of three gaps (∆1, ∆2, and ∆3). In the low-temperature and low-field region, the azimuthal-angle (φ) dependence of C shows a fourfold symmetric oscillation with sign change. On the other hand, the polar-angle (θ) dependence manifests as an anisotropy-inverted two-fold symmetry with unusual shoulder behavior. Combining the angle-resolved results and the theoretical calculation, the smaller gap ∆1 is proved to have two vertical-line nodes or gap minima along the kz direction, and is determined to reside on the electron-type ε band. ∆2 is found to be related to the electron-type δ band, and is isotropic in the ab-plane but largely anisotropic out of the plane. ∆3 residing on the hole-type α band shows a small out-of-plane anisotropy with a strong Pauli-paramagnetic effect.Superconducting (SC) gap structures are intimately related to the pairing mechanism, which is pivotal for high-temperature superconductors. This issue is crucial for FeSe because of the unexpectedly high superconducting temperature, T c , in this system. Although the initial T c is below 10 K [1], it can be easily enhanced to 37 K under pressure, [2,3] and to over 40 K by intercalating spacer layers [4]. Recently, a monolayer of FeSe grown on SrTiO 3 has even shown a sign of T c over 100 K [5,6]. FeSe manifests some intriguing properties, including a nematic state without long-range magnetic order [7], crossover from Bardeen-Cooper-Schrieffer (BCS) to Bose-Einstein-condensation (BEC) [8], and a Diraccone-like state [9][10][11], which are crucial to understanding high-T c superconductivity.Efforts have been made to elucidate FeSe's gap structure, and the presence of nodes [8,12] or deep minima [13][14][15][16] have been proposed. Even if the multi-gap structure is established [17,18], gap nodes or minima must still be located in the corresponding bands. Unfortunately, there have been few reports on this issue except for the recent Bogoliubov quasiparticle interference (BQPI) experiments, which found gap minima in both the α and ε bands [19]. However, there is still no bulk evidence of the locations of nodes or gap minima, and also no information about the gap from the δ band. More importantly, details of gap structure, including the three-dimensional (3D) locations of gap nodes or minima, remain unexplored. To solve these issues, a bulk technique capable of probing quasiparticle (QP) excitations with 3D angular resolution is needed. Field-angle-resolved specific heat (ARSH) measurement is an ideal tool for probing the density of states (DOS) of QPs without interference from surface effects; it is angle resolved because the lowlying QP excitations near the gap nodes (minima) are field-orientation dependent [20]. ARSH measurements have been well applied to investigating the locations and types of nod...