  We use N-body chemo-dynamic simulations to study the coupling between
morphology, kinematics and metallicity of the bar/bulge region of our Galaxy.
We make qualitative comparisons of our results with available observations and
find very good agreement. We conclude that this region is complex, since it
comprises several stellar components with different properties -- i.e. a
boxy/peanut bulge, thin and thick disc components, and, to lesser extents, a
disky pseudobulge, a stellar halo and a small classical bulge -- all cohabiting
in dynamical equilibrium. Our models show strong links between kinematics and
metallicity, or morphology and metallicity, as already suggested by a number of
recent observations. We discuss and explain these links.
