You can create a virtual field in the Field calculator for the grid layer and fill it with
aggregate( layer:='Your pointlayername', aggregate:=' concatenate', expression:=to_string(attributes()['attributenameofpointlayer']), filter:=intersects($geometry, geometry(@parent)), concatenator:=',')