We consider quasinormal modes and Hawking radiation of four dimensional asymptotically flat black holes in the most general up to-cubic-order-in-curvature dimension-independent theory of gravity that shares its graviton spectrum with the Einstein theory on constant curvature backgrounds. We show that damping rate and real oscillation frequencies of quasinormal modes for scalar, electromagnetic and Dirac fields are suppressed once the coupling with the cubic term is on. The intensity of Hawking radiation is suppressed as well, leading to, roughly, one order longer lifetime at a sufficiently large coupling constant.

Recently a $D$-dimensional regularization approach leading to the non-trivial $(3+1)$-dimensional Einstein-Gauss-Bonnet (EGB) effective description of gravity was formulated in [D. Glavan and C. Lin, Phys. Rev. Lett. 124, 081301 (2020)] which was claimed to bypasses the Lovelock's theorem and avoids Ostrogradsky instability. Later Aoki, Gorji and Mukohyama [arXiv:2005.03859] formulated a well-defined EGB theory which is four-dimensional owing to introduction of additional (scalar) degrees of freedom. A number solutions of the first naive approach are also solutions of the well-defined theory. Here we calculate quasinormal modes of scalar, electromagnetic and gravitational perturbations and find the radius of shadow for spherically symmetric and asymptotically flat black holes in this theory. We show that the black hole is gravitationally stable when ($-16 M^2<\alpha \lessapprox 0.6 M^2$). The instability in the outer range is the eikonal one and it develops at high multipole numbers. The radius of the shadow $R_{Sh}$ obeys the linear law with a remarkable accuracy.

Classical (quasinormal) and quantum (Hawking) radiations are investigated for test fields in the background of a four dimensional, spherically symmetric and asymptotically flat black hole in the Einstein-dilaton-Gauss-Bonnet (EdGB) theory. The geometry of the EdGB black hole deviates from the Schwarzschild geometry only slightly. Therefore, here we observe that the quasinormal spectrum also deviates from its Schwarzschild limit at most moderately, allowing for a $9\%$ decrease in the damping rate and up to a $6\%$ decrease in the real oscillation frequency. However, the intensity of Hawking radiation of an electromagnetic field turned out to be much more sensitive characteristic than its quasinormal spectrum, allowing for a $54\%$ increase of the energy emission rate. The analytical formula for the eikonal regime of quasinormal modes is derived for test fields and it is shown that the correspondence between the eikonal quasinormal modes and null geodesics is indeed fulfilled for test fields, but is not expected for the gravitational one.