  The LHCb collaboration has found that the production rate of $X(3872)$ in
proton-proton collisions decreases as final state particle multiplicity
increases. Moreover, the ALICE experiment at CERN has observed that the number
of deuterons produced increases with multiplicity, a behavior that is
qualitatively different from that of the $X(3872)$. These experimental findings
may point to a compact structure for the $X(3872)$ or, at least, that its
hadronization could proceed through a charm-anticharm core. We have recently
used a diffusion Monte Carlo method to solve the many-body Schr\"odinger
equation that describes the $X(3872)$ as a $c \bar c q \bar q$ tetraquark
system with quantum numbers $I^G(J^{PC})=0^+(1^{++})$ and $1^-(1^{++})$.
According to our structural analysis, the quark--(anti-)quark correlations
resemble light-meson--heavy-meson molecules of type $\omega J/\psi$ and $\rho
J/\psi$, rather than the most extended $D\bar D^{\ast}$ interpretation. It was
argued that this fact may be the key to make compatible the molecular features
of the $X(3872)$ with its production observables. The same formalism allows us
to compute the first color excited $c \bar c q \bar q$ tetraquark state with
either $I^G(J^{PC})=0^+(1^{++})$ or $1^-(1^{++})$. A bound-state is found in
each channel, their masses are around 4.0 GeV which is an energy region where
many new exotic candidates have been collected by the Particle Data Group.
Concerning their structural properties, these states cluster in a compact
diquark-antidiquark arrangement which matches perfectly with a so-called
Born-Oppenheimer tetraquark configuration. The promptly production rates of
these states in proton-proton, proton-nucleus and nucleus-nucleus collisions
should fall off equal to or even faster than those of the $X(3872)$.
