The phase diagram of two-dimensional continuous particle systems is studied using the event-chain Monte Carlo algorithm. For soft disks with repulsive power-law interactions ∝r^{-n} with n≳6, the recently established hard-disk melting scenario (n→∞) holds: a first-order liquid-hexatic and a continuous hexatic-solid transition are identified. Close to n=6, the coexisting liquid exhibits very long orientational correlations, and positional correlations in the hexatic are extremely short. For n≲6, the liquid-hexatic transition is continuous, with correlations consistent with the Kosterlitz-Thouless-Halperin-Nelson-Young (KTHNY) scenario. To illustrate the generality of these results, we demonstrate that Yukawa particles likewise may follow either the KTHNY or the hard-disk melting scenario, depending on the Debye-Hückel screening length as well as on the temperature.

Local structure characterization with the bond-orientational order parameters q 4 , q 6 ,. . . introduced by Steinhardt et al. has become a standard tool in condensed matter physics, with applications including glass, jamming, melting or crystallization transitions and cluster formation. Here we discuss two fundamental flaws in the definition of these parameters that significantly affect their interpretation for studies of disordered systems, and offer a remedy. First, the definition of the bond-orientational order parameters considers the geometrical arrangement of a set of neighboring spheres NN(p) around a given central particle p; we show that procedure to select the spheres constituting the neighborhood NN(p) can have greater influence on both the numerical values and qualitative trend of q l than a change of the physical parameters, such as packing fraction. Second, the discrete nature of neighborhood implies that NN(p) is not a continuous function of the particle coordinates; this discontinuity, inherited by q l , leads to a lack of robustness of the q l as structure metrics. Both issues can be avoided by a morphometric approach leading to the robust Minkowski structure metrics q ′ l . These q ′ l are of a similar mathematical form as the conventional bond-orientational order parameters and are mathematically equivalent to the recently introduced Minkowski tensors [Europhys. Lett. 90, 34001 (2010); Phys. Rev. E. 85, 030301 (2012)].