Ultracold neutral plasmas, formed by photoionizing laser-cooled atoms near the ionization threshold, have electron temperatures in the 1-1000 kelvin range and ion temperatures from tens of millikelvin to a few kelvin. They represent a new frontier in the study of neutral plasmas, which traditionally deals with much hotter systems, but they also blur the boundaries of plasma, atomic, condensed matter, and low temperature physics. Modelling these plasmas challenges computational techniques and theories of non-equilibrium systems, so the field has attracted great interest from the theoretical and computational physics communities. By varying laser intensities and wavelengths it is possible to accurately set the initial plasma density and energy, and charged-particle-detection and optical diagnostics allow precise measurements for comparison with theoretical predictions. Recent experiments using optical probes demonstrated that ions in the plasma equilibrate in a strongly coupled fluid phase. Strongly coupled plasmas, in which the electrical interaction energy between charged particles exceeds the average kinetic energy, reverse the traditional energy hierarchy underlying basic plasma concepts such as Debye screening and hydrodynamics. Equilibration in this regime is of particular interest because it involves the establishment of spatial correlations between particles, and it connects to the physics of the interiors of gas-giant planets and inertial confinement fusion devices.