How to sketch the Bode diagram of the output filter of a Buck converter?
Background
Starting with your result and proceeding:
$$\begin{align*}
H\left(s\right)&amp;=\frac{R}{R\,L\,C\,s^2+L\,s+R}\\\\
&amp;=\frac{1}{L\,C\,s^2+\frac LR\,s+1}\\\\
&amp;=\frac{\frac 1 {L\,C}}{s^2+\frac 1{R\,C}s +\frac 1{L\,C}}\tag{1}
\end{align*}$$
The denominator you see in equation (1) is sometimes called the characteristic equation. There are two relatively obvious time constants present there, \$\tau_{_0}=\sqrt{L\,C}\$ and \$\tau_{_1}={R\,C}\$, and therefore also two obvious angular frequencies, \$\omega_{_0}=\frac 1{\sqrt{L\,C}}\$ and \$\omega_{_1}=\frac 1{R\,C}\$.
So this would suggest the following:
$$\begin{align*}
H\left(s\right)
&amp;=\frac{\omega_{_0}^2}{s^2+\omega_{_1}\,s +\omega_{_0}^2}\tag{2}
\end{align*}$$
Standardized Form
Equation (2) isn't in standardized form, yet. And there's a good reason to do some more work. Equation (2) has two different angular frequencies and it's not clear how they might interact. Perhaps there's another way to express it that will help clarify things. Some kind of &quot;standard form&quot; that will be more readily meaningful?
The most interesting behaviors for equation (2) will be nearer where the denominator is zero and this is where focused attention would make sense. So let's solve the quadratic equation, \$\frac{-b\pm\sqrt{b^2-4\:a\:c}}{2\:a}\$. A quick glance suggests setting \$\alpha=\frac 12\omega_{_1}\$. Then we have \$s_1=-\alpha+\sqrt{\alpha^2-\omega_{_0}^2}\$ and \$s_2=-\alpha-\sqrt{\alpha^2-\omega_{_0}^2}\$ or else, written another way, \$s_1=-\alpha+j\sqrt{\omega_{_0}^2-\alpha^2}\$ and \$s_2=-\alpha-j\sqrt{\omega_{_0}^2-\alpha^2}\$. If we now set the damped angular frequency as \$\omega_d=\sqrt{\omega_{_0}^2-\alpha^2}\$ then \$s_1=-\alpha+j\,\omega_d\$ and \$s_2=-\alpha-j\,\omega_d\$.
If \$\alpha=\omega_{_0}\$ then \$\omega_d=0\$ and the system is critically-damped, with \$s_1=s_2=-\alpha\$, and both obviously the same real value. If \$\alpha\gt\omega_{_0}\$ then \$\omega_d\$ is imaginary and the system is over-damped, with \$s_1\ne s_2\$ but both real-valued. If \$\alpha\lt\omega_{_0}\$ then \$\omega_d\$ is real and the system is under-damped, \$s_1\$ and \$s_2\$ being complex conjugates.
Before proceeding further, it's worth looking again at \$\omega_d=\sqrt{\omega_{_0}^2-\alpha^2}\$. This can be expressed instead as \$\omega_d=\omega_{_0}\sqrt{1-\left[\frac{\alpha}{\omega_{_0}}\right]^2}\$. If we now set \$\zeta=\frac{\alpha}{\omega_{_0}}\$ (a dimensionless value) then \$\omega_d=\omega_{_0}\sqrt{1-\zeta^2}\$. It also then follows that \$\alpha=\zeta\,\omega_{_0}\$.
To zero in on the poles, the interesting areas surrounding the location where the denominator is zero, let's reformulate equation (2) above:
$$\begin{align*}
H\left(s\right)
&amp;=\frac{\omega_{_0}^2}{\left(s-s_1\right)\cdot\left(s-s_2\right)}\\\\
&amp;=\frac{\omega_{_0}^2}{\left(s-\left[-\alpha+j\,\omega_d\right]\right)\cdot\left(s-\left[-\alpha-j\,\omega_d\right]\right)}\\\\
&amp;=\frac{\omega_{_0}^2}{s^2+2\,\alpha\,s+\alpha^2+\omega_d^2}=\frac{\omega_{_0}^2}{\left(s+\alpha\right)^2+\omega_d^2}\\\\
&amp;=\frac{\omega_{_0}^2}{s^2+2\,\zeta\,\omega_{_0}\,s+\left(\zeta\,\omega_{_0}\right)^2+\omega_{_0}^2\left(1-\zeta^2\right)}\\\\
&amp;=\frac{\omega_{_0}^2}{s^2+2\,\zeta\,\omega_{_0}\,s+\omega_{_0}^2}\tag{3}
\end{align*}$$
This is the standardized form (except it is missing the gain, \$K\$) and something significant has been achieved here. We now have a single angular frequency, \$\omega_{_0}\$, and a new dimensionless damping factor, \$\zeta\$. If \$\zeta=1\$ then the system is critically damped. If \$\zeta \gt 1\$ then the system is over-damped. If \$\zeta\lt 1\$ then the system is under-damped.

Note 1: There is something else called a damping value (or in Sallen
and Key's writings it is called a dissipation factor), which has
units. It is designated by the letter \$d\$ and not \$\zeta\$. So be
cautious when reading. Earlier writing will often use the damping value
or dissipation factor, \$d\$. This will be particularly true for earlier
authors who are familiar with Sallen and Key's, &quot;A Practical Method of
Designing RC Active Filters&quot;, Technical Report No. 50 (now unclassified.)
Newer books have evolved to use the above approach and therefore have now
tended to standardize on using the unitless damping factor, \$\zeta\$, together with \$\omega_{_0}\$, instead.



Note 2: If you want to study the filter behavior without getting
bogged down by the angular frequency itself, just set
\$\omega_{_0}=1\:\frac{\text{rad}}{\text{s}}\$. The equation is a lot
simpler now, \$H\left(s\right)=\frac{1}{s^2+2\,\zeta\,s+1}\$ and
there's only one parameter to play with, \$\zeta\$. The amplitude is
then:
\$\frac{e_\text{out}}{e_\text{in}}=-20\operatorname{log}_{10}\left(\sqrt{\omega^4+\left(4\,\zeta^2-2\right)\omega^2+1}\right)\$
If you plug in \$\omega=1\:\frac{\text{rad}}{\text{s}}\$ you'll
find that the under-damped peak value is
\$\frac{e_\text{out}}{e_\text{in}}=-20\operatorname{log}_{10}\left(2\,\zeta\right)\$.
The phase response is also now just
\$\phi=\operatorname{tan}^{-1}\left(\frac{2\,\zeta\,\omega}{1-\omega^2}\right)\$.
So these two things make it very easy to study the amplitude and phase
response in this special case. But once you know the case where
\$\omega_{_0}=1\:\frac{\text{rad}}{\text{s}}\$, you have already
studied every possible other value for \$\omega_{_0}\$. All you have
to do is &quot;re-normalize&quot; things to different values of \$\omega_{_0}\$.
What you learned from \$\omega_{_0}=1\:\frac{\text{rad}}{\text{s}}\$
covers every other case. So just keep it simple while studying.

Your Question
In your case, the unitless damping factor is \$\zeta=\frac \alpha{\omega_0}\approx 0.158\$ and this means it's under-damped. From this, you can work out the peaking to be \$-20\operatorname{log}_{10}\left(2\zeta\right)\approx +10\:\text{dB}\$. You also know this is of \$2^\text{nd}\$ order, so the behavior at angular frequencies far higher than \$\omega_{_0}\approx 31.623\:\text{k}\frac{\text{rad}}{\text{s}}\$ drops off at \$-40\:\frac{\text{dB}}{\text{decade}}\$. At angular frequencies far less than \$\omega_{_0}\$, the magnitude will be \$0\:\text{dB}\$. And in general, far less means \$\lt \frac1{10}\,\omega_{_0}\$ and far more means \$\gt 10\,\omega_{_0}\$. Right at \$\omega_{_0}\$, the peak will be at \$+10\:\text{dB}\$ (as already mentioned.)
Given this is \$2^\text{nd}\$ order, you'll expect to see the phase flat and close to \$0^\circ\$ at \$\frac1{10}\,\omega_{_0}\$, curving rapidly towards \$-90^\circ\$ at approximately \$\omega_{_0}\$ (where the inflection point is at), and then curving back towards flat again at \$-180^\circ\$ at \$10\,\omega_{_0}\$.
I've hand-drawn a starting blue line to show the basic case without peaking. Just straight lines to the corner frequency (about \$5\:\text{kHz}\$.) The green bit is a hand-drawn bit of peaking to the computed peak value. The red line started out as three separate straight lines: one on the left to about \$500\:\text{Hz}\$ at \$0^\circ\$, another on the right going from about \$50\:\text{kHz}\$ and higher, and a third showing a rapid dive (almost but not quite vertical) going through \$-90^\circ\$ at \$5\:\text{kHz}\$. The rest was sketched in to connect those up. All such \$2^\text{nd}\$ order plots look the same. It's just positioning with values that makes one different from another. (The sharpness of the phase change is greater with higher peaking. Lower peaking will have a more gradual phase change. I'll leave the details for you to work out over time.)

None of this addresses your buck converter as a system. But it covers your output filter.