We show that a class of exactly solvable quantum Ising models, including the transverse-field Ising model and anisotropic XY model, can be characterized as the loops in a two-dimensional auxiliary space. The transverse-field Ising model corresponds to a circle and the XY model corresponds to an ellipse, while other models yield cardioid, limacon, hypocycloid, and Lissajous curves etc. It is shown that the variation of the ground state energy density, which is a function of the loop, experiences a nonanalytical point when the winding number of the corresponding loop changes. The winding number can serve as a topological quantum number of the quantum phases in the extended quantum Ising model, which sheds some light upon the relation between quantum phase transition and the geometrical order parameter characterizing the phase diagram.PACS numbers: 03.65. Vf, 75.10.Jm, 05.70.Fh, Introduction. Characterizing the quantum phase transitions (QPTs) is of central significance to both condensed matter physics and quantum information science. QPTs occur only at zero temperature due to the competition between different parameters describing the interactions of the system. A quantitative understanding of the second-order QPT is that the ground state undergoes qualitative changes when an external parameter passes through quantum critical points (QCPs).There are two prototypical models, Bose-Hubbard model and transverse-field Ising model, based on which the concept and characteristic of QPTs can be well demonstrated. However, among the two, only the transverse-field Ising model is exactly solvable [1], so as to be a unique paradigm for understanding the QPTs. Recently, more attention has been paid to theoretical studies of exactly solvable quantum spin models involving nearest-, next-nearest-neighbor interactions, and multiple spin exchange models, etc [3][4][5][6][7][8][9][10]. Those models are closer to real quasi-one-dimensional magnets [11][12][13] comparing to standard ones with only nearest-neighbor couplings. Furthermore, it has been shown that quantum spin models can be simulated in artificial quantum system with controllable parameters. Quantum simulation of spin chain can be experimentally realized through neutral atoms stored in an optical lattice [14,15], trapped ions [16][17][18][19][20][21][22][23][24] and NMR simulator [25]. This system often serves as a test bed for applying new ideas and methods to quantum phase transitions.A fundamental question is whether QPTs in Ising model can have a connection to some topological characterizations. It is interesting to note in this context that some simple Ising models have been found to exhibit topological characterization [26][27][28][29]. The purpose of the present work is to shed some light upon the relation between QPTs and a geometrical parameter characterizing the phase diagram, through the investigation of a class of quantum Ising models.In this work, we present an extended quantum Ising