We present a general construction for exact analytic Taylor states in axisymmetric toroidal geometries. In this construction, the Taylor equilibria are fully determined by specifying the aspect ratio, elongation, and triangularity of the desired plasma geometry. For equilibria with a magnetic X-point, the location of the X-point must also be specified. The flexibility and simplicity of these solutions make them useful for verifying the accuracy of numerical solvers and for theoretical studies of Taylor states in laboratory experiments.Plasmas in both astrophysical and laboratory settings have a strong tendency to relax to minimum energy states known as Taylor states or Woltjer-Taylor states [1][2][3][4][5][6][7][8][9][10][11][12] in which the magnetic fields are force-free fields given by the equationwhere λ is a global constant. A well-known analytic solution to equation (1) is often used for theoretical studies and to interpret experiments [5][6][7]13 . One of its main advantage is its simplicity, but it lacks the degrees of freedom necessary to describe the large variety of configurations observed in laboratory experiments. We present a new family of exact solutions to equation (1) and a general construction for the solutions that address this need. The new solutions, while still simple, have the flexibility to describe configurations within a wide range of aspect ratios, elongations, and triangularities. The plasma boundary can have a magnetic separatrix, if desired, and the location of the separatrix can be specified. The equilibria we describe in this article can thus be useful for a variety of applications, including the study of non-solenoidal current start-up in low aspect ratio toroidal devices 10 and plasma dynamics in spheromaks. They can also be used to verify the accuracy of numerical schemes developed to solve equation (1) in fusionrelevant geometries 14,15 . Efficient solvers for force-free magnetic fields have recently become particularly attractive as a building block in a promising formulation for three-dimensional equilibria in fusion devices 16,17 . The exact solutions we present in this article can in that sense be thought of as the equivalent of Solov'ev solutions used to benchmark Grad-Shafranov solvers which are designed to compute more general equilibria 18 . The ability to construct exact equilibria with magnetic X-points is very desirable, since X-points are usually a source of difficulty in both theoretical studies and in numerical solvers. Our construction of analytic solutions works as follows. We first turn equation (1) into its associated GradShafranov equation for the poloidal flux function ψ. We a) Electronic mail: cerfon@cims.nyu.edu b) Electronic mail: oneil@cims.nyu.edu then express the solution ψ as a finite sum of functions satisfying the Grad-Shafranov equation. Finally, in order to have the ψ contours conform with shaped plasmas relevant to laboratory experiments, we determine the free constants appearing in the finite sum of functions such that the edge of the plasma, g...