The neutron-rich zirconium isotopes are a key testing ground for nuclear models due to their sensitivity to shape changes, and because they cross the r-process path of explosive nucleosynthesis. Furthermore, recent experimental data have revealed a high-spin isomer in 108 Zr. Here we report the results of configuration-constrained potential-energy-surface calculations of ground states and high-K states in [104][105][106][107][108][109][110] Zr, restricted to quadrupole and hexadecapole shapes. These calculations enable the existing isomer data to be understood, and predictions are made for 110 Zr. The zirconium (Z = 40) isotopes are well known [1] for their variety of shapes, with considerable sensitivity to the neutron number, which is at least partly a consequence of the Z = 40 subshell gap in single-particle energies. On the neutron-rich side of stability, the shape sensitivity extends into the landscape of rapid-neutron-capture, r-process nucleosynthesis relevant to supernova detonations, giving an added incentive to understand the nuclear structure. In the past decade there have been many theoretical studies related to this theme (see for example Refs. [2-9]), but only Xu et al. [2], albeit with a more phenomenological approach, have given detailed predictions of isomeric states in the neutron-rich even-even zirconium isotopes. Now, with the success of the RIKEN RIBF accelerator complex, major experimental advances have been possible, extending knowledge of these isotopes [10,11]. One of the notable nuclear structure features was the discovery [11] of an isomer in 108 Zr, but not in 106 Zr or 104 Zr. The isomer in 108 Zr was tentatively interpreted as being due a tetrahedral shape, which has been predicted to be most favored in 110 Zr [12], though the excited-state structure of that isotope has not yet itself been studied experimentally. In the present work, we extend the earlier configuration-constrained potential-energysurface (PES) calculations for this region [2], now focused on even-even [104][105][106][107][108][109][110] Zr. We consider the quadrupole (β 2 and γ ) and hexadecapole (β 4 ) deformations as the relevant degrees of freedom, and show that the existing data can be understood on this basis. We also predict the low-lying two-quasiparticle structure of 110 Zr. The configuration-constrained PES model [13] has been widely used in various mass regions, including the superheavy [14,15] and drip-line [16] regions, to calculate the energies and shapes of high-K states, where K is the projection of the angular momentum on the nuclear axis of symmetry. High-K states at low excitation energy are often associated with isomerism, due to the K-forbidden nature of their decays * p.walker@surrey.ac.uk [17]. In the model, single-particle levels are obtained from the Woods-Saxon potential with the set of universal parameters [18]. For the pairing correlations, particle-number projection is approximated by the Lipkin-Nogami technique [19], with the pairing strength determined by the average-gap method [20]. ...