Spin-s wavefunctions with algebraic order

We generalize the Gutzwiller wavefunction for s = 1/2 spin chains to construct a family of wavefunctions for all s > 1/2. Through numerical simulations, we demonstrate that the spin spin correlation functions for all s decay as a power law with logarithmic corrections. This is done by mapping the model to a classical statistical mechanical model, which has coupled Ising spin chains with long range interactions. The power law exponents are those of the Wess Zumino Witten models with k = 2s. Thus these simple wavefunctions reproduce the spin correlations of the family of Hamiltonians obtained by the Algebraic Bethe Ansatz.