Put all possible broom::glance statistics of lm() model combinations with 9 variables into a dataframe in R

Question

As I am just learning R, I am not sure how to solve this. I am trying to get a data frame that shows me the following:

Model Number | adj.r.squared | sigma