The nuclear many-body problem for medium-mass systems is commonly addressed using wave-function expansion methods that build upon a second-quantized representation of operators with respect to chosen computational basis. While various options the basis are available, perturbatively constructed natural orbitals recently have been shown lead significant improvement in applications yielding faster...