The error message states that the weights of shape [16, 1, 5, 5] ([nb_filters, in_channels, kernel_h, kernel_w]) do not match the input size of your data.
You would have to fix both issues:
images = images.unsqueeze(1).float()
outputs = model(images)

Could you try this code again?