“…The knot matrix model is discussed through Harer-Zagier (HZ) formula defined below [50,51], which becomes a generating function of Gaussian mean. The knot polynomial of HOMFLYPT for a torus knot is realized as an average of a character [51]. For instance, the trefoil T (2, 3), which is a torus knot of (2,3) type, and the normalized HOMFLYPT polynomial P , by the multiplication of (v − v −1 )/z with the replacement of v → q N and z → q − q −1 , is expressed as P (T (2, 3) : v, z) = (q 2 + q −2 )q 5N − (q 2 + 1 + q −2 )q 3N + q N q − q −1 (9.5)…”