  In this paper, we obtain optimal time-decay rates in $L^r(\mathbb R^3_+)$ for
$r\ge 1$ of global strong solutions to the nematic liquid crystal flows in
$\mathbb R^3_+$, provided the initial data has small $L^3(\mathbb R^3_+)$-norm.
