Monte Carlo simulations in high-energy physics experiments face the non-trivial task of simulating realistically the response of individual detector components, while keeping the costs in terms of CPU time reasonably low. Such simulation procedures are called digitization and have to incorporate detector physics in as much detail as possible, while performing fast enough. Here, we present our approach to the simulation of the Belle II vertex detector (VXD) using the ILC software framework. We simulate the response of DEPFET (DEPleted Field Effect Transistor) pixel detectors (PXD) and double-sided silicon micro-strip detectors (SVD) in the presence of a magnetic field. To achieve sufficient performance, we use a combination of fast numerical techniques and reasonable simplification of in-detector physics. The simulation itself is divided into 3 steps: particle propagation through matter (using Geant4 with a detailed implementation of the VXD geometry), charge collection in the silicon detectors (digitization) and clustering. The second and third steps are performed in separate reconstruction modules -Marlin modules: SiPxlDigi (pixel detectors) and SiStripDigi (micro-strip detectors). The following physics processes are considered: continuous energy loss fluctuations, generation of e-h pairs, drift in electric field, diffusion, Lorentz shift, mutual crosstalk (strips), read-out/geometric pitch effect (strips), electronics/digital noise. The clustering procedure is based on the center-of-gravity and analog head-tail algorithms. All effects have been studied and compared to test beam data in order to validate the algorithms and to determine the relative importance of individual processes.