Abstract. We investigate theoretically the nonlinear optical response of a two-dimensional supercrystal comprized of semiconductor quantum dots. An isolated quantum dot is modeled as a three-level ladder-like system with ground, one-exciton, and biexciton states. It is shown that the optical response of supercrystal demonstrate a rich nonlinear dynamics, including bistability, self-oscillations, and dynamical chaos. Supercrystals (SCs) comprising semiconductor quantum dots (SQDs) represent a class of new materials not existing in nature. Modern nanotechnology has in its disposal a variety of methods to design such systems [1]. Optical properties of SCs depend on the SQD's size, shape, and chemical composition, as well as on the lattice geometry and can be easily controlled [2].We conduct a theoretical study of the optical response of a two-dimensional SC. Due to a high density of SQDs in the SC, the SQD-SQD dipole-dipole interaction plays an important role in the SC's optical response, both linear and nonlinear. This interaction provides a positive feedback which, together with the SQD's nonlinearity, results in a rich optical dynamics of SC, including bistability, self-oscillations, and dynamical chaos.It is assumed that the SC undergoes an action of an external field of an amplitude Е0 and frequency ω0, which is quasiresonant with the SQD allowed transitions. A single SQD is modelled as a three-level ladder-like quantum system with the ground |1⟩, one-exciton |2⟩, and bi-exction |3⟩ states. Only the optical transitions |1⟩↔ |2⟩ and |2⟩↔|3⟩ are dipoleallowed. They are characterized by the transition dipole moments d21 and d32, transition frequencies ω21 and ω32, and the spontaneous decay constants γ21 and γ32. The frequency ω32 is down-shifted with respect to ω21 by an amount ΔB being the biexciton binding energy.The optical dynamics of an isolated SQD is governed by 3x3 density matrix. The field acting on a given SQD in the SC represents a sum of the external field and the field produced by all others SQDs in place of the given one. In this way, the total (retarded) dipole-dipole interaction is taken into account. As the mean dipole moment of an SQD depends on how strong the SQD is excited, the SQD-SQD interaction appears to depend on the current SQD's condition as well. Similar to a one-dimensional case [3], the near-zone part of the retarded SQD-SQD dipole-dipole interaction gives rise to a dynamic shift of the transition frequencies