Modeling co-deposited samples

[m2patino]

I’m trying to reproduce the experimental results from Baldwin et al. [NME 23 (2020) 100743], where the trap energies and concentrations were used as input for FESTIM.

For codeposited samples, is it possible to set the D concentration in traps as an initial condition (and just model the TDS phase) or must the implantation phase also be modeled? Please note that I get all concentrations = 0 if I try to use set initial concentrations in traps, and the following TDS curves if I set the concentrations as a volumetric source or by implantation.

[remidm]

[posted here after private conversations]

The issue with your code is that you had:

my_model.initial_condition = [.....]

instead of

my_model.initial_conditions = [.....]

When using the wrong (singular) syntax, the initial conditions are simply ignored, which explain why the concentration are all zero.