TMAP - add 'heat map' style layer based on plot density

I have the below code which creates a map with two layers, 1) area borders and 2) points of interest on top.

Map <- tm_shape(Addresses)+
  tm_dots(col = "red", size = 0.2)+
  tm_shape(Boundary)+ 
  tm_borders(
    col = "black",
    lwd = 1,
    lty = "solid",
    alpha = NA,
    zindex = NA,
    group = NA)

I would like to create a heatmap layer broadly similar to the below style, i.e. illustrating where point density is highest and lowest - can anyone advise how to add this to the above code to achieve this effect? Thank you.