  The crystal structure, magnetic ground state, and the temperature dependent
microscopic spin-spin correlations of the frustrated honeycomb lattice
antiferromagnet Na2Co2TeO6 have been investigated by powder neutron
diffraction. A three dimensional (3D) long-range antiferromagnetic ordering has
been found below TN ~ 24.8 K. The magnetic ground state, determined to be
zigzag antiferromagnet [k = (1/2 0 0)], occurs due to the competing exchange
interactions up to third nearest neighbors within the honeycomb lattice. The
exceptional existence of a limited magnetic correlation length along the c axis
(perpendicular to the honeycomb layers in the ab planes) has been found within
the 3D long-range ordered state, even at 1.8 K, well below TN. The observed
limited correlation along the c axis is explained by the disorder distribution
of the Na ions within the intermediate layers. The reduced ordered moments
mCo(1) = 2.77(3) muB/Co2+ and mCo(2) = 2.45(2) muB/Co2+ at 1.8 K reflect the
persistence of spin fluctuations in the ordered state. Above TN, the presence
of short-range magnetic correlations, manifested by broad diffuse magnetic
peaks, has been found. Reverse Monte Carlo analysis of the experimental diffuse
magnetic scattering data reveals that the spin correlations are mainly confined
within the 2D honeycomb layers (ab plane) with a correlation length of ~ 12
{\AA} at 25 K. The nature of spin arrangements is found to be similar in both
the short-range and long-range ordered magnetic states. This implies that the
short-range correlation grows with decreasing temperature and leads to the
zigzag AFM ordering at T < TN. The present study provides a comprehensive
picture of the magnetic correlations over the temperature range above and below
TN and their relation to the crystal structure. The role of intermediate soft
Na-layers on the magnetic coupling between honeycomb planes is discussed.
