Active regions (AR) appearing on the surface of the Sun are classified into α, β, γ, and δ by the rules of the Mount Wilson Observatory, California on the basis of their topological complexity. Amongst these, the δ-sunspots are known to be super-active and produce the most X-ray flares. Here, we present results from a simulation of the Sun by mimicking the upper layers and the corona, but starting at a more primitive stage than any earlier treatment. We find that this initial state consisting of only a thin sub-photospheric magnetic sheet breaks into multiple flux-tubes which evolve into a colliding-merging system of spots of opposite polarity upon surface emergence, similar to those often seen on the Sun. The simulation goes on to produce many exotic δ-sunspot associated phenomena: repeated flaring in the range of typical solar flare energy release and ejective helical flux ropes with embedded cool-dense plasma filaments resembling solar coronal mass ejections.Delta-sunspots are formed when two sunspots of opposite polarity magnetic field appear very close to each other and reside in the same penumbra, the radial filamentary structure outside the umbral region of the strongest magnetic fields. Strong shear and horizontal magnetic fields often exist at the polarity-inversion line separating the two polarities [1]. The subsurface processes which form the δ-sunspots are still debated. Early observational studies [2,3] propose that δ-sunspots form from collision-merging of topologically separate dipoles, while numerical simulations by [4,5] show that kink unstable magnetic flux-tube -helical field lines winding around a central axis -emerging from the subsurface can have a δ-sunspot like structure. More recently, attempts to model the δ-spot in the NOAA AR 11158 utilized a uniformly twisted sub-surface flux-tube initially buoyant in two adjacent regions along its length [6,7]. Also, [8] found a magnetic flux concentration resembling a δ sunspot in their stratified helical dynamo simulation. These studies did not report any flaring activity. On the other hand [9] initialized their simulation with two parallel flux-tubes each lying at a different depth from the surface and with a different value of the initial magnetic twist which later evolved into a δ-sunspot like structure and powered multiple reconnection events. Delta-sunspots are highly flare-productive -95% of the strongest (X-class) X-ray flares originate from these regions [3]. Using a realistic numerical simulation [10] showed that interaction between adjacent expanding magnetic bipoles pressing against each other can lead to the formation of strong current layers in the atmosphere which in turn lead to repeated flaring in the region. Here, we report on a three-dimensional magneto-hydrodynamic (MHD) simulation of the formation of a δ-sunpot like region as a result of the break-up of a cool magnetic layer embedded in the upper solar convection zone into several flux-tubes due to the growth of three-dimensional unstable modes excited in the layer. This instab...