The crystal structure and crystal chemistry of wardite, ideally NaAl 3 (PO 4 ) 2 (OH) 4 ·2H 2 O, was investigated by single-crystal neutron diffraction (data collected at 20 K) and electron microprobe analysis in wavelength-dispersive mode. The empirical formula of the sample used in this study is: (Na 0.91 Ca 0.01 ) Σ=0.92 (Al 2.97 Fe 3+ 0.05 Ti 0.01 ) Σ=3.03 (P 2.10 O 8 )(OH) 4 •1.74H 2 O. The neutron diffraction data confirm that the crystal structure wardite can be described with a tetragonal symmetry (space groupP4 1 2 1 2, a = b =7.0577(5) and c = 19.0559(5) Å at 20 K) and consists of sheets made by edge-sharing Na-polyhedra and Al-octahedra along with vertex-sharing Al-octahedra, parallel to (001), connected by P-tetrahedra and H-bonds to form a (001) layer-type structure, which well explains the pronounced {001} cleavage of the wardite crystals. The present data show that four crystallographically independent H sites occur into the structure of wardite, two belonging to a H 2 O molecule (i.e., H1-O6-H2) and two forming hydroxyl groups (i.e., O5-H3 and O7-H4). The location of the hydrogen atoms allow us to define the extensive network of H-bonds: the H-atoms belonging to the H 2 O molecule form strong H-bonds, whereas both the H-atoms belonging to the two independent hydroxyl groups form weak interactions with bifurcated bonding schemes. As shown by the root-mean-square components of the displacement ellipsoids, oxygen and hydrogen atoms have slightly larger anisotropic displacement parameters if compared to the other sites (populated by P, Al and Na). The maximum ratio of the max and min root-mean-square components of the displacement ellipsoids is observed for the protons of the hydroxyl groups, which experience bifurcated H-bonding schemes. A comparative analysis of the crystal structure of wardite and fluorowardite is also provided.