A critical comparison of three polarization based approaches with the fields in AlN/GaN multiple quantum wells (MQWs) systems proved that they give identical results. The direct density functional theory (DFT) results, i.e., the fields, are in qualitative agreement with data obtained within the polarization theory. The results of DFT calculations of an AlN/GaN MQW system were used in the projection method to obtain a spatial distribution of the bands in the structure with atomic resolution. In parallel, the plane averaged and c-smoothed potential profiles obtained from the solution of the Poisson equation were used to determine the electric field in the multiquantum well structures and the magnitude of dipole layers at the AlN/GaN heterostructures. The dipole layers cause potential jumps of about 2.4 V that seriously affects the band offsets. The presence of the dipole layer is in good agreement with the potential measurements by electron holography. It was shown that the wells of the width up to 4 Ga layers behave as potential minima, but the wider layers behave as standard quantum wells. The barriers up to 3 Al layers do not localize the carriers. It is shown that the Quantum Confined Stark Effect causes a huge decrease of their energies and oscillator strengths of the optical transitions, especially for wider structures. For wider wells, the strengths fall much faster for perpendicular polarization which indicates the important role of the anisotropic band offsets. A direct simulation shows that the band offset for the valence band crystal field split off hole states, i.e., p z states are different from heavy and light hole (i.e., p ? ¼ p x p y ) states being equal to valence band offset ðVBOÞ ? ¼ 0:85 eV and rough estimate of VBO II ffi 0:5 eV, respectively. These values are in good agreement with the recently reported measurement of AlN/ GaN offsets. V C 2013 AIP Publishing LLC. [http://dx.