  We study the first derivative estimates for solutions to Monge-Amp\`ere
equations in terms of modulus of continuity. As a result, we establish the
optimal global log-Lipschitz continuity for the gradient of solutions to the
Monge-Amp\`ere equation with natural boundary conditions.
