We study the gradient flow, generated by the Landau-De Gennes energy functional, in the physically relevant spatial dimensions d= 2 , 3. We establish global well-posedness and global exponential time decay bounds for large H1 data in the 2D case, and uniform bounds for large data in 3D. This is indeed the best possible outcome for unrestricted coefficients in 3D, given that steady states do exist, at least for some coefficient configurations. We also establish leading order terms and in particular sharp asymptotics for the said dynamics in 2D. In 3D, we similarly isolate the leading order term, under the necessary assumption that a given, possibly large, solution converges to zero as t→ ∞. As a corollary, we prove an asymptotic formula for the correlation functional, c(y,t)=∫Rdtr(Q(x+y,t)Q(x,t))dx∫Rdtr(Q2(x,t))dx=e-|y|28t+OLy∞(t-12),d=2,3 for (potentially large) solutions Q(t) obeying a natural asymptotic condition. Such a formula was established in Kirr (J Stat Phys 155:625–657 (2014) for d= 3 and small initial data Q, subject to the non-degeneracy condition ∫R3Q0(x)dx≠0.