Nonuniqueness of nematic liquid crystal flows in dimension three

For suitable initial and boundary data, we construct infinitely many weak solutions to the nematic liquid crystal flows in dimension three. These solutions are in the axisymmetric class with bounded energy and “backward bubbling” at a large time.