Beam-beam collision example follow-up (ECP-WarpX#4578)

* Update species * Update checksum --------- Co-authored-by: Remi Lehe <remi.lehe@normalesup.org>
ax3l · Jan 5, 2024 · 920b3e8 · 920b3e8
1 parent 685ae3f
commit 920b3e8
Show file tree

Hide file tree

Showing 2 changed files with 126 additions and 168 deletions.
diff --git a/Examples/Physics_applications/beam-beam_collision/inputs b/Examples/Physics_applications/beam-beam_collision/inputs
@@ -71,8 +71,8 @@ algo.particle_pusher = vay
 #################################
 ########### PARTICLES ###########
 #################################
-particles.species_names = beam1 beam2 pho1 pho2 ele_nlbw1 pos_nlbw1 ele_nlbw2 pos_nlbw2 pho ele pos
-particles.photon_species = pho1 pho2 pho
+particles.species_names = beam1 beam2 pho1 pho2 ele1 pos1 ele2 pos2
+particles.photon_species = pho1 pho2
 
 beam1.species_type = electron
 beam1.injection_style = NUniformPerCell
@@ -117,68 +117,58 @@ beam2.do_classical_radiation_reaction = 0
 pho1.species_type = photon
 pho1.injection_style = none
 pho1.do_qed_breit_wheeler = 1
-pho1.qed_breit_wheeler_ele_product_species = ele_nlbw1
-pho1.qed_breit_wheeler_pos_product_species = pos_nlbw1
+pho1.qed_breit_wheeler_ele_product_species = ele1
+pho1.qed_breit_wheeler_pos_product_species = pos1
 
 pho2.species_type = photon
 pho2.injection_style = none
 pho2.do_qed_breit_wheeler = 1
-pho2.qed_breit_wheeler_ele_product_species = ele_nlbw2
-pho2.qed_breit_wheeler_pos_product_species = pos_nlbw2
-
-ele_nlbw1.species_type = electron
-ele_nlbw1.injection_style = none
-ele_nlbw1.self_fields_required_precision = 1e-11
-ele_nlbw1.self_fields_max_iters = 10000
-ele_nlbw1.do_qed_quantum_sync = 1
-ele_nlbw1.qed_quantum_sync_phot_product_species = pho
-ele_nlbw1.do_classical_radiation_reaction = 0
-
-pos_nlbw1.species_type = positron
-pos_nlbw1.injection_style = none
-pos_nlbw1.self_fields_required_precision = 1e-11
-pos_nlbw1.self_fields_max_iters = 10000
-pos_nlbw1.do_qed_quantum_sync = 1
-pos_nlbw1.qed_quantum_sync_phot_product_species = pho
-pos_nlbw1.do_classical_radiation_reaction = 0
-
-ele_nlbw2.species_type = electron
-ele_nlbw2.injection_style = none
-ele_nlbw2.self_fields_required_precision = 1e-11
-ele_nlbw2.self_fields_max_iters = 10000
-ele_nlbw2.do_qed_quantum_sync = 1
-ele_nlbw2.qed_quantum_sync_phot_product_species = pho
-ele_nlbw2.do_classical_radiation_reaction = 0
-
-pos_nlbw2.species_type = positron
-pos_nlbw2.injection_style = none
-pos_nlbw2.self_fields_required_precision = 1e-11
-pos_nlbw2.self_fields_max_iters = 10000
-pos_nlbw2.do_qed_quantum_sync = 1
-pos_nlbw2.qed_quantum_sync_phot_product_species = pho
-pos_nlbw2.do_classical_radiation_reaction = 0
-
-pho.species_type = photon
-pho.injection_style = none
-pho.do_qed_breit_wheeler = 1
-pho.qed_breit_wheeler_ele_product_species = ele
-pho.qed_breit_wheeler_pos_product_species = pos
-
-ele.species_type = electron
-ele.injection_style = none
-ele.self_fields_required_precision = 1e-11
-ele.self_fields_max_iters = 10000
-ele.do_qed_quantum_sync = 1
-ele.qed_quantum_sync_phot_product_species = pho
-ele.do_classical_radiation_reaction = 0
-
-pos.species_type = positron
-pos.injection_style = none
-pos.self_fields_required_precision = 1e-11
-pos.self_fields_max_iters = 10000
-pos.do_qed_quantum_sync = 1
-pos.qed_quantum_sync_phot_product_species = pho
-pos.do_classical_radiation_reaction = 0
+pho2.qed_breit_wheeler_ele_product_species = ele2
+pho2.qed_breit_wheeler_pos_product_species = pos2
+
+ele1.species_type = electron
+ele1.injection_style = none
+ele1.self_fields_required_precision = 1e-11
+ele1.self_fields_max_iters = 10000
+ele1.do_qed_quantum_sync = 1
+ele1.qed_quantum_sync_phot_product_species = pho1
+ele1.do_classical_radiation_reaction = 0
+
+pos1.species_type = positron
+pos1.injection_style = none
+pos1.self_fields_required_precision = 1e-11
+pos1.self_fields_max_iters = 10000
+pos1.do_qed_quantum_sync = 1
+pos1.qed_quantum_sync_phot_product_species = pho1
+pos1.do_classical_radiation_reaction = 0
+
+ele2.species_type = electron
+ele2.injection_style = none
+ele2.self_fields_required_precision = 1e-11
+ele2.self_fields_max_iters = 10000
+ele2.do_qed_quantum_sync = 1
+ele2.qed_quantum_sync_phot_product_species = pho2
+ele2.do_classical_radiation_reaction = 0
+
+pos2.species_type = positron
+pos2.injection_style = none
+pos2.self_fields_required_precision = 1e-11
+pos2.self_fields_max_iters = 10000
+pos2.do_qed_quantum_sync = 1
+pos2.qed_quantum_sync_phot_product_species = pho2
+pos2.do_classical_radiation_reaction = 0
+
+pho1.species_type = photon
+pho1.injection_style = none
+pho1.do_qed_breit_wheeler = 1
+pho1.qed_breit_wheeler_ele_product_species = ele1
+pho1.qed_breit_wheeler_pos_product_species = pos1
+
+pho2.species_type = photon
+pho2.injection_style = none
+pho2.do_qed_breit_wheeler = 1
+pho2.qed_breit_wheeler_ele_product_species = ele2
+pho2.qed_breit_wheeler_pos_product_species = pos2
 
 #################################
 ############# QED ###############
@@ -225,10 +215,10 @@ diagnostics.diags_names = diag1
 diag1.intervals = 0
 diag1.diag_type = Full
 diag1.write_species = 1
-diag1.fields_to_plot = Ex Ey Ez Bx By Bz rho_beam1 rho_beam2 rho_ele_nlbw1 rho_pos_nlbw1 rho_ele_nlbw2 rho_pos_nlbw2 rho_ele rho_pos
+diag1.fields_to_plot = Ex Ey Ez Bx By Bz rho_beam1 rho_beam2 rho_ele1 rho_pos1 rho_ele2 rho_pos2
 diag1.format = openpmd
 diag1.dump_last_timestep = 1
-diag1.species = pho1 pho2 pho ele_nlbw1 pos_nlbw1 ele_nlbw2 pos_nlbw2 ele pos beam1 beam2
+diag1.species = pho1 pho2 ele1 pos1 ele2 pos2 beam1 beam2
 
 # REDUCED
 warpx.reduced_diags_names = ParticleNumber ColliderRelevant_beam1_beam2

diff --git a/Regression/Checksum/benchmarks_json/BeamBeamCollision.json b/Regression/Checksum/benchmarks_json/BeamBeamCollision.json
@@ -1,128 +1,96 @@
 {
   "lev=0": {
-    "Bx": 970235841795.1099,
-    "By": 970175630167.5084,
-    "Bz": 51789226.68307851,
-    "Ex": 2.90874251210462e+20,
-    "Ey": 2.908934769466038e+20,
-    "Ez": 1.70819200999153e+17,
-    "rho_beam1": 7.969271809437626e+16,
-    "rho_beam2": 7.969216746820773e+16,
-    "rho_ele": 0.0,
-    "rho_ele_nlbw1": 299758753127495.56,
-    "rho_ele_nlbw2": 289860995480960.4,
-    "rho_pos": 613963689581.3057,
-    "rho_pos_nlbw1": 287891375746148.75,
-    "rho_pos_nlbw2": 295214367940598.1
+    "Bx": 971135657171.612,
+    "By": 971078812454.5405,
+    "Bz": 20140193.235893946,
+    "Ex": 2.9111756162943966e+20,
+    "Ey": 2.9113725115697712e+20,
+    "Ez": 1.0213536367191107e+17,
+    "rho_beam1": 7.970337929028706e+16,
+    "rho_beam2": 7.969213804851568e+16,
+    "rho_ele1": 343600677331163.7,
+    "rho_ele2": 302746939366837.25,
+    "rho_pos1": 333855946230626.06,
+    "rho_pos2": 310879461124837.44
   },
   "beam1": {
-    "particle_opticalDepthQSR": 104868.583557123,
-    "particle_position_x": 0.0015001423941487055,
-    "particle_position_y": 0.001500178877546805,
-    "particle_position_z": 0.004965525054740317,
-    "particle_momentum_x": 6.20719343794678e-15,
-    "particle_momentum_y": 6.1639864036830354e-15,
-    "particle_momentum_z": 6.807872002761295e-12,
-    "particle_weight": 635864683.4991333
+    "particle_opticalDepthQSR": 104909.59461909423,
+    "particle_position_x": 0.001500222221634118,
+    "particle_position_y": 0.0015002445303035634,
+    "particle_position_z": 0.0049656251227976015,
+    "particle_momentum_x": 6.205341799808723e-15,
+    "particle_momentum_y": 6.1592257603817594e-15,
+    "particle_momentum_z": 6.806886719670214e-12,
+    "particle_weight": 635949610.5135971
   },
   "beam2": {
-    "particle_opticalDepthQSR": 104166.96636860794,
-    "particle_position_x": 0.001500093343835549,
-    "particle_position_y": 0.0015001802698114421,
-    "particle_position_z": 0.00496555337137147,
-    "particle_momentum_x": 6.1976310171196696e-15,
-    "particle_momentum_y": 6.1844493116189144e-15,
-    "particle_momentum_z": 6.796731846337045e-12,
-    "particle_weight": 635863931.615053
+    "particle_opticalDepthQSR": 104164.848014815,
+    "particle_position_x": 0.0015001011957527532,
+    "particle_position_y": 0.001500139975740741,
+    "particle_position_z": 0.004965479176845744,
+    "particle_momentum_x": 6.200690794877584e-15,
+    "particle_momentum_y": 6.186048913459861e-15,
+    "particle_momentum_z": 6.7990490255176515e-12,
+    "particle_weight": 635863144.251134
   },
-  "ele": {
-    "particle_opticalDepthQSR": 0.0,
-    "particle_position_x": 0.0,
-    "particle_position_y": 0.0,
-    "particle_position_z": 0.0,
-    "particle_momentum_x": 0.0,
-    "particle_momentum_y": 0.0,
-    "particle_momentum_z": 0.0,
-    "particle_weight": 0.0
+  "ele1": {
+    "particle_opticalDepthQSR": 435.73003117907257,
+    "particle_position_x": 4.882183530045367e-06,
+    "particle_position_y": 4.841391483672882e-06,
+    "particle_position_z": 1.8449175055560687e-05,
+    "particle_momentum_x": 5.6656608489971696e-18,
+    "particle_momentum_y": 5.724425295258085e-18,
+    "particle_momentum_z": 2.6277553331470036e-15,
+    "particle_weight": 2696555.9200674472
   },
-  "ele_nlbw1": {
-    "particle_opticalDepthQSR": 405.16545575146574,
-    "particle_position_x": 5.119886480641899e-06,
-    "particle_position_y": 5.266875857997872e-06,
-    "particle_position_z": 1.855770689791573e-05,
-    "particle_momentum_x": 5.744620543080891e-18,
-    "particle_momentum_y": 5.463133845399354e-18,
-    "particle_momentum_z": 2.562592228511742e-15,
-    "particle_weight": 2359597.3527386785
-  },
-  "ele_nlbw2": {
-    "particle_opticalDepthQSR": 354.8472581769572,
-    "particle_position_x": 5.235909912617143e-06,
-    "particle_position_y": 5.068613695636293e-06,
-    "particle_position_z": 1.6665840940849017e-05,
-    "particle_momentum_x": 4.47916023711595e-18,
-    "particle_momentum_y": 4.092494634513803e-18,
-    "particle_momentum_z": 2.2277004588526963e-15,
-    "particle_weight": 2354612.144268994
-  },
-  "pho": {
-    "particle_opticalDepthBW": 462.76829758918615,
-    "particle_position_x": 6.163512168886876e-06,
-    "particle_position_y": 5.4614856451855956e-06,
-    "particle_position_z": 2.0316590334121627e-05,
-    "particle_momentum_x": 0.0,
-    "particle_momentum_y": 0.0,
-    "particle_momentum_z": 0.0,
-    "particle_weight": 2923343.7532989895
+  "ele2": {
+    "particle_opticalDepthQSR": 340.82229684726735,
+    "particle_position_x": 5.233059654483856e-06,
+    "particle_position_y": 4.781569085220371e-06,
+    "particle_position_z": 1.6293559425324337e-05,
+    "particle_momentum_x": 4.802611971470525e-18,
+    "particle_momentum_y": 4.3556825243407754e-18,
+    "particle_momentum_z": 2.41587659230925e-15,
+    "particle_weight": 2481137.860727036
   },
   "pho1": {
-    "particle_opticalDepthBW": 9769.560457892825,
-    "particle_position_x": 0.00013930069652083352,
-    "particle_position_y": 0.00014046343589264903,
-    "particle_position_z": 0.00046498427057467117,
+    "particle_opticalDepthBW": 9894.64959724129,
+    "particle_position_x": 0.00014191369823397644,
+    "particle_position_y": 0.00014347545392717968,
+    "particle_position_z": 0.00047442826029322116,
     "particle_momentum_x": 0.0,
     "particle_momentum_y": 0.0,
     "particle_momentum_z": 0.0,
-    "particle_weight": 59881947.934441805
+    "particle_weight": 61063948.68610941
   },
   "pho2": {
-    "particle_opticalDepthBW": 9986.841892753177,
-    "particle_position_x": 0.000143478808886347,
-    "particle_position_y": 0.00014235781306129165,
-    "particle_position_z": 0.0004760406640177353,
+    "particle_opticalDepthBW": 10292.1955840901,
+    "particle_position_x": 0.00014731892710321073,
+    "particle_position_y": 0.00014515617182809124,
+    "particle_position_z": 0.00048756513452074315,
     "particle_momentum_x": 0.0,
     "particle_momentum_y": 0.0,
     "particle_momentum_z": 0.0,
-    "particle_weight": 60412035.77539255
-  },
-  "pos": {
-    "particle_opticalDepthQSR": 1.0913353138611541,
-    "particle_position_x": 1.297022467502343e-08,
-    "particle_position_y": 6.533761995664651e-09,
-    "particle_position_z": 5.270388297177834e-08,
-    "particle_momentum_x": 9.474214808949084e-21,
-    "particle_momentum_y": 3.664691658234474e-21,
-    "particle_momentum_z": 5.559027877774493e-20,
-    "particle_weight": 4777.4993625771185
+    "particle_weight": 62299636.622087024
   },
-  "pos_nlbw1": {
-    "particle_opticalDepthQSR": 353.1921715878159,
-    "particle_position_x": 5.413239275792192e-06,
-    "particle_position_y": 5.599083716468022e-06,
-    "particle_position_z": 1.729829030925741e-05,
-    "particle_momentum_x": 5.130129014453181e-18,
-    "particle_momentum_y": 5.151236001734455e-18,
-    "particle_momentum_z": 2.7592435303859743e-15,
-    "particle_weight": 2298018.034763882
+  "pos1": {
+    "particle_opticalDepthQSR": 387.7441392212553,
+    "particle_position_x": 5.1462880118803425e-06,
+    "particle_position_y": 5.2613832293016684e-06,
+    "particle_position_z": 1.7054223425917483e-05,
+    "particle_momentum_x": 4.6437665862693495e-18,
+    "particle_momentum_y": 4.761862836969051e-18,
+    "particle_momentum_z": 2.3776996599289627e-15,
+    "particle_weight": 2625121.7841375084
   },
-  "pos_nlbw2": {
-    "particle_opticalDepthQSR": 350.113008933412,
-    "particle_position_x": 4.925666132330728e-06,
-    "particle_position_y": 4.677928261139806e-06,
-    "particle_position_z": 1.7528432992559524e-05,
-    "particle_momentum_x": 5.3410973872204896e-18,
-    "particle_momentum_y": 5.0088980554153576e-18,
-    "particle_momentum_z": 2.4803849492900593e-15,
-    "particle_weight": 2379194.375363683
+  "pos2": {
+    "particle_opticalDepthQSR": 361.943365907597,
+    "particle_position_x": 4.969019565031149e-06,
+    "particle_position_y": 4.361394970806125e-06,
+    "particle_position_z": 1.7413304358612675e-05,
+    "particle_momentum_x": 5.6348322786528905e-18,
+    "particle_momentum_y": 4.8171439953214205e-18,
+    "particle_momentum_z": 2.1937254860708963e-15,
+    "particle_weight": 2529794.7740602638
   }
 }