Abstract Borexino is a 300 ton sub-MeV liquid scintillator solar neutrino detector which has been running at the Laboratori Nazionali del Gran Sasso (Italy) since 2007. Thanks to its unprecedented radio-purity, it was able to measure the flux of ⁷Be, ⁸B, pp, and pep solar neutrinos and to detect geo-neutrinos. A reliable simulation of the detector is an invaluable tool for all Borexino physics analyses. The simulation accounts for the energy loss of particles in all the detector components, the generation of the scintillation photons, their propagation within the liquid scintillator volume, and a detailed simulation of the electronics chain. A novel efficient method for simulating the external background which survives the Borexino passive shield was developed. This technique allows to reliably predict the effect of the contamination in the peripheral construction materials. The techniques developed to simulate the Borexino detector and their level of refinement are of possible interest to the neutrino and dark matter communities, especially for current and future large-volume liquid scintillator experiments.