Come calcolare la distribuzione di Maxwell-Boltzmann delle velocità molecolari: distribuzione di probabilità per le velocità delle molecole in un gas ideale all equilibrio termico, formula fondamentale della teoria cinetica dei gas e fisica statistica. Formula PDF: f(v) = 4 pi (m/(2 pi k T))^(3/2) v^2 exp(-m v^2 / (2 k T)) dove v velocità molecolare (m/s), m massa molecola (kg), k costante Boltzmann 1.380649e-23 J/K, T temperatura (Kelvin).

Tre velocità caratteristiche: v_p most probable = sqrt(2 k T / m) (massimo PDF), v_mean = sqrt(8 k T / (pi m)) (media aritmetica), v_rms = sqrt(3 k T / m) (radice quadrata media energia cinetica). Sempre vale: v_p < v_mean < v_rms (rapporti fissi v_p:v_mean:v_rms = 1 : 1.128 : 1.225). Energia cinetica media per molecola = (3/2) k T; per mole = (3/2) R T = 3.72 kJ/mol a 298 K. Storia: James Clerk Maxwell (Scozia 1831-1879) Cambridge Trinity College 1860 paper Philosophical Magazine 19:19 'Illustrations of the Dynamical Theory of Gases' derivazione original molecular velocity distribution argumento simmetria probabilità indipendente direzioni; Ludwig Boltzmann (Austria 1844-1906) Vienna Graz Leipzig München 1872 paper Wiener Berichte 66:275 'Weitere Studien über das Wärmegleichgewicht unter Gasmolekülen' H-theorem derivazione rigorosa kinetic theory entropy + 1877 paper Wiener Berichte 76:373 entropy = k log W relation engraved tomba Vienna Zentralfriedhof.

Strumento essenziale per: fisica statistica Italia ricerca Sapienza Roma DIET Mathematica Statistica; Bologna alma mater 1088 prima Europa Faculty Physics; Federico II Napoli DIETI Computer Science; Padova DEI Galileo professor 1592-1610 Bioengineering; Pisa SNS Scuola Normale Superiore 1810 Mathematics Physics; Politecnico Milano DEIB Bioengineering; Politecnico Torino DET; Pavia 1361.

Italia INFN Istituto Nazionale Fisica Nucleare Frascati 1951 fondato 5000 ricercatori 28 sezioni Nuclear Physics; INFN Cosenza Catania Lecce sezioni meridionali; CNR Consiglio Nazionale Ricerche Roma 1923 fondato 8000 ricercatori 200 institutes; ENEA Roma energia ambiente; ASI Agenzia Spaziale Italiana Roma 1988; INAF Istituto Nazionale Astrofisica Roma 24 osservatori.

Italia fisici nobel storia: Guglielmo Marconi (Italia-Irlanda 1874-1937) Bologna Liverpool Premio Nobel Fisica 1909 telegrafia senza fili; Enrico Fermi (Italia-USA 1901-1954) Sapienza Roma 1922-1938 Nobel Fisica 1938 reazioni nucleari neutroni lenti, emigrato USA leggi razziali fascismo, Manhattan Project Los Alamos NM atomic bomb 1945, Chicago Pile-1 prima reazione nucleare controllata 1942; Emilio Segre (Italia-USA 1905-1989) Sapienza Roma + Berkeley Nobel Fisica 1959 antiprotone; Carlo Rubbia (Italia 1934-) CERN Premio Nobel Fisica 1984 W Z bosoni unificazione elettrodeboli; Riccardo Giacconi (Italia-USA 1931-2018) Premio Nobel Fisica 2002 raggi X astrofisica; Giorgio Parisi (Italia 1948-) Sapienza Roma Premio Nobel Fisica 2021 sistemi complessi spin glass disorder.

Italia industria gas e high-tech: SAES Getters Lainate 1940 getter materiali ultra-vuoto 1500 employees worldwide; Eni Milano 1953 petrolchimica gas naturale 32000 employees; Snam Milano 1941 gas pipeline; Italgas Milano 2017 gas distribution; Edison Milano 1884 energia gas; Saipem Milano oil&gas. Sangue umano composizione gas: O2 21% atmospheric, N2 78%, CO2 0.04%; respirazione capacità polmonare 6L; arterioso pO2 100 mmHg, venoso 40 mmHg; PaCO2 arterioso 35-45 mmHg; emoglobina P50 26 mmHg O2.

Formula Maxwell-Boltzmann

PDF velocità scalare: f(v) = 4 pi (m/(2 pi k T))^(3/2) v^2 exp(-m v^2 / (2 k T)) per v >= 0. Distribuzione modulo velocità (speed) ottenuta integrando velocity distribution f(vx, vy, vz) coordinate cartesiane su sfera 4 pi v^2 dv (Jacobian sferiche). m massa molecola singola (kg) = M / N_A dove M massa molare (kg/mol), N_A Avogadro 6.022e23. Costante Boltzmann k = 1.380649e-23 J/K (definita esattamente dal 2019 SI revision).

Tre velocità caratteristiche: v_p most probable = sqrt(2 k T / m) massimo PDF (df/dv=0); v_mean = sqrt(8 k T / (pi m)) integral v f(v) dv; v_rms = sqrt(3 k T / m) sqrt(integral v^2 f(v) dv). Rapporti universali v_p:v_mean:v_rms = 1 : sqrt(4/pi) : sqrt(3/2) = 1 : 1.1284 : 1.2247 indipendenti da T e m.

Energia cinetica: E_k = (1/2) m v_rms^2 = (3/2) k T per molecola. Per mole E_k = (3/2) R T dove R = 8.314 J/(mol K) costante gas. A T=298 K (25°C ambiente): E_k = 3.72 kJ/mol = 0.0257 eV per molecola. Equipartizione 3 gradi libertà traslazionali ciascuno (1/2) k T.

CDF cumulativa: F(v) = integral_0^v f(v') dv' = erf(v / v_p) - (2/sqrt(pi)) (v/v_p) exp(-(v/v_p)^2) dove erf error function. Frazione molecole con velocità v = 1 - F(v).

Effetto temperatura: aumento T sposta distribuzione verso velocità più alte e amplia spread. v_p, v_mean, v_rms scalano con sqrt(T). Distribuzione bassa T picco stretto velocità basse; alta T picco largo velocità alte. Esempi numerici N2 aria: T=200K v_rms=421 m/s; T=298K v_rms=515 m/s; T=1000K v_rms=943 m/s.

Storia e contesto

James Clerk Maxwell (Scozia 1831-1879) Edinburgh nato, Cambridge Trinity College mathematics natural philosophy 1850-1854 graduated Senior Wrangler Smith's Prize, King's College London 1860-1865, Cambridge Trinity Cavendish Laboratory founded 1871 first director. Maxwell 1860 papers Philosophical Magazine 19:19 'Illustrations of the Dynamical Theory of Gases Part I' + 19:21 Part II: derivazione originale Maxwell-Boltzmann distribution argumento simmetria probabilità velocità componenti indipendenti distributi normalmente, prodotto distribuzioni gaussiane in 3D dà PDF velocità scalare 4 pi v^2 exp(-m v^2 / (2 k T)). Maxwell electromagnetic theory 1865 Treatise 1873 'A Treatise on Electricity and Magnetism' equazioni Maxwell unificazione elettricità magnetismo luce 4 leggi fondamentali. Ludwig Boltzmann (Austria 1844-1906) Vienna Graz Universität fisica matematica, Leipzig Universität 1900-1902, Vienna Universität 1902-1906. Boltzmann 1872 paper Wiener Berichte 66:275 'Weitere Studien über das Wärmegleichgewicht unter Gasmolekülen' H-theorem derivazione kinetic theory rigorosa equilibrium distribution monotonic decrease H function entropy increase second law thermodynamics. Boltzmann 1877 paper Wiener Berichte 76:373 'Über die Beziehung zwischen dem zweiten Hauptsatze der mechanischen Wärmetheorie und der Wahrscheinlichkeitsrechnung' celebrated equation S = k log W entropy Boltzmann constant k probability microstates W engraved tombstone Vienna Zentralfriedhof.

Boltzmann statistical mechanics ergodic hypothesis 1871 ergodische Hypothese time average = ensemble average. Boltzmann tragic suicide Duino Italy near Trieste September 5 1906 mental illness depression. Erwin Schrödinger 1944 'What is Life?' Cambridge famous quote Boltzmann work fundamental quantum mechanics.

Continuazione fisica statistica: Albert Einstein (Germany-USA 1879-1955) Bern Patent Office 1905 'Annus Mirabilis' year photoelectric effect light quanta hv (Nobel 1921) Brownian motion paper 'Über die von der molekularkinetischen Theorie der Wärme geforderte Bewegung von in ruhenden Flüssigkeiten suspendierten Teilchen' Annalen der Physik 17:549 confirmed Boltzmann molecular kinetic theory experimental molecule reality. Jean Perrin (France 1870-1942) Sorbonne Paris Nobel Physics 1926 mesure Avogadro constant Brownian motion confirmation atomic theory.

Italia fisica statistica: Enrico Fermi (Italia-USA 1901-1954) Sapienza Roma 1922-1938 + Chicago + Los Alamos. Fermi nato Roma figlio capostazione FFSS, Pisa SNS Scuola Normale Superiore PhD 1922, Roma La Sapienza Faculty Physics 1922-1938 Faculty Director, Premio Nobel Fisica 1938 'demonstrations of the existence of new radioactive elements produced by neutron irradiation' neutron lenti reazioni nucleari ricerca Roma 1934-1938 fellow Roma. Fermi emigrazione USA 1938 Stockholm Nobel ceremony cancellation rientro Italia leggi razziali fasciste antisemitismo moglie Laura Capon ebrea, Columbia University NYC 1939, University Chicago 1942 Chicago Pile-1 prima reazione nucleare controllata 2 dicembre 1942 squash court below Stagg Field, Manhattan Project Los Alamos NM Robert Oppenheimer collaboratore atomic bomb Trinity test 1945. Fermi statistics fermion Dirac-Fermi 1926 + Bose-Einstein bosoni 1924; Fermi golden rule perturbation theory; Fermi paradox extraterrestrial life equation Drake. Italia fisici Roma scuola Fermi Boys Via Panisperna: Edoardo Amaldi (1908-1989), Bruno Pontecorvo (1913-1993) emigrato URSS Dubna, Emilio Segre (1905-1989) Berkeley Nobel 1959 antiprotone, Franco Rasetti (1901-2001) Quebec, Ettore Majorana (1906-1938 sparito) Naples Univ.

Italia istituzioni fisica: INFN Istituto Nazionale Fisica Nucleare Frascati 1951 fondato Edoardo Amaldi 5000 ricercatori 28 sezioni LNF Frascati Laboratori Nazionali Frascati DAFNE accelerator KLOE detector; INFN Cosenza Catania Lecce sezioni meridionali; INFN LNGS Laboratori Nazionali Gran Sasso 1987 fondato OPERA neutrino + Borexino solar neutrinos + LUNA underground dark matter; CERN Geneva Switzerland 1954 contributi italiani forti; INAF Istituto Nazionale Astrofisica Roma 1999 fondato 24 osservatori Brera Milano Padova Capodimonte Napoli; OAR Osservatorio Astronomico Roma Monte Porzio Catone; ASI Agenzia Spaziale Italiana Roma 1988; ENEA Centro Ricerche Frascati Roma + Bologna + Brasimone Bologna fusione nucleare ITER.

Software fisica statistica: Python NumPy SciPy + matplotlib Anaconda; Wolfram Mathematica Champaign IL Stephen Wolfram founder; MATLAB Statistics Toolbox MathWorks Natick MA; Origin OriginLab Northampton MA scientific graphing; LAMMPS Sandia National Laboratories molecular dynamics; GROMACS molecular dynamics University Groningen Netherlands; Quantum ESPRESSO Italia Trieste SISSA + CNR DFT Density Functional Theory ab initio simulations; CASTEP Cambridge UK; VASP Vienna Austria.

Tabella valori distribuzione

• Costante Boltzmann k: 1.380649e-23 J/K esatto SI 2019
• Costante gas R: 8.31446262 J/(mol K)
• Avogadro N_A: 6.02214076e23 mol^-1
• v_p most probable: sqrt(2 k T / m)
• v_mean media: sqrt(8 k T / (pi m))
• v_rms RMS: sqrt(3 k T / m)
• Rapporto v_p:v_mean:v_rms: 1 : 1.128 : 1.225
• E_k per molecola: (3/2) k T = 6.17e-21 J a 298K
• E_k per mole: (3/2) R T = 3.72 kJ/mol a 298K
• N2 aria 25°C v_p: 421 m/s
• N2 aria 25°C v_mean: 475 m/s
• N2 aria 25°C v_rms: 515 m/s
• O2 aria 25°C v_rms: 482 m/s
• H2 25°C v_rms: 1929 m/s (leggera)
• He 25°C v_rms: 1364 m/s
• Ar 25°C v_rms: 431 m/s
• CO2 25°C v_rms: 411 m/s
• H2O 25°C v_rms: 642 m/s
• CH4 25°C v_rms: 681 m/s
• N2 200K v_rms: 421 m/s (atmosfera alta)
• N2 1000K v_rms: 943 m/s (combustione)
• N2 2000K v_rms: 1334 m/s (plasma)
• Maxwell 1860 paper: Phil Mag 19:19
• Boltzmann 1872 paper: Wiener Berichte 66:275
• Boltzmann 1877 paper: Wiener Berichte 76:373 S=k log W
• Italia INFN Frascati 1951: 5000 ricercatori
• Italia LNGS Gran Sasso 1987: underground neutrini
• Italia INAF 1999: 24 osservatori
• Italia CNR 1923: 8000 ricercatori
• Italia Marconi 1909: Nobel Fisica radiotelegrafia
• Italia Fermi 1938: Nobel Fisica neutroni lenti
• Italia Segre 1959: Nobel Fisica antiprotone
• Italia Rubbia 1984: Nobel Fisica W Z bosoni
• Italia Giacconi 2002: Nobel Fisica raggi X
• Italia Parisi 2021: Nobel Fisica spin glass

Esempi numerici dettagliati

Esempio 1 - N2 aria temperatura ambiente. Aria 78% N2 a T=298K (25°C). M=28 g/mol N2. m=4.65e-26 kg. k T=4.11e-21 J. v_p=sqrt(2*4.11e-21/4.65e-26)=421 m/s; v_mean=sqrt(8*kT/(pi*m))=475 m/s; v_rms=sqrt(3*kT/m)=515 m/s. PDF a v=500 m/s: f(500)=2.0e-3 (s/m). Frazione molecole con velocità >1000 m/s: ~5%. Italia: ricerca chimica fisica Sapienza Roma DIET; Politecnico Milano DEIB; CNR ICCOM Florence.

Esempio 2 - Elio He gas leggero. He gas nobile T=298K. M=4 g/mol monoatomico inerte. m=6.65e-27 kg. v_p=sqrt(2*4.11e-21/6.65e-27)=1113 m/s; v_mean=1257 m/s; v_rms=1364 m/s. He molecole più veloci aria 2.7x dovuto massa ridotta. Atmosfera Terra: He scappa gravitazionalmente velocità >escape velocity 11.2 km/s. SAES Getters Lainate 1940 produttore He ultra-puro 1500 employees ricerca criogenia.

Esempio 3 - O2 alta temperatura forno. O2 ossigeno T=500K (forno). M=32 g/mol. m=5.31e-26 kg. v_p=494 m/s; v_mean=557 m/s; v_rms=605 m/s. PDF a v=600 m/s: 1.6e-3 (s/m). Aumento T da 298K a 500K (factor 1.68): v_rms aumenta sqrt(1.68)=1.30. Italia: industria petrolchimica Eni Milano 1953; Saipem oil&gas; Polynt Scanzorosciate chimica intermedi.

Esempio 4 - Idrogeno H2 atmosfera fredda. H2 gas leggero T=100K (atmosfera alta troposfera). M=2 g/mol. m=3.32e-27 kg. v_p=911 m/s; v_mean=1029 m/s; v_rms=1117 m/s. H2 escape Terra mass minor leggera + v alta = perdita atmosferica geologica (~ no H2 atmosfera attuale Terra). Marte H2 escape ancora più veloce gravità minore. Italia: ASI Agenzia Spaziale Italiana Roma 1988; INAF planetary science; INGV Roma sismologia.

Esempio 5 - CO2 combustione motori. CO2 motori interni T=2000K combustione. M=44 g/mol. m=7.31e-26 kg. v_p=775 m/s; v_mean=874 m/s; v_rms=950 m/s. Velocità alte propagazione fiamma deflagrante 1-100 m/s; detonante 1000-3000 m/s. Combustione scoppio motore Otto cycle ciclo termodinamico. Italia: Politecnico Torino DET motori interni Lancia FIAT Iveco; Eni Milano combustibili.

Limiti e considerazioni

1. Solo gas ideali equilibrio termico: distribuzione applicabile gas diluiti basse densità no interazioni intermolecolari (Lennard-Jones gas reali deviazioni).
2. Equilibrio termico richiesto: non valida sistemi non-equilibrio (jet supersonici Mach>1, plasmi non thermal, raggi cosmici).
3. Classica non-quantistica: a temperature molto basse (T<1K) o densità molto alte (degenerazione) Maxwell-Boltzmann sostituita da Bose-Einstein bosoni o Fermi-Dirac fermioni.
4. Solo traslazionale: distribuzione velocità molecole singolo grado libertà traslazionale 3D, non considera rotazione vibrazione (gas poliatomici).
5. Particelle classiche: assume particelle distinguibili classiche, fallisce per particelle identiche quantistiche (fermioni boson).
6. Massa fissa: assume tutte molecole stessa massa m, miscele gas distribuzione media pesata.
7. Volume infinito: sistema infinito non considera effetti finite size pareti contenitore, importanti micro/nanofluidica.
8. Stati stazionari: distribuzione equilibrio non considera transitori velocità di rilassamento ~ tempo collisione 10^-9 s.
9. Symmetry breaking: campi esterni (gravità, magnetico, elettrico) possono rompere simmetria, distribuzione anisotropa.
10. Quantum effects: a T molto basse de Broglie wavelength comparabile spaziatura intermolecolare, quantum corrections necessarie (Bose-Einstein condensation 1995 JILA Boulder Colorado).

Errori comuni

Confondere k_B costante Boltzmann e R costante gas: k_B = 1.381e-23 J/K per molecola singola; R = 8.314 J/(mol K) per mole. Relazione R = k_B * N_A. Energia cinetica media (3/2) k_B T per molecola; (3/2) R T per mole.

Mass molare g/mol vs kg: M tipicamente espressa g/mol; convertire kg dividendo 1000 + dividere per N_A per massa molecola singola: m = M / 1000 / N_A. Errore comune dimenticare conversione.

Confondere v_p, v_mean, v_rms: tre velocità distinte v_p<v_mean<v_rms (rapporto 1:1.128:1.225). v_rms usata energia cinetica E_k = (1/2) m v_rms²; v_mean usata flusso effusione attraverso fori; v_p picco distribuzione.

Temperatura Celsius vs Kelvin: T deve essere in Kelvin (assoluto). T(K) = T(°C) + 273.15. Errore frequente usare °C dà valori non fisici (negativi o T=0).

PDF dimensions: f(v) ha dimensioni 1/(m/s) = s/m (probabilità per unità velocità). Integral f(v) dv = 1 senza dimensioni.

Applicazioni pratiche

1. Teoria cinetica gas: pressione gas P = nkT (legge gas ideali); temperatura cinetica E_k=(3/2)kT; viscosità diffusione coefficienti.
2. Effusione Graham law: rate effusione attraverso foro proporzionale 1/sqrt(M) molar mass. UF6 isotope separation uranium enrichment Manhattan Project Oak Ridge K-25 plant.
3. Combustione motori: distribuzione velocità reagenti determina rate reazione Arrhenius E_a / k_B T frazione molecole con energia >E_a.
4. Astronomia stellare: temperatura plasma stellare 10^4-10^7 K; velocità ioni MeV-GeV in solar wind; nebulae HII T=10^4 K.
5. Atmosfere planetarie: escape velocity vs v_rms determines retention atmosfera. Earth retains N2 O2 (v_rms<<v_escape); Moon Mars lose H2 He.
6. Plasmi laboratorio: tokamak fusion T=10^7-10^8 K v_rms ioni 10^6 m/s. Italia: ENEA Frascati FTU Frascati Tokamak Upgrade; ITER Cadarache France contribution Italia.
7. Spettroscopia molecolare: Doppler broadening linee spettrali f_D ~ f_0 v/c proporzionale v_rms; risoluzione spettrometri.
8. Industria chimica: cracking petrolio pirolisi alta T; sintesi ammoniaca Haber-Bosch 400-500°C 200 atm; reattori chemical engineering.
9. Refrigerazione magnetica: laser cooling atoms 10^-6 K Bose-Einstein condensation 1995 JILA Boulder Colorado Wieman Cornell Nobel 2001.
10. Ingegneria aerospaziale: atmosphere entry shock waves Mach numbers; rocket combustion exhaust 3000 K 4000 m/s velocity. Italia: ASI Roma 1988; Avio Colleferro propulsione spaziale Vega rocket.

Aggiornato a luglio 2026