We consider here Musielak-Orlicz Sobolev (MOS) spaces $W^{k,\Phi}(\Omega)$, where $\Omega$ is an open subset of $\mathbb{R}^d$, $k\in\mathbb{N}$ and $\Phi$ is a Musielak-Orlicz function. The main outcomes consist of the results on density of the space of compactly supported smooth functions $C_C^\infty(\Omega)$ in $W^{k,\Phi}(\Omega)$. One section is devoted to compare the various conditions on $\Phi$ appearing in the literature in the context of maximal operator and density theorems in MOS spaces. The assumptions on $\Phi$ we apply here are substantially weaker than in the earlier papers on the topics of approximation by smooth functions \cite{Ahmida, Has1, Hasbook}. One of the reasons is that in the process of proving density theorems we do not use the fact that the Hardy-Littlewood maximal operator on Musielak-Orlicz space $L^\Phi(\Omega)$ is bounded, a standard tool employed in density results for different types of Sobolev spaces. We show in particular that under some regularity assumptions on $\Phi$, (A1) and $\Delta_2$ that are not sufficient for the maximal operator to be bounded, the space of $C_C^\infty(\mathbb{R}^d)$ is dense in $W^{k,\Phi}(\Omega)$. In the case of variable exponent Sobolev space $W^{k,p(\cdot)}(\mathbb{R}^d)$, we obtain the similar density result under assumption that $\Phi(x,t) = t^{p(x)}$, $p(x)\ge 1$, $t\ge 0$, $x\in \mathbb{R}^d$, satisfies the log-H\"older condition and the exponent function $p$ is essentially bounded.
Comment: new version will be posted