Radiated energy calculations with Urad

[jjl364]

Description

I am using the non-linear Compton model in SMILEI to determine the yield of photon radiation in a laser-plasma interaction. I only have one species (electrons) that are allowed to radiate, and I have been tracking the total radiated energy using the U_rad parameter from the scalar file.

I am finding that this U_rad number is underestimating the actual radiation energy by a significant margin. The discrepancy grows when the total laser energy increases. I get the numbers I expect by manually integrating the photon spectra for each output of photon RAW particles. Manual integration still might lose some photons that are created and lost to simulation walls in-between outputs, but it does set a lower bound for conversion efficiency and is much higher than U_rad suggests. Is there an easy fix for this?

Thanks for the help!

Parameters

• Smilei version 5.0
• Input file & command script - attached
• Computing resources
  • C++ compiler 10.3.1
  • MPI version 2021.6.0
  • HDF5 version 1.14.0
  • Python version 3.9.12

[mccoys]

I think this is a question for @MickaelGrechX ?

[MickaelGrechX]

I'm surprised.
Can you please provide the input file and post-processing scripts you are using?
You mentioned they are attached but I can't find them!

[jjl364]

Attached are the files needed to recreate the issue I am having. The wire_radiation file is the main script, and the troubleshoot_* files are for analysis in MATLAB. The analysis scripts should recreate the plots shown below, where the blue curve (left) is the total energy in photons based on summing the energies of all photons ever created in the simulation (sorted by ID to avoid double-counts). The red curve (right) is what I get from the scalar.txt file under the column for Urad_eonWire - which is the only species allowed to produce radiation. All other species are not allowed to radiate. Thanks for the help! Josh  [scalars.txt](https://github.com/user-attachments/files/18804658/scalars.txt) [troubleshoot_photon_energy.txt](https://github.com/user-attachments/files/18804680/troubleshoot_photon_energy.txt) [troubleshoot_scalars.txt](https://github.com/user-attachments/files/18804682/troubleshoot_scalars.txt)

[mccoys]

Please reply on the github website. Not from your email. Attachments won't work