diff --git a/results/soret/param_D.txt b/results/soret/param_D.txt
new file mode 100644
index 000000000..458f769a6
--- /dev/null
+++ b/results/soret/param_D.txt
@@ -0,0 +1,19 @@
+netParam=0 1 0 0 0
+networkFile=network6.h5
+dimensions=1
+vizHandler=dummy
+process=diff
+material=W100
+fluxDepthProfileFilePath=tridyn_D.dat
+petscArgs=-helium_retention -ts_adapt_clip 0.1,1.1 -ts_max_steps 100000000 -start_stop 100.0 -ts_final_time 25000.0 -ts_adapt_dt_max 5.0e-3 -ts_dt 1e-10 -ts_exact_final_time matchstep -ts_monitor -fieldsplit_0_pc_type sor -ts_max_snes_failures -1 -pc_fieldsplit_detect_coupling -snes_force_iteration -fieldsplit_1_pc_type redundant -ts_adapt_wnormtype INFINITY -pc_type fieldsplit
+flux=1.0e3
+gridType=geometric
+gridParam=500 1.0292315
+tempGridPower=2.5
+tempHandler=heat
+tempParam=6.0e-12 400
+voidPortion=0
+boundary=1 1
+xBCType=robin
+perfHandler=dummy
+interfaceLoc=-190000
diff --git a/results/soret/param_He.txt b/results/soret/param_He.txt
new file mode 100644
index 000000000..82c30fe50
--- /dev/null
+++ b/results/soret/param_He.txt
@@ -0,0 +1,18 @@
+netParam=1 0 0 0 0
+dimensions=1
+vizHandler=dummy
+process=diff
+material=W100
+fluxDepthProfileFilePath=tridyn_He.dat
+petscArgs=-helium_retention -ts_adapt_clip 0.1,1.1 -ts_max_steps 100000000 -start_stop 100.0 -ts_final_time 10000.0 -ts_adapt_dt_max 5.0e-3 -ts_dt 1e-10 -ts_exact_final_time matchstep -ts_monitor -fieldsplit_0_pc_type sor -ts_max_snes_failures -1 -pc_fieldsplit_detect_coupling -snes_force_iteration -fieldsplit_1_pc_type redundant -ts_adapt_wnormtype INFINITY -pc_type fieldsplit
+flux=1.0e3
+gridType=geometric
+gridParam=500 1.0292315
+tempGridPower=2.5
+tempHandler=heat
+tempParam=6.0e-12 400
+voidPortion=0
+boundary=1 1
+xBCType=robin
+perfHandler=dummy
+interfaceLoc=-15000000
diff --git a/results/soret/param_T.txt b/results/soret/param_T.txt
new file mode 100644
index 000000000..530702cc3
--- /dev/null
+++ b/results/soret/param_T.txt
@@ -0,0 +1,18 @@
+netParam=0 0 1 0 0
+dimensions=1
+vizHandler=dummy
+process=diff
+material=W100
+fluxDepthProfileFilePath=tridyn_T.dat
+petscArgs=-helium_retention -ts_adapt_clip 0.1,1.1 -ts_max_steps 100000000 -start_stop 100.0 -ts_final_time 25000.0 -ts_adapt_dt_max 5.0e-3 -ts_dt 1e-10 -ts_exact_final_time matchstep -ts_monitor -fieldsplit_0_pc_type sor -ts_max_snes_failures -1 -pc_fieldsplit_detect_coupling -snes_force_iteration -fieldsplit_1_pc_type redundant -ts_adapt_wnormtype INFINITY -pc_type fieldsplit
+flux=1.0e3
+gridType=geometric
+gridParam=500 1.0292315
+tempGridPower=2.5
+tempHandler=heat
+tempParam=6.0e-12 400
+voidPortion=0
+boundary=1 1
+xBCType=robin
+perfHandler=dummy
+interfaceLoc=-190000
diff --git a/results/soret/readme.txt b/results/soret/readme.txt
new file mode 100644
index 000000000..3d49de5c1
--- /dev/null
+++ b/results/soret/readme.txt
@@ -0,0 +1,8 @@
+readme.txt for Xolotl results
+20230623 Sophie Blondel
+
+This folder contains parameter files to reproduce results from https://doi.org/10.1088/1741-4326/acda4a . 
+- To use a damaged material change interfaceLoc=190
+- To add the Soret effect change process=diff soret
+
+The Xolotl commit is ed73cc62907257795e54480e8a8807d7fdc8ef41 .
diff --git a/results/soret/tridyn_D.dat b/results/soret/tridyn_D.dat
new file mode 100644
index 000000000..59a9bcb95
--- /dev/null
+++ b/results/soret/tridyn_D.dat
@@ -0,0 +1,2 @@
+D 1 1.0
+1059679.71284 247566.23629 3923948.05995 -3604314.08878 1568101.0087 -423999.229795 78064.4196404 -10206.2571617 968.141178442 -67.2494306824 3.41738425279 -0.125504952707 0.00324062967035 -5.57747749278e-05 5.7423497213e-07 -2.67442715008e-09 27.6567764334
diff --git a/results/soret/tridyn_He.dat b/results/soret/tridyn_He.dat
new file mode 100644
index 000000000..3600f7f78
--- /dev/null
+++ b/results/soret/tridyn_He.dat
@@ -0,0 +1,2 @@
+He 1 1.0
+-1365.18158572 24224.8878954 34639.344914 -72648.076295 57905.2884887 -27303.8598703 8483.78463815 -1825.31802219 279.6091276 -30.9339783622 2.47857639273 -0.142474884768 0.00572702746647 -0.000152834171315 2.43250628305e-06 -1.74741304962e-08 18.0400194055
diff --git a/results/soret/tridyn_T.dat b/results/soret/tridyn_T.dat
new file mode 100644
index 000000000..316db1d47
--- /dev/null
+++ b/results/soret/tridyn_T.dat
@@ -0,0 +1,2 @@
+T 1 1.0
+995770.337156 437447.115495 3451066.19453 -3215275.25638 1388081.01183 -370284.251218 67162.3853584 -8653.74428773 809.872638634 -55.5742523656 2.79346715487 -0.101598044627 0.00260063637513 -4.4412644642e-05 4.54064037876e-07 -2.10141070855e-09 27.6681164923