Dynamical properties of Néel and valence-bond phases in the J1-J2 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 (J1) and second-neighbor (J2) couplings, for J2< 0.5 J1. The ground state of the system shows long-range antiferromagnetic order for J2/J1< 0.23 (Néel phase), plaquette valence-bond order for 0.23 < J2/J1< 0.36, and columnar dimer order for J2/J1> 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 J2=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 J2/J1is increased and becomes gapless at the transition point, J2/J1≈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 J2/J1, 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.