Simulation of charge transfer in a tunnel junction: approaching the sub-e scale