In this paper, we consider the Beris-Edwards system for incompressible nematic liquid crystal flows. The system under investigation consists of the Navier-Stokes equations for the fluid velocity u coupled with an evolution equation for the order parameter Q-tensor. One important feature of the system is that its elastic free energy takes a general form and in particular, it contains a cubic term that possibly makes it unbounded from below. In the two dimensional periodic setting, we prove that if the initial L ∞ -norm of the Q-tensor is properly small, then the system admits a unique global weak solution. The proof is based on the construction of a specific approximating system that preserves the L ∞ -norm of the Q-tensor along the time evolution.Here, a, b, c ∈ R are material-dependent and temperature-dependent constants. On the other hand, the elastic free energy characterizes the distortion effect of the liquid crystal and its density function F elastic gives the strain energy density due to spatial variations in the Q-tensor:( 1.4) for 1 i, j, k, l d. Throughout this paper, we use the Einstein summation convention over repeated indices. The coefficients L 1 , L 2 , L 3 , L 4 are material-dependent elastic constants. We note that F elastic consists of three independent terms associated with L 1 , L 2 , L 3 that are quadratic in the first order partial derivatives of the Q-tensor, plus a cubic term associated with the coefficient L 4 . In the literature, the case L 2 = L 3 = L 4 = 0 is usually called isotropic, otherwise anisotropic if at least one of L 2 , L 3 , L 4 does not vanish. In particular, the retention of the cubic term (i.e., L 4 = 0) is due to the physically relevant consideration that it allows a complete reduction of the Landau-de Gennes energy E(Q) to the classical Oseen-Frank energy for nematic liquid crystals [6,31] (see also [24, Appendix B]). There exists a vast recent literature on the mathematical study of the Landau-de Gennes theory, and we refer interested readers to [2-7, 9-15, 17-24, 30-32, 34, 35, 38-40] as well as the references cited therein.In this paper, we study a basic model for the evolution of an incompressible nematic liquid crystal flow, which was first proposed by Beris-Edwards [7]. The resulting PDE system consists of the Navier-Stokes equations for the fluid velocity u and nonlinear convection-diffusion equations of parabolic type for the Q-tensor (see, e.g., [2,10,34]). To simplify the mathematical setting, throughout this paper we confine ourselves to the two dimensional periodic case such that d = 2 and Ω = T 2 , where T 2 stands for the periodic box with period ℓ i in the i-th direction with O = (0, ℓ 1 ) × (0, ℓ 2 ) being the periodic cell. Without loss of generality, we simply set O = (0, 1) 2 . Then the coupled system we are going to study takes the following form: