AbstractThe spectral transition from Very Low Mass stars (VLMs) to brown dwarfs (BDs) and planetary mass objects (Planemos) requires model atmospheres that can treat line, molecule, and dust-cloud formation with completeness and accuracy. One of the essential problems is the determination of the surface velocity field throughout the main sequence down to the BD and planemo mass regimes. We present local 2D and 3D radiation hydrodynamic simulations using the CO5BOLD code with binned Phoenix gas opacities, forsterite dust formation (and opacities) and rotation. The resulting velocity field vs depth and Teff has been used in the general purpose model atmosphere code Phoenix, adapted in static 1D spherical symmetry for these cool atmospheres. The result is a better understanding of the spectral transition from the stellar to substellar regimes. However, problems remain in reproducing the colors of the dustiest brown dwarfs. The global properties of rotation can change the averaged spectral properties of these objects. Our project for the period 2011-2015 is therefore to develop scaled down global 3D simulations of convection, cloud formation and rotation thanks to funding by the Agence Nationale de la Recherche in France.