Measurements of Higher Order Flow Harmonics in<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mi>Au</mml:mi><mml:mo>+</mml:mo><mml:mi>Au</mml:mi></mml:math>Collisions at<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:msqrt><mml:msub><mml:mi>s</mml:mi><mml:mrow><mml:mi>N</mml:mi><mml:mi>N</mml:mi></mml:mrow></mml:msub></mml:msqrt><mml:mo>=</mml:mo><mml:mn>200</mml:mn><mml:mtext> </mml:mtext><mml:mtext> </mml:mtext><mml:mi>GeV</mml:mi></mml:math>

Flow coefficients ν(n) for n=2, 3, 4, characterizing the anisotropic collective flow in Au+Au collisions at √s(NN)=200 GeV, are measured relative to event planes Ψ(n), determined at large rapidity. We report ν(n) as a function of transverse momentum and collision centrality, and study the correlations among the event planes of different order n. The ν(n) are well described by hydrodynamic models which employ a Glauber Monte Carlo initial state geometry with fluctuations, providing additional constraining power on the interplay between initial conditions and the effects of viscosity as the system evolves. This new constraint can serve to improve the precision of the extracted shear viscosity to entropy density ratio η/s.