The effect of the conformational disorder in the electronic ground state may be modelled either by performing ground-state molecular dynamics (GSMD) simulations (usually applying a thermostat to mimic temperature of environment) or by sampling Wigner function corresponding to the vibrational ground state of the electronic ground state. The vertical absoprtion spectra and exciton localization are then computed for a set of selected geometries. In the case of large (multichromophoric) molecules or molecular aggregates, the GSMD approach is usually used, since the Wigner function is computed in practice in the harmonic approximation and, hence, the Wigner sampling may be problematic for anharmonic modes, which are present in the complex systems. Recently, we have investigated the effect of the conformational disorder on the exciton states of a model azobenzene tetramer (with transazobenzene units assembled in an H-aggregate fashion) by time-dependent density functional theory (TD-DFT) calculations perfromed at nuclear configurations generated by the GSMD DFT approach. We have found that the ground-state conformational disorder leads to partial localization of the ππ * excitons (with a localization degree being dependent on temperature) and strong localization (to a single monomer) of the nπ * excitons. Here we should recall that the absorption spectrum of trans-azobenzene includes a weak lower-energy nπ * band and an intense higher-energy ππ * band. Further, illumination of molecular assemblies with light excites them from the electronic ground state to the electronically excited state(s), inducing subsequent exciton dynamics. These dynamics require, in general, nonadiabatic treatment, i.e., coupled electron-nuclear dynamics involving many potential energy surfaces should be described. Over the past years the exciton dynamics in multichromophoric systems have been modelled using various ap-proaches to nonadiabatic dynamics, including surface hopping (SH), multiconfiguration time-dependent Hartree (MCTDH) and multi-layer MCTDH (ML-MCTDH), Ehrenfest and multiconfigurational Ehrenfest (MCE), and symmetrical quasi-classical windowing model applied to the classical Meyer-Miller vibronic Hamiltonian (SQC/MM). However, only few studies have been devoted to the exciton dynamics in azobenzene aggregates. Sangiogo Gil, Persico, and Granucci used surface hopping combined with an exciton model to simulate Frenkel exciton dynamics in an azobenzene dimer (azobenzenophane) and an azobiphenil monolayer. Recently, we have employed a supermolecule approach to investigate exciton dynamics in several azobenzene tetramers (both free and embedded in a SAM-like environment) after ππ * excitation. The exciton evolution was explored using a transition density matrix analysis, allowing one to judge on spatial localization of excitons during surface hopping dynamics. In this work, we study the exciton dynamics after nπ * excitation, i.e., induced with visible light excitation, for the models introduced in our previous work on the photodynamics after ππ * excitation. While the ππ * absorption band of azobenzene is more intense than the nπ * band, the latter is located in the visible region of the electromagnetic spectrum, which is of preference for applications requiring lower energy photons, e.g., for use in biological systems. Indeed, one of the central goals of contemporary research in the field of photoswitches is to enable efficient isomerization with low-energy photons. Here, using nonadiabatic dynamics simulations, we study how the corresponding low-energy excitons formed in azobenzene aggregates evolve in time.