We present results of dynamical simulations of N-f = 2 degenerate Wilson twisted mass quarks at maximal twist in the range of pseudo scalar masses 300 MeV less than or similar to mps; 550 MeV less than or similar to Reaching such small masses was made possible owing to a recently developed variant of the HNIC algorithm. The simulations are performed at one value of the lattice spacing a less than or similar to 0. 1 fun. In order to have O(a) improvement and aiming at small residual O(a(2)) cutoff effects, the theory is tuned to maximal twist by requiring the vanishing of the untwisted quark mass. Precise results for the pseudo scalar decay constant and the pseudo scalar mass are confronted with chiral perturbation theory predictions and the low energy constants F, (l) over bar (3) and (l) over bar (4) are evaluated with small statistical errors. (c) 2007 Elsevier B.V. All rights reserved,