Add rough and dirty order of magnitude estimate of second order corre…

…ctions of QED processes
NeilZaim · Sep 11, 2023 · d794662 · d794662
1 parent 73c9529
commit d794662
Show file tree

Hide file tree

Showing 2 changed files with 10 additions and 0 deletions.
diff --git a/Source/Particles/ElementaryProcess/QEDInternals/BreitWheelerEngineWrapper.H b/Source/Particles/ElementaryProcess/QEDInternals/BreitWheelerEngineWrapper.H
@@ -173,12 +173,17 @@ public:
             chi_phot < m_bw_minimum_chi_phot)
             return 0;
 
+        const amrex::ParticleReal opt_depth_old = opt_depth;
+
         const auto is_out = pxr_bw::evolve_optical_depth<
             amrex::ParticleReal,
             BW_dndt_table_view,
             pxr_p::unit_system::SI>(
                 energy, chi_phot, dt, opt_depth, m_table_view);
 
+        constexpr amrex::ParticleReal two_thirds = amrex::ParticleReal(2.)/amrex::ParticleReal(3.);
+        opt_depth = opt_depth - (opt_depth_old-opt_depth)*PhysConst::alpha*std::pow(chi_phot,two_thirds);
+
         return is_out;
     }
 

diff --git a/Source/Particles/ElementaryProcess/QEDInternals/QuantumSyncEngineWrapper.H b/Source/Particles/ElementaryProcess/QEDInternals/QuantumSyncEngineWrapper.H
@@ -156,12 +156,17 @@ public:
         if (chi_part < m_qs_minimum_chi_part)
             return 0;
 
+        const amrex::ParticleReal opt_depth_old = opt_depth;
+
         const auto is_out = pxr_qs::evolve_optical_depth<
             amrex::ParticleReal,
             QS_dndt_table_view,
             pxr_p::unit_system::SI>(
                 energy, chi_part, dt, opt_depth, m_table_view);
 
+        constexpr amrex::ParticleReal two_thirds = amrex::ParticleReal(2.)/amrex::ParticleReal(3.);
+        opt_depth = opt_depth - (opt_depth_old-opt_depth)*PhysConst::alpha*std::pow(chi_part,two_thirds);
+
         return is_out;
     }