We study the transition between N\'eel and columnar valence-bond solid ordering in two-dimensional
S=3/2 square lattice quantum antiferromagnets with SO(3) symmetry. According to the deconfined criticality scenario, this transition can be direct and continuous like the well-studied
S=1/2 case. To study the global phase diagram, we work with four multi-spin couplings with full rotational symmetry, that are free of the sign-problem of quantum Monte Carlo. Exploring the phase diagram with quantum Monte Carlo simulations, we find that the phase transition between N\'eel and valence-bond solid is strongly first-order in the parts of the phase diagram that we have accessed.
We use quantum Monte Carlo methods to study the ground-state phase diagram of a
S=1/2 honeycomb lattice magnet in which a nearest-neighbor antiferromagnetic exchange
J (favoring N\'eel order) competes with two different multispin interaction terms: a six-spin interaction
Q3 that favors columnar valence-bond solid (VBS) order, and a four-spin interaction
Q2 that favors staggered VBS order. For ${Q}_{3}\ensuremath{\sim}{Q}_{2}\ensuremath{\gg}J$, we establish that the competition between the two different VBS orders stabilizes N\'eel order in a large swath of the phase diagram even when
J is the smallest energy scale in the Hamiltonian. When ${Q}_{3}\ensuremath{\gg}({Q}_{2},J)$ $\mathrm{[}{Q}_{2}\ensuremath{\gg}({Q}_{3},J)\text{]}$, this model exhibits at zero temperature phase transition from the N\'eel state to a columnar (staggered) VBS state. We establish that the N\'eel-columnar VBS transition is continuous for all values of
Q2, and that critical properties along the entire phase boundary are well characterized by critical exponents and amplitudes of the noncompact
CP1 (NCCP1) theory of deconfined criticality, similar to what is observed on a square lattice. However, a surprising threefold anisotropy of the phase of the VBS order parameter at criticality, whose presence was recently noted at the
Q2=0 deconfined critical point, is seen to persist all along this phase boundary. We use a classical analogy to explore this by studying the critical point of a three-dimensional
XY model with a fourfold anisotropy field which is known to be weakly irrelevant at the three-dimensional
XY critical point. In this case, we again find that the critical anisotropy appears to saturate to a nonzero value over the range of sizes accessible to our simulations.