# Modeling co-deposited samples

**URL:** <https://festim.discourse.group/t/modeling-co-deposited-samples/28>\
**Category:** User Support\
**Created:** [January 25, 2024, 8:32pm UTC](https://festim.discourse.group/t/modeling-co-deposited-samples/28 "2024-01-25T20:32:47Z")\
**Posts on this page:** 2\
**Page:** 1

<div class="post-metadata">

**Author:** ![m2patino](https://avatars.discourse-cdn.com/v4/letter/m/9e8a1a/32.png) [@m2patino](https://festim.discourse.group/u/m2patino)\
**Post date:** [January 25, 2024, 8:32pm UTC](https://festim.discourse.group/t/modeling-co-deposited-samples/28/1 "2024-01-25T20:32:47Z")

</div>

I’m trying to reproduce the experimental results from [Baldwin et al. [NME 23 (2020) 100743]](https://www.sciencedirect.com/science/article/pii/S2352179120300193), 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.

 ![FESTIM_Baldwin_implantation](https://global.discourse-cdn.com/free1/uploads/festim/original/1X/e7a0b4ed64b2cf8bac223c90ae4f7c00833be6be.png)  
 ![FESTIM_Baldwin_sources](https://global.discourse-cdn.com/free1/uploads/festim/original/1X/4e73c9c7340ee3cd4761db6c30fde3d363975292.png)

---

<div class="post-metadata">

**Author:** ![remidm](https://yyz2.discourse-cdn.com/free1/user_avatar/festim.discourse.group/remidm/32/4_2.png) [@remidm](https://festim.discourse.group/u/remidm)\
**Post date:** [January 25, 2024, 8:54pm UTC](https://festim.discourse.group/t/modeling-co-deposited-samples/28/2 "2024-01-25T20:54:38Z")

</div>

[posted here after private conversations]

The issue with your code is that you had:

```python
my_model.initial_condition = [.....]

```

instead of

```python
my_model.initial_conditions = [.....]

```

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