Yes, one output is fine. You could use a nn.Sigmoid layer as the model’s output and BCELoss.
In case you would like to use logits, you could skip the sigmoid and just use BCEWithLogitsLoss. The second approach should be numerically more stable.
In the post before I said logsigmoid which was wrong regarding the BCELoss.