moment d a . In the linear case, i. e., when the laser intensity is significantly lower than the saturation intensity, the results for the two-level model coincide exactly with those for the classical harmonic oscillator, see [2]. We further only consider lattices of atoms interacting via quasistatic near-field dipole forces. In general, this assumption is valid if, on one hand, the minimum separation l a between atoms is not too small, so that their atomic orbitals do not overlap, and, on the other hand, the lattice is finite and its overall dimensions are smaller than the laser wavelength λ . However, in most cases, in particular within the NNA, which is the most common one and used throughout this paper, the sufficient condition is much less restrictive, only requiring that the interatomic separation l a ≪ λ , which is the same limit as in the standard LLT [3]. The LF acting upon a given atom located at a point r is a superposition of the incident (“external”) field E$_{in}$ and the sum of the fields E$_{dp}$ ( r , r $^{′}$) from surrounding dipoles at all other lattice positions r ′ over the entire lattice: E$_{L}$ ( r ) = E$_{in}$ ( r ) + $^{$_{r}$′}$̸ = r ∑ lattice E$_{dp}$ ( r , r $^{′}$) , (1) where the near fields of the surrounding atomic dipoles p ( r $^{′}$) are dominated by the non- radiative (quasi-static) components [14], Here u ≡ ( r − r $^{′}$) / | r − r $^{′}$| is a unit vector along the line connecting the two atoms, ϵ is the background dielectric constant ( ϵ = 1 in vacuum and ϵ = 1 if a host medium is present). The atomic dipoles p ( r $^{′}$) are, in turn, induced by the LF acting upon them. In the case of linear optical response of a two-level atom, its induced dipole moment is where δ ≡ T ( ω − ω$_{0}$ ) is a dimensionless laser frequency detuning and T = 2 / Γ is the transverse relaxation time of the atom with a resonant homogeneous linewidth Γ. The condition that the electron orbitals of neighboring atoms do not overlap implies that l a ≫ | d a | /e and justifies our use of a semiclassical approach. Nonlinear effects in systems of resonant atoms could be included into this picture via the saturation nonlinearity of the two- level system. As we have recently demonstrated [1, 2], the nonlinearity enables interesting E$_{dp}$ ( r , r $^{′}$) = 3 u [ p ( r $^{′}$) · u ] − p ( r $^{′}$) ϵ | r ′ − r | 3 . (2) p ( r ) = − E$_{L}$ ( r ) 2 | d a | 2 ¯ h Γ( δ + i ) , (3) 6