What is the magnetic field induced near an antiferromagnet? A ferromagnet has a nonvanishing total magnetization which results in a finite magnetic field induced outside the magnet. For an Heisenberg antiferromagnet instead, the total magnetization is zero and therefore the "macroscopic" magnetic field induced is also zero in first approximation. 
However, since an Heisenberg antiferromagnet can be thought in first approximation as a discrete arrangement of antiparallel magnetic moments (e.g., the Néel state), it is reasonable to imagine that on a very short scale, very close to the system (assume a planar surface), the magnetic field induced is non zero. 
Ferromagnet: $\quad\quad\quad\quad\quad\quad\quad\uparrow\,\uparrow\,\uparrow\,\uparrow\,\uparrow\,\uparrow\,\uparrow\,\uparrow\,\uparrow\,\uparrow\,\uparrow\,\uparrow\,\uparrow\,\uparrow\,\uparrow\,\uparrow$
Antiferromagnet (Néel state): $\quad\uparrow \, \downarrow \,\uparrow \, \downarrow \,\uparrow \, \downarrow \,\uparrow \, \downarrow \,\uparrow \, \downarrow \,\uparrow \, \downarrow \,\uparrow \, \downarrow \,\uparrow \, \downarrow \,$
So my questions are: 
1) Is this reasoning correct? It is safe to assume a Néel state to calculate the field outside the antiferromagnet? 
2) Is it possible to have a rough estimate of how big is this field and how fast it decays with the distance from the surface?
Edit regarding 1) The Neel is an approximate groundstate for an AFM, not the exact (which can be obtained using the Bethe ansatz). Thus, if you calculate the AFM field from the exact groundstate (suppose you could do it), you will get an "exact" field which is probably different from the approximate field calculated using the Néel state. 
Is the difference between two solutions relevant in condensed matter systems?
How good is this approximation?
 A: *

*It is accurate to assume that an antiferromagnetic material has the Néel form as long as it is experimentally verified (e.g. with Magnetic Neutron Scattering). On the theoretical side, whether the Neel state is the ground state depends on the Hamiltonian and dimensionality under consideration. For the Ising model it is an eigenstate, but for the Heisenberg model it is not. If you want an estimate to how well the Neel state works, you can approach it via Mean-field theory and estimate the rms fluctuations. Additionally, you can examine the effects of the excited states perturbatively on the Neel order parameter. Much of this is outlined in this pdf.  


Note: However, as is the case for ferromagnets, most antiferromagnets have microscopic domains. Because there is no net dipole moment, getting these domains to line up is much more difficult than for a ferromagnet. Additionally, many antiferromagnets have structures that are more complicated than a simple cubic arrangement of spins such as in the classical Néel state. For example, see the image below from here:



*The overall dipole moment of an antiferromagnetic domain is zero, so the next order term is the quadrupole term. The quadrupole magnetic field goes as ${r^{-4}}$, as compared to the dipole field which goes as ${r^{-3}}$. To get the exact form you would need to calculate the full quadrupole matrix, something that is described in the section on multipole expansion in Griffiths. In any case, the fields are even weaker when taking domain sizes into consideration. Nonetheless, although the fields from antiferromagnets are very weak, you can still measure them. The idea technique for this is scanning SQUID microscopy. Which put simply, is to put a superconducting junction near the material to detect very weak magnetic fields. This was done for example [here] (http://journals.aps.org/prl/abstract/10.1103/PhysRevLett.42.257).


Edit: Noticed OP wanted a quantum mechanical explanation. The Néel state is the (approximate) ground state of many quantum mechanical models, but in particular the Ising, Heisenberg, and XY models for the correct sign of the magnetic coupling term $J$ (see below). 
$$H=J\sum_{i,j}\mathbf{\vec{S}}_i\cdot\mathbf{\vec{S}}_j$$
However, dimensionality, range of interaction, and finite system sizes are key. In 1D you run into the Mermin-Wagner theorem, and the Néel state is almost never the correct ground state; you can get more details in Ashcroft and Mermin and the Oxford series textbook on Magnetism in Solids.
With regards to how to describe the field coming from a quantum mechanical magnetic moment (either from spin or angular momentum), the dipole formula for the field is exact for spins. That is, a quantum mechanical spin is a pure dipole. With regards to angular momentum, you technically would need to know how to compute the exact current operator $\mathbf{j}=$ together with the Biot-Savart law to get the exact form of the magnetic field. However, because electrons with large angular momenta are typically localized, it is appropriate to use a multipole expansion and keep just the dipole term.  
