the Gauss-Bonnet (GB) invariant,
a choice which allows for the scalarization of vacuum Schwarzschild BH (we recall that
\(L_{GB}=-48 M^2/r^6\) 
in this case, with \(M\)  the BH mass).
This model has been extensively studied, starting with
[1]}, [2]}, [3]}, [4]}
where the first examples of scalarized BHs resulting from this type of mechanism have been
reported.
Further work includes the study of scalarized BHs in various extensions
of the initial framework
[5]}, [6]}, [7]}, [8]}, [9]}, [10]}, [11]}, [12]}, [13]}, [14]}, [15]}, [16]}, [17]}, [18]}, [19]}, [20]}, [21]}
and
the investigation of solutions' stability
[22]}, [23]}, [24]}, [25]};
furthermore,
partial analytical results are reported in Refs.
[26]}, [27]}, [28]}, [29]},
while scalarized, rotating BHs
are studied in
[30]}, [31]}, [32]}, [33]}, [34]}, [35]}, [36]}.
