{"id":1897,"date":"2022-12-30T23:30:00","date_gmt":"2022-12-31T04:30:00","guid":{"rendered":"http:\/\/underthehood.blogwyrm.com\/?p=1897"},"modified":"2022-11-27T16:18:16","modified_gmt":"2022-11-27T21:18:16","slug":"kinetic-theory-3-exploring-the-maxwell-boltzmann-distribution","status":"publish","type":"post","link":"https:\/\/underthehood.blogwyrm.com\/?p=1897","title":{"rendered":"Kinetic Theory 3 &#8211; Exploring the Maxwell-Boltzmann Distribution"},"content":{"rendered":"\n<p>In this post, we explore some of the physical implications of the Maxwell-Boltzmann speed distribution<\/p>\n\n\n\n<p>\\[ f(v) = 4 \\pi \\left( \\frac{m}{2 \\pi k_B T} \\right)^{3\/2} v^2 \\exp \\left( &#8211; \\frac{m v^2}{2 k_B T} \\right) \\; \\]<\/p>\n\n\n\n<p>derived in the previous post.<\/p>\n\n\n\n<p>The first thing to note is that the decaying exponential favors lower speed (in the limit as $v \\rightarrow \\infty$, $f(v) = 0$) while the factor of $v^4$ favors higher ones ($f(v=0) = 0 $).\u00a0 As a result, we expect that there is a peak in the distribution somewhere between these two extremes.\u00a0 To verify this, we can plot $f(v)$ for the case where $m$ is the mass of a nitrogen gas molecule ($m_{N_2} \\approx 28 Da$, where $1 Da \\approx 1.66 \\times 10^{-27} \\, kg$) and the temperature is $300 K$, corresponding, roughly, to room temperature.<\/p>\n\n\n\n<figure class=\"wp-block-image size-large is-resized\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/underthehood.blogwyrm.com\/wp-content\/uploads\/2022\/11\/UTH_12Dec_N2_room_temperature.png\" alt=\"\" class=\"wp-image-1896\" width=\"656\" height=\"378\" srcset=\"https:\/\/underthehood.blogwyrm.com\/wp-content\/uploads\/2022\/11\/UTH_12Dec_N2_room_temperature.png 832w, https:\/\/underthehood.blogwyrm.com\/wp-content\/uploads\/2022\/11\/UTH_12Dec_N2_room_temperature-300x173.png 300w, https:\/\/underthehood.blogwyrm.com\/wp-content\/uploads\/2022\/11\/UTH_12Dec_N2_room_temperature-768x443.png 768w, https:\/\/underthehood.blogwyrm.com\/wp-content\/uploads\/2022\/11\/UTH_12Dec_N2_room_temperature-810x467.png 810w\" sizes=\"auto, (max-width: 656px) 100vw, 656px\" \/><\/figure>\n\n\n\n<p>There is a distinct peak, which can be estimated by eye, at approximately $420 \\, km\/s$.&nbsp; This value is of the order-of-magnitude of the RMS speed given by<\/p>\n\n\n\n<p>\\[ v_{RMS} = \\sqrt{ \\frac{3 k_B T }{m_{N_2}}}&nbsp; = 510.77 \\, km\/s \\; ,\\]<\/p>\n\n\n\n<p>but is substantially smaller by about $20 \\, %$.&nbsp; The exact value at the peak, which corresponds to the most probable value, comes from finding the maximum of the distribution in the usual way.&nbsp; First we differentiate the distribution with respect to $v$ to get<\/p>\n\n\n\n<p>\\[ \\frac{d f(v)}{d v} = 4 \\pi \\left( \\frac{m}{2 \\pi k_B T} \\right)^{3\/2} \\left( 2 v &#8211; \\frac{ m v^3}{k_B T} \\right) \\exp \\left( &#8211; \\frac{v^2}{2 k_B T} \\right) \\; . \\]<\/p>\n\n\n\n<p>Setting this expression to zero and solving for $v$ yields the most probable speed<\/p>\n\n\n\n<p>\\[ v_{mp} = \\sqrt{\\frac{2 k_B T}{m} } \\; .\\]<\/p>\n\n\n\n<p>Plugging in the numerical values used above gives $v_{mp} = 417.04 \\, km\/s$, which is consistent with eyeball estimate of $420 \\, km\/s$.<\/p>\n\n\n\n<p>The final speed of interest is the average speed defined by<\/p>\n\n\n\n<p>\\[ v_{ave} = \\int_0^{\\infty}\\, dv \\,&nbsp; v f(v) \\; .\\]<\/p>\n\n\n\n<p>Ordinarily, the odd moment of any Gaussian-like distribution would be zero (e.g., the average of any component of the velocity would be zero &#8211; see last post) but since the speed is confined to the interval $[0,\\infty)$, the integral yields a finite value.&nbsp; To get that moment, we start with<\/p>\n\n\n\n<p>\\[ J_1 = \\int_0^{\\infty} \\, v e^{-qv^2} dv = \\frac{1}{2} \\int_0^{\\infty} \\, d(v^2) e^{-q v^2} \\; .\\]<\/p>\n\n\n\n<p>Substituting&nbsp; $w = v^2$ yields a simple integral<\/p>\n\n\n\n<p>\\[ J_1 = \\frac{1}{2} \\int_0^{\\infty} \\, dw e^{-q w} = \\frac{1}{2q} \\; . \\]<\/p>\n\n\n\n<p>Higher order moments come by differentiation with<\/p>\n\n\n\n<p>\\[ J_3 = \\int_0^{\\infty} v^3 e^{-q v^2} dv = -\\frac{d}{dq} \\int_0^{\\infty} v e^{-q v^2} = -\\frac{d}{dq} \\frac{1}{2q} = \\frac{1}{2q^2} \\; ,&nbsp; \\]<\/p>\n\n\n\n<p>which, up to some constants, is the desired result.<\/p>\n\n\n\n<p>Using this result, we arrive at the expression<\/p>\n\n\n\n<p>\\[ v_{ave} = 4 \\pi \\left(\\frac{m}{2 \\pi k_B T}\\right)^{3\/2} \\frac{k_B T}{m} = \\sqrt{ \\frac{8 k_B T}{\\pi m}} \\; .\\]<\/p>\n\n\n\n<p>Plugging in the numerical values used above gives $v_{ave} = 470.58 \\, km\/s$, which falls between the most probable and the RMS speeds.\u00a0 We can annotate the graph with these three lines to draw out the distinctions<\/p>\n\n\n\n<figure class=\"wp-block-image size-large is-resized\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/underthehood.blogwyrm.com\/wp-content\/uploads\/2022\/11\/UTH_12Dec_N2_room_temperature_anotated-1024x482.png\" alt=\"\" class=\"wp-image-1895\" width=\"637\" height=\"300\" srcset=\"https:\/\/underthehood.blogwyrm.com\/wp-content\/uploads\/2022\/11\/UTH_12Dec_N2_room_temperature_anotated-1024x482.png 1024w, https:\/\/underthehood.blogwyrm.com\/wp-content\/uploads\/2022\/11\/UTH_12Dec_N2_room_temperature_anotated-300x141.png 300w, https:\/\/underthehood.blogwyrm.com\/wp-content\/uploads\/2022\/11\/UTH_12Dec_N2_room_temperature_anotated-768x361.png 768w, https:\/\/underthehood.blogwyrm.com\/wp-content\/uploads\/2022\/11\/UTH_12Dec_N2_room_temperature_anotated-810x381.png 810w, https:\/\/underthehood.blogwyrm.com\/wp-content\/uploads\/2022\/11\/UTH_12Dec_N2_room_temperature_anotated.png 1146w\" sizes=\"auto, (max-width: 637px) 100vw, 637px\" \/><\/figure>\n\n\n\n<p>and we can note that the average speed is about $12.8 \\, %$ higher than the most probable speed while the RMS speed is $22.5 \\, %$ higher.<\/p>\n\n\n\n<p>The final point to be explored is how the curve shifts as a function of temperature.\u00a0 Since all the speeds have the same functional form, differing only in the numerical coefficient, it is straightforward to see that each speed scales as $\\sqrt{T}$.\u00a0 However, the overall shape of the curve can be a little surprising, which the following plot illustrates by looking at broad range of temperatures.<\/p>\n\n\n\n<figure class=\"wp-block-image size-large\"><img loading=\"lazy\" decoding=\"async\" width=\"731\" height=\"558\" src=\"http:\/\/underthehood.blogwyrm.com\/wp-content\/uploads\/2022\/11\/UTH_12Dec_N2_varying_temperature.png\" alt=\"\" class=\"wp-image-1894\" srcset=\"https:\/\/underthehood.blogwyrm.com\/wp-content\/uploads\/2022\/11\/UTH_12Dec_N2_varying_temperature.png 731w, https:\/\/underthehood.blogwyrm.com\/wp-content\/uploads\/2022\/11\/UTH_12Dec_N2_varying_temperature-300x229.png 300w\" sizes=\"auto, (max-width: 731px) 100vw, 731px\" \/><\/figure>\n\n\n\n<p>As the temperature increases, the distribution tends to being more symmetric by shifting right, effectively eating into the long tail.&nbsp;<\/p>\n\n\n\n<p>Of course, the Maxwell-Boltzmann distribution is physically unrealizable as there is always a finite probability for having a speed equal or greater to the speed of light.&nbsp; But this is of little concern as very little of the distribution is found at these higher speeds.&nbsp; For example, the vast majority of the distribution is found below $2000 \\, km\/s$, which is $ &lt; 0.01 c$, even at $T = 1000 K$.&nbsp;<\/p>\n","protected":false},"excerpt":{"rendered":"<p>In this post, we explore some of the physical implications of the Maxwell-Boltzmann speed distribution \\[ f(v) = 4 \\pi \\left( \\frac{m}{2 \\pi k_B T} \\right)^{3\/2} v^2 \\exp \\left( &#8211;&#8230; <a class=\"read-more-button\" href=\"https:\/\/underthehood.blogwyrm.com\/?p=1897\">Read more &gt;<\/a><\/p>\n","protected":false},"author":2,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[1],"tags":[],"class_list":["post-1897","post","type-post","status-publish","format-standard","hentry","category-uncategorized"],"_links":{"self":[{"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=\/wp\/v2\/posts\/1897","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=1897"}],"version-history":[{"count":3,"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=\/wp\/v2\/posts\/1897\/revisions"}],"predecessor-version":[{"id":1900,"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=\/wp\/v2\/posts\/1897\/revisions\/1900"}],"wp:attachment":[{"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=1897"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=1897"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/underthehood.blogwyrm.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=1897"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}