Ionized electron momentum using Tomassini's model

src/particles/plasma/PlasmaParticleContainer.cpp

@@ -629,7 +629,7 @@ LaserIonization (const int islice,
 
         // This kernel supports multiple deposition orders (0, 1, 2, 3) at compile time.
         // It calculates the momentum of ionized electrons based on equations (B8) and (B9)
-        // from the Massimo (2020) article and equation (14) from the C. Schroeder (2014) article.
+        // from the Massimo (2020) article and equation (12) from the P. Tomassini (2021) article.
         // Additionally, it computes the energy of emitted electrons and assigns their properties
         // (momentum, position, etc.) to newly created electrons in the plasma container.
         amrex::AnyCTO(
@@ -675,8 +675,8 @@ LaserIonization (const int islice,
                     Ep *= E0;
                     amrex::Real delta = std::sqrt(Ep) * laser_dp_prefactor[ion_lev_loc];
                     amrex::Real delta2 = delta * delta;
-                    amrex::Real width_p = amrex::abs(A) * delta * (1._rt - (3._rt/4._rt) * delta2
-                        - (3._rt/2._rt) * delta2 + laser_dp_second_prefactor[ion_lev_loc] * delta2);
+                    amrex::Real width_p = amrex::abs(A) * delta * (1._rt - (7._rt/4._rt) * delta2
+                        + laser_dp_second_prefactor[ion_lev_loc] * delta2);
                     amrex::Real p_pol = amrex::RandomNormal(0.0, width_p, engine);
                     ux = p_pol; // linear polarization is assumed along x.
                     uz = (amrex::abs(A * A) * 0.25_rt + p_pol * p_pol * 0.5_rt);

src/particles/plasma/PlasmaParticleContainerInit.cpp

@@ -453,7 +453,7 @@ InitIonizationModule (const amrex::Geometry& geom, const amrex::Real background_
         h_adk_exp_prefactor[i] = -2./3. * std::pow( Uion/UH,3./2.) * Ea;
         h_laser_adk_prefactor[i] = (3./MathConst::pi) * std::pow(Uion/UH, -3./2.) / Ea;
         h_laser_dp_prefactor[i] = std::sqrt(3./2./Ea) * std::pow(UH/Uion, 3./4.);
-        h_laser_dp_second_prefactor[i] = 2 * ion_atomic_number * std::sqrt(UH/Uion) - 1./2.;
+        h_laser_dp_second_prefactor[i] = ion_atomic_number * std::sqrt(UH/Uion) - 1./2.;
     }
 
     amrex::Gpu::copy(amrex::Gpu::hostToDevice,

tests/checksum/benchmarks_json/laser_ionization.1Rank.json

@@ -1,23 +1,23 @@
 {
   "lev=0": {
-    "Bx": 13.149192792269,
-    "By": 9.6319466943698,
-    "Bz": 2214.7586533665,
+    "Bx": 12.895401136926,
+    "By": 9.7408143220707,
+    "Bz": 2091.4259740328,
     "ExmBy": 0.0,
     "EypBx": 0.0,
-    "Ez": 779215198153.58,
+    "Ez": 735839981397.02,
     "Psi": 0.0,
-    "Sx": 1114892448780.4,
-    "Sy": 1516523075637.6,
+    "Sx": 1120131620270.7,
+    "Sy": 1484350002437.1,
     "aabs": 0.46877495720552,
-    "chi": 265993671876340.0,
-    "jx": 704455285849490.0,
+    "chi": 265993671878070.0,
+    "jx": 665650594373930.0,
     "jx_beam": 0.0,
     "jy": 0.0,
     "jy_beam": 0.0,
     "jz_beam": 0.0,
     "laserEnvelope": 74.898701643489,
-    "rho_elec": 1202832873.9999,
+    "rho_elec": 1202832466.3739,
     "rho_ion": 1202806069.8845,
     "rhomjz": 0.0
   }