320 likes | 488 Views
Classical N Body Dynamics? Sous titre à la grammaire incertaine et à l’orthographe approximative!. A.Bonasera, V.Latora J.Randrup, C.Dorso, A.Strachan, A.Chernomoretz J.Bondorf et al. S.Pratt et al. X.Campi, H.Krivine, N.Sator, E.Plagnol, F.Lavaud T.J.Schlagel, V.R.Pandharipande
E N D
Classical N Body Dynamics? Sous titre à la grammaire incertaine et à l’orthographe approximative! A.Bonasera, V.Latora J.Randrup, C.Dorso, A.Strachan, A.Chernomoretz J.Bondorf et al. S.Pratt et al. X.Campi, H.Krivine, N.Sator, E.Plagnol, F.Lavaud T.J.Schlagel, V.R.Pandharipande M.Ison, F.Gulminelli Classical N Body Dynamics D.Cussol D.C. PRC 65 (2002) 054614 D.C. PRC 68 (2003) 014602 Thèse F.Morisseau Séminaire Informel 16 mai 2007
Modeles Modèles Dynamiques Formation de noyaux excités + émissions PE Modèles Statistiques Désexcitation des noyaux chauds LV AMD BUU • GEMINI QMD(s) tcouplage? • SIMON • SMM • MMMC • EES E*/A,T J A,Z temps Séminaire Informel 16 mai 2007
Test QI Rubriques A BracTome 1 “Testons notre Q.I.” M.Gotlib Données Modèles Séminaire Informel 16 mai 2007
CNBD Classical NBody Dynamics • Contient • N corps en interaction forte • Equations classiques du mouvement • interaction à 4 paramètres • Lois de conservation (1 paramètre de simulation) • Ne contient pas • pas d’interaction Coulombienne • pas de corrections quantiques • pas de désexcitation secondaire N’est pas fait pour reproduire les données! Ce n’est pas de la physique nucléaire! Ce n’est pas de la théorie! (expériences numériques) C’est juste "pour voir"! D.Cussol PRC 65 (2002) 054614 Séminaire Informel 16 mai 2007
Un peu de vocabulaire C’est Navrant de Bétise Daniel Un peu de vocabulaire • Unités arbitaires • Distance Simulation Units • Energy Simulation Units • Time Simulation Units • Energies • Ecm:énergie cinétique du projectile et de la cible dans le centre de masses du système. • Edisp: énergie disponible dans le centre de masse du système • E* énergie d’excitation Séminaire Informel 16 mai 2007
statique Propriétés statiques Eleast_bound Ebind Mini cristaux! "Equation d’état" r’=r eff=(1/)3 Séminaire Informel 16 mai 2007
collisions Neck Multifragmentation Fusion Participant/Spectator Binary collision Deep Inelastic Collision Particle(s) exchange Quelques collisions les animations sont visibles à http://caeinfo.in2p3.fr/theorie/CNBD/node15.html Séminaire Informel 16 mai 2007
Mimf vs b Comment Nous Barrer Discrètement? Collisions binaires Fusion M(N > 2) Multifragmentation b/bmax Multifragmentation Energy close to Eleast_bound Low Energy Deux particules i et j appartiennent au même fragment si rij < r1 Energy close to Ebind High Energy Séminaire Informel 16 mai 2007
Dalitz Fusion/ Evaporation "Fission" ELeastBound (68) D1 D2 D3 Multifragmentation Dalitz plots CNBD, 34+34, bred < 0.1 N1 > N2 > N3 > ... Di=Ni/(N1+N2+N3) Séminaire Informel 16 mai 2007
E vs V-N 50+50 ELeastBound(N) ELeastBound(100) EBind (N)/N E*/N = f(V//) Limitation de l’énergie d’excitation E*/N < ELeastBound D.Cussol PRC68 (2003) 014602 Séminaire Informel 16 mai 2007
Thermalisation? Multifragmentation Fusion Thermalisation? Multiplicité Taux de production Séminaire Informel 16 mai 2007
Freeze-outs Freeze-out clusters Freeze-out Vz Freeze-out et Freeze-out 50 + 50 Ecm/N=120 E.S.U. 0,3 < b/bmax < 0,4 D.Cussol, soumis à PRC. Séminaire Informel 16 mai 2007
Calcul Nul Bidouillé par Daniel INDRA DATA Xe + Sn 45 MeV/A CNBD Nproj=50 Ntarg=50 Ecm/N=120 E.S.U. Cartes taille-Vz Z-Vz Séminaire Informel 16 mai 2007
Ta + Au + lourd Nproj=50 Ntarg=50, Ecm/N=120 E.S.U. Ta+Au at 39 MeV/A INDRA data CNBD ZMAX Effet de hiérarchie Le fragment le plus lourd est le plus rapide et sa direction d’émission est proche de celle du QP. J.Colin et al (INDRA collaboration) PRC67 (2003) 064603 J.Normand, PhD Thesis, Université de Caen (2001) F.Bocage et al (INDRA collaboration) NPA676 (2000) 391 B.Davin et al, PRC65 (2002) 064614 L.Gingras et al, PRC65 (2002) 054613 Z MAX-1 Z MAX-2 Z MAX-3
Vr vs cos(Tetaprox) Crazy Naughty Boy Daniel CNBD Nproj=50 Ntarg=50 Ecm/N=120 E.S.U. INDRA Data Ta + Au 39 MeV/u Corrélations internes J.Colin et al (INDRA collaboration) PRC67 (2003) 064603 D.Cussol, soumis à PRC.
Cneg Data constant P constant V M.D’Agostino et al, PLB473 (2000) 219 (Multics + miniball, Au + Au) M.D’Agostino et al, NPA699 (2002) 795 (INDRA, Xe + Sn) C<0 : Données expérimentales F.Gulminelli, P.Chomaz, PRL82 (1999) 1402 Séminaire Informel 16 mai 2007
Determination 50+50 bred < 0.1 ELeastBound EBind/N Détermination de C Séminaire Informel 16 mai 2007
INDRA Au + Au 80 MeV/u Comment Nous Bullons Délicieusement Bimodalité: données expérimentales (INDRA) CNBD and bimodality Z1 > Z2 > Z3 > … S’il y a une transition de phase, la distribution d’un paramètre d’ordre doit avoir deux "bosses". P.Chomaz, F.Gulminelli, Physica A330 (2003) 451, NPA749 (2005) 3c M.Pichon, B.Tamain, R.Bougault, O.Lopez and INDRA NPA749 (2005) 93c Séminaire Informel 16 mai 2007
Bimodalités ELeastBound Bimodalité dans CNBD Collisions centrales b/bmax < 0.1 Séminaire Informel 16 mai 2007
Conclusion 1 • Pour la multifragmentation: • collisions centrales • ELeastBoundest le seuil de fragmentation • Saturation de l’énergie d’excitation: E*/N ELeastBound • Fragmentation de systèmes incomplètement thermalisés • Signaux de transition de phase Etransition/N ELeastBound • Collisions périphériques • Les effets de hiérarchie et les corrélations internes sont présents sans ajouter de paramètres • “Gel” de la distribution des vitesses avant que les fragments soient séparés On a un accord qualitatif surprenant (troublant?) de CNBD avec les données expérimentales de collisions noyau-noyau et ceci dans un cadre UNIQUE. Constats Séminaire Informel 16 mai 2007
Three Physical configurations Systèmes pré-thermalisés E* = Epot + Ekin -Ebind(N) Rconf E* = Epot + Ekin -Ebind(N) Dans le centre de masses Collisions centrales b/bmax < 0.1 Cussol Nous Barbe Drôlement Deux configurations physiques ou ? Séminaire Informel 16 mai 2007
t < tconf t = tconf t > tconf Pré-thermalisation des systèmes dans CNBD Etape 1 Etape 2 Etape 3 Rconf Deux particules i et j appartiennent au même fragment si rij < r1 (portée) ttherm < tconf Séminaire Informel 16 mai 2007
50+50 bred < 0.1 ELeastBound Ebind/N C < 0 dans CNBD N=100 Rconf =16.41 D.S.U. N=100 Rconf =17.68 D.S.U. • La transition dépend du volume de confinement. • Collisions centrales compatibles avec un système pré-thermalisé très comprimé N=100 Rconf =42.94 D.S.U. N=100 Rconf =59.55 D.S.U. Séminaire Informel 16 mai 2007
Bimodalité(s) • L’énergie de transition diminue lorsque le rayon de confinement diminue. • Collisions centrales compatibles avec un système pré-thermalisé très comprimé Séminaire Informel 16 mai 2007
Comparaison des distributions de taille s CNBD Ntot=100 ≈ 0 Séminaire Informel 16 mai 2007
Comparaison des distributions de taille CNBD Ntot=100 ≈ 0/8 Séminaire Informel 16 mai 2007
Comparaison des distributions de taille s CNBD Ntot=100 ≈ 2,5 0 Séminaire Informel 16 mai 2007
Chérissons Notre Bon Daniel Distributions angulaires cos() Séminaire Informel 16 mai 2007
Distributions angulaires cos() Séminaire Informel 16 mai 2007
Emoinslié Centrales b/bmax < 0.1 Processus hors équilibre Processus thermalisés 50+50 E*/N ELeastBound ELeastBound Collisions Systèmes pré-thermalisés Collisions Systèmes pré-thermalisés dilués Systèmes pré-thermalisés comprimés Ebind/N Rôle de ELeastBound dans CNBD N ELeastBound(N) Séminaire Informel 16 mai 2007
? Récapitulons Récapitulation ? Séminaire Informel 16 mai 2007
Commandes Neuf Bières Daniel! Questions Questions • Pourquoi existe-t-il une telle similitude entre CNBD et les collisions noyau-noyau? • Effets quantiques faibles? • Comportements communs à tous les systèmes à N corps? • Prédominance des corrélations à N corps? • Le concept d’équation d’état de la matière infinie est-il adapté pour comprendre les mécanismes de collisions entre systèmes à N corps? • Effets de nombres finis déterminants • Géométrie de la collision importante • Est-il toujours pertinent d’utiliser des modèles "thermiques" pour calculer les taux de production des particules et des fragments? • Energie d’excitation limitée • Effets de voie d’entrée déterminants: ils DOIVENTêtre pris en compte • Les distributions des vitesses sont figées bien avant que les fragments soient séparés. Les modèles dynamiques sont notre seule issue... Séminaire Informel 16 mai 2007