Calcolatore effetto tunnel quantistico barriera

Il calcolatore effetto tunnel quantistico determina la probabilità di trasmissione T di una particella attraverso una barriera di potenziale rettangolare di altezza V₀ e larghezza L con energia E inferiore alla barriera (E < V₀), fenomeno puramente quantistico privo di analogo classico previsto da Friedrich Hund 1927 nello studio dell'inversione della molecola di ammoniaca NH₃ + sviluppato indipendentemente da George Gamow 1928 Z. Phys. 51 e Ronald Gurney + Edward Condon 1928 Nature 122 per spiegare il decadimento alfa nuclei pesanti (U-238 t₁/₂ = 4,5 Gyr). La formula classica barriera rettangolare: T ≈ 16·E·(V₀-E)/V₀² · exp(-2·κ·L), dove κ = √(2m(V₀-E))/ℏ costante di decadimento esponenziale dentro la barriera, m massa particella in kg, V₀ altezza barriera in J, E energia incidente in J, L larghezza barriera in m, ℏ = 1,0546×10⁻³⁴ J·s costante di Planck ridotta. Il prefattore 16·E·(V₀-E)/V₀² rappresenta il termine di accoppiamento d'onda alle interfacce, mentre l'esponenziale exp(-2κL) domina la soppressione tunneling quando κL ≫ 1.

In approssimazione WKB Wentzel-Kramers-Brillouin 1926 per barriera generale V(x) la formula diventa T ≈ exp(-2·∫√(2m(V(x)-E))/ℏ · dx) integrale tra punti di inversione classica. La profondità di penetrazione 1/κ indica la scala spaziale di decadimento esponenziale dell'ampiezza di probabilità nella regione classicamente proibita. Le costanti fondamentali: ℏ = 1,0546×10⁻³⁴ J·s, m_e = 9,109×10⁻³¹ kg massa elettrone, e = 1,602×10⁻¹⁹ C carica elementare, 1 eV = 1,602×10⁻¹⁹ J. Le 10 applicazioni principali:

• decadimento alfa nuclei pesanti Gamow 1928 U-238 Th-232 Ra-226 (curva Geiger-Nuttall log T vs Q-value)
• STM Scanning Tunneling Microscope Gerd Binnig + Heinrich Rohrer 1981 IBM Zürich Nobel 1986 (corrente exp(-2κL) sensibile 0,01 nm)
• diodo tunnel Esaki Leo Esaki 1957 Sony Nobel 1973 (giunzione p-n degenerata, resistenza differenziale negativa NDR)
• giunzione Josephson Brian Josephson 1962 Nobel 1973 (Cooper pair tunneling, SQUID magnetometri 10⁻¹⁵ T)
• Flash memory NAND Floating Gate Fowler-Nordheim tunneling SiO₂ 8 nm
• Field emission displays + emissione fredda elettroni
• fusione nucleare stellare D + T → ⁴He + n (Sole pp-chain 4·H → ⁴He, Gamow factor exp(-2π·Z₁Z₂e²/ℏv))
• decadimento beta W boson tunneling electroweak
• STMicroelectronics SGS-Thomson Catania 1987 + Agrate Brianza fab semiconduttori 28-12 nm Flash memory tunnel ossido
• computing quantistico qubit transmon Josephson junction superconduttore Al/AlO_x

Tabella applicazioni tunnel quantistico italiana

(Fonte: Hund 1927 + Gamow 1928 Z. Phys. + Gurney-Condon 1928 Nature + Esaki 1957 + Binnig-Rohrer 1981 + INFN LNGS Gran Sasso).

Esempi calcolo italiano effetto tunnel

• Decadimento alfa U-238 INFN LNGS Gamow 1928: m = 4·u = 6,64×10⁻²⁷ kg = 7.344·m_e, V₀ = 27 MeV (barriera coulombiana 2·Z·e²/(4πε₀·R_nuc)), E = 4,2 MeV (Q-value), L ≈ 2·R_nuc = 15 fm. κ = √(2·6,64×10⁻²⁷·(27-4,2)·10⁶·1,602×10⁻¹⁹)/1,0546×10⁻³⁴ ≈ 1,5×10¹⁵ /m. exp(-2κL) ≈ exp(-45) ≈ 10⁻²⁰ probabilità per attempt. Frequenza vibrazione 10²¹ Hz → λ = 10⁻²⁰ · 10²¹ = 10 /s → t₁/₂ = ln2/λ. Gamow factor riproduce Geiger-Nuttall curva log T vs Q^(-1/2) lineare 5 ordini grandezza t₁/₂ da μs a 10¹⁷ s (Th-232 14 Gyr)
• STM Au-Au tip Binnig-Rohrer 1981 INFM-CNR Trieste: m = m_e, V₀ = 5 eV (work function Au φ), E = 4 eV, L = 1 nm gap. κ = √(2·9,109×10⁻³¹·1·1,602×10⁻¹⁹)/1,0546×10⁻³⁴ = 5,1×10⁹ /m. T = 16·4·1/25 · exp(-10,2) ≈ 2,56 · 3,7×10⁻⁵ ≈ 10⁻⁴. Sensibilità Δz = 0,01 nm corrente cambia ×e ≈ 2,7 (risoluzione atomica 0,1 nm laterale)
• Esaki diode Si STMicro Catania 1987: m = 0,07·m_e (massa effettiva Si), V₀ = 1 eV gap, L = 10 nm (zona svuotamento Si degenerato p+/n+ doping 10²⁰ /cm³). κ = √(2·0,07·9,1×10⁻³¹·0,3·1,6×10⁻¹⁹)/1,05×10⁻³⁴ ≈ 4×10⁸ /m. T ≈ 16·0,7·0,3·exp(-8) ≈ 10⁻³. Resistenza differenziale negativa per V_bias 0,1-0,3 V (oscillatori microonde)
• Josephson SQUID INRIM Torino metrologia: 2 superconduttori Nb separati ossido AlO_x 1 nm. Cooper pair tunneling I_c = (πΔ/2eR_n) tanh(Δ/2kT). T ≈ 1 a zero-bias (resistenza zero). SQUID sensibilità campo magnetico 10⁻¹⁵ T (MEG magnetoencefalografia, riferimento volt SI Φ₀ = h/2e = 2,068×10⁻¹⁵ Wb)
• Fusione D+T ENEA Frascati DTT 2028: V₀ = 1 MeV (Coulomb d-t Z₁Z₂=1·1, R = ℏc/(αm_pc²) ~ fm), E = kT = 10 keV (plasma 100 milioni K), m_red = m_d·m_t/(m_d+m_t) = 1,2·u. Gamow factor exp(-2π·Z₁Z₂α·c/v_rel) = exp(-2π·1/137·c/v) v=√(2E/m)=10⁶ m/s → exp(-31)/(...) ≈ 10⁻⁴. ITER France 2025 first plasma + DTT Frascati 2028 confinamento Q>1 break-even
• Flash NAND STMicro Agrate Brianza fab: ossido SiO₂ tunnel 8 nm, V₀ = 3,2 eV (banda conduzione SiO₂), E = 0,1 eV (elettroni canale Si), Vg=15 V Fowler-Nordheim regime. T ≈ 10⁻⁹ ma 10¹³ elettroni/s in 100 μs Vg → carica gate flottante 1 fC retention 10 anni. STMicro 28 nm + Micron Agrate 3D NAND 232-layer 2024
• Transmon qubit IIT Genova 2003: junction Al/AlO_x 2 nm, Cooper pair tunneling E_J/E_C ≈ 50 (regime transmon Koch 2007). Anharmonicity α/2π = -300 MHz per gate qubit ω_01 ≈ 5 GHz (T_1 coerenza 100 μs IBM Quantum 2024)

Ricerca quantistica italiana

La ricerca quantistica italiana:

• INFN Istituto Nazionale Fisica Nucleare 1951 22 sezioni 4 Laboratori Nazionali (LNF Frascati 1957 DAFNE 1999 + KLOE-2 rare decays B-meson tunneling; LNGS Gran Sasso 1987 1.400 m sotto roccia: GERDA + LEGEND-200 200 kg Ge-76 doppio decadimento beta senza neutrini 0νββ tunneling W-boson, CRESST + DAMA-LIBRA dark matter direct detection annual modulation, Borexino solar neutrino CNO cycle 2020 first detection, XENONnT 8,5 t Xe WIMP scattering 6×10⁻⁴⁸ cm² sensibilità; LNS Sud Catania 1976 ciclotrone superconduttore + INFN-Lab Legnaro 1968 Tandem)
• ENEA Frascati DTT 2028 + FTU 1989: tokamak fusione plasma D+T, divertore liquido Ga DTT €600 mln PNRR; ITER France 2025 first plasma ingegneria italiana Ansaldo Nucleare Genova 1853 (cryostat 30 m), ICAS Torino (toroidal field coils Nb₃Sn), Walter Tosto Chieti (vacuum vessel 8.000 t), Simic Camerana (CdC) tubature criogeniche
• STMicroelectronics SGS-Thomson Catania 1987 + Agrate Brianza fab semiconduttori SiC GaN BCD 28 nm 12 nm; FinFET 7 nm + Flash memory tunnel ossido NAND 3D 232-layer 2024 fabbisogno PNRR Chips Act €5 mld
• Università e centri ricerca: Politecnico Milano 1863 + Politecnico Torino 1859 ingegneria elettronica; Università Bologna 1088 + Pisa Scuola Normale 1810 (Andrea Ferrara cosmic dawn); Sapienza Roma 1303 Marcello Conversi Pancini muon 1947; ICTP Trieste 1964 + SISSA Trieste 1978 meccanica quantistica; IIT Genova 2003 quantum technologies + transmon qubit; INRIM Torino 1923 metrologia Josephson SI volt + lattice Cs orologio atomico 10⁻¹⁶
• HPC supercalcolo CINECA Casalecchio Bologna 1969: Marconi100 30 PFlops + Leonardo pre-exascale 250 PFlops 2023 (5° TOP500 mondiale) simulazioni quantum chemistry DFT VASP Quantum ESPRESSO Gaussian
• Premi Nobel italiani: Marconi 1909 wireless tunneling diode anticipo, Fermi 1938 fissione, Segrè 1959 antiprotone, Carlo Rubbia 1984 CERN W/Z bosons electroweak Geneva (UA1 esperimento 2.300 collaboratori), Riccardo Giacconi 2002 X-ray astronomy (Uhuru 1970, Einstein 1978, Chandra 1999), Giorgio Parisi 2021 spin glass complex systems. Leo Esaki 1973 Nobel tunnel diode legato Italia via STMicro Catania 1987 ex SGS-Thomson licenza
• Spazio e satelliti italiani ASI Roma 1988: AGILE 2007 astronomy gamma-ray INAF; AMS-02 ISS 2011 raggi cosmici antimateria positroni; BepiColombo ESA-JAXA 2018 ASI 35% magnetometro Mercurio; Athena ESA 2035 X-IFU INAF Palermo 25%; Euclid ESA 2023 ASI 25% materia oscura
• Particle physics LHC CERN: ATLAS + CMS 2012 Higgs boson scoperta INFN-INAF + Politecnico Milano-Torino-Pisa; ICFA International Linear Collider INFN Pisa partecipazione; Fermilab MicroBooNE/DUNE neutrino INFN. Bonus fiscali italiani 2025: Credito d'Imposta R&S art. 1 c. 200 L. 160/2019 fino 25%; Patent Box L. 178/2020 50%; Credito d'Imposta Industria/Transizione 5.0 art. 38 DL 39/2024 fino 35%; Iperammortamento 250%; PNRR Missione 4 Componente 2 6,9 mld € ricerca quantum + Missione 1 Componente 2 idrogeno 3,2 mld €. Disclaimer educativo: i risultati indicativi presuppongono barriera rettangolare ideale; per applicazioni reali (decadimento alfa, dispositivi semiconduttori) usare WKB integrale completo o Schrödinger numerica + parametri sperimentali (bibliografia INFN, Eur. Phys. J., Phys. Rev.)

Aggiornato a luglio 2026