Add edges to scale-free network igraph

This is my code:

network2 <- sample_pa(1000, m = 3, directed=F)
ecount(network2)

My goal is a network with 1000 vertices (done) and 3000 edges (I only get to 2994). Is there a way to add edges according to the Barabasi/Altert model and not just random edges?