How can the 1 dimensional Maxwell-Boltzmann Distribution be derived? In terms of energy, the 3 dimensional MB Distribution is giving the probability for a particle to have an energy $E \geq E + dE$ is:
$$f(E) = \frac{2}{\sqrt \pi} \cdot \bigg(\frac{1}{k_BT}\bigg)^{\frac{3}{2}} \cdot e^{-\frac{E}{k_BT}} \cdot \sqrt{E} \cdot dE$$
It is said that the 1 dimensional MB Distribution, giving the probability for a particle to have a certain energy in 1 degree of freedom, say $E_x \geq E_x + dE_x$, is:
$$f(E_x) = \sqrt{\frac{1}{\pi E_x k_B T}} \cdot e^{-\frac{E_x}{k_BT}} \cdot dE_x$$
How is this 1D MB-Distribution derived from scratch and from the 3D MB Distribution?
 A: The MB distribution for the one degree of freedom in terms of momentum is:
$$
w(p_x) dp_x = \frac1{\sqrt{2\pi m k_BT}} e^{-\frac{p_x^2}{2mk_BT}} dp_x
$$
The relation between $p_x$ and $E$ is
$$
E = \frac{p_x^2}{2m}\ \longleftrightarrow \ p_x(E) = \pm\sqrt{2mE}.
$$
Probability "conservation" requires
$$
f(E) dE = 2 w(p_x(E)) dp_x(E),
$$
where $p_x = \sqrt{2mE}$ is chosen and $2$ factor accounts for two possibilities. From the last equation follows
$$
f(E) = 2\ w(p_x(E))\ |p_x'(E)|.
$$
This formula leads to the needed expression.
A: Consider a canonical ensemble of a single particle in a 1D container of length $L$. According to the (classical) Boltzmann probability distribution, the probability density function for the state of the system being the microstate with position $x'$ and momentum $p'$ is
$$ f_{x,p}(x',p') = \frac{e^{-\beta H(x', p')}}{\int\limits_0^L dx'' \int\limits_{-\infty}^\infty dp'' e^{-\beta H(x'',p'')}}$$ 
where $ H(x, p) = p^2/2m $ is the energy of the particle with position $x$ and momentum $p$, and $\beta = 1/k_BT$. The marginal probablilty density function for the momentum is
$$ f_p(p') = \int\limits_0^L dx' f_{x,p}(x',p') = \frac{\int\limits_0^L dx'e^{-\beta p'^2/2m}}{\int\limits_0^L dx'' \int\limits_{-\infty}^\infty dp'' e^{-\beta p''^2/2m}} = \frac{L e^{-\beta p'^2/2m}}{L\int\limits_{-\infty}^\infty dp'' e^{-\beta p''^2/2m}} = \sqrt{\frac{\beta}{2\pi m}} e^{-\beta p'^2/2m}. $$
To calculate the probability distribution function for the energy, note that the probability that the energy is less than $E'$ is given by
$$ \int\limits_0^E' dE'' f_E(E'') = \int\limits_{-\sqrt{2mE'}}^\sqrt{2mE'} dp' f_p(p'') = 2\int\limits_0^\sqrt{2mE'} dp' f_p(p'') . $$
Differentiating both sides and using the fundamental theorem of calculus,
$$f_E(E') = 2 \sqrt{\frac{m}{2E'}} f_p(\sqrt{2mE'}) = \sqrt{\frac{2m}{E'}} \sqrt{\frac{\beta}{2\pi m}} e^{-\beta E'} = \sqrt{\frac{\beta}{\pi E'}} e^{-\beta E'}. $$
If you don't know where the Boltzmann probability distribution that I started with comes from, it is a standard result of statistical mechanics that should be found in any statistical mechanics textbook. If you want to ask about that, I think that probably deserves its own question.
If you wish to start with a particle free to move in a 3D container of volume $V$,
$$ f_{\vec{x},\vec{p}}(\vec{x}',\vec{p}') = \frac{e^{-\beta H(\vec{x}', \vec{p}')}}{\int\limits_V d^3\vec{x}'' \int\limits_{-\infty}^\infty d^3\vec{p}'' e^{-\beta H(\vec{x}'',\vec{p}'')}}$$
$$ f_\vec{p}(\vec{p}') = \int\limits_V d^3\vec{x}' f_{\vec{x},\vec{p}}(\vec{x}',\vec{p}') = \frac{\int\limits_V d^3\vec{x}'e^{-\beta p'^2/2m}}{\int\limits_V d^3\vec{x}'' \int\limits_{-\infty}^\infty d^3\vec{p}'' e^{-\beta p''^2/2m}} = \frac{V e^{-\beta p'^2/2m}}{V\int\limits_{-\infty}^\infty d^3\vec{p}'' e^{-\beta p''^2/2m}} = \frac{e^{-\beta p'^2/2m}}{\left( \int\limits_{-\infty}^\infty dp_x'' e^{-\beta p_x''^2/2m} \right) \left( \int\limits_{-\infty}^\infty dp_y'' e^{-\beta p_y''^2/2m} \right) \left( \int\limits_{-\infty}^\infty dp_z'' e^{-\beta p_z''^2/2m} \right)} = \left( \frac{\beta}{2\pi m} \right)^{3/2} e^{-\beta p'^2/2m} = \left( \frac{\beta}{2\pi m} \right)^{3/2} e^{-\beta (p_x'^2 + p_y'^2 + p_z'^2)^2/2m}. $$
This is the joint probability density function for all three components of the momentum. If you want the marginal probability density function for the $x$-component only,
$$ f_{p_x}(p_x') = \int\limits_{-\infty}^\infty dp_y' \int\limits_{-\infty}^\infty dp_z' f_\vec{p}(\vec{p}') = \left( \frac{\beta}{2\pi m} \right)^{3/2} e^{-\beta p_x'^2/2m} \int\limits_{-\infty}^\infty dp_y' e^{-\beta p_y'^2/2m} \int\limits_{-\infty}^\infty dp_z' e^{-\beta p_z'^2/2m} = \sqrt{\frac{\beta}{2\pi m}} e^{-\beta p_x'^2/2m} $$
so you end up with the same distribution as the 1D case. If you want the distribution of the energy due to the motion along the $x$ direction (i.e. $E_x = p_x^2/2m$), the math to derive it is identical to the 1D case above, and you end up with
$$f_{E_x}(E_x') = \sqrt{\frac{\beta}{\pi E_x'}} e^{-\beta E_x'}. $$
