The mechanisms underlying the deformation of nanovoids in Ta single crystals are analyzed when they are subjected to cyclic uniaxial deformation using numerical simulations. Boundary and cell-size effects have been mitigated by means of the Quasicontinuum (QC) method. We have considered ≈1 billion-atom systems containing 10.9 nm voids. Two kinds of simulations have been performed, each characterized by a different boundary condition. First, we compress the material along the nominal [0 0 1] direction, resulting in a highly symmetric configuration that results in high stresses. Second, we load the material along the high-index direction to confine plasticity to a single slip system and break the symmetry. We find that the plastic response under these two conditions is strikingly different, the former governed by dislocation loop emission and dipole formation, while the latter is dominated by twinning. We calculate the irreversible plastic work budget derived from a loading–unloading cycle and identify the most relevant yield points. These calculations represent the first fully three-dimensional, fully non-local simulations of any body-centered cubic metal using QC.