“…where Ω ⊂ R 3 , u ∈ W 1,p (Ω, R 3 ) for some p > 3, n ∈ H 1 (u(Ω), R 2 ), and W mec (F, n) = W (α −1 n ⊗ n + √ α(I − n ⊗ n))F (1.2) for a certain α > 0 and some polyconvex energy function W . Functionals with a similar structure appear also in models describing the nematic mesogens with the Landau-de Gennes theory, and in magnetoelasticity and plasticity, see, e.g., [6,12,18,28,5]. The major difficulties are that I depends on the composition of the two unknowns and that the nematic director n is defined in the domain u(Ω) which is also determined only as a part of the solution of the variational problem.…”