In our analysis, we have relied on the effective field theory expansion in the Einstein frame, which is expected to be meaningful during preheating.
Gravitational dark matter production can also be encoded in the DM non–minimal coupling to gravity [1]}, which corresponds to a specific choice of higher dimensional operators in our approach.
A related option is provided by
gravity–induced inflaton decay in Starobinsky–like models [2]}.
