As a much-enriched supplement to the previous review paper entitled the “Effective work functions for ionic and electronic emissions from mono- and polycrystalline surfaces” Prog. Surf. Sci. 83 (2008) 1–165, the present monograph summarizes a comprehensive and up-to-date database in Table 1, which includes more than ten thousands of experimental and theoretical data accumulated mainly during the last half century on the work functions (ϕ+, ϕe and ϕ−) effective for positive-ionic, electronic and negative-ionic emissions from mono- and polycrystalline surfaces of 88 kinds of chemical elements (1H–99Es), and also which includes the main experimental condition and method employed for each sample specimen (bulk or film) together with 490 footnotes. From the above database originating from 4461 references published to date in the fields of both physics and chemistry, the most probable values of ϕ+, ϕe and ϕ− for substantially clean surfaces are statistically estimated for about 600 surface species of mono- and polycrystals. The values recommended for ϕe together with ϕ+ and ϕ− in Table 2 are much more abundant in both surface species and data amount, and also they may be more reliable and convenient than those in popular handbooks and reviews consulted widely still today by great many workers, because the latter is based on less-plentiful data on ϕe published generally before ∼1980 and also because it covers no value recommended for ϕ+ and ϕ−. Consequently, Table 1 may be more advantageous as the latest and most abundant database on work functions (especially ϕe) for quickly referring to a variety of data obtained under specified conditions. Comparison of the most probable values of ϕe recommended for each surface species between this article and other literatures listed in Tables 2 and 3 indicates that consideration of the recent work function data accumulated particularly during the last ∼40 years is very important for correct analysis of these surface phenomena or processes concerned with either work function or its changes. On the basis of our simple model about the work function of polycrystal consisting of a number of patchy faces (1–i) having each a fractional area (Fi) and a local work function (ϕi), its values of both ϕ+ and ϕe are theoretically calculated and also critically compared with a plenty of experimental data. In addition, the “polycrystalline thermionic work function contrast” (Δϕ∗≡ϕ+−ϕe) well-known as the thermionic peculiarity inherent in every polycrystal is carefully analyzed as a function of the degree of monocrystallization (δm) corresponding to the largest (Fm) among Fi’s (Tables 4–6 and Fig. 1), thereby yielding the conclusions as follows: (1) Δϕ∗≃ const (>0) holds for the generally called “polycrystalline” surfaces (usually δm < 50%), (2) Δϕ∗ ranges from ∼0.3 eV (Pt) to 0.7 eV (Nb) depending upon the polycrystalline surface species, (3) in the case of the “submonocrystal” (50 <δm < 100%) tentatively named here, Δϕ∗ decreases parabolically down to zero as δm increases from ∼50% up to 100% (monocrystal), (4) Δϕ∗=0.0eV applies to a clean and smooth monocrystalline surface (δm≈ 100%) alone, (5) regarding negative ion emission, on the other hand, our theoretical prediction of Δϕ∗∗≡ϕ−−ϕe=0.0eV is experimentally verified to hold for any surface species under any surface conditions (Table 7), (6) every polycrystal (usually, δm < 50%) may be concluded in general to have a unique value of ϕe characteristic of its species with little dependence upon δm, (7) this conclusion affords us first a sound basis for supporting theoretically the experimental fact (Table 2) that every species of polycrystal has a nearly constant value of ϕe as well as ϕ+ (usually within the uncertainty of ±0.1 eV) depending little upon the difference in the surface components (Fi and ϕi) among specimens so long as δm < 50%, (8) on the contrary to polycrystal (δm < 50%), any submonocrystal (50 <δm < 100%) has such an anomaly that it does not possess the unique value of work function characteristic of the surface species itself, because its ϕe as well as ϕ+ changes considerably depending upon δm, (9) consequently, submonocrystal must be taken as another type (category) different from both poly- and monocrystals, (10) in this way, δm acts as the key factor mainly governing the work functions in the different mode between poly- and submonocrystals with δm lower and higher than the “critical point” of 50%, respectively, (11) on the contrary to δm, ϕm belonging to δm has a differential effect on both ϕ+ and ϕe, but their values remain nearly constant so long as δm < 50% and, thus interestingly, (12) the complicate governance of ϕ+ and ϕe by both δm and ϕm and also the anomaly of submonocrystal (cf. (8) above) observed first by our theoretical analysis may be considered as a new contribution to the work function studies developed to date. Together with brief comments and experimental conditions, typical data on ϕe and/or ϕ+ are summarized from the various aspects of (1) examination of the work function dependence upon the surface atom density of low-Miller-index monocrystals of typical metals such as Al, Ni, W and Re (Table 8), (2) demonstration of the above dependence usually called the “anisotropic work function dependence sequences” of both ϕe(110) >ϕe(100) >ϕe(111) and ϕ+(110) >ϕ+(100) >ϕ+(111) for various bcc-metals (e.g., Nb, Mo, Ta and W) exactly obeying the Smoluchowski rule (Table 9), (3) substantiation of both ϕe(111) >ϕe(100) >ϕe(110) for a variety of fcc-metals (except Al and Pb) and ϕ+(111) >ϕ+(100) >ϕ+(110) for Ni strictly following the above rule (Table 10), (4) verification of the quantitative relations between work function and surface energy and also melting point of the three low index planes of several metals (typically, Ni), (5) examination of the work function change (Δϕe) due to allotropic transformation from α to β or β to γ phase (Table 11) together with a concise outline of the Burgers orientation relationship, (6) evaluation of Δϕe due to liquefying (Table 12), (7) estimation of Δϕe due to transformation from ferro- to paramagnetic state (Table 13) in addition to a brief description of the Curie point dependence upon metastable metal film thickness above one monolayer, (8) estimation of Δϕe due to transition from normal to superconducting state (Table 14), (9) study of the work function dependence on the Wigner–Seitz radius and also comparison between its theoretical values (by Kohn) and experimental data (Fig. 2), (10) inspection of the annealing effect on work function for layers or films, (11) verification of the coincidence of work function values among different experimental methods, and (12) inquisition of the work function dependence upon the size of fine particles (∼20–100 Å in radius) studied by theory and experiment. •Ten thousands of work function data are summarized together with methods and conditions.•Most probable work function values are listed for 600 surface species of 88 elements.•The above values are discussed comparatively with data recommended previously.•Three kinds of effective work functions are examined by theory and experiment.•Work function changes at various critical temperatures are analyzed critically.