Des neutrons atmosphériques à la capture du bore : jusqu’où va l’idée ?

En 2011, nous avons posé une question inhabituelle : les neutrons présents à l’altitude des avions peuvent-ils déclencher des captures dans une tumeur enrichie en bore-10 ? Le mécanisme existe. Mais notre calcul trouve une dose moyenne trop faible pour traiter une tumeur macroscopique. L’intérêt de cette étude est de quantifier cette limite, avant de promettre une application.

Explorer le modèle · Lire le résultat · Article et équipe

Pourquoi cette question appartient à mon parcours

La physique médicale et la prévision solaire partagent ici une méthode : partir du phénomène, expliciter la géométrie et les hypothèses, simuler, puis vérifier si l’ordre de grandeur répond au besoin. L’IA n’est pas nécessaire pour poser cette question. Le modèle sert aussi à reconnaître une piste qui ne suffit pas, et à préciser ce qui resterait à démontrer.

Une cascade dans l’atmosphère, puis une réaction à l’échelle cellulaire

Les rayons cosmiques galactiques, principalement des protons et des noyaux énergétiques, produisent des particules secondaires en rencontrant l’atmosphère, dont des neutrons. Le champ magnétique terrestre, l’altitude et le spectre des particules comptent. Le vent solaire module l’arrivée des rayons cosmiques ; il ne faut pas le confondre avec les événements de particules solaires énergétiques. Le fond galactique est généralement plus élevé au minimum solaire. Le schéma ci-dessous est qualitatif : il ne calcule pas une météorologie de l’espace.

Le résultat : une limite quantitative, pas un hôpital volant validé

Dans le cas de référence de l’article : altitude 12 km, latitude 45°, flux neutronique de 14,7 cm⁻²·s⁻¹, couche d’eau de 10 cm, crâne de 0,5 cm et tumeur d’environ 1 g contenant 20 µg de bore-10 par gramme. Le modèle hybride associe transport Monte-Carlo et approximations analytiques. Il estime une dose tumorale moyenne de l’ordre de 1 µGy pour une heure. C’est insuffisant pour une irradiation thérapeutique macroscopique.

Pour la concentration de référence uniquement, l’outil illustre D̄ ≈ 1 µGy × (t / 1 h). Un autre panneau représente I = 100 % × (t / 1 h) × (c₁₀B / 20 µg·g⁻¹) : un indice relatif du nombre d’interactions sur le bore, normalisé à 100 % pour une heure à 20 µg/g. Il ne donne ni une probabilité de capture ni une dose totale à une autre concentration. Spectre, géométrie et masse restent fixés ; le transport de l’article n’est pas réexécuté. La concentration est celle d’un fantôme numérique, pas une quantité à administrer à une personne.

Ce que l’étude apporte à la modélisation

Le travail relie trois niveaux habituellement séparés : les réactions dans les molécules porteuses du bore, le transport à travers un modérateur et les tissus, puis les ordres de grandeur de dose. Sa seconde contribution est de remplacer certaines étapes coûteuses du transport par des expressions analytiques de transmission et de libre parcours moyen. L’article rapporte 42 secondes au lieu de 15 heures pour 6 000 neutrons, dans les conditions de calcul de l’époque. Ce n’est pas un benchmark de ce navigateur : c’est un exemple concret de simplification physique au service d’une question scientifique.

Monte-Carlo : suivre beaucoup d’histoires de neutrons

Une section efficace microscopique σ(E), exprimée en unités de surface, caractérise l’interaction d’un neutron d’énergie E avec un noyau. Dans un matériau, Σ = Σᵢ Nᵢσᵢ est la section macroscopique (inverse d’une longueur), où Nᵢ est la densité de noyaux de l’espèce i. Elle donne un libre parcours moyen 1/Σ. Un tirage uniforme U entre 0 et 1 permet, dans un milieu homogène, de tirer une distance ℓ = −ln(U)/Σ. Le code choisit ensuite le noyau et le canal à partir des sections efficaces, puis suit l’énergie et la direction après collision. Il répète ces histoires et additionne les dépôts d’énergie.

Les résultats statistiques convergent avec le nombre d’histoires ; sous les hypothèses usuelles d’indépendance et de variance finie, l’erreur-type diminue comme 1/√N. Cela ne corrige pas une géométrie ou des données physiques inadéquates. Notre étude de 2011 associe ces tirages à des expressions analytiques ; l’animation du navigateur illustre le mécanisme, sans réexécuter ce transport.

Los Alamos, ENDF/B et le code de l’article

L’article indique avoir utilisé la bibliothèque ENDF/B, avec un accès alors référencé sur t2.lanl.gov, à Los Alamos. Ces évaluations fournissent les sections efficaces et les distributions des réactions. ENDF/B est produite par le groupe CSEWG et distribuée par le NNDC de Brookhaven ; elle ne doit pas être attribuée exclusivement à Los Alamos. Notre code était écrit en Visual Basic. Le papier cite MCNP comme exemple de code de transport utilisant ces évaluations ; cela ne signifie pas que notre simulation utilisait MCNP ni qu’il existait un partenariat institutionnel avec Los Alamos.

Pourquoi le tableau affiche-t-il (n,p) ≈ 0 % ?

Il s’agit du résultat arrondi de la comparaison neutron–molécule BSH/BPA dans le spectre étudié (tableau 2), pas d’une impossibilité physique. La méthode de l’article inclut les canaux (n,n), (n,γ), (n,p) et (n,α). Dans un tissu, la réaction ¹⁴N(n,p)¹⁴C constitue aussi une composante de la dose BNCT, distincte de la capture sur le bore, comme l’explique l’AIEA, Advances in Boron Neutron Capture Therapy, tableau 7. Ce panneau ne calcule pas ces différentes composantes de dose. Sur le dessin, les épaisseurs sont indiquées par les nombres ; la largeur de l’os est exagérée pour rester visible.

BSH ou BPA : ce que la comparaison dit vraiment

La BNCT (Boron Neutron Capture Therapy, thérapie par capture neutronique du bore) utilise des vecteurs du bore tels que le BSH et la BPA. Le tableau 2 compare les réactions lors de l’interaction neutron–molécule dans le spectre étudié : le canal (n,α) représente 94 % pour le BSH et 88 % pour la BPA. Ces distributions conditionnelles ne donnent ni la probabilité de capture de chaque neutron du vol, ni la sélectivité tumorale, ni un taux de succès clinique. À concentration de bore-10 identique, nous ne leur attribuons donc pas un multiplicateur de dose différent.

Très local ne veut pas dire partout

Après capture, ¹⁰B + n → ⁷Li + α : les particules chargées parcourent quelques micromètres (environ 4 µm pour le lithium et 8 µm pour l’alpha dans l’article). Le texte rapporte aussi une dose locale d’environ 0,5 Gy par impact à l’échelle cellulaire. Cette valeur conditionnelle ne signifie pas que chaque cellule reçoit cette dose. Le présent outil ne reproduit pas cette microdosimétrie : le volume de score et la fraction d’énergie déposée devraient être explicités pour la recalculer.

La dose du vol est une autre grandeur

Pour un cylindre simplifié de 56,5 kg, l’article trouve 4,37 µSv/h, comparés à une référence historique SIEVERT Paris–New York de 6,1 µSv/h : environ 28 % de moins, avec des périodes solaires et des approximations différentes. Cette comparaison d’ordre de grandeur n’est pas un calculateur de vol actuel. Le sievert (Sv) pondère l’exposition radiologique ; le gray (Gy) mesure l’énergie absorbée par unité de masse. Les deux panneaux ne se convertissent pas directement l’un dans l’autre.

Ce que l’on pourrait rechercher ensuite

L’article envisage les micrométastases et une piste adjuvante, tout en demandant des expériences biologiques. Il faudrait notamment établir la distribution du bore, la sélectivité, la réponse aux faibles doses, la microdosimétrie et le bénéfice par rapport aux risques. « Transformer les avions en salles de radiothérapie » est une perspective exploratoire de la conclusion, pas une efficacité démontrée. Cette expérience n’est ni une prescription ni un outil de planification clinique.

L’article et toute l’équipe

Cyril Voyant, Rudy Roustit, Jennifer Tatje, Katia Biffi, Delphine Leschi, Jérome Briançon, Céline Lantieri Marcovici. Therapeutic potential of atmospheric neutrons. Reports of Practical Oncology and Radiotherapy, 16(1), 21–31 (2011). DOI: 10.1016/j.rpor.2010.11.002 · Texte intégral chez PMC.

Origine des particules et modulation solaire : NOAA, Galactic Cosmic Rays · FAA, What aircrews should know about their occupational exposure.

Poursuivre : LQ–LQL : calculer lorsque la relation est connue · Le fil de mes recherches

Dans un autre domaine : le modèle ne fait pas tout, variabilité et mémoire en prévision solaire, caractériser le problème avant de choisir l’outil.