“…One is the energy-minimization based approach (Cohen, Hardt, Kinderlehrer, Lin and Luskin 1987, Alouges 1997, Adler, Atherton, Emerson and MacLachlan 2015, Nochetto, Walker and Zhang 2017, Majumdar and Wang 2018, Gartland Jr, Palffy-Muhoray and Varga 1991, which is often numerically solved by Newton-type or quasi-Newton method. The other approach is to follow the gradient flow dynamics driven by the free energy corresponding to individual model of LCs (Fukuda, Stark, Yoneya and Yokoyama 2004, Ravnik and Žumer 2009, Canevari, Majumdar and Spicer 2017, Wang, Canevari and Majumdar 2019, MacDonald, Mackenzie and Ramage 2020. Various efficient numerical methods have been developed to solve the gradient flow equations, including energy stable numerical schemes such as convex splitting method (Elliott and Stuart 1993), invariant energy quadratization method (Yang 2016), scalar auxiliary variable method , etc.…”