Dynamical properties of Néel and valence-bond phases in the J -J model on the honeycomb lattice.

By using a variational Monte Carlo technique based upon Gutzwiller-projected fermionic states, we investigate the dynamical structure factor of the antiferromagnetic S = 1/2 Heisenberg model on the honeycomb lattice, in presence of first-neighbor (J ) and second-neighbor (J ) couplings, for J < 0.5J . The ground state of the system shows long-range antiferromagnetic order for J /J ≲ 0.23 (Néel phase), plaquette valence-bond order for 0.23 ≲ J /J ≲ 0.36, and columnar dimer order for J /J ≳ 0.36. Within the Néel phase, a well-defined magnon mode is observed, whose dispersion is in relatively good agreement with linear spin-wave approximation for J = 0. When a nonzero second-neighbor super-exchange is included, a roton-like mode develops around the K point (i.e., the corner of the Brillouin zone). This mode softens when J /J is increased and becomes gapless at the transition point, J /J ≈ 0.23. Here, a broad continuum of states is clearly visible in the dynamical spectrum, suggesting that nearly-deconfined spinon excitations could exist, at least at relatively high energies. For larger values of J /J , valence-bond order is detected and the spectrum of the system becomes clearly gapped, with a triplon mode at low energies. This is particularly evident for the spectrum of the dimer valence-bond phase, in which the triplon mode is rather well separated from the continuum of excitations that appears at higher energies.